The rearing environment persistently modulates mouse phenotypes from the molecular to the behavioural level.

The phenotype of an organism results from its genotype and the influence of the environment throughout development. Even when using animals of the same genotype, independent studies may test animals of different phenotypes, resulting in poor replicability due to genotype-by-environment interactions. Thus, genetically defined strains of mice may respond differently to experimental treatments depending on their rearing environment. However, the extent of such phenotypic plasticity and its implications for the replicability of research findings have remained unknown. Here, we examined the extent to which common environmental differences between animal facilities modulate the phenotype of genetically homogeneous (inbred) mice. We conducted a comprehensive multicentre study, whereby inbred C57BL/6J mice from a single breeding cohort were allocated to and reared in 5 different animal facilities throughout early life and adolescence, before being transported to a single test laboratory. We found persistent effects of the rearing facility on the composition and heterogeneity of the gut microbial community. These effects were paralleled by persistent differences in body weight and in the behavioural phenotype of the mice. Furthermore, we show that environmental variation among animal facilities is strong enough to influence epigenetic patterns in neurons at the level of chromatin organisation. We detected changes in chromatin organisation in the regulatory regions of genes involved in nucleosome assembly, neuronal differentiation, synaptic plasticity, and regulation of behaviour. Our findings demonstrate that common environmental differences between animal facilities may produce facility-specific phenotypes, from the molecular to the behavioural level. Furthermore, they highlight an important limitation of inferences from single-laboratory studies and thus argue that study designs should take environmental background into account to increase the robustness and replicability of findings.

BAZ2A safeguards genome architecture of ground-state pluripotent stem cells.

Chromosomes have an intrinsic tendency to segregate into compartments, forming long-distance contacts between loci of similar chromatin states. How genome compartmentalization is regulated remains elusive. Here, comparison of mouse ground-state embryonic stem cells (ESCs) characterized by open and active chromatin, and advanced serum ESCs with a more closed and repressed genome, reveals distinct regulation of their genome organization due to differential dependency on BAZ2A/TIP5, a component of the chromatin remodeling complex NoRC. On ESC chromatin, BAZ2A interacts with SNF2H, DNA topoisomerase 2A (TOP2A) and cohesin. BAZ2A associates with chromatin sub-domains within the active A compartment, which intersect through long-range contacts. We found that ground-state chromatin selectively requires BAZ2A to limit the invasion of active domains into repressive compartments. BAZ2A depletion increases chromatin accessibility at B compartments. Furthermore, BAZ2A regulates H3K27me3 genome occupancy in a TOP2A-dependent manner. Finally, ground-state ESCs require BAZ2A for growth, differentiation, and correct expression of developmental genes. Our results uncover the propensity of open chromatin domains to invade repressive domains, which is counteracted by chromatin remodeling to establish genome partitioning and preserve cell identity.

Dissection of Besnoitia besnoiti intermediate host life cycle stages: From morphology to gene expression.

Cyst-forming Apicomplexa (CFA) of the Sarcocystidae have a ubiquitous presence as pathogens of humans and farm animals transmitted through the food chain between hosts with few notable exceptions. The defining hallmark of this family of obligate intracellular protists consists of their ability to remain for very long periods as infectious tissue cysts in chronically infected intermediate hosts. Nevertheless, each closely related species has evolved unique strategies to maintain distinct reservoirs on global scales and ensuring efficient transmission to definitive hosts as well as between intermediate hosts. Here, we present an in-depth comparative mRNA expression analysis of the tachyzoite and bradyzoite stages of Besnoitia besnoiti strain Lisbon14 isolated from an infected farm animal based on its annotated genome sequence. The B. besnoiti genome is highly syntenic with that of other CFA and also retains the capacity to encode a large majority of known and inferred factors essential for completing a sexual cycle in a yet unknown definitive host. This work introduces Besnoitia besnoiti as a new model for comparative biology of coccidian tissue cysts which can be readily obtained in high purity. This model provides a framework for addressing fundamental questions about the evolution of tissue cysts and the biology of this pharmacologically intractable infectious parasite stage.

Ecosystem type drives soil eukaryotic diversity and composition in Europe.

Soil eukaryotes play a crucial role in maintaining ecosystem functions and services, yet the factors driving their diversity and distribution remain poorly understood. While many studies focus on some eukaryotic groups (mostly fungi), they are limited in their spatial scale. Here, we analyzed an unprecedented amount of observational data of soil eukaryomes at continental scale (787 sites across Europe) to gain further insights into the impact of a wide range of environmental conditions (climatic and edaphic) on their community composition and structure. We found that the diversity of fungi, protists, rotifers, tardigrades, nematodes, arthropods, and annelids was predominantly shaped by ecosystem type (annual and permanent croplands, managed and unmanaged grasslands, coniferous and broadleaved woodlands), and higher diversity of fungi, protists, nematodes, arthropods, and annelids was observed in croplands than in less intensively managed systems, such as coniferous and broadleaved woodlands. Also in croplands, we found more specialized eukaryotes, while the composition between croplands was more homogeneous compared to the composition of other ecosystems. The observed high proportion of overlapping taxa between ecosystems also indicates that DNA has accumulated from previous land uses, hence mimicking the land transformations occurring in Europe in the last decades. This strong ecosystem-type influence was linked to soil properties, and particularly, soil pH was driving the richness of fungi, rotifers, and annelids, while plant-available phosphorus drove the richness of protists, tardigrades, and nematodes. Furthermore, the soil organic carbon to total nitrogen ratio crucially explained the richness of fungi, protists, nematodes, and arthropods, possibly linked to decades of agricultural inputs. Our results highlighted the importance of long-term environmental variables rather than variables measured at the time of the sampling in shaping soil eukaryotic communities, which reinforces the need to include those variables in addition to ecosystem type in future monitoring programs and conservation efforts.

Alterations of Stress-Related Glucocorticoids and Endocannabinoids in Hair of Chronic Cocaine Users.

BACKGROUND: Previous research in animals and humans has demonstrated a potential role of stress regulatory systems, such as the hypothalamic-pituitary-adrenal (HPA) axis and the endocannabinoid (eCB) system, in the development of substance use disorders. We thus investigated alterations of HPA and eCB markers in individuals with chronic cocaine use disorder by using an advanced hair analysis technique. METHODS: We compared hair concentrations of glucocorticoids (cortisone, cortisol) and the eCBs 2-arachidonylglycerol, anandamide (AEA), oleoylethanolamide (OEA), and palmitoylethanolamide (PEA) between 48 recreational cocaine users (RCU), 25 dependent cocaine users (DCU), and 67 stimulant-naïve controls. Self-reported substance use and hair concentrations of substances were also assessed. RESULTS: Significantly higher concentrations of hair cortisone were found in RCU and DCU compared with controls. Hair concentrations of OEA and PEA were significantly lower in DCU compared with RCU and controls. Additionally, within cocaine users, elevated cocaine hair concentration was a significant predictor for increased glucocorticoid and decreased OEA hair levels. Moreover, higher 3,4-methyl​enedioxymethamphetamine hair concentration was correlated with elevated cortisone and AEA, OEA, and PEA levels in hair within cocaine users, whereas more self-reported cannabis use was associated with lower eCBs levels in hair across the total sample. CONCLUSION: Our findings support the hypothesis that the HPA axis and eCB system might be important regulators for substance use disorders. The mechanistic understanding of changes in glucocorticoid and eCB levels in future research might be a promising pharmacological target to reduce stress-induced craving and relapse specifically in cocaine use disorder.

Hi-C analysis in Arabidopsis identifies the KNOT, a structure with similarities to the flamenco locus of Drosophila.

Chromosomes are folded, spatially organized, and regulated by epigenetic marks. How chromosomal architecture is connected to the epigenome is not well understood. We show that chromosomal architecture of Arabidopsis is tightly linked to the epigenetic state. Furthermore, we show how physical constraints, such as nuclear size, correlate with the folding principles of chromatin. We also describe a nuclear structure, termed KNOT, in which genomic regions of all five Arabidopsis chromosomes interact. These KNOT ENGAGED ELEMENT (KEE) regions represent heterochromatic islands within euchromatin. Similar to PIWI-interacting RNA clusters, such as flamenco in Drosophila, KEEs represent preferred landing sites for transposable elements, which may be part of a transposon defense mechanism in the Arabidopsis nucleus.

Inheritance of DNA methylation differences in the mangrove Rhizophora mangle.

The capacity to respond to environmental challenges ultimately relies on phenotypic variation which manifests from complex interactions of genetic and nongenetic mechanisms through development. While we know something about genetic variation and structure of many species of conservation importance, we know very little about the nongenetic contributions to variation. Rhizophora mangle is a foundation species that occurs in coastal estuarine habitats throughout the neotropics where it provides critical ecosystem functions and is potentially threatened by anthropogenic environmental changes. Several studies have documented landscape-level patterns of genetic variation in this species, but we know virtually nothing about the inheritance of nongenetic variation. To assess one type of nongenetic variation, we examined the patterns of DNA sequence and DNA methylation in maternal plants and offspring from natural populations of R. mangle from the Gulf Coast of Florida. We used a reduced representation bisulfite sequencing approach (epi-genotyping by sequencing; epiGBS) to address the following questions: (a) What are the levels of genetic and epigenetic diversity in natural populations of R. mangle? (b) How are genetic and epigenetic variation structured within and among populations? (c) How faithfully is epigenetic variation inherited? We found low genetic diversity but high epigenetic diversity from natural populations of maternal plants in the field. In addition, a large portion (up to ~25%) of epigenetic differences among offspring grown in common garden was explained by maternal family. Therefore, epigenetic variation could be an important source of response to challenging environments in the genetically depauperate populations of this foundation species.

Abscisic acid is a substrate of the ABC transporter encoded by the durable wheat disease resistance gene Lr34.

The wheat Lr34res allele, coding for an ATP-binding cassette transporter, confers durable resistance against multiple fungal pathogens. The Lr34sus allele, differing from Lr34res by two critical nucleotide polymorphisms, is found in susceptible wheat cultivars. Lr34res is functionally transferrable as a transgene into all major cereals, including rice, barley, maize, and sorghum. Here, we used transcriptomics, physiology, genetics, and in vitro and in vivo transport assays to study the molecular function of Lr34. We report that Lr34res results in a constitutive induction of transcripts reminiscent of an abscisic acid (ABA)-regulated response in transgenic rice. Lr34-expressing rice was altered in biological processes that are controlled by this phytohormone, including dehydration tolerance, transpiration and seedling growth. In planta seedling and in vitro yeast accumulation assays revealed that both LR34res and LR34sus act as ABA transporters. However, whereas the LR34res protein was detected in planta the LR34sus version was not, suggesting a post-transcriptional regulatory mechanism. Our results identify ABA as a substrate of the LR34 ABC transporter. We conclude that LR34res-mediated ABA redistribution has a major effect on the transcriptional response and physiology of Lr34res-expressing plants and that ABA is a candidate molecule that contributes to Lr34res-mediated disease resistance.

Evidence for rapid evolution in a grassland biodiversity experiment.

In long-term grassland experiments, positive biodiversity effects on plant productivity commonly increase with time. Subsequent glasshouse experiments showed that these strengthened positive biodiversity effects persist not only in the local environment but also when plants are transferred into a common environment. Thus, we hypothesized that community diversity had acted as a selective agent, resulting in the emergence of plant monoculture and mixture types with differing genetic composition. To test our hypothesis, we grew offspring from plants that were grown for eleven years in monoculture or mixture environments in a biodiversity experiment (Jena Experiment) under controlled glasshouse conditions in monocultures or two-species mixtures. We used epiGBS, a genotyping-by-sequencing approach combined with bisulphite conversion, to provide integrative genetic and epigenetic (i.e., DNA methylation) data. We observed significant divergence in genetic and DNA methylation data according to selection history in three out of five perennial grassland species, namely Galium mollugo, Prunella vulgaris and Veronica chamaedrys, with DNA methylation differences mostly reflecting the genetic differences. In addition, current diversity levels in the glasshouse had weak effects on epigenetic variation. However, given the limited genome coverage of the reference-free bisulphite method epiGBS, it remains unclear how much of the differences in DNA methylation was independent of underlying genetic differences. Our results thus suggest that selection of genetic variants, and possibly epigenetic variants, caused the rapid emergence of monoculture and mixture types within plant species in the Jena Experiment.

Feedbacks of plant identity and diversity on the diversity and community composition of rhizosphere microbiomes from a long-term biodiversity experiment.

Soil microbes are known to be key drivers of several essential ecosystem processes such as nutrient cycling, plant productivity and the maintenance of plant species diversity. However, how plant species diversity and identity affect soil microbial diversity and community composition in the rhizosphere is largely unknown. We tested whether, over the course of 11 years, distinct soil bacterial communities developed under plant monocultures and mixtures, and if over this time frame plants with a monoculture or mixture history changed in the bacterial communities they associated with. For eight species, we grew offspring of plants that had been grown for 11 years in the same field monocultures or mixtures (plant history in monoculture vs. mixture) in pots inoculated with microbes extracted from the field monoculture and mixture soils attached to the roots of the host plants (soil legacy). After 5 months of growth in the glasshouse, we collected rhizosphere soil from each plant and used 16S rRNA gene sequencing to determine the community composition and diversity of the bacterial communities. Bacterial community structure in the plant rhizosphere was primarily determined by soil legacy and by plant species identity, but not by plant history. In seven of the eight plant species the number of individual operational taxonomic units with increased abundance was larger when inoculated with microbes from mixture soil. We conclude that plant species richness can affect below-ground community composition and diversity, feeding back to the assemblage of rhizosphere bacterial communities in newly establishing plants via the legacy in soil.

Transcriptome analysis of Plasmodium berghei during exo-erythrocytic development.

BACKGROUND: The complex life cycle of malaria parasites requires well-orchestrated stage specific gene expression. In the vertebrate host the parasites grow and multiply by schizogony in two different environments: within erythrocytes and within hepatocytes. Whereas erythrocytic parasites are well-studied in this respect, relatively little is known about the exo-erythrocytic stages. METHODS: In an attempt to fill this gap, genome wide RNA-seq analyses of various exo-erythrocytic stages of Plasmodium berghei including sporozoites, samples from a time-course of liver stage development and detached cells were performed. These latter contain infectious merozoites and represent the final step in exo-erythrocytic development. RESULTS: The analysis represents the complete transcriptome of the entire life cycle of P. berghei parasites with temporal detailed analysis of the liver stage allowing comparison of gene expression across the progression of the life cycle. These RNA-seq data from different developmental stages were used to cluster genes with similar expression profiles, in order to infer their functions. A comparison with published data from other parasite stages confirmed stage-specific gene expression and revealed numerous genes that are expressed differentially in blood and exo-erythrocytic stages. One of the most exo-erythrocytic stage-specific genes was PBANKA_1003900, which has previously been annotated as a "gametocyte specific protein". The promoter of this gene drove high GFP expression in exo-erythrocytic stages, confirming its expression profile seen by RNA-seq. CONCLUSIONS: The comparative analysis of the genome wide mRNA expression profiles of erythrocytic and different exo-erythrocytic stages could be used to improve the understanding of gene regulation in Plasmodium parasites and can be used to model exo-erythrocytic stage metabolic networks toward the identification of differences in metabolic processes during schizogony in erythrocytes and hepatocytes.

Plant germline formation: common concepts and developmental flexibility in sexual and asexual reproduction.

The life cycle of flowering plants alternates between two heteromorphic generations: a diploid sporophytic generation and a haploid gametophytic generation. During the development of the plant reproductive lineages - the germlines - typically, single sporophytic (somatic) cells in the flower become committed to undergo meiosis. The resulting spores subsequently develop into highly polarized and differentiated haploid gametophytes that harbour the gametes. Recent studies have provided insights into the genetic basis and regulatory programs underlying cell specification and the acquisition of reproductive fate during both sexual reproduction and asexual (apomictic) reproduction. As we review here, these recent advances emphasize the importance of transcriptional, translational and post-transcriptional regulation, and the role of epigenetic regulatory pathways and hormonal activity.

Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes.

Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.

Starch Turnover and Metabolism during Flower and Early Embryo Development.

The accumulation of starch within photosynthetic tissues and within dedicated storage organs has been characterized extensively in many species, and a function in buffering carbon availability or in fueling later growth phases, respectively, has been proposed. However, developmentally regulated starch turnover within heterotrophic tissues other than dedicated storage organs is poorly characterized, and its function is not well understood. Here, we report on the characterization of starch turnover during flower, early embryo, and silique development in Arabidopsis (Arabidopsis thaliana) using a combined clearing-staining technique on whole-mount tissue. Besides the two previously documented waves of transient starch accumulation in the stamen envelope, occurring during meiosis and pollen mitosis I, we identified a novel, third wave of starch amylogenesis/amylolysis during the last stages of stamen development. To gain insights into the underlying molecular mechanisms, we analyzed publicly available microarray data, which revealed a developmentally coordinated expression of carbohydrate transport and metabolism genes during these waves of transient starch accumulation. Based on this analysis, we characterized starch dynamics in mutants affecting hexose phosphate metabolism and translocation, and identified the Glc-6-phosphate/phosphate antiporter GPT1 as the putative translocator of Glc-6-phosphate for starch biosynthesis in reproductive tissues. Based on these results, we propose a model of starch synthesis within the pollen grain and discuss the nutrient transport route feeding the embryo within the developing seed.

Apomictic and sexual germline development differ with respect to cell cycle, transcriptional, hormonal and epigenetic regulation.

Seeds of flowering plants can be formed sexually or asexually through apomixis. Apomixis occurs in about 400 species and is of great interest for agriculture as it produces clonal offspring. It differs from sexual reproduction in three major aspects: (1) While the sexual megaspore mother cell (MMC) undergoes meiosis, the apomictic initial cell (AIC) omits or aborts meiosis (apomeiosis); (2) the unreduced egg cell of apomicts forms an embryo without fertilization (parthenogenesis); and (3) the formation of functional endosperm requires specific developmental adaptations. Currently, our knowledge about the gene regulatory programs underlying apomixis is scarce. We used the apomict Boechera gunnisoniana, a close relative of Arabidopsis thaliana, to investigate the transcriptional basis underlying apomeiosis and parthenogenesis. Here, we present the first comprehensive reference transcriptome for reproductive development in an apomict. To compare sexual and apomictic development at the cellular level, we used laser-assisted microdissection combined with microarray and RNA-Seq analyses. Conservation of enriched gene ontologies between the AIC and the MMC likely reflects functions of importance to germline initiation, illustrating the close developmental relationship of sexuality and apomixis. However, several regulatory pathways differ between sexual and apomictic germlines, including cell cycle control, hormonal pathways, epigenetic and transcriptional regulation. Enrichment of specific signal transduction pathways are a feature of the apomictic germline, as is spermidine metabolism, which is associated with somatic embryogenesis in various plants. Our study provides a comprehensive reference dataset for apomictic development and yields important new insights into the transcriptional basis underlying apomixis in relation to sexual reproduction.

Rcount: simple and flexible RNA-Seq read counting.

SUMMARY: Analysis of differential gene expression by RNA sequencing (RNA-Seq) is frequently done using feature counts, i.e. the number of reads mapping to a gene. However, commonly used count algorithms (e.g. HTSeq) do not address the problem of reads aligning with multiple locations in the genome (multireads) or reads aligning with positions where two or more genes overlap (ambiguous reads). Rcount specifically addresses these issues. Furthermore, Rcount allows the user to assign priorities to certain feature types (e.g. higher priority for protein-coding genes compared to rRNA-coding genes) or to add flanking regions. AVAILABILITY AND IMPLEMENTATION: Rcount provides a fast and easy-to-use graphical user interface requiring no command line or programming skills. It is implemented in C++ using the SeqAn (www.seqan.de) and the Qt libraries (qt-project.org). Source code and 64 bit binaries for (Ubuntu) Linux, Windows (7) and MacOSX are released under the GPLv3 license and are freely available on github.com/MWSchmid/Rcount. CONTACT: marcschmid@gmx.ch SUPPLEMENTARY INFORMATION: Test data, genome annotation files, useful Python and R scripts and a step-by-step user guide (including run-time and memory usage tests) are available on github.com/MWSchmid/Rcount.

HiCdat: a fast and easy-to-use Hi-C data analysis tool.

BACKGROUND: The study of nuclear architecture using Chromosome Conformation Capture (3C) technologies is a novel frontier in biology. With further reduction in sequencing costs, the potential of Hi-C in describing nuclear architecture as a phenotype is only about to unfold. To use Hi-C for phenotypic comparisons among different cell types, conditions, or genetic backgrounds, Hi-C data processing needs to be more accessible to biologists. RESULTS: HiCdat provides a simple graphical user interface for data pre-processing and a collection of higher-level data analysis tools implemented in R. Data pre-processing also supports a wide range of additional data types required for in-depth analysis of the Hi-C data (e.g. RNA-Seq, ChIP-Seq, and BS-Seq). CONCLUSIONS: HiCdat is easy-to-use and provides solutions starting from aligned reads up to in-depth analyses. Importantly, HiCdat is focussed on the analysis of larger structural features of chromosomes, their correlation to genomic and epigenomic features, and on comparative studies. It uses simple input and output formats and can therefore easily be integrated into existing workflows or combined with alternative tools.

Selection-driven evolution of sex-biased genes is consistent with sexual selection in Arabidopsis thaliana.

Sex-biased genes are genes with a preferential or specific expression in one sex and tend to show an accelerated rate of evolution in animals. Various hypotheses--which are not mutually exclusive--have been put forth to explain observed patterns of rapid evolution. One possible explanation is positive selection, but this has been shown only in few animal species and mostly for male-specific genes. Here, we present a large-scale study that investigates evolutionary patterns of sex-biased genes in the predominantly self-fertilizing plant Arabidopsis thaliana. Unlike most animal species, A. thaliana does not possess sex chromosomes, its flowers develop both male and female sexual organs, and it is characterized by low outcrossing rates. Using cell-specific gene expression data, we identified genes whose expression is enriched in comparison with all other tissues in the male and female gametes (sperm, egg, and central cell), as well as in synergids, pollen, and pollen tubes, which also play an important role in reproduction. Genes specifically expressed in gametes and synergids show higher rates of protein evolution compared with the genome-wide average and no evidence for positive selection. In contrast, pollen- and pollen tube-specific genes not only have lower rates of protein evolution but also exhibit a higher proportion of adaptive amino acid substitutions. We show that this is the result of increased levels of purifying and positive selection among genes with pollen- and pollen tube-specific expression. The increased proportion of adaptive substitutions cannot be explained by the fact that pollen- and pollen tube-expressed genes are enriched in segmental duplications, are on average older, or have a larger effective population size. Our observations are consistent with prezygotic sexual selection as a result of interactions during pollination and pollen tube growth such as pollen tube competition.

Evolutionary trajectory of pattern recognition receptors in plants.

Cell-surface receptors play pivotal roles in many biological processes, including immunity, development, and reproduction, across diverse organisms. How cell-surface receptors evolve to become specialised in different biological processes remains elusive. To shed light on the immune-specificity of cell-surface receptors, we analyzed more than 200,000 genes encoding cell-surface receptors from 350 genomes and traced the evolutionary origin of immune-specific leucine-rich repeat receptor-like proteins (LRR-RLPs) in plants. Surprisingly, we discovered that the motifs crucial for co-receptor interaction in LRR-RLPs are closely related to those of the LRR-receptor-like kinase (RLK) subgroup Xb, which perceives phytohormones and primarily governs growth and development. Functional characterisation further reveals that LRR-RLPs initiate immune responses through their juxtamembrane and transmembrane regions, while LRR-RLK-Xb members regulate development through their cytosolic kinase domains. Our data suggest that the cell-surface receptors involved in immunity and development share a common origin. After diversification, their ectodomains, juxtamembrane, transmembrane, and cytosolic regions have either diversified or stabilised to recognise diverse ligands and activate differential downstream responses. Our work reveals a mechanism by which plants evolve to perceive diverse signals to activate the appropriate responses in a rapidly changing environment.

Covering the bases: Population genomic structure of Lemna minor and the cryptic species L. japonica in Switzerland.

Duckweeds, including the common duckweed Lemna minor, are increasingly used to test eco-evolutionary theories. Yet, despite its popularity and near-global distribution, the understanding of its population structure (and genetic variation therein) is still limited. It is essential that this is resolved, because of the impact genetic diversity has on experimental responses and scientific understanding. Through whole-genome sequencing, we assessed the genetic diversity and population genomic structure of 23 natural Lemna spp. populations from their natural range in Switzerland. We used two distinct analytical approaches, a reference-free kmer approach and the classical reference-based one. Two genetic clusters were identified across the described species distribution of L. minor, surprisingly corresponding to species-level divisions. The first cluster contained the targeted L. minor individuals and the second contained individuals from a cryptic species: Lemna japonica. Within the L. minor cluster, we identified a well-defined population structure with little intra-population genetic diversity (i.e., within ponds) but high inter-population diversity (i.e., between ponds). In L. japonica, the population structure was significantly weaker and genetic variation between a subset of populations was as low as within populations. This study revealed that L. japonica is more widespread than previously thought. Our findings signify that thorough genotype-to-phenotype analyses are needed in duckweed experimental ecology and evolution.