Experimental evolution links post-transcriptional regulation to Leishmania fitness gain.

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Genome instability drives epistatic adaptation in the human pathogen Leishmania.

How genome instability is harnessed for fitness gain despite its potential deleterious effects is largely elusive. An ideal system to address this important open question is provided by the protozoan pathogen Leishmania, which exploits frequent variations in chromosome and gene copy number to regulate expression levels. Using ecological genomics and experimental evolution approaches, we provide evidence that Leishmania adaptation relies on epistatic interactions between functionally associated gene copy number variations in pathways driving fitness gain in a given environment. We further uncover posttranscriptional regulation as a key mechanism that compensates for deleterious gene dosage effects and provides phenotypic robustness to genetically heterogenous parasite populations. Finally, we correlate dynamic variations in small nucleolar RNA (snoRNA) gene dosage with changes in ribosomal RNA 2'-O-methylation and pseudouridylation, suggesting translational control as an additional layer of parasite adaptation. Leishmania genome instability is thus harnessed for fitness gain by genome-dependent variations in gene expression and genome-independent compensatory mechanisms. This allows for polyclonal adaptation and maintenance of genetic heterogeneity despite strong selective pressure. The epistatic adaptation described here needs to be considered in Leishmania epidemiology and biomarker discovery and may be relevant to other fast-evolving eukaryotic cells that exploit genome instability for adaptation, such as fungal pathogens or cancer.

Identification and functional implications of pseudouridine RNA modification on small noncoding RNAs in the mammalian pathogen Trypanosoma brucei.

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vault RNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq that the Ψ on C/D snoRNA guiding 2'-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these noncoding RNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of noncoding RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

Effects of Vineyard Agro-management Practices on Soil Bacterial Community Composition, and Diversity.

Changes in land use strongly affect soil biological and physico-chemical structure and characteristics, which are strongly related to agricultural conversion of natural habitats to man-made usage. These are among the most important and not always beneficial changes, affecting loss of habitats. In Golan Heights basaltic soils, vineyards are currently a driving force in land-use change. Such changes could have an important effect on soil microbial community that play an important role in maintaining stable functioning of soil ecosystems. This study investigated the microbial communities in five different agro-managements using molecular tools that can clarify the differences in microbial community structure and function. Significant differences in soil microbial community composition were found. However, no differences in alpha diversity or functionality were found between the treatments. To the best of our knowledge, this is the first report indicating that the bacterial community in different agro-managements provide an insight into the potential function of a vineyard system.

Plant Gender Affects Soil Fungal Microbiota Associated with Welwitschia mirabilis, an Unusual Desert Gymnosperm.

In a recent study, we found a distinct soil bacterial community associated with male and female plants of the desert gymnosperm Welwitschia mirabilis. In this subsequent study, we also found that the soil fungal community associated with Welwitschia differs between male and female plants, and between unvegetated areas and the soil under plants. Site location, pH, and soil moisture also had an important influence on the composition of the fungal community. A number of Ascomycota and Chytrid species were found to be distinct indicators of male and female plants, respectively, but there was no overall difference at the phylum level or in terms of diversity. The unvegetated areas between plants also differed in terms of several Ascomycota OTUs. Network connectivity of the fungal communities was found to be higher under both male and female Welwitschia plants than in unvegetated control areas. As with the bacterial community, it is unclear what processes produce the gender-distinct fungal community, and also the more general plant-associated community, and also what the effects on the biology of the plants are. One possibility behind the gender-related difference in fungal community is that there are differences in the production of pollen or nectar between the two plant genders, affecting the below-ground soil community.

Pseudouridines on Trypanosoma brucei mRNAs are developmentally regulated: Implications to mRNA stability and protein binding.

The parasite Trypanosoma brucei cycles between an insect and a mammalian host and is the causative agent of sleeping sickness. Here, we performed high-throughput mapping of pseudouridines (Ψs) on mRNA from two life stages of the parasite. The analysis revealed ~273 Ψs, including developmentally regulated Ψs that are guided by homologs of pseudouridine synthases (PUS1, 3, 5, and 7). Mutating the U that undergoes pseudouridylation in the 3' UTR of valyl-tRNA synthetase destabilized the mRNA level. To investigate the mechanism by which Ψ affects the stability of this mRNA, proteins that bind to the 3' UTR were identified, including the RNA binding protein RBSR1. The binding of RBSR1 protein to the 3' UTR was stronger when lacking Ψ compared to transcripts carrying the modification, suggesting that Ψ can inhibit the binding of proteins to their target and thus affect the stability of mRNAs. Consequently, Ψ modification on mRNA adds an additional level of regulation to the dominant post-transcriptional control in these parasites.

Drosophila TRF2 is a preferential core promoter regulator.

Transcription of protein-coding genes is highly dependent on the RNA polymerase II core promoter. Core promoters, generally defined as the regions that direct transcription initiation, consist of functional core promoter motifs (such as the TATA-box, initiator [Inr], and downstream core promoter element [DPE]) that confer specific properties to the core promoter. The known basal transcription factors that support TATA-dependent transcription are insufficient for in vitro transcription of DPE-dependent promoters. In search of a transcription factor that supports DPE-dependent transcription, we used a biochemical complementation approach and identified the Drosophila TBP (TATA-box-binding protein)-related factor 2 (TRF2) as an enriched factor in the fractions that support DPE-dependent transcription. We demonstrate that the short TRF2 isoform preferentially activates DPE-dependent promoters. DNA microarray analysis reveals the enrichment of DPE promoters among short TRF2 up-regulated genes. Using primer extension analysis and reporter assays, we show the importance of the DPE in transcriptional regulation of TRF2 target genes. It was previously shown that, unlike TBP, TRF2 fails to bind DNA containing TATA-boxes. Using microfluidic affinity analysis, we discovered that short TRF2-bound DNA oligos are enriched for Inr and DPE motifs. Taken together, our findings highlight the role of short TRF2 as a preferential core promoter regulator.

Small nucleolar RNAs controlling rRNA processing in Trypanosoma brucei.

In trypanosomes, in contrast to most eukaryotes, the large subunit (LSU) ribosomal RNA is fragmented into two large and four small ribosomal RNAs (srRNAs) pieces, and this additional processing likely requires trypanosome-specific factors. Here, we examined the role of 10 abundant small nucleolar RNAs (snoRNAs) involved in rRNA processing. We show that each snoRNA involved in LSU processing associates with factors engaged in either early or late biogenesis steps. Five of these snoRNAs interact with the intervening sequences of rRNA precursor, whereas the others only guide rRNA modifications. The function of the snoRNAs was explored by silencing snoRNAs. The data suggest that the LSU rRNA processing events do not correspond to the order of rRNA transcription, and that srRNAs 2, 4 and 6 which are part of LSU are processed before srRNA1. Interestingly, the 6 snoRNAs that affect srRNA1 processing guide modifications on rRNA positions that span locations from the protein exit tunnel to the srRNA1, suggesting that these modifications may serve as check-points preceding the liberation of srRNA1. This study identifies the highest number of snoRNAs so far described that are involved in rRNA processing and/or rRNA folding and highlights their function in the unique trypanosome rRNA maturation events.

Pseudouridines on Trypanosoma brucei spliceosomal small nuclear RNAs and their implication for RNA and protein interactions.

The parasite Trypanosoma brucei, the causative agent of sleeping sickness, cycles between an insect and a mammalian host. Here, we investigated the presence of pseudouridines (Ψs) on the spliceosomal small nuclear RNAs (snRNAs), which may enable growth at the very different temperatures characterizing the two hosts. To this end, we performed the first high-throughput mapping of spliceosomal snRNA Ψs by small RNA Ψ-seq. The analysis revealed 42 Ψs on T. brucei snRNAs, which is the highest number reported so far. We show that a trypanosome protein analogous to human protein WDR79, is essential for guiding Ψ on snRNAs but not on rRNAs. snoRNA species implicated in snRNA pseudouridylation were identified by a genome-wide approach based on ligation of RNAs following in vivo UV cross-linking. snRNA Ψs are guided by single hairpin snoRNAs, also implicated in rRNA modification. Depletion of such guiding snoRNA by RNAi compromised the guided modification on snRNA and reduced parasite growth at elevated temperatures. We further demonstrate that Ψ strengthens U4/U6 RNA-RNA and U2B"/U2A' proteins-U2 snRNA interaction at elevated temperatures. The existence of single hairpin RNAs that modify both the spliceosome and ribosome RNAs is unique for these parasites, and may be related to their ability to cycle between their two hosts that differ in temperature.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA.

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Sex differences in testosterone reactivity and sensitivity in a non-model gerbil.

Although testosterone (T) is a key regulator in vertebrate development, physiology, and behaviour in both sexes, studies suggest that its regulation may be sex-specific. We measured circulating T levels in Baluchistan gerbils (Gerbillus nanus) in the field and in the lab all year round and found no significant sex differences. However, we observed sex differences in circulating T levels following gonadotropin-releasing hormone (GnRH) challenge and T implants in this non-model species. Whereas only males elevated T following a GnRH challenge, females had higher serum T concentrations following T implant insertion. These differences may be a result of different points of regulation along the hypothalamic-pituitary-gonadal (HPG) axis. Consequently, we examined sex differences in the mRNA expression of the androgen receptor (AR) in multiple brain regions. We identified AR and ß-actin sequences in assembled genomic sequences of members of the Gerbillinae, which were analogous to rat sequences, and designed primers for them. The distribution of the AR in G. nanus brain regions was similar to documented expression profiles in rodents. We found lower AR mRNA levels in females in the striatum. Additionally, G. nanus that experienced housing in mixed-sex pairs had higher adrenal AR expression than G. nanus that were housed alone. Regulation of the gerbil HPG axis may reflect evolutionary sex differences in life-history strategies, with males ready to reproduce when receptive females are available, while the possible reproductive costs associated with female T direct its regulation upstream.

Cucumber ovaries inhibited by dominant fruit express a dynamic developmental program, distinct from either senescence-determined or fruit-setting ovaries.

Cucurbits represent an attractive model to explore the dynamics of fruit set, whose regulation is not fully understood, despite its importance for yield determination. A fertilized ovary must integrate signals from distant plant parts and 'decide' whether to set fruit, or remain inhibited and later senesce. Here, we set out to characterize first-fruit inhibition (FFI), that is, the inhibitory effect of the first fruit on subsequent development of younger ovaries during pollination-induced and parthenocarpic fruit set. After the first fertilized ovaries set fruit, younger fertilized ovaries remained in a temporary state of inhibition. Such ovaries preserved their size and green color, and if the older fruit were removed within a 1-week reversibility window, they set fruit. The FFI effect was documented in both fertilized and parthenocarpic ovaries. We compared the gene expression profiles of pollinated ovaries (committed to set fruit) with respect to those affected by FFI, and to non-pollinated ovaries (undergoing senescence). The three fates of the ovaries were characterized by wide changes in gene expression, with several specific transcripts being up- or down-regulated in response to pollination, and to the presence of inhibitory fruit. Metabolic profiling was undertaken and integrated with the transcriptomic data in order to characterize early physiological changes that occur in post-anthesis ovaries in parthenocarpic and non-parthenocarpic genotypes. The combined results are discussed with respect to current models of fruit set and specifically with regard to FFI. Moreover, these metabolome and transcriptome data provide a valuable resource for studying ovary development and fruit set.

Sedimentary Marl mudstone as a substrate in a xeric environment revealed by microbiome analysis.

The sedimentary Marl mudstone soil is composed primarily of CaCO3, and is an important pedologic and geomorphologic element known as Marl, extensively dispersed in slopes and ridges in the northern Negev Desert, Israel. The wide Marl soil-layer areas are barren, with well-developed streamsides and no perennial vegetation cover. Soil systems in the Negev Desert have been widely studied, yet very little information was collected on Marl soils, and even less on the microbiome present in the Negev. Thus, an evaluation of the microbial-community inhabitants in a Marl soil layer was conducted in an attempt to distinguish between Marl with surface green mat and bare Marl soil layer. Our objective was to investigate the microbiome and abiotic components of the upper layer (0-5 cm) of Marl and Marl-with-green-mat soil collected in the Negev Desert. Plate-counting enabled the estimation of fungal and bacterial population size, while nested polymerase chain reaction (nPCR) and Ion Torrent sequencing were used to analyze biological diversity. The results indicate significant differences in microbial biomass and microbial-community diversity between Marl and Marl-green mat, despite similar pH levels. Further study is needed to enhance understanding of the activity of the biotic components and their contribution to slope stabilization.

Desert Perennial Shrubs Shape the Microbial-Community Miscellany in Laimosphere and Phyllosphere Space.

Microbial function, composition, and distribution play a fundamental role in ecosystem ecology. The interaction between desert plants and their associated microbes is expected to greatly affect their response to changes in this harsh environment. Using comparative analyses, we studied the impact of three desert shrubs, Atriplex halimus (A), Artemisia herba-alba (AHA), and Hammada scoparia (HS), on soil- and leaf-associated microbial communities. DNA extracted from the leaf surface and soil samples collected beneath the shrubs were used to study associated microbial diversity using a sequencing survey of variable regions of bacterial 16S rRNA and fungal ribosomal internal transcribed spacer (ITS1). We found that the composition of bacterial and fungal orders is plant-type-specific, indicating that each plant type provides a suitable and unique microenvironment. The different adaptive ecophysiological properties of the three plant species and the differential effect on their associated microbial composition point to the role of adaptation in the shaping of microbial diversity. Overall, our findings suggest a link between plant ecophysiological adaptation as a "temporary host" and the biotic-community parameters in extreme xeric environments.

Genome-wide analysis of small nucleolar RNAs of Leishmania major reveals a rich repertoire of RNAs involved in modification and processing of rRNA.

Trypanosomatids are protozoan parasites and the causative agent of infamous infectious diseases. These organisms regulate their gene expression mainly at the post-transcriptional level and possess characteristic RNA processing mechanisms. In this study, we analyzed the complete repertoire of Leishmania major small nucleolar (snoRNA) RNAs by performing RNA-seq analysis on RNAs that were affinity-purified using the C/D snoRNA core protein, SNU13, and the H/ACA core protein, NHP2. This study revealed a large collection of C/D and H/ACA snoRNAs, organized in gene clusters generally containing both snoRNA types. Abundant snoRNAs were identified and predicted to guide trypanosome-specific rRNA cleavages. The repertoire of snoRNAs was compared to that of the closely related Trypanosoma brucei, and 80% of both C/D and H/ACA molecules were found to have functional homologues. The comparative analyses elucidated how snoRNAs evolved to generate molecules with analogous functions in both species. Interestingly, H/ACA RNAs have great flexibility in their ability to guide modifications, and several of the RNA species can guide more than one modification, compensating for the presence of single hairpin H/ACA snoRNA in these organisms. Placing the predicted modifications on the rRNA secondary structure revealed hypermodification regions mostly in domains which are modified in other eukaryotes, in addition to trypanosome-specific modifications.

RNA-seq analysis of small RNPs in Trypanosoma brucei reveals a rich repertoire of non-coding RNAs.

The discovery of a plethora of small non-coding RNAs (ncRNAs) has fundamentally changed our understanding of how genes are regulated. In this study, we employed the power of deep sequencing of RNA (RNA-seq) to examine the repertoire of ncRNAs present in small ribonucleoprotein particles (RNPs) of Trypanosoma brucei, an important protozoan parasite. We identified new C/D and H/ACA small nucleolar RNAs (snoRNAs), as well as tens of putative novel non-coding RNAs; several of these are processed from trans-spliced and polyadenylated transcripts. The RNA-seq analysis provided information on the relative abundance of the RNAs, and their 5'- and 3'-termini. The study demonstrated that three highly abundant snoRNAs are involved in rRNA processing and highlight the unique trypanosome-specific repertoire of these RNAs. Novel RNAs were studied using in situ hybridization, association in RNP complexes, and 'RNA walk' to detect interaction with their target RNAs. Finally, we showed that the abundance of certain ncRNAs varies between the two stages of the parasite, suggesting that ncRNAs may contribute to gene regulation during the complex parasite's life cycle. This is the first study to provide a whole-genome analysis of the large repertoire of small RNPs in trypanosomes.

CLOCK evolved in cnidaria to synchronize internal rhythms with diel environmental cues.

The circadian clock enables anticipation of the day/night cycle in animals ranging from cnidarians to mammals. Circadian rhythms are generated through a transcription-translation feedback loop (TTFL or pacemaker) with CLOCK as a conserved positive factor in animals. However, CLOCK's functional evolutionary origin and mechanism of action in basal animals are unknown. In the cnidarian Nematostella vectensis, pacemaker gene transcript levels, including NvClk (the Clock ortholog), appear arrhythmic under constant darkness, questioning the role of NvCLK. Utilizing CRISPR/Cas9, we generated a NvClk allele mutant (NvClkΔ), revealing circadian behavior loss under constant dark (DD) or light (LL), while maintaining a 24 hr rhythm under light-dark condition (LD). Transcriptomics analysis revealed distinct rhythmic genes in wild-type (WT) polypsunder LD compared to DD conditions. In LD, NvClkΔ/Δ polyps exhibited comparable numbers of rhythmic genes, but were reduced in DD. Furthermore, under LD, the NvClkΔ/Δ polyps showed alterations in temporal pacemaker gene expression, impacting their potential interactions. Additionally, differential expression of non-rhythmic genes associated with cell division and neuronal differentiation was observed. These findings revealed that a light-responsive pathway can partially compensate for circadian clock disruption, and that the Clock gene has evolved in cnidarians to synchronize rhythmic physiology and behavior with the diel rhythm of the earth's biosphere.

Structural and mechanistic insights into the function of Leishmania ribosome lacking a single pseudouridine modification.

Leishmania is the causative agent of cutaneous and visceral diseases affecting millions of individuals worldwide. Pseudouridine (Ψ), the most abundant modification on rRNA, changes during the parasite life cycle. Alterations in the level of a specific Ψ in helix 69 (H69) affected ribosome function. To decipher the molecular mechanism of this phenotype, we determine the structure of ribosomes lacking the single Ψ and its parental strain at ∼2.4-3 Å resolution using cryo-EM. Our findings demonstrate the significance of a single Ψ on H69 to its structure and the importance for its interactions with helix 44 and specific tRNAs. Our study suggests that rRNA modification affects translation of mRNAs carrying codon bias due to selective accommodation of tRNAs by the ribosome. Based on the high-resolution structures, we propose a mechanism explaining how the ribosome selects specific tRNAs.

The long non-coding RNA ERIC is regulated by E2F and modulates the cellular response to DNA damage.

BACKGROUND: The human genome encodes thousands of unique long non-coding RNAs (lncRNAs), and these transcripts are emerging as critical regulators of gene expression and cell fate. However, the transcriptional regulation of their expression is not fully understood. The pivotal transcription factor E2F1 which can induce both proliferation and cell death, is a critical downstream target of the tumor suppressor, RB. The retinoblastoma pathway is often inactivated in human tumors resulting in deregulated E2F activity. RESULTS: Here, we report that lncRNA XLOC 006942, which we named ERIC, is regulated by E2F1 and, most probably, also E2F3. We show that expression levels of ERIC were elevated upon activation of exogenous E2F1, E2F3 or endogenous E2Fs. Moreover, knockdown of either E2F1 or E2F3 reduced ERIC levels and endogenous E2F1 binds ERIC's promoter. Expression of ERIC was cell cycle regulated and peaked in G1 in an E2F1-dependent manner. Inhibition of ERIC expression increased E2F1-mediated apoptosis, suggesting that E2F1 and ERIC constitute a negative feedback loop that modulates E2F1 activity. Furthermore, ERIC levels were increased following DNA damage by the chemotherapeutic drug Etoposide, and inhibition of ERIC expression enhanced Etoposide -induced apoptosis. CONCLUSIONS: Our data identify ERIC as a novel lncRNA that is transcriptionally regulated by E2Fs, and restricts apoptosis induced by E2F1, as well as by DNA damage.

A single DPE core promoter motif contributes to in vivo transcriptional regulation and affects cardiac function.

Transcription is initiated at the core promoter, which confers specific functions depending on the unique combination of core promoter elements. The downstream core promoter element (DPE) is found in many genes related to heart and mesodermal development. However, the function of these core promoter elements has thus far been studied primarily in isolated, in vitro or reporter gene settings. tinman (tin) encodes a key transcription factor that regulates the formation of the dorsal musculature and heart. Pioneering a novel approach utilizing both CRISPR and nascent transcriptomics, we show that a substitution mutation of the functional tin DPE motif within the natural context of the core promoter results in a massive perturbation of Tinman's regulatory network orchestrating dorsal musculature and heart formation. Mutation of endogenous tin DPE reduced the expression of tin and distinct target genes, resulting in significantly reduced viability and an overall decrease in adult heart function. We demonstrate the feasibility and importance of characterizing DNA sequence elements in vivo in their natural context, and accentuate the critical impact a single DPE motif has during Drosophila embryogenesis and functional heart formation.