Multiple LTR-retrotransposon families in the asexual yeast Candida albicans

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.

Intraspecies Transcriptional Profiling Reveals Key Regulators of Candida albicans Pathogenic Traits

The human commensal and opportunistic fungal pathogen Candida albicans displays extensive genetic and phenotypic variation across clinical isolates. Here, we performed RNA sequencing on 21 well-characterized isolates to examine how genetic variation contributes to gene expression differences and to link these differences to phenotypic traits. C. albicans adapts primarily through clonal evolution, and yet hierarchical clustering of gene expression profiles in this set of isolates did not reproduce their phylogenetic relationship. Strikingly, strain-specific gene expression was prevalent in some strain backgrounds. Association of gene expression with phenotypic data by differential analysis, linear correlation, and assembly of gene networks connected both previously characterized and novel genes with 23 C. albicans traits. Construction of de novo gene modules produced a gene atlas incorporating 67% of C. albicans genes and revealed correlations between expression modules and important phenotypes such as systemic virulence. Furthermore, targeted investigation of two modules that have novel roles in growth and filamentation supported our bioinformatic predictions. Together, these studies reveal widespread transcriptional variation across C. albicans isolates and identify genetic and epigenetic links to phenotypic variation based on coexpression network analysis.IMPORTANCE Infectious fungal species are often treated uniformly despite clear evidence of genotypic and phenotypic heterogeneity being widespread across strains. Identifying the genetic basis for this phenotypic diversity is extremely challenging because of the tens or hundreds of thousands of variants that may distinguish two strains. Here, we use transcriptional profiling to determine differences in gene expression that can be linked to phenotypic variation among a set of 21 Candida albicans isolates. Analysis of this transcriptional data set uncovered clear trends in gene expression characteristics for this species and new genes and pathways that were associated with variation in pathogenic processes. Direct investigation confirmed functional predictions for a number of new regulators associated with growth and filamentation, demonstrating the utility of these approaches in linking genes to important phenotypes.

Twenty years of transposable element analysis in the Arabidopsis thaliana genome

Transposable elements (TEs) are mobile repetitive DNA sequences shown to be major drivers of genome evolution. As the first plant to have its genome sequenced and analyzed at the genomic scale, Arabidopsis thaliana has largely contributed to our TE knowledge.

The present report describes 20 years of accumulated TE knowledge gained through the study of the Arabidopsis genome and covers the known TE families, their relative abundance, and their genomic distribution. It presents our knowledge of the different TE family activities, mobility, population and long-term evolutionary dynamics. Finally, the role of TE as substrates for new genes and their impact on gene expression is illustrated through a few selected demonstrative cases. Promising future directions for TE studies in this species conclude the review.

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.

To Repeat or Not to Repeat: Repetitive Sequences Regulate Genome Stability in Candida albicans

Genome instability often leads to cell death but can also give rise to innovative genotypic and phenotypic variation through mutation and structural rearrangements. Repetitive sequences and chromatin architecture in particular are critical modulators of recombination and mutability. In Candida albicans, four major classes of repeats exist in the genome: telomeres, subtelomeres, the major repeat sequence (MRS), and the ribosomal DNA (rDNA) locus. Characterization of these loci has revealed how their structure contributes to recombination and either promotes or restricts sequence evolution. The mechanisms of recombination that give rise to genome instability are known for some of these regions, whereas others are generally unexplored. More recent work has revealed additional repetitive elements, including expanded gene families and centromeric repeats that facilitate recombination and genetic innovation. Together, the repeats facilitate C. albicans evolution through construction of novel genotypes that underlie C. albicans adaptive potential and promote persistence across its human host.

Chromatin-Mediated Regulation of Genome Plasticity in Human Fungal Pathogens

Human fungal pathogens, such as Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans, are a public health problem, causing millions of infections and killing almost half a million people annually. The ability of these pathogens to colonise almost every organ in the human body and cause life-threating infections relies on their capacity to adapt and thrive in diverse hostile host-niche environments. Stress-induced genome instability is a key adaptive strategy used by human fungal pathogens as it increases genetic diversity, thereby allowing selection of genotype(s) better adapted to a new environment. Heterochromatin represses gene expression and deleterious recombination and could play a key role in modulating genome stability in response to environmental changes. However, very little is known about heterochromatin structure and function in human fungal pathogens. In this review, I use our knowledge of heterochromatin structure and function in fungal model systems as a road map to review the role of heterochromatin in regulating genome plasticity in the most common human fungal pathogens: Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans.

Chromatin Profiling of the Repetitive and Nonrepetitive Genomes of the Human Fungal Pathogen Candida albicans

The fungus Candida albicans is an opportunistic pathogen that normally lives on the human body without causing any harm. However, C. albicans is also a dangerous pathogen responsible for millions of infections annually. C. albicans is such a successful pathogen because it can adapt to and thrive in different environments. Chemical modifications of chromatin, the structure that packages DNA into cells, can allow environmental adaptation by regulating gene expression and genome organization. Surprisingly, the contribution of chromatin modification to C. albicans biology is still largely unknown. For the first time, we analyzed C. albicans chromatin modifications on a genome-wide basis. We demonstrate that specific chromatin states are associated with distinct regions of the C. albicans genome and identify the roles of the chromatin modifiers Sir2 and Set1 in shaping C. albicans chromatin and gene expression.

Eukaryotic genomes are packaged into chromatin structures that play pivotal roles in regulating all DNA-associated processes. Histone posttranslational modifications modulate chromatin structure and function, leading to rapid regulation of gene expression and genome stability, key steps in environmental adaptation. Candida albicans, a prevalent fungal pathogen in humans, can rapidly adapt and thrive in diverse host niches. The contribution of chromatin to C. albicans biology is largely unexplored. Here, we generated the first comprehensive chromatin profile of histone modifications (histone H3 trimethylated on lysine 4 [H3K4me³], histone H3 acetylated on lysine 9 [H3K9Ac], acetylated lysine 16 on histone H4 [H4K16Ac], and γH2A) across the C. albicans genome and investigated its relationship to gene expression by harnessing genome-wide sequencing approaches. We demonstrated that gene-rich nonrepetitive regions are packaged into canonical euchromatin in association with histone modifications that mirror their transcriptional activity. In contrast, repetitive regions are assembled into distinct chromatin states; subtelomeric regions and the ribosomal DNA (rDNA) locus are assembled into heterochromatin, while major repeat sequences and transposons are packaged in chromatin that bears features of euchromatin and heterochromatin. Genome-wide mapping of γH2A, a marker of genome instability, identified potential recombination-prone genomic loci. Finally, we present the first quantitative chromatin profiling in C. albicans to delineate the role of the chromatin modifiers Sir2 and Set1 in controlling chromatin structure and gene expression. This report presents the first genome-wide chromatin profiling of histone modifications associated with the C. albicans genome. These epigenomic maps provide an invaluable resource to understand the contribution of chromatin to C. albicans biology and identify aspects of C. albicans chromatin organization that differ from that of other yeasts.

Genome plasticity in Candida albicans is driven by long repeat sequences

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.

Activation of transposable elements and genetic instability during long-term culture of the human fungal pathogen Candida albicans

It has been repeatedly reported that transposable elements (TE) become active and/or mobile in the genomes of replicatively and stress-induced senescent mammalian cells. However, the biological role of senescence-associated transposon activation and its occurrence and relevance in other eukaryotic cells remain to be elucidated. In the present study, Candida albicans, a prevalent opportunistic fungal pathogen in humans, was used to analyze changes in gene copy number of selected TE, namely Cirt2, Moa and Cmut1 during long-term culture (up to 90 days). The effects of stress stimuli (fluconazole, hydrogen peroxide, hypochlorite) and ploidy state (haploid, diploid, tetraploid cells) were also considered. An increase in copy number of Cirt2 and Moa was the most accented in tetraploid cells after 90 days of culture that was accompanied by changes in karyotype patterns and slightly more limited growth rate compared to haploid and diploid cells. Stress stimuli did not potentiate TE activity. Elevation in chromosomal DNA breaks was also observed during long-term culture of cells of different ploidy, however this was not correlated with increased TE activity. Our results suggest that increased TE activity may promote genomic diversity and plasticity, and cellular heterogeneity during long-term culture of C. albicans cells.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

A single nucleotide polymorphism uncovers a novel function for the transcription factor Ace2 during Candida albicans hyphal development

Candida albicans is a major invasive fungal pathogen in humans. An important virulence factor is its ability to switch between the yeast and hyphal forms, and these filamentous forms are important in tissue penetration and invasion. A common feature for filamentous growth is the ability to inhibit cell separation after cytokinesis, although it is poorly understood how this process is regulated developmentally. In C. albicans, the formation of filaments during hyphal growth requires changes in septin ring dynamics. In this work, we studied the functional relationship between septins and the transcription factor Ace2, which controls the expression of enzymes that catalyze septum degradation. We found that alternative translation initiation produces two Ace2 isoforms. While full-length Ace2, Ace2L, influences septin dynamics in a transcription-independent manner in hyphal cells but not in yeast cells, the use of methionine-55 as the initiation codon gives rise to Ace2S, which functions as the nuclear transcription factor required for the expression of cell separation genes. Genetic evidence indicates that Ace2L influences the incorporation of the Sep7 septin to hyphal septin rings in order to avoid inappropriate activation of cell separation during filamentous growth. Interestingly, a natural single nucleotide polymorphism (SNP) present in the C. albicans WO-1 background and other C. albicans commensal and clinical isolates generates a stop codon in the ninth codon of Ace2L that mimics the phenotype of cells lacking Ace2L. Finally, we report that Ace2L and Ace2S interact with the NDR kinase Cbk1 and that impairing activity of this kinase results in a defect in septin dynamics similar to that of hyphal cells lacking Ace2L. Together, our findings identify Ace2L and the NDR kinase Cbk1 as new elements of the signaling system that modify septin ring dynamics in hyphae to allow cell-chain formation, a feature that appears to have evolved in specific C. albicans lineages.

Evolutionary dynamics of hAT DNA transposon families in Saccharomycetaceae

Transposable elements (TEs) are widespread in eukaryotes but uncommon in yeasts of the Saccharomycotina subphylum, in terms of both host species and genome fraction. The class II elements are especially scarce, but the hAT element Rover is a noteworthy exception that deserves further investigation. Here, we conducted a genome-wide analysis of hAT elements in 40 ascomycota. A novel family, Roamer, was found in three species, whereas Rover was detected in 15 preduplicated species from Kluyveromyces, Eremothecium, and Lachancea genera, with up to 41 copies per genome. Rover acquisition seems to have occurred by horizontal transfer in a common ancestor of these genera. The detection of remote Rover copies in Naumovozyma dairenensis and in the sole Saccharomyces cerevisiae strain AWRI1631, without synteny, suggests that two additional independent horizontal transfers took place toward these genomes. Such patchy distribution of elements prevents any anticipation of TE presence in incoming sequenced genomes, even closely related ones. The presence of both putative autonomous and defective Rover copies, as well as their diversification into five families, indicate particular dynamics of Rover elements in the Lachancea genus. Especially, we discovered the first miniature inverted-repeat transposable elements (MITEs) to be described in yeasts, together with their parental autonomous copies. Evidence of MITE insertion polymorphism among Lachancea waltii strains suggests their recent activity. Moreover, 40% of Rover copies appeared to be involved in chromosome rearrangements, showing the large structural impact of TEs on yeast genome and opening the door to further investigations to understand their functional and evolutionary consequences.

The structure and retrotransposition mechanism of LTR-retrotransposons in the asexual yeast Candida albicans

Retrotransposons constitute a major part of the genome in a number of eukaryotes. Long-terminal repeat (LTR) retrotransposons are one type of the retrotransposons. Candida albicans have 34 distinct LTR-retrotransposon families. They respectively belong to the Ty1/copia and Ty3/gypsy groups which have been extensively studied in the model yeast Saccharomyces cerevisiae. LTR-retrotransposons carry two LTRs flanking a long internal protein-coding domain, open reading frames. LTR-retrotransposons use RNA as intermediate to synthesize double-stranded DNA copies. In this article, we describe the structure feature, retrotransposition mechanism and the influence on organism diversity of LTR retrotransposons in C. albicans. We also discuss the relationship between pathogenicity and LTR retrotransposons in C. albicans.

Structural features and mechanism of translocation of non-LTR retrotransposons in Candida albicans

A number of abundant mobile genetic elements called retrotransposons reverse transcribe RNA to generate DNA for insertion into eukaryotic genomes. Non-long-terminal repeat (non-LTR) retrotransposons represent a major class of retrotransposons, and transposons that move by target-primed reverse transcription lack LTRs characteristic of retroviruses and retroviral-like transposons. Yeast model systems in Candida albicans and Saccharomyces cerevisiae have been developed for the study of non-LTR retrotransposons. Non-LTR retrotransposons are divided into LINEs (long interspersed nuclear elements), SINEs (short interspersed nuclear elements), and SVA (SINE, VNTR, and Alu). LINE-1 elements have been described in fungi, and several families called Zorro elements have been detected from C. albicans. They are all members of L1 clades. Through a mechanism named target-primed reverse transcription (TPRT), LINEs translocate the new copy into the target site to initiate DNA synthesis primed by the 3′ OH of the broken strand. In this article, we describe some advances in the research on structural features and origin of non-LTR retrotransposons in C. albicans, and discuss mechanisms underlying their reverse transcription and integration of the donor copy into the target site.

Molecular epidemiology, phylogeny and evolution of Candida albicans

A small number of Candida species form part of the normal microbial flora of mucosal surfaces in humans and may give rise to opportunistic infections when host defences are impaired. Candida albicans is by far the most prevalent commensal and pathogenic Candida species. Several different molecular typing approaches including multilocus sequence typing, multilocus microsatellite typing and DNA fingerprinting using C. albicans-specific repetitive sequence-containing DNA probes have yielded a wealth of information regarding the epidemiology and population structure of this species. Such studies revealed that the C. albicans population structure consists of multiple major and minor clades, some of which exhibit geographical or phenotypic enrichment and that C. albicans reproduction is predominantly clonal. Despite this, losses of heterozygosity by recombination, the existence of a parasexual cycle, toleration of a wide range of aneuploidies and the recent description of viable haploid strains have all demonstrated the extensive plasticity of the C. albicans genome. Recombination and gross chromosomal rearrangements are more common under stressful environmental conditions, and have played a significant role in the evolution of this opportunistic pathogen. Surprisingly, Candida dubliniensis, the closest relative of C. albicans exhibits more karyotype variability than C. albicans, but is significantly less adaptable to unfavourable environments. This disparity most likely reflects the evolutionary processes that occurred during or soon after the divergence of both species from their common ancestor. Whilst C. dubliniensis underwent significant gene loss and pseudogenisation, C. albicans expanded gene families considered to be important in virulence. It is likely that technological developments in whole genome sequencing and data analysis in coming years will facilitate its routine use for population structure, epidemiological investigations, and phylogenetic analyses of Candida species. These are likely to reveal more minor C. albicans clades and to enhance our understanding of the population biology of this versatile organism.

Functional and structural divergence of an unusual LTR retrotransposon family in plants

Retrotransposons with long terminal repeats (LTRs) more than 3 kb are not frequent in most eukaryotic genomes. Rice LTR retrotransposon, Retrosat2, has LTRs greater than 3.2 kb and two open reading frames (ORF): ORF1 encodes enzymes for retrotransposition whereas no function can be assigned to ORF0 as it is not found in any other organism. A variety of experimental and in silico approaches were used to determine the origin of Retrosat2 and putative function of ORF0. Our data show that not only is Retrosat2 highly abundant in the Oryza genus, it may yet be active in rice. Homologs of Retrosat2 were identified in maize, sorghum, Arabidopsis and other plant genomes suggesting that the Retrosat2 family is of ancient origin. Several putatively cis-acting elements, some multicopy, that regulate retrotransposon replication or responsiveness to environmental factors were found in the LTRs of Retrosat2. Unlike the ORF1, the ORF0 sequences from Retrosat2 and homologs are divergent at the sequence level, 3D-structures and predicted biological functions. In contrast to other retrotransposon families, Retrosat2 and its homologs are dispersed throughout genomes and not concentrated in the specific chromosomal regions, such as centromeres. The genomic distribution of Retrosat2 homologs varies across species which likely reflects the differing evolutionary trajectories of this retrotransposon family across diverse species.

DNA methylation regulates phenotype-dependent transcriptional activity in Candida albicans

DNA methylation is a common epigenetic signaling mechanism associated with silencing of repeated DNA and transcriptional regulation in eukaryotes. Here we report that DNA methylation in the human fungal pathogen Candida albicans is primarily localized within structural genes and modulates transcriptional activity. Major repeat sequences and multigene families are largely free of DNA methylation. Among the genes subject to DNA methylation are those associated with dimorphic transition between yeast and hyphal forms, switching between white and opaque cells, and iron metabolism. Transcriptionally repressed methylated loci showed increased frequency of C-to-T transitions during asexual growth, an evolutionarily stable pattern of repression associated mutation that could bring about genetic alterations under changing environmental or host conditions. Dynamic differential DNA methylation of structural genes may be one factor contributing to morphological plasticity that is cued by nutrition and host interaction.

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.

Fungal sex and pathogenesis

Human fungal pathogens are associated with diseases ranging from dandruff and skin colonization to invasive bloodstream infections. The major human pathogens belong to the Candida, Aspergillus, and Cryptococcus clades, and infections have high and increasing morbidity and mortality. Many human fungal pathogens were originally assumed to be asexual. However, recent advances in genome sequencing, which revealed that many species have retained the genes required for the sexual machinery, have dramatically influenced our understanding of the biology of these organisms. Predictions of a rare or cryptic sexual cycle have been supported experimentally for some species. Here, I examine the evidence that human pathogens reproduce sexually. The evolution of the mating-type locus in ascomycetes (including Candida and Aspergillus species) and basidiomycetes (Malassezia and Cryptococcus) is discussed. I provide an overview of how sex is suppressed in different species and discuss the potential associations with pathogenesis.

Exploring repetitive DNA landscapes using REPCLASS, a tool that automates the classification of transposable elements in eukaryotic genomes

Eukaryotic genomes contain large amount of repetitive DNA, most of which is derived from transposable elements (TEs). Progress has been made to develop computational tools for ab initio identification of repeat families, but there is an urgent need to develop tools to automate the annotation of TEs in genome sequences. Here we introduce REPCLASS, a tool that automates the classification of TE sequences. Using control repeat libraries, we show that the program can classify accurately virtually any known TE types. Combining REPCLASS to ab initio repeat finding in the genomes of Caenorhabditis elegans and Drosophila melanogaster allowed us to recover the contrasting TE landscape characteristic of these species. Unexpectedly, REPCLASS also uncovered several novel TE families in both genomes, augmenting the TE repertoire of these model species. When applied to the genomes of distant Caenorhabditis and Drosophila species, the approach revealed a remarkable conservation of TE composition profile within each genus, despite substantial interspecific covariations in genome size and in the number of TEs and TE families. Lastly, we applied REPCLASS to analyze 10 fungal genomes from a wide taxonomic range, most of which have not been analyzed for TE content previously. The results showed that TE diversity varies widely across the fungi “kingdom” and appears to positively correlate with genome size, in particular for DNA transposons. Together, these data validate REPCLASS as a powerful tool to explore the repetitive DNA landscapes of eukaryotes and to shed light onto the evolutionary forces shaping TE diversity and genome architecture.

Loss of heterozygosity in commensal isolates of the asexual diploid yeast Candida albicans

Candida albicans is a commensal and the most frequent fungal pathogen of humans. One mechanism of genetic variation in this diploid asexual yeast involves loss of heterozygosity (LOH). LOH events occur upon infection and contribute to the acquisition of antifungal resistance in patients. In contrast, little is known about the nature and extent of LOH events during commensalism. Using a combination of single nucleotide polymorphism typing, positional transcript profiling and karyotyping, we have characterized related C. albicans commensal isolates that differ by LOH events. Most of these LOH events encompassed the entirety of the chromosome or a large region extending to the telomere, suggesting chromosome loss or mitotic recombination/break-induced replication events, respectively. They were frequently accompanied by karyotype alterations such as chromosome length polymorphism and copy number variations at other chromosomes. These results demonstrate the high plasticity of the C. albicans genome during commensalism.

Phylogenetic profiles reveal evolutionary relationships within the "twilight zone" of sequence similarity

Inferring evolutionary relationships among highly divergent protein sequences is a daunting task. In particular, when pairwise sequence alignments between protein sequences fall <25% identity, the phylogenetic relationships among sequences cannot be estimated with statistical certainty. Here, we show that phylogenetic profiles generated with the Gestalt Domain Detection Algorithm-Basic Local Alignment Tool (GDDA-BLAST) are capable of deriving, ab initio, phylogenetic relationships for highly divergent proteins in a quantifiable and robust manner. Notably, the results from our computational case study of the highly divergent family of retroelements accord with previous estimates of their evolutionary relationships. Taken together, these data demonstrate that GDDA-BLAST provides an independent and powerful measure of evolutionary relationships that does not rely on potentially subjective sequence alignment. We demonstrate that evolutionary relationships can be measured with phylogenetic profiles, and therefore propose that these measurements can provide key insights into relationships among distantly related and/or rapidly evolving proteins.

The role of alternative genetic codes in viral evolution and emergence

Although the 'universal' genetic code is widespread among life-forms, a number of diverse lineages have evolved unique codon reassignments. The proteomes of these organisms and organelles must, by necessity, use the same codon assignments. Likewise, for an exogenous genetic element, such as an infecting viral genome, to be accurately and completely expressed with the host's translation system, it must employ the same genetic code. This raises a number of intriguing questions regarding the origin and evolution of viruses. In particular, it is extremely unlikely that viruses of hosts utilizing the universal genetic code would emerge, via cross-species transmission, in hosts utilizing alternative codes, and vice versa. Consequently, more parsimonious scenarios for the origins of such viruses include the prolonged co-evolution of viruses with cellular life, or the escape of genetic material from host genomes. Further, we raise the possibility that emerging viruses provide the selection pressure favoring the use of alternative codes in potential hosts, such that the evolution of a variant genetic code acts as a unique and powerful antiviral strategy. As such, in the face of new emerging viruses, hosts with codon reassignments would have a significant selective advantage compared to hosts utilizing the universal code.

Assembly of the Candida albicans genome into sixteen supercontigs aligned on the eight chromosomes

For Assembly 20 of the Candida albicans genome, the sequence of each of the eight chromosomes was determined, revealing new insights into gene family creation and dispersion, subtelomere organization, and chromosome evolution.

Background

The 10.9× genomic sequence of Candida albicans, the most important human fungal pathogen, was published in 2004. Assembly 19 consisted of 412 supercontigs, of which 266 were a haploid set, since this fungus is diploid and contains an extensive degree of heterozygosity but lacks a complete sexual cycle. However, sequences of specific chromosomes were not determined.

Results

Supercontigs from Assembly 19 (183, representing 98.4% of the sequence) were assigned to individual chromosomes purified by pulse-field gel electrophoresis and hybridized to DNA microarrays. Nine Assembly 19 supercontigs were found to contain markers from two different chromosomes. Assembly 21 contains the sequence of each of the eight chromosomes and was determined using a synteny analysis with preliminary versions of the Candida dubliniensis genome assembly, bioinformatics, a sequence tagged site (STS) map of overlapping fosmid clones, and an optical map. The orientation and order of the contigs on each chromosome, repeat regions too large to be covered by a sequence run, such as the ribosomal DNA cluster and the major repeat sequence, and telomere placement were determined using the STS map. Sequence gaps were closed by PCR and sequencing of the products. The overall assembly was compared to an optical map; this identified some misassembled contigs and gave a size estimate for each chromosome.

Conclusion

Assembly 21 reveals an ancient chromosome fusion, a number of small internal duplications followed by inversions, and a subtelomeric arrangement, including a new gene family, the TLO genes. Correlations of position with relatedness of gene families imply a novel method of dispersion. The sequence of the individual chromosomes of C. albicans raises interesting biological questions about gene family creation and dispersion, subtelomere organization, and chromosome evolution.

Centromere size and position in Candida albicans are evolutionarily conserved independent of DNA sequence heterogeneity

The centromere regions (CEN) of all eight chromosomes in Candida albicans have been characterized in terms of nucleotide sequence and size. The boundaries of each of the eight CEN DNA regions were mapped by chromatin immunoprecipitation-PCR using polyclonal rabbit antibodies generated against C. albicans centromere-specific protein CaCse4p (CENP-A homolog). A single 3-4.5 kb unique DNA sequence on each chromosome was found to be bound to CaCse4p. Sequence analysis revealed that the eight CEN regions in C. albicans lack any conserved DNA sequence motifs common to the group; all are quite different in overall DNA sequence. In contrast to centromeres in many organisms, the C. albicans centromeres are generally free of repeated DNA elements and transposons. However, a few small inverted repeats and long terminal repeats do occur in the centromeric and pericentric regions on a few chromosomes. We also characterized the CEN DNAs in four groups of phylogenetically divergent C. albicans strains, estimated to be separated from each other by 1-3 million years. The same eight different and unique 3-4.5 kb DNA sequences are utilized as centromeres in all of these strains. The chromosomal locations and the sizes of CEN DNAs have remained conserved, in agreement with the idea that CEN function in C. albicans is templated by heritable epigenetic mechanisms.

Structure and transcriptional impact of divergent repetitive elements inserted within Phanerochaete chrysosporium strain RP-78 genes

We describe the structure, organization, and transcriptional impact of repetitive elements within the lignin-degrading basidiomycete, Phanerochaete chrysosporium. Searches of the P. chrysosporium genome revealed five copies of pce1, a ~1,750-nt non-autonomous, class II element. Alleles encoding a putative glucosyltransferase and a cytochrome P450 harbor pce insertions and produce incomplete transcripts. Class I elements included pcret1, an intact 8.14-kb gypsy-like retrotransposon inserted within a member of the multicopper oxidase gene family. Additionally, we describe a complex insertion of nested transposons within another putative cytochrome P450 gene. The disrupted allele lies within a cluster of >14 genes, all of which encode family 64 cytochrome P450s. Components of the insertion include a disjoint copia-like element, pcret3, the pol domain of a second retroelement, pcret2, and a duplication of an extended ORF of unknown function. As in the case of the pce elements, pcret1 and pcret2/3 insertions are confined to single alleles, transcripts of which are truncated. The corresponding wild-type alleles are apparently unaffected. In aggregate, P. chrysosporium harbors a complex array of repetitive elements, at least five of which directly influence expression of genes within families of structurally related sequences.

The long hard road to a completed Candida albicans genome

After almost a decade of work, the sequencing, assembly, and annotation of the genome of the fungal pathogen Candida albicans is finally close at hand. This review covers the early history of the C. albicans genome project, from the release of early assemblies that provided the impetus for an explosion in functional genomics research, to a community-based annotation and a preview of the work that was necessary for the production of a final genome assembly.

[Structural and functional analyses of MRS (major repeated sequnece) in Candida albicans. Application for genotyping and the eternal way to the complete genome sequence]

There are several different types of repeated sequences in the genome of Candida albicans, including the MRS (Major repeated sequence). In 2004, the whole genome sequence of C. albicans was published. Assembly of the sequences to chromosomal length contigs was not achieved, mainly due to interruption of the sequences by MRS. However, MRS including Ca3, 27A and RPS have been playing important roles in a number of epidemiological studies and basic biological investigations into C. albicans chromosome loss events and associated phenotypic changes. Here we summarize structural analyses from subrepeat sequences to the chromosome level, and functional analyses of MRS.

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.

Genome-wide transcription profiling of the early phase of biofilm formation by Candida albicans

The ability to adhere to surfaces and develop as a multicellular community is an adaptation used by most microorganisms to survive in changing environments. Biofilm formation proceeds through distinct developmental phases and impacts not only medicine but also industry and evolution. In organisms such as the opportunistic pathogen Candida albicans, the ability to grow as biofilms is also an important mechanism for persistence, facilitating its growth on different tissues and a broad range of abiotic surfaces used in medical devices. The early stage of C. albicans biofilm is characterized by the adhesion of single cells to the substratum, followed by the formation of an intricate network of hyphae and the beginning of a dense structure. Changes in the transcriptome begin within 30 min of contact with the substrate and include expression of genes related to sulfur metabolism, in particular MET3, and the equivalent gene homologues of the Ribi regulon in Saccharomyces cerevisiae. Some of these changes are initiated early and maintained throughout the process; others are restricted to the earliest stages of biofilm formation. We identify here a potential alternative pathway for cysteine metabolism and the biofilm-associated expression of genes involved in glutathione production in C. albicans.

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.

APE-type non-LTR retrotransposons: determinants involved in target site recognition

Non-long terminal repeat (Non-LTR) retrotransposons represent a diverse and widely distributed group of transposable elements and an almost ubiquitous component of eukaryotic genomes that has a major impact on evolution. Their copy number can range from a few to several million and they often make up a significant fraction of the genomes. The members of the dominating subtype of non-LTR retrotransposons code for an endonuclease with homology to apurinic/apyrimidinic endonucleases (APE), and are thus termed APE-type non-LTR retrotransposons. In the last decade both the number of identified non-LTR retrotransposons and our knowledge of biology and evolution of APE-type non-LTR retrotransposons has increased tremendously.

Genome evolution mediated by Ty elements in Saccharomyces

How mobile genetic elements molded eukaryotic genomes is a key evolutionary question that gained wider popularity when mobile DNA sequences were shown to comprise about half of the human genome. Although Saccharomyces cerevisiae does not suffer such "genome obesity", five families of LTR-retrotransposons, Ty1, Ty2, Ty3, Ty4, and Ty5 elements, comprise about 3% of its genome. The availability of complete genome sequences from several Saccharomyces species, including members of the closely related sensu stricto group, present new opportunities for analyzing molecular mechanisms for chromosome evolution, speciation, and reproductive isolation. In this review I present key experiments from both the pre- and current genomic sequencing eras suggesting how Ty elements mediate genome evolution.

Sequence finishing and gene mapping for Candida albicans chromosome 7 and syntenic analysis against the Saccharomyces cerevisiae genome

The size of the genome in the opportunistic fungus Candida albicans is 15.6 Mb. Whole-genome shotgun sequencing was carried out at Stanford University where the sequences were assembled into 412 contigs. C. albicans is a diploid basically, and analysis of the sequence is complicated due to repeated sequences and to sequence polymorphism between homologous chromosomes. Chromosome 7 is 1 Mb in size and the best characterized of the 8 chromosomes in C. albicans. We assigned 16 of the contigs, ranging in length from 7309 to 267,590 bp, to chromosome 7 and determined sequences of 16 regions. These regions included four gaps, a misassembled sequence, and two major repeat sequences (MRS) of >16 kb. The length of the continuous sequence attained was 949,626 bp and provided complete coverage of chromosome 7 except for telomeric regions. Sequence analysis was carried out and predicted 404 genes, 11 of which included at least one intron. A 7-kb indel, which might be caused by a retrotransposon, was identified as the largest difference between the homologous chromosomes. Synteny analysis revealed that the degree of synteny between C. albicans and Saccharomyces cerevisiae is too weak to use for completion of the genomic sequence in C. albicans.

The genome of the social amoeba Dictyostelium discoideum

The social amoebae are exceptional in their ability to alternate between unicellular and multicellular forms. Here we describe the genome of the best-studied member of this group, Dictyostelium discoideum. The gene-dense chromosomes of this organism encode approximately 12,500 predicted proteins, a high proportion of which have long, repetitive amino acid tracts. There are many genes for polyketide synthases and ABC transporters, suggesting an extensive secondary metabolism for producing and exporting small molecules. The genome is rich in complex repeats, one class of which is clustered and may serve as centromeres. Partial copies of the extrachromosomal ribosomal DNA (rDNA) element are found at the ends of each chromosome, suggesting a novel telomere structure and the use of a common mechanism to maintain both the rDNA and chromosomal termini. A proteome-based phylogeny shows that the amoebozoa diverged from the animal-fungal lineage after the plant-animal split, but Dictyostelium seems to have retained more of the diversity of the ancestral genome than have plants, animals or fungi.

The functional basis of mycophenolic acid resistance in Candida albicans IMP dehydrogenase

Candida albicans is an important fungal pathogen of immunocompromised patients. In cell culture, C. albicans is sensitive to mycophenolic acid (MPA) and mizoribine, both natural product inhibitors of IMP dehydrogenase (IMPDH). These drugs have opposing interactions with the enzyme. MPA prevents formation of the closed enzyme conformation by binding to the same site as a mobile flap. In contrast, mizoribine monophosphate, the active metabolite of mizoribine, induces the closed conformation. Here, we report the characterization of IMPDH from wild-type and MPA-resistant strains of C. albicans. The wild-type enzyme displays significant differences from human IMPDHs, suggesting that selective inhibitors that could be novel antifungal agents may be developed. IMPDH from the MPA-resistant strain contains a single substitution (A251T) that is far from the MPA-binding site. The A251T variant was 4-fold less sensitive to MPA as expected. This substitution did not affect the k(cat) value, but did decrease the K(m) values for both substrates, so the mutant enzyme is more catalytically efficient as measured by the value of k(cat)/K(m). These simple criteria suggest that the A251T variant would be the evolutionarily superior enzyme. However, the A251T substitution caused the enzyme to be 40-fold more sensitive to mizoribine monophosphate. This result suggests that A251T stabilizes the closed conformation, and this hypothesis is supported by further inhibitor analysis. Likewise, the MPA-resistant strain was more sensitive to mizoribine in cell culture. These observations illustrate the evolutionary challenge posed by the gauntlet of chemical warfare at the microbial level.

The diploid genome sequence of Candida albicans

We present the diploid genome sequence of the fungal pathogen Candida albicans. Because C. albicans has no known haploid or homozygous form, sequencing was performed as a whole-genome shotgun of the heterozygous diploid genome in strain SC5314, a clinical isolate that is the parent of strains widely used for molecular analysis. We developed computational methods to assemble a diploid genome sequence in good agreement with available physical mapping data. We provide a whole-genome description of heterozygosity in the organism. Comparative genomic analyses provide important clues about the evolution of the species and its mechanisms of pathogenesis.

Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78

White rot fungi efficiently degrade lignin, a complex aromatic polymer in wood that is among the most abundant natural materials on earth. These fungi use extracellular oxidative enzymes that are also able to transform related aromatic compounds found in explosive contaminants, pesticides and toxic waste. We have sequenced the 30-million base-pair genome of Phanerochaete chrysosporium strain RP78 using a whole genome shotgun approach. The P. chrysosporium genome reveals an impressive array of genes encoding secreted oxidases, peroxidases and hydrolytic enzymes that cooperate in wood decay. Analysis of the genome data will enhance our understanding of lignocellulose degradation, a pivotal process in the global carbon cycle, and provide a framework for further development of bioprocesses for biomass utilization, organopollutant degradation and fiber bleaching. This genome provides a high quality draft sequence of a basidiomycete, a major fungal phylum that includes important plant and animal pathogens.

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.

Retrotransposons and their recognition of pol II promoters: a comprehensive survey of the transposable elements from the complete genome sequence of Schizosaccharomyces pombe

The complete DNA sequence of the genome of Schizosaccharomyces pombe provides the opportunity to investigate the entire complement of transposable elements (TEs), their association with specific sequences, their chromosomal distribution, and their evolution. Using homology-based sequence identification, we found that the sequenced strain of S. pombe contained only one family of full-length transposons. This family, Tf2, consisted of 13 full-length copies of a long terminal repeat (LTR) retrotransposon. We found that LTR-LTR recombination of previously existing transposons had resulted in extensive populations of solo LTRs. These included 35 solo LTRs of Tf2, as well as 139 solo LTRs from other Tf families. Phylogenetic analysis of solo Tf LTRs reveals that Tf1 and Tf2 were the most recently active elements within the genome. The solo LTRs also served as footprints for previous insertion events by the Tf retrotransposons. Analysis of 186 genomic insertion events revealed a close association with RNA polymerase II promoters. These insertions clustered in the promoter-proximal regions of genes, upstream of protein coding regions by 100 to 400 nucleotides. The association of Tf insertions with pol II promoters was very similar to the preference previously observed for Tf1 integration. We found that the recently active Tf elements were absent from centromeres and pericentromeric regions of the genome containing tandem tRNA gene clusters. In addition, our analysis revealed that chromosome III has twice the density of insertion events compared to the other two chromosomes. Finally we describe a novel repetitive sequence, wtf, which was also preferentially located on chromosome III, and was often located near solo LTRs of Tf elements.

Ty3/gypsy-like retrotransposons in Candida albicans and Candida dubliniensis: Tca3 and Tcd3

Ty3/gypsy retrotransposons are a widespread group of eukaryote mobile genetic elements. They are similar in structure to, and may be ancestors of, the vertebrate retroviruses. Here we describe the first Ty3/gypsy retrotransposons from the pathogenic yeasts Candida albicans and Candida dubliniensis, which we refer to as Tca3 and Tcd3, respectively. Tca3 was first identified in a variety of strains as an element lacking a large part of its coding region. Comparative analyses between C. albicans and C. dubliniensis allowed us to identify the closely related full-length Tcd3 element, and, subsequently, the full-length Tca3 elements. The full-length versions of Tca3 and Tcd3 are broadly similar in structure to other Ty3/gypsy elements, but have several features of special interest, e.g. both elements appear to have a novel mechanism for priming minus-strand DNA synthesis, probably involving conserved secondary structures adjacent to the 5' LTRs. Also, while closely related to each other, the two elements appear to be fairly distantly related to other known Ty3/gypsy-like elements. Finally, the occurrence of the internally deleted forms of Tca3 in many strains raises interesting questions concerning the evolution of these transposable elements in Candida and the evolution of Candida itself. The sequences reported in this paper have been assigned GenBank Accession Nos AF499463, AF499464 and AF510498.

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.

Three retrotransposon families in the genome of Giardia lamblia: two telomeric, one dead

Transposable elements inhabiting eukaryotic genomes are generally regarded either as selfish DNA, which is selectively neutral to the host organism, or as parasitic DNA, deleterious to the host. Thus far, the only agreed-upon example of beneficial eukaryotic transposons is provided by Drosophila telomere-associated retrotransposons, which transpose directly to the chromosome ends and thereby protect them from degradation. This article reports the transposon content of the genome of the protozoan Giardia lamblia, one of the earliest-branching eukaryotes. A total of three non-long terminal repeat retrotransposon families have been identified, two of which are located at the ends of chromosomes, and the third one contains exclusively dead copies with multiple internal deletions, nucleotide substitutions, and frame shifts. No other reverse transcriptase- or transposase-related sequences were found. Thus, the entire genome of this protozoan, which is not known to reproduce sexually, contains only retrotransposons that are either confined to telomeric regions and possibly beneficial, or inactivated and completely nonfunctional.

Evolutionary history of Cer elements and their impact on the C. elegans genome

We report the results of sequence analysis and chromosomal distribution of all distinguishable long terminal repeat (LTR) retrotransposons (Cer elements) in the Caenorhabditis elegans genome. Included in this analysis are all readily recognizable full-length and fragmented elements, as well as solo LTRs. Our results indicate that there are 19 families of Cer elements, some of which display significant subfamily structure. Cer elements can be clustered based on their tRNA primer binding sites (PBSs). These clusters are in concordance with our reverse transcriptase- and LTR-based phylogenies. Although we find that most Cer elements are located in the gene depauperate chromosome ends, some elements are located in or near putative genes and may contribute to gene structure and function. The results of RT-PCR analyses are consistent with this prediction.