Non-coding RNA prediction and verification in Saccharomyces cerevisiae.

Non-coding RNA (ncRNA) play an important and varied role in cellular function. A significant amount of research has been devoted to computational prediction of these genes from genomic sequence, but the ability to do so has remained elusive due to a lack of apparent genomic features. In this work, thermodynamic stability of ncRNA structural elements, as summarized in a Z-score, is used to predict ncRNA in the yeast Saccharomyces cerevisiae. This analysis was coupled with comparative genomics to search for ncRNA genes on chromosome six of S. cerevisiae and S. bayanus. Sets of positive and negative control genes were evaluated to determine the efficacy of thermodynamic stability for discriminating ncRNA from background sequence. The effect of window sizes and step sizes on the sensitivity of ncRNA identification was also explored. Non-coding RNA gene candidates, common to both S. cerevisiae and S. bayanus, were verified using northern blot analysis, rapid amplification of cDNA ends (RACE), and publicly available cDNA library data. Four ncRNA transcripts are well supported by experimental data (RUF10, RUF11, RUF12, RUF13), while one additional putative ncRNA transcript is well supported but the data are not entirely conclusive. Six candidates appear to be structural elements in 5' or 3' untranslated regions of annotated protein-coding genes. This work shows that thermodynamic stability, coupled with comparative genomics, can be used to predict ncRNA with significant structural elements.

Recent evolution of the human pathogen Cryptococcus neoformans by intervarietal transfer of a 14-gene fragment.

The availability of the whole-genome sequence from the 2 known varieties of the human pathogenic fungus Cryptococcus neoformans provides an opportunity to study the relative contribution of divergence and introgression during the process of speciation in a genetically tractable organism. At the genomic level, these varieties are nearly completely syntenic, share approximately 85-90% nucleotide identity, and are believed to have diverged approximately 18 MYA. Via a comparative genomic approach, we identified a 14-gene region (approximately 40 kb) that is nearly identical between the 2 varieties that resulted from a nonreciprocal transfer event from var. grubii to var. neoformans approximately 2 MYA. The majority of clinical and environmental var. neoformans strains from around the world contain this sequence obtained from var. grubii. This introgression event likely occurred via an incomplete intervarietal sexual cycle, creating a hybrid intermediate where mobile elements common to both lineages mediated the exchange. The subsequent duplication in laboratory strains of a fragment of this same genomic region supports evolutionary theories that instabilities in subtelomeric regions promote adaptive evolution through gene amplification and subsequent adaptation. Along with a more ancient predicted transfer event in C. neoformans and a recently reported example from Saccharomyces cerevisiae, these data indicate that DNA exchange between closely related sympatric varieties or species may be a recurrent theme in the evolution of fungal species. It further suggests that although evolutionary divergence is the primary force driving speciation, rare introgression events also play a potentially important role.