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

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

DNA Origami-Enabled Gene Localization of Repetitive Sequences

Repetitive sequences, which make up over 50% of human DNA, have diverse applications in disease diagnosis, forensic identification, paternity testing, and population genetic analysis due to their crucial functions for gene regulation. However, representative detection technologies such as sequencing and fluorescence imaging suffer from time-consuming protocols, high cost, and inaccuracy of the position and order of repetitive sequences. Here, we develop a precise and cost-effective strategy that combines the high resolution of atomic force microscopy with the shape customizability of DNA origami for repetitive sequence-specific gene localization. "Tri-block" DNA structures were specifically designed to connect repetitive sequences to DNA origami tags, thereby revealing precise genetic information in terms of position and sequence for high-resolution and high-precision visualization of repetitive sequences. More importantly, we achieved the results of simultaneous detection of different DNA repetitive sequences on the gene template with a resolution of ∼6.5 nm (19 nt). This strategy is characterized by high efficiency, high precision, low operational complexity, and low labor/time costs, providing a powerful complement to sequencing technologies for gene localization of repetitive sequences.

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.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

S-Locus Genotyping in Japanese Plum by High Throughput Sequencing Using a Synthetic S-Loci Reference Sequence

Self-incompatibility in Prunus species is governed by a single locus consisting of two highly multi-allelic and tightly linked genes, one coding for an F-box protein-i.e., SFB in Prunus- controlling the pollen specificity and one coding for an S-RNase gene controlling the pistil specificity. Genotyping the allelic combination in a fruit tree species is an essential procedure both for cross-based breeding and for establishing pollination requirements. Gel-based PCR techniques using primer pairs designed from conserved regions and spanning polymorphic intronic regions are traditionally used for this task. However, with the great advance of massive sequencing techniques and the lowering of sequencing costs, new genotyping-by-sequencing procedures are emerging. The alignment of resequenced individuals to reference genomes, commonly used for polymorphism detection, yields little or no coverage in the S-locus region due to high polymorphism between different alleles within the same species, and cannot be used for this purpose. Using the available sequences of Japanese plum S-loci concatenated in a rosary-like structure as synthetic reference sequence, we describe a procedure to accurately genotype resequenced individuals that allowed the analysis of the S-genotype in 88 Japanese plum cultivars, 74 of them are reported for the first time. In addition to unraveling two new S-alleles from published reference genomes, we identified at least two S-alleles in 74 cultivars. According to their S-allele composition, they were assigned to 22 incompatibility groups, including nine new incompatibility groups reported here for the first time (XXVII-XXXV).

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

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.

The Scope for Postmating Sexual Selection in Plants

Sexual selection is known to shape plant traits that affect access to mates during the pollination phase, but it is less well understood to what extent it affects traits relevant to interactions between pollen and pistils after pollination. This is surprising, because both of the two key modes of sexual selection, male-male competition and female choice, could plausibly operate during pollen-pistil interactions where physical male-female contact occurs. Here, we consider how the key processes of sexual selection might affect traits involved in pollen-pistil interactions, including 'Fisherian runaway' and 'good-genes' models. We review aspects of the molecular and cellular biology of pollen-pistil interactions on which sexual selection could act and point to research that is needed to investigate them.

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.

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.