Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Temporal and spatial dynamics of amphioxus population (Branchiostoma belcheri tsingtaneuse) and its influential factors in Luan River Estuary,China

Population distribution of amphioxus (Branchiostoma belcheri tsingtauense) was investigated in the Luan River Estuary from 1999 to 2011, to describe its trends and discover the factors that influence its decline. Amphioxus are distributed on the seabed at 5–10 m depth, between the Xinkai Estuary and the Dapu River Estuary, where the primary sand components of the sediment are medium and fine, with particle size ranging from 0.001 to 1 mm. In recent years, the population density and biomass of amphioxus have sharply decreased. Changes in sediment granularity composition may significantly influence amphioxus distribution. Our data showed that the highest density region of amphioxus occurred where >90% of the sediment particles were between 0.063 and 0.5 mm in size. A reduction in sediment discharge from Luan River and the expansion of raft‐breeding mariculture may be causing the decline in amphioxus through habitat destruction.     

Whole‐genome sequencing of two North American Drosophila melanogaster populations reveals genetic differentiation and positive selection

The prevailing demographic model for Drosophila melanogaster suggests that the colonization of North America occurred very recently from a subset of European flies that rapidly expanded across the continent. This model implies a sudden population growth and range expansion consistent with very low or no population subdivision. As flies adapt to new environments, local adaptation events may be expected. To describe demographic and selective events during North American colonization, we have generated a data set of 35 individual whole‐genome sequences from inbred lines of D. melanogaster from a west coast US population (Winters, California, USA) and compared them with a public genome data set from Raleigh (Raleigh, North Carolina, USA). We analysed nuclear and mitochondrial genomes and described levels of variation and divergence within and between these two North American D. melanogaster populations. Both populations exhibit negative values of Tajima's D across the genome, a common signature of demographic expansion. We also detected a low but significant level of genome‐wide differentiation between the two populations, as well as multiple allele surfing events, which can be the result of gene drift in local subpopulations on the edge of an expansion wave. In contrast to this genome‐wide pattern, we uncovered a 50‐kilobase segment in chromosome arm 3L that showed all the hallmarks of a soft selective sweep in both populations. A comparison of allele frequencies within this divergent region among six populations from three continents allowed us to cluster these populations in two differentiated groups, providing evidence for the action of natural selection on a global scale.     

Demographic inference from whole‐genome and RAD sequencing data suggests alternating human impacts on goose populations since the last ice age

We investigated how population changes and fluctuations in the pink‐footed goose might have been affected by climatic and anthropogenic factors. First, genomic data confirmed the existence of two separate populations: western (Iceland) and eastern (Svalbard/Denmark). Second, demographic inference suggests that the species survived the last glacial period as a single ancestral population with a low population size (100–1,000 individuals) that split into the current populations at the end of the last glacial maximum with Iceland being the most plausible glacial refuge. While population changes during the last glaciation were clearly environmental, we hypothesize that more recent demographic changes are human‐related: (1) the inferred population increase in the Neolithic is due to deforestation to establish new lands for agriculture, increasing available habitat for pink‐footed geese, (2) the decline inferred during the Middle Ages is due to human persecution, and (3) improved protection explains the increasing demographic trends during the 20th century. Our results suggest both environmental (during glacial cycles) and anthropogenic effects (more recent) can be a threat to species survival.     

Whole‐genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea

Anopheles mosquitoes are the vectors of several human diseases including malaria. In many malaria endemic areas, several species of Anopheles coexist, sometimes in the form of related sibling species that are morphologically indistinguishable. Determining the size and organization of Anopheles populations, and possible ongoing gene flow among them is important for malaria control and, in particular, for monitoring the spread of insecticide resistance alleles. However, these parameters have been difficult to evaluate in most Anopheles species due to the paucity of genetic data available. Here, we assess the extent of contemporary gene flow and historical variations in population size by sequencing and de novo assembling the genomes of wild‐caught mosquitoes from four species of the Anopheles punctulatus group of Papua New Guinea. Our analysis of more than 50 Mb of orthologous DNA sequences revealed no evidence of contemporary gene flow among these mosquitoes. In addition, investigation of the demography of two of the An. punctulatus species revealed distinct population histories. Overall, our analyses suggest that, despite their similarities in morphology, behaviour and ecology, contemporary sympatric populations of An. punctulatus are evolving independently.     

Whole‐exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high‐throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.     

A genetic demographic analysis of Lake Malawi rock‐dwelling cichlids using spatio‐temporal sampling

We estimated the effective population sizes (N_e) and tested for short‐term temporal demographic stability of populations of two Lake Malawi cichlids: Maylandia benetos, a micro‐endemic, and Maylandia zebra, a widespread species found across the lake. We sampled a total of 351 individuals, genotyped them at 13 microsatellite loci and sequenced their mitochondrial D‐loop to estimate genetic diversity, population structure, demographic history and effective population sizes. At the microsatellite loci, genetic diversity was high in all populations. Yet, genetic diversity was relatively low for the sequence data. Microsatellites yielded mean N_e estimates of 481 individuals (±99 SD) for M. benetos and between 597 (±106.3 SD) and 1524 (±483.9 SD) individuals for local populations of M. zebra. The microsatellite data indicated no deviations from mutation–drift equilibrium. Maylandia zebra was further found to be in migration–drift equilibrium. Temporal fluctuations in allele frequencies were limited across the sampling period for both species. Bayesian Skyline analyses suggested a recent expansion of M. zebra populations in line with lake‐level fluctuations, whereas the demographic history of M. benetos could only be estimated for the very recent past. Divergence time estimates placed the origin of M. benetos within the last 100 ka after the refilling of the lake and suggested that it split off the sympatric M. zebra population. Overall, our data indicate that micro‐endemics and populations in less favourable habitats have smaller N_e, indicating that drift may play an important role driving their divergence. Yet, despite small population sizes, high genetic variation can be maintained.     

Effective size of density‐dependent two‐sex populations: the effect of mating systems

Density dependence in vital rates is a key feature affecting temporal fluctuations of natural populations. This has important implications for the rate of random genetic drift. Mating systems also greatly affect effective population sizes, but knowledge of how mating system and density regulation interact to affect random genetic drift is poor. Using theoretical models and simulations, we compare N_e in short‐lived, density‐dependent animal populations with different mating systems. We study the impact of a fluctuating, density‐dependent sex ratio and consider both a stable and a fluctuating environment. We find a negative relationship between annual N_e/N and adult population size N due to density dependence, suggesting that loss of genetic variation is reduced at small densities. The magnitude of this decrease was affected by mating system and life history. A male‐biased, density‐dependent sex ratio reduces the rate of genetic drift compared to an equal, density‐independent sex ratio, but a stochastic change towards male bias reduces the N_e/N ratio. Environmental stochasticity amplifies temporal fluctuations in population size and is thus vital to consider in estimation of effective population sizes over longer time periods. Our results on the reduced loss of genetic variation at small densities, particularly in polygamous populations, indicate that density regulation may facilitate adaptive evolution at small population sizes.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

Short‐term temporal instability in fine‐scale genetic structure of masu salmon

相似文献   

Long‐term mark‐and‐recapture study of a freshwater mussel reveals patterns of habitat use and an association between survival and river discharge

相似文献   

Spatial and temporal genetic structure of Symbiodinium populations within a common reef‐building coral on the Great Barrier Reef

The dinoflagellate photosymbiont Symbiodinium plays a fundamental role in defining the physiological tolerances of coral holobionts, but little is known about the dynamics of these endosymbiotic populations on coral reefs. Sparse data indicate that Symbiodinium populations show limited spatial connectivity; however, no studies have investigated temporal dynamics for in hospite Symbiodinium populations following significant mortality and recruitment events in coral populations. We investigated the combined influences of spatial isolation and disturbance on the population dynamics of the generalist Symbiodinium type C2 (ITS1 rDNA) hosted by the scleractinian coral Acropora millepora in the central Great Barrier Reef. Using eight microsatellite markers, we genotyped Symbiodinium in a total of 401 coral colonies, which were sampled from seven sites across a 12‐year period including during flood plume–induced coral bleaching. Genetic differentiation of Symbiodinium was greatest within sites, explaining 70–86% of the total genetic variation. An additional 9–27% of variation was explained by significant differentiation of populations among sites separated by 0.4–13 km, which is consistent with low levels of dispersal via water movement and historical disturbance regimes. Sampling year accounted for 6–7% of total genetic variation and was related to significant coral mortality following severe bleaching in 1998 and a cyclone in 2006. Only 3% of the total genetic variation was related to coral bleaching status, reflecting generally small (8%) reductions in allelic diversity within bleached corals. This reduction probably reflected a loss of genotypes in hospite during bleaching, although no site‐wide changes in genetic diversity were observed. Combined, our results indicate the importance of disturbance regimes acting together with limited oceanographic transport to determine the genetic composition of Symbiodinium types within reefs.     

Demographic history and local adaptation of Myripnois dioica (Asteraceae) provide insight on plant evolution in northern China flora

The flora of northern China forms the main part of the Sino‐Japanese floristic region and is located in a south–north vegetative transect in East Asia. Phylogeographic studies have demonstrated that an arid belt in this region has promoted divergence of plants in East Asia. However, little is known about how plants that are restricted to the arid belt of flora in northern China respond to climatic oscillation and environmental change. Here, we used genomic‐level data of Myripnois dioica across its distribution as a representative of northern China flora to reconstruct plant demographic history, examine local adaptation related to environmental disequilibrium, and investigate the factors related to effective population size change. Our results indicate M. dioica originated from the northern area and expanded to the southern area, with the Taihang Mountains serving as a physical barrier promoting population divergence. Genome‐wide evidence found strong correlation between genomic variation and environmental factors, specifically signatures associated with local adaptation to drought stress in heterogeneous environments. Multiple linear regression analyses revealed joint effects of population age, mean temperature of coldest quarter, and precipitation of wettest month on effective population size (Ne). Our current study uses M. dioica as a case for providing new insights into the evolutionary history and local adaptation of northern China flora and provides qualitative strategies for plant conservation.     

Whole‐genome sequencing reveals a large deletion in the MITF gene in horses with white spotted coat colour and increased risk of deafness

White spotting phenotypes in horses are highly valued in some breeds. They are quite variable and may range from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for white spotting phenotypes in horses. For the present study, we investigated an American Paint Horse family segregating a phenotype involving white spotting and blue eyes. Six of eight horses with the white‐spotting phenotype were deaf. We obtained whole‐genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~63‐kb deletion spanning exons 6–9 of the MITF gene (chr16:21 503 211–21 566 617). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. PCR‐based genotyping revealed that all eight available affected horses from the family carried the deletion. The finding of an MITF variant fits well with the syndromic phenotype involving both depigmentation and an increased risk for deafness and corresponds to human Waardenburg syndrome type 2A. Our findings will enable more precise genetic testing for depigmentation phenotypes in horses.     

Genetic assignment of recruits reveals short‐ and long‐distance larval dispersal in Pocillopora damicornis on the Great Barrier Reef

Understanding connectivity of coral populations among and within reefs over ecologically significant timescales is essential for developing evidence‐based management strategies, including the design of marineprotected areas. Here, we present the first assessment of contemporary connectivity among populations of two Molecular Operational Taxonomic Units (MOTUs) of the brooding coral Pocillopora damicornis. We used individual‐based genetic assignment methods to identify the proportions of philopatric and migrant larval recruits, settling over 12 months at sites around Lizard Island (northern Great Barrier Reef [GBR]) and over 24 months at sites around the Palms Islands (central GBR). Overall, we found spatially and temporally variable rates of self‐recruitment and dispersal, demonstrating the importance of variation in local physical characteristics in driving dispersal processes. Recruitment patterns and inferred dispersal distances differed between the two P. damicornis MOTUs, with type α recruits exhibiting predominantly philopatric recruitment, while the majority of type β recruits were either migrants from identified putative source populations or assumed migrants based on genetic exclusion from all known populations. While P. damicornis invests much energy into brooding clonal larvae, we found that only 15% and 7% of type α and type β recruits, respectively, were clones of sampled adult colonies or other recruits, challenging the hypothesis that reproduction is predominantly asexual in this species on the GBR. We explain high rates of self‐recruitment and low rates of clonality in these MOTUs by suggesting that locally retained larvae originate predominantly from spawned gametes, while brooded larvae are mainly vagabonds.     

Effective size of density‐dependent populations in fluctuating environments

Reliable estimates of effective population size are of central importance in population genetics and evolutionary biology. For populations that fluctuate in size, harmonic mean population size is commonly used as a proxy for (multi‐) generational effective size. This assumes no effects of density dependence on the ratio between effective and actual population size, which limits its potential application. Here, we introduce density dependence on vital rates in a demographic model of variance effective size. We derive an expression for the ratio in a density‐regulated population in a fluctuating environment. We show by simulations that yearly genetic drift is accurately predicted by our model, and not proportional to as assumed by the harmonic mean model, where N is the total population size of mature individuals. We find a negative relationship between and N. For a given N, the ratio depends on variance in reproductive success and the degree of resource limitation acting on the population growth rate. Finally, our model indicate that environmental stochasticity may affect not only through fluctuations in N, but also for a given N at a given time. Our results show that estimates of effective population size must include effects of density dependence and environmental stochasticity.     

Diploidization and genome size change in allopolyploids is associated with differential dynamics of low‐ and high‐copy sequences

Recent advances have highlighted the ubiquity of whole‐genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid‐like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph‐based clustering of next‐generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low‐copy repeats. The loss of high‐copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus‐like Ty3/Gypsy retroelements and some low‐complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low‐ and high‐copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low‐copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high‐copy‐number sequences.     

Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Genetic substructure and admixture as important factors in linkage disequilibrium‐based estimation of effective number of breeders in recovering wildlife populations

The number of effective breeders (N_b) and effective population size (N_e) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium‐based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993–2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on N_b and N_e, by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for N_e were larger than N_b and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of N_b and N_e. Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on N_b and N_e.