Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

Interactions between U and V sex chromosomes during the life cycle of Ectocarpus.

In many animals and flowering plants, sex determination occurs in the diploid phase of the life cycle with XX/XY or ZW/ZZ sex chromosomes. However, in early diverging plants and most macroalgae, sex is determined by female (U) or male (V) sex chromosomes in a haploid phase called the gametophyte. Once the U and V chromosomes unite at fertilization to produce a diploid sporophyte, sex determination no longer occurs, raising key questions about the fate of the U and V sex chromosomes in the sporophyte phase. Here, we investigate genetic and molecular interactions of the UV sex chromosomes in both the haploid and diploid phases of the brown alga Ectocarpus. We reveal extensive developmental regulation of sex chromosome genes across its life cycle and implicate the TALE-HD transcription factor OUROBOROS in suppressing sex determination in the diploid phase. Small RNAs may also play a role in the repression of a female sex-linked gene, and transition to the diploid sporophyte coincides with major reconfiguration of histone H3K79me2, suggesting a more intricate role for this histone mark in Ectocarpus development than previously appreciated.

JGI Plant Gene Atlas: an updateable transcriptome resource to improve functional gene descriptions across the plant kingdom.

Gene functional descriptions offer a crucial line of evidence for candidate genes underlying trait variation. Conversely, plant responses to environmental cues represent important resources to decipher gene function and subsequently provide molecular targets for plant improvement through gene editing. However, biological roles of large proportions of genes across the plant phylogeny are poorly annotated. Here we describe the Joint Genome Institute (JGI) Plant Gene Atlas, an updateable data resource consisting of transcript abundance assays spanning 18 diverse species. To integrate across these diverse genotypes, we analyzed expression profiles, built gene clusters that exhibited tissue/condition specific expression, and tested for transcriptional response to environmental queues. We discovered extensive phylogenetically constrained and condition-specific expression profiles for genes without any previously documented functional annotation. Such conserved expression patterns and tightly co-expressed gene clusters let us assign expression derived additional biological information to 64 495 genes with otherwise unknown functions. The ever-expanding Gene Atlas resource is available at JGI Plant Gene Atlas (https://plantgeneatlas.jgi.doe.gov) and Phytozome (https://phytozome.jgi.doe.gov/), providing bulk access to data and user-specified queries of gene sets. Combined, these web interfaces let users access differentially expressed genes, track orthologs across the Gene Atlas plants, graphically represent co-expressed genes, and visualize gene ontology and pathway enrichments.

Co-action of COP1, SPA and cryptochrome in light signal transduction and photomorphogenesis of the moss Physcomitrium patens.

The Arabidopsis COP1/SPA ubiquitin ligase suppresses photomorphogenesis in darkness. In the light, photoreceptors inactivate COP1/SPA to allow a light response. While SPA genes are specific to the green lineage, COP1 also exists in humans. This raises the question of when in evolution plant COP1 acquired the need for SPA accessory proteins. We addressed this question by generating Physcomitrium Ppcop1 mutants and comparing their visible and molecular phenotypes with those of Physcomitrium Ppspa mutants. The phenotype of Ppcop1 nonuple mutants resembles that of Ppspa mutants. Most importantly, both mutants produce green chloroplasts in complete darkness. They also exhibit dwarfed gametophores, disturbed branching of protonemata and absent gravitropism. RNA-sequencing analysis indicates that both mutants undergo weak constitutive light signaling in darkness. PpCOP1 and PpSPA proteins form a complex and they interact via their WD repeat domains with the VP motif of the cryptochrome CCE domain in a blue light-dependent manner. This resembles the interaction of Arabidopsis SPA proteins with Arabidopsis CRY1, and is different from that with Arabidopsis CRY2. Taken together, the data indicate that PpCOP1 and PpSPA act together to regulate growth and development of Physcomitrium. However, in contrast to their Arabidopsis orthologs, PpCOP1 and PpSPA proteins execute only partial suppression of light signaling in darkness. Hence, additional repressors may exist that contribute to the repression of a light response in dark-exposed Physcomitrium.

Evolutionary Conserved and Divergent Responses to Copper Zinc Superoxide Dismutase Inhibition in Plants.

After an initial evolution in a reducing environment, life got successively challenged by reactive oxygen species (ROS), especially during the great oxidation event (GOE) that followed the development of photosynthesis. Therefore, ROS are deeply intertwined into the physiological, morphological and transcriptional responses of most present-day organisms. Copper-zinc superoxide dismutases (CuZnSODs) evolved during the GOE and are present in charophytes and extant land plants, but nearly absent from chlorophytes. The chemical inhibitor of CuZnSOD, lung cancer screen 1 (LCS-1), could greatly facilitate the study of SODs in diverse plants. Here, we determined the impact of chemical inhibition of plant CuZnSOD activity, on plant growth, transcription and metabolism. We followed a comparative approach by using different plant species, including Marchantia Polymorpha and Physcomitrium patens, representing bryophytes, the sister lineage to vascular plants, and Arabidopsis thaliana. We show that LCS-1 causes oxidative stress in plants and that the inhibition of CuZnSODs provoked a similar core response that mainly impacted glutathione homoeostasis in all plant species analysed. That said, Physcomitrium and Arabidopsis, which contain multiple CuZnSOD isoforms showed a more complex and exacerbated response. In addition, an untargeted metabolomics approach revealed a specific metabolic signature for each plant species. Our comparative analysis exposes a conserved core response at the physiological and transcriptional level towards LCS-1, while the metabolic response largely varies. These differences correlate with the number and localization of the CuZnSOD isoforms present in each species.

A vertically transmitted amalgavirus is present in certain accessions of the bryophyte Physcomitrium patens.

In the last few years, next-generation sequencing techniques have started to be used to identify new viruses infecting plants. This has allowed to rapidly increase our knowledge on viruses other than those causing symptoms in economically important crops. Here we used this approach to identify a virus infecting Physcomitrium patens that has the typical structure of the double-stranded RNA endogenous viruses of the Amalgaviridae family, which we named Physcomitrium patens amalgavirus 1, or PHPAV1. PHPAV1 is present only in certain accessions of P. patens, where its RNA can be detected throughout the cell cycle of the plant. Our analysis demonstrates that PHPAV1 can be vertically transmitted through both paternal and maternal germlines, in crosses between accessions that contain the virus with accessions that do not contain it. This work suggests that PHPAV1 can replicate in genomic backgrounds different from those that actually contain the virus and opens the door for future studies on virus-host coevolution.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes.

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens.

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Podocyte endowment and the impact of adult body size on kidney health.

Low birth weight is a risk factor for chronic kidney disease, whereas adult podocyte depletion is a key event in the pathogenesis of glomerulosclerosis. However, whether low birth weight due to poor maternal nutrition is associated with low podocyte endowment and glomerulosclerosis in later life is not known. Female Sprague-Dawley rats were fed a normal-protein diet (NPD; 20%) or low-protein diet (LPD; 8%), to induce low birth weight, from 3 wk before mating until postnatal day 21 (PN21), when kidneys from some male offspring were taken for quantitation of podocyte number and density in whole glomeruli using immunolabeling, tissue clearing, and confocal microscopy. The remaining offspring were fed a normal- or high-fat diet until 6 mo to induce catch-up growth and excessive weight gain, respectively. At PN21, podocyte number per glomerulus was 15% lower in low birth weight (LPD) than normal birth weight (NPD) offspring, with this deficit greater in outer glomeruli. Surprisingly, podocyte number in LPD offspring increased in outer glomeruli between PN21 and 6 mo, although an overall 9% podocyte deficit persisted. Postnatal fat feeding to LPD offspring did not alter podometric indexes or result in glomerular pathology at 6 mo, whereas fat feeding in NPD offspring was associated with far greater body and fat mass as well as podocyte loss, reduced podocyte density, albuminuria, and glomerulosclerosis. This is the first report that maternal diet can influence podocyte endowment. Our findings provide new insights into the impact of low birth weight, podocyte endowment, and postnatal weight on podometrics and kidney health in adulthood.NEW & NOTEWORTHY The present study shows, for the first time, that low birth weight as a result of maternal nutrition is associated with low podocyte endowment. However, a mild podocyte deficit at birth did not result in glomerular pathology in adulthood. In contrast, postnatal podocyte loss in combination with excessive body weight led to albuminuria and glomerulosclerosis. Taken together, these findings provide new insights into the associations between birth weight, podocyte indexes, postnatal weight, and glomerular pathology.

Coregulation of gene expression by White collar 1 and phytochrome in Ustilago maydis.

Ustilago maydis encodes ten predicted light-sensing proteins. The biological functions of only a few of them are elucidated. Among the characterized ones are two DNA-photolyases and two rhodopsins that act as DNA-repair enzymes or green light-driven proton pumps, respectively. Here we report on the role of two other photoreceptors in U. maydis, namely White collar 1 (Wco1) and Phytochrome 1 (Phy1). We show that they bind flavins or biliverdin as chromophores, respectively. Both photoreceptors undergo a photocycle in vitro. Wco1 is the dominant blue light receptor in the saprophytic phase, controlling all of the 324 differentially expressed genes in blue light. U. maydis also responds to red and far-red light. However, the number of red or far-red light-controlled genes is less compared to blue light-regulated ones. Moreover, most of the red and far-red light-controlled genes not only depend on Phy1 but also on Wco1, indicating partial coregulation of gene expression by both photoreceptors. GFP-fused Wco1 is preferentially located in the nucleus, Phy1 in the cytosol, thus providing no hint that these photoreceptors directly interact or operate within the same complex. This is the first report on a functional characterization and coaction of White collar 1 and phytochrome orthologs in basidiomycetes.

The Biotrophic Development of Ustilago maydis Studied by RNA-Seq Analysis.

The maize smut fungus Ustilago maydis is a model organism for elucidating host colonization strategies of biotrophic fungi. Here, we performed an in depth transcriptional profiling of the entire plant-associated development of U. maydis wild-type strains. In our analysis, we focused on fungal metabolism, nutritional strategies, secreted effectors, and regulatory networks. Secreted proteins were enriched in three distinct expression modules corresponding to stages on the plant surface, establishment of biotrophy, and induction of tumors. These modules are likely the key determinants for U. maydis virulence. With respect to nutrient utilization, we observed that expression of several nutrient transporters was tied to these virulence modules rather than being controlled by nutrient availability. We show that oligopeptide transporters likely involved in nitrogen assimilation are important virulence factors. By measuring the intramodular connectivity of transcription factors, we identified the potential drivers for the virulence modules. While known components of the b-mating type cascade emerged as inducers for the plant surface and biotrophy module, we identified a set of yet uncharacterized transcription factors as likely responsible for expression of the tumor module. We demonstrate a crucial role for leaf tumor formation and effector gene expression for one of these transcription factors.

Characterisation of evolutionarily conserved key players affecting eukaryotic flagellar motility and fertility using a moss model.

Defects in flagella/cilia are often associated with infertility and disease. Motile male gametes (sperm cells) are an ancestral eukaryotic trait that has been lost in several lineages like flowering plants. Here, we made use of a phenotypic male fertility difference between two moss (Physcomitrella patens) ecotypes to explore spermatozoid function. We compare genetic and epigenetic variation as well as expression profiles between the Gransden and Reute ecotype to identify a set of candidate genes associated with moss male infertility. We generated a loss-of-function mutant of a coiled-coil domain containing 39 (ccdc39) gene that is part of the flagellar hydin network. Defects in mammal and algal homologues of this gene coincide with a loss of fertility, demonstrating the evolutionary conservation of flagellar function related to male fertility across kingdoms. The Ppccdc39 mutant resembles the Gransden phenotype in terms of male fertility. Potentially, several somatic (epi-)mutations occurred during prolonged vegetative propagation of Gransden, causing regulatory differences of for example the homeodomain transcription factor BELL1. Probably these somatic changes are causative for the observed male fertility defect. We propose that moss spermatozoids might be employed as an easily accessible system to study male infertility of humans and animals in terms of flagellar structure and movement.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data.

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Sexual reproduction, sporophyte development and molecular variation in the model moss Physcomitrella patens: introducing the ecotype Reute.

Rich ecotype collections are used for several plant models to unravel the molecular causes of phenotypic differences, and to investigate the effects of environmental adaption and acclimation. For the model moss Physcomitrella patens collections of accessions are available, and have been used for phylogenetic and taxonomic studies, for example, but few have been investigated further for phenotypic differences. Here, we focus on the Reute accession and provide expression profiling and comparative developmental data for several stages of sporophyte development, as well as information on genetic variation via genomic sequencing. We analysed cross-technology and cross-laboratory data to define a confident set of 15 mature sporophyte-specific genes. We find that the standard laboratory strain Gransden produces fewer sporophytes than Reute or Villersexel, although gametangia develop with the same time course and do not show evident morphological differences. Reute exhibits less genetic variation relative to Gransden than Villersexel, yet we found variation between Gransden and Reute in the expression profiles of several genes, as well as variation hot spots and genes that appear to evolve under positive Darwinian selection. We analyzed expression differences between the ecotypes for selected candidate genes in the GRAS transcription factor family, the chalcone synthase family and in genes involved in cell wall modification that are potentially related to phenotypic differences. We confirm that Reute is a P. patens ecotype, and suggest its use for reverse-genetics studies that involve progression through the life cycle and multiple generations.

Electron Interference in Hall Effect Measurements on GaAs/InAs Core/Shell Nanowires.

We present low-temperature magnetotransport measurements on GaAs/InAs core/shell nanowires contacted by regular source-drain leads as well as laterally attached Hall contacts, which only touch parts of the nanowire sidewalls. Low-temperature measurements between source and drain contacts show typical phase coherent effects, such as universal conductance fluctuations in a magnetic field aligned perpendicularly to the nanowire axis as well as Aharonov-Bohm-type oscillations in a parallel aligned magnetic field. However, the signal between the Hall contacts shows a Hall voltage buildup, when the magnetic field is turned perpendicular to the nanowire axis while current is driven through the wire using the source-drain contacts. At low temperatures, the phase coherent effects measured between source and drain leads are superimposed on the Hall voltage, which can be explained by nonlocal probing of large segments of the nanowire. In addition, the Aharonov-Bohm-type oscillations are also observed in the magnetoconductance at magnetic fields aligned parallel to the nanowire axis, using the laterally contacted leads. This measurement geometry hereby directly corresponds to classical Aharonov-Bohm experiments using planar quantum rings. In addition, the Hall voltage is used to characterize the nanowires in terms of charge carrier concentration and mobility, using temperature- and gate-dependent measurements as well as measurements in tilted magnetic fields. The GaAs/InAs core/shell nanowire used in combination with laterally attached contacts is therefore the ideal system to three-dimensionally combine quantum ring experiments using the cross-sectional plane and Hall experiments using the axial nanowire plane.

Anisotropic phase coherence in GaAs/InAs core/shell nanowires.

Low-temperature transport in nanowires is accompanied by phase-coherent effects, which are observed as modulation of the conductance in an external magnetic field. In the GaAs/InAs core/shell nanowires investigated here, these are h/e flux periodic oscillations in a magnetic field aligned parallel to the nanowire axis and aperiodic universal conductance fluctuations in a field aligned perpendicularly to the nanowire axis. Both electron interference effects are used to analyse the phase coherence of the system. Temperature-dependent measurements are carried out, in order to derive the phase coherence lengths in the cross-sectional plane as well as along the nanowire sidewalls. It is found that these values show a strong anisotropy, which can be explained by the crystal structure of the GaAs/InAs core/shell nanowire. For nanowires with a radius as low as 45 nm, flux periodic oscillations were observed up to a temperature of 55 K.

Dynamics of transcriptional start site selection during nitrogen stress-induced cell differentiation in Anabaena sp. PCC7120.

The fixation of atmospheric N(2) by cyanobacteria is a major source of nitrogen in the biosphere. In Nostocales, such as Anabaena, this process is spatially separated from oxygenic photosynthesis and occurs in heterocysts. Upon nitrogen step-down, these specialized cells differentiate from vegetative cells in a process controlled by two major regulators: NtcA and HetR. However, the regulon controlled by these two factors is only partially defined, and several aspects of the differentiation process have remained enigmatic. Using differential RNA-seq, we experimentally define a genome-wide map of >10,000 transcriptional start sites (TSS) of Anabaena sp. PCC7120, a model organism for the study of prokaryotic cell differentiation and N(2) fixation. By analyzing the adaptation to nitrogen stress, our global TSS map provides insight into the dynamic changes that modify the transcriptional organization at a critical step of the differentiation process. We identify >900 TSS with minimum fold change in response to nitrogen deficiency of eight. From these TSS, at least 209 were under control of HetR, whereas at least 158 other TSS were potentially directly controlled by NtcA. Our analysis of the promoters activated during the switch to N(2) fixation adds hundreds of protein-coding genes and noncoding transcripts to the list of potentially involved factors. These data experimentally define the NtcA regulon and the DIF(+) motif, a palindrome at or close to position -35 that seems essential for heterocyst-specific expression of certain genes.