A surprisingly large RNase P RNA in Candida glabrata

We have found an extremely large ribonuclease P (RNase P) RNA (RPR1) in the human pathogen Candida glabrata and verified that this molecule is expressed and present in the active enzyme complex of this hemiascomycete yeast. A structural alignment of the C. glabrata sequence with 36 other hemiascomycete RNase P RNAs (abbreviated as P RNAs) allows us to characterize the types of insertions. In addition, 15 P RNA sequences were newly characterized by searching in the recently sequenced genomes Candida albicans, C. glabrata, Debaryomyces hansenii, Eremothecium gossypii, Kluyveromyces lactis, Kluyveromyces waltii, Naumovia castellii, Saccharomyces kudriavzevii, Saccharomyces mikatae, and Yarrowia lipolytica; and by PCR amplification for other Candida species (Candida guilliermondii, Candida krusei, Candida parapsilosis, Candida stellatoidea, and Candida tropicalis). The phylogenetic comparative analysis identifies a hemiascomycete secondary structure consensus that presents a conserved core in all species with variable insertions or deletions. The most significant variability is found in C. glabrata P RNA in which three insertions exceeding in total 700 nt are present in the Specificity domain. This P RNA is more than twice the length of any other homologous P RNAs known in the three domains of life and is eight times the size of the smallest. RNase P RNA, therefore, represents one of the most diversified noncoding RNAs in terms of size variation and structural diversity.     

Preferential recombination between GC clusters in yeast mitochondrial DNA.

Yeast mitochondrial DNA molecules have long, AT-rich intergenic spacers punctuated by short GC clusters. GC-rich elements have previously been characterized by others as preferred sites for intramolecular recombination leading to the formation of subgenomic petite molecules. In the present study we show that GC clusters are favored sites for intermolecular recombination between a petite and the wild-type grande genome. The petite studied retains 6.5 kb of mitochondrial DNA reiterated tandemly to form molecules consisting of repeated units. Genetic selection for integration of tandem 6.5 kb repeats of the petite into the grande genome yielded a novel recombination event. One of two crossovers in a double exchange event occurred as expected in the 6.5 kb of matching sequence between the genomes, whereas the second exchange involved a 44 bp GC cluster in the petite and another 44 bp GC cluster in the grande genome 700 bp proximal to the region of homology. Creation of a mitochondrial DNA molecule with a repetitive region led to secondary recombination events that generated a family of molecules with zero to several petite units. The finding that 44 bp GC clusters are preferred as sites for intermolecular exchange adds to the data on petite excision implicating these elements as recombinational hotspots in the yeast mitochondrial genome.     

A bunch of fun-guys: the whole-genome view of yeast evolution

The resolution of the complete sequences of several hemiascomycete genomes provides new insights into the ways that yeast genomes change in size and in gene contents. These genomes provide evidence of whole-genome duplication occurring before the divergence of Saccharomyces cerevisiae and Candida glabrata, followed by massive gene loss that restored diploidy. The pattern of genome evolution in yeast differs from that in bacteria apparently as a result of stronger selective constraints on bacterial chromosomes.     

Why do some yeast species require niacin for growth? Different modes of NAD synthesis

NAD holds a key position in metabolism and cellular regulatory events as a major redox carrier and a signalling molecule. NAD biosynthesis pathways have been reconstructed and compared in seven yeast species with completely sequenced genomes, including Saccharomyces cerevisiae, Kluyveromyces lactis, Candida glabrata, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and Schizosaccharomyces pombe. Both amino acid and nucleotide sequence similarity analysis in silico indicated that de novo NAD biosynthesis might not exist in K. lactis, C. glabrata and Schiz. pombe, while other species have the kynurenine pathway. It also showed that the NAD salvage pathway via nicotinic acid and nicotinic acid mononucleotide is conserved in all of these yeasts. Deletion of KlNPT1 (the gene for nicotinate phosphoribosyl-transferase) is lethal, which demonstrates that this salvage pathway, utilizing exogenous nicotinic acid, is the unique route to synthesize NAD in K. lactis. The results suggested that the basis of the variation of niacin requirements in yeasts lies in their different combinations of NAD biosynthesis pathways. The de novo pathway is absent but the salvage pathway is conserved in niacin-negative yeasts, while both pathways coexist in niacin-positive yeasts.     

Comparative Analysis of the Mitochondrial Genomes of Callitettixini Spittlebugs (Hemiptera: Cercopidae) Confirms the Overall High Evolutionary Speed of the AT-Rich Region but Reveals the Presence of Short Conservative Elements at the Tribal Level

The present study compares the mitochondrial genomes of five species of the spittlebug tribe Callitettixini (Hemiptera: Cercopoidea: Cercopidae) from eastern Asia. All genomes of the five species sequenced are circular double-stranded DNA molecules and range from 15,222 to 15,637 bp in length. They contain 22 tRNA genes, 13 protein coding genes (PCGs) and 2 rRNA genes and share the putative ancestral gene arrangement of insects. The PCGs show an extreme bias of nucleotide and amino acid composition. Significant differences of the substitution rates among the different genes as well as the different codon position of each PCG are revealed by the comparative evolutionary analyses. The substitution speeds of the first and second codon position of different PCGs are negatively correlated with their GC content. Among the five species, the AT-rich region features great differences in length and pattern and generally shows a 2–5 times higher substitution rate than the fastest PCG in the mitochondrial genome, atp8. Despite the significant variability in length, short conservative segments were identified in the AT-rich region within Callitettixini, although absent from the other groups of the spittlebug superfamily Cercopoidea.     

双翅目昆虫线粒体基因组结构特点及其测序通用引物的设计和应用

研究双翅目昆虫线粒体基因组的结构特点, 并设计其测序的通用引物, 为今后双翅目昆虫线粒体基因组的研究提供参考和依据。利用比较基因组学和生物信息学方法, 分析了已经完全测序的26个双翅目昆虫线粒体基因组的结构特点、 碱基组成和保守区, 并据此设计了双翅目昆虫基因组测序的通用引物。结果表明： 双翅目昆虫线粒体基因组长14 503~19 517 bp, 其结构保守, 含有37个编码基因, 包括13个蛋白质编码基因, 22个tRNA编码基因和2个rRNA编码基因, 此外还包含一段长度差异很大的非编码区(AT富含区)。基因组内基因排列次序稳定, 除个别基因外, 其余都与黑腹果蝇Drosophila melanogaster基因排列次序一致。基因组的碱基组成不均衡, AT含量在72.59%~85.15%之间, 碱基使用存在偏向性, 偏好使用AC碱基。全基因组的核苷酸和氨基酸序列保守, 共鉴定了11个保守区。在保守区内共设计了26对双翅目线粒体基因组测序通用引物, 扩增的目标片段都在1 200 bp以内。将该套通用引物用于葱蝇Delia antiqua线粒体全基因组测序, 结果证明其高效、 合用。     

Biased distributions and decay of long interspersed nuclear elements in the chicken genome

The genomes of birds are much smaller than mammalian genomes, and transposable elements (TEs) make up only 10% of the chicken genome, compared with the 45% of the human genome. To study the mechanisms that constrain the copy numbers of TEs, and as a consequence the genome size of birds, we analyzed the distributions of LINEs (CR1's) and SINEs (MIRs) on the chicken autosomes and Z chromosome. We show that (1) CR1 repeats are longest on the Z chromosome and their length is negatively correlated with the local GC content; (2) the decay of CR1 elements is highly biased, and the 5'-ends of the insertions are lost much faster than their 3'-ends; (3) the GC distribution of CR1 repeats shows a bimodal pattern with repeats enriched in both AT-rich and GC-rich regions of the genome, but the CR1 families show large differences in their GC distribution; and (4) the few MIRs in the chicken are most abundant in regions with intermediate GC content. Our results indicate that the primary mechanism that removes repeats from the chicken genome is ectopic exchange and that the low abundance of repeats in avian genomes is likely to be the consequence of their high recombination rates.     

Exon-intron structure of genes in complete fungal genomes

Computer analysis of genes was performed for lower fungi Aspergillus fumigatus, Candida glabrata, Cryptococcus neoformans, Debaryomyces hansenii, Encephalitozoon cuniculi, Eremothecium gossypii, Kluyveromyces lactis, Magnaporthe grisea, Neurospora crassa, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Ustilago maydis, and Yarrowia lipolytica. The content of genes with an exon-intron structure in their genomes varied from 0.7 to 97.0%. The exon-intron structure substantially changes with an increasing portion of intron-containing genes. Gene size and total exon length proved to linearly depend on the intron number in the A. fumigatus, C. neoformans, M. grisea, N. crassa, S. pombe, and U. maydis genomes.     

Evolutionary relationships among pathogenic Candida species and relatives.

Small subunit rRNA sequences have been determined for 10 of the most clinically important pathogenic species of the yeast genus Candida (including Torulopsis [Candida] glabrata and Yarrowia [Candida] lipolytica) and for Hansenula polymorpha. Phylogenetic analyses of these sequences and those of Saccharomyces cerevisiae, Kluyveromyces marxianus var. lactis, and Aspergillus fumigatus indicate that Candida albicans, C. tropicalis, C. parapsilosis, and C. viswanathii form a subgroup within the genus. The remaining significant pathogen, T. glabrata, falls into a second, distinct subgroup and is specifically related to S. cerevisiae and more distantly related to C. kefyr (psuedotropicalis) and K. marxianus var. lactis. The 18S rRNA sequence of Y. lipolytica has evolved rapidly in relation to the other Candida sequences examined and appears to be only distantly related to them. As anticipated, species of several other genera appear to bear specific relationships to members of the genus Candida.     

A mitochondrial molecular marker, ori-rep-tra, for differentiation of yeast species.

Yeasts exhibit various mechanisms for the inheritance of their mitochondrial genomes. Differences among these mechanisms are based on variations within nuclear as well as mitochondrial genetic elements. Here we report diagnostic differences in the presence of biologically active mitochondrial intergenic sequences, ori-reptra, among related yeasts in the genera Saccharomyces, Arxiozyma, Debaryomyces, Kluyveromyces, Pachytichospora, Torulaspora, and Zygosaccharomyces. A molecular probe containing ori-rep-tra can be employed specifically for the differentiation and identification of isolates belonging to the species complex Saccharomyces sensu stricto.     

Combined phylogenetic and neighbourhood analysis of the hexose transporters and glucose sensors in yeasts

The sugar porter family in yeasts encompasses a wide variety of transporters including the hexose transporters and glucose sensors. We analysed a total of 75 members from both groups in nine hemiascomycetous species, with complete and well-annotated genomes: Saccharomyces cerevisiae, Candida glabrata, Zygosaccharomyces rouxii, Kluyveromyces thermotolerans, Saccharomyces kluyverii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii and Yarrowia lipolytica . We present a model for the evolution of the hexose transporters and glucose sensors, supported by two types of complementary evidences: phylogeny and neighbourhood analysis. Five lineages of evolution were identified and discussed according to different mechanisms of gene evolution: lineage A for HXT1, HXT3 , HXT4, HXT5 , HXT6 and HXT7 ; lineage B for HXT2 and HXT10 ; lineage C for HXT8 ; lineage D for HXT14 ; and lineage E for SNF3 and RGT2 .     

Megasatellites: a peculiar class of giant minisatellites in genes involved in cell adhesion and pathogenicity in Candida glabrata

Minisatellites are DNA tandem repeats that are found in all sequenced genomes. In the yeast Saccharomyces cerevisiae, they are frequently encountered in genes encoding cell wall proteins. Minisatellites present in the completely sequenced genome of the pathogenic yeast Candida glabrata were similarly analyzed, and two new types of minisatellites were discovered: minisatellites that are composed of two different intermingled repeats (called compound minisatellites), and minisatellites containing unusually long repeated motifs (126-429 bp). These long repeat minisatellites may reach unusual length for such elements (up to 10 kb). Due to these peculiar properties, they have been named 'megasatellites'. They are found essentially in genes involved in cell-cell adhesion, and could therefore be involved in the ability of this opportunistic pathogen to colonize the human host. In addition to megasatellites, found in large paralogous gene families, there are 93 minisatellites with simple shorter motifs, comparable to those found in S. cerevisiae. Most of the time, these minisatellites are not conserved between C. glabrata and S. cerevisiae, although their host genes are well conserved, raising the question of an active mechanism creating minisatellites de novo in hemiascomycetes.     

Detection of nucleolar organizer and mitochondrial DNA insertion regions based on the isochore map of Arabidopsis thaliana

Eukaryotic genomes are composed of isochores, i.e. long sequences relatively homogeneous in GC content. In this paper, the isochore structure of Arabidopsis thaliana genome has been studied using a windowless technique based on the Z curve method and intuitive curves are drawn for all the five chromosomes. Using these curves, we can calculate the GC content at any resolution, even at the base level. It is observed that all the five chromosomes are composed of several GC-rich and AT-rich regions alternatively. Usually, these regions, named 'isochore-like regions', have large fluctuations in the GC content. Five isochores with little fluctuations are also observed. Detailed analyses have been performed for these isochores. A GC-rich 'isochore-like region' and a GC-isochore in chromosome II and IV, respectively, are the nucleolar organizer regions (NORs), and genes located in the two regions prefer to use GC-ending codons. Another GC-isochore located in chromosome II is a mitochondrial DNA insertion region, the position and size of this region is precisely predicted by the current method. The amino acid usage and codon preference of genes in this organellar-to-nuclear transfer region show significant difference from other regions. Moreover, the centromeres are located in GC-rich 'isochore-like regions' in all the five chromosomes. The current method can provide a useful tool for analyzing whole genomic sequences of eukaryotes.     

Complete sequences of mitochondria genomes of Aedes aegypti and Culex quinquefasciatus and comparative analysis of mitochondrial DNA fragments inserted in the nuclear genomes

We present complete sequences of the mitochondrial genomes for two important mosquitoes, Aedes aegypti and Culex quinquefasciatus, that are major vectors of dengue virus and lymphatic filariasis, respectively. The A. aegypti mitochondrial genome is 16,655 bp in length and that of C. quinquefasciatus is 15,587 bp, yet both contain 13 protein coding genes, 22 transfer RNA (tRNA) genes, one 12S ribosomal RNA (rRNA) gene, one 16S rRNA gene and a control region (CR) in the same order. The difference in the genome size is due to the difference in the length of the control region. We also analyzed insertions of nuclear copies of mtDNA-like sequences (NUMTs) in a comparative manner between the two mosquitoes. The NUMT sequences occupy ~0.008% of the A. aegypti genome and ~0.001% of the C. quinquefasciatus genome. Several NUMTs were found localized in the introns of predicted protein coding genes in both genomes (32 genes in A. aegypti but only four in C. quinquefasciatus). None of these NUMT-containing genes had an ortholog between the two species or had paralogous copies within a genome that was also NUMT-containing. It was further observed that the NUMT-containing genes were relatively longer but had lower GC content compared to the NUMT-less paralogous copies. Moreover, stretches of homologies are present among the genic and non-genic NUMTs that may play important roles in genomic rearrangement of NUMTs in these genomes. Our study provides new insights on understanding the roles of nuclear mtDNA sequences in genome complexities of these mosquitoes.     

var1 Gene on the mitochondrial genome of Torulopsis glabrata

We have cloned and sequenced a region of the Torulopsis glabrata mitochondrial genome homologous to the Saccharomyces cerevisiae var1 gene (var1Sc). An open reading frame that could encode a protein of 339 amino acids was found with 72.7% amino acid and 85.3% nucleotide sequence homology to the S. cerevisiae var1 gene. The T. glabrata gene (var1Tg) is transcribed yielding two stable RNAs, a more abundant 13.5 S RNA and a less abundant 18 S species. We have also identified a candidate for a T. glabrata var1 protein among mitochondrial translation products labeled in isolated mitochondria. The var1Tg gene is even more A + T-rich (93%) than var1Sc (89.6%) and has conserved the strong codon bias of var1Sc. Major differences between the two sequences were found. Significant among these are that no GC clusters are found in var1Tg and the sequences surrounding each of the sites where known polymorphisms exist in var1Sc have deletions at the corresponding sites in var1Tg. These data are discussed with respect to possible origins of these var1 genes and translocation of GC clusters in S. cerevisiae mitochondrial DNA.     

Promiscuous DNA in the nuclear genomes of hemiascomycetous yeasts

Transfer of fragments of mtDNA to the nuclear genome is a general phenomenon that gives rise to NUMTs (NUclear sequences of MiTochondrial origin). We present here the first comparative analysis of the NUMT content of entirely sequenced species belonging to a monophyletic group, the hemiascomycetous yeasts ( Candida glabrata, Kluyveromyces lactis, Kluyveromyces thermotolerans, Debaryomyces hansenii and Yarrowia lipolytica , along with the updated NUMT content of Saccharomyces cerevisiae ). This study revealed a huge diversity in NUMT number and organization across the six species. Debaryomyces hansenii harbors the highest number of NUMTs (145), half of which are distributed in numerous large mosaics of up to eight NUMTs arising from multiple noncontiguous mtDNA fragments inserted at the same chromosomal locus. Most NUMTs, in all species, are found within intergenic regions including seven NUMTs in pseudogenes. However, five NUMTs overlap a gene, suggesting a positive impact of NUMTs on protein evolution. Contrary to the other species, K. lactis and K. thermotolerans harbor only a few diverged NUMTs, suggesting that mitochondrial transfer to the nuclear genome has decreased or ceased in these phylogenetic branches. The dynamics of NUMT acquisition and loss are illustrated here by their species-specific distribution.     

Asymmetry Indices for Analysis and Prediction of Replication Origins in Eukaryotic Genomes

DNA replication was recently shown to induce the formation of compositional skews in the genomes of the yeasts Saccharomyces cerevisiae and Kluyveromyces lactis. In this work, I have characterized further GC and TA skew variations in the vicinity of S. cerevisiae replication origins and termination sites, and defined asymmetry indices for origin analysis and prediction. The presence of skew jumps at some termination sites in the S. cerevisiae genome was established. The majority of S. cerevisiae replication origins are marked by an oriented consensus sequence called ACS, but no evidence could be found for asymmetric origin firing that would be linked to ACS orientation. Asymmetry indices related to GC and TA skews were defined, and a global asymmetry index I_GC,TA was described. I_GC,TA was found to strongly correlate with origin efficiency in S. cerevisiae and to allow the determination of sets of intergenes significantly enriched in origin loci. The generalized use of asymmetry indices for origin prediction in naive genomes implies the determination of the direction of the skews, i.e. the identification of which strand, leading or lagging, is enriched in G and which one is enriched in T. Recent work indicates that in Candida albicans and in several related species, centromeres contain early and efficient replication origins. It has been proposed that the skew jumps observed at these positions would reflect the activity of these origins, thus allowing to determine the direction of the skews in these genomes. However, I show here that the skew jumps at C. albicans centromeres are not related to replication and that replication-associated GC and TA skews in C. albicans have in fact the opposite directions of what was proposed.