Recombination and speciation: loci near centromeres are more differentiated than loci near telomeres between subspecies of the European rabbit (Oryctolagus cuniculus)

Conformation-associated health in pet rabbits in the UK: A VetCompass cohort study

BACKGROUND

Domestic rabbit breeds vary substantially from the wild rabbit body type. However, little is known about how the conformation of pet rabbits influences their health.

METHODS

Data were extracted from VetCompass anonymised clinical records of rabbits under UK primary veterinary care during 2019.

RESULTS

The study included 162,107 rabbits. Based on 88,693 rabbits with relevant breed information recorded, skull shape was classified as brachycephalic (79.69%), mesaticephalic (16.80%) and dolichocephalic (3.51%). Based on 83,821 rabbits with relevant breed information recorded, ear carriage was classified as lop-eared (57.05%) and erect-eared (42.95%). From a random sample of 3933 rabbits, the most prevalent disorders recorded overall were overgrown nail(s) (28.19%), overgrown molar(s) (14.90%) and obesity (8.82%). Compared to those with a mesaticephalic skull shape, brachycephalic rabbits had lower odds of obesity, anorexia and gastrointestinal stasis and higher odds of perineal faecal impaction, tear duct abnormality and haircoat disorder. Compared to erect-eared rabbits, lop-eared rabbits had higher odds of perineal faecal impaction and tear duct abnormality.

LIMITATION

A large proportion of records with incomplete breed information hindered full analysis for breed-related and conformation-related attributes.

CONCLUSION

Limited evidence for major links between skull shape or ear carriage conformations and overall disorder risk suggests that factors such as husbandry or even just living life as a domesticated species may be bigger drivers of common health issues in pet rabbits in the UK.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

The roles of recombination and selection in shaping genomic divergence in an incipient ecological species complex

Speciation genomic studies have revealed that genomes of diverging lineages are shaped jointly by the actions of gene flow and selection. These evolutionary forces acting in concert with processes such as recombination and genome features such as gene density shape a mosaic landscape of divergence. We investigated the roles of recombination and gene density in shaping the patterns of differentiation and divergence between the cyclically parthenogenetic ecological sister-taxa, Daphnia pulicaria and Daphnia pulex. First, we assembled a phased chromosome-scale genome assembly using trio-binning for D. pulicaria and constructed a genetic map using an F2-intercross panel to understand sex-specific recombination rate heterogeneity. Finally, we used a ddRADseq data set with broad geographic sampling of D. pulicaria, D. pulex, and their hybrids to understand the patterns of genome-scale divergence and demographic parameters. Our study provides the first sex-specific estimates of recombination rates for a cyclical parthenogen, and unlike other eukaryotic species, we observed male-biased heterochiasmy in D. pulicaria, which may be related to this somewhat unique breeding mode. Additionally, regions of high gene density and recombination are generally more divergent than regions of suppressed recombination. Outlier analysis indicated that divergent genomic regions are probably driven by selection on D. pulicaria, the derived lineage colonizing a novel lake habitat. Together, our study supports a scenario of selection acting on genes related to local adaptation shaping genome-wide patterns of differentiation despite high local recombination rates in this species complex. Finally, we discuss the limitations of our data in light of demographic uncertainty.

Centromere drive: chromatin conflict in meiosis

Centromeres are essential loci in eukaryotes that are necessary for the faithful segregation of chromosomes in mitosis and meiosis. Centromeres organize the kinetochore, the protein machine that attaches sister chromatids or homologous chromosomes to spindle microtubules and regulates their disjunction. Centromeres have both genetic and epigenetic determinants, which can come into conflict in asymmetric female meiosis in seed plants and animals. The centromere drive model was proposed to describe this conflict and explain how it leads to the rapid evolution of both centromeres and kinetochores. Recent studies confirm key aspects of the centromere drive model, clarify its mechanisms, and implicate rapid centromere/kinetochore evolution in hybrid inviability between species.

Satellitome analysis illuminates the evolution of ZW sex chromosomes of Triportheidae fishes (Teleostei: Characiformes)

Satellites are an abundant source of repetitive DNAs that play an essential role in the chromosomal organization and are tightly linked with the evolution of sex chromosomes. Among fishes, Triportheidae stands out as the only family where almost all species have a homeologous ZZ/ZW sex chromosomes system. While the Z chromosome is typically conserved, the W is always smaller, with variations in size and morphology between species. Here, we report an analysis of the satellitome of Triportheus auritus (TauSat) by integrating genomic and chromosomal data, with a special focus on the highly abundant and female-biased satDNAs. In addition, we investigated the evolutionary trajectories of the ZW sex chromosomes in the Triportheidae family by mapping satDNAs in selected representative species of this family. The satellitome of T. auritus comprised 53 satDNA families of which 24 were also hybridized by FISH. Most satDNAs differed significantly between sexes, with 19 out of 24 being enriched on the W chromosome of T. auritus. The number of satDNAs hybridized into the W chromosomes of T. signatus and T. albus decreased to six and four, respectively, in accordance with the size of their W chromosomes. No TauSat probes produced FISH signals on the chromosomes of Agoniates halecinus. Despite its apparent conservation, our results indicate that each species differs in the satDNA accumulation on the Z chromosome. Minimum spanning trees (MSTs), generated for three satDNA families with different patterns of FISH mapping data, revealed different homogenization rates between the Z and W chromosomes. These results were linked to different levels of recombination between them. The most abundant satDNA family (TauSat01) was exclusively hybridized in the centromeres of all 52 chromosomes of T. auritus, and its putative role in the centromere evolution was also highlighted. Our results identified a high differentiation of both ZW chromosomes regarding satellites composition, highlighting their dynamic role in the sex chromosomes evolution.

Bidirectional introgression between Mus musculus domesticus and Mus spretus

Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles—including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Keeping an eye on the use of eye-lens weight as a universal indicator of age for European wild rabbits

Accurate methods for age determination are critical to the knowledge of wildlife populations' age structure and, therefore, to their successful management. The reliability of age estimation may have profound economic and ecological consequences on the management of the European wild rabbits, Oryctolagus cuniculus, in its native and introduced range, where it is a keystone species and a major pest, respectively. As in other mammal species, European rabbits' age is often estimated using the Gompertz relationship between age and lens' weight. The growth rate formula has been developed based on data collected from European rabbits introduced in Australia, where a single subspecies (O. cuniculus cuniculus, Occ) is present. However, this curve has never been validated in the species native range, the Iberian Peninsula, where two subspecies (Occ, and O. c. algirus, Oca) coexist naturally. In this study, we tested the relationship between age and lens' weight using 173 Occ and 112 Oca wild rabbits that were surveyed in two experimental facilities in Spain. Our findings show that, in the native range, the published growth curve formula fits well Occ but not Oca data. Therefore, we recommend using the formula reported in this study to estimate the age of Oca (Lens dry weight = 240 × 10^{(-64.9/(Age+32))}). This study supports Oca rabbits' distinctiveness revealed by previous studies, which suggests that management interventions should be applied to protect this subspecies whose distribution range is very narrow and whose populations seem to be declining. More broadly, our findings point to the importance of testing the suitability of growth curves defined for other species with different genetic forms as occurs in the European wild rabbit case.

The Genomic Landscape of Divergence Across the Speciation Continuum in Island-Colonising Silvereyes (Zosterops lateralis)

Inferring the evolutionary dynamics at play during the process of speciation by analyzing the genomic landscape of divergence is a major pursuit in population genomics. However, empirical assessments of genomic landscapes under varying evolutionary scenarios that are known a priori are few, thereby limiting our ability to achieve this goal. Here we combine RAD-sequencing and individual-based simulations to evaluate the genomic landscape of divergence in the silvereye (Zosterops lateralis). Using pairwise comparisons that differ in divergence timeframe and the presence or absence of gene flow, we document how genomic patterns accumulate along the speciation continuum. In contrast to previous predictions, our results provide limited support for the idea that divergence accumulates around loci under divergent selection or that genomic islands widen with time. While a small number of genomic islands were found in populations diverging with and without gene flow, in few cases were SNPs putatively under selection tightly associated with genomic islands. The transition from localized to genome-wide levels of divergence was captured using individual-based simulations that considered only neutral processes. Our results challenge the ubiquity of existing verbal models that explain the accumulation of genomic differences across the speciation continuum and instead support the idea that divergence both within and outside of genomic islands is important during the speciation process.

Staphylococci among Wild European Rabbits from the Azores: A Potential Zoonotic Issue?

ABSTRACT

The prevalence and diversity of Staphylococcus species from wild European rabbits (Oryctolagus cuniculus) in the Azores were investigated, and the antibiotic resistance phenotype and genotype of the isolates were determined. Nasal samples from 77 wild European rabbits from São Jorge and São Miguel islands in Azores were examined. Antibiotic susceptibility of the isolates was determined using the Kirby-Bauer disk diffusion method, and the presence of antimicrobial resistance genes and virulence factors was determined by PCR. The genetic lineages of S. aureus isolates were characterized by spa typing and multilocus sequence typing. A total of 49 staphylococci were obtained from 35 of the 77 wild rabbits. Both coagulase-positive (8.2%) and coagulase-negative (91.8%) staphylococci were detected: 4 S. aureus, 17 S. fleurettii, 13 S. sciuri, 7 S. xylosus, 4 S. epidermidis, and 1 each of S. simulans, S. saprophyticus, S. succinus, and S. equorum. The four S. aureus isolates showed methicillin susceptibility and were characterized as spa type t272/CC121, Panton-Valentine leukocidin negative, and hlB positive. Most of the coagulase-negative staphylococci showed resistance to fusidic acid and beta-lactams, and multidrug resistance was identified especially among S. epidermidis isolates. The mecA gene was detected in 20 isolates of the species S. fleurettii and S. epidermidis, associated with the blaZ gene in one S. epidermidis isolate. Five antimicrobial resistance genes were detected in one S. epidermidis isolate (mecA,dfrA,dfrG,aac6'-aph2'', and ant4). Our results highlight that wild rabbits are reservoirs or "temporary hosts" of Staphylococcus species with zoonotic potential, some of them carrying relevant antimicrobial resistances.

HIGHLIGHTS

Evolutionary history of two cryptic species of northern African jerboas

BACKGROUND

Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence.

RESULTS

The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization.

CONCLUSIONS

Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.

What makes a centromere?

Centromeres are the eukaryotic chromosomal sites at which the kinetochore forms and attaches to spindle microtubules to orchestrate chromosomal segregation in mitosis and meiosis. Although centromeres are essential for cell division, their sequences are not conserved and evolve rapidly. Centromeres vary dramatically in size and organization. Here we categorize their diversity and explore the evolutionary forces shaping them. Nearly all centromeres favor AT-rich DNA that is gene-free and transcribed at a very low level. Repair of frequent centromere-proximal breaks probably contributes to their rapid sequence evolution. Point centromeres are only ~125 bp and are specified by common protein-binding motifs, whereas short regional centromeres are 1-5 kb, typically have unique sequences, and may have pericentromeric repeats adapted to facilitate centromere clustering. Transposon-rich centromeres are often ~100-300 kb and are favored by RNAi machinery that silences transposons, by suppression of meiotic crossovers at centromeres, and by the ability of some transposons to target centromeres. Megabase-length satellite centromeres arise in plants and animals with asymmetric female meiosis that creates centromere competition, and favors satellite monomers one or two nucleosomes in length that position and stabilize centromeric nucleosomes. Holocentromeres encompass the length of a chromosome and may differ dramatically between mitosis and meiosis. We propose a model in which low level transcription of centromeres facilitates the formation of non-B DNA that specifies centromeres and promotes loading of centromeric nucleosomes.

Genomic Prediction Including SNP-Specific Variance Predictors

The increasing amount of available biological information on the markers can be used to inform the models applied for genomic selection to improve predictions. The objective of this study was to propose a general model for genomic selection using a link function approach within the hierarchical generalized linear model framework (hglm) that can include external information on the markers. These models can be fitted using the well-established hglm package in R. We also present an R package (CodataGS) to fit these models, which is significantly faster than the hglm package. Simulated data were used to validate the proposed model. We tested categorical, continuous and combination models where the external information on the markers was related to 1) the location of the QTL on the genome with varying degree of uncertainty, 2) the relationship of the markers with the QTL calculated as the LD between them, and 3) a combination of both. The proposed models showed improved accuracies from 3.8% up to 23.2% compared to the SNP-BLUP method in a simulated population derived from a base population with 100 individuals. Moreover, the proposed categorical model was tested on a dairy cattle dataset for two traits (Milk Yield and Fat Percentage). These results also showed improved accuracy compared to SNP-BLUP, especially for the Fat% trait. The performance of the proposed models depended on the genetic architecture of the trait, as traits that deviate from the infinitesimal model benefited more from the external information. Also, the gain in accuracy depended on the degree of uncertainty of the external information provided to the model. The usefulness of these type of models is expected to increase with time as more accurate information on the markers becomes available.

Selection, Linkage, and Population Structure Interact To Shape Genetic Variation Among Threespine Stickleback Genomes

The outcome of selection on genetic variation depends on the geographic organization of individuals and populations as well as the organization of loci within the genome. Spatially variable selection between marine and freshwater habitats has had a significant and heterogeneous impact on patterns of genetic variation across the genome of threespine stickleback fish. When marine stickleback invade freshwater habitats, more than a quarter of the genome can respond to divergent selection, even in as little as 50 years. This process largely uses standing genetic variation that can be found ubiquitously at low frequency in marine populations, can be millions of years old, and is likely maintained by significant bidirectional gene flow. Here, we combine population genomic data of marine and freshwater stickleback from Cook Inlet, Alaska, with genetic maps of stickleback fish derived from those same populations to examine how linkage to loci under selection affects genetic variation across the stickleback genome. Divergent selection has had opposing effects on linked genetic variation on chromosomes from marine and freshwater stickleback populations: near loci under selection, marine chromosomes are depauperate of variation, while these same regions among freshwater genomes are the most genetically diverse. Forward genetic simulations recapitulate this pattern when different selective environments also differ in population structure. Lastly, dense genetic maps demonstrate that the interaction between selection and population structure may impact large stretches of the stickleback genome. These findings advance our understanding of how the structuring of populations across geography influences the outcomes of selection, and how the recombination landscape broadens the genomic reach of selection.

Parallel adaptation of rabbit populations to myxoma virus

In the 1950s the myxoma virus was released into European rabbit populations in Australia and Europe, decimating populations and resulting in the rapid evolution of resistance. We investigated the genetic basis of resistance by comparing the exomes of rabbits collected before and after the pandemic. We found a strong pattern of parallel evolution, with selection on standing genetic variation favoring the same alleles in Australia, France, and the United Kingdom. Many of these changes occurred in immunity-related genes, supporting a polygenic basis of resistance. We experimentally validated the role of several genes in viral replication and showed that selection acting on an interferon protein has increased the protein's antiviral effect.

Incomplete lineage sorting impacts the inference of macroevolutionary regimes from molecular phylogenies when concatenation is employed: An analysis based on Cetacea

Interest in methods that estimate speciation and extinction rates from molecular phylogenies has increased over the last decade. The application of such methods requires reliable estimates of tree topology and node ages, which are frequently obtained using standard phylogenetic inference combining concatenated loci and molecular dating. However, this practice disregards population-level processes that generate gene tree/species tree discordance. We evaluated the impact of employing concatenation and coalescent-based phylogeny inference in recovering the correct macroevolutionary regime using simulated data based on the well-established diversification rate shift of delphinids in Cetacea. We found that under scenarios of strong incomplete lineage sorting, macroevolutionary analysis of phylogenies inferred by concatenating loci failed to recover the delphinid diversification shift, while the coalescent-based tree consistently retrieved the correct rate regime. We suggest that ignoring microevolutionary processes reduces the power of methods that estimate macroevolutionary regimes from molecular data.

Natural selection interacts with recombination to shape the evolution of hybrid genomes

To investigate the consequences of hybridization between species, we studied three replicate hybrid populations that formed naturally between two swordtail fish species, estimating their fine-scale genetic map and inferring ancestry along the genomes of 690 individuals. In all three populations, ancestry from the "minor" parental species is more common in regions of high recombination and where there is linkage to fewer putative targets of selection. The same patterns are apparent in a reanalysis of human and archaic admixture. These results support models in which ancestry from the minor parental species is more likely to persist when rapidly uncoupled from alleles that are deleterious in hybrids. Our analyses further indicate that selection on swordtail hybrids stems predominantly from deleterious combinations of epistatically interacting alleles.

Genomic islands of differentiation in two songbird species reveal candidate genes for hybrid female sterility

Hybrid sterility is a common first step in the evolution of postzygotic reproductive isolation. According to Haldane's Rule, it affects predominantly the heterogametic sex. While the genetic basis of hybrid male sterility in organisms with heterogametic males has been studied for decades, the genetic basis of hybrid female sterility in organisms with heterogametic females has received much less attention. We investigated the genetic basis of reproductive isolation in two closely related avian species, the common nightingale (Luscinia megarhynchos) and the thrush nightingale (L. luscinia), that hybridize in a secondary contact zone and produce viable hybrid progeny. In accordance with Haldane's Rule, hybrid females are sterile, while hybrid males are fertile, allowing gene flow to occur between the species. Using transcriptomic data from multiple individuals of both nightingale species, we identified genomic islands of high differentiation (F_ST ) and of high divergence (D_xy ), and we analysed gene content and patterns of molecular evolution within these islands. Interestingly, we found that these islands were enriched for genes related to female meiosis and metabolism. The islands of high differentiation and divergence were also characterized by higher levels of linkage disequilibrium than the rest of the genome in both species indicating that they might be situated in genomic regions of low recombination. This study provides one of the first insights into genetic basis of hybrid female sterility in organisms with heterogametic females.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Elevated Proportions of Deleterious Genetic Variation in Domestic Animals and Plants

A fraction of genetic variants segregating in any population are deleterious, which negatively impacts individual fitness. The domestication of animals and plants is associated with population bottlenecks and artificial selection, which are predicted to increase the proportion of deleterious variants. However, the extent to which this is a general feature of domestic species is unclear. Here, we examine the effects of domestication on the prevalence of deleterious variation using pooled whole-genome resequencing data from five domestic animal species (dog, pig, rabbit, chicken, and silkworm) and two domestic plant species (rice and soybean) compared with their wild ancestors. We find significantly reduced genetic variation and increased proportion of nonsynonymous amino acid changes in all but one of the domestic species. These differences are observable across a range of allele frequencies, both common and rare. We find proportionally more single nucleotide polymorphisms in highly conserved elements in domestic species and a tendency for domestic species to harbor a higher proportion of changes classified as damaging. Our findings most likely reflect an increased incidence of deleterious variants in domestic species, which is most likely attributable to population bottlenecks that lead to a reduction in the efficacy of selection. An exception to this pattern is displayed by European domestic pigs, which do not show traces of a strong population bottleneck and probably continued to exchange genes with wild boar populations after domestication. The results presented here indicate that an elevated proportion of deleterious variants is a common, but not ubiquitous, feature of domestic species.

The evolution of recombination rates in finite populations during ecological speciation

Recombination can impede ecological speciation with gene flow by mixing locally adapted genotypes with maladapted migrant genotypes from a divergent population. In such a scenario, suppression of recombination can be selectively favoured. However, in finite populations evolving under the influence of random genetic drift, recombination can also facilitate adaptation by reducing Hill-Robertson interference between loci under selection. In this case, increased recombination rates can be favoured. Although these two major effects on recombination have been studied individually, their joint effect on ecological speciation with gene flow remains unexplored. Using a mathematical model, we investigated the evolution of recombination rates in two finite populations that exchange migrants while adapting to contrasting environments. Our results indicate a two-step dynamic where increased recombination is first favoured (in response to the Hill-Robertson effect), and then disfavoured, as the cost of recombining locally with maladapted migrant genotypes increases over time (the maladaptive gene flow effect). In larger populations, a stronger initial benefit for recombination was observed, whereas high migration rates intensify the long-term cost of recombination. These dynamics may have important implications for our understanding of the conditions that facilitate incipient speciation with gene flow and the evolution of recombination in finite populations.

Hybridization patterns in two contact zones of grass snakes reveal a new Central European snake species

Recent studies found major conflicts between traditional taxonomy and genetic differentiation of grass snakes and identified previously unknown secondary contact zones. Until now, little is known about gene flow across these contact zones. Using two mitochondrial markers and 13 microsatellite loci, we examined two contact zones. One, largely corresponding to the Rhine region, involves the western subspecies Natrix natrix helvetica and the eastern subspecies N. n. natrix, whereas in the other, more easterly, contact zone two lineages meet that are currently identified with N. n. natrix and N. n. persa. This second contact zone runs across Central Europe to the southern Balkans. Our analyses reveal that the western contact zone is narrow, with parapatrically distributed mitochondrial lineages and limited, largely unidirectional nuclear gene flow. In contrast, the eastern contact zone is very wide, with massive nuclear admixture and broadly overlapping mitochondrial lineages. In combination with additional lines of evidence (morphology, phylogeny, divergence times), we conclude that these differences reflect different stages in the speciation process and that Natrix helvetica should be regarded as a distinct species. We suggest a nomenclatural framework for presently recognized grass snake taxa and highlight the need for reconciling the conflicts between genetics and taxonomy.

Genomewide patterns of variation in genetic diversity are shared among populations, species and higher-order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST , PBS, d_xy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

Mapping centromeres of microchromosomes in the zebra finch (Taeniopygia guttata) using half-tetrad analysis

Centromeres usually consist of hundreds of kilobases of repetitive sequence which renders them difficult to assemble. As a consequence, centromeres are often missing from assembled genomes and their locations on physical chromosome maps have to be inferred from flanking sequences via fluorescence in situ hybridization (FISH). Alternatively, centromere positions can be mapped using linkage analyses in accidentally triploid individuals formed by half-tetrads (resulting from the inheritance of two chromatids from a single meiosis). The current genome assembly of the zebra finch (Taeniopygia guttata) comprises 32 chromosomes, but only for the ten largest chromosomes centromere positions have been mapped using FISH. We here map the positions of most of the remaining centromeres using half-tetrad analyses. For this purpose, we genotyped 37 zebra finches that were triploid or tetraploid due to inheritance errors (and mostly died as embryos) together with their parents at 64 microsatellite markers (at least two per chromosome). Using the information on centromere positions on the ten largest chromosomes, we were able to identify 12 cases of non-disjunction in maternal meiosis I and 10 cases of non-disjunction in maternal meiosis II. These 22 informative cases allowed us to infer centromere positions on additional 19 microchromosomes in reference to the current genome assembly. This knowledge will be valuable for studies of chromosome evolution, meiotic drive and species divergence in the avian lineage.

Linkage mapping of a polymorphic plumage locus associated with intermorph incompatibility in the Gouldian finch (Erythrura gouldiae)

Colour polymorphism is known to facilitate speciation but the genetic basis of animal pigmentation and how colour polymorphisms contribute to speciation is poorly understood. Restricted recombination may promote linkage disequilibrium between the colour locus and incompatibility genes. Genomic rearrangement and the position of relevant loci within a chromosome are important factors that influence the frequency of recombination. Therefore, it is important to know the position of the colour locus, gene order and recombination landscape of the chromosome to understand the mechanism that generates incompatibilities between morphs. Recent studies showed remarkable pre- and postzygotic incompatibilities between sympatric colour morphs of the Gouldian finch (Erythrura gouldiae), in which head feather colour is genetically determined by a single sex-linked locus, Red. We constructed a genetic map for the Z chromosome of the Gouldian finch (male-specific map distance=131 cM), using 618 captive-bred birds and 34 microsatellite markers, to investigate the extent of inter- and intraspecific genomic rearrangements and variation in recombination rate within the Z chromosome. We refined the location of the Red locus to a ~7.2-cM interval in a region with a moderate recombination rate but outside the least-recombining, putative centromeric region. There was no evidence of chromosome-wide genomic rearrangements between the chromosomes carrying the red or black alleles with the current marker resolution. This work will contribute to identifying the causal gene, which will in turn enable alternative explanations for the association between incompatibility and colouration, such as fine-scale linkage disequilibrium, genomic rearrangements and pleiotropy, to be tested.

Inbreeding drives maize centromere evolution

Functional centromeres, the chromosomal sites of spindle attachment during cell division, are marked epigenetically by the centromere-specific histone H3 variant cenH3 and typically contain long stretches of centromere-specific tandem DNA repeats (∼1.8 Mb in maize). In 23 inbreds of domesticated maize chosen to represent the genetic diversity of maize germplasm, partial or nearly complete loss of the tandem DNA repeat CentC precedes 57 independent cenH3 relocation events that result in neocentromere formation. Chromosomal regions with newly acquired cenH3 are colonized by the centromere-specific retrotransposon CR2 at a rate that would result in centromere-sized CR2 clusters in 20,000-95,000 y. Three lines of evidence indicate that CentC loss is linked to inbreeding, including (i) CEN10 of temperate lineages, presumed to have experienced a genetic bottleneck, contain less CentC than their tropical relatives; (ii) strong selection for centromere-linked genes in domesticated maize reduced diversity at seven of the ten maize centromeres to only one or two postdomestication haplotypes; and (iii) the centromere with the largest number of haplotypes in domesticated maize (CEN7) has the highest CentC levels in nearly all domesticated lines. Rare recombinations introduced one (CEN2) or more (CEN5) alternate CEN haplotypes while retaining a single haplotype at domestication loci linked to these centromeres. Taken together, this evidence strongly suggests that inbreeding, favored by postdomestication selection for centromere-linked genes affecting key domestication or agricultural traits, drives replacement of the tandem centromere repeats in maize and other crop plants. Similar forces may act during speciation in natural systems.

Inbreeding drives maize centromere evolution

Functional centromeres, the chromosomal sites of spindle attachment during cell division, are marked epigenetically by the centromere-specific histone H3 variant cenH3 and typically contain long stretches of centromere-specific tandem DNA repeats (∼1.8 Mb in maize). In 23 inbreds of domesticated maize chosen to represent the genetic diversity of maize germplasm, partial or nearly complete loss of the tandem DNA repeat CentC precedes 57 independent cenH3 relocation events that result in neocentromere formation. Chromosomal regions with newly acquired cenH3 are colonized by the centromere-specific retrotransposon CR2 at a rate that would result in centromere-sized CR2 clusters in 20,000-95,000 y. Three lines of evidence indicate that CentC loss is linked to inbreeding, including (i) CEN10 of temperate lineages, presumed to have experienced a genetic bottleneck, contain less CentC than their tropical relatives; (ii) strong selection for centromere-linked genes in domesticated maize reduced diversity at seven of the ten maize centromeres to only one or two postdomestication haplotypes; and (iii) the centromere with the largest number of haplotypes in domesticated maize (CEN7) has the highest CentC levels in nearly all domesticated lines. Rare recombinations introduced one (CEN2) or more (CEN5) alternate CEN haplotypes while retaining a single haplotype at domestication loci linked to these centromeres. Taken together, this evidence strongly suggests that inbreeding, favored by postdomestication selection for centromere-linked genes affecting key domestication or agricultural traits, drives replacement of the tandem centromere repeats in maize and other crop plants. Similar forces may act during speciation in natural systems.

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres - an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half-tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination - a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high-density maps in species with strong recombination interference and will enrich many existing and future mapping resources.

Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae)

BACKGROUND

The evolution of species boundaries and the relative impact of selection and gene flow on genomic divergence are best studied in populations and species pairs exhibiting various levels of divergence along the speciation continuum. We studied species boundaries in Iberian barbels, Barbus and Luciobarbus, a system of populations and species spanning a wide degree of genetic relatedness, as well as geographic distribution and range overlap. We jointly analyze multiple types of molecular markers and morphological traits to gain a comprehensive perspective on the nature of species boundaries in these cyprinid fishes.

RESULTS

Intraspecific molecular and morphological differentiation is visible among many populations. Genomes of all sympatric species studied are porous to gene flow, even if they are not sister species. Compared to their allopatric counterparts, sympatric representatives of different species share alleles and show an increase in all measures of nucleotide polymorphism (S, Hd, K, π and θ). High molecular diversity is particularly striking in L. steindachneri from the Tejo and Guadiana rivers, which co-varies with other sympatric species. Interestingly, different nuclear markers introgress across species boundaries at various levels, with distinct impacts on population trees. As such, some loci exhibit limited introgression and population trees resemble the presumed species tree, while alleles at other loci introgress more freely and population trees reflect geographic affinities and interspecific gene flow. Additionally, extent of introgression decreases with increasing genetic divergence in hybridizing species pairs.

CONCLUSIONS

We show that reproductive isolation in Iberian Barbus and Luciobarbus is not complete and species boundaries are semi-permeable to (some) gene flow, as different species (including non-sister) are exchanging genes in areas of sympatry. Our results support a speciation-with-gene-flow scenario with heterogeneous barriers to gene flow across the genome, strengthening with genetic divergence. This is consistent with observations coming from other systems and supports the notion that speciation is not instantaneous but a gradual process, during which different species are still able to exchange some genes, while selection prevents gene flow at other loci. We also provide evidence for a hybrid origin of a barbel ecotype, L. steindachneri, suggesting that ecology plays a key role in species coexistence and hybridization in Iberian barbels. This ecotype with intermediate, yet variable, molecular, morphological, trophic and ecological characteristics is the local product of introgressive hybridization of L. comizo with up to three different species (with L. bocagei in the Tejo, with L. microcephalus and L. sclateri in the Guadiana). In spite of the homogenizing effects of ongoing gene flow, species can still be discriminated using a combination of morphological and molecular markers. Iberian barbels are thus an ideal system for the study of species boundaries, since they span a wide range of genetic divergences, with diverse ecologies and degrees of sympatry.

Spatial genetic structure across a hybrid zone between European rabbit subspecies

The Iberian Peninsula is the only region in the world where the two existing subspecies of the European rabbit (Oryctolagus cuniculus) naturally occur and hybridize. In this study we explore the relative roles of historical and contemporary processes in shaping the spatial genetic structure of the rabbit across its native distribution range, and how they differently affect each subspecies and the hybrid zone. For that purpose we obtained multilocus genotypes and mitochondrial DNA data from 771 rabbits across most of the distribution range of the European rabbit in Spain. Based on the nuclear markers we observed a hierarchical genetic structure firstly comprised by two genetic groups, largely congruent with the mitochondrial lineages and subspecies distributions (O. c. algirus and O. c. cuniculus), which were subsequently subdivided into seven genetic groups. Geographic distance alone emerged as an important factor explaining genetic differentiation across the whole range, without the need to invoke for the effect for geographical barriers. Additionally, the significantly positive spatial correlation up to a distance of only 100 km supported the idea that differentiation at a local level is of greater importance when considering the species overall genetic structure. When looking at the subspecies, northern populations of O. c. cuniculus showed more spatial genetic structure and differentiation than O. c. algirus. This could be due to local geographic barriers, limited resources, soil type and/or social behavior limiting dispersal. The hybrid zone showed similar genetic structure to the southern populations but a larger introgression from the northern lineage genome. These differences have been attributed to selection against the hybrids rather than to behavioral differences between subspecies. Ultimately, the genetic structure of the rabbit in its native distribution range is the result of an ensemble of factors, from geographical and ecological, to behavioral and molecular, that hierarchically interact through time and space.

Effective population size does not predict codon usage bias in mammals

Synonymous codons are not used at equal frequency throughout the genome, a phenomenon termed codon usage bias (CUB). It is often assumed that interspecific variation in the intensity of CUB is related to species differences in effective population sizes (N_e), with selection on CUB operating less efficiently in species with small N_e. Here, we specifically ask whether variation in N_e predicts differences in CUB in mammals and report two main findings. First, across 41 mammalian genomes, CUB was not correlated with two indirect proxies of N_e (body mass and generation time), even though there was statistically significant evidence of selection shaping CUB across all species. Interestingly, autosomal genes showed higher codon usage bias compared to X-linked genes, and high-recombination genes showed higher codon usage bias compared to low recombination genes, suggesting intraspecific variation in N_e predicts variation in CUB. Second, across six mammalian species with genetic estimates of N_e (human, chimpanzee, rabbit, and three mouse species: Mus musculus, M. domesticus, and M. castaneus), N_e and CUB were weakly and inconsistently correlated. At least in mammals, interspecific divergence in N_e does not strongly predict variation in CUB. One hypothesis is that each species responds to a unique distribution of selection coefficients, confounding any straightforward link between N_e and CUB.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

The genomic architecture of population divergence between subspecies of the European rabbit

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

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. FST 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.

A role for migration-linked genes and genomic islands in divergence of a songbird

Next-generation sequencing has made it possible to begin asking questions about the process of divergence at the level of the genome. For example, recently, there has been a debate around the role of 'genomic islands of divergence' (i.e. blocks of outlier loci) in facilitating the process of speciation-with-gene-flow. The Swainson's thrush, Catharus ustulatus, is a migratory songbird with two genetically distinct subspecies that differ in a number of traits known to be involved in reproductive isolation in birds (plumage coloration, song and migratory behaviour), despite contemporary gene flow along a secondary contact zone. Here, we use RAD-PE sequencing to test emerging hypotheses about the process of divergence at the level of the genome and identify genes and gene regions involved in differentiation in this migratory songbird. Our analyses revealed distinct genomic islands on 15 of the 23 chromosomes and an accelerated rate of divergence on the Z chromosome, one of the avian sex chromosomes. Further, an analysis of loci linked to traits known to be involved in reproductive isolation in songbirds showed that genes linked to migration are significantly more differentiated than expected by chance, but that these genes lie primarily outside the genomic islands. Overall, our analysis supports the idea that genes linked to migration play an important role in divergence in migratory songbirds, but we find no compelling evidence that the observed genomic islands are facilitating adaptive divergence in migratory behaviour.

Adaptive introgression between Anopheles sibling species eliminates a major genomic island but not reproductive isolation

Adaptive introgression can provide novel genetic variation to fuel rapid evolutionary responses, though it may be counterbalanced by potential for detrimental disruption of the recipient genomic background. We examine the extent and impact of recent introgression of a strongly selected insecticide-resistance mutation (Vgsc-1014F) located within one of two exceptionally large genomic islands of divergence separating the Anopheles gambiae species pair. Here we show that transfer of the Vgsc mutation results in homogenization of the entire genomic island region (~1.5% of the genome) between species. Despite this massive disruption, introgression is clearly adaptive with a dramatic rise in frequency of Vgsc-1014F and no discernable impact on subsequent reproductive isolation between species. Our results show (1) how resilience of genomes to massive introgression can permit rapid adaptive response to anthropogenic selection and (2) that even extreme prominence of genomic islands of divergence can be an unreliable indicator of importance in speciation.

Highly divergent genomic islands segregate between a species pair of the mosquito, Anopheles gambiae. Here Clarkson et al. show that loss of one of the largest islands, driven by adaptive introgression of an insecticide-resistance mutation, had no impact on reproductive isolation.

Speciation and introgression between Mimulus nasutus and Mimulus guttatus

Mimulus guttatus and M. nasutus are an evolutionary and ecological model sister species pair differentiated by ecology, mating system, and partial reproductive isolation. Despite extensive research on this system, the history of divergence and differentiation in this sister pair is unclear. We present and analyze a population genomic data set which shows that M. nasutus budded from a central Californian M. guttatus population within the last 200 to 500 thousand years. In this time, the M. nasutus genome has accrued genomic signatures of the transition to predominant selfing, including an elevated proportion of nonsynonymous variants, an accumulation of premature stop codons, and extended levels of linkage disequilibrium. Despite clear biological differentiation, we document genomic signatures of ongoing, bidirectional introgression. We observe a negative relationship between the recombination rate and divergence between M. nasutus and sympatric M. guttatus samples, suggesting that selection acts against M. nasutus ancestry in M. guttatus.

Recovering power in association mapping panels with variable levels of linkage disequilibrium

Association mapping has permitted the discovery of major QTL in many species. It can be applied to existing populations and, as a consequence, it is generally necessary to take into account structure and relatedness among individuals in the statistical model to control false positives. We analytically studied power in association studies by computing noncentrality parameter of the tests and its relationship with parameters characterizing diversity (genetic differentiation between groups and allele frequencies) and kinship between individuals. Investigation of three different maize diversity panels genotyped with the 50k SNPs array highlighted contrasted average power among panels and revealed gaps of power of classical mixed models in regions with high linkage disequilibrium (LD). These gaps could be related to the fact that markers are used for both testing association and estimating relatedness. We thus considered two alternative approaches to estimating the kinship matrix to recover power in regions of high LD. In the first one, we estimated the kinship with all the markers that are not located on the same chromosome than the tested SNP. In the second one, correlation between markers was taken into account to weight the contribution of each marker to the kinship. Simulations revealed that these two approaches were efficient to control false positives and were more powerful than classical models.

Genome-wide dissection of hybrid sterility in Drosophila confirms a polygenic threshold architecture

To date, different studies about the genetic basis of hybrid male sterility (HMS), a postzygotic reproductive barrier thoroughly investigated using Drosophila species, have demonstrated that no single major gene can produce hybrid sterility without the cooperation of several genetic factors. Early work using hybrids between Drosophila koepferae (Dk) and Drosophila buzzatii (Db) was consistent with the idea that HMS requires the cooperation of several genetic factors, supporting a polygenic threshold (PT) model. Here we present a genome-wide mapping strategy to test the PT model, analyzing serially backcrossed fertile and sterile males in which the Dk genome was introgressed into the Db background. We identified 32 Dk-specific markers significantly associated with hybrid sterility. Our results demonstrate 1) a strong correlation between the number of segregated sterility markers and males' degree of sterility, 2) the exchangeability among markers, 3) their tendency to cluster into low-recombining chromosomal regions, and 4) the requirement for a minimum number (threshold) of markers to elicit sterility. Although our findings do not contradict a role for occasional major hybrid-sterility genes, they conform more to the view that HMS primarily evolves by the cumulative action of many interacting genes of minor effect in a complex PT architecture.