Conservation of divergent transcription in fungi

Regulatory mechanisms ensuring coordinated expression of functionally related genes

Coordinated spatiotemporal expression of large sets of genes is required for the development and homeostasis of organisms. To achieve this goal, organisms use myriad strategies where they form operons, utilize bidirectional promoters, cluster genes, share enhancers among genes by DNA looping, and form topologically associated domains and transcriptional condensates. Coexpression achieved by these different strategies is hypothesized to have functional importance in minimizing gene expression variability, establishing dosage balance to ensure stoichiometry of protein complexes, and minimizing accumulation of toxic intermediate metabolites. By combining gene-editing tools with computational modeling, recent studies tested the advantages of adjacent genes located in pairs and clusters. We propose that with the advancement of gene editing, single-cell sequencing, and imaging tools, one could readily test the functional importance of different coexpression strategies in a variety of biological processes.

Genomic analyses reveal evolutionary and geologic context for the plateau fungus Ophiocordyceps sinensis

Background

Ophiocordyceps sinensis, which is only naturally found in the high-elevation extreme environment of the Tibetan Plateau, has been used in traditional Chinese medicine. Information concerning the evolutionary and geologic context of O. sinensis remains limited, however.

Methods

We constructed the high-quality genome of O. sinensis and provided insight into the evolution and ecology of O. sinensis using comparative genomics.

Results

We mapped the whole genome of the anamorph/asexual form Hirsutella of O. sinensis using Illumina and PacBio sequencing technologies and obtained a well assembled genome of 119.2 Mbp size. Long-read Single Molecule Real Time (SMRT) sequencing technology generated an assembly with more accurate representation of repeat sequence abundances and placement. Evolutionary analyses indicated that O. sinensis diverged from other fungi 65.9 Mya in the Upper Cretaceous, during the uplift of the Tibetan Plateau. Gene family expansions and contractions in addition to genome inflation via long terminal repeat (LTR) retrotransposon insertions were implicated as an important driver of O. sinensis divergence. The insertion rate of LTR sequences into the O. sinensis genome peaked ~ 30-40 Mya, when the Tibetan Plateau rose rapidly. Gene Ontology (GO) enrichment analysis suggested that O. sinensis contained more genes related to ice binding compared to other closely related fungi, which may aid in their adaptability to the cold Tibetan Plateau. Further, heavy metal resistance genes were in low abundance in the O. sinensis genome, which may help to explain previous observations that O. sinensis tissues contain high levels of heavy metals.

Conclusions

Our results reveal the evolutionary, geological, and ecological context for the evolution of the O. sinensis genome and the factors that have contributed to the environmental adaptability of this valuable fungus. These findings suggest that genome inflation via LTR retrotransposon insertions in O. sinensis coincided with the uplift of the Tibetan Plateau. LTRs and the specific genetic mechanisms of O. sinensis contributed to its adaptation to the environment on the plateau.

Additional Layer of Regulation via Convergent Gene Orientation in Yeasts

Convergent gene pairs can produce transcripts with complementary sequences. We had shown that mRNA duplexes form in vivo in Saccharomyces cerevisiae via interactions of mRNA overlapping 3′-ends and can lead to posttranscriptional regulatory events. Here we show that mRNA duplex formation is restricted to convergent genes separated by short intergenic distance, independently of their 3′-untranslated region (UTR) length. We disclose an enrichment in genes involved in biological processes related to stress among these convergent genes. They are markedly conserved in convergent orientation in budding yeasts, meaning that this mode of posttranscriptional regulation could be shared in these organisms, conferring an additional level for modulating stress response. We thus investigated the mechanistic advantages potentially conferred by 3′-UTR mRNA interactions. Analysis of genome-wide transcriptome data revealed that Pat1 and Lsm1 factors, having 3′-UTR binding preference and participating to the remodeling of messenger ribonucleoprotein particles, bind differently these messenger-interacting mRNAs forming duplexes in comparison to mRNAs that do not interact (solo mRNAs). Functionally, messenger-interacting mRNAs show limited translational repression upon stress. We thus propose that mRNA duplex formation modulates the regulation of mRNA expression by limiting their access to translational repressors. Our results thus show that posttranscriptional regulation is an additional factor that determines the order of coding genes.

The Glycerol Phosphatase Gpp2: A Link to Osmotic Stress, Sulfur Assimilation and Virulence in Cryptococcus neoformans

Cryptococcus neoformans is an opportunist fungal pathogen that causes meningoencephalitis in immunocompromised patients. During infection, this basidiomycete yeast has to adapt to several adverse conditions, especially nutrient availability. The interruption on various amino acid biosynthetic pathways and on amino acid uptake causes reduced viability, inability to cope with various stresses, failure in virulence factors expression and avirulence in animal model of infection. The sulfur amino acid biosynthesis and uptake is an important feature for pathogen survival in vivo and in vitro. Our previous work demonstrates that C. neoformans Cys3 BZip transcription factor controls the gene expression in several steps of the sulfur assimilation and sulfur amino acid biosynthesis. Also, we have shown that Gpp2 phosphatase modulates Cys3 activity. In Saccharomyces cerevisiae Gpp2 is induced in response to hyper osmotic or oxidative stress and during diauxic shift. In this work, we will show that, in C. neoformans, Gpp2 is required to respond to stresses, mainly osmotic stress; also its transcription is induced during exposure to NaCl. Global transcriptional profile of gpp2Δ by RNAseq shows that CYS3 and other genes in the sulfur assimilation pathway are up regulated, which is consistent with our previous report, in which Gpp2 acts by avoiding Cys3 accumulation and nuclear localization. In addition, several transporters genes, especially amino acid permeases and oxidative stress genes are induced in the gpp2Δ strain; on the contrary, genes involved in glucose and tricarboxylic acid metabolism are down regulated. gpp2Δ strain fails to express virulence factors, as melanin, phospholipase, urease and has virulence attenuation in Galleria mellonella. Our data suggest that Gpp2 is an important factor for general pathogen adaptation to various stresses and also to the host, and perhaps it could be an interesting target for therapeutic use.

The birth, evolution and death of metabolic gene clusters in fungi

Fungi contain a remarkable diversity of both primary and secondary metabolic pathways involved in ecologically specialized or accessory functions. Genes in these pathways are frequently physically linked on fungal chromosomes, forming metabolic gene clusters (MGCs). In this Review, we describe the diversity in the structure and content of fungal MGCs, their population-level and species-level variation, the evolutionary mechanisms that underlie their formation, maintenance and decay, and their ecological and evolutionary impact on fungal populations. We also discuss MGCs from other eukaryotes and the reasons for their preponderance in fungi. Improved knowledge of the evolutionary life cycle of MGCs will advance our understanding of the ecology of specialized metabolism and of the interplay between the lifestyle of an organism and genome architecture.

A Robust Phylogenomic Time Tree for Biotechnologically and Medically Important Fungi in the Genera Aspergillus and Penicillium

Understanding the evolution of traits across technologically and medically significant fungi requires a robust phylogeny. Even though species in the Aspergillus and Penicillium genera (family Aspergillaceae, class Eurotiomycetes) are some of the most significant technologically and medically relevant fungi, we still lack a genome-scale phylogeny of the lineage or knowledge of the parts of the phylogeny that exhibit conflict among analyses. Here, we used a phylogenomic approach to infer evolutionary relationships among 81 genomes that span the diversity of Aspergillus and Penicillium species, to identify conflicts in the phylogeny, and to determine the likely underlying factors of the observed conflicts. Using a data matrix comprised of 1,668 genes, we found that while most branches of the phylogeny of the Aspergillaceae are robustly supported and recovered irrespective of method of analysis, a few exhibit various degrees of conflict among our analyses. Further examination of the observed conflict revealed that it largely stems from incomplete lineage sorting and hybridization or introgression. Our analyses provide a robust and comprehensive evolutionary genomic roadmap for this important lineage, which will facilitate the examination of the diverse technologically and medically relevant traits of these fungi in an evolutionary context.

The filamentous fungal family Aspergillaceae contains >1,000 known species, mostly in the genera Aspergillus and Penicillium. Several species are used in the food, biotechnology, and drug industries (e.g., Aspergillus oryzae and Penicillium camemberti), while others are dangerous human and plant pathogens (e.g., Aspergillus fumigatus and Penicillium digitatum). To infer a robust phylogeny and pinpoint poorly resolved branches and their likely underlying contributors, we used 81 genomes spanning the diversity of Aspergillus and Penicillium to construct a 1,668-gene data matrix. Phylogenies of the nucleotide and amino acid versions of this full data matrix as well as of several additional data matrices were generated using three different maximum likelihood schemes (i.e., gene-partitioned, unpartitioned, and coalescence) and using both site-homogenous and site-heterogeneous models (total of 64 species-level phylogenies). Examination of the topological agreement among these phylogenies and measures of internode certainty identified 11/78 (14.1%) bipartitions that were incongruent and pinpointed the likely underlying contributing factors, which included incomplete lineage sorting, hidden paralogy, hybridization or introgression, and reconstruction artifacts associated with poor taxon sampling. Relaxed molecular clock analyses suggest that Aspergillaceae likely originated in the lower Cretaceous and that the Aspergillus and Penicillium genera originated in the upper Cretaceous. Our results shed light on the ongoing debate on Aspergillus systematics and taxonomy and provide a robust evolutionary and temporal framework for comparative genomic analyses in Aspergillaceae. More broadly, our approach provides a general template for phylogenomic identification of resolved and contentious branches in densely genome-sequenced lineages across the tree of life.

Neighboring genes are closely related to whole genome duplications after their separation

BACKGROUND

The gene order in a eukaryotic genome is not random. Some neighboring genes show specific similarities, while others become separated during evolution. Whole genome duplication events (WGDs) have been recognized as an important evolutionary force. The potential relationship between the separation of neighboring genes and WGDs needs to be investigated. In this study, we investigated whether there is a potential relationship between separated neighboring gene pairs and WGDs, and the mechanism by which neighboring genes are separated. Additionally, we studied whether neighboring genes tend to show intrachromosomal colocalization after their neighborhood was disrupted and the factors facilitating the intrachromosomal colocalization of separated neighboring genes.

RESULTS

The separation of neighboring gene pairs is closely related to whole genome duplication events. Furthermore, we found that there is a double linear relationship between separated neighboring genes, total genes, and WGDs. The process of separation of neighboring genes caused by WGDs is also not random but abides by the double linear model. Separated neighboring gene pairs tend to show intrachromosomal colocalization. The conservativism of separated neighboring genes and histone modification facilitate the intrachromosomal colocalization of neighboring genes after their separation.

CONCLUSIONS

These results provide new insight into the understanding of evolutionary roles of locations and the relationship of neighboring gene pairs with whole genome duplications. Furthermore, understanding the proposed mechanism for intrachromosomal colocalization of separated genes benefits our knowledge of chromosomal interactions in the nucleus.

A Molecular Portrait of De Novo Genes in Yeasts

New genes, with novel protein functions, can evolve "from scratch" out of intergenic sequences. These de novo genes can integrate the cell's genetic network and drive important phenotypic innovations. Therefore, identifying de novo genes and understanding how the transition from noncoding to coding occurs are key problems in evolutionary biology. However, identifying de novo genes is a difficult task, hampered by the presence of remote homologs, fast evolving sequences and erroneously annotated protein coding genes. To overcome these limitations, we developed a procedure that handles the usual pitfalls in de novo gene identification and predicted the emergence of 703 de novo gene candidates in 15 yeast species from 2 genera whose phylogeny spans at least 100 million years of evolution. We validated 85 candidates by proteomic data, providing new translation evidence for 25 of them through mass spectrometry experiments. We also unambiguously identified the mutations that enabled the transition from noncoding to coding for 30 Saccharomyces de novo genes. We established that de novo gene origination is a widespread phenomenon in yeasts, only a few being ultimately maintained by selection. We also found that de novo genes preferentially emerge next to divergent promoters in GC-rich intergenic regions where the probability of finding a fortuitous and transcribed ORF is the highest. Finally, we found a more than 3-fold enrichment of de novo genes at recombination hot spots, which are GC-rich and nucleosome-free regions, suggesting that meiotic recombination contributes to de novo gene emergence in yeasts.

Kynurenine 3-Monooxygenase Gene Associated With Nicotine Initiation and Addiction: Analysis of Novel Regulatory Features at 5' and 3'-Regions

Tobacco smoking is widespread behavior in Qatar and worldwide and is considered one of the major preventable causes of ill health and death. Nicotine is part of tobacco smoke that causes numerous health risks and is incredibly addictive; it binds to the α7 nicotinic acetylcholine receptor (α7nAChR) in the brain. Recent studies showed α7nAChR involvement in the initiation and addiction of smoking. Kynurenic acid (KA), a significant tryptophan metabolite, is an antagonist of α7nAChR. Inhibition of kynurenine 3-monooxygenase enzyme encoded by KMO enhances the KA levels. Modulating KMO gene expression could be a useful tactic for the treatment of tobacco initiation and dependence. Since KMO regulation is still poorly understood, we aimed to investigate the 5' and 3'-regulatory factors of KMO gene to advance our knowledge to modulate KMO gene expression. In this study, bioinformatics methods were used to identify the regulatory sequences associated with expression of KMO. The displayed differential expression of KMO mRNA in the same tissue and different tissues suggested the specific usage of the KMO multiple alternative promoters. Eleven KMO alternative promoters identified at 5'-regulatory region contain TATA-Box, lack CpG Island (CGI) and showed dinucleotide base-stacking energy values specific to transcription factor binding sites (TFBSs). The structural features of regulatory sequences can influence the transcription process and cell type-specific expression. The uncharacterized LOC105373233 locus coding for non-coding RNA (ncRNA) located on the reverse strand in a convergent manner at the 3'-side of KMO locus. The two genes likely expressed by a promoter that lacks TATA-Box harbor CGI and two TFBSs linked to the bidirectional transcription, the NRF1, and ZNF14 motifs. We identified two types of microRNA (miR) in the uncharacterized LOC105373233 ncRNA, which are like hsa-miR-5096 and hsa-miR-1285-3p and can target the miR recognition element (MRE) in the KMO mRNA. Pairwise sequence alignment identified 52 nucleotides sequence hosting MRE in the KMO 3' UTR untranslated region complementary to the ncRNA LOC105373233 sequence. We speculate that the identified miRs can modulate the KMO expression and together with alternative promoters at the 5'-regulatory region of KMO might contribute to the development of novel diagnostic and therapeutic algorithm for tobacco smoking.

Genomic Organization and Expression of Iron Metabolism Genes in the Emerging Pathogenic Mold Scedosporium apiospermum

The ubiquitous mold Scedosporium apiospermum is increasingly recognized as an emerging pathogen, especially among patients with underlying disorders such as immunodeficiency or cystic fibrosis (CF). Indeed, it ranks the second among the filamentous fungi colonizing the respiratory tract of CF patients. However, our knowledge about virulence factors of this fungus is still limited. The role of iron-uptake systems may be critical for establishment of Scedosporium infections, notably in the iron-rich environment of the CF lung. Two main strategies are employed by fungi to efficiently acquire iron from their host or from their ecological niche: siderophore production and reductive iron assimilation (RIA) systems. The aim of this study was to assess the existence of orthologous genes involved in iron metabolism in the recently sequenced genome of S. apiospermum. At first, a tBLASTn analysis using A. fumigatus iron-related proteins as query revealed orthologs of almost all relevant loci in the S. apiospermum genome. Whereas the genes putatively involved in RIA were randomly distributed, siderophore biosynthesis and transport genes were organized in two clusters, each containing a non-ribosomal peptide synthetase (NRPS) whose orthologs in A. fumigatus have been described to catalyze hydroxamate siderophore synthesis. Nevertheless, comparative genomic analysis of siderophore-related clusters showed greater similarity between S. apiospermum and phylogenetically close molds than with Aspergillus species. The expression level of these genes was then evaluated by exposing conidia to iron starvation and iron excess. The expression of several orthologs of A. fumigatus genes involved in siderophore-based iron uptake or RIA was significantly induced during iron starvation, and conversely repressed in iron excess conditions. Altogether, these results indicate that S. apiospermum possesses the genetic information required for efficient and competitive iron uptake. They also suggest an important role of the siderophore production system in iron uptake by S. apiospermum.

Major sulfonate transporter Soa1 in Saccharomyces cerevisiae and considerable substrate diversity in its fungal family

Sulfate is a well-established sulfur source for fungi; however, in soils sulfonates and sulfate esters, especially choline sulfate, are often much more prominent. Here we show that Saccharomyces cerevisiae YIL166C(SOA1) encodes an inorganic sulfur (sulfate, sulfite and thiosulfate) transporter that also catalyses sulfonate and choline sulfate uptake. Phylogenetic analysis of fungal SOA1 orthologues and expression of 20 members in the sul1Δ sul2Δ soa1Δ strain, which is deficient in inorganic and organic sulfur compound uptake, reveals that these transporters have diverse substrate preferences for sulfur compounds. We further show that SOA2, a S. cerevisiae SOA1 paralogue found in S. uvarum, S. eubayanus and S. arboricola is likely to be an evolutionary remnant of the uncharacterized open reading frames YOL163W and YOL162W. Our work highlights the importance of sulfonates and choline sulfate as sulfur sources in the natural environment of S. cerevisiae and other fungi by identifying fungal transporters for these compounds.

Sulfonates are a major source of sulphur for soil microbes but their cellular uptake is still not fully understood. Here the authors show that Saccharomyces cerevisiae YIL166C(SOA1) encodes for an inorganic sulphur transporter that can also function as a sulfonate and choline sulphate transporter.

A low-temperature-responsive element involved in the regulation of the Arabidopsis thaliana At1g71850/At1g71860 divergent gene pair

The bidirectional promoter of the Arabidopsis thaliana gene pair At1g71850/At1g71860 harbors low-temperature-responsive elements, which participate in anti-correlated transcription regulation of the driving genes in response to environmental low temperature. A divergent gene pair is defined as two adjacent genes organized head to head in opposite orientation, sharing a common promoter region. Divergent gene pairs are mainly coexpressed, but some display opposite regulation. The mechanistic basis of such anti-correlated regulation is not well understood. Here, the regulation of the Arabidopsis thaliana gene pair At1g71850/At1g71860 was investigated. Semi-quantitative RT-PCR and Genevestigator analyses showed that while one of the pair was upregulated by exposure to low temperature, the same treatment downregulated the other. Promoter::GUS fusion transgenes were used to show that this behavior was driven by a bidirectional promoter, which harbored an as-1 motif, associated with the low-temperature response; mutation of this sequence produced a significant decrease in cold-responsive expression. With regard to the as-1 motif in the native orientation repressing the promoter's low-temperature responsiveness, the same as-1 motif introduced in the reverse direction showed a slight enhancement in the promoter's responsiveness to low-temperature exposure, indicating that the orientation of the motif was important for the promoter's activity. These findings provide new insights into the complex transcriptional regulation of bidirectional gene pairs as well as plant stress response.

Regulation of cell-to-cell variability in divergent gene expression

Cell-to-cell variability (noise) is an important feature of gene expression that impacts cell fitness and development. The regulatory mechanism of this variability is not fully understood. Here we investigate the effect on gene expression noise in divergent gene pairs (DGPs). We generated reporters driven by divergent promoters, rearranged their gene order, and probed their expressions using time-lapse fluorescence microscopy and single-molecule fluorescence in situ hybridization (smFISH). We show that two genes in a co-regulated DGP have higher expression covariance compared with the separate, tandem and convergent configurations, and this higher covariance is caused by more synchronized firing of the divergent transcriptions. For differentially regulated DGPs, the regulatory signal of one gene can stochastically ‘leak' to the other, causing increased gene expression noise. We propose that the DGPs' function in limiting or promoting gene expression noise may enhance or compromise cell fitness, providing an explanation for the conservation pattern of DGPs.

Gene expression noise affects cell fitness and development. Here, Yan et al. show that co-regulated divergent gene pairs (DGPs) suppress uncorrelated gene expression noise due to more synchronized transcription firing, and differentially regulated DGPs enhance gene expression noise due to transcription leakage.

The AVR2-SIX5 gene pair is required to activate I-2-mediated immunity in tomato

Plant-invading microbes betray their presence to a plant by exposure of antigenic molecules such as small, secreted proteins called 'effectors'. In Fusarium oxysporum f. sp. lycopersici (Fol) we identified a pair of effector gene candidates, AVR2-SIX5, whose expression is controlled by a shared promoter. The pathogenicity of AVR2 and SIX5 Fol knockouts was assessed on susceptible and resistant tomato (Solanum lycopersicum) plants carrying I-2. The I-2 NB-LRR protein confers resistance to Fol races carrying AVR2. Like Avr2, Six5 was found to be required for full virulence on susceptible plants. Unexpectedly, each knockout could breach I-2-mediated disease resistance. So whereas Avr2 is sufficient to induce I-2-mediated cell death, Avr2 and Six5 are both required for resistance. Avr2 and Six5 interact in yeast two-hybrid assays as well as in planta. Six5 and Avr2 accumulate in xylem sap of plants infected with the reciprocal knockouts, showing that lack of I-2 activation is not due to a lack of Avr2 accumulation in the SIX5 mutant. The effector repertoire of a pathogen determines its host specificity and its ability to manipulate plant immunity. Our findings challenge an oversimplified interpretation of the gene-for-gene model by showing requirement of two fungal genes for immunity conferred by one resistance gene.

Intersecting transcription networks constrain gene regulatory evolution

Epistasis—the non-additive interactions between different genetic loci—constrains evolutionary pathways, blocking some and permitting others^1–8. For biological networks such as transcription circuits, the nature of these constraints and their consequences are largely unknown. Here we describe the evolutionary pathways of a transcription network that controls the response to mating pheromone in yeasts⁹. A component of this network, the transcription regulator Ste12, has evolved two different modes of binding to a set of its target genes. In one group of species, Ste12 binds to specific DNA binding sites, while in another lineage it occupies DNA indirectly, relying on a second transcription regulator to recognize DNA. We show, through the construction of various possible evolutionary intermediates, that evolution of the direct mode of DNA binding was not directly accessible to the ancestor. Instead, it was contingent on a lineage-specific change to an overlapping transcription network with a different function, the specification of cell type. These results show that analyzing and predicting the evolution of cis-regulatory regions requires an understanding of their positions in overlapping networks, as this placement constrains the available evolutionary pathways.

Uncovering the functional constraints underlying the genomic organization of the odorant-binding protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

A highly abundant bacteriophage discovered in the unknown sequences of human faecal metagenomes

Metagenomics, or sequencing of the genetic material from a complete microbial community, is a promising tool to discover novel microbes and viruses. Viral metagenomes typically contain many unknown sequences. Here we describe the discovery of a previously unidentified bacteriophage present in the majority of published human faecal metagenomes, which we refer to as crAssphage. Its ~97 kbp genome is six times more abundant in publicly available metagenomes than all other known phages together; it comprises up to 90% and 22% of all reads in virus-like particle (VLP)-derived metagenomes and total community metagenomes, respectively; and it totals 1.68% of all human faecal metagenomic sequencing reads in the public databases. The majority of crAssphage-encoded proteins match no known sequences in the database, which is why it was not detected before. Using a new co-occurrence profiling approach, we predict a Bacteroides host for this phage, consistent with Bacteroides-related protein homologues and a unique carbohydrate-binding domain encoded in the phage genome.

Metagenomic studies of microbial communities often report DNA sequences from unidentified viruses. Here, Dutilh et al. analyse metagenomic data to reveal the complete genome of an abundant, ubiquitous virus from human faeces, and predict that the virus infects bacteria of the Bacteroides group.

Neighboring genes show interchromosomal colocalization after their separation

The order of genes on eukaryotic chromosomes is nonrandom. Some neighboring genes show order conservation among species, while some neighboring genes separate during evolution. Here, we investigated whether neighboring genes show interactions after their separation. We found that neighboring gene pairs tend to show interchromosomal colocalization (i.e., nuclear colocalization) in the species in which they separate. These nuclear colocalized separated neighboring gene pairs 1) show neighborhood conservation in more species, 2) tend to be regulated by the same transcription factor, and 3) tend to be regulated by the same histone modification. These results suggest a mechanism by which neighboring genes could retain nuclear proximity after their separation.

Physical linkage of metabolic genes in fungi is an adaptation against the accumulation of toxic intermediate compounds

Genomic analyses have proliferated without being tied to tangible phenotypes. For example, although coordination of both gene expression and genetic linkage have been offered as genetic mechanisms for the frequently observed clustering of genes participating in fungal metabolic pathways, elucidation of the phenotype(s) favored by selection, resulting in cluster formation and maintenance, has not been forthcoming. We noted that the cause of certain well-studied human metabolic disorders is the accumulation of toxic intermediate compounds (ICs), which occurs when the product of an enzyme is not used as a substrate by a downstream neighbor in the metabolic network. This raises the hypothesis that the phenotype favored by selection to drive gene clustering is the mitigation of IC toxicity. To test this, we examined 100 diverse fungal genomes for the simplest type of cluster, gene pairs that are both metabolic neighbors and chromosomal neighbors immediately adjacent to each other, which we refer to as "double neighbor gene pairs" (DNGPs). Examination of the toxicity of their corresponding ICs shows that, compared with chromosomally nonadjacent metabolic neighbors, DNGPs are enriched for ICs that have acutely toxic LD50 doses or reactive functional groups. Furthermore, DNGPs are significantly more likely to be divergently oriented on the chromosome; remarkably, ∼40% of these DNGPs have ICs known to be toxic. We submit that the structure of synteny in metabolic pathways of fungi is a signature of selection for protection against the accumulation of toxic metabolic intermediates.

Explaining microbial phenotypes on a genomic scale: GWAS for microbes

There is an increasing availability of complete or draft genome sequences for microbial organisms. These data form a potentially valuable resource for genotype-phenotype association and gene function prediction, provided that phenotypes are consistently annotated for all the sequenced strains. In this review, we address the requirements for successful gene-trait matching. We outline a basic protocol for microbial functional genomics, including genome assembly, annotation of genotypes (including single nucleotide polymorphisms, orthologous groups and prophages), data pre-processing, genotype-phenotype association, visualization and interpretation of results. The methodologies for association described herein can be applied to other data types, opening up possibilities to analyze transcriptome-phenotype associations, and correlate microbial population structure or activity, as measured by metagenomics, to environmental parameters.

Impact of chromosomal inversions on the yeast DAL cluster

Chromosomal rearrangements occur readily in nature and are a major reshaping force during genome evolution. Such large scale modifications are usually deleterious causing several fitness defects, but sometimes can confer an advantage and become adaptive. For example the DAL metabolic cluster in yeast was assembled in recent evolutionary times in the Hemiascomycetes lineage, through a set of rearrangements that brought together the genes involved in the allantoin degradation pathway. In eukaryotes, the existence of physical clustering of genes with related functions supports the notion that neighbouring ORFs tend to be co-expressed and that the order of genes along the chromosomes may have biological significance, rather than being random as previously believed. In this study, we investigate the phenotypic effect that inversions have on the DAL gene cluster, expressed during nitrogen starvation. In all Saccharomyces "sensu stricto" species the order of the DAL cluster is conserved, while in the "sensu lato" species Naumovia castellii, which grows significantly worse than S. cerevisiae on allantoin, the cluster includes two nested inversions encompassing three DAL genes. We constructed several inverted and non-inverted S. cerevisiae strains possessing different inversions including those to mimic the configuration of the N. castellii DAL cluster. We showed that the inversion of DAL2 lower its own expression and reduces yeast fitness during nitrogen starvation. This rearrangement also altered the expression of the neighbouring genes DAL1 and DAL4. Moreover, we showed that the expression of the DAL4 anti-sense transcript (SUT614) does not change upon inversions of DAL2 and therefore is unlikely to be involved in its regulation. These results show that the order of the DAL cluster has an impact on the phenotype and gene expression, suggesting that these rearrangements may have been adaptive in the "sensu stricto" group in relation to the low availability of nitrogen in the environment.

Long noncoding RNAs in C. elegans

Thousands of long noncoding RNAs (lncRNAs) have been found in vertebrate animals, a few of which have known biological roles. To better understand the genomics and features of lncRNAs in invertebrates, we used available RNA-seq, poly(A)-site, and ribosome-mapping data to identify lncRNAs of Caenorhabditis elegans. We found 170 long intervening ncRNAs (lincRNAs), which had single- or multiexonic structures that did not overlap protein-coding transcripts, and about sixty antisense lncRNAs (ancRNAs), which were complementary to protein-coding transcripts. Compared to protein-coding genes, the lncRNA genes tended to be expressed in a stage-dependent manner. Approximately 25% of the newly identified lincRNAs showed little signal for sequence conservation and mapped antisense to clusters of endogenous siRNAs, as would be expected if they serve as templates and targets for these siRNAs. The other 75% tended to be more conserved and included lincRNAs with intriguing expression and sequence features associating them with processes such as dauer formation, male identity, sperm formation, and interaction with sperm-specific mRNAs. Our study provides a glimpse into the lncRNA content of a nonvertebrate animal and a resource for future studies of lncRNA function.

Coexpression of linked gene pairs persists long after their separation

In many eukaryotes, physically linked gene pairs tend to be coexpressed. However, it is still controversial to what extent this neighbor coexpression is maintained by selection and to what extent it is nonselective, purely mechanistic "leaky expression." Here, we analyze expression patterns of gene pairs that have lost their linkage in the evolution of Saccharomyces cerevisiae since its last common ancestor with Kluyveromyces waltii or that were never linked in the S. cerevisiae lineage but became neighbors in a related yeast. We demonstrate that coexpression of many linked genes is retained long after their separation and is thus likely to be functionally important. In addition, unlinked gene pairs that recently became neighbors in other yeast species tend to be coexpressed in S. cerevisiae. This suggests that natural selection often favors chromosomal rearrangements in which coexpressed genes become neighbors. Contrary to previous suggestions, selectively favorable coexpression appears not to be restricted to bidirectional promoters.

Minimal regulatory spaces in yeast genomes

BACKGROUND

The regulatory information encoded in the DNA of promoter regions usually enforces a minimal, non-zero distance between the coding regions of neighboring genes. However, the size of this minimal regulatory space is not generally known. In particular, it is unclear if minimal promoter size differs between species and between uni- and bi-directionally acting regulatory regions.

RESULTS

Analyzing the genomes of 11 yeasts, we show that the lower size limit on promoter-containing regions is species-specific within a relatively narrow range (80-255 bp). This size limit applies equally to regions that initiate transcription on one or both strands, indicating that bi-directional promoters and uni-directional promoters are constrained similarly. We further find that young, species-specific regions are on average much longer than older regions, suggesting either a bias towards deletions or selection for genome compactness in yeasts. While the length evolution of promoter-less intergenic regions is well described by a simplistic, purely neutral model, regions containing promoters typically show an excess of unusually long regions. Regions flanked by divergently transcribed genes have a bi-modal length distribution, with short lengths found preferentially among older regions. These old, short regions likely harbor evolutionarily conserved bi-directionally active promoters. Surprisingly, some of the evolutionarily youngest regions in two of the eleven species (S. cerevisiae and K. waltii) are shorter than the lower limit observed in older regions.

CONCLUSIONS

The minimal chromosomal space required for transcriptional regulation appears to be relatively similar across yeast species, and is the same for uni-directional and bi-directional promoters. New intergenic regions created by genome rearrangements tend to evolve towards the more narrow size distribution found among older regions.

Transcriptional coupling of neighboring genes and gene expression noise: evidence that gene orientation and noncoding transcripts are modulators of noise

How is noise in gene expression modulated? Do mechanisms of noise control impact genome organization? In yeast, the expression of one gene can affect that of a very close neighbor. As the effect is highly regionalized, we hypothesize that genes in different orientations will have differing degrees of coupled expression and, in turn, different noise levels. Divergently organized gene pairs, in particular those with bidirectional promoters, have close promoters, maximizing the likelihood that expression of one gene affects the neighbor. With more distant promoters, the same is less likely to hold for gene pairs in nondivergent orientation. Stochastic models suggest that coupled chromatin dynamics will typically result in low abundance-corrected noise (ACN). Transcription of noncoding RNA (ncRNA) from a bidirectional promoter, we thus hypothesize to be a noise-reduction, expression-priming, mechanism. The hypothesis correctly predicts that protein-coding genes with a bidirectional promoter, including those with a ncRNA partner, have lower ACN than other genes and divergent gene pairs uniquely have correlated ACN. Moreover, as predicted, ACN increases with the distance between promoters. The model also correctly predicts ncRNA transcripts to be often divergently transcribed from genes that a priori would be under selection for low noise (essential genes, protein complex genes) and that the latter genes should commonly reside in divergent orientation. Likewise, that genes with bidirectional promoters are rare subtelomerically, cluster together, and are enriched in essential gene clusters is expected and observed. We conclude that gene orientation and transcription of ncRNAs are candidate modulators of noise.

Conservation, rearrangement, and deletion of gene pairs during the evolution of four grass genomes

Gene order and content differ among homologous regions of closely related genomes. Similarities in the expression profiles of physically adjacent genes suggest that the proper functioning of these genes depends on maintaining a specific position relative to each other. To better understand the results of the interaction of these two genomic forces, convergent, divergent, and tandem gene pairs in rice and sorghum, as well as their homologs in rice, sorghum, maize, and Brachypodium were analyzed. The status of each pair in all four species: whether it was conserved, inverted, rearranged, or missing homologs was determined. We observed that divergent gene pairs had lower rates of conservation than convergent or tandem pairs, but higher rates of rearranged pairs and missing homologs in maize than in any other species. We also discovered species-specific gene pairs in rice and sorghum. In rice, gene pairs with strongly correlated expression levels were conserved significantly more often than those with little or no correlation. We assigned three types of gene pair to one of 14 possible evolutionary history categories to uncover their evolutionary dynamics during the evolution of grass genomes.

Genome-wide analysis of the cis-regulatory modules of divergent gene pairs in yeast

In budding yeast, approximately a quarter of adjacent genes are divergently transcribed (divergent gene pairs). Whether genes in a divergent pair share the same regulatory system is still unknown. By examining transcription factor (TF) knockout experiments, we found that most TF knockout only altered the expression of one gene in a divergent pair. This prompted us to conduct a comprehensive analysis in silico to estimate how many divergent pairs are regulated by common sets of TFs (cis-regulatory modules, CRMs) using TF binding sites and expression data. Analyses of ten expression datasets show that only a limited number of divergent gene pairs share CRMs in any single dataset. However, around half of divergent pairs do share a regulatory system in at least one dataset. Our analysis suggests that genes in a divergent pair tend to be co-regulated in at least one condition; however, in most conditions, they may not be co-regulated.

ATP-citrate lyase is required for production of cytosolic acetyl coenzyme A and development in Aspergillus nidulans

Acetyl coenzyme A (CoA) is a central metabolite in carbon and energy metabolism and in the biosynthesis of cellular molecules. A source of cytoplasmic acetyl-CoA is essential for the production of fatty acids and sterols and for protein acetylation, including histone acetylation in the nucleus. In Saccharomyces cerevisiae and Candida albicans acetyl-CoA is produced from acetate by cytoplasmic acetyl-CoA synthetase, while in plants and animals acetyl-CoA is derived from citrate via ATP-citrate lyase. In the filamentous ascomycete Aspergillus nidulans, tandem divergently transcribed genes (aclA and aclB) encode the subunits of ATP-citrate lyase, and we have deleted these genes. Growth is greatly diminished on carbon sources that do not result in cytoplasmic acetyl-CoA, such as glucose and proline, while growth is not affected on carbon sources that result in the production of cytoplasmic acetyl-CoA, such as acetate and ethanol. Addition of acetate restores growth on glucose or proline, and this is dependent on facA, which encodes cytoplasmic acetyl-CoA synthetase, but not on the regulatory gene facB. Transcription of aclA and aclB is repressed by growth on acetate or ethanol. Loss of ATP-citrate lyase results in severe developmental effects, with the production of asexual spores (conidia) being greatly reduced and a complete absence of sexual development. This is in contrast to Sordaria macrospora, in which fruiting body formation is initiated but maturation is defective in an ATP-citrate lyase mutant. Addition of acetate does not repair these defects, indicating a specific requirement for high levels of cytoplasmic acetyl-CoA during differentiation. Complementation in heterokaryons between aclA and aclB deletions for all phenotypes indicates that the tandem gene arrangement is not essential.