Benjamin Laenen

Rapid Evolution of Genomic Imprinting in Two Species of the Brassicaceae

Genomic imprinting is a rapidly evolving phenomenon in plants; however, for those genes where imprinting conveys an advantage, imprinted expression is maintained over long evolutionary timescales. Genomic imprinting is an epigenetic phenomenon occurring in mammals and flowering plants that causes genes to adopt a parent-of-origin-specific mode of expression. While the imprinting status of genes is well conserved in mammals, clear estimates for the degree of conservation were lacking in plants. We therefore analyzed the genome-wide imprinting status of Capsella rubella, which shared a common recent ancestor with Arabidopsis thaliana ∼10 to 14 million years ago. However, only ∼14% of maternally expressed genes (MEGs) and ∼29% of paternally expressed genes (PEGs) in C. rubella were commonly imprinted in both species, revealing that genomic imprinting is a rapidly evolving phenomenon in plants. Nevertheless, conserved PEGs exhibited signs of selection, suggesting that a subset of imprinted genes play an important functional role and are therefore maintained in plants. Like in Arabidopsis, PEGs in C. rubella are frequently associated with the presence of transposable elements that preferentially belong to helitron and MuDR families. Our data further reveal that MEGs and PEGs differ in their targeting by 24-nucleotide small RNAs and asymmetric DNA methylation, suggesting different mechanisms establishing DNA methylation at MEGs and PEGs.

Phylogeny, classification and species delimitation in the liverwort genus Odontoschisma (Cephaloziaceae)

A species-level phylogeny of Odontoschisma (Jungermanniidae: Cephaloziaceae) was produced to revisit the infrageneric classification and delimit the species within the genus. New methods of species delimitation have been used to explicitly contrast taxonomic hypotheses and test the relevance of the morphological traits traditionally used in this group. The results confirm previous evidence suggesting that the circumscription of Odontoschisma needs to be enlarged to include Iwatsukia and Cladopodiella, and further indicate that a third genus, Anomoclada, is nested within it. Twentythree molecular entities were recognized based the results of generalized mixed Yule-coalescent (GMYC) analyses. These entities partly conflicted with traditionally defined species that were shown to belong to different, not necessarily closely related entities, adding to the growing body of evidence calling for an extensive revision of species delimitations in taxa with reduced morphologies like leafy liverworts, using an integrative taxonomic approach. A fully revised classification of Odontoschisma into five sections, twenty species, and three subspecies is presented. While the species reported from Europe, Asia and North America were of polyphyletic origins, all Neotropical species were resolved as monophyletic, which could result from a combination of fast speciation rates and reduced dispersal in the Neotropics, and potential extinction in other areas, especially sub-Saharan Africa.

Genomic analysis reveals major determinants of cis-regulatory variation in Capsella grandiflora

Significance Despite long-standing interest in the contribution of cis-regulatory changes to adaptation, we still have a limited understanding of the selective importance and genomic determinants of cis-regulatory variation in natural populations. We use a combination of analyses of allele-specific expression and population genomic analyses to investigate the selective forces and genomic determinants of cis-regulatory variation in the outcrossing plant species Capsella grandiflora. We conclude that gene-specific functional constraints shape cis-regulatory variation and that genes with cis-regulatory variation are under relaxed purifying selection compared with other genes. Finally, we identify a link between gene body methylation and the extent of cis-regulatory constraint in natural populations. Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora. We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Early burst in body size evolution is uncoupled from species diversification in diving beetles (Dytiscidae)

Changes in morphology are often thought to be linked to changes in species diversification, which is expected to leave a signal of early burst (EB) in phenotypic traits. However, such signal is rarely recovered in empirical phylogenies, even for groups with well‐known adaptive radiation. Using a comprehensive phylogenetic approach in Dytiscidae, which harbours ~4,300 species with as much as 50‐fold variation in body size among them, we ask whether pattern of species diversification correlates with morphological evolution. Additionally, we test whether the large variation in body size is linked to habitat preference and whether the latter influences species turnover. We found, in sharp contrast to most animal groups, that Dytiscidae body size evolution follows an early‐burst model with subsequent high phylogenetic conservatism. However, we found no evidence for associated shifts in species diversification, which point to an uncoupled evolution of morphology and species diversification. We recovered the ancestral habitat of Dytiscidae as lentic (standing water), with many transitions to lotic habitat (running water) that are concomitant to a decrease in body size. Finally, we found no evidence for difference in net diversification rates between habitats nor difference in turnover in lentic and lotic species. This result, together with recent findings in dragonflies, contrasts with some theoretical expectations of the habitat stability hypothesis. Thus, a thorough reassessment of the impact of dispersal, gene flow and range size on the speciation process is needed to fully encompass the evolutionary consequences of the lentic–lotic divide for freshwater fauna.

Range size heritability and diversification patterns in the liverwort genus Radula

Why some species exhibit larger geographical ranges than others, and to what extent does variation in range size affect diversification rates, remains a fundamental, but largely unanswered question in ecology and evolution. Here, we implement phylogenetic comparative analyses and ancestral area estimations in Radula, a liverwort genus of Cretaceous origin, to investigate the mechanisms that explain differences in geographical range size and diversification rates among lineages. Range size was phylogenetically constrained in the two sub-genera characterized by their almost complete Australasian and Neotropical endemicity, respectively. The congruence between the divergence time of these lineages and continental split suggests that plate tectonics could have played a major role in their present distribution, suggesting that a strong imprint of vicariance can still be found in extant distribution patterns in these highly mobile organisms. Amentuloradula, Volutoradula and Metaradula species did not appear to exhibit losses of dispersal capacities in terms of dispersal life-history traits, but evidence for significant phylogenetic signal in macroecological niche traits suggests that niche conservatism accounts for their restricted geographic ranges. Despite their greatly restricted distribution to Australasia and Neotropics respectively, Amentuloradula and Volutoradula did not exhibit significantly lower diversification rates than more widespread lineages, in contrast with the hypothesis that the probability of speciation increases with range size by promoting geographic isolation and increasing the rate at which novel habitats are encountered. We suggest that stochastic long-distance dispersal events may balance allele frequencies across large spatial scales, leading to low genetic structure among geographically distant areas or even continents, ultimately decreasing the diversification rates in highly mobile, widespread lineages.

High migration rates shape the postglacial history of amphi‐Atlantic bryophytes

Paleontological evidence and current patterns of angiosperm species richness suggest that European biota experienced more severe bottlenecks than North American ones during the last glacial maximum. How well this pattern fits other plant species is less clear. Bryophytes offer a unique opportunity to contrast the impact of the last glacial maximum in North America and Europe because about 60% of the European bryoflora is shared with North America. Here, we use population genetic analyses based on approximate Bayesian computation on eight amphi‐Atlantic species to test the hypothesis that North American populations were less impacted by the last glacial maximum, exhibiting higher levels of genetic diversity than European ones and ultimately serving as a refugium for the postglacial recolonization of Europe. In contrast with this hypothesis, the best‐fit demographic model involved similar patterns of population size contractions, comparable levels of genetic diversity and balanced migration rates between European and North American populations. Our results thus suggest that bryophytes have experienced comparable demographic glacial histories on both sides of the Atlantic. Although a weak, but significant genetic structure was systematically recovered between European and North American populations, evidence for migration from and towards both continents suggests that amphi‐Atlantic bryophyte population may function as a metapopulation network. Reconstructing the biogeographic history of either North American or European bryophyte populations therefore requires a large, trans‐Atlantic geographic framework.

Targeted Long-Read Sequencing of a Locus Under Long-Term Balancing Selection in Capsella

Rapid advances in short-read DNA sequencing technologies have revolutionized population genomic studies, but there are genomic regions where this technology reaches its limits. Limitations mostly arise due to the difficulties in assembly or alignment to genomic regions of high sequence divergence and high repeat content, which are typical characteristics for loci under strong long-term balancing selection. Studying genetic diversity at such loci therefore remains challenging. Here, we investigate the feasibility and error rates associated with targeted long-read sequencing of a locus under balancing selection. For this purpose, we generated bacterial artificial chromosomes (BACs) containing the Brassicaceae S-locus, a region under strong negative frequency-dependent selection which has previously proven difficult to assemble in its entirety using short reads. We sequence S-locus BACs with single-molecule long-read sequencing technology and conduct de novo assembly of these S-locus haplotypes. By comparing repeated assemblies resulting from independent long-read sequencing runs on the same BAC clone we do not detect any structural errors, suggesting that reliable assemblies are generated, but we estimate an indel error rate of 5.7×10−5. A similar error rate was estimated based on comparison of Illumina short-read sequences and BAC assemblies. Our results show that, until de novo assembly of multiple individuals using long-read sequencing becomes feasible, targeted long-read sequencing of loci under balancing selection is a viable option with low error rates for single nucleotide polymorphisms or structural variation. We further find that short-read sequencing is a valuable complement, allowing correction of the relatively high rate of indel errors that result from this approach.

Integrative taxonomy reveals too extensive lumping and a new species in the moss genus Amphidium (Bryophyta)

An integrative taxonomic approach, including molecular phylogenetic reconstructions based on plastid rps4-trnF and nuclear ITS sequences, statistical analysis of morphological-anatomical characters, and classical taxonomy, indicates that the reduction of 13 Amphidium species to three in a recent morphological revision represents a case of too extensive lumping. Instead, six Amphidium species can be distinguished based on molecular and morphological data, the widespread Amphidium lapponicum, A. mougeotii, and A. tortuosum, as well as the Macaronesian endemic A. curvipes, the North American endemic A. californicum, and a newly discovered species from Central Asia (southern Siberia and northern Mongolia), A. asiaticum sp. nov. Diagnostic morphological characters for all six species are discussed. The present data confirm that species diversity of Amphidium is highest in the Holarctic, where all six species occur.

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

Shifts from outcrossing to self-fertilisation have occurred repeatedly in many different lineages of flowering plants, and often involve the breakdown of genetic outcrossing mechanisms. In the Brassicaceae, self-incompatibility (SI) allows plants to ensure outcrossing by recognition and rejection of self-pollen on the stigma. This occurs through the interaction of female and male specificity components, consisting of a pistil based receptor and a pollen-coat protein, both of which are encoded by tightly linked genes at the S-locus. When benefits of selfing are higher than costs of inbreeding, theory predicts that loss-of-function mutations in the male (pollen) SI component should be favoured, especially if they are dominant. However, it remains unclear whether mutations in the male component of SI are predominantly responsible for shifts to self-compatibility, and testing this prediction has been difficult due to the challenges of sequencing the highly polymorphic and repetitive ~100 kbp S-locus. The crucifer genus Capsella offers an excellent opportunity to study multiple transitions from outcrossing to self-fertilization, but so far, little is known about the genetic basis and timing of loss of SI in the self-fertilizing diploid Capsella orientalis. Here, we show that loss of SI in C. orientalis occurred within the past 2.6 Mya and maps as a dominant trait to the S-locus. Using targeted long-read sequencing of multiple complete S-haplotypes, we identify a frameshift deletion in the male specificity gene SCR that is fixed in C. orientalis, and we confirm loss of male SI specificity. We further analyze RNA sequencing data to identify a conserved, S-linked small RNA (sRNA) that is predicted to cause dominance of self-compatibility. Our results suggest that degeneration of pollen SI specificity in dominant S-alleles is important for shifts to self-fertilization in the Brassicaceae. Author Summary Already Darwin was fascinated by the widely varying modes of plant reproduction. The shift from outcrossing to self-fertilization is considered one of the most frequent evolutionary transitions in flowering plants, yet we still know little about the genetic basis of these shifts. In the Brassicaceae, outcrossing is enforced by a self-incompatibility (SI) system that enables the recognition and rejection of self pollen. This occurs through the action of two tightly linked genes at the S-locus, that encode a receptor protein located on the stigma (female component) and a pollen ligand protein (male component), respectively. Nevertheless, SI has frequently been lost, and theory predicts that mutations in the male component should have an advantage during the loss of SI, especially if they are dominant. To test this hypothesis, we mapped the loss of SI in a selfing species from the genus Capsella, a model system for evolutionary genomics. We found that loss of SI mapped to the S-locus, which harbored a dominant loss-of-function mutation in the male SI protein, and as expected, we found that male specificity was indeed lost in C. orientalis. Our results suggest that transitions to selfing often involve parallel genetic changes.

Pervasive population genomic consequences of genome duplication in Arabidopsis arenosa

Ploidy-variable species allow direct inference of the effects of chromosome copy number on fundamental evolutionary processes. While an abundance of theoretical work suggests polyploidy should leave distinct population genomic signatures, empirical data remains sparse. We sequenced ∼300 individuals from 39 populations of Arabidopsis arenosa, a naturally diploid-autotetraploid species. We find the impacts of polyploidy on population genomic processes are subtle yet pervasive, including reduced efficiency on linked and purifying selection as well as rampant gene flow from diploids. Initial masking of deleterious mutations, faster rates of nucleotide substitution, and interploidy introgression all conspire to shape the evolutionary potential of polyploids.