The diversity of retrotransposons in the yeast Cryptococcus neoformans

YaliCMulti and YaliHMulti: Stable, efficient multi-copy integration tools for engineering Yarrowia lipolytica

Yarrowia lipolytica is widely used in biotechnology to produce recombinant proteins, food ingredients and diverse natural products. However, unstable expression of plasmids, difficult and time-consuming integration of single and low-copy-number plasmids hampers the construction of efficient production pathways and application to industrial production. Here, by exploiting sequence diversity in the long terminal repeats (LTRs) of retrotransposons and ribosomal DNA (rDNA) sequences, a set of vectors and methods that can recycle multiple and high-copy-number plasmids was developed that can achieve stable integration of long-pathway genes in Y. lipolytica. By combining these sequences, amino acids and antibiotic tags with the Cre-LoxP system, a series of multi-copy site integration recyclable vectors were constructed and assessed using the green fluorescent protein (HrGFP) reporter system. Furthermore, by combining the consensus sequence with the vector backbone of a rapidly degrading selective marker and a weak promoter, multiple integrated high-copy-number vectors were obtained and high levels of stable HrGFP expression were achieved. To validate the universality of the tools, simple integration of essential biosynthesis modules was explored, and 7.3 g/L of L-ergothioneine and 8.3 g/L of (2S)-naringenin were achieved in a 5 L fermenter, the highest titres reported to date for Y. lipolytica. These novel multi-copy genome integration strategies provide convenient and effective tools for further metabolic engineering of Y. lipolytica.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Uncontrolled transposition following RNAi loss causes hypermutation and antifungal drug resistance in clinical isolates of Cryptococcus neoformans

Cryptococcus neoformans infections cause approximately 15% of AIDS-related deaths owing to a combination of limited antifungal therapies and drug resistance. A collection of clinical and environmental C. neoformans isolates were assayed for increased mutation rates via fluctuation analysis, and we identified two hypermutator C. neoformans clinical isolates with increased mutation rates when exposed to the combination of rapamycin and FK506. Sequencing of drug target genes found that Cnl1 transposon insertions conferred the majority of resistance to rapamycin and FK506 and could also independently cause resistance to 5-fluoroorotic acid and the clinically relevant antifungal 5-flucytosine. Whole-genome sequencing revealed both hypermutator genomes harbour a nonsense mutation in the RNA-interference component ZNF3 and hundreds of Cnl1 elements organized into massive subtelomeric arrays on each of the fourteen chromosomes. Quantitative trait locus mapping in 28 progeny derived from a cross between a hypermutator and wild-type identified a locus associated with hypermutation that included znf3. CRISPR editing of the znf3 nonsense mutation abolished hypermutation and restored small-interfering-RNA production. We conclude that hypermutation and drug resistance in these clinical isolates result from RNA-interference loss and accumulation of Cnl1 elements.

Obligate sexual reproduction of a homothallic fungus closely related to the Cryptococcus pathogenic species complex

eLife digest

.

Fungi are enigmatic organisms that flourish in soil, on decaying plants, or during infection of animals or plants. Growing in myriad forms, from single-celled yeast to multicellular molds and mushrooms, fungi have also evolved a variety of strategies to reproduce. Normally, fungi reproduce in one of two ways: either they reproduce asexually, with one individual producing a new individual identical to itself, or they reproduce sexually, with two individuals of different 'mating types' contributing to produce a new individual. However, individuals of some species exhibit 'homothallism' or self-fertility: these individuals can produce reproductive cells that are universally compatible, and therefore can reproduce sexually with themselves or with any other cell in the population.

Homothallism has evolved multiple times throughout the fungal kingdom, suggesting it confers advantage when population numbers are low or mates are hard to find. Yet some homothallic fungi been overlooked compared to heterothallic species, whose mating types have been well characterised. Understanding the genetic basis of homothallism and how it evolved in different species can provide insights into pathogenic species that cause fungal disease.

With that in mind, Passer, Clancey et al. explored the genetic basis of homothallism in Cryptococcus depauperatus, a close relative of C. neoformans, a species that causes fungal infections in humans. A combination of genetic sequencing techniques and experiments were applied to analyse, compare, and manipulate C. depauperatus' genome to see how this species evolved self-fertility.

Passer, Clancey et al. showed that C. depauperatus evolved the ability to reproduce sexually by itself via a unique evolutionary pathway. The result is a form of homothallism never reported in fungi before. C. depauperatus lost some of the genes that control mating in other species of fungi, and acquired genes from the opposing mating types of a heterothallic ancestor to become self-fertile.

Passer, Clancey et al. also found that, unlike other Cryptococcus species that switch between asexual and sexual reproduction, C. depauperatus grows only as long, branching filaments called hyphae, a sexual form. The species reproduces sexually with itself throughout its life cycle and is unable to produce a yeast (asexual) form, in contrast to other closely related species.

This work offers new insights into how different modes of sexual reproduction have evolved in fungi. It also provides another interesting case of how genome plasticity and evolutionary pressures can produce similar outcomes, homothallism, via different evolutionary paths. Lastly, assembling the complete genome of C. depauperatus will foster comparative studies between pathogenic and non-pathogenic Cryptococcus species.

A comprehensive annotation dataset of intact LTR retrotransposons of 300 plant genomes

LTR retrotransposons (LTR-RTs) are ubiquitous and represent the dominant repeat element in plant genomes, playing important roles in functional variation, genome plasticity and evolution. With the advent of new sequencing technologies, a growing number of whole-genome sequences have been made publicly available, making it possible to carry out systematic analyses of LTR-RTs. However, a comprehensive and unified annotation of LTR-RTs in plant groups is still lacking. Here, we constructed a plant intact LTR-RTs dataset, which is designed to classify and annotate intact LTR-RTs with a standardized procedure. The dataset currently comprises a total of 2,593,685 intact LTR-RTs from genomes of 300 plant species representing 93 families of 46 orders. The dataset is accompanied by sequence, diverse structural and functional annotation, age determination and classification information associated with the LTR-RTs. This dataset will contribute valuable resources for investigating the evolutionary dynamics and functional implications of LTR-RTs in plant genomes.

Factors enforcing the species boundary between the human pathogens Cryptococcus neoformans and Cryptococcus deneoformans

Hybridization has resulted in the origin and variation in extant species, and hybrids continue to arise despite pre- and post-zygotic barriers that limit their formation and evolutionary success. One important system that maintains species boundaries in prokaryotes and eukaryotes is the mismatch repair pathway, which blocks recombination between divergent DNA sequences. Previous studies illuminated the role of the mismatch repair component Msh2 in blocking genetic recombination between divergent DNA during meiosis. Loss of Msh2 results in increased interspecific genetic recombination in bacterial and yeast models, and increased viability of progeny derived from yeast hybrid crosses. Hybrid isolates of two pathogenic fungal Cryptococcus species, Cryptococcus neoformans and Cryptococcus deneoformans, are isolated regularly from both clinical and environmental sources. In the present study, we sought to determine if loss of Msh2 would relax the species boundary between C. neoformans and C. deneoformans. We found that crosses between these two species in which both parents lack Msh2 produced hybrid progeny with increased viability and high levels of aneuploidy. Whole-genome sequencing revealed few instances of recombination among hybrid progeny and did not identify increased levels of recombination in progeny derived from parents lacking Msh2. Several hybrid progeny produced structures associated with sexual reproduction when incubated alone on nutrient-rich medium in light, a novel phenotype in Cryptococcus. These findings represent a unique, unexpected case where rendering the mismatch repair system defective did not result in increased meiotic recombination across a species boundary. This suggests that alternative pathways or other mismatch repair components limit meiotic recombination between homeologous DNA and enforce species boundaries in the basidiomycete Cryptococcus species.

Several mechanisms enforce species boundaries by either preventing the formation of hybrid zygotes, known as pre-zygotic barriers, or preventing the viability and fecundity of hybrids, known as post-zygotic barriers. Despite these barriers, interspecific hybrids form at an appreciable frequency, such as hybrid isolates of the human fungal pathogenic species, Cryptococcus neoformans and Cryptococcus deneoformans, which are regularly isolated from both clinical and environmental sources. C. neoformans x C. deneoformans hybrids are typically highly aneuploid, sterile, and display phenotypes intermediate to those of either parent, although self-fertile isolates and transgressive phenotypes have been observed. One important mechanism known to enforce species boundaries or lead to incipient speciation is the DNA mismatch repair system, which blocks recombination between divergent DNA sequences during meiosis. The aim of this study was to determine if genetically deleting the DNA mismatch repair component Msh2 would relax the species boundary between C. neoformans and C. deneoformans. Progeny derived from C. neoformans x C. deneoformans crosses in which both parental strains lacked Msh2 had higher viability, and unlike previous studies in Saccharomyces, these Cryptococcus hybrid progeny had higher levels of aneuploidy and no observable increase in meiotic recombination at the whole-genome level.

Diverse transposable element landscapes in pathogenic and nonpathogenic yeast models: the value of a comparative perspective

Genomics and other large-scale analyses have drawn increasing attention to the potential impacts of transposable elements (TEs) on their host genomes. However, it remains challenging to transition from identifying potential roles to clearly demonstrating the level of impact TEs have on genome evolution and possible functions that they contribute to their host organisms. I summarize TE content and distribution in four well-characterized yeast model systems in this review: the pathogens Candida albicans and Cryptococcus neoformans, and the nonpathogenic species Saccharomyces cerevisiae and Schizosaccharomyces pombe. I compare and contrast their TE landscapes to their lifecycles, genomic features, as well as the presence and nature of RNA interference pathways in each species to highlight the valuable diversity represented by these models for functional studies of TEs. I then review the regulation and impacts of the Ty1 and Ty3 retrotransposons from Saccharomyces cerevisiae and Tf1 and Tf2 retrotransposons from Schizosaccharomyces pombe to emphasize parallels and distinctions between these well-studied elements. I propose that further characterization of TEs in the pathogenic yeasts would enable this set of four yeast species to become an excellent set of models for comparative functional studies to address outstanding questions about TE-host relationships.

Genetic and Genomic Analyses Reveal Boundaries between Species Closely Related to Cryptococcus Pathogens

Speciation is a central mechanism of biological diversification. While speciation is well studied in plants and animals, in comparison, relatively little is known about speciation in fungi. One fungal model is the Cryptococcus genus, which is best known for the pathogenic Cryptococcus neoformans/Cryptococcus gattii species complex that causes >200,000 new human infections annually. Elucidation of how these species evolved into important human-pathogenic species remains challenging and can be advanced by studying the most closely related nonpathogenic species, Cryptococcus amylolentus and Tsuchiyaea wingfieldii However, these species have only four known isolates, and available data were insufficient to determine species boundaries within this group. By analyzing full-length chromosome assemblies, we reappraised the phylogenetic relationships of the four available strains, confirmed the genetic separation of C. amylolentus and T. wingfieldii (now Cryptococcus wingfieldii), and revealed an additional cryptic species, for which the name Cryptococcus floricola is proposed. The genomes of the three species are ∼6% divergent and exhibit significant chromosomal rearrangements, including inversions and a reciprocal translocation that involved intercentromeric ectopic recombination, which together likely impose significant barriers to genetic exchange. Using genetic crosses, we show that while C. wingfieldii cannot interbreed with any of the other strains, C. floricola can still undergo sexual reproduction with C. amylolentus However, most of the resulting spores were inviable or sterile or showed reduced recombination during meiosis, indicating that intrinsic postzygotic barriers had been established. Our study and genomic data will foster additional studies addressing fungal speciation and transitions between nonpathogenic and pathogenic Cryptococcus lineages.IMPORTANCE The evolutionary drivers of speciation are critical to our understanding of how new pathogens arise from nonpathogenic lineages and adapt to new environments. Here we focus on the Cryptococcus amylolentus species complex, a nonpathogenic fungal lineage closely related to the human-pathogenic Cryptococcus neoformans/Cryptococcus gattii complex. Using genetic and genomic analyses, we reexamined the species boundaries of four available isolates within the C. amylolentus complex and revealed three genetically isolated species. Their genomes are ∼6% divergent and exhibit chromosome rearrangements, including translocations and small-scale inversions. Although two of the species (C. amylolentus and newly described C. floricola) were still able to interbreed, the resulting hybrid progeny were usually inviable or sterile, indicating that barriers to reproduction had already been established. These results advance our understanding of speciation in fungi and highlight the power of genomics in assisting our ability to correctly identify and discriminate fungal species.

Systematic survey of plant LTR-retrotransposons elucidates phylogenetic relationships of their polyprotein domains and provides a reference for element classification

BACKGROUND

Plant LTR-retrotransposons are classified into two superfamilies, Ty1/copia and Ty3/gypsy. They are further divided into an enormous number of families which are, due to the high diversity of their nucleotide sequences, usually specific to a single or a group of closely related species. Previous attempts to group these families into broader categories reflecting their phylogenetic relationships were limited either to analyzing a narrow range of plant species or to analyzing a small numbers of elements. Furthermore, there is no reference database that allows for similarity based classification of LTR-retrotransposons.

RESULTS

We have assembled a database of retrotransposon encoded polyprotein domains sequences extracted from 5410 Ty1/copia elements and 8453 Ty3/gypsy elements sampled from 80 species representing major groups of green plants (Viridiplantae). Phylogenetic analysis of the three most conserved polyprotein domains (RT, RH and INT) led to dividing Ty1/copia and Ty3/gypsy retrotransposons into 16 and 14 lineages respectively. We also characterized various features of LTR-retrotransposon sequences including additional polyprotein domains, extra open reading frames and primer binding sites, and found that the occurrence and/or type of these features correlates with phylogenies inferred from the three protein domains.

CONCLUSIONS

We have established an improved classification system applicable to LTR-retrotransposons from a wide range of plant species. This system reflects phylogenetic relationships as well as distinct sequence and structural features of the elements. A comprehensive database of retrotransposon protein domains (REXdb) that reflects this classification provides a reference for efficient and unified annotation of LTR-retrotransposons in plant genomes. Access to REXdb related tools is implemented in the RepeatExplorer web server (https://repeatexplorer-elixir.cerit-sc.cz/) or using a standalone version of REXdb that can be downloaded seaparately from RepeatExplorer web page (http://repeatexplorer.org/).

RNAi is a critical determinant of centromere evolution in closely related fungi

The “centromere paradox” refers to rapidly evolving and highly diverse centromere DNA sequences even in closely related eukaryotes. However, factors contributing to this rapid divergence are largely unknown. Here, we identified large regional, LTR retrotransposon-rich centromeres in a group of human fungal pathogens belonging to the Cryptococcus species complex. We provide evidence that loss-of-functional RNAi machinery and possibly cytosine DNA methylation trigger instability of the genome by activation of centromeric retrotransposons presumably suppressed by RNAi. We propose that RNAi, together with cytosine DNA methylation, serves as a critical determinant that maintains repetitive transposon-rich centromere structures. This study explores the direct link between RNAi and centromere structure evolution.

The centromere DNA locus on a eukaryotic chromosome facilitates faithful chromosome segregation. Despite performing such a conserved function, centromere DNA sequence as well as the organization of sequence elements is rapidly evolving in all forms of eukaryotes. The driving force that facilitates centromere evolution remains an enigma. Here, we studied the evolution of centromeres in closely related species in the fungal phylum of Basidiomycota. Using ChIP-seq analysis of conserved inner kinetochore proteins, we identified centromeres in three closely related Cryptococcus species: two of which are RNAi-proficient, while the other lost functional RNAi. We find that the centromeres in the RNAi-deficient species are significantly shorter than those of the two RNAi-proficient species. While centromeres are LTR retrotransposon-rich in all cases, the RNAi-deficient species lost all full-length retroelements from its centromeres. In addition, centromeres in RNAi-proficient species are associated with a significantly higher level of cytosine DNA modifications compared with those of RNAi-deficient species. Furthermore, when an RNAi-proficient Cryptococcus species and its RNAi-deficient mutants were passaged under similar conditions, the centromere length was found to be occasionally shortened in RNAi mutants. In silico analysis of predicted centromeres in a group of closely related Ustilago species, also belonging to the Basidiomycota, were found to have undergone a similar transition in the centromere length in an RNAi-dependent fashion. Based on the correlation found in two independent basidiomycetous species complexes, we present evidence suggesting that the loss of RNAi and cytosine DNA methylation triggered transposon attrition, which resulted in shortening of centromere length during evolution.

Comparative genomics of Cryptococcus neoformans var. grubii associated with meningitis in HIV infected and uninfected patients in Vietnam

The vast burden of cryptococcal meningitis occurs in immunosuppressed patients, driven by HIV, and is caused by Cryptococcus neoformansvar. grubii. We previously reported cryptococcal meningitis in Vietnam arising atypically in HIV uninfected, apparently immunocompetent patients, caused by a single amplified fragment length polymorphism (AFLP) cluster of C. neoformansvar. grubii (VNIγ). This variant was less common in HIV infected individuals; it remains unclear why this lineage is associated with apparently immunocompetent patients. To study this host tropism we aimed to further our understanding of clinical phenotype and genomic variation within Vietnamese C. neoformansvar. grubii. After performing MLST on C. neoformans clinical isolates we identified 14 sequence types (STs); ST5 correlated with the VNIγ cluster. We next compared clinical phenotype by lineage and found HIV infected patients with cryptococcal meningitis caused by ST5 organisms were significantly more likely to have lymphadenopathy (11% vs. 4%, p = 0.05 Fisher’s exact test) and higher blood lymphocyte count (median 0.76 versus 0.55 X10⁹ cells/L, p = 0.001, Kruskal-Wallis test). Furthermore, survivors of ST5 infections had evidence of worse disability outcomes at 70 days (72.7% (40/55) in ST5 infections versus 57.1% (52/91) non-ST5 infections (OR 2.11, 95%CI 1.01 to 4.41), p = 0.046). To further investigate the relationship between strain and disease phenotype we performed genome sequencing on eight Vietnamese C. neoformansvar. grubii. Eight genome assemblies exhibited >99% nucleotide sequence identity and we identified 165 kbp of lineage specific to Vietnamese isolates. ST5 genomes harbored several strain specific regions, incorporating 19 annotated coding sequences and eight hypothetical proteins. These regions included a phenolic acid decarboxylase, a DEAD-box ATP-dependent RNA helicase 26, oxoprolinases, a taurine catabolism dioxygenase, a zinc finger protein, membrane transport proteins and various drug transporters. Our work outlines the complexity of genomic pathogenicity in cryptococcal infections and identifies a number of gene candidates that may aid the disaggregation of the pathways associated with the pathogenesis of Cryptococcus neoformansvar. grubii.

Cryptococcal meningitis is a brain infection caused by a yeast, Cryptococcus neoformans, and results in an estimated 600 000 deaths each year. Disease usually only occurs in patients who have some problem with their immune systems—most commonly Human Immunodeficiency Virus (HIV) infection. However, it is increasingly recognized that disease can occur, particularly in southeast and east Asia, in patients with apparently normal immune systems (‘immunocompetent’). We previously showed that almost all infections in immunocompetent patients in Vietnam are due to just one small ‘family’ (or lineage) of Cryptococcus neoformans var. grubii, which we called VNIγ. This is in contrast to disease in HIV infected patients, which can be caused by a number of different families. This suggests that VNIγ strains have an increased ability to cause disease. Here, we define the pattern of disease caused by VNIγ infections compared with other strains in HIV infected patients, and use whole genome sequencing—comparing the entire genetic codes from different strains—to try and understand which genes give the VNIγ family this special ability to cause disease in immunocompetent patients.

sRNAs as possible regulators of retrotransposon activity in Cryptococcus gattii VGII

Background

The absence of Argonaute genes in the fungal pathogen Cryptococcus gattii R265 and other VGII strains indicates that yeasts of this genotype cannot have a functional RNAi pathway, an evolutionarily conserved gene silencing mechanism performed by small RNAs. The success of the R265 strain as a pathogen that caused the Pacific Northwest and Vancouver Island outbreaks may imply that RNAi machinery loss could be beneficial under certain circumstances during evolution. As a result, a hypermutant phenotype would be created with high rates of genome retrotransposition, for instance. This study therefore aimed to evaluate in silicio the effect of retrotransposons and their control mechanisms by small RNAs on genomic stability and synteny loss of C. gattii R265 through retrotransposons sequence comparison and orthology analysis with other 16 C. gattii genomic sequences available.

Results

Retrotransposon mining identified a higher sequence count to VGI genotype compared to VGII, VGIII, and VGIV. However, despite the lower retrotransposon number, VGII exhibited increased synteny loss and genome rearrangement events. RNA-Seq analysis indicated highly expressed retrotransposons as well as sRNA production.

Conclusions

Genome rearrangement and synteny loss may suggest a greater retrotransposon mobilization caused by RNAi pathway absence, but the effective presence of sRNAs that matches retrotransposon sequences means that an alternative retrotransposon silencing mechanism could be active in genomic integrity maintenance of C. gattii VGII strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3688-4) contains supplementary material, which is available to authorized users.

Trichoderma reesei complete genome sequence, repeat-induced point mutation, and partitioning of CAZyme gene clusters

BACKGROUND

Trichoderma reesei (Ascomycota, Pezizomycotina) QM6a is a model fungus for a broad spectrum of physiological phenomena, including plant cell wall degradation, industrial production of enzymes, light responses, conidiation, sexual development, polyketide biosynthesis, and plant-fungal interactions. The genomes of QM6a and its high enzyme-producing mutants have been sequenced by second-generation-sequencing methods and are publicly available from the Joint Genome Institute. While these genome sequences have offered useful information for genomic and transcriptomic studies, their limitations and especially their short read lengths make them poorly suited for some particular biological problems, including assembly, genome-wide determination of chromosome architecture, and genetic modification or engineering.

RESULTS

We integrated Pacific Biosciences and Illumina sequencing platforms for the highest-quality genome assembly yet achieved, revealing seven telomere-to-telomere chromosomes (34,922,528 bp; 10877 genes) with 1630 newly predicted genes and >1.5 Mb of new sequences. Most new sequences are located on AT-rich blocks, including 7 centromeres, 14 subtelomeres, and 2329 interspersed AT-rich blocks. The seven QM6a centromeres separately consist of 24 conserved repeats and 37 putative centromere-encoded genes. These findings open up a new perspective for future centromere and chromosome architecture studies. Next, we demonstrate that sexual crossing readily induced cytosine-to-thymine point mutations on both tandem and unlinked duplicated sequences. We also show by bioinformatic analysis that T. reesei has evolved a robust repeat-induced point mutation (RIP) system to accumulate AT-rich sequences, with longer AT-rich blocks having more RIP mutations. The widespread distribution of AT-rich blocks correlates genome-wide partitions with gene clusters, explaining why clustering of genes has been reported to not influence gene expression in T. reesei.

CONCLUSION

Compartmentation of ancestral gene clusters by AT-rich blocks might promote flexibilities that are evolutionarily advantageous in this fungus' soil habitats and other natural environments. Our analyses, together with the complete genome sequence, provide a better blueprint for biotechnological and industrial applications.

Reconstructing the evolutionary history of gypsy retrotransposons in the Périgord black truffle (Tuber melanosporum Vittad.)

Truffles are ascomycete fungi belonging to genus Tuber, and they form ectomycorrhizal associations with trees and shrubs. Transposable elements constitute more than 50 % of the black Périgord truffle (Tuber melanosporum) genome, which are mainly class 1 gypsy retrotransposons, but their impact on its genome is unknown. The aims of this study are to investigate the diversity of gypsy retrotransposons in this species and their evolutionary history by analysing the reference genome and six resequenced genomes of different geographic accessions. Using the reverse transcriptase sequences, six different gypsy retrotransposon clades were identified. Tmt1 and Tmt6 are the most abundant transposable elements, representing 14 and 13 % of the T. melanosporum genome, respectively. Tmt6 showed a major burst of proliferation between 1 and 4 million years ago, but evidence of more recent transposition was observed. Except for Tmt2, the other clades tend to aggregate, and their mode of transposition excluded the master copy model. This suggests that each new copy has the same probability of transposing as other copies. This study provides a better view of the diversity and dynamic nature of gypsy retrotransposons in T. melanosporum. Even if the major gypsy retrotransposon bursts are old, some elements seem to have transposed recently, suggesting that they may continue to model the truffle genomes.

Gene Network Polymorphism Illuminates Loss and Retention of Novel RNAi Silencing Components in the Cryptococcus Pathogenic Species Complex

RNAi is a ubiquitous pathway that serves central functions throughout eukaryotes, including maintenance of genome stability and repression of transposon expression and movement. However, a number of organisms have lost their RNAi pathways, including the model yeast Saccharomyces cerevisiae, the maize pathogen Ustilago maydis, the human pathogen Cryptococcus deuterogattii, and some human parasite pathogens, suggesting there may be adaptive benefits associated with both retention and loss of RNAi. By comparing the RNAi-deficient genome of the Pacific Northwest Outbreak C. deuterogattii strain R265 with the RNAi-proficient genomes of the Cryptococcus pathogenic species complex, we identified a set of conserved genes that were lost in R265 and all other C. deuterogattii isolates examined. Genetic and molecular analyses reveal several of these lost genes play roles in RNAi pathways. Four novel components were examined further. Znf3 (a zinc finger protein) and Qip1 (a homolog of N. crassa Qip) were found to be essential for RNAi, while Cpr2 (a constitutive pheromone receptor) and Fzc28 (a transcription factor) are involved in sex-induced but not mitosis-induced silencing. Our results demonstrate that the mitotic and sex-induced RNAi pathways rely on the same core components, but sex-induced silencing may be a more specific, highly induced variant that involves additional specialized or regulatory components. Our studies further illustrate how gene network polymorphisms involving known components of key cellular pathways can inform identification of novel elements and suggest that RNAi loss may have been a core event in the speciation of C. deuterogattii and possibly contributed to its pathogenic trajectory.

Genome instability and mutations provoked by transposon movement are counteracted by novel defense mechanisms in organisms as diverse as fungi, plants, and mammals. In the human fungal pathogen Cryptococcus neoformans, an RNAi silencing pathway operates to defend the genome against mobile elements and transgene repeats. RNAi silencing pathways are conserved in the Cryptococcus pathogenic species complex and are mediated by canonical RNAi components. Surprisingly, several of these components are missing from all analyzed C. deuterogattii VGII strains, the molecular type responsible for the North American Pacific Northwest outbreak. To identify novel components of the RNAi pathways, we surveyed the reference genomes of C. deuterogattii, C. gattii, C. neoformans, and C. deneoformans. We identified 14 otherwise conserved genes missing in R265, including the RDP1, AGO1, and DCR1 canonical RNAi components, and focused on four potentially novel RNAi components: ZNF3, QIP1, CPR2, and FZC28. We found that Znf3 and Qip1 are both required for mitotic- and sex-induced silencing, while Cpr2 and Fzc28 contribute to sex-induced but not mitosis-induced silencing. Our studies reveal elements of RNAi pathways that operate to defend the genome during sexual development and vegetative growth and illustrate the power of network polymorphisms to illuminate novel components of biological pathways.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.

Comparative Genomics of Serial Isolates of Cryptococcus neoformans Reveals Gene Associated With Carbon Utilization and Virulence

The opportunistic fungal pathogen Cryptococcus neoformans is a leading cause of mortality among the human immunodeficiency virus/acquired immunodeficiency syndrome population and is known for frequently causing life-threatening relapses. To investigate the potential contribution of in-host microevolution to persistence and relapse, we have analyzed two serial isolates obtained from a patient with acquired immunodeficiency syndrome who suffered an initial and relapse episode of cryptococcal meningoencephalitis. Despite being identical by multilocus sequence typing, the isolates differ phenotypically, exhibiting changes in key virulence factors, nutrient acquisition, metabolic profiles, and the ability to disseminate in an animal model. Whole-genome sequencing uncovered a clonal relationship, with only a few unique differences. Of these, two key changes are expected to explain the phenotypic differences observed in the relapse isolate: loss of a predicted AT-rich interaction domain protein and changes in copy number of the left and right arms of chromosome 12. Gene deletion of the predicted transcriptional regulator produced changes in melanin, capsule, carbon source use, and dissemination in the host, consistent with the phenotype of the relapse isolate. In addition, the deletion mutant displayed altered virulence in the murine model. The observed differences suggest the relapse isolate evolved subsequent to penetration of the central nervous system and may have gained dominance following the administration of antifungal therapy. These data reveal the first molecular insights into how the Cryptococcus neoformans genome changes during infection of humans and the manner in which microevolution progresses in this deadly fungal pathogen.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

Characterizing the role of RNA silencing components in Cryptococcus neoformans

The RNA interference (RNAi) mediated by homology-dependent degradation of the target mRNA with small RNA molecules plays a key role in controlling transcription and translation processes in a number of eukaryotic organisms. The RNAi machinery is also evolutionarily conserved in a wide variety of fungal species, including pathogenic fungi. To elucidate the physiological functions of the RNAi pathway in Cryptococcus neoformans that causes fungal meningitis, here we performed genetic analyses for genes encoding Argonaute (AGO1 and AGO2), RNA-dependent RNA polymerase (RDP1), and Dicers (DCR1 and DCR2) in both serotype A and D C. neoformans. The present study shows that Ago1, Rdp1, and Dcr2 are the major components of the RNAi process occurring in C. neoformans. However, the RNAi machinery is not involved in regulation of production of two virulence factors (capsule and melanin), sexual differentiation, and diverse stress response. Comparative transcriptome analysis of the serotype A and D RNAi mutants revealed that only modest changes occur in the genome-wide transcriptome profiles when the RNAi process was perturbed. Notably, the serotype D rdp1Δ mutants showed an increase in transcript abundance of active retrotransposons and transposons, such as T2 and T3, the latter of which is a novel serotype D-specific transposon of C. neoformans. In a wild type background both T2 and T3 were found to be weakly active mobile elements, although we found no evidence of Cnl1 retrotransposon mobility. In contrast, all three transposable elements exhibited enhanced mobility in the rdp1Δ mutant strain. In conclusion, the RNAi pathway plays an important role in controlling transposon activity and genome integrity of C. neoformans.

Evolutionary genomics revealed interkingdom distribution of Tcn1-like chromodomain-containing Gypsy LTR retrotransposons among fungi and plants

Background

Chromodomain-containing Gypsy LTR retrotransposons or chromoviruses are widely distributed among eukaryotes and have been found in plants, fungi and vertebrates. The previous comprehensive survey of chromoviruses from mosses (Bryophyta) suggested that genomes of non-seed plants contain the clade which is closely related to the retrotransposons from fungi. The origin, distribution and evolutionary history of this clade remained unclear mainly due to the absence of information concerning the diversity and distribution of LTR retrotransposons in other groups of non-seed plants as well as in fungal genomes.

Results

In present study we preformed in silico analysis of chromodomain-containing LTR retrotransposons in 25 diverse fungi and a number of plant species including spikemoss Selaginella moellendorffii (Lycopodiophyta) coupled with an experimental survey of chromodomain-containing Gypsy LTR retrotransposons from diverse non-seed vascular plants (lycophytes, ferns, and horsetails). Our mining of Gypsy LTR retrotransposons in genomic sequences allowed identification of numerous families which have not been described previously in fungi. Two new well-supported clades, Galahad and Mordred, as well as several other previously unknown lineages of chromodomain-containing Gypsy LTR retrotransposons were described based on the results of PCR-mediated survey of LTR retrotransposon fragments from ferns, horsetails and lycophytes. It appeared that one of the clades, namely Tcn1 clade, was present in basidiomycetes and non-seed plants including mosses (Bryophyta) and lycophytes (genus Selaginella).

Conclusions

The interkingdom distribution is not typical for chromodomain-containing LTR retrotransposons clades which are usually very specific for a particular taxonomic group. Tcn1-like LTR retrotransposons from fungi and non-seed plants demonstrated high similarity to each other which can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission.

Most environmental isolates of Cryptococcus neoformans var. grubii (serotype A) are not lethal for mice

Most cases of cryptococcosis are caused by Cryptococcus neoformans var. grubii (serotype A), which is widespread in the environment, where it is primarily associated with pigeon excreta. A number of molecular epidemiological studies indicate that many environmental and clinical isolates of serotype A are indistinguishable. However, the murine virulence of environmental strains of C. neoformans has not been thoroughly evaluated. We used the murine intranasal model of cryptococcosis to compare the lethality of clinical and environmental strains of serotype A that possessed identical genotypes as determined by amplified fragment length polymorphisms (AFLP) and multilocus sequence typing (MLST). Eleven environmental strains were tested, and only one caused disease within 60 days postinfection, at which time the experiments were terminated. Conversely, 7 of 10 clinical isolates were lethal for mice at median times of 19 to 40 days. Passing environmental isolates in mice (up to three times) did not significantly increase their lethality. In follow-up studies, we developed a new genotyping technique based on hybridization with TCN2 and TCN4 retrotransposon-specific probes. Although the retrotransposon banding patterns were unstable after prolonged incubation in the laboratory, this method was able to differentiate clinical and environmental strains that had the same AFLP/MLST genotypes.

Chromosomal rearrangements between serotype A and D strains in Cryptococcus neoformans

Cryptococcus neoformans is a major human pathogenic fungus that can cause meningoencephalitis in immunocompromised hosts. It contains two divergent varieties, var. grubii (serotype A) and var. neoformans (serotype D), as well as hybrids (serotype AD) between these two varieties. In this study, we investigated the extent of chromosomal rearrangements between the two varieties, estimated the effects of chromosomal rearrangements on recombination frequencies, and surveyed the potential polymorphisms of the rearrangements among natural strains of the three serotypes. Through the analyses of two sequenced genomes from strains H99 (representing var. grubii) and JEC21 (representing var. neoformans), we revealed a total of 32 unambiguous chromosome rearrangements, including five translocations, nine simple inversions, and 18 complex rearrangements. Our analyses identified that overall, rearranged regions had recombination frequencies about half of those around syntenic regions. Using a direct PCR screening strategy, we examined the potential polymorphisms of 11 rearrangements among 64 natural C. neoformans strains from five countries. We found no polymorphism within var. neoformans and very limited polymorphism within var. grubii. However, strains of serotype AD showed significant polymorphism, consistent with their hybrid origins coupled with differential loss of heterozygosity. We discuss the implications of these results on the genome structure, ecology, and evolution of C. neoformans.

LINE-like retrotransposition in Saccharomyces cerevisiae

Over one-third of human genome sequence is a product of non-LTR retrotransposition. The retrotransposon that currently drives this process in humans is the highly abundant LINE-1 (L1) element. Despite the ubiquitous nature of L1's in mammals, we still lack a complete mechanistic understanding of the L1 replication cycle and how it is regulated. To generate a genetically amenable model for non-LTR retrotransposition, we have reengineered the Zorro3 retrotransposon, an L1 homolog from Candida albicans, for use in the budding yeast Saccharomyces cerevisiae. We found that S. cerevisiae, which has no endogenous L1 homologs or remnants, can still support Zorro3 retrotransposition. Analysis of Zorro3 mutants and insertion structures suggest that this is authentic L1-like retrotransposition with remarkable resemblance to mammalian L1-mediated events. This suggests that S. cerevisiae has unexpectedly retained the basal host machinery required for L1 retrotransposition. This model will also serve as a powerful system to study the cell biology of L1 elements and for the genetic identification and characterization of cellular factors involved in L1 retrotransposition.

A retrotransposon-derived probe for discriminating strains of Cryptococcus neoformans

Hybridization of digested DNA to probes derived from repeated sequences has proven to be an extremely powerful epidemiologic tool for studying the relatedness of fungi. The dispersed nature of these sequences throughout the genome provides the discriminatory power for distinguishing two independent isolates from each other based on banding pattern. The genome of Cryptococcus neoformans contains a number of classes of transposable elements, which are often present in multiple copies. We characterized a probe related to the Ty3/gypsy class of transposable elements called TCN1 and used it to screen multiple isolates from all four serotypes of C. neoformans. DNA with TCN1 homology could be amplified from each isolate of serotypes A and D and all isolates hybridized to a probe derived from TCN1. Isolates from serotype B and C were also tested for the presence of a TCN1 homolog, however, only some of these isolates yielded both a TCN1-specific PCR product or hybridization signal. Comparison of the TCN1 hybridization patterns of serotypes A and D to multiple RAPD patterns of the same isolates suggested that TCN1 was more discriminating and therefore, a useful epidemiological tool.

Characterization of AFLAV, a Tf1/Sushi retrotransposon from Aspergillus flavus

The plasmid, pAF28, a genomic clone from Aspergillus flavus NRRL 6541, has been used as a hybridization probe to fingerprint A. flavus strains isolated in corn and peanut fields. The insert of pAF28 contains a 4.5 kb region which encodes a truncated retrotransposon (AfRTL-1). In search for a full-length and intact copy of retrotransposon, we exploited a novel PCR cloning strategy by amplifying a 3.4 kb region from the genomic DNA of A. flavus NRRL 6541. The fragment was cloned into pCR 4-TOPO. Sequence analysis confirmed that this region encoded putative domains of partial reverse transcriptase, RNase H, and integrase of the predicted retrotransposon. The two flanking long terminal repeats (LTRs) and the sequence between them comprise a putative full-length LTR retrotransposon of 7799 bp in length. This intact retrotransposon sequence is named AFLAV (A. flavus Retrotransposon). The order of the predicted catalytic domains in the polyprotein (Pol) placed AFLAV in the Tf1/sushi subgroup of the Ty3/gypsy retrotransposon family. Primers derived from AFLAV sequence were used to screen this retrotransposon in other strains of A. flavus. More than fifty strains of A. flavus isolated from different geological origins were surveyed and the results show that many strains have extensive deletions in the regions encoding the capsid (Gag) and Pol.

The biology of the Cryptococcus neoformans species complex

Cryptococcus neoformans is a major cause of fungal meningoencephalitis in immunocompromised patients. Despite recent advances in the genetics and molecular biology of C. neoformans, and improved techniques for molecular epidemiology, aspects of the ecology, population structure, and mode of reproduction of this environmental pathogen remain to be established. Application of recent insights into the life cycle of C. neoformans and its different ways of engaging in sexual reproduction under laboratory conditions has just begun to affect research on the ecology and epidemiology of this human pathogenic fungus. The melding of these disparate disciplines should yield rich dividends in our understanding of the evolution of microbial pathogens, providing insights relevant to diagnosis, treatment, and prevention.

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.

Tyl6, a novel Ty3/gypsy-like retrotransposon in the genome of the dimorphic fungus Yarrowia lipolytica

The novel LTR retrotransposon Tyl6 was detected in the genome of the dimorphic fungus Yarrowia lipolytica. Sequence analysis revealed that this element is related to the well-known Ty3 element of Saccharomyces cerevisiae and, especially, to the recently described Tse3 retrotransposon of Saccharomyces exiguus and to the del1-like plant retrotransposons. Tyl6 is 5108 bp long, is flanked by two identical long terminal repeats (LTR), each of 276 bp, and its ORFs are separated by a -1 frameshift. Both ORFs are intact and deduced translation products display a significant similarity with those of previously described Ty3/gypsy retrotransposons. Distribution of Tyl6 among Y. lipolytica strains of different origins was also analysed. A single copy of the novel retrotransposon is present in some commonly used laboratory strains, which are derivatives of the wild-type isolate YB423-12, whereas other strains of independent origin are devoid of Ty16. No solo LTR of Tyl6 was detected in the analysed strains.

Deciphering the model pathogenic fungus Cryptococcus neoformans

Cryptococcus neoformans is a basidiomycete fungal pathogen of humans that has diverged considerably from other model fungi such as Neurospora crassa, Aspergillus nidulans, Saccharomyces cerevisiae and the common human fungal pathogen Candida albicans. The recent completion of the genome sequences of two related C. neoformans strains and the ongoing genome sequencing of three other divergent Cryptococcus strains with different virulence phenotypes and environmental distributions should improve our understanding of this important pathogen. We discuss the biology of C. neoformans in light of this genomic data, with a special emphasis on the role that evolution and sexual reproduction have in the complex relationships of the fungus with the environment and the host.

Mobile genetic elements and sexual reproduction

Transposable elements (TE) are prominent components of most eukaryotic genomes. In addition to their possible participation in the origin of sexual reproduction in eukaryotes, they may be also involved in its maintenance as important contributors to the deleterious mutation load. Comparative analyses of transposon content in the genomes of sexually reproducing and anciently asexual species may help to understand the contribution of different TE classes to the deleterious load. The apparent absence of deleterious retrotransposons from the genomes of ancient asexuals is in agreement with the hypothesis that they may play a special role in the maintenance of sexual reproduction and in early extinction for which most species are destined upon the abandonment of sex.

Phylogenomic analysis of chromoviruses

Genome sequences of model organisms provide a unique opportunity to obtain insight into the complete diversity of any transposable element (TE) group. A limited number of chromoviruses, the chromodomain containing genus of Metaviridae, is known from plant, fungal and vertebrate genomes. By searching diverse eukaryotic genome databases, we have found a surprisingly large number of new, structurally intact and highly conserved chromoviral elements, greatly exceeding the number of previously known chromoviruses. In this study, we examined the diversity, origin and evolution of chromoviruses in Eukaryota. Chromoviral diversity in plants, fungi and vertebrates, as shown by phylogenetic analyses, was found to be much greater than previously expected. A novel centromere-specific chromoviral lineage was found to be widespread and highly conserved in all seed plants. The age of chromoviruses has been significantly extended by finding their representatives in the most basal plant lineages (green and red algae), in Heterokonta (oomycetes) and in Cercozoa (plasmodiophorids). The evolutionary origin of chromoviruses has been found to be no earlier than in Cercozoa, since none can be found in the basal eukaryotic lineages, despite the extensive genome data. The evolutionary dynamics of chromoviruses can be explained by a strict vertical transmission in plants and fungi, while in Metazoa it is more complex. The currently available genome data clearly show that chromoviruses are the most widespread and one of the oldest Metaviridae clade.

Repeat-induced point mutation and the population structure of transposable elements in Microbotryum violaceum

Repeat-induced point mutation (RIP) is a genome defense in fungi that hypermutates repetitive DNA and is suggested to limit the accumulation of transposable elements. The genome of Microbotryum violaceum has a high density of transposable elements compared to other fungi, but there is also evidence of RIP activity. This is the first report of RIP in a basidiomycete and was obtained by sequencing multiple copies of the integrase gene of a copia-type transposable element and the helicase gene of a Helitron-type element. In M. violaceum, the targets for RIP mutations are the cytosine residues of TCG trinucleotide combinations. Although RIP is a linkage-dependent process that tends to increase the variation among repetitive sequences, a chromosome-specific substructuring was observed in the transposable element population. The observed chromosome-specific patterns are not consistent with RIP, but rather suggest an effect of gene conversion, which is also a linkage-dependent process but results in a homogenization of repeated sequences. Particular sequences were found more widely distributed within the genome than expected by chance and may reflect the recently active variants. Therefore, sequence variation of transposable elements in M. violaceum appears to be driven by selection for transposition ability in combination with the context-specific forces of the RIP and gene conversion.

Repetitive DNA in the automictic fungus Microbotryum violaceum

The small genomes of fungi are expected to have little repetitive content other than rDNA genes. Moreover, among asexual or highly selfing lineages, the diversity of repetitive elements is also expected to be very low. However, in the automictic fungus Microbotryum violaceum, a very large proportion of random DNA fragments from the autosomes and the fungal sex chromosomes are repetitive in nature, either as retrotransposon or helicase sequences. Among the retrotransposon sequences, examples were found from each major kind of elements, including copia, gypsy, and non-LTR sequences. The most numerous were copia-like elements, which are believed to be rare in fungi, particularly among basidiomycetes. The many helicase sequences appear to belong to the recently discovered Helitron type of transposable elements. Also, sequences that could not be identified as a known type of gene were also very repetitive within the database of random fragments from M. violaceum. The differentiated pair of fungal sex chromosomes and suppression of recombination may be the major forces determining the highly repetitive content in the small genome of M. violaceum.

Mutator-like element in the yeast Yarrowia lipolytica displays multiple alternative splicings

A new type of DNA transposon, Mutyl, has been identified in the sequenced genome of the yeast Yarrowia lipolytica. This transposon is 7,413 bp long and carries two open reading frames (ORFs) which potentially encode proteins of 459 and 1,178 amino acids, respectively. Whereas the first ORF shows no significant homology to previously described proteins, the second ORF shows sequence similarities with various Mutator-like element (MULE)-encoded transposases, including the bacterial transposase signature sequence. Other MULE features shared by Mutyl include a zinc finger motif in the putative transposase, a 22-bp-long imperfect inverted repeat at each end, and a 9- to 10-bp duplication of its target site in the chromosome. Of the five copies of Mutyl present in the genome, one has a deletion of the first 8 bases, and the others are full length with a single base change in one element. The first potential gene of Mutyl, mutB, was shown to be expressed in exponentially growing cells. Its sequence contains a predicted intron with two 5' splice sites, a single branch point, and two 3' splice sites. Its mRNA is alternatively spliced, as judged by reverse transcription-PCR, and generates four mRNAs corresponding to protein-coding sequences of 128, 156, 161, and 190 amino acids. Of the three distinct lineages characterized in Y. lipolytica, strains from the German lineage and the French lineage do not carry Mutyl. A study of the distribution of Mutyl in strains of the French lineage evidenced a recent transposition event. Taken together, these results indicate that Mutyl is still active.

Chromosomal translocation and segmental duplication in Cryptococcus neoformans

Large chromosomal events such as translocations and segmental duplications enable rapid adaptation to new environments. Here we marshal genomic, genetic, meiotic mapping, and physical evidence to demonstrate that a chromosomal translocation and segmental duplication occurred during construction of a congenic strain pair in the fungal human pathogen Cryptococcus neoformans. Two chromosomes underwent telomere-telomere fusion, generating a dicentric chromosome that broke to produce a chromosomal translocation, forming two novel chromosomes sharing a large segmental duplication. The duplication spans 62,872 identical nucleotides and generated a second copy of 22 predicted genes, and we hypothesize that this event may have occurred during meiosis. Gene disruption studies of one embedded gene (SMG1) corroborate that this region is duplicated in an otherwise haploid genome. These findings resolve a genome project assembly anomaly and illustrate an example of rapid genome evolution in a fungal genome rich in repetitive elements.

Cryptococcus neoformans virulence gene discovery through insertional mutagenesis

Insertional mutagenesis was applied to Cryptococcus neoformans to identify genes associated with virulence attributes. Using biolistic transformation, we generated 4,300 nourseothricin (NAT)-resistant strains, of which 590 exhibited stable resistance. We focused on mutants with defects in established virulence factors and identified two with reduced growth at 37 degrees C, four with reduced production of the antioxidant pigment melanin, and two with an increased sensitivity to nitric oxide (NO). The NAT insertion and mutant phenotypes were genetically linked in five of eight mutants, and the DNA flanking the insertions was characterized. For the strains with altered growth at 37 degrees C and altered melanin production, mutations were in previously uncharacterized genes, while the two NO-sensitive strains bore insertions in the flavohemoglobin gene FHB1, whose product counters NO stress. Because of the frequent instability of nourseothricin resistance associated with biolistic transformation, Agrobacterium-mediated transformation was tested. This transkingdom DNA delivery approach produced 100% stable nourseothricin-resistant transformants, and three melanin-defective strains were identified from 576 transformants, of which 2 were linked to NAT in segregation analysis. One of these mutants contained a T-DNA insertion in the promoter of the LAC1 (laccase) gene, which encodes a key enzyme required for melanin production, while the second contained an insertion in the promoter of the CLC1 gene, encoding a voltage-gated chloride channel. Clc1 and its homologs are required for ion homeostasis, and in their absence Cu+ transport into the secretory pathway is compromised, depriving laccase and other Cu(+)-dependent proteins of their essential cofactor. The NAT resistance cassette was optimized for cryptococcal codon usage and GC content and was then used to disrupt a mitogen-activated protein kinase gene, a predicted gene, and two putative chloride channel genes to analyze their contributions to fungal physiology. Our findings demonstrate that both insertional mutagenesis methods can be applied to gene identification, but Agrobacterium-mediated transformation is more efficient and generates exclusively stable insertion mutations.

Convergent evolution of chromosomal sex-determining regions in the animal and fungal kingdoms

Sexual identity is governed by sex chromosomes in plants and animals, and by mating type (MAT) loci in fungi. Comparative analysis of the MAT locus from a species cluster of the human fungal pathogen Cryptococcus revealed sequential evolutionary events that fashioned this large, highly unusual region. We hypothesize that MAT evolved via four main steps, beginning with acquisition of genes into two unlinked sex-determining regions, forming independent gene clusters that then fused via chromosomal translocation. A transitional tripolar intermediate state then converted to a bipolar system via gene conversion or recombination between the linked and unlinked sex-determining regions. MAT was subsequently subjected to intra- and interallelic gene conversion and inversions that suppress recombination. These events resemble those that shaped mammalian sex chromosomes, illustrating convergent evolution in sex-determining structures in the animal and fungal kingdoms.

A comparative genomic analysis of the sex determining region in fungi reveals a remarkable similarity between its evolution and the events which shaped mammalian sex chromosomes

A comparison of the nature and abundance of microsatellites in 14 fungal genomes

An overview of the character of microsatellites in 14 fungal genomes was obtained by analyzing databases containing complete or nearly complete genome sequences. Low GC content, rather than genome size, was the best predictor of high microsatellite density, although very long iterations of tandem repeats were less common in small genomes. Motif type correlated with %GC in that low-GC genomes were more likely to be dominated by A/T-rich motifs, and vice versa, although some exceptions were noted. The experimentally useful dinucleotide and trinucleotide arrays were analyzed in greater detail. Although these varied in sequence and length among fungal species, some that are likely to be universally useful were identified. This information will be useful for researchers wanting to identify the most useful microsatellites to analyze for the fungi included in this survey and provides a platform for choosing microsatellites to target in fungi that are not yet sequenced.

Cryptons: a group of tyrosine-recombinase-encoding DNA transposons from pathogenic fungi

A new group of transposable elements, which the authors have named cryptons, was detected in several pathogenic fungi, including the basidiomyceteCryptococcus neoformans, and the ascomycetesCoccidioides posadasiiandHistoplasma capsulatum. These elements are unlike any previously described transposons. An archetypal member of the group, cryptonCn1, is 4 kb in length and is present at a low but variable copy number in a variety ofC. neoformansstrains. It displays interstrain variations in its insertion sites, suggesting recent mobility. The internal region contains a long gene, interrupted by several introns. The product of this gene contains a putative tyrosine recombinase near its middle, and a region similar in sequence to the DNA-binding domains of several fungal transcription factors near its C-terminus. The element contains no long repeat sequences, but is bordered by short direct repeats which may have been produced by its insertion into the host genome by recombination. Many of the structural features of cryptonCn1are conserved in the other known cryptons, suggesting that these elements represent the functional forms. The presence of cryptons in ascomycetes and basidiomycetes suggests that this is an ancient group of elements (>400 million years old). Sequence comparisons suggest that cryptons may be related to the DIRS1 andNgaro1groups of tyrosine-recombinase-encoding retrotransposons.

Mating-type locus of Cryptococcus neoformans: a step in the evolution of sex chromosomes

The sexual development and virulence of the fungal pathogen Cryptococcus neoformans is controlled by a bipolar mating system determined by a single locus that exists in two alleles, alpha and a. The alpha and a mating-type alleles from two divergent varieties were cloned and sequenced. The C. neoformans mating-type locus is unique, spans >100 kb, and contains more than 20 genes. MAT-encoded products include homologs of regulators of sexual development in other fungi, pheromone and pheromone receptors, divergent components of a MAP kinase cascade, and other proteins with no obvious function in mating. The alpha and a alleles of the mating-type locus have extensively rearranged during evolution and strain divergence but are stable during genetic crosses and in the population. The C. neoformans mating-type locus is strikingly different from the other known fungal mating-type loci, sharing features with the self-incompatibility systems and sex chromosomes of algae, plants, and animals. Our study establishes a new paradigm for mating-type loci in fungi with implications for the evolution of cell identity and self/nonself recognition.

SINEs and LINEs: the art of biting the hand that feeds you

SINEs and LINEs are short and long interspersed retrotransposable elements, respectively, that invade new genomic sites using RNA intermediates. SINEs and LINEs are found in almost all eukaryotes (although not in Saccharomyces cerevisiae) and together account for at least 34% of the human genome. The noncoding SINEs depend on reverse transcriptase and endonuclease functions encoded by partner LINEs. With the completion of many genome sequences, including our own, the database of SINEs and LINEs has taken a great leap forward. The new data pose new questions that can only be answered by detailed studies of the mechanism of retroposition. Current work ranges from the biochemistry of reverse transcription and integration invitro, target site selection in vivo, nucleocytoplasmic transport of the RNA and ribonucleoprotein intermediates, and mechanisms of genomic turnover. Two particularly exciting new ideas are that SINEs may help cells survive physiological stress, and that the evolution of SINEs and LINEs has been shaped by the forces of RNA interference. Taken together, these studies promise to explain the birth and death of SINEs and LINEs, and the contribution of these repetitive sequence families to the evolution of genomes.

Genomic evolution of the long terminal repeat retrotransposons in hemiascomycetous yeasts

We identified putative long terminal repeat- (LTR) retrotransposon sequences among the 50,000 random sequence tags (RSTs) obtained by the Génolevures project from genomic libraries of 13 Hemiascomycetes species. In most cases additional sequencing enabled us to assemble the whole sequences of these retrotransposons. These approaches identified 17 distinct families, 10 of which are defined by full-length elements. We also identified five families of solo LTRs that were not associated with retrotransposons. Ty1-like retrotransposons were found in four of five species that are phylogenetically related to Saccharomyces cerevisiae (S. uvarum, S. exiguus, S. servazzii, and S. kluyveri but not Zygosaccharomyces rouxii), and in two of three Kluyveromyces species (K. lactis and K. marxianus but not K. thermotolerans). Only multiply crippled elements could be identified in the K. lactis and S. servazzii strains analyzed, and only solo LTRs could be identified in S. uvarum. Ty4-like elements were only detected in S. uvarum, indicating that these elements appeared recently before speciation of the Saccharomyces sensu stricto species. Ty5-like elements were detected in S. exiguus, Pichia angusta, and Debaryomyces hansenii. A retrotransposon homologous with Tca2 from Candida albicans, an element absent from S. cerevisiae, was detected in the closely related species D. hansenii. A complete Ty3/gypsy element was present in S. exiguus, whereas only partial, often degenerate, sequences resembling this element were found in S. servazzii, Z. rouxii, S. kluyveri, C. tropicalis, and Yarrowica lipolytica. P. farinosa (syn. P. sorbitophila) is currently the only yeast species in which no LTR retrotransposons or remnants have been found. Thorough analysis of protein sequences, structural characteristics of the elements, and phylogenetic relationships deduced from these data allowed us to propose a classification for the Ty1/copia elements of hemiascomycetous yeasts and a model of LTR-retrotransposon evolution in yeasts.

Ylli, a non-LTR retrotransposon L1 family in the dimorphic yeast Yarrowia lipolytica

During the course of a random sequencing project of the genome of the dimorphic yeast Yarrowia lipolytica, we have identified sequences that were repeated in the genome and that matched the reverse transcriptase (RT) sequence of non-long terminal repeat (non-LTR) retrotransposons. Extension of sequencing on each side of this zone of homology allowed the definition of an element over 6 kb long. The conceptual translation of this sequence revealed two open reading frames (ORFs) that displayed several characteristics of non-LTR retrotransposons: a Cys-rich motif in the ORF1, an N-terminal endonuclease, a central RT, and a C-terminal zinc finger domain in the ORF2. We called this element Ylli (for Y. lipolytica LINE). A total of 19 distinct repeats carrying the 3' untranslated region (UTR) and all ending with a poly-A tail were detected. Most of them were very short, 17 being 134 bp long or less. The number of copies of Ylli was estimated to be around 100 if these short repeats are 5' truncations. No 5' UTR was clearly identified, indicating that entire and therefore active elements might be very rare in the Y. lipolytica strain tested. Ylli does not seem to have any insertion specificity. Phylogenetic analysis of the RT domain unambiguously placed Ylli within the L1 clade. It forms a monophyletic group with the Zorro non-LTR retrotransposons discovered in another dimorphic yeast Candida albicans. BLAST comparisons showed that ORF2 of Ylli is closely related to that of the slime mold Dictyostelium discoideum L1 family, TRE.

Ancient lineages of non-LTR retrotransposons in the primitive eukaryote, Giardia lamblia

Mobile elements that use reverse transcriptase to make new copies of themselves are found in all major lineages of eukaryotes. The non-long terminal repeat (non-LTR) retrotransposons have been suggested to be the oldest of these eukaryotic elements. Phylogenetic analysis of non-LTR elements suggests that they have predominantly undergone vertical transmission, as opposed to the frequent horizontal transmissions found for other mobile elements. One prediction of this vertical model of inheritance is that the oldest lineages of eukaryotes should exclusively harbor the oldest lineages of non-LTR retrotransposons. Here we characterize the non-LTR retrotransposons present in one of the most primitive eukaryotes, the diplomonad Giardia lamblia. Two families of elements were detected in the WB isolate of G. lamblia currently being used for the genome sequencing project. These elements are clearly distinct from all other previously described non-LTR lineages. Phylogenetic analysis indicates that these Genie elements (for Giardia early non-LTR insertion element) are among the oldest known lineages of non-LTR elements consistent with strict vertical descent. Genie elements encode a single open reading frame with a carboxyl terminal endonuclease domain. Genie 1 is site specific, as seven to eight copies are present in a single tandem array of a 771-bp repeat near the telomere of one chromosome. The function of this repeat is not known. One additional, highly divergent, element within the Genie 1 lineage is not located in this tandem array but is near a second telomere. Four different telomere addition sites could be identified within or near the Genie elements on each of these chromosomes. The second lineage of non-LTR elements, Genie 2, is composed of about 10 degenerate copies. Genie 2 elements do not appear to be site specific in their insertion. An unusual aspect of Genie 2 is that all copies contain inverted repeats up to 172 bp in length.