An active retrotransposon in Candida albicans

Candida albicans isolates contain frequent heterozygous structural variants and transposable elements within genes and centromeres

The human fungal pathogen Candida albicans poses a significant burden on global health, causing high rates of mortality and antifungal drug resistance. C. albicans is a heterozygous diploid organism that reproduces asexually. Structural variants (SVs) are an important source of genomic rearrangement, particularly in species that lack sexual recombination. To comprehensively investigate SVs across clinical isolates of C. albicans, we conducted long-read sequencing and genome-wide SV analysis in three distantly related clinical isolates. Our work includes a new, comprehensive analysis of transposable element (TE) composition, location, and diversity. SVs and TEs are frequently close to coding sequences and many SVs are heterozygous, suggesting that SVs might impact gene and allele-specific expression. Most SVs are uniquely present in only one clinical isolate, indicating that SVs represent a significant source of intraspecies genetic variation. We identify multiple, distinct SVs at the centromeres of Chromosome 4 and Chromosome 5, including inversions and transposon polymorphisms. These two chromosomes are often aneuploid in drug-resistant clinical isolates and can form isochromosome structures with breakpoints near the centromere. Further screening of 100 clinical isolates confirms the widespread presence of centromeric SVs in C. albicans, often appearing in a heterozygous state, indicating that SVs are contributing to centromere evolution in C. albicans Together, these findings highlight that SVs and TEs are common across diverse clinical isolates of C. albicans and that the centromeres of this organism are important sites of genome rearrangement.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

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.

Comparative metatranscriptome analysis revealed broad response of microbial communities in two soil types, agriculture versus organic soil

BACKGROUND

Studying expression of genes by direct sequencing and analysis of metatranscriptomes at a particular time and space can disclose structural and functional insights about microbial communities. The present study reports comparative analysis of metatranscriptome from two distinct soil ecosystems referred as M1 (agriculture soil) and O1 (organic soil).

RESULTS

Analysis of sequencing reads revealed Proteobacteria as major dominant phyla in both soil types. The order of the top 3 abundant phyla in M1 sample was Proteobacteria > Ascomycota > Firmicutes, whereas in sample O1, the order was Proteobacteria > Cyanobacteria > Actinobacteria. Analysis of differentially expressed genes demonstrated high expression of transcripts related to copper-binding proteins, proteins involved in electron carrier activity, DNA integration, endonuclease activity, MFS transportation, and other uncharacterized proteins in M1 compared to O1. Of the particular interests, several transcripts related to nitrification, ammonification, stress response, and alternate carbon fixation pathways were highly expressed in M1. In-depth analysis of the sequencing data revealed that transcripts of archaeal origin had high expression in M1 compared to O1 indicating the active role of Archaea in metal- and pesticide-contaminated environment. In addition, transcripts encoding 4-hydroxyphenylpyruvate dioxygenase, glyoxalase/bleomycin resistance protein/dioxygenase, metapyrocatechase, and ring hydroxylating dioxygenases of aromatic hydrocarbon degradation pathways had high expression in M1. Altogether, this study provided important insights about the transcripts and pathways upregulating in the presence of pesticides and herbicides.

CONCLUSION

Altogether, this study claims a high expression of microbial transcripts in two ecosystems with a wide range of functions. It further provided clue about several molecular markers which could be a strong indicator of metal and pesticide contamination in soils. Interestingly, our study revealed that Archaea are playing a significant role in nitrification process as compared to bacteria in metal- and pesticide-contaminated soil. In particular, high expression of transcripts related to aromatic hydrocarbon degradation in M1 soil indicates their important role in biodegradation of pollutants, and therefore, further investigation is needed.

Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies

Transposable Elements (TEs) are mobile DNA sequences that make up significant fractions of amniote genomes. However, they are difficult to detect and annotate ab initio because of their variable features, lengths and clade-specific variants. We have addressed this problem by refining and developing a Comprehensive ab initio Repeat Pipeline (CARP) to identify and cluster TEs and other repetitive sequences in genome assemblies. The pipeline begins with a pairwise alignment using krishna, a custom aligner. Single linkage clustering is then carried out to produce families of repetitive elements. Consensus sequences are then filtered for protein coding genes and then annotated using Repbase and a custom library of retrovirus and reverse transcriptase sequences. This process yields three types of family: fully annotated, partially annotated and unannotated. Fully annotated families reflect recently diverged/young known TEs present in Repbase. The remaining two types of families contain a mixture of novel TEs and segmental duplications. These can be resolved by aligning these consensus sequences back to the genome to assess copy number vs. length distribution. Our pipeline has three significant advantages compared to other methods for ab initio repeat identification: 1) we generate not only consensus sequences, but keep the genomic intervals for the original aligned sequences, allowing straightforward analysis of evolutionary dynamics, 2) consensus sequences represent low-divergence, recently/currently active TE families, 3) segmental duplications are annotated as a useful by-product. We have compared our ab initio repeat annotations for 7 genome assemblies to other methods and demonstrate that CARP compares favourably with RepeatModeler, the most widely used repeat annotation package.

Sequence Assembly of Yarrowia lipolytica Strain W29/CLIB89 Shows Transposable Element Diversity

Yarrowia lipolytica, an oleaginous yeast, is capable of accumulating significant cellular mass in lipid making it an important source of biosustainable hydrocarbon-based chemicals. In spite of a similar number of protein-coding genes to that in other Hemiascomycetes, the Y. lipolytica genome is almost double that of model yeasts. Despite its economic importance and several distinct strains in common use, an independent genome assembly exists for only one strain. We report here a de novo annotated assembly of the chromosomal genome of an industrially-relevant strain, W29/CLIB89, determined by hybrid next-generation sequencing. For the first time, each Y. lipolytica chromosome is represented by a single contig. The telomeric rDNA repeats were localized by Irys long-range genome mapping and one complete copy of the rDNA sequence is reported. Two large structural variants and retroelement differences with reference strain CLIB122 including a full-length, novel Ty3/Gypsy long terminal repeat (LTR) retrotransposon and multiple LTR-like sequences are described. Strikingly, several of these are adjacent to RNA polymerase III-transcribed genes, which are almost double in number in Y. lipolytica compared to other Hemiascomycetes. In addition to previously-reported dimeric RNA polymerase III-transcribed genes, tRNA pseudogenes were identified. Multiple full-length and truncated LINE elements are also present. Therefore, although identified transposons do not constitute a significant fraction of the Y. lipolytica genome, they could have played an active role in its evolution. Differences between the sequence of this strain and of the existing reference strain underscore the utility of an additional independent genome assembly for this economically important organism.

Genetic and phenotypic intra-species variation in Candida albicans

Candida albicans is a commensal fungus of the human gastrointestinal tract and a prevalent opportunistic pathogen. To examine diversity within this species, extensive genomic and phenotypic analyses were performed on 21 clinical C. albicans isolates. Genomic variation was evident in the form of polymorphisms, copy number variations, chromosomal inversions, subtelomeric hypervariation, loss of heterozygosity (LOH), and whole or partial chromosome aneuploidies. All 21 strains were diploid, although karyotypic changes were present in eight of the 21 isolates, with multiple strains being trisomic for Chromosome 4 or Chromosome 7. Aneuploid strains exhibited a general fitness defect relative to euploid strains when grown under replete conditions. All strains were also heterozygous, yet multiple, distinct LOH tracts were present in each isolate. Higher overall levels of genome heterozygosity correlated with faster growth rates, consistent with increased overall fitness. Genes with the highest rates of amino acid substitutions included many cell wall proteins, implicating fast evolving changes in cell adhesion and host interactions. One clinical isolate, P94015, presented several striking properties including a novel cellular phenotype, an inability to filament, drug resistance, and decreased virulence. Several of these properties were shown to be due to a homozygous nonsense mutation in the EFG1 gene. Furthermore, loss of EFG1 function resulted in increased fitness of P94015 in a commensal model of infection. Our analysis therefore reveals intra-species genetic and phenotypic differences in C. albicans and delineates a natural mutation that alters the balance between commensalism and pathogenicity.

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.

A core filamentation response network in Candida albicans is restricted to eight genes

Although morphological plasticity is a central virulence trait of Candida albicans, the number of filament-associated genes and the interplay of mechanisms regulating their expression remain unknown. By correlation-based network modeling of the transcriptional response to different defined external stimuli for morphogenesis we identified a set of eight genes with highly correlated expression patterns, forming a core filamentation response. This group of genes included ALS3, ECE1, HGT2, HWP1, IHD1 and RBT1 which are known or supposed to encode for cell- wall associated proteins as well as the Rac1 guanine nucleotide exchange factor encoding gene DCK1 and the unknown function open reading frame orf19.2457. The validity of network modeling was confirmed using a dataset of advanced complexity that describes the transcriptional response of C. albicans during epithelial invasion as well as comparing our results with other previously published transcriptome studies. Although the set of core filamentation response genes was quite small, several transcriptional regulators are involved in the control of their expression, depending on the environmental condition.

High-frequency genetic contents variations in clinical Candida albicans isolates

Genome plasticity is a hallmark of Candida albicans and is believed to be an adaptation strategy. But the extent of such genomic variability is not well investigated. In this study, genetic contents of clinical C. albicans isolates were investigated at whole-genome level with array-based comparative genomic hybridization (array CGH) technology. It was revealed that C. albicans possessed variations of genetic contents, as well as aneuploidy. The variable genes were scattered across the chromosomes, as well clustered in particular regions, including sub-telomeric regions, retrotransposon-insertion sites and a variable region on chromosome 6.

Transposable elements in yeasts

With the development of new sequencing technologies in the past decade, yeast genomes have been extensively sequenced and their structures investigated. Transposable elements (TEs) are ubiquitous in eukaryotes and constitute a limited part of yeast genomes. However, due to their ability to move in genomes and generate dispersed repeated sequences, they contribute to modeling yeast genomes and thereby induce plasticity. This review assesses the TE contents of yeast genomes investigated so far. Their diversity and abundance at the inter- and intraspecific levels are presented, and their effects on gene expression and genome stability is considered. Recent results concerning TE-host interactions are also analyzed.

Phenotypic screening, transcriptional profiling, and comparative genomic analysis of an invasive and non-invasive strain of Candida albicans

Background

Invasion of host tissue by the human fungal pathogen Candida albicans is an important step during the development of candidosis. However, not all C. albicans strains possess the same invasive and virulence properties. For example, the two clinical isolates SC5314 and ATCC10231 differ in their ability to invade host tissue and cause experimental infections. Strain SC5314 is invasive whereas strain ATCC10231 is non-invasive and strongly attenuated in virulence compared to SC5314. In this study we compare the in vitro phenotypic, transcriptional and genomic profiles of these two widely used laboratory strains in order to determine the principal biological and genetic properties responsible for their differential virulence.

Results

In all media tested, the two strains showed the same metabolic flexibility, stress resistance, adhesion properties and hydrolytic enzyme secretion in vitro. However, differences were observed in response to cell-surface disturbing agents and alkaline pH. Furthermore, reduced hyphal formation in strain ATCC10231 under certain conditions correlated with reduced invasive properties in an in vitro invasion assay and a reduced ability to invade epithelial tissue. Despite these diverse phenotypic properties, no substantial genomic differences were detected by comparative genome hybridisation within the open reading frames. However, in vitro transcriptional profiling displayed major differences in the gene expression of these two strains, even under normal in vitro growth conditions.

Conclusion

Our data suggest that the reason for differential virulence of C. albicans strains is not due to the absence of specific genes, but rather due to differences in the expression, function or activity of common genes.

Control of gag-pol gene expression in the Candida albicans retrotransposon Tca2

Background

In the C. albicans retrotransposon Tca2, the gag and pol ORFs are separated by a UGA stop codon, 3' of which is a potential RNA pseudoknot. It is unclear how the Tca2 gag UGA codon is bypassed to allow pol expression. However, in other retroelements, translational readthrough of the gag stop codon can be directed by its flanking sequence, including a 3' pseudoknot.

Results

The hypothesis was tested that in Tca2, gag stop codon flanking sequences direct translational readthrough and synthesis of a gag-pol fusion protein. Sequence from the Tca2 gag-UGA-pol junction (300 nt) was inserted between fused lacZ and luciferase (luc) genes in a Saccharomyces cerevisiae dual reporter construct. Although downstream of UGA, luc was expressed, but its expression was unaffected by inserting additional stop codons at the 3' end of lacZ. Luc expression was instead being driven by a previously unknown minor promoter activity within the gag-pol junction region. Evidence together indicated that junction sequence alone cannot direct UGA readthrough. Using reporter genes in C. albicans, the activities of this gag-pol junction promoter and the Tca2 long terminal repeat (LTR) promoter were compared. Of the two promoters, only the LTR promoter was induced by heat-shock, which also triggers retrotransposition. Tca2 pol protein, epitope-tagged in C. albicans to allow detection, was also heat-shock induced, indicating that pol proteins were expressed from a gag-UGA-pol RNA.

Conclusion

This is the first demonstration that the LTR promoter directs Tca2 pol protein expression, and that pol proteins are translated from a gag-pol RNA, which thus requires a mechanism for stop codon bypass. However, in contrast to most other retroelement and viral readthrough signals, immediate gag UGA-flanking sequences were insufficient to direct stop readthrough in S. cerevisiae, indicating non-canonical mechanisms direct gag UGA bypass in Tca2.

A yeast model for target-primed (non-LTR) retrotransposition

Background

Target-primed (non-LTR) retrotransposons, such as the human L1 element, are mobile genetic elements found in many eukaryotic genomes. They are often present in large numbers and their retrotransposition can cause mutations and genomic rearrangements. Despite their importance, many aspects of their replication are not well understood.

Results

We have developed a yeast model system for studying target-primed retrotransposons. This system uses the Zorro3 element from Candida albicans. A cloned copy of Zorro3, tagged with a retrotransposition indicator gene, retrotransposes at a high frequency when introduced into an appropriate C. albicans host strain. Retrotransposed copies of the tagged element exhibit similar features to the native copies, indicating that the natural retrotransposition pathway is being used. Retrotransposition is dependent on the products of the tagged element's own genes and is highly temperature-regulated. The new assay permits the analysis of the effects of specific mutations introduced into the cloned element.

Conclusion

This Zorro3 retrotransposition assay system complements previously available target-primed retrotransposition assays. Due to the relative simplicity of the growth, manipulation and analysis of yeast cells, the system should advance our understanding of target-primed retrotransposition.

Loss and fragmentation of chromosome 5 are major events linked to the adaptation of rad52-DeltaDelta strains of Candida albicans to sorbose

Candida albicans can adapt and grow on sorbose plates by losing one copy of Chr5. Since rad52 mutants of Saccharomyces cerevisiae lose chromosomes at a high rate, we have investigated the ability of C. albicans rad52 to adapt to sorbose. Carad52-DeltaDelta mutants generate Sou(+) strains earlier than wild-type but the final yield is lower, probably because they die at a higher rate in sorbose. As other strains of C. albicans, CAF2 and rad52-DeltaDelta derivatives generate Sou(+) strains by a loss of one copy of Chr5 about 75% of the time. In addition, rad52 strains were able to produce Sou(+) strains by a fragmentation/deletion event in one copy of Chr5, consisting of loss of a region adjacent to the right telomere. Finally, both CAF2 and rad52-DeltaDelta produced Sou(+) strains with two apparent full copies of Chr5, suggesting that additional genomic changes may also regulate adaptation to sorbose.

Developmental regulation of an adhesin gene during cellular morphogenesis in the fungal pathogen Candida albicans

Candida albicans expresses specific virulence traits that promote disease establishment and progression. These traits include morphological transitions between yeast and hyphal growth forms that are thought to contribute to dissemination and invasion and cell surface adhesins that promote attachment to the host. Here, we describe the regulation of the adhesin gene ALS3, which is expressed specifically during hyphal development in C. albicans. Using a combination of reporter constructs and regulatory mutants, we show that this regulation is mediated by multiple factors at the transcriptional level. The analysis of ALS3 promoter deletions revealed that this promoter contains two activation regions: one is essential for activation during hyphal development, while the second increases the amplitude of this activation. Further deletion analyses using the Renilla reniformis luciferase reporter delineate the essential activation region between positions -471 and -321 of the promoter. Further 5' or 3' deletions block activation. ALS3 transcription is repressed mainly by Nrg1 and Tup1, but Rfg1 contributes to this repression. Efg1, Tec1, and Bcr1 are essential for the transcriptional activation of ALS3, with Tec1 mediating its effects indirectly through Bcr1 rather than through the putative Tec1 sites in the ALS3 promoter. ALS3 transcription is not affected by Cph2, but Cph1 contributes to full ALS3 activation. The data suggest that multiple morphogenetic signaling pathways operate through the promoter of this adhesin gene to mediate its developmental regulation in this major fungal pathogen.

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.

Automated genotyping of Saccharomyces cerevisiae using the RiboPrinter®

This research note addresses the development of an automated molecular typing system for yeast. Specifically, our objectives were to generate specific probes for genotyping yeast with an automated fingerprinting system. We have adapted the RiboPrinter microbial characterization system for use with Saccharomyces cerevisiae yeast using alternative probes based on specific multi-copy gene families. Manual construction and labeling of probes proved successful in initial experiments. Results indicate that this method could be applied to food or clinical environments if the appropriate probes are developed.

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.

Functional genomics approaches for the identification and validation of antifungal drug targets

So far, antifungal drug discovery seems to have benefited little from the enormous advances in the field of genomics in the last decade. Although it has become clear that traditional drug screening is not delivering the long-awaited novel potent antifungals, little has been reported on efforts to use novel genome-based methodologies in the quest for new drugs acting on human pathogenic fungi. Although the market for a novel systemic and even topical broad-spectrum antifungal appears considerable, many large pharmaceutical companies have decided to scale back their activities in antifungal drug discovery. Here we report on some of the recent advances in genomics-based technologies that will allow us not only to identify and validate novel drug targets but hopefully also to discover active therapeutic agents. Novel drug targets have already been found by 'en masse' gene inactivation strategies (e.g. using antisense RNA inhibition). In addition, genome expression profiling using DNA microarrays helps to assign gene function but also to understand better the mechanism of action of known drugs (e.g. itraconazole) and to elucidate how new drug candidates work. No doubt, we have a long way to go just to catch up with the advances made in other therapeutic areas, but all tools are at hand to derive practical benefits from the genomics revolution. The next few years should prove a very exciting time in the history of antifungal drug discovery.

Current awareness on yeast

In order to keep subscribers up-to-date with the latest developments in their field, this current awareness service is provided by John Wiley & Sons and contains newly-published material on yeasts. Each bibliography is divided into 10 sections. 1 Books, Reviews & Symposia; 2 General; 3 Biochemistry; 4 Biotechnology; 5 Cell Biology; 6 Gene Expression; 7 Genetics; 8 Physiology; 9 Medical Mycology; 10 Recombinant DNA Technology. Within each section, articles are listed in alphabetical order with respect to author. If, in the preceding period, no publications are located relevant to any one of these headings, that section will be omitted. (3 weeks journals - search completed 5th. Dec. 2001)