A novel TFL1 gene induces flowering in the mast seeding alpine snow tussock, Chionochloa pallens (Poaceae).

Masting, the synchronous, highly variable flowering across years by a population of perennial plants, has been reported to be precipitated by various factors including nitrogen levels, drought conditions, and spring and summer temperatures. However, the molecular mechanism leading to the initiation of flowering in masting plants in particular years remains largely unknown, despite the potential impact of climate change on masting phenology. We studied genes controlling flowering in the alpine snow tussock Chionochloa pallens (Poaceae), a strongly masting perennial grass. We used a range of in situ and manipulated plants to obtain leaf samples from tillers (shoots) which subsequently remained vegetative or flowered. Here, we show that a novel orthologue of TERMINAL FLOWER 1 (TFL1; normally a repressor of flowering in other species) promotes the induction of flowering in C. pallens (hence Anti-TFL1), a conclusion supported by structural, functional and expression analyses. Global transcriptomic analysis indicated differential expression of CpTPS1, CpGA20ox1, CpREF6 and CpHDA6, emphasizing the role of endogenous cues and epigenetic regulation in terms of responsiveness of plants to initiate flowering. Our molecular-based study provides insights into the cellular mechanism of flowering in masting plants and will supplement ecological and statistical models to predict how masting will respond to global climate change.

Structural Phylogenetics with Confidence.

For evaluating the deepest evolutionary relationships among proteins, sequence similarity is too low for application of sequence-based homology search or phylogenetic methods. In such cases, comparison of protein structures, which are often better conserved than sequences, may provide an alternative means of uncovering deep evolutionary signal. Although major protein structure databases such as SCOP and CATH hierarchically group protein structures, they do not describe the specific evolutionary relationships within a hierarchical level. Structural phylogenies have the potential to fill this gap. However, it is difficult to assess evolutionary relationships derived from structural phylogenies without some means of assessing confidence in such trees. We therefore address two shortcomings in the application of structural data to deep phylogeny. First, we examine whether phylogenies derived from pairwise structural comparisons are sensitive to differences in protein length and shape. We find that structural phylogenetics is best employed where structures have very similar lengths, and that shape fluctuations generated during molecular dynamics simulations impact pairwise comparisons, but not so drastically as to eliminate evolutionary signal. Second, we address the absence of statistical support for structural phylogeny. We present a method for assessing confidence in a structural phylogeny using shape fluctuations generated via molecular dynamics or Monte Carlo simulations of proteins. Our approach will aid the evolutionary reconstruction of relationships across structurally defined protein superfamilies. With the Protein Data Bank now containing in excess of 158,000 entries (December 2019), we predict that structural phylogenetics will become a useful tool for ordering the protein universe.

Molecular control of the floral transition in the mast seeding plant Celmisia lyallii (Asteraceae).

Mast flowering (or masting) is synchronous, highly variable flowering among years in populations of perennial plants. Despite having widespread consequences for seed consumers, endangered fauna and human health, masting is hard to predict. While observational studies show links to various weather patterns in different plant species, the mechanism(s) underpinning the regulation of masting is still not fully explained. We studied floral induction in Celmisia lyallii (Asteraceae), a mast flowering herbaceous alpine perennial, comparing gene expression in flowering and nonflowering plants. We performed translocation experiments to induce the floral transition in C. lyallii plants followed by both global and targeted expression analysis of flowering-pathway genes. Differential expression analysis showed elevated expression of ClSOC1 and ClmiR172 (promoters of flowering) in leaves of plants that subsequently flowered, in contrast to elevated expression of ClAFT and ClTOE1 (repressors of flowering) in leaves of plants that did not flower. The warm summer conditions that promoted flowering led to differential regulation of age and hormonal pathway genes, including ClmiR172 and ClGA20ox2, known to repress the expression of floral repressors and permit flowering. Upregulated expression of epigenetic modifiers of floral promoters also suggests that plants may maintain a novel "summer memory" across years to induce flowering. These results provide a basic mechanistic understanding of floral induction in masting plants and evidence of their ability to imprint various environmental cues to synchronize flowering, allowing us to better predict masting events under climate change.

An Evaluation of Function of Multicopy Noncoding RNAs in Mammals Using ENCODE/FANTOM Data and Comparative Genomics.

Mammalian diversification has coincided with a rapid proliferation of various types of noncoding RNAs, including members of both snRNAs and snoRNAs. The significance of this expansion however remains obscure. While some ncRNA copy-number expansions have been linked to functionally tractable effects, such events may equally likely be neutral, perhaps as a result of random retrotransposition. Hindering progress in our understanding of such observations is the difficulty in establishing function for the diverse features that have been identified in our own genome. Projects such as ENCODE and FANTOM have revealed a hidden world of genomic expression patterns, as well as a host of other potential indicators of biological function. However, such projects have been criticized, particularly from practitioners in the field of molecular evolution, where many suspect these data provide limited insight into biological function. The molecular evolution community has largely taken a skeptical view, thus it is important to establish tests of function. We use a range of data, including data drawn from ENCODE and FANTOM, to examine the case for function for the recent copy number expansion in mammals of six evolutionarily ancient RNA families involved in splicing and rRNA maturation. We use several criteria to assess evidence for function: conservation of sequence and structure, genomic synteny, evidence for transposition, and evidence for species-specific expression. Applying these criteria, we find that only a minority of loci show strong evidence for function and that, for the majority, we cannot reject the null hypothesis of no function.

The fitness challenge of studying molecular adaptation.

Advances in bioinformatics and high-throughput genetic analysis increasingly allow us to predict the genetic basis of adaptive traits. These predictions can be tested and confirmed, but the molecular-level changes - i.e. the molecular adaptation - that link genetic differences to organism fitness remain generally unknown. In recent years, a series of studies have started to unpick the mechanisms of adaptation at the molecular level. In particular, this work has examined how changes in protein function, activity, and regulation cause improved organismal fitness. Key to addressing molecular adaptations is identifying systems and designing experiments that integrate changes in the genome, protein chemistry (molecular phenotype), and fitness. Knowledge of the molecular changes underpinning adaptations allow new insight into the constraints on, and repeatability of adaptations, and of the basis of non-additive interactions between adaptive mutations. Here we critically discuss a series of studies that examine the molecular-level adaptations that connect genetic changes and fitness.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Does the Ribosome Challenge our Understanding of the RNA World?

In a recent article published in these pages, Bowman and colleagues propose that the ribosome represents a challenge to the RNA world model, a long-standing framework to explain the origin of DNA and genetically encoded proteins from a hypothetical RNA-based system. Specifically, they outline a scenario for the emergence and subsequent coevolution of the peptidyl transferase centre (PTC) of the ribosome with non-templated peptide products of this RNA through chemical evolution. They also propose that the PTC would have predated the emergence of enzymatic RNA replication, and that this in turn indicates that the RNA world never existed. We and others have previously incorporated non-templated peptide production as an early stage in the evolution of protein synthesis, which we would count as a chemical process, in agreement with Bowman and colleagues' model. However, their model raises an important question: to what extent could early protein synthesis and its products have evolved in the absence of Darwinian processes? We argue that evolution of the early ribosome requires Darwinian evolution, and that, while chemical evolution could give rise to peptidyl transferase activity, it is insufficient for subsequent improvement of a proto-PTC, or for ongoing coevolution of the proto-PTC with its early non-templated peptide products. We conclude that it is difficult to preclude the involvement of replicative processes, themselves subject to Darwinian evolution, from the evolution of the PTC. Finally, Bowman et al. call into question current models for the RNA to protein transition. We show that the difficulty that Bowman et al. have with this scenario is down to a misreading of our previous work.

Bacterial natural transformation by highly fragmented and damaged DNA.

DNA molecules are continuously released through decomposition of organic matter and are ubiquitous in most environments. Such DNA becomes fragmented and damaged (often <100 bp) and may persist in the environment for more than half a million years. Fragmented DNA is recognized as nutrient source for microbes, but not as potential substrate for bacterial evolution. Here, we show that fragmented DNA molecules (≥ 20 bp) that additionally may contain abasic sites, cross-links, or miscoding lesions are acquired by the environmental bacterium Acinetobacter baylyi through natural transformation. With uptake of DNA from a 43,000-y-old woolly mammoth bone, we further demonstrate that such natural transformation events include ancient DNA molecules. We find that the DNA recombination is RecA recombinase independent and is directly linked to DNA replication. We show that the adjacent nucleotide variations generated by uptake of short DNA fragments escape mismatch repair. Moreover, double-nucleotide polymorphisms appear more common among genomes of transformable than nontransformable bacteria. Our findings reveal that short and damaged, including truly ancient, DNA molecules, which are present in large quantities in the environment, can be acquired by bacteria through natural transformation. Our findings open for the possibility that natural genetic exchange can occur with DNA up to several hundreds of thousands years old.

Robust identification of noncoding RNA from transcriptomes requires phylogenetically-informed sampling.

Noncoding RNAs are integral to a wide range of biological processes, including translation, gene regulation, host-pathogen interactions and environmental sensing. While genomics is now a mature field, our capacity to identify noncoding RNA elements in bacterial and archaeal genomes is hampered by the difficulty of de novo identification. The emergence of new technologies for characterizing transcriptome outputs, notably RNA-seq, are improving noncoding RNA identification and expression quantification. However, a major challenge is to robustly distinguish functional outputs from transcriptional noise. To establish whether annotation of existing transcriptome data has effectively captured all functional outputs, we analysed over 400 publicly available RNA-seq datasets spanning 37 different Archaea and Bacteria. Using comparative tools, we identify close to a thousand highly-expressed candidate noncoding RNAs. However, our analyses reveal that capacity to identify noncoding RNA outputs is strongly dependent on phylogenetic sampling. Surprisingly, and in stark contrast to protein-coding genes, the phylogenetic window for effective use of comparative methods is perversely narrow: aggregating public datasets only produced one phylogenetic cluster where these tools could be used to robustly separate unannotated noncoding RNAs from a null hypothesis of transcriptional noise. Our results show that for the full potential of transcriptomics data to be realized, a change in experimental design is paramount: effective transcriptomics requires phylogeny-aware sampling.

Evolutionary history of the TBP-domain superfamily.

The TATA binding protein (TBP) is an essential transcription initiation factor in Archaea and Eucarya. Bacteria lack TBP, and instead use sigma factors for transcription initiation. TBP has a symmetric structure comprising two repeated TBP domains. Using sequence, structural and phylogenetic analyses, we examine the distribution and evolutionary history of the TBP domain, a member of the helix-grip fold family. Our analyses reveal a broader distribution than for TBP, with TBP-domains being present across all three domains of life. In contrast to TBP, all other characterized examples of the TBP domain are present as single copies, primarily within multidomain proteins. The presence of the TBP domain in the ubiquitous DNA glycosylases suggests that this fold traces back to the ancestor of all three domains of life. The TBP domain is also found in RNase HIII, and phylogenetic analyses show that RNase HIII has evolved from bacterial RNase HII via TBP-domain fusion. Finally, our comparative genomic screens confirm and extend earlier reports of proteins consisting of a single TBP domain among some Archaea. These monopartite TBP-domain proteins suggest that this domain is functional in its own right, and that the TBP domain could have first evolved as an independent protein, which was later recruited in different contexts.

The case for an early biological origin of DNA.

All life generates deoxyribonucleotides, the building blocks of DNA, via ribonucleotide reductases (RNRs). The complexity of this reaction suggests it did not evolve until well after the advent of templated protein synthesis, which in turn suggests DNA evolved later than both RNA and templated protein synthesis. However, deoxyribonucleotides may have first been synthesised via an alternative, chemically simpler route--the reversal of the deoxyriboaldolase (DERA) step in deoxyribonucleotide salvage. In light of recent work demonstrating that this reaction can drive synthesis of deoxyribonucleosides, we consider what pressures early adoption of this pathway would have placed on cell metabolism. This in turn provides a rationale for the replacement of DERA-dependent DNA production by RNR-dependent production.

The Constructive Black Queen hypothesis: new functions can evolve under conditions favouring gene loss.

Duplication is a major route for the emergence of new gene functions. However, the emergence of new gene functions via this route may be reduced in prokaryotes, as redundant genes are often rapidly purged. In lineages with compact, streamlined genomes, it thus appears challenging for novel function to emerge via duplication and divergence. A further pressure contributing to gene loss occurs under Black Queen dynamics, as cheaters that lose the capacity to produce a public good can instead acquire it from neighbouring producers. We propose that Black Queen dynamics can favour the emergence of new function because, under an emerging Black Queen dynamic, there is high gene redundancy spread across a community of interacting cells. Using computational modelling, we demonstrate that new gene functions can emerge under Black Queen dynamics. This result holds even if there is deletion bias due to low duplication rates and selection against redundant gene copies resulting from the high cost associated with carrying a locus. However, when the public good production costs are high, Black Queen dynamics impede the fixation of new functions. Our results expand the mechanisms by which new gene functions can emerge in prokaryotic systems.

Use of structural phylogenetic networks for classification of the ferritin-like superfamily.

In the postgenomic era, bioinformatic analysis of sequence similarity is an immensely powerful tool to gain insight into evolution and protein function. Over long evolutionary distances, however, sequence-based methods fail as the similarities become too low for phylogenetic analysis. Macromolecular structure generally appears better conserved than sequence, but clear models for how structure evolves over time are lacking. The exponential growth of three-dimensional structural information may allow novel structure-based methods to drastically extend the evolutionary time scales amenable to phylogenetics and functional classification of proteins. To this end, we analyzed 80 structures from the functionally diverse ferritin-like superfamily. Using evolutionary networks, we demonstrate that structural comparisons can delineate and discover groups of proteins beyond the "twilight zone" where sequence similarity does not allow evolutionary analysis, suggesting that considerable and useful evolutionary signal is preserved in three-dimensional structures.

Comparative analysis of RNA families reveals distinct repertoires for each domain of life.

The RNA world hypothesis, that RNA genomes and catalysts preceded DNA genomes and genetically-encoded protein catalysts, has been central to models for the early evolution of life on Earth. A key part of such models is continuity between the earliest stages in the evolution of life and the RNA repertoires of extant lineages. Some assessments seem consistent with a diverse RNA world, yet direct continuity between modern RNAs and an RNA world has not been demonstrated for the majority of RNA families, and, anecdotally, many RNA functions appear restricted in their distribution. Despite much discussion of the possible antiquity of RNA families, no systematic analyses of RNA family distribution have been performed. To chart the broad evolutionary history of known RNA families, we performed comparative genomic analysis of over 3 million RNA annotations spanning 1446 families from the Rfam 10 database. We report that 99% of known RNA families are restricted to a single domain of life, revealing discrete repertoires for each domain. For the 1% of RNA families/clans present in more than one domain, over half show evidence of horizontal gene transfer (HGT), and the rest show a vertical trace, indicating the presence of a complex protein synthesis machinery in the Last Universal Common Ancestor (LUCA) and consistent with the evolutionary history of the most ancient protein-coding genes. However, with limited interdomain transfer and few RNA families exhibiting demonstrable antiquity as predicted under RNA world continuity, our results indicate that the majority of modern cellular RNA repertoires have primarily evolved in a domain-specific manner.

Characterisation of the trans-membrane nucleoporins GP210 and NDC1 in Arabidopsis thaliana.

The nuclear pore is structurally conserved across eukaryotes as are many of the pore's constituent proteins. The transmembrane nuclear pore proteins GP210 and NDC1 span the nuclear envelope holding the nuclear pore in place. Orthologues of GP210 and NDC1 in Arabidopsis were investigated through characterisation of T-DNA insertional mutants. While the T-DNA insert into GP210 reduced expression of the gene, the insert in the NDC1 gene resulted in increased expression in both the ndc1 mutant as well as the ndc1/gp210 double mutant. The ndc1 and gp210 individual mutants showed little phenotypic difference from wild-type plants, but the ndc1/gp210 mutant showed a range of phenotypic effects. As with many plant nuclear pore protein mutants, these effects included non-nuclear phenotypes such as reduced pollen viability, reduced growth and glabrous leaves in mature plants. Importantly, however, ndc1/gp210 exhibited nuclear-specific effects including modifications to nuclear shape in different cell types. We also observed functional changes to nuclear transport in ndc1/gp210 plants, with low levels of cytoplasmic fluorescence observed in cells expressing nuclear-targeted GFP. The lack of phenotypes in individual insertional lines, and the relatively mild phenotype suggests that additional transmembrane nucleoporins, such as the recently-discovered CPR5, likely compensate for their loss.

Structome: a tool for the rapid assembly of datasets for structural phylogenetics.

Summary: Protein structures carry signal of common ancestry and can therefore aid in reconstructing their evolutionary histories. To expedite the structure-informed inference process, a web server, Structome, has been developed that allows users to rapidly identify protein structures similar to a query protein and to assemble datasets useful for structure-based phylogenetics. Structome was created by clustering ∼94% of the structures in RCSB PDB using 90% sequence identity and representing each cluster by a centroid structure. Structure similarity between centroid proteins was calculated, and annotations from PDB, SCOP, and CATH were integrated. To illustrate utility, an H3 histone was used as a query, and results show that the protein structures returned by Structome span both sequence and structural diversity of the histone fold. Additionally, the pre-computed nexus-formatted distance matrix, provided by Structome, enables analysis of evolutionary relationships between proteins not identifiable using searches based on sequence similarity alone. Our results demonstrate that, beginning with a single structure, Structome can be used to rapidly generate a dataset of structural neighbours and allows deep evolutionary history of proteins to be studied. Availability and Implementation: Structome is available at: https://structome.bii.a-star.edu.sg.

Characterisation of an Escherichia coli line that completely lacks ribonucleotide reduction yields insights into the evolution of parasitism and endosymbiosis.

Life requires ribonucleotide reduction for de novo synthesis of deoxyribonucleotides. As ribonucleotide reduction has on occasion been lost in parasites and endosymbionts, which are instead dependent on their host for deoxyribonucleotide synthesis, it should in principle be possible to knock this process out if growth media are supplemented with deoxyribonucleosides. We report the creation of a strain of Escherichia coli where all three ribonucleotide reductase operons have been deleted following introduction of a broad spectrum deoxyribonucleoside kinase from Mycoplasma mycoides. Our strain shows slowed but substantial growth in the presence of deoxyribonucleosides. Under limiting deoxyribonucleoside levels, we observe a distinctive filamentous cell morphology, where cells grow but do not appear to divide regularly. Finally, we examined whether our lines can adapt to limited supplies of deoxyribonucleosides, as might occur in the switch from de novo synthesis to dependence on host production during the evolution of parasitism or endosymbiosis. Over the course of an evolution experiment, we observe a 25-fold reduction in the minimum concentration of exogenous deoxyribonucleosides necessary for growth. Genome analysis reveals that several replicate lines carry mutations in deoB and cdd. deoB codes for phosphopentomutase, a key part of the deoxyriboaldolase pathway, which has been hypothesised as an alternative to ribonucleotide reduction for deoxyribonucleotide synthesis. Rather than complementing the loss of ribonucleotide reduction, our experiments reveal that mutations appear that reduce or eliminate the capacity for this pathway to catabolise deoxyribonucleotides, thus preventing their loss via central metabolism. Mutational inactivation of both deoB and cdd is also observed in a number of obligate intracellular bacteria that have lost ribonucleotide reduction. We conclude that our experiments recapitulate key evolutionary steps in the adaptation to life without ribonucleotide reduction.

Comparative genomics of eukaryotic small nucleolar RNAs reveals deep evolutionary ancestry amidst ongoing intragenomic mobility.

BACKGROUND: Small nucleolar (sno)RNAs are required for posttranscriptional processing and modification of ribosomal, spliceosomal and messenger RNAs. Their presence in both eukaryotes and archaea indicates that snoRNAs are evolutionarily ancient. The location of some snoRNAs within the introns of ribosomal protein genes has been suggested to belie an RNA world origin, with the exons of the earliest protein-coding genes having evolved around snoRNAs after the advent of templated protein synthesis. Alternatively, this intronic location may reflect more recent selection for coexpression of snoRNAs and ribosomal components, ensuring rRNA modification by snoRNAs during ribosome synthesis. To gain insight into the evolutionary origins of this genetic organization, we examined the antiquity of snoRNA families and the stability of their genomic location across 44 eukaryote genomes. RESULTS: We report that dozens of snoRNA families are traceable to the Last Eukaryotic Common Ancestor (LECA), but find only weak similarities between the oldest eukaryotic snoRNAs and archaeal snoRNA-like genes. Moreover, many of these LECA snoRNAs are located within the introns of host genes independently traceable to the LECA. Comparative genomic analyses reveal the intronic location of LECA snoRNAs is not ancestral however, suggesting the pattern we observe is the result of ongoing intragenomic mobility. Analysis of human transcriptome data indicates that the primary requirement for hosting intronic snoRNAs is a broad expression profile. Consistent with ongoing mobility across broadly-expressed genes, we report a case of recent migration of a non-LECA snoRNA from the intron of a ubiquitously expressed non-LECA host gene into the introns of two LECA genes during the evolution of primates. CONCLUSIONS: Our analyses show that snoRNAs were a well-established family of RNAs at the time when eukaryotes began to diversify. While many are intronic, this association is not evolutionarily stable across the eukaryote tree; ongoing intragenomic mobility has erased signal of their ancestral gene organization, and neither introns-first nor evolved co-expression adequately explain our results. We therefore present a third model - constrained drift - whereby individual snoRNAs are intragenomically mobile and may occupy any genomic location from which expression satisfies phenotype.

Endosymbiont gene functions impaired and rescued by polymerase infidelity at poly(A) tracts.

Among host-dependent bacteria that have evolved by extreme reductive genome evolution, long-term bacterial endosymbionts of insects have the smallest (160-790 kb) and most A + T-rich (>70%) bacterial genomes known to date. These genomes are riddled with poly(A) tracts, and 5-50% of genes contain tracts of 10 As or more. Here, we demonstrate transcriptional slippage at poly(A) tracts within genes of Buchnera aphidicola associated with aphids and Blochmannia pennsylvanicus associated with ants. Several tracts contain single frameshift deletions; these apparent pseudogenes showed patterns of constraint consistent with purifying selection on the encoded proteins. Transcriptional slippage yielded a heterogeneous population of transcripts with variable numbers of As in the tract. Across several frameshifted genes, including B. aphidicola cell wall biosynthesis genes and a B. pennsylvanicus histidine biosynthesis gene, 12-50% of transcripts contained corrected reading frames that could potentially yield full-length proteins. In situ immunostaining confirmed the production of the cell wall biosynthetic enzyme UDP-N-acetylmuramyl pentapeptide synthase encoded by the frameshifted murF gene. Simulation studies indicated an overrepresentation of poly(A) tracts in endosymbiont genomes relative to other A + T-rich bacterial genomes. Polymerase infidelity at poly(A) tracts rescues the functionality of genes with frameshift mutations and, conversely, reduces the efficiency of expression for in-frame genes carrying poly(A) regions. These features of homopolymeric tracts could be exploited to manipulate gene expression in small synthetic genomes.