Overlapping coding regions and trancriptional units of two essential chromosomal genes (CCT8, TRP1)in the fungal pathogen Candida albicans

HexSE: Simulating evolution in overlapping reading frames

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may provide a mechanism to increase the information content of compact genomes. The presence of overlapping reading frames (OvRFs) can skew estimates of selection based on the rates of non-synonymous and synonymous substitutions, since a substitution that is synonymous in one reading frame may be non-synonymous in another and vice versa. To understand the impact of OvRFs on molecular evolution, we implemented a versatile simulation model of nucleotide sequence evolution along a phylogeny with any distribution of open reading frames in linear or circular genomes. We use a custom data structure to track the substitution rates at every nucleotide site, which is determined by the stationary nucleotide frequencies, transition bias and the distribution of selection biases (dN/dS) in the respective reading frames. Our simulation model is implemented in the Python scripting language. All source code is released under the GNU General Public License version 3 and are available at https://github.com/PoonLab/HexSE.

Overlapping genes: a new strategy of thermophilic stress tolerance in prokaryotes

Overlapping genes (OGs) draw the focus of recent day's research. However, the significance of OGs in prokaryotic genomes remained unexplored. As an adaptation to high temperature, thermophiles were shown to eliminate their intergenic regions. Therefore, it could be possible that prokaryotes would increase their OG content to adapt to high temperature. To test this hypothesis, we carried out a comparative study on OG frequency of 256 prokaryotic genomes comprising both thermophiles and non-thermophiles. It was found that thermophiles exhibit higher frequency of overlapping genes than non-thermophiles. Moreover, overlap frequency was found to correlate with optimal growth temperature (OGT) in prokaryotes. Long overlap frequency was found to hold a positive correlation with OGT resulting in an abundance of long overlaps in thermophiles compared to non-thermophiles. On the other hand, short overlap (1-4 nucleotides) frequency (SOF) did not yield any direct correlation with OGT. However, the correlation of SOF with CAIavg (extent of variation of codon usage bias measured as the mean of codon adaptation index of all genes in a given genome) and IG% (proportion of intergenic regions) indicate that they might upregulate the aforementioned factors (CAIavg and IG%) which are already known to be vital forces for thermophilic adaptation. From these evidences, we propose that the OG content bears a strong link to thermophily. Long overlaps are important for their genome compaction and short overlaps are important to uphold high CAIavg. Our findings will surely help in better understanding of the significance of overlapping gene content in prokaryotic genomes.

PAP1 [poly(A) polymerase 1] homozygosity and hyperadenylation are major determinants of increased mRNA stability of CDR1 in azole-resistant clinical isolates of Candida albicans

Using genetically matched azole-susceptible (AS) and azole-resistant (AR) clinical isolates of Candida albicans, we recently demonstrated that CDR1 overexpression in AR isolates is due to its enhanced transcriptional activation and mRNA stability. This study examines the molecular mechanisms underlying enhanced CDR1 mRNA stability in AR isolates. Mapping of the 3' untranslated region (3' UTR) of CDR1 revealed that it was rich in adenylate/uridylate (AU) elements, possessed heterogeneous polyadenylation sites, and had putative consensus sequences for RNA-binding proteins. Swapping of heterologous and chimeric lacZ-CDR1 3' UTR transcriptional reporter fusion constructs did not alter the reporter activity in AS and AR isolates, indicating that cis-acting sequences within the CDR1 3' UTR itself are not sufficient to confer the observed differential mRNA decay. Interestingly, the poly(A) tail of the CDR1 mRNA of AR isolates was approximately 35-50 % hyperadenylated as compared with AS isolates. C. albicans poly(A) polymerase (PAP1), responsible for mRNA adenylation, resides on chromosome 5 in close proximity to the mating type-like (MTL) locus. Two different PAP1 alleles, PAP1-a/PAP1-alpha, were recovered from AS (MTL-a/MTL-alpha), while a single type of PAP1 allele (PAP1-alpha) was recovered from AR isolates (MTL-alpha/MTL-alpha). Among the heterozygous deletions of PAP1-a (Deltapap1-a/PAP1-alpha) and PAP1-alpha (PAP1-a/Deltapap1-alpha), only the former led to relatively enhanced drug resistance, to polyadenylation and to transcript stability of CDR1 in the AS isolate. This suggests a dominant negative role of PAP1-a in CDR1 transcript polyadenylation and stability. Taken together, our study provides the first evidence, to our knowledge, that loss of heterozygosity at the PAP1 locus is linked to hyperadenylation and subsequent increased stability of CDR1 transcripts, thus contributing to enhanced drug resistance.

MFS transportome of the human pathogenic yeast Candida albicans

BACKGROUND

The major facilitator superfamily (MFS) is one of the two largest superfamilies of membrane transporters present ubiquitously in bacteria, archaea, and eukarya and includes members that function as uniporters, symporters or antiporters. We report here the complete transportome of MFS proteins of a human pathogenic yeast Candida albicans.

RESULTS

Computational analysis of C. albicans genome enabled us to identify 95 potential MFS proteins which clustered into 17 families using Saier's Transport Commission (TC) system. Among these SP, DHA1, DHA2 and ACS represented major families consisting of 22, 22, 9 and 16 members, respectively. Family designations in C. albicans were validated by subjecting Saccharomyces cerevisiae genome to TC system. Based on the published available genomics/proteomics data, 87 of the putative MFS genes of C. albicans were found to express either at mRNA or protein levels. We checked the expression of the remaining 8 genes by using RT-PCR and observed that they are not expressed under basal growth conditions implying that either these 8 genes are expressed under specific growth conditions or they may be candidates for pseudogenes.

CONCLUSION

The in silico characterisation of MFS transporters in Candida albicans genome revealed a large complement of MFS transporters with most of them showing expression. Considering the clinical relevance of C. albicans and role of MFS members in antifungal resistance and nutrient transport, this analysis would pave way for identifying their physiological relevance.

mRNA processing in Antonospora locustae spores

Microsporidia are a group of intracellular parasites characterized by highly reduced and compact genomes. The presence of a high gene density had several consequences for microsporidian genomes, including a high frequency of overlap between transcripts of adjacent genes. This phenomenon is apparently widespread in microsporidia, and strongly correlated with gene density. However, all analyses to date have focused on one or a few transcripts from many loci, so it is unclear how diverse the pool of transcripts at a given locus may be. To address this question, we characterized initiation and termination points from 62 transcripts in gene-dense regions in Antonospora locustae spores using both conventional and fluorescence-based RACE-PCR procedures. In parallel, we investigated the abundance and nature of transcripts along a 6 kb region surrounding the actin locus of A. locustae using northern blotting, RACE-PCR and previously characterised EST sequences. Overall, we confirmed previous suggestions that most transcripts in A. locustae spores overlap with the downstream gene, but that at the 5' end untranslated regions are very short and overlap is rare. From fluorescence-based RACE-PCR we show that transcription of most genes (31 out of 34) initiates at a single position, whereas 35% of loci analyzed with 3' RACE polyadenylate mRNA at several sites. Finally, we identified the presence of previously unsuspected and very large transcripts in A. locustae spores. Those transcripts were found to overlap up to four open reading frames in different strands, adding a novel layer of complexity in the mRNA transcription of this microsporidian species.

Comparative profiling of overlapping transcription in the compacted genomes of microsporidia Antonospora locustae and Encephalitozoon cuniculi

Microsporidia are highly adapted parasites related to fungi with compact, gene-dense genomes. It has previously been shown in the microsporidian Antonospora locustae that transcripts from any given gene overlap with adjacent genes at a high frequency, perhaps due to the compact nature of its genome. However, it is still not known if this phenomenon is widespread among microsporidia or conserved between species, or even whether it is strictly correlated with compaction. To address these questions, we performed a comparison of transcription profiles in two microsporidian species, A. locustae and Encephalitozoon cuniculi. Transcription overlap was characterized at many A. locustae loci representing a range of gene densities, to determine if overlapping transcription correlates with the length of intergenic spacers. In parallel, we examined the first cases of transcription overlap in E. cuniculi. Using regions of the genome where the order of genes is conserved between A. locustae and E. cuniculi, we identified the transcriptional processing points in both species to determine how the process changes through evolutionary time. We show that there is little conservation of processing points between species and indeed that the process differs in important ways in the two genomes. Overall, A. locustae transcripts generally start just upstream of the start codon, but terminate well within or beyond downstream genes. In contrast, E. cuniculi transcripts often initiate within upstream genes, but more frequently terminate prior to the downstream gene. This process appears to have predictable characteristics within a given genome, but to be relatively flexible between species, presenting further challenges to the study of gene expression in these obligately intracellular parasites.

A high frequency of overlapping gene expression in compacted eukaryotic genomes

The gene density of eukaryotic nuclear genomes is generally low relative to prokaryotes, but several eukaryotic lineages (many parasites or endosymbionts) have independently evolved highly compacted, gene-dense genomes. The best studied of these are the microsporidia, highly adapted fungal parasites, and the nucleomorphs, relict nuclei of endosymbiotic algae found in cryptomonads and chlorarachniophytes. These systems are now models for the effects of compaction on the form and dynamics of the nuclear genome. Here we report a large-scale investigation of gene expression from compacted eukaryotic genomes. We have conducted EST surveys of the microsporidian Antonospora locustae and nucleomorphs of the cryptomonad Guillardia theta and the chlorarachniophyte Bigelowiella natans. In all three systems we find a high frequency of mRNA molecules that encode sequence from more than one gene. There is no bias for these genes to be on the same strand, so it is unlikely that these mRNAs represent operons. Instead, compaction appears to have reduced the intergenic regions to such an extent that control elements like promoters and terminators have been forced into or beyond adjacent genes, resulting in long untranslated regions that encode other genes. Normally, transcriptional overlap can interfere with expression of a gene, but these genomes cope with high frequencies of overlap and with termination signals within expressed genes. These findings also point to serious practical difficulties in studying expression in compacted genomes, because many techniques, such as arrays or serial analysis of gene expression will be misleading.

Reinvestigation of the Saccharomyces cerevisiae genome annotation by comparison to the genome of a related fungus: Ashbya gossypii

BACKGROUND

The recently sequenced genome of the filamentous fungus Ashbya gossypii revealed remarkable similarities to that of the budding yeast Saccharomyces cerevisiae both at the level of homology and synteny (conservation of gene order). Thus, it became possible to reinvestigate the S. cerevisiae genome in the syntenic regions leading to an improved annotation.

RESULTS

We have identified 23 novel S. cerevisiae open reading frames (ORFs) as syntenic homologs of A. gossypii genes; for all but one, homologs are present in other eukaryotes including humans. Other comparisons identified 13 overlooked introns and suggested 69 potential sequence corrections resulting in ORF extensions or ORF fusions with improved homology to the syntenic A. gossypii homologs. Of the proposed corrections, 25 were tested and confirmed by resequencing. In addition, homologs of nearly 1,000 S. cerevisiae ORFs, presently annotated as hypothetical, were found in A. gossypii at syntenic positions and can therefore be considered as authentic genes. Finally, we suggest that over 400 S. cerevisiae ORFs that overlap other ORFs in S. cerevisiae and for which no homolog can be detected in A. gossypii should be regarded as spurious.

CONCLUSIONS

Although, the S. cerevisiae genome is rightly considered as one of the most accurately sequenced and annotated eukaryotic genomes, we have shown that it still benefits substantially from comparison to the completed sequence and syntenic gene map of A. gossypii, an evolutionarily related fungus. This type of approach will strongly support the annotation of more complex genomes such as the human and murine genomes.

The POLARIS gene of Arabidopsis encodes a predicted peptide required for correct root growth and leaf vascular patterning

The POLARIS (PLS) gene of Arabidopsis was identified as a promoter trap transgenic line, showing beta-glucuronidase fusion gene expression predominantly in the embryonic and seedling root, with low expression in aerial parts. Cloning of the PLS locus revealed that the promoter trap T-DNA had inserted into a short open reading frame (ORF). Rapid amplification of cDNA ends PCR, RNA gel blot analysis, and RNase protection assays showed that the PLS ORF is located within a short ( approximately 500 nucleotides) auxin-inducible transcript and encodes a predicted polypeptide of 36 amino acid residues. pls mutants exhibit a short-root phenotype and reduced vascularization of leaves. pls roots are hyperresponsive to exogenous cytokinins and show increased expression of the cytokinin-inducible gene ARR5/IBC6 compared with the wild type. pls seedlings also are less responsive to the growth-inhibitory effects of exogenous auxin and show reduced expression of the auxin-inducible gene IAA1 compared with the wild type. The PLS peptide-encoding region of the cDNA partially complements the pls mutation and requires the PLS ORF ATG for activity, demonstrating the functionality of the peptide-encoding ORF. Ectopic expression of the PLS ORF reduces root growth inhibition by exogenous cytokinins and increases leaf vascularization. We propose that PLS is required for correct auxin-cytokinin homeostasis to modulate root growth and leaf vascular patterning.

Transcription initiation of genes associated with azole resistance in Candida albicans

Oral infections with the opportunistic pathogenic yeast Candida albicans are one of the earliest and most frequent infections in immunosuppressed individuals. In these populations, drug-resistant isolates have emerged with the widespread use of antifungal azole drugs. Many molecular mechanisms of resistance have been identified, including overexpression of two types of efflux pumps, the major facilitator MDR1 and the ABC-transporters (CDR1 and CDR2), and the overexpression or mutation of the target enzyme, ERG11. With overexpression of these four genes implicated in multidrug resistance, identification of regulatory regions of the promoters is important. 5' rapid amplification of cDNA ends (RACE) was used to identify transcription initiation sites for genes associated with multidrug resistance (CDR1, CDR2, MDR1 and ERG11). These results were confirmed by cloning and sequencing of 5' RACE products and by primer extension. This research will allow further analysis of the regulation of transcription for these genes.

Overlapping genes in parasitic protist Giardia lamblia

The parasitic protist Giardia lamblia lacks mitochondria and peroxisomes, as well as many typical membrane-bound organella characteristics of higher eukaryotic cells, together with extremely economized usage of DNA sequence, as demonstrated by the lack of introns. We describe here the presence of overlapping genes in G. lamblia, in which a part of the protein coding sequence of one mRNA exists in a region corresponding to the 3'-noncoding region of another mRNA transcribed from a gene on the opposite strand. Recently we isolated 13 kinesin-related cDNAs from G. lamblia. Nine of these cDNAs contain long 3'-noncoding sequences in which long open reading frames (ORFs) exist (in the remaining four cDNAs, the lengths of the 3'-noncoding sequences are very short). The predicted amino acid sequences of these ORFs were subjected to a search for homologies with sequences in databases. The amino acid sequences of the six ORFs exhibited significant sequence similarities with known sequences. These lines of evidence suggest the frequent occurrence of gene overlap in Giardial genome.

OTC and AUL1, two convergent and overlapping genes in the nuclear genome of Arabidopsis thaliana

In contrast to bacterial, fungal and vertebrate ornithine transcarbamylases (OTCs; EC 2.1.3.3), very little is known about the enzyme in plants. We report here the isolation of a T-DNA-tagged mutant displaying sensitivity to ornithine, whose characterization has allowed for the identification of several complementary and genomic DNA clones encoding the OTC and auxilin-like 1 (AUL1) proteins of the crucifer Arabidopsis thaliana. Transcript mapping revealed that at least 22 bp within the OTC-AUL1 intercoding region are transcribed from both strands, which makes this one of the rarely described cases of convergent and overlapping transcription units in the nuclear genome of a multicellular eukaryote. Transcription of the OTC gene was shown to be ubiquitous in aerial organs of adult plants, whereas that of AUL1 was obscured by the existence of a putative second copy of the gene. The OTC-AUL1 locus maps at the bottom of chromosome 1.