Genome-wide computational prediction of tandem gene arrays: application in yeasts

Colocality to Cofunctionality: Eukaryotic Gene Neighborhoods as a Resource for Function Discovery

Diverging from the classic paradigm of random gene order in eukaryotes, gene proximity can be leveraged to systematically identify functionally related gene neighborhoods in eukaryotes, utilizing techniques pioneered in bacteria. Current methods of identifying gene neighborhoods typically rely on sequence similarity to characterized gene products. However, this approach is not robust for nonmodel organisms like algae, which are evolutionarily distant from well-characterized model organisms. Here, we utilize a comparative genomic approach to identify evolutionarily conserved proximal orthologous gene pairs conserved across at least two taxonomic classes of green algae. A total of 317 gene neighborhoods were identified. In some cases, gene proximity appears to have been conserved since before the streptophyte–chlorophyte split, 1,000 Ma. Using functional inferences derived from reconstructed evolutionary relationships, we identified several novel functional clusters. A putative mycosporine-like amino acid, “sunscreen,” neighborhood contains genes similar to either vertebrate or cyanobacterial pathways, suggesting a novel mosaic biosynthetic pathway in green algae. One of two putative arsenic-detoxification neighborhoods includes an organoarsenical transporter (ArsJ), a glyceraldehyde 3-phosphate dehydrogenase-like gene, homologs of which are involved in arsenic detoxification in bacteria, and a novel algal-specific phosphoglycerate kinase-like gene. Mutants of the ArsJ-like transporter and phosphoglycerate kinase-like genes in Chlamydomonas reinhardtii were found to be sensitive to arsenate, providing experimental support for the role of these identified neighbors in resistance to arsenate. Potential evolutionary origins of neighborhoods are discussed, and updated annotations for formerly poorly annotated genes are presented, highlighting the potential of this strategy for functional annotation.

Expression attenuation as a mechanism of robustness against gene duplication

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein-protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Life cycle adapted upstream open reading frames (uORFs) in Trypanosoma congolense: A post-transcriptional approach to accurate gene regulation

The presented work explores the regulatory influence of upstream open reading frames (uORFs) on gene expression in Trypanosoma congolense. More than 31,000 uORFs in total were identified and characterized here. We found evidence for the uORFs’ appearance in the transcriptome to be correlated with proteomic expression data, clearly indicating their repressive potential in T. congolense, which has to rely on post-transcriptional gene expression regulation due to its unique genomic organization. Our data show that uORF’s translation repressive potential does not only correlate with elemental sequence features such as length, position and quantity, but involves more subtle components, in particular the codon and amino acid profiles. This corresponds with the popular mechanistic model of a ribosome shedding initiation factors during the translation of a uORF, which can prevent reinitiation at the downstream start codon of the actual protein-coding sequence, due to the former extensive consumption of crucial translation components. We suggest that uORFs with uncommon codon and amino acid usage can slow down the translation elongation process in T. congolense, systematically deplete the limited factors, and restrict downstream reinitiation, setting up a bottleneck for subsequent translation of the protein-coding sequence. Additionally we conclude that uORFs dynamically influence the T. congolense life cycle. We found evidence that transition to epimastigote form could be supported by gain of uORFs due to alternative trans-splicing, which down-regulate housekeeping genes’ expression and render the trypanosome in a metabolically reduced state of endurance.

Genetic basis of priority effects: insights from nectar yeast

Priority effects, in which the order of species arrival dictates community assembly, can have a major influence on species diversity, but the genetic basis of priority effects remains unknown. Here, we suggest that nitrogen scavenging genes previously considered responsible for starvation avoidance may drive priority effects by causing rapid resource depletion. Using single-molecule sequencing, we de novo assembled the genome of the nectar-colonizing yeast, Metschnikowia reukaufii, across eight scaffolds and complete mitochondrion, with gap-free coverage over gene spaces. We found a high rate of tandem gene duplication in this genome, enriched for nitrogen metabolism and transport. Both high-capacity amino acid importers, GAP1 and PUT4, present as tandem gene arrays, were highly expressed in synthetic nectar and regulated by the availability and quality of amino acids. In experiments with competitive nectar yeast, Candida rancensis, amino acid addition alleviated suppression of C. rancensis by early arrival of M. reukaufii, corroborating that amino acid scavenging may contribute to priority effects. Because niche pre-emption via rapid resource depletion may underlie priority effects in a broad range of microbial, plant and animal communities, nutrient scavenging genes like the ones we considered here may be broadly relevant to understanding priority effects.

Massive Amplification at an Unselected Locus Accompanies Complex Chromosomal Rearrangements in Yeast

Gene amplification has been observed in different organisms in response to environmental constraints, such as limited nutrients or exposure to a variety of toxic compounds, conferring them with specific phenotypic adaptations via increased expression levels. However, the presence of multiple gene copies in natural genomes has generally not been found in the absence of specific functional selection. Here, we show that the massive amplification of a chromosomal locus (up to 880 copies per cell) occurs in the absence of any direct selection, and is associated with low-order amplifications of flanking segments in complex chromosomal alterations. These results were obtained from mutants with restored phenotypes that spontaneously appeared from genetically engineered strains of the yeast Saccharomyces cerevisiae suffering from severe fitness reduction. Grossly extended chromosomes (macrotene) were formed, with complex structural alterations but sufficient stability to propagate unchanged over successive generations. Their detailed molecular analysis, including complete genome sequencing, identification of sequence breakpoints, and comparisons between mutants, revealed novel mechanisms causing their formation, whose combined action underlies the astonishing dynamics of eukaryotic chromosomes and their consequences.

Rapid evolution of recombinant Saccharomyces cerevisiae for Xylose fermentation through formation of extra-chromosomal circular DNA

Circular DNA elements are involved in genome plasticity, particularly of tandem repeats. However, amplifications of DNA segments in Saccharomyces cerevisiae reported so far involve pre-existing repetitive sequences such as ribosomal DNA, Ty elements and Long Terminal Repeats (LTRs). Here, we report the generation of an eccDNA, (extrachromosomal circular DNA element) in a region without any repetitive sequences during an adaptive evolution experiment. We performed whole genome sequence comparison between an efficient D-xylose fermenting yeast strain developed by metabolic and evolutionary engineering, and its parent industrial strain. We found that the heterologous gene XylA that had been inserted close to an ARS sequence in the parent strain has been amplified about 9 fold in both alleles of the chromosomal locus of the evolved strain compared to its parent. Analysis of the amplification process during the adaptive evolution revealed formation of a XylA-carrying eccDNA, pXI2-6, followed by chromosomal integration in tandem arrays over the course of the evolutionary adaptation. Formation of the eccDNA occurred in the absence of any repetitive DNA elements, probably using a micro-homology sequence of 8 nucleotides flanking the amplified sequence. We isolated the pXI2-6 eccDNA from an intermediate strain of the evolutionary adaptation process, sequenced it completely and showed that it confers high xylose fermentation capacity when it is transferred to a new strain. In this way, we have provided clear evidence that gene amplification can occur through generation of eccDNA without the presence of flanking repetitive sequences and can serve as a rapid means of adaptation to selection pressure.

Xylose is an important component of lignocellulose hydrolysates used for the production of bioethanol, but the yeast Saccharomyces cerevisiae is unable to utilize xylose. Insertion of a bacterial xylose isomerase gene and improvement of growth on xylose by evolutionary adaptation resulted in amplification of this gene and efficient xylose fermentation capacity. Further analysis of the final and intermediate strains from the evolutionary adaptation process revealed interesting features about the mechanisms involved in gene amplification events, which have occurred frequently in natural evolution. We now show that a circular DNA element was spontaneously created by the yeast, encompassing the xylose isomerase gene and an ARS element, present by coincidence adjacent of the inserted xylose isomerase gene. ARS elements are the sites where DNA polymerase initiates duplication of DNA. Interestingly, this has revealed for the first time in yeast that circular DNA plasmids can be created from genomic DNA in the absence of flanking repetitive sequences.

Detecting long tandem duplications in genomic sequences

BACKGROUND

Detecting duplication segments within completely sequenced genomes provides valuable information to address genome evolution and in particular the important question of the emergence of novel functions. The usual approach to gene duplication detection, based on all-pairs protein gene comparisons, provides only a restricted view of duplication.

RESULTS

In this paper, we introduce ReD Tandem, a software using a flow based chaining algorithm targeted at detecting tandem duplication arrays of moderate to longer length regions, with possibly locally weak similarities, directly at the DNA level. On the A. thaliana genome, using a reference set of tandem duplicated genes built using TAIR,(a) we show that ReD Tandem is able to predict a large fraction of recently duplicated genes (dS  <  1) and that it is also able to predict tandem duplications involving non coding elements such as pseudo-genes or RNA genes.

CONCLUSIONS

ReD Tandem allows to identify large tandem duplications without any annotation, leading to agnostic identification of tandem duplications. This approach nicely complements the usual protein gene based which ignores duplications involving non coding regions. It is however inherently restricted to relatively recent duplications. By recovering otherwise ignored events, ReD Tandem gives a more comprehensive view of existing evolutionary processes and may also allow to improve existing annotations.

Clusters of ancestrally related genes that show paralogy in whole or in part are a major feature of the genomes of humans and other species

Arrangements of genes along chromosomes are a product of evolutionary processes, and we can expect that preferable arrangements will prevail over the span of evolutionary time, often being reflected in the non-random clustering of structurally and/or functionally related genes. Such non-random arrangements can arise by two distinct evolutionary processes: duplications of DNA sequences that give rise to clusters of genes sharing both sequence similarity and common sequence features and the migration together of genes related by function, but not by common descent [1], [2], [3]. To provide a background for distinguishing between the two, which is important for future efforts to unravel the evolutionary processes involved, we here provide a description of the extent to which ancestrally related genes are found in proximity.

Towards this purpose, we combined information from five genomic datasets, InterPro, SCOP, PANTHER, Ensembl protein families, and Ensembl gene paralogs. The results are provided in publicly available datasets (http://cgd.jax.org/datasets/clustering/paraclustering.shtml) describing the extent to which ancestrally related genes are in proximity beyond what is expected by chance (i.e. form paraclusters) in the human and nine other vertebrate genomes, as well as the D. melanogaster, C. elegans, A. thaliana, and S. cerevisiae genomes.

With the exception of Saccharomyces, paraclusters are a common feature of the genomes we examined. In the human genome they are estimated to include at least 22% of all protein coding genes. Paraclusters are far more prevalent among some gene families than others, are highly species or clade specific and can evolve rapidly, sometimes in response to environmental cues. Altogether, they account for a large portion of the functional clustering previously reported in several genomes.