Genome-wide differentiation in closely related populations: the roles of selection and geographic isolation.

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as 'Isolation by Adaptation' (IBA). IBA can be complementary to the well-known pattern of 'Isolation by Distance' (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome-wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole genome sequence for H. rustica, to which we aligned a set of 9,493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behavior and sexual signaling, as well as geographic distance together, explain over 70% of genome-wide divergence among populations. Controlling for IBD, we find 42% of genome-wide divergence is attributable to IBA through pairwise differences in traits related to migratory behavior and sexual signaling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population-level differentiation in this group of closely related populations. This article is protected by copyright. All rights reserved

Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations

Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application

Genome-culture coevolution promotes rapid divergence of killer whale ecotypes.

Analysing population genomic data from killer whale ecotypes, which we estimate have globally radiated within less than 250,000 years, we show that genetic structuring including the segregation of potentially functional alleles is associated with socially inherited ecological niche. Reconstruction of ancestral demographic history revealed bottlenecks during founder events, likely promoting ecological divergence and genetic drift resulting in a wide range of genome-wide differentiation between pairs of allopatric and sympatric ecotypes. Functional enrichment analyses provided evidence for regional genomic divergence associated with habitat, dietary preferences and post-zygotic reproductive isolation. Our findings are consistent with expansion of small founder groups into novel niches by an initial plastic behavioural response, perpetuated by social learning imposing an altered natural selection regime. The study constitutes an important step towards an understanding of the complex interaction between demographic history, culture, ecological adaptation and evolution at the genomic level

Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map

Ajuts: Funding support is provided by Gus R. Douglass Institute (Evans Allen Project to Nimmakayala) and USDA-NIFA (2010-02247 and 2012-02511).Melon (Cucumis melo L.) is a phenotypically diverse eudicot diploid (2n = 2x = 24) has climacteric and non-climacteric morphotypes and show wide variation for fruit firmness, an important trait for transportation and shelf life. We generated 13,789 SNP markers using genotyping-by-sequencing (GBS) and anchored them to chromosomes to understand genome-wide fixation indices (Fst) between various melon morphotypes and genomewide linkage disequilibrium (LD) decay. The FST between accessions of cantalupensis and inodorus was 0.23. The FST between cantalupensis and various agrestis accessions was in a range of 0.19-0.53 and between inodorus and agrestis accessions was in a range of 0.21-0.59 indicating sporadic to wide ranging introgression. The EM (Expectation Maximization) algorithm was used for estimation of 1436 haplotypes. Average genome-wide LD decay for the melon genome was noted to be 9.27 Kb. In the current research, we focused on the genome-wide divergence underlying diverse melon horticultural groups. A high-resolution genetic map with 7153 loci was constructed. Genome-wide segregation distortion and recombination rate across various chromosomes were characterized. Melon has climacteric and non-climacteric morphotypes and wide variation for fruit firmness, a very important trait for transportation and shelf life. Various levels of QTLs were identified with high to moderate stringency and linked to fruit firmness using both genome-wide association study (GWAS) and biparental mapping. Gene annotation revealed some of the SNPs are located in β-D-xylosidase, glyoxysomal malate synthase, chloroplastic anthranilate phosphoribosyltransferase, and histidine kinase, the genes that were previously characterized for fruit ripening and softening in other crops

Contrasting signatures of genomic divergence during sympatric speciation

Population genomic analyses of Midas cichlid fishes in young Nicaraguan crater lakes suggest that sympatric speciation is promoted by polygenic architectures. The transition from 'well-marked varieties' of a single species into 'well-defined species'-especially in the absence of geographic barriers to gene flow (sympatric speciation)-has puzzled evolutionary biologists ever since Darwin(1,2). Gene flow counteracts the buildup of genome-wide differentiation, which is a hallmark of speciation and increases the likelihood of the evolution of irreversible reproductive barriers (incompatibilities) that complete the speciation process(3). Theory predicts that the genetic architecture of divergently selected traits can influence whether sympatric speciation occurs(4), but empirical tests of this theory are scant because comprehensive data are difficult to collect and synthesize across species, owing to their unique biologies and evolutionary histories(5). Here, within a young species complex of neotropical cichlid fishes (Amphilophus spp.), we analysed genomic divergence among populations and species. By generating a new genome assembly and re-sequencing 453 genomes, we uncovered the genetic architecture of traits that have been suggested to be important for divergence. Species that differ in monogenic or oligogenic traits that affect ecological performance and/or mate choice show remarkably localized genomic differentiation. By contrast, differentiation among species that have diverged in polygenic traits is genomically widespread and much higher overall, consistent with the evolution of effective and stable genome-wide barriers to gene flow. Thus, we conclude that simple trait architectures are not always as conducive to speciation with gene flow as previously suggested, whereas polygenic architectures can promote rapid and stable speciation in sympatry.Peer reviewe

History cleans up messes: The impact of time in driving divergence and introgression in a tropical suture zone

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/137755/1/evo13278-sup-0001-SuppMat.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/137755/2/evo13278.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/137755/3/evo13278_am.pd

RAD sequencing sheds new light on the genetic structure and local adaptation of European scallops and resolves their demographic histories.

Recent developments in genomics are advancing our understanding of the processes shaping population structure in wild organisms. In particular, reduced representation sequencing has facilitated the generation of dense genetic marker datasets that provide greater power for resolving population structure, investigating the role of selection and reconstructing demographic histories. We therefore used RAD sequencing to study the great scallop Pecten maximus and its sister species P. jacobeus along a latitudinal cline in Europe. Analysis of 219 samples genotyped at 82,439 single nucleotide polymorphisms clearly resolved an Atlantic and a Norwegian group within P. maximus as well as P. jacobeus, in support of previous studies. Fine-scale structure was also detected, including pronounced differences involving Mulroy Bay in Ireland, where scallops are commercially cultured. Furthermore, we identified a suite of 279 environmentally associated loci that resolved a contrasting phylogenetic pattern to the remaining neutral loci, consistent with ecologically mediated divergence. Finally, demographic inference provided support for the two P. maximus groups having diverged during the last glacial maximum and subsequently expanded, whereas P. jacobeus diverged around 95,000 generations ago and experienced less pronounced expansion. Our results provide an integrative perspective on the factors shaping genome-wide differentiation in a commercially important marine invertebrate

The Effect of Balancing Selection on Population Differentiation: A Study with HLA Genes.

Balancing selection is defined as a class of selective regimes that maintain polymorphism above what is expected under neutrality. Theory predicts that balancing selection reduces population differentiation, as measured by FST. However, balancing selection regimes in which different sets of alleles are maintained in different populations could increase population differentiation. To tackle the connection between balancing selection and population differentiation, we investigated population differentiation at the HLA genes, which constitute the most striking example of balancing selection in humans. We found that population differentiation of single nucleotide polymorphisms (SNPs) at the HLA genes is on average lower than that of SNPs in other genomic regions. We show that these results require using a computation that accounts for the dependence of FST on allele frequencies. However, in pairs of closely related populations, where genome-wide differentiation is low, differentiation at HLA is higher than in other genomic regions. Such increased population differentiation at HLA genes for recently diverged population pairs was reproduced in simulations of overdominant selection, as long as the fitness of the homozygotes differs between the diverging populations. The results give insight into a possible "divergent overdominance" mechanism for the nature of balancing selection on HLA genes across human populations

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

MM received funding from the Max Planck innovation funds for this project. PGDF was supported by a Marie Curie European Reintegration Grant (proposal nr 270891). CE was supported by German Science Foundation grants (DFG, EI 841/4-1 and EI 841/6-1)