skip to main content


Title: Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome
Abstract

A whole‐genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid‐origin lineages. However, little is known about general consequences of aWGDbecause gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid‐origin species. We demonstrate a new method that enables genome‐wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus‐specific copy number into account. We apply this method toRADsequence data from different ecotypes of a polyploid‐origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences ofWGDs and local segmental gene duplications.

 
more » « less
NSF-PAR ID:
10039732
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
26
Issue:
17
ISSN:
0962-1083
Format(s):
Medium: X Size: p. 4509-4522
Size(s):
["p. 4509-4522"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identifySNPmarkers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein‐coding and nearby 5′ and 3′ untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis) exon capture data and directly from the domestic sheep genome (Ovis ariesv. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR‐basedSNPchip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositanandbayescan), we detected 28SNPloci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease‐regulating functions (e.g. Ovar‐DRA,APC,BATF2,MAGEB18), cell regulation signalling pathways (e.g.KRIT1,PI3K,ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene‐targetedSNPdiscovery and subsequentSNPchip genotyping using low‐quality samples in a nonmodel species.

     
    more » « less
  2. Abstract

    Identifying genomic signatures of natural selection can be challenging against a background of demographic changes such as bottlenecks and population expansions. Here, we disentangle the effects of demography from selection in the House Finch (Haemorhous mexicanus) using samples collected before and after a pathogen‐induced selection event. Using ddRADseq, we genotyped over 18,000SNPs across the genome in native pre‐epizootic westernUSbirds, introduced birds from Hawaii and the eastern United States, post‐epizootic eastern birds, and western birds sampled across a similar time span. We found 14% and 7% reductions in nucleotide diversity, respectively, in Hawaiian and pre‐epizootic eastern birds relative to pre‐epizootic western birds, as well as elevated levels of linkage disequilibrium and other signatures of founder events. Despite finding numerous significant frequency shifts (outlier loci) between pre‐epizootic native and introduced populations, we found no signal of reduced genetic diversity, elevated linkage disequilibrium, or outlier loci as a result of the epizootic. Simulations demonstrate that the proportion of outliers associated with founder events could be explained by genetic drift. This rare view of genetic evolution across time in an invasive species provides direct evidence that demographic shifts like founder events have genetic consequences more widespread across the genome than natural selection.

     
    more » « less
  3. Premise

    Male gametophytes of most seed plants deliver sperm to eggs via a pollen tube. Pollen tube growth rates (PTGRs) of angiosperms are exceptionally rapid, a pattern attributed to more effective haploid selection under stronger pollen competition. Paradoxically, whole genome duplication (WGD) has been common in angiosperms but rare in gymnosperms. Pollen tube polyploidy should initially acceleratePTGRbecause increased heterozygosity and gene dosage should increase metabolic rates. However, polyploidy should also independently increase tube cell size, causing more work which should decelerate growth. We asked how genome size changes have affected the evolution of seed plantPTGRs.

    Methods

    We assembled a phylogenetic tree of 451 species with knownPTGRs. We then used comparative phylogenetic methods to detect effects of neo‐polyploidy (within‐genus origins),DNAcontent, andWGDhistory onPTGR, and correlated evolution ofPTGRandDNAcontent.

    Results

    Gymnosperms had significantly higherDNAcontent and slowerPTGRoptima than angiosperms, and theirPTGRandDNAcontent were negatively correlated. For angiosperms, 89% of model weight favored Ornstein‐Uhlenbeck models with a fasterPTGRoptimum for neo‐polyploids, whereasPTGRandDNAcontent were not correlated. For within‐genus and intraspecific‐cytotype pairs,PTGRs of neo‐polyploids < paleo‐polyploids.

    Conclusions

    Genome size increases should negatively affectPTGRwhen genetic consequences ofWGDs are minimized, as found in intra‐specific autopolyploids (low heterosis) and gymnosperms (fewWGDs). But in angiosperms, the higherPTGRoptimum of neo‐polyploids and non‐negativePTGRDNAcontent correlation suggest that recurrentWGDs have caused substantialPTGRevolution in a non‐haploid state.

     
    more » « less
  4. Abstract

    A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome‐wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome‐wide patterns of divergence. Aided by a first‐draft genome assembly, we combined anFCTscan and geographic cline analysis to characterize patterns of genome‐wide divergence between divergent pollination ecotypes ofMimulus aurantiacus.FCTanalysis of 58 872SNPs generated viaRAD‐seq revealed little ecotypic differentiation (meanFCT = 0.041), although a small number of loci were moderately‐to‐highly diverged. Consistent with our previous results from the geneMaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under‐utilized framework for studying genome‐wide patterns of divergence.

     
    more » « less
  5. Abstract

    Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria,Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i ‘amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations whereP. relictumexists, and can sustain infection without major fitness consequences. High‐elevation, unexposed populations of ‘amakihi display little to no tolerance. To explore the genomic basis of adaptation toP. relictumin low‐elevation ‘amakihi, we genotyped 125 ‘amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containingSNPs and used the reference ‘amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low‐ and high‐elevation population pairs and identified loci with signatures of selection within low‐elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta‐defensin, glycoproteins and interleukin‐related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

     
    more » « less