A growing variety of "genotype-by-sequencing" (GBS) methods use restriction enzymes and high throughput DNA sequencing to generate data for a subset of genomic loci, allowing the simultaneous ...discovery and genotyping of thousands of polymorphisms in a set of multiplexed samples. We evaluated a "double-digest" restriction-site associated DNA sequencing (ddRAD-seq) protocol by 1) comparing results for a zebra finch (Taeniopygia guttata) sample with in silico predictions from the zebra finch reference genome; 2) assessing data quality for a population sample of indigobirds (Vidua spp.); and 3) testing for consistent recovery of loci across multiple samples and sequencing runs. Comparison with in silico predictions revealed that 1) over 90% of predicted, single-copy loci in our targeted size range (178-328 bp) were recovered; 2) short restriction fragments (38-178 bp) were carried through the size selection step and sequenced at appreciable depth, generating unexpected but nonetheless useful data; 3) amplification bias favored shorter, GC-rich fragments, contributing to among locus variation in sequencing depth that was strongly correlated across samples; 4) our use of restriction enzymes with a GC-rich recognition sequence resulted in an up to four-fold overrepresentation of GC-rich portions of the genome; and 5) star activity (i.e., non-specific cutting) resulted in thousands of "extra" loci sequenced at low depth. Results for three species of indigobirds show that a common set of thousands of loci can be consistently recovered across both individual samples and sequencing runs. In a run with 46 samples, we genotyped 5,996 loci in all individuals and 9,833 loci in 42 or more individuals, resulting in <1% missing data for the larger data set. We compare our approach to similar methods and discuss the range of factors (fragment library preparation, natural genetic variation, bioinformatics) influencing the recovery of a consistent set of loci among samples.
Serpentine barrens represent extreme hazards for plant colonists. These sites are characterized by high porosity leading to drought, lack of essential mineral nutrients, and phytotoxic levels of ...metals. Nevertheless, nature forged populations adapted to these challenges. Here, we use a population-based evolutionary genomic approach coupled with elemental profiling to assess how autotetraploid Arabidopsis arenosa adapted to a multichallenge serpentine habitat in the Austrian Alps. We first demonstrate that serpentine-adapted plants exhibit dramatically altered elemental accumulation levels in common conditions, and then resequence 24 autotetraploid individuals from three populations to perform a genome scan. We find evidence for highly localized selective sweeps that point to a polygenic, multitrait basis for serpentine adaptation. Comparing our results to a previous study of independent serpentine colonizations in the closely related diploid Arabidopsis lyrata in the United Kingdom and United States, we find the highest levels of differentiation in 11 of the same loci, providing candidate alleles for mediating convergent evolution. This overlap between independent colonizations in different species suggests that a limited number of evolutionary strategies are suited to overcome the multiple challenges of serpentine adaptation. Interestingly, we detect footprints of selection in A. arenosa in the context of substantial gene flow from nearby off-serpentine populations of A. arenosa, as well as from A. lyrata. In several cases, quantitative tests of introgression indicate that some alleles exhibiting strong selective sweep signatures appear to have been introgressed from A. lyrata. This finding suggests that migrant alleles may have facilitated adaptation of A. arenosa to this multihazard environment.
Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci ...and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi‐locus data requires substantial genomic coverage. In this study, we used ddRAD‐seq to sample the genomes of five recently diverged, New World “mallards” (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD‐seq loci and 177 loci from the Z‐chromosome. In contrast to previous studies, the ddRAD‐seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD‐seq data were sufficient to identify previously unknown outlier regions across the Z‐chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.
Display omitted
•ddRAD-seq data yield several distinct categories of informative polymorphism data.•The utility of ddRAD-seq for phylogenetic analysis declines with genetic distance.•Concatenation ...and newer species tree methods produce largely consistent results.•Rapid speciation in brood parasitic Vidua results in lower gene tree concordance.
Genotype-by-sequencing (GBS) methods have revolutionized the field of molecular ecology, but their application in molecular phylogenetics remains somewhat limited. In addition, most phylogenetic studies based on large GBS data sets have relied on analyses of concatenated data rather than species tree methods that explicitly account for genealogical stochasticity among loci. We explored the utility of “double-digest” restriction site-associated DNA sequencing (ddRAD-seq) for phylogenetic analyses of the Lagonosticta firefinches (family Estrildidae) and the Vidua brood parasitic finches (family Viduidae). As expected, the number of homologous loci shared among samples was negatively correlated with genetic distance due to the accumulation of restriction site polymorphisms. Nonetheless, for each genus, we obtained data sets of ∼3000 loci shared in common among all samples, including a more distantly related outgroup taxon. For all samples combined, we obtained >1000 homologous loci despite ∼20my divergence between estrildid and parasitic finches. In addition to nucleotide polymorphisms, the ddRAD-seq data yielded large sets of indel and locus presence–absence polymorphisms, all of which had higher consistency indices than mtDNA sequence data in the context of concatenated parsimony analyses. Species tree methods, using individual gene trees or single nucleotide polymorphisms as input, generated results broadly consistent with analyses of concatenated data, particularly for Lagonosticta, which appears to have a well resolved, bifurcating history. Results for Vidua were also generally consistent across methods and data sets, although nodal support and results from different species tree methods were more variable. Lower gene tree congruence in Vidua is likely the result of its unique evolutionary history, which includes rapid speciation by host shift and occasional hybridization and introgression due to incomplete reproductive isolation. We conclude that ddRAD-seq is a cost-effective method for generating robust phylogenetic data sets, particularly for analyses of closely related species and genera.
Smoking negatively affects B cell function and immunoglobulin levels, but it is unclear if this immune dysfunction contributes to the risk of severe COVID-19 in smokers. We evaluated binding IgM, IgA ...and IgG antibodies to spike and receptor binding domain antigens, and used a pseudovirus assay to quantify neutralization titers in a set of 27 patients with severe COVID-19. We found no significant differences between binding and neutralization antibody responses for people with a smoking history and people who never smoked. High plasma viral load, but not antibody titers, was linked to an increased risk of death. Humoral immune dysfunction was not a major driver of severe COVID-19 in smokers.
Interspecific hybridization is recognized as an important process in the evolutionary dynamics of both speciation and the reversal of speciation. However, our understanding of the spatial and ...temporal patterns of hybridization that erode versus promote species boundaries is incomplete. The endangered, endemic koloa maoli (or Hawaiian duck, Anas wyvilliana) is thought to be threatened with genetic extinction through ongoing hybridization with an introduced congener, the feral mallard (A. platyrhynchos). We investigated spatial and temporal variation in hybrid prevalence in populations throughout the main Hawaiian Islands, using genomic data to characterize population structure of koloa, quantify the extent of hybridization, and compare hybrid proportions over time. To accomplish this, we genotyped 3,308 double‐digest restriction‐site‐associated DNA (ddRAD) loci in 425 putative koloa, mallards, and hybrids from populations across the main Hawaiian Islands. We found that despite a population decline in the last century, koloa genetic diversity is high. There were few hybrids on the island of Kauaʻi, home to the largest population of koloa. By contrast, we report that sampled populations outside of Kauaʻi can now be characterized as hybrid swarms, in that all individuals sampled were of mixed koloa × mallard ancestry. Further, there is some evidence that these swarms are stable over time. These findings demonstrate spatial variation in the extent and consequences of interspecific hybridization, and highlight how islands or island‐like systems with small population sizes may be especially prone to genetic extinction when met with a congener that is not reproductively isolated.
For more than 225 million y, all seed plants were woody trees, shrubs, or vines. Shortly after the origin of angiosperms ∼140 million y ago (MYA), the Nymphaeales (water lilies) became one of the ...first lineages to deviate from their ancestral, woody habit by losing the vascular cambium, the meristematic population of cells that produces secondary xylem (wood) and phloem. Many of the genes and gene families that regulate differentiation of secondary tissues also regulate the differentiation of primary xylem and phloem, which are produced by apical meristems and retained in nearly all seed plants. Here, we sequenced and assembled a draft genome of the water lily Nymphaea thermarum, an emerging system for the study of early flowering plant evolution, and compared it to genomes from other cambium-bearing and cambium-less lineages (e.g., monocots and Nelumbo). This revealed lineage-specific patterns of gene loss and divergence. Nymphaea is characterized by a significant contraction of the HD-ZIP III transcription factors, specifically loss of REVOLUTA, which influences cambial activity in other angiosperms. We also found the Nymphaea and monocot copies of cambium-associated CLE signaling peptides display unique substitutions at otherwise highly conserved amino acids. Nelumbo displays no obvious divergence in cambium-associated genes. The divergent genomic signatures of convergent loss of vascular cambium reveals that even pleiotropic genes can exhibit unique divergence patterns in association with independent events of trait loss. Our results shed light on the evolution of herbaceousness—one of the key biological innovations associated with the earliest phases of angiosperm evolution.
Speciation is a continuous and dynamic process, and studying organisms during the early stages of this process can aid in identifying speciation mechanisms. The mallard (Anas platyrhynchos) and ...Mexican duck (A. p. diazi) are two recently diverged taxa with a history of hybridization and controversial taxonomy. To understand their evolutionary history, we conducted genomic scans to characterize patterns of genetic diversity and divergence across the mitochondrial DNA (mtDNA) control region, 3523 autosomal loci and 172 Z‐linked sex chromosome loci. Between the two taxa, Z‐linked loci (ΦST = 0.088) were 5.2 times more differentiated than autosomal DNA (ΦST = 0.017) but comparable to mtDNA (ΦST = 0.092). This elevated Z differentiation deviated from neutral expectations inferred from simulated data that incorporated demographic history and differences in effective population sizes between marker types. Furthermore, 3% of Z‐linked loci, compared to <0.1% of autosomal loci, were detected as outlier loci under divergent selection with elevated relative (ΦST) and absolute (dXY) estimates of divergence. In contrast, the ratio of Z‐linked and autosomal differentiation among the seven Mexican duck sampling locations was close to 1:1 (ΦST = 0.018 for both markers). We conclude that between mallards and Mexican ducks, divergence at autosomal markers is largely neutral, whereas greater divergence on the Z chromosome (or some portions thereof) is likely the product of selection that has been important in speciation. Our results contribute to a growing body of literature indicating elevated divergence on the Z chromosome and its likely importance in avian speciation.
Large multi-locus data sets provide an unprecedented opportunity for detecting fine-scale population structure and identifying populations that warrant independent conservation and management. In ...this study, we used 3161 polymorphic loci obtained from double-digest restriction-site associated DNA sequencing (ddRAD-seq) to examine genetic differentiation between populations of mottled ducks (Anas fulvigula) and their close relative, the mallard (A. platyrhynchos). Florida (FL) and western gulf coast (WGC) mottled ducks group into discrete populations in ordination analyses and can be assigned unambiguously to their population of origin based on multi-locus genotypes. Moreover, we find limited evidence of recent admixture between populations, and the range gap separating them is a significant barrier to gene flow. We found no evidence of diversifying selection in the sample of loci we analyzed; instead, the genome-wide pattern of differentiation between mottled duck populations is consistent with historical isolation and divergence by neutral genetic drift. Thus, our results indicate that FL and WGC mottled ducks are discrete populations on separate evolutionary pathways and should be managed as independent conservation units, and perhaps recognized as distinct subspecies given observed differences in plumage and morphology. In contrast, the recently established South Atlantic coastal zone population is indistinguishable from the WGC population but shows evidence of a small amount of admixture from FL mottled ducks. Our study demonstrates a quantitative, operational approach for delineating discrete populations for conservation using ddRAD-seq data; multi-locus genotypes allow for a robust test of historical isolation and evolutionary independence on much shorter evolutionary timescales than the overly conservative standard of reciprocal monophyly in mtDNA.
Display omitted
•Florida (FL) and west gulf coast (WGC) mottled ducks are genetically discrete.•A sharp transition in genetic similarity is coincident with the FL-WGC range gap.•Estimates of genetic admixture and gene flow between FL and WGC are low.•Mallard introgression is likely higher than gene flow between populations.•FL and WGC mottled ducks should be considered as distinct conservation units.
In coevolutionary arms races, interacting species impose selection on each other, generating reciprocal adaptations and counter adaptations. This process is typically enhanced by genetic ...recombination and heterozygosity, but these sources of evolutionary novelty may be secondarily lost when uniparental inheritance evolves to ensure the integrity of sex-linked adaptations. We demonstrate that host-specific egg mimicry in the African cuckoo finch Anomalospiza imberbis is maternally inherited, confirming the validity of an almost century-old hypothesis. We further show that maternal inheritance not only underpins the mimicry of different host species but also additional mimetic diversification that approximates the range of polymorphic egg “signatures” that have evolved within host species as an escalated defense against parasitism. Thus, maternal inheritance has enabled the evolution and maintenance of nested levels of mimetic specialization in a single parasitic species. However, maternal inheritance and the lack of sexual recombination likely disadvantage cuckoo finches by stifling further adaptation in the ongoing arms races with their individual hosts, which we show have retained biparental inheritance of egg phenotypes. The inability to generate novel genetic combinations likely prevents cuckoo finches from mimicking certain host phenotypes that are currently favored by selection (e.g., the olive-green colored eggs laid by some tawny-flanked prinia, Prinia subflava, females). This illustrates an important cost of coding coevolved adaptations on the nonrecombining sex chromosome, which may impede further coevolutionary change by effectively reversing the advantages of sexual reproduction in antagonistic coevolution proposed by the Red Queen hypothesis.