Genomics and museum specimens

Patterns of genetic variation reflect interactions among microevolutionary forces that vary in strength with changing demography. Here, patterns of variation within and among samples of the mouthbrooding gafftopsail catfish (Bagre marinus, Family Ariidae) captured in the U.S. Atlantic and throughout the Gulf of Mexico were analyzed using genomics to generate neutral and non-neutral SNP data sets. Because genomic resources are lacking for ariids, linkage disequilibrium network analysis was used to examine patterns of putatively adaptive variation. Finally, historical demographic parameters were estimated from site frequency spectra. The results show four differentiated groups, corresponding to the (1) U.S. Atlantic, and the (2) northeastern, (3) northwestern, and (4) southern Gulf of Mexico. The non-neutral data presented two contrasting signals of structure, one due to increases in diversity moving west to east and north to south, and another to increased heterozygosity in the Atlantic. Demographic analysis suggested that recently reduced long-term effective population size in the Atlantic is likely an important driver of patterns of genetic variation and is consistent with a known reduction in population size potentially due to an epizootic. Overall, patterns of genetic variation resemble that of other fishes that use the same estuarine habitats as nurseries, regardless of the presence/absence of a larval phase, supporting the idea that adult/juvenile behavior and habitat are important predictors of contemporary patterns of genetic structure.

Response to Cavalli-Sforza Interview: [Human Biology 82(3):245-266 (June 2010)

Key insights into the evolutionary history of recently extinct or critically endangered species can be obtained through analysis of genomic data collected using high-throughput sequencing and ancient DNA from museum specimens, particularly where specimens are rare. For instance, the evolutionary history of the critically endangered Puebla deer mouse, Peromyscus mekisturus, remains unclear due to discordance between morphological and molecular phylogenetic analyses. However, previous molecular analyses were based on PCR and Sanger sequencing of only a few mitochondrial genes. Here, we used ancient DNA from historical museum specimens followed by target enrichment and high-throughput sequencing of several thousand nuclear ultraconserved elements and whole mitochondrial genomes to test the validity of the previous phylogenetic placement of P. mekisturus. Based on UCEs and mitogenomes, our results revealed that P. mekisturus forms a well-supported distinct lineage outside the clade containing all other members of the Peromyscus melanophrys group. Additionally, the mitogenome phylogeny further supports the placement of P. mekisturus as the sister species of the genus Reithrodontomys. This conflicts with the previous mtDNA phylogenetic reconstruction, in which P. mekisturus was nested within the species P. melanophrys. Our study demonstrates that high-throughput sequencing of ancient DNA, appropriately controlling for contamination and degradation, can provide a robust resolution of the evolutionary history and taxonomic status of species for which few or no modern genetic samples exist. In light of our results and pending further analysis with denser taxon sampling and the addition of morphological data, a re-evaluation of the taxonomy and conservation management plans of P. mekisturus is needed to ensure that the evolutionary distinctiveness of this species is recognized in future conservation efforts.

Laboratory techniques for high‐throughput sequencing have enhanced our ability to generate DNA sequence data from millions of natural history specimens collected prior to the molecular era, but remain poorly tested at shallower evolutionary time scales. Hybridization capture using restriction site‐associated DNA probes (hyRAD) is a recently developed method for population genomics with museum specimens. The hyRAD method employs fragments produced in a restriction site‐associated double digestion as the basis for probes that capture orthologous loci in samples of interest. While promising in that it does not require a reference genome, hyRAD has yet to be applied across study systems in independent laboratories. Here, we provide an independent assessment of the effectiveness of hyRAD on both fresh avian tissue and dried tissue from museum specimens up to 140 years old and investigate how variable quantities of input DNA affect sequencing, assembly, and population genetic inference. We present a modified bench protocol and bioinformatics pipeline, including three steps for detection and removal of microbial and mitochondrial DNA contaminants. We confirm that hyRAD is an effective tool for sampling thousands of orthologous SNPs from historic museum specimens to describe phylogeographic patterns. We find that modern DNA performs significantly better than historical DNA better during sequencing but that assembly performance is largely equivalent. We also find that the quantity of input DNA predicts %GC content of assembled contiguous sequences, suggesting PCR bias. We caution against sampling schemes that include taxonomic or geographic autocorrelation across modern and historic samples.

Past climate change is known to have strongly impacted current patterns of genetic variation of animals and plants in Europe. However, ecological factors also have the potential to influence demographic history and thus patterns of genetic variation. In this study, we investigated the impact of past climate, and also the potential impact of host plant species abundance, on intraspecific genetic variation in three codistributed and related specialized solitary bees of the genus Melitta with very similar life history traits and dispersal capacities. We sequenced five independent loci in samples collected from the three species. Our analyses revealed that the species associated with the most abundant host plant species (Melitta leporina) displays unusually high genetic variation, to an extent that is seldom reported in phylogeographic studies of animals and plants. This suggests a potential role of food resource abundance in determining current patterns of genetic variation in specialized herbivorous insects. Patterns of genetic variation in the two other species indicated lower overall levels of diversity, and that M.nigricans could have experienced a recent range expansion. Ecological niche modelling of the three Melitta species and their main host plant species suggested a strong reduction in range size during the last glacial maximum. Comparing observed sequence data with data simulated using spatially explicit models of coalescence suggests that M.leporina recovered a range and population size close to their current levels at the end of the last glaciation, and confirms recent range expansion as the most likely scenario for M.nigricans. Overall, this study illustrates that both demographic history and ecological factors may have contributed to shape current phylogeographic patterns.

BackgroundPleistocene climatic oscillations and historical geological events may both influence current patterns of genetic variation, and the species in southern China that faced unique climatic and topographical events have complex evolutionary histories. However, the relative contributions of climatic oscillations and geographical events to the genetic variation of these species remain undetermined. To investigate patterns of genetic variation and to test the hypotheses about the factors that shaped the distribution of this genetic variation in species of southern China, mitochondrial genes (cytochrome b and NADH dehydrogenase subunit 2) and nine microsatellite loci of the Omei tree frog (Rhacophorus omeimontis) were amplified in this study.ResultsThe genetic diversity in the populations of R. omeimontis was high. The phylogenetic trees reconstructed from the mitochondrial DNA (mtDNA) haplotypes and the Bayesian genetic clustering analysis based on microsatellite data both revealed that all populations were divided into three lineages (SC, HG and YN). The two most recent splitting events among the lineages coincided with recent geological events (including the intense uplift of the Qinghai-Tibet Plateau, QTP and the subsequent movements of the Yun-Gui Plateau, YGP) and the Pleistocene glaciations. Significant expansion signals were not detected in mismatch analyses or neutrality tests. And the effective population size of each lineage was stable during the Pleistocene.ConclusionsBased on the results of this study, complex geological events (the recent dramatic uplift of the QTP and the subsequent movements of the YGP) and the Pleistocene glaciations were apparent drivers of the rapid divergence of the R. omeimontis lineages. Each diverged lineages survived in situ with limited gene exchanges, and the stable demographics of lineages indicate that the Pleistocene climatic oscillations were inconsequential for this species. The analysis of genetic variation in populations of R. omeimontis contributes to the understanding of the effects of changes in climate and of geographical events on the dynamic development of contemporary patterns of genetic variation in the species of southern China.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0572-1) contains supplementary material, which is available to authorized users.

The relation between populations structure variables and regional diversity was investigated among the thirteen endogamous populations of Puria division, Bihar, who belonged to different castes and ethnic groups with various migratory origins The phenotypic and gene frequencies of four common loci (ABO, Rh, PTC and Color Blindness) show wide differences between populations. In general B is more common than A, but in four populations, where group O is predominant. Rh negative frequency is the highest among Brahmin and non-taster gene frequency is the least among Chamars. Among the castes, Brahmins show a maximum value of the average heterozygosity per locus (37.6%). The extent of genetic differentiation between populations (as measured by Gst) is 1.45%. The genetic distance analysis and the cluster analysis based on Nei's and Harpending and Jenkins' measure, show agreement between the observed levels of genetic diversity in line with the ethno-historical migration. This is also in agreement with the caste hierarchical structure of the regional populations. Principal component analysis based on similarity index ahows the same pattern. The apportionment of genetic diversity analysis supports the hypothesis of differential migratory origins for these populations. The Harpending and Ward model used to investigate the association between population structure and genetic diversity fits the data in predicting the heterozygosity levels in relation to the long range gene flow, especially from the eastern migrants. Introduction Understanding the extent and pattern of genetic variation and explaining the plausible causative factors among human populations are the major objectives of anthropological genetics. Previous studies have approached the problem by model-free methods, which attempt to infer possible influence of population structure variables by distance measures and dendrograms by comparing the variation between populations. Recently model-bound methods have been developed for using polygenic quantitative traits to test the hypothesis of population structure variables (especially for gene flow) on the between-population variation and genetic distance measures (Chakroborty, 1990). The application of model-based studies in some populations have revealed the extent of the influence of migration and other factors on the genetic differentiation between populations, as for example, among the Jirels of Nepal (Williams-- Blangero, 1990). In this study the model-bound approach has been considered among regional populations of known historical migration background to detect the influence of migration on the genetic differentiation. Studies among Indian populations are mostly based on model-- free methods. And the results of these studies have shown wide genetic diversity characterized by gradients of gene frequencies with uneven distribution of very high or Low frequency of a few characters in certain populations (Malhotra and Vasulu, 1993: Majumder and Mukherjee, 1993 and etc.). Attempts to explain the causative factors of this uneven pattern of distribution have shown that geographical contiguity and ethno-historical factors are the primary explanatory variables, although other associated population structure variables also play an important role. However, the relative importance of these explanatory variables differ between populations, so that to discern a general trend, studies covering many populaitons over a region and small scale local studies, especially specific populations, are necessary. In this context, regional level studies aiming to discern the genetic diversity at the level of local populations (e.g. tribal, ethnic or isolated) is expected to provide a better explanation of the causative factors involved. Indeed such regional studies in some specific populations have revealed the role of population structure in explaining the pattern of genetic variation. For example, the pattern of affinity between the two neighboring populations of Himachal Pradesh, especially, Kamte and Kotis, are due to admixture with the neighboring Tibetans (Papiha, 1985). …

The resolution of the Tree of Life has accelerated with advances in DNA sequencing technology. To achieve dense taxon sampling, it is often necessary to obtain DNA from historical museum specimens to supplement modern genetic samples. However, DNA from historical material is generally degraded, which presents various challenges. In this study, we evaluated how the coverage at variant sites and missing data among historical and modern samples impacts phylogenomic inference. We explored these patterns in the brush-tongued parrots (lories and lorikeets) of Australasia by sampling ultraconserved elements in 105 taxa. Trees estimated with low coverage characters had several clades where relationships appeared to be influenced by whether the sample came from historical or modern specimens, which were not observed when more stringent filtering was applied. To assess if the topologies were affected by missing data, we performed an outlier analysis of sites and loci, and a data reduction approach where we excluded sites based on data completeness. Depending on the outlier test, 0.15% of total sites or 38% of loci were driving the topological differences among trees, and at these sites, historical samples had 10.9× more missing data than modern ones. In contrast, 70% data completeness was necessary to avoid spurious relationships. Predictive modeling found that outlier analysis scores were correlated with parsimony informative sites in the clades whose topologies changed the most by filtering. After accounting for biased loci and understanding the stability of relationships, we inferred a more robust phylogenetic hypothesis for lories and lorikeets.

Over the past 50 years, geneticists have made great strides in understanding how our species' evolutionary history gave rise to current patterns of human genetic diversity classically summarized by Lewontin in his 1972 paper, ‘The Apportionment of Human Diversity’. One evolutionary process that requires special attention in both population genetics and statistical genetics is admixture: gene flow between two or more previously separated source populations to form a new admixed population. The admixture process introduces ancestry-based structure into patterns of genetic variation within and between populations, which in turn influences the inference of demographic histories, identification of genetic targets of selection and prediction of complex traits. In this review, we outline some challenges for admixture population genetics, including limitations of applying methods designed for populations without recent admixture to the study of admixed populations. We highlight recent studies and methodological advances that aim to overcome such challenges, leveraging genomic signatures of admixture that occurred in the past tens of generations to gain insights into human history, natural selection and complex trait architecture.This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Recent developments in molecular genetics have given us much information about DNA markers and a detailed look at evolutionary causes of genetic variation in humans. Patterns of genetic variation within populations and genetic differences among populations is best understood in light of our species’ African origin 200,000 years ago and the subsequent dispersion of human populations throughout the world 70,000 to 100,000 years ago. Although an evolutionary focus on genetic variation using methods from the field of population genetics provides us with considerable insight regarding our history, older ideas based on traditional definitions of biological race have no explanatory value.

Conifers are an ecologically and economically important seed plant group that can provide significant insights into the evolution of land plants. Molecular phylogenetics has developed as an important approach in evolutionary studies, although there have been relatively few studies of conifers that employ large-scale data sourced from multiple nuclear genes. Target enrichment sequencing (target capture, exon capture, or Hyb-Seq) has developed as a key approach in modern phylogenomic studies. However, until now, there has been no bait set that specifically targets the entire conifer clade. REMcon is a target sequence capture probe set intended for family- and species-level phylogenetic studies of conifers that target c. 100 single-copy nuclear loci. We tested the REMcon probe set using 69 species, including 44 conifer genera across six families and four other gymnosperm taxa, to evaluate the efficiency of target capture to efficiently generate comparable DNA sequence data across conifers. The recovery of target loci was high, with, on average, 94% of the targeted regions recovered across samples with high read coverage. A phylogenetic analysis of these data produced a well-supported topology that is consistent with the current understanding of relationships among conifers. The REMcon bait set will be useful in generating relatively large-scale nuclear data sets consistently for any conifer lineage.

Ejaculates function as an integrated unit to ensure male fertility and paternity, can have a complex structure, and can experience multiple episodes of selection. Current studies on the evolution of ejaculates typically focus on phenotypic variation in sperm number, size, or related traits such as testes size as adaptations to postcopulatory male-male competition. However, the evolution of the integrated nature of ejaculate structure and function depends on genetic variation in and covariation between the component parts. Here we report a quantitative genetic study of the components of the ejaculate of the cockroach Nauphoeta cinerea, including those we know to experience postcopulatory sexual selection, in the context of functional integration of ejaculate characters. We use the patterns of genetic variation and covariation to infer how the integration of the functions of the ejaculate constrain and shape its evolution. Ejaculate components were highly variable, showed significant additive genetic variance, and moderate to high evolvability. The level of genetic variation in these characters, despite strong directional or truncating selection, may reflect the integration of multiple episodes of selection that occur in N. cinerea. There were few significant phenotypic correlations, but all the genetic correlations among ejaculate characters were significantly different from zero. The patterns of genetic variation and covariation suggest that there are important trade-offs among individual traits of the ejaculate and that evolution of ejaculate characteristics will not proceed unconstrained. Fully describing the genetic relationships among traits that perform as an integrated unit helps us understand how functional relationships constrain or facilitate the evolution of the complex structure that is the ejaculate.

Historic DNA data have the potential to identify phenotypic information otherwise invisible in the historical, archaeological and palaeontological record. In order to determine whether a single nucleotide polymorphism typing protocol based on single based extension (SNaPshot™) could produce reliable phenotypic data from historic samples, we genotyped three coat colour markers for a sample of historic Thoroughbred horses for which both phenotypic and correct genotypic information were known from pedigree information in the General Stud Book. Experimental results were consistent with the pedigrees in all cases. Thus we demonstrate that historic DNA techniques can produce reliable phenotypic information from museum specimens.

While fluctuating asymmetry (FA; small, random deviations from perfect symmetry in bilaterally symmetrical traits) is widely regarded as a proxy for environmental and genetic stress effects, empirical associations between FA and stress are often weak or heterogeneous among traits. A conceptually important source of heterogeneity in relationships with FA is variation in the selection history of the trait(s) under study, i.e. traits that experienced a (recent) history of directional change are predicted to be developmentally less stable, potentially through the loss of canalizing modifiers. Here we applied X-ray photography on museum specimens and live captures to test to what extent the magnitude of FA and FA-stress relationships covary with directional shifts in traits related to the flight apparatus of four East-African rainforest birds that underwent recent shifts in habitat quality and landscape connectivity. Both the magnitude and direction of phenotypic change varied among species, with some traits increasing in size while others decreased or maintained their original size. In three of the four species, traits that underwent larger directional changes were less strongly buffered against random perturbations during their development, and traits that increased in size over time developed more asymmetrically than those that decreased. As we believe that spurious relationships due to biased comparisons of historic (museum specimens) and current (field captures) samples can be ruled out, these results support the largely untested hypothesis that directional shifts may increase the sensitivity of developing traits to random perturbations of environmental or genetic origin.

We examined spatio-temporal genetic variation at 53 single nucleotide polymorphisms in anadromous Arctic char populations from Western Greenland, a region experiencing pronounced climate change. The study was based on contemporary and historical samples, the latter represented by DNA extracted from otoliths and scales from the 1950s–1960s. We investigated whether genetic population structure was temporarily stable or unstable, the latter due to relatively small spawning and nursery areas combined with a harsh Arctic environment. Furthermore, in order to evaluate the potential for adaptive responses and local adaptation we estimated effective population size (Ne) and migration rate (m). Temporal stability of genetic population structure was suggested, based on a hierarchical analysis of genetic differentiation showing much higher differentiation among samples from different populations (FCT = 0.091) than among temporal samples from the same populations (FSC = 0.01). This was further supported by a neighbor-joining tree and assignment of individuals that showed high contingency between historical and contemporary samples. Estimates of Ne were high (> 500) in three out of four populations, with a lower estimate in one population potentially reflecting fishing pressure or suboptimal environmental conditions. Estimates of m were in most cases low, ≤ 0.01. Ne and m estimates suggest a potential for adaptive responses and local adaptation. However, long generation time may also cause adaptive responses by microevolution to be unable to track climate change, especially considering the low migration rates that reduce potential evolutionary rescue by gene flow from populations better adapted to the altered environments.