Lack of gene flow: Narrow and dispersed differentiation islands in a triplet of Leptidea butterfly species

Genome scans in recently separated species can inform on molecular mechanisms and evolutionary processes driving divergence. Large‐scale polymorphism data from multiple species pairs are also key to investigate the repeatability of divergence—whether radiations tend to show parallel responses to similar selection pressures and/or underlying molecular forces. Here, we used whole‐genome resequencing data from six wood white (Leptidea sp.) butterfly populations, representing three closely related species with karyomorph variation, to infer the species' demographic history and characterize patterns of genomic diversity and differentiation. The analyses supported previously established species relationships, and there was no evidence for postdivergence gene flow. We identified significant intraspecific genetic structure, in particular between karyomorph extremes in the wood white (L. sinapis)—a species with a remarkable chromosome number cline across the distribution range. The genomic landscapes of differentiation were erratic, and outlier regions were narrow and dispersed. Highly differentiated (F_ST) regions generally had low genetic diversity (θ_π), but increased absolute divergence (D_XY) and excess of rare frequency variants (low Tajima's D). A minority of differentiation peaks were shared across species and population comparisons. However, highly differentiated regions contained genes with overrepresented functions related to metabolism, response to stimulus and cellular processes, indicating recurrent directional selection on a specific set of traits in all comparisons. In contrast to the majority of genome scans in recently diverged lineages, our data suggest that divergence landscapes in Leptidea have been shaped by directional selection and genetic drift rather than stable recombination landscapes and/or introgression.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single‐nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south‐western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.     

Exome analysis and functional classification of identified variants in racing Quarter Horses

The main objectives of this study were to identify and functionally classify SNPs and indels by exome sequencing of animals of the racing line of Quarter Horses. Based on the individual genomic estimated breeding values (GEBVs) for maximum speed index (SImax) obtained for 349 animals, two groups of 20 extreme animals were formed. Of these individuals, 20 animals with high GEBVs for SImax and 19 with low GEBVs for SImax had their exons and 5′ and 3′ UTRs sequenced. Considering SNPs and indels, 105 182 variants were identified in the expressed regions of the Quarter Horse genome. Of these, 72 166 variants were already known and 33 016 are new variants and were deposited in a database. The analysis of the set of gene variants significantly related (P_adjusted < 0.05) to extreme animals in conjunction with the predicted impact of the changes and the physiological role of protein product pointed to two candidate genes potentially related to racing performance: SLC3A1 on ECA15 and CCN6 on ECA10.     

High Variation in Single Nucleotide Polymorphisms (SNPs) and Insertions/Deletions (Indels) in the Highly Invasive Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Middle East-Asia Minor 1 (MEAM1)

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Middle East-Asia Minor 1 (MEAM1) is invasive and adaptive to varied environments throughout the world. The adaptability is closely related to genomic variation such as single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). In order to elucidate the feature of SNPs and indels in MEAM1, and reveal the association between SNPs/indels and adaptive capacity to various environments, a computational approach with QualitySNP was used to identify reliable SNPs and indels on the basis of 9110-expressed sequence tags of MEAM1 present in the NCBI database. There were 575 SNPs detected with a density of 10.1 SNPs/kb and 6.4 SNPs/contig. Also, 237 transitions (39.3%) and 366 transversions (60.7%) were obtained, where the ratio of transitions to transversions was 0.65:1. In addition, 581 indels with a density of 14.1 indels/kb and 9.2 indels/contig were detected. Collectively, it showed that invasive MEAM1 has high SNPs density, and higher SNPs percentage than non-invasive B. tabaci species. A high SNPs density/percentage in MEAM1 yielded a high genomic variation that might have allowed it to adapt to varied environments, which provides some support to understand the invasive nature of MEAM1 at the genomic level. High levels of genomic variation are implicated in the level of adaptive capacity and invasive species are thought to exhibit higher levels of adaptive capacity than non-invasive species.     

BGVD:An Integrated Database for Bovine Sequencing Variations and Selective Signatures

Next-generation sequencing has yielded a vast amount of cattle genomic data for global characterization of population genetic diversity and identification of genomic regions under natural and artificial selection. However, efficient storage, querying, and visualization of such large datasets remain challenging. Here, we developed a comprehensive database, the Bovine Genome Variation Database (BGVD). It provides six main functionalities: gene search, variation search, genomic signature search, Genome Browser, alignment search tools, and the genome coordinate conversion tool. BGVD contains information on genomic variations comprising ~60.44 M SNPs, ~6.86 M indels, 76,634 CNV regions, and signatures of selective sweeps in 432 samples from modern cattle worldwide. Users can quickly retrieve distribution patterns of these variations for 54 cattle breeds through an interactive source of breed origin map, using a given gene symbol or genomic region for any of the three versions of the bovine reference genomes (ARS-UCD1.2, UMD3.1.1, and Btau 5.0.1). Signals of selection sweep are displayed as Manhattan plots and Genome Browser tracks. To further investigate and visualize the relationships between variants and signatures of selection, the Genome Browser integrates all variations, selection data, and resources, from NCBI, the UCSC Genome Browser, and Animal QTLdb. Collectively, all these features make the BGVD a useful archive for in-depth data mining and analyses of cattle biology and cattle breeding on a global scale. BGVD is publicly available at http://animal.nwsuaf.edu.cn/BosVar.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

Developing genomic resources in two Linum species via 454 pyrosequencing and genomic reduction

Recent advances in next-generation DNA sequencing (NGS) have enhanced the development of genomic resources such as contigs or single-nucleotide polymorphisms (SNPs) for evolutionary studies of a nonmodel species with a complex and unsequenced genome. This study presents an application of a NGS technique in combination with genomic reduction and advanced bioinformatics tools to identify contigs and SNPs from multiple samples of two Linum species. A full Roche 454 GS FLX run of 16 diverse Linum samples representing cultivated flax (Linum usitatissimum L.) and its wild progenitor (Linum bienne Mill.) generated approximately 1.6 million sequence reads with a total length of 498 Mbp. Application of the computational pipeline de novo identification of alleles identified 713 contigs and 1067 SNPs. A blast search revealed alignments of all 713 contigs with 491 existing Linum scaffolds and gene annotations associated with 512 contigs. Sanger sequencing confirmed 95% of 79 selected contigs and 94% of 272 SNPs and identified 211 new SNPs and 19 new indels. The scored 454 SNP data were highly imbalanced for assayed samples. These findings not only are useful for evolutionary studies of Linum species but also help to illustrate the utility of NGS technologies in SNP discovery for nonmodel organisms.     

A survey of genome‐wide single nucleotide polymorphisms through genome resequencing in the Périgord black truffle (Tuber melanosporum Vittad.)

The Périgord black truffle (Tuber melanosporum Vittad.), considered a gastronomic delicacy worldwide, is an ectomycorrhizal filamentous fungus that is ecologically important in Mediterranean French, Italian and Spanish woodlands. In this study, we developed a novel resource of single nucleotide polymorphisms (SNPs) for T. melanosporum using Illumina high‐throughput resequencing. The genome from six T. melanosporum geographical accessions was sequenced to a depth of approximately 20×. These geographical accessions were selected from different populations within the northern and southern regions of the geographical species distribution. Approximately 80% of the reads for each of the six resequenced geographical accessions mapped against the reference T. melanosporum genome assembly, estimating the core genome size of this organism to be approximately 110 Mbp. A total of 442 326 SNPs corresponding to 3540 SNPs/Mbps were identified as being included in all seven genomes. The SNPs occurred more frequently in repeated sequences (85%), although 4501 SNPs were also identified in the coding regions of 2587 genes. Using the ratio of nonsynonymous mutations per nonsynonymous site (pN) to synonymous mutations per synonymous site (pS) and Tajima's D index scanning the whole genome, we were able to identify genomic regions and genes potentially subjected to positive or purifying selection. The SNPs identified represent a valuable resource for future population genetics and genomics studies.     

GCAT|Panel,a comprehensive structural variant haplotype map of the Iberian population from high-coverage whole-genome sequencing

The combined analysis of haplotype panels with phenotype clinical cohorts is a common approach to explore the genetic architecture of human diseases. However, genetic studies are mainly based on single nucleotide variants (SNVs) and small insertions and deletions (indels). Here, we contribute to fill this gap by generating a dense haplotype map focused on the identification, characterization, and phasing of structural variants (SVs). By integrating multiple variant identification methods and Logistic Regression Models (LRMs), we present a catalogue of 35 431 441 variants, including 89 178 SVs (≥50 bp), 30 325 064 SNVs and 5 017 199 indels, across 785 Illumina high coverage (30x) whole-genomes from the Iberian GCAT Cohort, containing a median of 3.52M SNVs, 606 336 indels and 6393 SVs per individual. The haplotype panel is able to impute up to 14 360 728 SNVs/indels and 23 179 SVs, showing a 2.7-fold increase for SVs compared with available genetic variation panels. The value of this panel for SVs analysis is shown through an imputed rare Alu element located in a new locus associated with Mononeuritis of lower limb, a rare neuromuscular disease. This study represents the first deep characterization of genetic variation within the Iberian population and the first operational haplotype panel to systematically include the SVs into genome-wide genetic studies.     

Functional constraint and small insertions and deletions in the ENCODE regions of the human genome

Background  

We describe the distribution of indels in the 44 Encyclopedia of DNA Elements (ENCODE) regions (about 1% of the human genome) and evaluate the potential contributions of small insertion and deletion polymorphisms (indels) to human genetic variation. We relate indels to known genomic annotation features and measures of evolutionary constraint.     

Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly

Here we use whole-genome de novo assembly of second-generation sequencing reads to map structural variation (SV) in an Asian genome and an African genome. Our approach identifies small- and intermediate-size homozygous variants (1-50 kb) including insertions, deletions, inversions and their precise breakpoints, and in contrast to other methods, can resolve complex rearrangements. In total, we identified 277,243 SVs ranging in length from 1-23 kb. Validation using computational and experimental methods suggests that we achieve overall <6% false-positive rate and <10% false-negative rate in genomic regions that can be assembled, which outperforms other methods. Analysis of the SVs in the genomes of 106 individuals sequenced as part of the 1000 Genomes Project suggests that SVs account for a greater fraction of the diversity between individuals than do single-nucleotide polymorphisms (SNPs). These findings demonstrate that whole-genome de novo assembly is a feasible approach to deriving more comprehensive maps of genetic variation.     

Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species

High‐throughput DNA sequencing facilitates the analysis of large portions of the genome in nonmodel organisms, ensuring high accuracy of population genetic parameters. However, empirical studies evaluating the appropriate sample size for these kinds of studies are still scarce. In this study, we use double‐digest restriction‐associated DNA sequencing (ddRADseq) to recover thousands of single nucleotide polymorphisms (SNPs) for two physically isolated populations of Amphirrhox longifolia (Violaceae), a nonmodel plant species for which no reference genome is available. We used resampling techniques to construct simulated populations with a random subset of individuals and SNPs to determine how many individuals and biallelic markers should be sampled for accurate estimates of intra‐ and interpopulation genetic diversity. We identified 3646 and 4900 polymorphic SNPs for the two populations of A. longifolia, respectively. Our simulations show that, overall, a sample size greater than eight individuals has little impact on estimates of genetic diversity within A. longifolia populations, when 1000 SNPs or higher are used. Our results also show that even at a very small sample size (i.e. two individuals), accurate estimates of F_ST can be obtained with a large number of SNPs (≥1500). These results highlight the potential of high‐throughput genomic sequencing approaches to address questions related to evolutionary biology in nonmodel organisms. Furthermore, our findings also provide insights into the optimization of sampling strategies in the era of population genomics.     

Screening for signatures of selection of Tianzhu white yak using genome‐wide re‐sequencing

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai‐Tibet Plateau. Specific ecological and long‐term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome‐wide re‐sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and F_ST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high‐altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.     

Combining genome‐wide association study and FST‐based approaches to identify targets of Borrelia‐mediated selection in natural rodent hosts

Recent advances in high‐throughput sequencing technologies provide opportunities to gain novel insights into the genetic basis of phenotypic trait variation. Yet to date, progress in our understanding of genotype–phenotype associations in nonmodel organisms in general and natural vertebrate populations in particular has been hampered by small sample sizes typically available for wildlife populations and a resulting lack of statistical power, as well as a limited ability to control for false‐positive signals. Here we propose to combine a genome‐wide association study (GWAS) and F_ST‐based approach with population‐level replication to partly overcome these limitations. We present a case study in which we used this approach in combination with genotyping‐by‐sequencing (GBS) single nucleotide polymorphism (SNP) data to identify genomic regions associated with Borrelia afzelii resistance or susceptibility in the natural rodent host of this Lyme disease‐causing spirochete, the bank vole (Myodes glareolus). Using this combined approach we identified four consensus SNPs located in exonic regions of the genes Slc26a4, Tns3, Wscd1 and Espnl, which were significantly associated with the voles’ Borrelia infectious status within and across populations. Functional links between host responses to bacterial infections and most of these genes have previously been demonstrated in other rodent systems, making them promising new candidates for the study of evolutionary host responses to Borrelia emergence. Our approach is applicable to other systems and may facilitate the identification of genetic variants underlying disease resistance or susceptibility, as well as other ecologically relevant traits, in wildlife populations.     

Time-ordering japonica/geng genomes analysis indicates the importance of large structural variants in rice breeding

Temperate japonica/geng (GJ) rice yield has significantly improved due to intensive breeding efforts, dramatically enhancing global food security. However, little is known about the underlying genomic structural variations (SVs) responsible for this improvement. We compared 58 long-read assemblies comprising cultivated and wild rice species in the present study, revealing 156 319 SVs. The phylogenomic analysis based on the SV dataset detected the putatively selected region of GJ sub-populations. A significant portion of the detected SVs overlapped with genic regions were found to influence the expression of involved genes inside GJ assemblies. Integrating the SVs and causal genetic variants underlying agronomic traits into the analysis enables the precise identification of breeding signatures resulting from complex breeding histories aimed at stress tolerance, yield potential and quality improvement. Further, the results demonstrated genomic and genetic evidence that the SV in the promoter of LTG1 is accounting for chilling sensitivity, and the increased copy numbers of GNP1 were associated with positive effects on grain number. In summary, the current study provides genomic resources for retracing the properties of SVs-shaped agronomic traits during previous breeding procedures, which will assist future genetic, genomic and breeding research on rice.     

PRIMER NOTE. Using the incomplete genome of the ectomycorrhizal fungus Amanita bisporigera to identify molecular polymorphisms in the related Amanita phalloides

We describe the cross‐genomic isolation of 13 single nucleotide polymorphisms (SNPs) and one variable microsatellite from five loci for the death cap mushroom Amanita phalloides. Microsatellite repeats were identified by searching the partial Amanita bisporigera genome. Flanking primers were designed for 25 of these microsatellite loci and tested for cross‐amplification in A. phalloides. One locus contained an interrupted, compound microsatellite, and four loci contained one to six SNPs. These results demonstrate the usefulness of even an incomplete genome to identify molecular markers for population studies in nonmodel organisms.     

Multiple capacitors for natural genetic variation in Drosophila melanogaster

Cryptic genetic variation (CGV) or a standing genetic variation that is not ordinarily expressed as a phenotype is released when the robustness of organisms is impaired under environmental or genetic perturbations. Evolutionary capacitors modulate the amount of genetic variation exposed to natural selection and hidden cryptically; they have a fundamental effect on the evolvability of traits on evolutionary timescales. In this study, I have demonstrated the effects of multiple genomic regions of Drosophila melanogaster on CGV in wing shape. I examined the effects of 61 genomic deficiencies on quantitative and qualitative natural genetic variation in the wing shape of D. melanogaster. I have identified 10 genomic deficiencies that do not encompass a known candidate evolutionary capacitor, Hsp90, exposing natural CGV differently depending on the location of the deficiencies in the genome. Furthermore, five genomic deficiencies uncovered qualitative CGV in wing morphology. These findings suggest that CGV in wing shape of wild‐type D. melanogaster is regulated by multiple capacitors with divergent functions. Future analysis of genes encompassed by these genomic regions would help elucidate novel capacitor genes and better understand the general features of capacitors regarding natural genetic variation.