Cryptic MHC Polymorphism Revealed but Not Explained by Selection on the Class IIB Peptide-Binding Region

The immune genes of the major histocompatibility complex (MHC) are characterized by extraordinarily high levels of nucleotide and haplotype diversity. This variation is maintained by pathogen-mediated balancing selection that is operating on the peptide-binding region (PBR). Several recent studies have found, however, that some populations possess large clusters of alleles that are translated into virtually identical proteins. Here, we address the question of how this nucleotide polymorphism is maintained with little or no functional variation for selection to operate on. We investigate circa 750–850 bp of MHC class II DAB genes in four wild populations of the guppy Poecilia reticulata. By sequencing an extended region, we uncovered 40.9% more sequences (alleles), which would have been missed if we had amplified the exon 2 alone. We found evidence of several gene conversion events that may have homogenized sequence variation. This reduces the visible copy number variation (CNV) and can result in a systematic underestimation of the CNV in studies of the MHC and perhaps other multigene families. We then focus on a single cluster, which comprises 27 (of a total of 66) sequences. These sequences are virtually identical and show no signal of selection. We use microsatellites to reconstruct the populations' demography and employ simulations to examine whether so many similar nucleotide sequences can be maintained in the populations. Simulations show that this variation does not behave neutrally. We propose that selection operates outside the PBR, for example, on linked immune genes or on the “sheltered load” that is thought to be associated to the MHC. Future studies on the MHC would benefit from extending the amplicon size to include polymorphisms outside the exon with the PBR. This may capture otherwise cryptic haplotype variation and CNV, and it may help detect other regions in the MHC that are under selection

Whole‐genome sequencing and genome regions of special interest : lessons from major histocompatibility complex, sex determination, and plant self‐incompatibility

Whole-genome sequencing of non-model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data-handling. Three such regions are the major histocompatibility complex (MHC), sex-determining regions (SDRs) and the plant self-incompatibility locus (S-locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short-read sequences against multiple alternative references improves genotyping comprehensiveness at the S-locus thereby contributing to more accurate assessments of allelic frequencies. Long-read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S-locus haplotypes. These advances are adding value to short-read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches

Ocean currents help explain population genetic structure

Management and conservation can be greatly informed by considering explicitly how environmental factors influence population genetic structure. Using simulated larval dispersal estimates based on ocean current observations, we demonstrate how explicit consideration of frequency of exchange of larvae among sites via ocean advection can fundamentally change the interpretation of empirical population genetic structuring as compared with conventional spatial genetic analyses. Both frequency of larval exchange and empirical genetic difference were uncorrelated with Euclidean distance between sites. When transformed into relative oceanographic distances and integrated into a genetic isolation-by-distance framework, however, the frequency of larval exchange explained nearly 50 per cent of the variance in empirical genetic differences among sites over scales of tens of kilometres. Explanatory power was strongest when we considered effects of multiple generations of larval dispersal via intermediary locations on the long-term probability of exchange between sites. Our results uncover meaningful spatial patterning to population genetic structuring that corresponds with ocean circulation. This study advances our ability to interpret population structure from complex genetic data characteristic of high gene flow species, validates recent advances in oceanographic approaches for assessing larval dispersal and represents a novel approach to characterize population connectivity at small spatial scales germane to conservation and fisheries management

Genetic basis and timing of a major mating system shift in Capsella

YesA crucial step in the transition from outcrossing to self-fertilization is the loss of genetic self-incompatibility (SI). In the Brassicaceae, SI involves the interaction of female and male speci-ficity components, encoded by the genesSRKandSCRat the self-incompatibility locus (S-lo-cus). Theory predicts thatS-linked mutations, and especially dominant mutations inSCR, arelikely to contribute to loss of SI. However, few studies have investigated the contribution ofdominant mutations to loss of SI in wild plant species. Here, we investigate the genetic basis of loss of SI in the self-fertilizing crucifer speciesCapsella orientalis, by combining genetic mapping, long-read sequencing of completeS-hap-lotypes, gene expression analyses and controlled crosses. We show that loss of SI inC. orientalisoccurred<2.6 Mya and maps as a dominant trait totheS-locus. We identify a fixed frameshift deletion in the male specificity geneSCRand con-firm loss of male SI specificity. We further identify anS-linked small RNA that is predicted tocause dominance of self-compatibility. Our results agree with predictions on the contribution of dominantS-linked mutations toloss of SI, and thus provide new insights into the molecular basis of mating system transitions.Work at Uppsala Genome Center is funded by 550 RFI / VR and Science for Life Laboratory, Sweden. The SNP&SEQ Platform is supported by 551 the Swedish Research Council and the Knut and Alice Wallenberg Foundation. V.C. 552 acknowledges support by a grant from the European Research Council (NOVEL project, 553 grant #648321). The authors thank the French Ministère de l’Enseignement Supérieur et de la 554 Recherche, the Hauts de France Region and the European Funds for Regional Economical 555 Development for their financial support to this project. This work was supported by a grant 556 from the Swedish Research Council (grant #D0432001) and by a grant from the Science for 557 Life Laboratory, Swedish Biodiversity Program to T.S. The Swedish Biodiversity Program is 558 supported by the Knut and Alice Wallenberg Foundation

Structural and Content Diversity of Mitochondrial Genome in Beet: A Comparative Genomic Analysis

Despite their monophyletic origin, mitochondrial (mt) genomes of plants and animals have developed contrasted evolutionary paths over time. Animal mt genomes are generally small, compact, and exhibit high mutation rates, whereas plant mt genomes exhibit low mutation rates, little compactness, larger sizes, and highly rearranged structures. We present the (nearly) whole sequences of five new mt genomes in the Beta genus: four from Beta vulgaris and one from B. macrocarpa, a sister species belonging to the same Beta section. We pooled our results with two previously sequenced genomes of B. vulgaris and studied genome diversity at the species level with an emphasis on cytoplasmic male-sterilizing (CMS) genomes. We showed that, contrary to what was previously assumed, all three CMS genomes belong to a single sterile lineage. In addition, the CMSs seem to have undergone an acceleration of the rates of substitution and rearrangement. This study suggests that male sterility emergence might have been favored by faster rates of evolution, unless CMS itself caused faster evolution

Contrasted Patterns of Molecular Evolution in Dominant and Recessive Self-Incompatibility Haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae

Fine-Scale Genetic Structure Arises during Range Expansion of an Invasive Gecko

Processes of range expansion are increasingly important in light of current concerns about invasive species and range shifts due to climate change. Theoretical studies suggest that genetic structuring may occur during range expansion. Ephemeral genetic structure can have important evolutionary implications, such as propagating genetic changes along the wave front of expansion, yet few studies have shown evidence of such structure. We tested the hypothesis that genetic structure arises during range expansion in Hemidactylus mabouia, a nocturnal African gecko recently introduced to Florida, USA. Twelve highly variable microsatellite loci were used to screen 418 individuals collected from 43 locations from four sampling sites across Florida, representing a gradient from earlier (∼1990s) to very recent colonization. We found earlier colonized locations had little detectable genetic structure and higher allelic richness than more recently colonized locations. Genetic structuring was pronounced among locations at spatial scales of tens to hundreds of meters near the leading edge of range expansion. Despite the rapid pace of range expansion in this introduced gecko, dispersal is limited among many suitable habitat patches. Fine-scale genetic structure is likely the result of founder effects during colonization of suitable habitat patches. It may be obscured over time and by scale-dependent modes of dispersal. Further studies are needed to determine if such genetic structure affects adaptation and trait evolution in range expansions and range shifts

What explains rare and conspicuous colours in a snail? A test of time-series data against models of drift, migration or selection

It is intriguing that conspicuous colour morphs of a prey species may be maintained at low frequencies alongside cryptic morphs. Negative frequency-dependent selection by predators using search images ('apostatic selection') is often suggested without rejecting alternative explanations. Using a maximum likelihood approach we fitted predictions from models of genetic drift, migration, constant selection, heterozygote advantage or negative frequency-dependent selection to time-series data of colour frequencies in isolated populations of a marine snail (Littorina saxatilis), re-established with perturbed colour morph frequencies and followed for >20 generations. Snails of conspicuous colours (white, red, banded) are naturally rare in the study area (usually <10%) but frequencies were manipulated to levels of ~50% (one colour per population) in 8 populations at the start of the experiment in 1992. In 2013, frequencies had declined to ~15-45%. Drift alone could not explain these changes. Migration could not be rejected in any population, but required rates much higher than those recorded. Directional selection was rejected in three populations in favour of balancing selection. Heterozygote advantage and negative frequency-dependent selection could not be distinguished statistically, although overall the results favoured the latter. Populations varied idiosyncratically as mild or variable colour selection (3-11%) interacted with demographic stochasticity, and the overall conclusion was that multiple mechanisms may contribute to maintaining the polymorphisms.Heredity advance online publication, 21 September 2016; doi:10.1038/hdy.2016.77

Adaptive Introgression across Species Boundaries in Heliconius Butterflies

It is widely documented that hybridisation occurs between many closely related species, but the importance of introgression in adaptive evolution remains unclear, especially in animals. Here, we have examined the role of introgressive hybridisation in transferring adaptations between mimetic Heliconius butterflies, taking advantage of the recent identification of a gene regulating red wing patterns in this genus. By sequencing regions both linked and unlinked to the red colour locus, we found a region that displays an almost perfect genotype by phenotype association across four species, H. melpomene, H. cydno, H. timareta, and H. heurippa. This particular segment is located 70 kb downstream of the red colour specification gene optix, and coalescent analysis indicates repeated introgression of adaptive alleles from H. melpomene into the H. cydno species clade. Our analytical methods complement recent genome scale data for the same region and suggest adaptive introgression has a crucial role in generating adaptive wing colour diversity in this group of butterflies

In Vivo Fluorescence Imaging of Bacteriogenic Cyanide in the Lungs of Live Mice Infected with Cystic Fibrosis Pathogens

10.1371/journal.pone.0021387PLoS ONE67