The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene.

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes.

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Long-Term Balancing Selection and the Genetic Load Linked to the Self-Incompatibility Locus in Arabidopsis halleri and A. lyrata.

Balancing selection is a form of natural selection maintaining diversity at the sites it targets and at linked nucleotide sites. Due to selection favoring heterozygosity, it has the potential to facilitate the accumulation of a "sheltered" load of tightly linked recessive deleterious mutations. However, precisely evaluating the extent of these effects has remained challenging. Taking advantage of plant self-incompatibility as one of the best-understood examples of long-term balancing selection, we provide a highly resolved picture of the genomic extent of balancing selection on the sheltered genetic load. We used targeted genome resequencing to reveal polymorphism of the genomic region flanking the self-incompatibility locus in three sample sets in each of the two closely related plant species Arabidopsis halleri and Arabidopsis lyrata, and used 100 control regions from throughout the genome to factor out differences in demographic histories and/or sample structure. Nucleotide polymorphism increased strongly around the S-locus in all sample sets, but only over a limited genomic region, as it became indistinguishable from the genomic background beyond the first 25-30 kb. Genes in this chromosomal interval exhibited no excess of mutations at 0-fold degenerated sites relative to putatively neutral sites, hence revealing no detectable weakening of the efficacy of purifying selection even for these most tightly linked genes. Overall, our results are consistent with the predictions of a narrow genomic influence of linkage to the S-locus and clarify how natural selection in one genomic region affects the evolution of the adjacent genomic regions.

Inferring the Demographic History and Inheritance Mode of Tetraploid Species Using ABC.

Genomic patterns of diversity and divergence are impacted by certain life history traits, reproductive systems, and demographic history. The latter is characterized by fluctuations in population sizes over time, as well as by temporal patterns of introgression. For a given organism, identifying a demographic history that deviates from the standard neutral model allows a better understanding of its evolution but also helps to reduce the risk of false positives when screening for molecular targets of natural selection. Tetraploid organisms and beyond have demographic histories that are complicated by the mode of polyploidization, the mode of inheritance, and different scenarios of gene flow between sub-genomes and diploid parental species. Here we provide guidelines for experimenters wishing to address these issues through a flexible statistical framework: approximate Bayesian computation (ABC). The emphasis is on the general philosophy of the approach to encourage future users to exploit the enormous flexibility of ABC beyond the limitations imposed by generalist data analysis pipelines.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms.

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

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.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population.

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

DILS: Demographic inferences with linked selection by using ABC.

We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single-population or two-population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.

Evolution of self-incompatibility in the Brassicaceae: Lessons from a textbook example of natural selection.

Self-incompatibility (SI) is a self-recognition genetic system enforcing outcrossing in hermaphroditic flowering plants and results in one of the arguably best understood forms of natural (balancing) selection maintaining genetic variation over long evolutionary times. A rich theoretical and empirical population genetics literature has considerably clarified how the distribution of SI phenotypes translates into fitness differences among individuals by a combination of inbreeding avoidance and rare-allele advantage. At the same time, the molecular mechanisms by which self-pollen is specifically recognized and rejected have been described in exquisite details in several model organisms, such that the genotype-to-phenotype map is also pretty well understood, notably in the Brassicaceae. Here, we review recent advances in these two fronts and illustrate how the joint availability of detailed characterization of genotype-to-phenotype and phenotype-to-fitness maps on a single genetic system (plant self-incompatibility) provides the opportunity to understand the evolutionary process in a unique perspective, bringing novel insight on general questions about the emergence, maintenance, and diversification of a complex genetic system.

Base-Pairing Requirements for Small RNA-Mediated Gene Silencing of Recessive Self-Incompatibility Alleles in Arabidopsis halleri.

Small noncoding RNAs are central regulators of genome activity and stability. Their regulatory function typically involves sequence similarity with their target sites, but understanding the criteria by which they specifically recognize and regulate their targets across the genome remains a major challenge in the field, especially in the face of the diversity of silencing pathways involved. The dominance hierarchy among self-incompatibility alleles in Brassicaceae is controlled by interactions between a highly diversified set of small noncoding RNAs produced by dominant S-alleles and their corresponding target sites on recessive S-alleles. By controlled crosses, we created numerous heterozygous combinations of S-alleles in Arabidopsis halleri and developed an real-time quantitative PCR assay to compare allele-specific transcript levels for the pollen determinant of self-incompatibility (SCR). This provides the unique opportunity to evaluate the precise base-pairing requirements for effective transcriptional regulation of this target gene. We found strong transcriptional silencing of recessive SCR alleles in all heterozygote combinations examined. A simple threshold model of base pairing for the small RNA-target interaction captures most of the variation in SCR transcript levels. For a subset of S-alleles, we also measured allele-specific transcript levels of the determinant of pistil specificity (SRK), and found sharply distinct expression dynamics throughout flower development between SCR and SRK In contrast to SCR, both SRK alleles were expressed at similar levels in the heterozygote genotypes examined, suggesting no transcriptional control of dominance for this gene. We discuss the implications for the evolutionary processes associated with the origin and maintenance of the dominance hierarchy among self-incompatibility alleles.

Diversification dynamics of freshwater bivalves (Unionidae: Parreysiinae: Coelaturini) indicate historic hydrographic connections throughout the East African Rift System.

Invertebrates are exceptionally diverse, but many are in decline because of anthropogenic changes to their habitat. This situation is particularly problematic for taxa that are not well monitored or taxonomically poorly understood, because the lack of knowledge hampers conservation. Despite their important functional role in freshwater ecosystems, African bivalves of the family Unionidae remain poorly studied compared to their highly threatened relatives in Europe, the U.S.A. and Canada. To resolve relationships and to study diversification dynamics in space and time, we performed time-calibrated phylogenetic studies and biogeographical modeling on the unionids from the East African Rift System and surroundings, including representatives of all currently recognized Afrotropical genera except for Brazzaea (and Unio from southern Africa). Our analyses indicate that all sampled Afrotropical unionids belong to the tribe Coelaturini (subfamily Parreysiinae), as does the genus Moncetia from Lake Tanganyika, which is currently attributed to the family Iridinidae. Colonization of Africa from Eurasia by Parreysiinae occurred ~17 Ma ago, and the subsequent diversification of Coelaturini in Africa continued at a steady pace, although net diversification decreased over time as more niches and ecoregions became occupied. Clades in Coelaturini largely reflect drainage basins, with the oldest lineages and highest regional diversity occurring in Lake Tanganyika, followed by the Congo Basin watershed in general. The species assemblage of Lake Tanganyika reflects multiple independent events of colonization and intralacustrine diversification since the Late Miocene or Early Pliocene. The clades of other regions, including that containing the species from Lake Malawi, are comparatively young. Biogeographical analyses indicate that the colonization history was mainly driven by cladogenesis in sympatry, whereas few anagenetic events contributed to the modern distribution of Coelaturini. Ancestral range estimations demonstrate that Coelaturini originated in the Victoria and/or Tanganyika ecoregions, and that the Congo Basin played an essential role in the colonization of Africa by Coelaturini.

Impact of whole genome triplication on the evolutionary history and the functional dynamics of regulatory genes involved in Brassica self-incompatibility signalling pathway.

Polyploidy or whole genome duplication is a frequent and recurrent phenomenon in flowering plants that has played a major role in their diversification, adaptation and speciation. The adaptive success of polyploids relates to the different evolutionary fates of duplicated genes. In this study, we explored the impact of the whole genome triplication (WGT) event in the Brassiceae tribe on the genes involved in the self-incompatibility (SI) signalling pathway, a mechanism allowing recognition and rejection of self-pollen in hermaphrodite plants. By taking advantage of the knowledge acquired on this pathway as well as of several reference genomes in Brassicaceae species, we determined copy number of the different genes involved in this pathway and investigated their structural and functional evolutionary dynamics. We could infer that whereas most genes involved in the SI signalling returned to single copies after the WGT event (i.e. ARC1, JDP1, THL1, THL2, Exo70A01) in diploid Brassica species, a few were retained in duplicated (GLO1 and PLDα) or triplicated copies (MLPK). We also carefully studied the gene structure of these latter duplicated genes (including the conservation of functional domains and active sites) and tested their transcription in the stigma to identify which copies seem to be involved in the SI signalling pathway. By taking advantage of these analyses, we then explored the putative origin of a contrasted SI phenotype between two Brassica rapa varieties that have been fully sequenced and shared the same S-allele (S60).

Adaptive divergence in shell morphology in an ongoing gastropod radiation from Lake Malawi.

BACKGROUND: Ecological speciation is a prominent mechanism of diversification but in many evolutionary radiations, particularly in invertebrates, it remains unclear whether supposedly critical ecological traits drove or facilitated diversification. As a result, we lack accurate knowledge on the drivers of diversification for most evolutionary radiations along the tree of life. Freshwater mollusks present an enigmatic example: Putatively adaptive radiations are being described in various families, typically from long-lived lakes, whereas other taxa represent celebrated model systems in the study of ecophenotypic plasticity. Here we examine determinants of shell-shape variation in three nominal species of an ongoing ampullariid radiation in the Malawi Basin (Lanistes nyassanus, L. solidus and Lanistes sp. (ovum-like)) with a common garden experiment and semi-landmark morphometrics. RESULTS: We found significant differences in survival and fecundity among these species in contrasting habitats. Morphological differences observed in the wild persisted in our experiments for L. nyassanus versus L. solidus and L. sp. (ovum-like), but differences between L. solidus and L. sp. (ovum-like) disappeared and re-emerged in the F1 and F2 generations, respectively. These results indicate that plasticity occurred, but that it is not solely responsible for the observed differences. Our experiments provide the first unambiguous evidence for genetic divergence in shell morphology in an ongoing freshwater gastropod radiation in association with marked fitness differences among species under controlled habitat conditions. CONCLUSIONS: Our results indicate that differences in shell morphology among Lanistes species occupying different habitats have an adaptive value. These results also facilitate an accurate reinterpretation of morphological variation in fossil Lanistes radiations, and thus macroevolutionary dynamics. Finally, our work testifies that the shells of freshwater gastropods may retain signatures of adaptation at low taxonomic levels, beyond representing an evolutionary novelty responsible for much of the diversity and disparity in mollusks altogether.

Genotyping and De Novo Discovery of Allelic Variants at the Brassicaceae Self-Incompatibility Locus from Short-Read Sequencing Data.

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Asymmetrical diversification of the receptor-ligand interaction controlling self-incompatibility in Arabidopsis.

How two-component genetic systems accumulate evolutionary novelty and diversify in the course of evolution is a fundamental problem in evolutionary systems biology. In the Brassicaceae, self-incompatibility (SI) is a spectacular example of a diversified allelic series in which numerous highly diverged receptor-ligand combinations are segregating in natural populations. However, the evolutionary mechanisms by which new SI specificities arise have remained elusive. Using in planta ancestral protein reconstruction, we demonstrate that two allelic variants segregating as distinct receptor-ligand combinations diverged through an asymmetrical process whereby one variant has retained the same recognition specificity as their (now extinct) putative ancestor, while the other has functionally diverged and now represents a novel specificity no longer recognized by the ancestor. Examination of the structural determinants of the shift in binding specificity suggests that qualitative rather than quantitative changes of the interaction are an important source of evolutionary novelty in this highly diversified receptor-ligand system.

Chromosome-scale assemblies of plant genomes using nanopore long reads and optical maps.

Plant genomes are often characterized by a high level of repetitiveness and polyploid nature. Consequently, creating genome assemblies for plant genomes is challenging. The introduction of short-read technologies 10 years ago substantially increased the number of available plant genomes. Generally, these assemblies are incomplete and fragmented, and only a few are at the chromosome scale. Recently, Pacific Biosciences and Oxford Nanopore sequencing technologies were commercialized that can sequence long DNA fragments (kilobases to megabase) and, using efficient algorithms, provide high-quality assemblies in terms of contiguity and completeness of repetitive regions1-4. However, even though genome assemblies based on long reads exhibit high contig N50s (>1 Mb), these methods are still insufficient to decipher genome organization at the chromosome level. Here, we describe a strategy based on long reads (MinION or PromethION sequencers) and optical maps (Saphyr system) that can produce chromosome-level assemblies and demonstrate applicability by generating high-quality genome sequences for two new dicotyledon morphotypes, Brassica rapa Z1 (yellow sarson) and Brassica oleracea HDEM (broccoli), and one new monocotyledon, Musa schizocarpa (banana). All three assemblies show contig N50s of >5 Mb and contain scaffolds that represent entire chromosomes or chromosome arms.

Immune failure reveals vulnerability of populations exposed to pollution in the bioindicator species Hediste diversicolor.

Human activities on the shoreline generate a growing pollution, creating deleterious habitats in coastal zones. Some species nevertheless succeed in such harsh milieus, raising the question of their tolerance to environmental stress. The annelid Hediste diversicolor lives buried in the sediments, directly exposed to contaminants trapped in the mud. After verifying the similarity of their genetic contexts, we compared reproductive output and individual immune resistance measures of populations living in polluted vs. 'clean' sediments, and related these assessments with measures of phthalates and metal pollution, and associated toxicity indices. Chemical analyses predicted no toxicity to the local infauna, and phenological studies evidenced no direct cost of living in noxious habitats. However, populations exposed to pollutants showed a significantly reduced survival upon infection with a local pathogen. Surprisingly, physiological studies evidenced a basal overinflammatory state in the most exposed populations. This over-activated baseline immune phenotype likely generates self-damage leading to enhanced immune cell death rate and immune failure. Monitoring the immune status of individual worms living in anthropic areas could thus be used as a reliable source of information regarding the actual health of wild populations.

Controlling for genetic identity of varieties, pollen contamination and stigma receptivity is essential to characterize the self-incompatibility system of Olea europaea L.

Bervillé et al. express concern about the existence of the diallelic self-incompatibility (DSI) system in Olea europaea, mainly because our model does not account for results from previous studies from their group that claimed to have documented asymmetry of the incompatibility response in reciprocal crosses. In this answer to their comment, we present original results based on reciprocal stigma tests that contradict conclusions from these studies. We show that, in our hands, not a single case of asymmetry was confirmed, endorsing that symmetry of incompatibility reactions seems to be the rule in Olive. We discuss three important aspects that were not taken into account in the studies cited in their comments and that can explain the discrepancy: (i) the vast uncertainty around the actual genetic identity of vernacular varieties, (ii) the risk of massive contamination associated with the pollination protocols that they used and (iii) the importance of checking for stigma receptivity in controlled crosses. These studies were thus poorly genetically controlled, and we stand by our original conclusion that Olive tree exhibits DSI.

Elucidation of the genetic architecture of self-incompatibility in olive: Evolutionary consequences and perspectives for orchard management.

The olive (Olea europaea L.) is a typical important perennial crop species for which the genetic determination and even functionality of self-incompatibility (SI) are still largely unresolved. It is still not known whether SI is under gametophytic or sporophytic genetic control, yet fruit production in orchards depends critically on successful ovule fertilization. We studied the genetic determination of SI in olive in light of recent discoveries in other genera of the Oleaceae family. Using intra- and interspecific stigma tests on 89 genotypes representative of species-wide olive diversity and the compatibility/incompatibility reactions of progeny plants from controlled crosses, we confirmed that O. europaea shares the same homomorphic diallelic self-incompatibility (DSI) system as the one recently identified in Phillyrea angustifolia and Fraxinus ornus. SI is sporophytic in olive. The incompatibility response differs between the two SI groups in terms of how far pollen tubes grow before growth is arrested within stigma tissues. As a consequence of this DSI system, the chance of cross-incompatibility between pairs of varieties in an orchard is high (50%) and fruit production may be limited by the availability of compatible pollen. The discovery of the DSI system in O. europaea will undoubtedly offer opportunities to optimize fruit production.