Genotyping‐by‐sequencing for estimating relatedness in nonmodel organisms: Avoiding the trap of precise bias

There has been remarkably little attention to using the high resolution provided by genotyping‐by‐sequencing (i.e., RADseq and similar methods) for assessing relatedness in wildlife populations. A major hurdle is the genotyping error, especially allelic dropout, often found in this type of data that could lead to downward‐biased, yet precise, estimates of relatedness. Here, we assess the applicability of genotyping‐by‐sequencing for relatedness inferences given its relatively high genotyping error rate. Individuals of known relatedness were simulated under genotyping error, allelic dropout and missing data scenarios based on an empirical ddRAD data set, and their true relatedness was compared to that estimated by seven relatedness estimators. We found that an estimator chosen through such analyses can circumvent the influence of genotyping error, with the estimator of Ritland (Genetics Research, 67, 175) shown to be unaffected by allelic dropout and to be the most accurate when there is genotyping error. We also found that the choice of estimator should not rely solely on the strength of correlation between estimated and true relatedness as a strong correlation does not necessarily mean estimates are close to true relatedness. We also demonstrated how even a large SNP data set with genotyping error (allelic dropout or otherwise) or missing data still performs better than a perfectly genotyped microsatellite data set of tens of markers. The simulation‐based approach used here can be easily implemented by others on their own genotyping‐by‐sequencing data sets to confirm the most appropriate and powerful estimator for their data.     

A chromosome‐scale reference genome and genome‐wide genetic variations elucidate adaptation in yak

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu Kush ranges of the Himalayas. The yak genome was sequenced in 2012, but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyse it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied towards high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein‐coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here, we identified breed‐specific sequences and genes by whole‐genome alignment, which may facilitate yak breed identification.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

中国人群参考基因组及基因组变异图谱资源库

随着人类基因组计划和国际千人基因组计划的实施,已公开数百个中国人个体的全基因组数据。建立高精度的中国人群参考基因组序列,发现并解析中国人群特有的序列变异,是我国未来精准医学研究的基础。为满足未来精准医学研究中国人基因组数据持续增长的科学管理和深入研究的需求,中国科学院北京基因组研究所发展并建立了基于中国人群全基因组测序数据的虚拟中国人基因组数据库(Virtual Chinese Genome Database, VCGDB)和中国人群基因组变异数据库(Genome Variation Map, GVM),面向国内外用户提供数据检索、共享、下载和在线分析服务。本文重点介绍了这两个数据库的特点和功能,以及未来发展与应用前景,以期为中国人群参考基因组及基因组变异图谱资源库的推广使用、发展完善提供有益信息。     

Plant intraspecific functional trait variation is related to within‐habitat heterogeneity and genetic diversity in Trifolium montanum L.

Intraspecific trait variation (ITV), based on available genetic diversity, is one of the major means plant populations can respond to environmental variability. The study of functional trait variation and diversity has become popular in ecological research, for example, as a proxy for plant performance influencing fitness. Up to now, it is unclear which aspects of intraspecific functional trait variation (iFD_CV) can be attributed to the environment or genetics under natural conditions. Here, we examined 260 individuals from 13 locations of the rare (semi‐)dry calcareous grassland species Trifolium montanum L. in terms of iFD_CV, within‐habitat heterogeneity, and genetic diversity. The iFD_CV was assessed by measuring functional traits (releasing height, biomass, leaf area, specific leaf area, leaf dry matter content, F_v/F_m, performance index, stomatal pore surface, and stomatal pore area index). Abiotic within‐habitat heterogeneity was derived from altitude, slope exposure, slope, leaf area index, soil depth, and further soil factors. Based on microsatellites, we calculated expected heterozygosity (H_e) because it best‐explained, among other indices, iFD_CV. We performed multiple linear regression models quantifying relationships among iFD_CV, abiotic within‐habitat heterogeneity and genetic diversity, and also between separate functional traits and abiotic within‐habitat heterogeneity or genetic diversity. We found that abiotic within‐habitat heterogeneity influenced iFD_CV twice as strong compared to genetic diversity. Both aspects together explained 77% of variation in iFD_CV ( = .77, F_{2, 10} = 21.66, p < .001). The majority of functional traits (releasing height, biomass, specific leaf area, leaf dry matter content, F_v/F_m, and performance index) were related to abiotic habitat conditions indicating responses to environmental heterogeneity. In contrast, only morphology‐related functional traits (releasing height, biomass, and leaf area) were related to genetics. Our results suggest that both within‐habitat heterogeneity and genetic diversity affect iFD_CV and are thus crucial to consider when aiming to understand or predict changes of plant species performance under changing environmental conditions.     

Genetic differentiation and reduced genetic diversity at the northern range edge of two species with different dispersal modes

Theory predicts that genetic variation should be reduced at range margins, but empirical support is equivocal. Here, we used genotyping‐by‐sequencing technology to investigate genetic variation in central and marginal populations of two species in the marine gastropod genus Crepidula. These two species have different development and dispersal types and might therefore show different spatial patterns of genetic variation. Both allelic richness and the proportion of private alleles were highest in the most central populations of both species, and lower at the margin. The species with low dispersal, Crepidula convexa, showed high degrees of structure throughout the range that conform to the pattern found in previous studies using other molecular markers. The northernmost populations of the high‐dispersing species, Crepidula fornicata, are distinct from more central populations, although this species has been previously observed to have little genetic structure over much of its range. Although genetic diversity was significantly lower at the range margin, the absolute reduction in diversity observed with these genomewide markers was slight, and it is not yet known whether there are functional consequences for the marginal populations.     

Exploration of Genetic Variations through Single‐cell Whole‐genome Sequencing in the Model Ciliate Tetrahymena thermophila

Ciliates are unicellular eukaryotes with separate germline and somatic genomes and diverse life cycles, which make them a unique model to improve our understanding of population genetics through the detection of genetic variations. However, traditional sequencing methods cannot be directly applied to ciliates because the majority are uncultivated. Single‐cell whole‐genome sequencing (WGS) is a powerful tool for studying genetic variation in microbes, but no studies have been performed in ciliates. We compared the use of single‐cell WGS and bulk DNA WGS to detect genetic variation, specifically single nucleotide polymorphisms (SNPs), in the model ciliate Tetrahymena thermophila. Our analyses showed that (i) single‐cell WGS has excellent performance regarding mapping rate and genome coverage but lower sequencing uniformity compared with bulk DNA WGS due to amplification bias (which was reproducible); (ii) false‐positive SNP sites detected by single‐cell WGS tend to occur in genomic regions with particularly high sequencing depth and high rate of C:G to T:A base changes; (iii) SNPs detected in three or more cells should be reliable (an detection efficiency of 83.4–97.4% was obtained for combined data from three cells). This analytical method could be adapted to measure genetic variation in other ciliates and broaden research into ciliate population genetics.     

Founder events predict changes in genetic diversity during human‐mediated range expansions

Intentional or accidental introduction of species to new locations is predicted to result in loss of genetic variation and increase the likelihood of inbreeding, thus reducing population viability and evolutionary potential. However, multiple introductions and large founder numbers can prevent loss of genetic diversity and may therefore facilitate establishment success and range expansion. Based on a meta‐analysis of 119 introductions of 85 species of plants and animals, we here show a quantitative effect of founding history on genetic diversity in introduced populations. Both introduction of large number of individuals and multiple introduction events significantly contribute to maintaining or even increasing genetic diversity in introduced populations. The most consistent loss of genetic diversity is seen in insects and mammals, whereas introduced plant populations tend to have higher genetic variation than native populations. However, loss or gain of genetic diversity does not explain variation in the extent to which plant or animal populations become invasive outside of their native range. These results provide strong support for predictions from population genetics theory with respect to patterns of genetic diversity in introduced populations, but suggest that invasiveness is not limited by genetic bottlenecks.     

Exploring a Pool‐seq‐only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool‐seq data to generate a de novo genome assembly for mining exons, upon which Pool‐seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual‐based single nucleotide polymorphisms [SNPs]) and from mapping the Pool‐seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F_ST) between the two introduced populations exceeds that of the naturally sympatric populations (F_ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( ≈ 0.002 and  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high‐quality reference assembly from a divergent species. We conclude that the Pool‐seq‐only approach can be suitable for detecting and quantifying genome‐wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.     

Population genetic connectivity of an endemic New Zealand passerine after large‐scale local extirpations: a model of re‐colonization potential

Little is known about how a 70% loss of native forests has affected the genetic connectivity of remnant bird populations in New Zealand. We use the common and widely distributed New Zealand Bellbird Anthornis melanura as an indicator species of population connectivity for well‐flighted birds. Using eight microsatellite loci, we identified five main genetic populations in the North Island, South Island, sub‐Antarctic Auckland Islands and two small remnant island populations adjacent to a large region of avian extirpations in northern North Island. Only one remnant island population, on a 30‐year‐old conservation reserve at Tiritiri Matangi, displayed a clear signature of recent genetic bottleneck. The 7% migration rate at Tiritiri Matangi indicates that bottlenecks can be maintained despite habitat rehabilitation, possibly through behavioural barriers to gene flow. Adjacent to the same extirpation zone, Bellbirds on the Poor Knights Islands were found to have low genetic diversity and low re‐colonization potential. Two gaps concordant with deforestation patterns separated the Kapiti Coast of southern North Island from populations to both the north and the south. In summary, we identified linked avian habitats, as well as isolated and inbred populations and suggest that Bellbirds are good re‐colonizers. We emphasize the importance of genetic studies that assess animal dispersal among newly rehabilitated habitat patches.     

Intraspecific trait variation across scales: implications for understanding global change responses

Recognition of the importance of intraspecific variation in ecological processes has been growing, but empirical studies and models of global change have only begun to address this issue in detail. This review discusses sources and patterns of intraspecific trait variation and their consequences for understanding how ecological processes and patterns will respond to global change. We examine how current ecological models and theories incorporate intraspecific variation, review existing data sources that could help parameterize models that account for intraspecific variation in global change predictions, and discuss new data that may be needed. We provide guidelines on when it is most important to consider intraspecific variation, such as when trait variation is heritable or when nonlinear relationships are involved. We also highlight benefits and limitations of different model types and argue that many common modeling approaches such as matrix population models or global dynamic vegetation models can allow a stronger consideration of intraspecific trait variation if the necessary data are available. We recommend that existing data need to be made more accessible, though in some cases, new experiments are needed to disentangle causes of variation.     

Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

Whole‐genome sequencing approaches for conservation biology: Advantages,limitations and practical recommendations

Whole‐genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled‐individual DNA (Pool‐seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.     

Fine‐scale genetic diversity and landscape‐scale genetic structuring in three foundational bulrush species: implications for wetland revegetation

The appropriate sourcing of seeds for restoration is critical for establishing foundational plant species that support ecosystem functions and services. Genetic analyses of such species can yield insights into patterns of genetic diversity and structuring to inform seed collections. Here we document, for three foundational bulrush species, distinct genetic patterns to guide restoration of wetlands along the iconic Great Salt Lake, the largest lake in western North America. Specifically, Schoenoplectus acutus and Schoenoplectus americanus had moderate levels of site‐scale genet richness and relatively low genet richness levels within 1‐m² plots. These patterns contrast with Bolboschoenus maritimus, which had higher levels of site‐ and plot‐level genet richness, and has therefore likely experienced more recent seedling establishment. At the landscape scale, we found some evidence for genetic isolation of individuals at more remote sites (namely Fish Springs National Wildlife Refuge in the West Desert of Utah), but all species are relatively well dispersed over hundreds of kilometers, a pattern most likely to occur via avian dispersal. In our mechanistic dispersal assessment, we found abundant bulrush seeds present in waterfowl gizzards and those seeds germinated readily despite (or because of) partial digestion. Migratory waterfowl likely facilitate the broad dispersal of all species and may aid in bulrush establishment by breaking seed dormancy. These findings suggest that seeds for restoration should be collected within and among seed source sites to ensure a diverse restoration seed lot that does not disrupt gene flow patterns.     

Range limits,large‐scale biogeographic variation,and localized evolutionary dynamics in a polymorphic damselfly

Studies of heritable colour polymorphisms allow investigators to track the genetic dynamics of natural populations. By comparing polymorphic populations over large geographic areas and across generations, issues about both morph stability and evolutionary dynamics can be addressed, increasing our understanding of the potential mechanisms maintaining genetic polymorphisms. In the present study, we investigated population morph frequencies in a sex‐limited heritable colour polymorphic damselfly (Ischnura elegans, Vander Linden), with three discrete female morphs. We compared the frequencies of these three female morphs in 120 different populations from ten European countries at differing latitudes and longitudes. There were pronounced differences in morph frequencies both across the entire European biogeographic range, as well as at a smaller scale within regions. We also found considerable between‐population variation at the local scale within regions, particularly at the edges of the range of this species. We discuss these findings in the context of recent models of adaptive population divergence along the range of a species. This polymorphism is thus highly dynamic, with stable morph frequencies at the core of the species range but fluctuating morph dynamics at the range limits. We finish with a discussion of how local interactions and climatic factors can be expected to have a strong influence on the biogeographic patterns in this species and other sexually selected polymorphisms. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 775–785.     

Fine‐scale spatiotemporal patterns of genetic variation reflect budding dispersal coupled with strong natal philopatry in a cooperatively breeding mammal

The relatedness structure of animal populations is thought to be a critically important factor underlying the evolution of mating systems and social behaviours. While previous work has shown that population structure is shaped by many biological processes, few studies have investigated how these factors vary over time. Consequently, we explored the fine‐scale spatiotemporal genetic structure of an intensively studied population of cooperatively breeding banded mongooses (Mungos mungo) over a 10‐year period. Overall population structure was strong (average F_ST = 0.129) but groups with spatially overlapping territories were not more genetically similar to one another than noncontiguous groups. Instead, genetic differentiation was associated with historical group‐fission (budding) events, with new groups diverging from their parent groups over time. Within groups, relatedness was high within but not between the sexes, although the latter increased over time since group formation due to group founders being replaced by philopatric young. This trend was not mirrored by a decrease in average offspring heterozygosity over time, suggesting that close inbreeding may often be avoided, even when immigration into established groups is virtually absent and opportunities for extra‐group matings are rare. Fine‐scale spatiotemporal population structure could have important implications in social species, where relatedness between interacting individuals is a vital component in the evolution of patterns of inbreeding avoidance, reproductive skew and kin‐selected helping and harming.     

湖北海棠的等位酶变异和遗传多样性研究

采用超薄平板微型聚丙烯酰胺等电聚焦电泳方法对湖北海棠（Malus hupehensis)的9个野生居群和2个人工栽培居群的等位酶变异和遗传多样性进行了初步研究。通过对12个酶系统29个酶位点的检测,结果表明湖北海棠有25个酶位点的等位基因频率分布差异,,有10个居群发现稀有等位基因,并有11个（37．9％）重复位点;湖北海棠的遗传多样性水平很高,等位基因平均数A＝2．127,多态位点百分率P＝74．927,平均预期杂合度He=0.376;居群间的基因分化系数GST＝0．224。与其他苹果属植物相比,湖北海棠具有中等丰富的遗传变异水平。居群间的基因流仅为Nm=0.866,表明遗传漂变是影响居群遗传变异和遗传结构的一个重要因素。     

Genotyping‐by‐sequencing approaches to characterize crop genomes: choosing the right tool for the right application

In the last decade, the revolution in sequencing technologies has deeply impacted crop genotyping practice. New methods allowing rapid, high‐throughput genotyping of entire crop populations have proliferated and opened the door to wider use of molecular tools in plant breeding. These new genotyping‐by‐sequencing (GBS) methods include over a dozen reduced‐representation sequencing (RRS) approaches and at least four whole‐genome resequencing (WGR) approaches. The diversity of methods available, each often producing different types of data at different cost, can make selection of the best‐suited method seem a daunting task. We review the most common genotyping methods used today and compare their suitability for linkage mapping, genomewide association studies (GWAS), marker‐assisted and genomic selection and genome assembly and improvement in crops with various genome sizes and complexity. Furthermore, we give an outline of bioinformatics tools for analysis of genotyping data. WGR is well suited to genotyping biparental cross populations with complex, small‐ to moderate‐sized genomes and provides the lowest cost per marker data point. RRS approaches differ in their suitability for various tasks, but demonstrate similar costs per marker data point. These approaches are generally better suited for de novo applications and more cost‐effective when genotyping populations with large genomes or high heterozygosity. We expect that although RRS approaches will remain the most cost‐effective for some time, WGR will become more widespread for crop genotyping as sequencing costs continue to decrease.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.