Search

Abstract Background This paper describes an efficient in silico method for detecting tandem gene arrays (TGAs) in fully sequenced and compact genomes such as those of prokaryotes or unicellular eukaryotes. The originality of this method lies in the search of protein sequence similarities in the vicinity of each coding sequence, which allows the prediction of tandem duplicated gene copies independently of their functionality. Results Applied to nine hemiascomycete yeast genomes, this method predicts that 2% of the genes are involved in TGAs and gene relics are present in 11% of TGAs. The frequency of TGAs with degenerated gene copies means that a significant fraction of tandem duplicated genes follows the birth-and-death model of evolution. A comparison of sequence identity distributions between sets of homologous gene pairs shows that the different copies of tandem arrayed paralogs are less divergent than copies of dispersed paralogs in yeast genomes. It suggests that paralogs included in tandem structures are more recent or more subject to the gene conversion mechanism than other paralogs. Conclusion The method reported here is a useful computational tool to provide a database of TGAs composed of functional or nonfunctional gene copies. Such a database has obvious applications in the fields of structural and comparative genomics. Notably, a detailed study of the TGA catalog will make it possible to tackle the fundamental questions of the origin and evolution of tandem gene clusters.

tRNAs are not only essential for decoding the genetic code, but their abundance also has a strong impact on the rate of protein production, folding, and on the stability of the translated messenger RNAs. Plasmodium expresses a unique surface protein called tRip, involved in the import of exogenous tRNAs into the parasite. Comparative proteomic analysis of the blood stage of wild-type and tRip-KO variant of P. berghei parasites revealed that downregulated proteins in the mutant parasite are distinguished by a bias in their asparagine content. Furthermore, the demonstration of the possibility of charging host tRNAs with Plasmodium aminoacyl-tRNA synthetases led us to propose that imported host tRNAs participate in parasite protein synthesis. These results also suggest a novel mechanism of translational control in which import of host tRNAs emerge as regulators of gene expression in the Plasmodium developmental cycle and pathogenesis, by enabling the synthesis of asparagine-rich regulatory proteins that efficiently and selectively control the parasite infectivity., (© 2024 International Union of Biochemistry and Molecular Biology.)

In eukaryotes, translation initiation is a highly regulated process, which combines cis- regulatory sequences located on the messenger RNA along with trans- acting factors like eukaryotic initiation factors (eIF). One critical step of translation initiation is the start codon recognition by the scanning 43S particle, which leads to ribosome assembly and protein synthesis. In this study, we investigated the involvement of secondary structures downstream the initiation codon in the so-called START (STructure-Assisted RNA translation) mechanism on AUG and non-AUG translation initiation. The results demonstrate that downstream secondary structures can efficiently promote non-AUG translation initiation if they are sufficiently stable to stall a scanning 43S particle and if they are located at an optimal distance from non-AUG codons to stabilize the codon-anticodon base pairing in the P site. The required stability of the downstream structure for efficient translation initiation varies in distinct cell types. We extended this study to genome-wide analysis of functionally characterized alternative translation initiation sites in Homo sapiens . This analysis revealed that about 25% of these sites have an optimally located downstream secondary structure of adequate stability which could elicit START, regardless of the start codon. We validated the impact of these structures on translation initiation for several selected uORFs., (© The Author(s) 2024. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

Translation initiation is a key step in the protein synthesis stage of the gene expression pathway of all living cells. In this important process, ribosomes have to accurately find the AUG start codon in order to ensure the integrity of the proteome. "Structure Assisted RNA Translation", or "START", has been proposed to use stable secondary structures located in the coding sequence to augment start site selection by steric hindrance of the progression of pre-initiation complex on messenger RNA. This implies that such structures have to be located downstream and at on optimal distance from the AUG start codon (i.e., downstream nucleotide +16). In order to assess the importance of the START mechanism in the overall mRNA translation process, we developed a bioinformatic tool to screen coding sequences for such stable structures in a 50 nucleotide-long window spanning the nucleotides from +16 to +65. We screened eight bacterial genomes and six eukaryotic genomes. We found stable structures in 0.6-2.5% of eukaryotic coding sequences. Among these, approximately half of them were structures predicted to form G-quadruplex structures. In humans, we selected 747 structures. In bacteria, the coding sequences from Gram-positive bacteria contained 2.6-4.2% stable structures, whereas the structures were less abundant in Gram-negative bacteria (0.2-2.7%). In contrast to eukaryotes, putative G-quadruplex structures are very rare in the coding sequence of bacteria. Altogether, our study reveals that the START mechanism seems to be an ancient strategy to facilitate the start codon recognition that is used in different kingdoms of life.

Decoding of the 61 sense codons of the genetic code requires a variable number of tRNAs that establish codon-anticodon interactions. Thanks to the wobble base pairing at the third codon position, less than 61 different tRNA isoacceptors are needed to decode the whole set of codons. On the tRNA, a subtle distribution of nucleoside modifications shapes the anticodon loop structure and participates to accurate decoding and reading frame maintenance. Interestingly, although the 61 anticodons should exist in tRNAs, a strict absence of some tRNAs decoders is found in several codon families. For instance, in Eukaryotes, G34-containing tRNAs translating 3-, 4- and 6-codon boxes are absent. This includes tRNA specific for Ala, Arg, Ile, Leu, Pro, Ser, Thr, and Val. tRNA Gly is the only exception for which in the three kingdoms, a G34-containing tRNA exists to decode C3 and U3-ending codons. To understand why G34-tRNA Gly exists, we analysed at the genome wide level the codon distribution in codon +1 relative to the four GGN Gly codons. When considering codon GGU, a bias was found towards an unusual high usage of codons starting with a G whatever the amino acid at +1 codon. It is expected that GGU codons are decoded by G34-containing tRNA Gly , decoding also GGC codons. Translation studies revealed that the presence of a G at the first position of the downstream codon reduces the +1 frameshift by stabilizing the G34•U3 wobble interaction. This result partially explains why G34-containing tRNA Gly exists in Eukaryotes whereas all the other G34-containing tRNAs for multiple codon boxes are absent., Competing Interests: The authors declare no conflict of interest.

Non-coding RNA regulatory elements are important for viral replication, making them promising targets for therapeutic intervention. However, regulatory RNA is challenging to detect and characterise using classical structure-function assays. Here, we present in cell Mutational Interference Mapping Experiment (in cell MIME) as a way to define RNA regulatory landscapes at single nucleotide resolution under native conditions. In cell MIME is based on (i) random mutation of an RNA target, (ii) expression of mutated RNA in cells, (iii) physical separation of RNA into functional and non-functional populations, and (iv) high-throughput sequencing to identify mutations affecting function. We used in cell MIME to define RNA elements within the 5' region of the HIV-1 genomic RNA (gRNA) that are important for viral replication in cells. We identified three distinct RNA motifs controlling intracellular gRNA production, and two distinct motifs required for gRNA packaging into virions. Our analysis reveals the 73AAUAAA78 polyadenylation motif within the 5' PolyA domain as a dual regulator of gRNA production and gRNA packaging, and demonstrates that a functional polyadenylation signal is required for viral packaging even though it negatively affects gRNA production.

RNA regulates many biological processes; however, identifying functional RNA sequences and structures is complex and time-consuming. We introduce a method, mutational interference mapping experiment (MIME), to identify, at single-nucleotide resolution, the primary sequence and secondary structures of an RNA molecule that are crucial for its function. MIME is based on random mutagenesis of the RNA target followed by functional selection and next-generation sequencing. Our analytical approach allows the recovery of quantitative binding parameters and permits the identification of base-pairing partners directly from the sequencing data. We used this method to map the binding site of the human immunodeficiency virus-1 (HIV-1) Pr55(Gag) protein on the viral genomic RNA in vitro, and showed that, by analyzing permitted base-pairing patterns, we could model RNA structure motifs that are crucial for protein binding.

Background: The industrially important yeast Blastobotrys (Arxula) adeninivorans is an asexual hemiascomycete phylogenetically very distant from Saccharomyces cerevisiae. Its unusual metabolic flexibility allows it to use a wide range of carbon and nitrogen sources, while being thermotolerant, xerotolerant and osmotolerant., Results: The sequencing of strain LS3 revealed that the nuclear genome of A. adeninivorans is 11.8 Mb long and consists of four chromosomes with regional centromeres. Its closest sequenced relative is Yarrowia lipolytica, although mean conservation of orthologs is low. With 914 introns within 6116 genes, A. adeninivorans is one of the most intron-rich hemiascomycetes sequenced to date. Several large species-specific families appear to result from multiple rounds of segmental duplications of tandem gene arrays, a novel mechanism not yet described in yeasts. An analysis of the genome and its transcriptome revealed enzymes with biotechnological potential, such as two extracellular tannases (Atan1p and Atan2p) of the tannic-acid catabolic route, and a new pathway for the assimilation of n-butanol via butyric aldehyde and butyric acid., Conclusions: The high-quality genome of this species that diverged early in Saccharomycotina will allow further fundamental studies on comparative genomics, evolution and phylogenetics. Protein components of different pathways for carbon and nitrogen source utilization were identified, which so far has remained unexplored in yeast, offering clues for further biotechnological developments. In the course of identifying alternative microorganisms for biotechnological interest, A. adeninivorans has already proved its strengthened competitiveness as a promising cell factory for many more applications.

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993(T)), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ~13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

Polyploidization is an important process in the evolution of eukaryotic genomes, but ensuing molecular mechanisms remain to be clarified. Autopolyploidization or whole-genome duplication events frequently are resolved in resulting lineages by the loss of single genes from most duplicated pairs, causing transient gene dosage imbalance and accelerating speciation through meiotic infertility. Allopolyploidization or formation of interspecies hybrids raises the problem of genetic incompatibility (Bateson-Dobzhansky-Muller effect) and may be resolved by the accumulation of mutational changes in resulting lineages. In this article, we show that an osmotolerant yeast species, Pichia sorbitophila, recently isolated in a concentrated sorbitol solution in industry, illustrates this last situation. Its genome is a mosaic of homologous and homeologous chromosomes, or parts thereof, that corresponds to a recently formed hybrid in the process of evolution. The respective parental contributions to this genome were characterized using existing variations in GC content. The genomic changes that occurred during the short period since hybrid formation were identified (e.g., loss of heterozygosity, unilateral loss of rDNA, reciprocal exchange) and distinguished from those undergone by the two parental genomes after separation from their common ancestor (i.e., NUMT (NUclear sequences of MiTochondrial origin) insertions, gene acquisitions, gene location movements, reciprocal translocation). We found that the physiological characteristics of this new yeast species are determined by specific but unequal contributions of its two parents, one of which could be identified as very closely related to an extant Pichia farinosa strain.

This short article presents an overview of tandem gene arrays (TGAs) in hemiascomycete yeasts. In silico and in vivo analyses are combined to address structural, functional and evolutionary aspects of these particular chromosomal structures. Genomic instability of TGAs is discussed. We conclude that TGAs are generally dynamic regions of the genome in that they are the seats of chromosomal rearrangement events. In addition, they are often breeding grounds of new genes for a rapid adaptation of cells to demands of the environment., (Copyright © 2011 Académie des sciences. Published by Elsevier SAS. All rights reserved.)

Background: This paper describes an efficient in silico method for detecting tandem gene arrays (TGAs) in fully sequenced and compact genomes such as those of prokaryotes or unicellular eukaryotes. The originality of this method lies in the search of protein sequence similarities in the vicinity of each coding sequence, which allows the prediction of tandem duplicated gene copies independently of their functionality., Results: Applied to nine hemiascomycete yeast genomes, this method predicts that 2% of the genes are involved in TGAs and gene relics are present in 11% of TGAs. The frequency of TGAs with degenerated gene copies means that a significant fraction of tandem duplicated genes follows the birth-and-death model of evolution. A comparison of sequence identity distributions between sets of homologous gene pairs shows that the different copies of tandem arrayed paralogs are less divergent than copies of dispersed paralogs in yeast genomes. It suggests that paralogs included in tandem structures are more recent or more subject to the gene conversion mechanism than other paralogs., Conclusion: The method reported here is a useful computational tool to provide a database of TGAs composed of functional or nonfunctional gene copies. Such a database has obvious applications in the fields of structural and comparative genomics. Notably, a detailed study of the TGA catalog will make it possible to tackle the fundamental questions of the origin and evolution of tandem gene clusters.

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

The DUP gene family of Saccharomyces cerevisiae comprises 23 members that can be divided into two subfamilies--DUP240 and DUP380. The location of the DUP loci suggests that at least three mechanisms were responsible for their genomic dispersion: nonreciprocal translocation at chromosomal ends, tandem duplication and Ty-associated duplication. The data we present here suggest that these nonessential genes encode proteins that facilitate membrane trafficking processes. Dup240 proteins have three conserved domains (C1, C2 and C3) and two predicted transmembrane segments (H1 and H2). A direct repetition of the C1-H1-H2-C2 module is observed in Dup380p sequences. In this article, we propose an evolutionary model to account for the emergence of the two gene subfamilies.

Duplication, resulting in gene redundancy, is well known to be a driving force of evolutionary change. Gene families are therefore useful targets for approaching genome evolution. To address the gene death process, we examined the fate of the 10-member-large S288C DUP240 family in 15 Saccharomyces cerevisiae strains. Using an original three-step method of analysis reported here, both slightly and highly degenerate DUP240 copies, called pseudo-open reading frames (ORFs) and relics, respectively, were detected in strain S288C. It was concluded that two previously annotated ORFs correspond, in fact, to pseudo-ORFs and three additional relics were identified in intergenic areas. Comparative intraspecies analysis of these degenerate DUP240 loci revealed that the two pseudo-ORFs are present in a nondegenerate state in some other strains. This suggests that within a given gene family different loci are the target of the gene erasure process, which is therefore strain dependent. Besides, the variable positions observed indicate that the relic sequence may diverge faster than the flanking regions. All in all, this study shows that short conserved protein motifs provide a useful tool for detecting and accurately mapping degenerate gene remnants. The present results also highlight the strong contribution of comparative genomics for gene relic detection because the possibility of finding short conserved protein motifs in intergenic regions (IRs) largely depends on the choice of the most closely related paralog or ortholog. By mapping new genetic components in previously annotated IRs, our study constitutes a further refinement step in the crucial stage of genome annotation and provides a strategy for retracing ancient chromosomal reshaping events and, hence, for deciphering genome history.

The influence of duplicated sequences on chromosomal stability is poorly understood. To characterize chromosomal rearrangements involving duplicated sequences, we compared the organization of tandem repeats of the DUP240 gene family in 15 Saccharomyces cerevisiae strains of various origins. The DUP240 gene family consists of 10 members of unknown function in the reference strain S288C. Five DUP240 paralogs on chromosome I and two on chromosome VII are arranged as tandem repeats that are highly polymorphic in copy number and sequence. We characterized DNA sequences that are likely involved in homologous or nonhomologous recombination events and are responsible for intra- and interchromosomal rearrangements that cause the creation and disappearance of DUP240 paralogs. The tandemly repeated DUP240 genes seem to be privileged sites of gene birth and death., (Copyright 2004 Genetics Society of America)

Identifying the mechanisms of eukaryotic genome evolution by comparative genomics is often complicated by the multiplicity of events that have taken place throughout the history of individual lineages, leaving only distorted and superimposed traces in the genome of each living organism. The hemiascomycete yeasts, with their compact genomes, similar lifestyle and distinct sexual and physiological properties, provide a unique opportunity to explore such mechanisms. We present here the complete, assembled genome sequences of four yeast species, selected to represent a broad evolutionary range within a single eukaryotic phylum, that after analysis proved to be molecularly as diverse as the entire phylum of chordates. A total of approximately 24,200 novel genes were identified, the translation products of which were classified together with Saccharomyces cerevisiae proteins into about 4,700 families, forming the basis for interspecific comparisons. Analysis of chromosome maps and genome redundancies reveal that the different yeast lineages have evolved through a marked interplay between several distinct molecular mechanisms, including tandem gene repeat formation, segmental duplication, a massive genome duplication and extensive gene loss.

The DUP240 gene family of Saccharomyces cerevisiae is composed of 10 members. They encode proteins of about 240 amino acids which contain two predicted transmembrane domains. Database searches identified only one homologue in the closely related species Saccharomyces bayanus, indicating that the DUP240 genes encode proteins specific to Saccharomyces sensu stricto. The short-flanking homology PCR gene-replacement strategy with a variety of selective markers for replacements, and classical genetic methods, were used to generate strains deleted for all 10 DUP240 genes. All of the knock-out strains were viable and had similar growth kinetics to the wild-type. Two-hybrid screens, hSos1p fusions and GFP fusions were carried out; the results indicated that the Dup240 proteins are membrane associated, and that some of them are concentrated around the plasma membrane.

The yeast Saccharomyces cerevisiae possesses two methionyl-tRNA synthetases (MetRS), one in the cytoplasm and the other in mitochondria. The cytoplasmic MetRS has a zinc-finger motif of the type Cys-X(2)-Cys-X(9)-Cys-X(2)-Cys in an insertion domain that divides the nucleotide-binding fold into two halves, whereas no such motif is present in the mitochondrial MetRS. Here, we show that tightly bound zinc atom is present in the cytoplasmic MetRS but not in the mitochondrial MetRS. To test whether the presence of a zinc-binding site is required for cytoplasmic functions of MetRS, we constructed a yeast strain in which cytoplasmic MetRS gene was inactivated and the mitochondrial MetRS gene was expressed in the cytoplasm. Provided that methionine-accepting tRNA is overexpressed, this strain was viable, indicating that mitochondrial MetRS was able to aminoacylate tRNA(Met) in the cytoplasm. Site-directed mutagenesis demonstrated that the zinc domain was required for the stability and consequently for the activity of cytoplasmic MetRS. Mitochondrial MetRS, like cytoplasmic MetRS, supported homocysteine editing in vivo in the yeast cytoplasm. Both MetRSs catalyzed homocysteine editing and aminoacylation of coenzyme A in vitro. Thus, identical synthetic and editing functions can be carried out in different structural frameworks of cytoplasmic and mitochondrial MetRSs., (Copyright 2001 Academic Press.)

As for Escherichia coli methionine tRNAs, the anticodon triplet of yeast tRNA(Met) plays an important role in the recognition by the yeast methionyl-tRNA synthetase (MetRS), indicating that this determinant for methionine identity is conserved in yeast. Efficient aminoacylation of the E. coli tRNA(Met) transcript by the heterologous yeast methionine enzyme also suggests conservation of the protein determinants that interact with the CAU anticodon sequence. We have analysed by site-directed mutagenesis the peptide region 655 to 663 of the yeast MetRS that is equivalent to the anticodon binding region of the E. coli methionine enzyme. Only one change, converting Leu658 into Ala significantly reduced tRNA aminoacylation. Semi-conservative substitutions of L658 allow a correlation to be drawn between side-chain volume of the hydrophobic residue at this site and activity. The analysis of the L658A mutant shows that Km is mainly affected. This suggests that the peptide region 655 to 663 contributes partially to the binding of the anticodon, since separate mutational analysis of the anticodon bases shows that kcat is the most critical parameter in the recognition of tRNA(Met) by the yeast synthetase. We have analysed the role of peptide region (583-GNLVNR-588) that is spatially close to the region 655 to 663. Replacements of residues N584 and R588 reduces significantly the kcat of aminoacylation. The peptide region 583-GNLVNR-588 is highly conserved in all MetRS so far sequenced. We therefore propose that the hydrogen donor/acceptor amino acid residues within this region are the most critical protein determinants for the positive selection of the methionine tRNAs.

Sequence comparisons among methionyl-tRNA synthetases from different organisms reveal only one block of homology beyond the last beta strand of the mononucleotide fold. We have introduced a series of semi-conservative amino acid replacements in the conserved motif of yeast methionyl-tRNA synthetase. The results indicate that replacements of two polar residues (Asn584 and Arg588) affected specifically the aminoacylation reaction. The location of these residues in the tertiary structure of the enzyme is compatible with a direct interaction of the amino acid side-chains with the tRNA anticodon.