Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Molecular Population Genetics of Rice Domestication

Domestication is a selection process that genetically modifies species to meet human needs. A most intriguing feature of domestication is the extreme phenotypic diversification among breeds. What could be the ultimate source of such genetic variations? Another notable outcome of artificial selection is the reduction in the fitness of domesticated species when they live in the wild without human assistance. The complete sequences of the two subspecies of rice cultivars provide an opportunity to address these questions. Between the two subspecies, we found much higher rates of non‐synonymous (N) than synonymous (S) substitutions and the N/S ratios are higher between cultivars than between wild species. Most interestingly, substitutions of highly dissimilar amino acids that are deleterious and uncommon between natural species are disproportionately common between the two subspecies of rice. We suggest strong selection in the absence of effective recombination may be the driving force, which we called the domestication‐associated Hill‐Robertson effect. These hitchhiking mutations may contribute to some fitness reduction in cultivars. Comparisons of the two genomes also reveal the existence of highly divergent regions in the genomes. Haplotypes in these regions often form highly polymorphic linkage blocks that are much older than speciation between wild species. Genes from such regions could contribute to the differences between indica and japonica and are likely to be involved in the diversifying selection under domestication. Their existence suggests that the amount of genetic variation within the single progenitor species Oryza rufipogon may be insufficient to account for the variation among rice cultivars, which may come from a more inclusive gene pool comprising most of the A‐genome wild species. Genes from the highly polymorphic regions also provide strong support for the independent domestication of the two subspecies. The genomic variation in rice has revealing implications for studying the genetic basis of indica‐japonica differentiation under rice domestication and subsequent improvement.     

Patterns of sequence divergence in 5' intergenic spacers and linked coding regions in 10 species of pathogenic bacteria reveal distinct recombinational histories

We compared the pattern of nucleotide difference in 8034 genes and in their 5' intergenic spacers between conspecific pairs of genomes from 10 species of pathogenic bacteria. Certain genes or spacers showed much greater sequence divergence between the genotypes compared to others; such divergent regions plausibly originated by recombinational events by which a gene and/or spacers was donated from a divergent genome. Different patterns of divergence in genes and spacers identified different recombinational patterns. For example, in Chlamydophila pneumoniae, there were examples of both unusually divergent spacers and unusually divergent genes, but there were no cases in which a gene and its spacer were both unusually divergent. This pattern suggests that, in C. pneumoniae, recombination events have broken up the linkage between genes and 5' spacers. By contrast, in Streptococcus agalactiae, there were a number of cases in which both spacer and gene were unusually divergent, indicating that a number of large-scale recombination events that included both genes and 5' spacers have occurred; there was evidence of at least two large-scale recombination events in the genomic region including the pur genes in S. agalactiae.     

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the ‘species trait'' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory.     

Sex Ratio Meiotic Drive as a Plausible Evolutionary Mechanism for Hybrid Male Sterility

Biological diversity on Earth depends on the multiplication of species or speciation, which is the evolution of reproductive isolation such as hybrid sterility between two new species. An unsolved puzzle is the exact mechanism(s) that causes two genomes to diverge from their common ancestor so that some divergent genes no longer function properly in the hybrids. Here we report genetic analyses of divergent genes controlling male fertility and sex ratio in two very young fruitfly species, Drosophila albomicans and D. nasuta. A majority of the genetic divergence for both traits is mapped to the same regions by quantitative trait loci mappings. With introgressions, six major loci are found to contribute to both traits. This genetic colocalization implicates that genes for hybrid male sterility have evolved primarily for controlling sex ratio. We propose that genetic conflicts over sex ratio may operate as a perpetual dynamo for genome divergence. This particular evolutionary mechanism may largely contribute to the rapid evolution of hybrid male sterility and the disproportionate enrichment of its underlying genes on the X chromosome – two patterns widely observed across animals.     

Patterns of genome duplication within the Brassica napus genome.

The progenitor diploid genomes (A and C) of the amphidiploid Brassica napus are extensively duplicated with 73% of genomic clones detecting two or more duplicate sequences within each of the diploid genomes. This comprehensive duplication of loci is to be expected in a species that has evolved through a polyploid ancestor. The majority of the duplicate loci within each of the diploid genomes were found in distinct linkage groups as collinear blocks of linked loci, some of which had undergone a variety of rearrangements subsequent to duplication, including inversions and translocations. A number of identical rearrangements were observed in the two diploid genomes, suggesting they had occurred before the divergence of the two species. A number of linkage groups displayed an organization consistent with centric fusion and (or) fission, suggesting this mechanism may have played a role in the evolution of Brassica genomes. For almost every genetically mapped locus detected in the A genome a homologous locus was found in the C genome; the collinear arrangement of these homologous markers allowed the primary regions of homoeology between the two genomes to be identified. At least 16 gross chromosomal rearrangements differentiated the two diploid genomes during their divergence from a common ancestor.     

Patterns of polymorphism and gene flow of gender-associated mitochondrial DNA lineages in European mussel populations

Mussels of the genus Mytilus have two types of mitochondrial DNA (mtDNA). The M type is transmitted paternally and the F type is transmitted maternally. RFLP analysis is used to assess phylogenetic relationships and nucleotide diversity and divergence for both mtDNA genomes in European populations of M. edulis and Atlantic and Mediterranean forms of M. galloprovincialis. Ten restriction endonucleases were used to assay variation in regions of the ND2 and COIII genes for a total of 77 individuals. F and M genomes show a concordant phylogenetic split into two major divergent clades, one specific to Mediterranean M. galloprovincialis and the other containing haplotypes from the three taxa. For both genomes, the geographical distribution of mtDNA variation suggests: (i) extensive levels of mtDNA introgression; (ii) asymmetric mtDNA gene flow from Atlantic to Mediterranean populations; and (iii) recurrent historical hybridization events. Significantly higher mtDNA diversity and divergence are observed for the M than F genome in all three Mytilus taxa, although the evolutionary forces responsible for these differences cannot be resolved. The extensive mtDNA gene flow between European Mytilus taxa conflicts with the restricted mtDNA introgression observed in American mussels , implying geographical variation in the nature of nuclear/mtDNA interactions regulating biparental inheritance.     

Breaking RAD: an evaluation of the utility of restriction site‐associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site‐associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome‐scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq‐based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD‐tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best‐case modelling scenario, we found that RADseq using six‐ or eight‐base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD‐tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.     

A test of genomic modularity among life‐history adaptations promoting speciation with gene flow

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or “modularity” among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple‐infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping‐by‐sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1–3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life‐history components.     

Divergence of Mammalian Higher Order Chromatin Structure Is Associated with Developmental Loci

Several recent studies have examined different aspects of mammalian higher order chromatin structure – replication timing, lamina association and Hi-C inter-locus interactions — and have suggested that most of these features of genome organisation are conserved over evolution. However, the extent of evolutionary divergence in higher order structure has not been rigorously measured across the mammalian genome, and until now little has been known about the characteristics of any divergent loci present. Here, we generate a dataset combining multiple measurements of chromatin structure and organisation over many embryonic cell types for both human and mouse that, for the first time, allows a comprehensive assessment of the extent of structural divergence between mammalian genomes. Comparison of orthologous regions confirms that all measurable facets of higher order structure are conserved between human and mouse, across the vast majority of the detectably orthologous genome. This broad similarity is observed in spite of many loci possessing cell type specific structures. However, we also identify hundreds of regions (from 100 Kb to 2.7 Mb in size) showing consistent evidence of divergence between these species, constituting at least 10% of the orthologous mammalian genome and encompassing many hundreds of human and mouse genes. These regions show unusual shifts in human GC content, are unevenly distributed across both genomes, and are enriched in human subtelomeric regions. Divergent regions are also relatively enriched for genes showing divergent expression patterns between human and mouse ES cells, implying these regions cause divergent regulation. Particular divergent loci are strikingly enriched in genes implicated in vertebrate development, suggesting important roles for structural divergence in the evolution of mammalian developmental programmes. These data suggest that, though relatively rare in the mammalian genome, divergence in higher order chromatin structure has played important roles during evolution.     

Corn and humans: recombination and linkage disequilibrium in two genomes of similar size

Two species with genomes of almost identical size, maize and human, have different evolutionary histories, and as a result their genomes differ greatly in their content of retroelements, average size of the genes and amount of genetic diversity. However, there are also significant similarities: they both have undergone bottlenecks during their recent history and seem to have non-uniform distribution of recombination events. The human genome has been shown to contain large linkage blocks characterized by a limited number of haplotypes. A similar linkage block structure is likely to exist in maize. Although highly diverse maize populations show rapid decline of linkage disequilibrium, as in humans, it is possible to define populations with strong linkage disequilibrium, suitable for whole-genome scan association mapping. The genetic diversity and lack of sequence homology found in maize influences recombinational properties and local linkage disequilibrium levels but also challenges our understanding of the relationship between the genome sequence and species definition.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Can AFLP genome scans detect small islands of differentiation? The case of shell sculpture variation in the periwinkle Echinolittorina hawaiiensis

Genome scans have identified candidate regions of the genome undergoing selection in a wide variety of organisms, yet have rarely been applied to broadly dispersing marine organisms experiencing divergent selection pressures, where high recombination rates can reduce the extent of linkage disequilibrium (LD) and the ability to detect genomic regions under selection. The broadly dispersing periwinkle Echinolittorina hawaiiensis exhibits a heritable shell sculpture polymorphism that is correlated with environmental variation. To elucidate the genetic basis of phenotypic variation, a genome scan using over 1000 AFLP loci was conducted on smooth and sculptured snails from divergent habitats at four replicate sites. Approximately 5% of loci were identified as outliers with Dfdist, whereas no outliers were identified by BayeScan. Closer examination of the Dfdist outliers supported the conclusion that these loci were false positives. These results highlight the importance of controlling for Type I error using multiple outlier detection approaches, multitest corrections and replicate population comparisons. Assuming shell phenotypes have a genetic basis, our failure to detect outliers suggests that the life history of the target species needs to be considered when designing a genome scan.     

Linkage disequilibria between pairs of loci within a highly polymorphic region of chromosome 2Q.

A recently described region on chromosome 2q contains seven restriction fragment length polymorphisms (RFLPs) revealed by single-copy probes isolated from a 20-kilobase (kb) segment of a single cosmid insert. Analysis of six of these loci demonstrates modest amounts of linkage disequilibrium. This reflects the presence of a substantial number of different haplotypes in this chromosome region and indicates that the region could be used as one highly polymorphic locus. No consistent relationship is found between the amount of linkage disequilibrium and the physical distance between pairs of loci. For seven of the 10 pairs of diallelic loci studied, the observed disequilibrium can be attributed primarily to the absence of the minor haplotype from the population. These results suggest that, for small regions of the genome, factors such as mutation, genetic drift, and population admixture may have effects that outweight those of recombination. In addition, results are reviewed which show that estimates of linkage disequilibrium coefficients for tightly linked loci are very imprecise. Thus, the inference of gene order from linkage disequilibrium values must be regarded with caution.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

Evolutionary dynamics of the Pgk1 gene in the polyploid genus Kengyilia (Triticeae: Poaceae) and its diploid relatives

The level and pattern of nucleotide variation in duplicate gene provide important information on the evolutionary history of polyploids and divergent process between homoeologous loci within lineages. Kengyilia is a group of allohexaploid species with the StYP genomic constitutions in the wheat tribe. To investigate the evolutionary dynamics of the Pgk1 gene in Kengyilia and its diploid relatives, three copies of Pgk1 homoeologues were isolated from all sampled hexaploid Kengyilia species and analyzed with the Pgk1 sequences from 47 diploid taxa representing 18 basic genomes in Triticeae. Sequence diversity patterns and genealogical analysis suggested that (1) Kengyilia species from the Central Asia and the Qinghai-Tibetan plateau have independent origins with geographically differentiated P genome donors and diverged levels of nucleotide diversity at Pgk1 locus; (2) a relatively long-time sweep event has allowed the Pgk1 gene within Agropyron to adapt to cold climate triggered by the recent uplifts of the Qinghai-Tibetan Plateau; (3) sweep event and population expansion might result in the difference in the d(N)/d(S) value of the Pgk1 gene in allopatric Agropyron populations, and this difference may be genetically transmitted to Kengyilia lineages via independent polyploidization events; (4) an 83 bp MITE element insertion has shaped the Pgk1 loci in the P genome lineage with different geographical regions; (5) the St and P genomes in Kengyilia were donated by Pseudoroegneria and Agropyron, respectively, and the Y genome is closely related to the Xp genome of Peridictyon sanctum. The interplay of evolutionary forces involving diverged natural selection, population expansion, and transposable events in geographically differentiated P genome donors could attribute to geographical differentiation of Kengyilia species via independent origins.     

基于大豆胞囊线虫病抗性候选基因的SNP位点遗传变异分析

以我国363份栽培和野生大豆资源为材料, 对大豆胞囊线虫抗性候选基因(rhg1和Rhg4)的SNP位点(8个)进行遗传变异分析, 以期阐明野生和栽培大豆间遗传多样性及连锁不平衡水平差异。结果表明, 与野生大豆相比, 代表我国栽培大豆总体资源多样性的微核心种质及其补充材料的连锁不平衡水平较高(R²值为0.216)。在栽培大豆群体内, 基因内和基因间分别有100％和16.6％的SNP位点对连锁不平衡显著, 形成两个基因特异的连锁不平衡区间(Block)。在所有供试材料中共检测到单倍型46个, 野生大豆的单倍型数目(27)少于栽培大豆(31), 但单倍型多样性(0.916)稍高于栽培大豆(0.816)。单倍型大多数(67.4％)为群体所特有(31个), 其中15个为野生大豆特有单倍型。野生大豆的两个主要优势单倍型(Hap_10和Hap_11)在栽培大豆中的发生频率也明显下降, 推测野生大豆向栽培大豆进化过程中, 一方面形成了新的单倍型, 另一方面因为瓶颈效应部分单倍型的频率降低甚至消失。