Extrachromosomal Circular DNA: An Emerging Potential Biomarker for Inflammatory Bowel Diseases?

Inflammatory bowel disease (IBD) comprising ulcerative colitis and Crohn's disease is a chronic immune-mediated disease which affects the gastrointestinal tract with a relapsing and remitting course, causing lifelong morbidity. IBD pathogenesis is determined by multiple factors including genetics, immune and microbial factors, and environmental factors. Although therapy options are expanding, remission rates are unsatisfiable, and together with the disease course, response to therapy remains unpredictable. Therefore, the identification of biomarkers that are predictive for the disease course and response to therapy is a significant challenge. Extrachromosomal circular DNA (eccDNA) fragments exist in all tissue tested so far. These fragments, ranging in length from a few hundreds of base pairs to mega base pairs, have recently gained more interest due to technological advances. Until now, eccDNA has mainly been studied in relation to cancer due to its ability to act as an amplification site for oncogenes and drug resistance genes. However, eccDNA could also play an important role in inflammation, expressed both locally in the- involved tissue and at distant sites. Here, we review the current evidence on the molecular mechanisms of eccDNA and its role in inflammation and IBD. Additionally, the potential of eccDNA as a tissue or plasma marker for disease severity and/or response to therapy is evaluated.

Methods for the purification and detection of single nucleotide KRAS mutations on extrachromosomal circular DNA in human plasma

BACKGROUNDS

Despite recent advances, many cancers are still detected too late for curative treatment. There is, therefore, a need for the development of new diagnostic methods and biomarkers. One approach may arise from the detection of extrachromosomal circular DNA (eccDNA), which is part of cell-free DNA in human plasma.

AIMS

First, we assessed and compared two methods for the purification of eccDNA from plasma. Second, we tested for an easy diagnostic application of eccDNA liquid biopsy-based assays.

MATERIALS & METHODS

For the comparison we tested a solid-phase silica purification method and a phenol/chloroform method with salt precipitation. For the diagnostic application of eccDNA we developed and tested a qPCR primer-based SNP detection system, for the detection of two well-established cancer-causing KRAS mutations (G12V and G12R) on circular DNA. This investigation was supported by purifying, sequencing, and analysing clinical plasma samples for eccDNAs containing KRAS mutant alleles in 0.5 mL plasma from 16 pancreatic ductal adenocarcinoma patients and 19 healthy controls.

RESULTS

In our method comparison we observed, that following exonuclease treatment a lower eccDNA yield was found for the phenol/chloroform method (15.7%-26.7%) compared with the solid-phase purification approach (47.8%-65.9%). For the diagnostic application of eccDNA tests, the sensitivity of the tested qPCR assay only reached ~10-3 in a background of 105 wild type (wt) KRAS circular entities, which was not improved by general amplification or primer-based inhibition of wt KRAS amplification. Furthermore, we did not detect eccDNA containing KRAS in any of the clinical samples.

DISCUSSION

A potential explanation for our inability to detect any KRAS mutations in the clinical samples may be related to the general low abundance of eccDNA in plasma.

CONCLUSION

Taken together our results provide a benchmark for eccDNA purification methods while raising the question of what is required for the optimal fast and sensitive detection of SNP mutations on eccDNA with greater sensitivity than primer-based qPCR detection.

Variation of extrachromosomal circular DNA in cancer cell lines

The presence of oncogene carrying eccDNAs is strongly associated with carcinogenesis and poor patient survival. Tumour biopsies and in vitro cancer cell lines are frequently utilized as models to investigate the role of eccDNA in cancer. However, eccDNAs are often lost during the in vitro growth of cancer cell lines, questioning the reproducibility of studies utilizing cancer cell line models. Here, we conducted a comprehensive analysis of eccDNA variability in seven cancer cell lines (MCA3D, PDV, HaCa4, CarC, MIA-PaCa-2, AsPC-1, and PC-3). We compared the content of unique eccDNAs between triplicates of each cell line and found that the number of unique eccDNA is specific to each cell line, while the eccDNA sequence content varied greatly among triplicates (∼ 0-1% eccDNA coordinate commonality). In the PC-3 cell line, we found that the large eccDNA (ecDNA) with MYC is present in high-copy number in an NCI cell line isolate but not present in ATCC isolates. Together, these results reveal that the sequence content of eccDNA is highly variable in cancer cell lines. This highlights the importance of testing cancer cell lines before use, and to enrich for subclones in cell lines with the desired eccDNA to get relatively pure population for studying the role of eccDNA in cancer.

Circular and Circulating DNA in Inflammatory Bowel Disease: From Pathogenesis to Potential Molecular Therapies

Inflammatory bowel diseases (IBD), including Crohn's Disease (CD) and Ulcerative Colitis (UC) are chronic multifactorial disorders which affect the gastrointestinal tract with variable extent. Despite extensive research, their etiology and exact pathogenesis are still unknown. Cell-free DNAs (cfDNAs) are defined as any DNA fragments which are free from the origin cell and able to circulate into the bloodstream with or without microvescicles. CfDNAs are now being increasingly studied in different human diseases, like cancer or inflammatory diseases. However, to date it is unclear how IBD etiology is linked to cfDNAs in plasma. Extrachromosomal circular DNA (eccDNA) are non-plasmidic, nuclear, circular and closed DNA molecules found in all eukaryotes tested. CfDNAs appear to play an important role in autoimmune diseases, inflammatory processes, and cancer; recently, interest has also grown in IBD, and their role in the pathogenesis of IBD has been suggested. We now suggest that eccDNAs also play a role in IBD. In this review, we have comprehensively collected available knowledge in literature regarding cfDNA, eccDNA, and structures involving them such as neutrophil extracellular traps and exosomes, and their role in IBD. Finally, we focused on old and novel potential molecular therapies and drug delivery systems, such as nanoparticles, for IBD treatment.

Did circular DNA shape the evolution of mammalian genomes?

Extrachromosomal circular DNA (eccDNA) can shape the genomes of somatic cells, but how it impacts genomes across generations is largely unexplored. We propose that genomes can rearrange via circular intermediates across generations and show that up to 6% of a mammalian genome can have changed gene order through eccDNA.

CReSIL: accurate identification of extrachromosomal circular DNA from long-read sequences

Extrachromosomal circular DNA (eccDNA) of chromosomal origin is found in many eukaryotic species and cell types, including cancer, where eccDNAs with oncogenes drive tumorigenesis. Most studies of eccDNA employ short-read sequencing for their identification. However, short-read sequencing cannot resolve the complexity of genomic repeats, which can lead to missing eccDNA products. Long-read sequencing technologies provide an alternative to constructing complete eccDNA maps. We present a software suite, Construction-based Rolling-circle-amplification for eccDNA Sequence Identification and Location (CReSIL), to identify and characterize eccDNA from long-read sequences. CReSIL's performance in identifying eccDNA, with a minimum F1 score of 0.98, is superior to the other bioinformatic tools based on simulated data. CReSIL provides many useful features for genomic annotation, which can be used to infer eccDNA function and Circos visualization for eccDNA architecture investigation. We demonstrated CReSIL's capability in several long-read sequencing datasets, including datasets enriched for eccDNA and whole genome datasets from cells containing large eccDNA products. In conclusion, the CReSIL suite software is a versatile tool for investigating complex and simple eccDNA in eukaryotic cells.

A unifying model for extrachromosomal circular DNA load in eukaryotic cells

Extrachromosomal circular DNA (eccDNA) with exons and whole genes are common features of eukaryotic cells. Work from especially tumours and the yeast Saccharomyces cerevisiae has revealed that eccDNA can provide large selective advantages and disadvantages. Besides the phenotypic effect due to expression of an eccDNA fragment, eccDNA is different from other mutations in that it is released from 1:1 segregation during cell division. This means that eccDNA can quickly change copy number, pickup secondary mutations and reintegrate into a chromosome to establish substantial genetic variation that could not have evolved via canonical mechanisms. We propose a unifying 5-factor model for conceptualizing the eccDNA load of a eukaryotic cell, emphasizing formation, replication, segregation, selection and elimination. We suggest that the magnitude of these sequential events and their interactions determine the copy number of eccDNA in mitotically dividing cells. We believe that our model will provide a coherent framework for eccDNA research, to understand its biology and the factors that can be manipulated to modulate eccDNA load in eukaryotic cells.

Circle-Seq reveals genomic and disease-specific hallmarks in urinary cell-free extrachromosomal circular DNAs

BACKGROUND

Extrachromosomal circular deoxyribonucleic acid (eccDNA) is evolving as a valuable biomarker, while little is known about its presence in urine.

METHODS

Here, we report the discovery and analysis of urinary cell-free eccDNAs (ucf-eccDNAs) in healthy controls and patients with advanced chronic kidney disease (CKD) by Circle-Seq.

RESULTS

Millions of unique ucf-eccDNAs were identified and comprehensively characterised. The ucf-eccDNAs are GC-rich. Most ucf-eccDNAs are less than 1000 bp and are enriched in four pronounced peaks at 207, 358, 553 and 732 bp. Analysis of the genomic distribution of ucf-eccDNAs shows that eccDNAs are found on all chromosomes but enriched on chromosomes 17, 19 and 20 with a high density of protein-coding genes, CpG islands, short interspersed transposable elements (SINEs) and simple repeat elements. Analysis of eccDNA junction sequences further suggests that microhomology and palindromic repeats might be involved in eccDNA formation. The ucf-eccDNAs in CKD patients are significantly higher than those in healthy controls. Moreover, eccDNA with miRNA genes is highly enriched in CKD ucf-eccDNA.

CONCLUSIONS

This work discovers and provides the first deep characterisation of ucf-eccDNAs and suggests ucf-eccDNA as a valuable noninnvasive biomarker for urogenital disorder diagnosis and monitoring.

Extrachromosomal circular DNA in cancer: history, current knowledge, and methods

Extrachromosomal circular DNA (eccDNA) is a closed-circle, nuclear, nonplasmid DNA molecule found in all tested eukaryotes. eccDNA plays important roles in cancer pathogenesis, evolution of tumor heterogeneity, and therapeutic resistance. It is known under many names, including very large cancer-specific circular extrachromosomal DNA (ecDNA), which carries oncogenes and is often amplified in cancer cells. Our understanding of eccDNA has historically been limited and fragmented. To provide better a context of new and previous research on eccDNA, in this review we give an overview of the various names given to eccDNA at different times. We describe the different mechanisms for formation of eccDNA and the methods used to study eccDNA thus far. Finally, we explore the potential clinical value of eccDNA.

Circular DNA in the human germline and its association with recombination

Extrachromosomal circular DNA (eccDNA) is common in somatic tissue, but its existence and effects in the human germline are unexplored. We used microscopy, long-read DNA sequencing, and new analytic methods to document thousands of eccDNAs from human sperm. EccDNAs derived from all genomic regions and mostly contained a single DNA fragment, although some consisted of multiple fragments. The generation of eccDNA inversely correlates with the meiotic recombination rate, and chromosomes with high coding-gene density and Alu element abundance form the least eccDNA. Analysis of insertions in human genomes further indicates that eccDNA can persist in the human germline when the circular molecules reinsert themselves into the chromosomes. Our results suggest that eccDNA has transient and permanent effects on the germline. They explain how differences in the physical and genetic map might arise and offer an explanation of how Alu elements coevolved with genes to protect genome integrity against deleterious mutations producing eccDNA.

Targeted removal of mitochondrial DNA from mouse and human extrachromosomal circular DNA with CRISPR-Cas9

Extrachromosomal circular DNA (eccDNA) of chromosomal origin is common in eukaryotic cells. Amplification of oncogenes on large eccDNA (ecDNA) can drive biological processes such as tumorigenesis, and identification of eccDNA by sequencing after removal of chromosomal DNA is therefore important for understanding their impact on the expressed phenotype. However, the circular mitochondrial DNA (mtDNA) might challenge the detection of eccDNA because the average somatic cell has hundreds of copies of mtDNA. Here we show that 61.2–99.5% of reads from eccDNA-enriched samples correspond to mtDNA in mouse tissues. We have developed a method to selectively remove mtDNA from total circular DNA by CRISPR/Cas9 guided cleavage of mtDNA with one single-guide RNA (sgRNA) or two sgRNAs followed by exonuclease degradation of the linearized mtDNA. Sequencing revealed that mtDNA reads were 85.9% ± 12.6% removed from eccDNA of 9 investigated mouse tissues. CRISPR/Cas9 cleavage also efficiently removed mtDNA from a human HeLa cell line and colorectal cancer samples. We identified up to 14 times more, and also larger eccDNA in CRISPR/Cas9 treated colorectal cancer samples than in untreated samples. We foresee that the method can be applied to effectively remove mtDNA from any eukaryotic species to obtain higher eccDNA yields.

Isolation, characterization, and genome assembly of Barnettozyma botsteinii sp. nov. and novel strains of Kurtzmaniella quercitrusa isolated from the intestinal tract of the termite Macrotermes bellicosus

Bioconversion of hemicelluloses into simpler sugars leads to the production of a significant amount of pentose sugars, such as d-xylose. However, efficient utilization of pentoses by conventional yeast production strains remains challenging. Wild yeast strains can provide new industrially relevant characteristics and efficiently utilize pentose sugars. To explore this strategy, we isolated gut-residing yeasts from the termite Macrotermes bellicosus collected in Comoé National Park, Côte d'Ivoire. The yeasts were classified through their Internal Transcribed Spacer/Large Subunit sequence, and their genomes were sequenced and annotated. We identified a novel yeast species, which we name Barnettozyma botsteinii sp. nov. 1118T (MycoBank: 833563, CBS 16679T and IBT 710) and two new strains of Kurtzmaniella quercitrusa: var. comoensis (CBS 16678, IBT 709) and var. filamentosus (CBS 16680, IBT 711). The two K. quercitrusa strains grow 15% faster on synthetic glucose medium than Saccharomyces cerevisiae CEN.PKT in acidic conditions (pH = 3.2) and both strains grow on d-xylose as the sole carbon source at a rate of 0.35 h-1. At neutral pH, the yeast form of K. quercitrusa var. filamentosus, but not var. comoensis, switched to filamentous growth in a carbon source-dependent manner. Their genomes are 11.0-13.2 Mb in size and contain between 4888 and 5475 predicted genes. Together with closely related species, we did not find any relationship between gene content and ability to grow on xylose. Besides its metabolism, K. quercitrusa var. filamentosus has a large potential as a production organism, because of its capacity to grow at low pH and to undergo a dimorphic shift.

Replicative aging is associated with loss of genetic heterogeneity from extrachromosomal circular DNA in Saccharomyces cerevisiae

Circular DNA can arise from all parts of eukaryotic chromosomes. In yeast, circular ribosomal DNA (rDNA) accumulates dramatically as cells age, however little is known about the accumulation of other chromosome-derived circles or the contribution of such circles to genetic variation in aged cells. We profiled circular DNA in Saccharomyces cerevisiae populations sampled when young and after extensive aging. Young cells possessed highly diverse circular DNA populations but 94% of the circular DNA were lost after ∼15 divisions, whereas rDNA circles underwent massive accumulation to >95% of circular DNA. Circles present in both young and old cells were characterized by replication origins including circles from unique regions of the genome and repetitive regions: rDNA and telomeric Y' regions. We further observed that circles can have flexible inheritance patterns: [HXT6/7circle] normally segregates to mother cells but in low glucose is present in up to 50% of cells, the majority of which must have inherited this circle from their mother. Interestingly, [HXT6/7circle] cells are eventually replaced by cells carrying stable chromosomal HXT6 HXT6/7 HXT7 amplifications, suggesting circular DNAs are intermediates in chromosomal amplifications. In conclusion, the heterogeneity of circular DNA offers flexibility in adaptation, but this heterogeneity is remarkably diminished with age.

Near-Random Distribution of Chromosome-Derived Circular DNA in the Condensed Genome of Pigeons and the Larger, More Repeat-Rich Human Genome

Extrachromosomal circular DNA (eccDNA) elements of chromosomal origin are known to be common in a number of eukaryotic species. However, it remains to be addressed whether genomic features such as genome size, the load of repetitive elements within a genome, and/or animal physiology affect the number of eccDNAs. Here, we investigate the distribution and numbers of eccDNAs in a condensed and less repeat-rich genome compared with the human genome, using Columba livia domestica (domestic rock pigeon) as a model organism. By sequencing eccDNA in blood and breast muscle from three pigeon breeds at various ages and with different flight behavior, we characterize 30,000 unique eccDNAs. We identify genomic regions that are likely hotspots for DNA circularization in breast muscle, including genes involved in muscle development. We find that although eccDNA counts do not correlate with the biological age in pigeons, the number of unique eccDNAs in a nonflying breed (king pigeons) is significantly higher (9-fold) than homing pigeons. Furthermore, a comparison between eccDNA from skeletal muscle in pigeons and humans reveals ∼9-10 times more unique eccDNAs per human nucleus. The fraction of eccDNA sequences, derived from repetitive elements, exist in proportions to genome content, that is, human 72.4% (expected 52.5%) and pigeon 8.7% (expected 5.5%).

Overall, our results support that eccDNAs are common in pigeons, that the amount of unique eccDNA types per nucleus can differ between species as well as subspecies, and suggest that eccDNAs from repeats are found in proportions relative to the content of repetitive elements in a genome.

Sensitive detection of circular DNAs at single-nucleotide resolution using guided realignment of partially aligned reads

BACKGROUND

Circular DNA has recently been identified across different species including human normal and cancerous tissue, but short-read mappers are unable to align many of the reads crossing circle junctions hence limiting their detection from short-read sequencing data.

RESULTS

Here, we propose a new method, Circle-Map that guides the realignment of partially aligned reads using information from discordantly mapped reads to map the short unaligned portions using a probabilistic model. We compared Circle-Map to similar up-to-date methods for circular DNA and RNA detection and we demonstrate how the approach implemented in Circle-Map dramatically increases sensitivity for detection of circular DNA on both simulated and real data while retaining high precision.

CONCLUSION

Circle-Map is an easy-to-use command line tool that implements the required pipeline to accurately detect circular DNA from circle enriched next generation sequencing experiments. Circle-Map is implemented in python3.6 and it is freely available at https://github.com/iprada/Circle-Map.

To Be or Not to Be: Circular RNAs or mRNAs From Circular DNAs?

In recent years, there has been a growing interest in circular RNAs (circRNAs) since they are involved in a wide spectrum of cellular functions that might have a large impact on phenotype and disease. CircRNAs are mainly recorded by RNA-Seq and computational methods focused on the detection of back-splicing junction sequences considered the diagnostic feature of circRNAs. While some protocols remove linear RNA prior to sequencing, many have characterized circRNAs by sorting through total RNA sequencing data without excluding the possibility that some linear RNA can provide the same signal as a circRNA. Recent studies have revealed that circular DNAs of chromosomal origin are common in eukaryotic genomes and that they can be transcribed. Transcription events across the junction of circular DNAs would result in a transcript with a junction similar to those present in circRNAs. Therefore, in this report, we want to draw attention to transcripts from such circular DNAs both as an interesting new player in the transcriptome and also as a confounding factor that must be taken into account when studying circRNAs.

Multicellular group formation in Saccharomyces cerevisiae

Understanding how and why cells cooperate to form multicellular organisms is a central aim of evolutionary biology. Multicellular groups can form through clonal development (where daughter cells stick to mother cells after division) or by aggregation (where cells aggregate to form groups). These different ways of forming groups directly affect relatedness between individual cells, which in turn can influence the degree of cooperation and conflict within the multicellular group. It is hard to study the evolution of multicellularity by focusing only on obligately multicellular organisms, like complex animals and plants, because the factors that favour multicellular cooperation cannot be disentangled, as cells cannot survive and reproduce independently. We support the use of Saccharomyces cerevisiae as an ideal model for studying the very first stages of the evolution of multicellularity. This is because it can form multicellular groups both clonally and through aggregation and uses a family of proteins called ‘flocculins’ that determine the way in which groups form, making it particularly amenable to laboratory experiments. We briefly review current knowledge about multicellularity in S. cerevisiae and then propose a framework for making predictions about the evolution of multicellular phenotypes in yeast based on social evolution theory. We finish by explaining how S. cerevisiae is a particularly useful experimental model for the analysis of open questions concerning multicellularity.

RNAi as a Tool to Study Virulence in the Pathogenic Yeast Candida glabrata

The yeast Candida glabrata is a major opportunistic pathogen causing mucosal and systemic infections in humans. Systemic infections caused by this yeast have high mortality rates and are difficult to treat due to this yeast’s intrinsic and frequently adapting antifungal resistance. To understand and treat C. glabrata infections, it is essential to investigate the molecular basis of C. glabrata virulence and resistance. We established an RNA interference (RNAi) system in C. glabrata by expressing the Dicer and Argonaute genes from Saccharomyces castellii (a budding yeast with natural RNAi). Our experiments with reporter genes and putative virulence genes showed that the introduction of RNAi resulted in 30 and 70% gene-knockdown for the construct-types antisense and hairpin, respectively. The resulting C. glabrata RNAi strain was used for the screening of a gene library for new virulence-related genes. Phenotypic profiling with a high-resolution quantification of growth identified genes involved in the maintenance of cell integrity, antifungal drugs, and ROS resistance. The genes identified by this approach are promising targets for the treatment of C. glabrata infections.

CRISPR-C: circularization of genes and chromosome by CRISPR in human cells

Extrachromosomal circular DNA (eccDNA) and ring chromosomes are genetic alterations found in humans with genetic disorders. However, there is a lack of genetic engineering tools to recapitulate and study the biogenesis of eccDNAs. Here, we created a dual-fluorescence biosensor cassette, which upon the delivery of pairs of CRISPR/Cas9 guide RNAs, CRISPR-C, allows us to study the biogenesis of a specific fluorophore expressing eccDNA in human cells. We show that CRISPR-C can generate functional eccDNA, using the novel eccDNA biosensor system. We further reveal that CRISPR-C also can generate eccDNAs from intergenic and genic loci in human embryonic kidney 293T cells and human mammary fibroblasts. EccDNAs mainly forms by end-joining mediated DNA-repair and we show that CRISPR-C is able to generate endogenous eccDNAs in sizes from a few hundred base pairs and ranging up to 207 kb. Even a 47.4 megabase-sized ring chromosome 18 can be created by CRISPR-C. Our study creates a new territory for CRISPR gene editing and highlights CRISPR-C as a useful tool for studying the cellular impact, persistence and function of eccDNAs.

Regulation of apoptosis and autophagy in mouse and human skeletal muscle with aging and lifelong exercise training

Exercise training has been reported to prevent the age-induced decline in muscle mass and fragmentation of mitochondria, as well as to affect autophagy and mitophagy. The interaction between these pathways during aging as well as the similarity between such changes in human and mouse skeletal muscle is however not fully understood. Therefore the aim of the present study was to test the hypothesis that cellular degradation pathways, including apoptosis, autophagy and mitophagy are coordinately regulated in mouse and human skeletal muscle during aging and lifelong exercise training through a PGC-1α-p53 axis. Muscle samples were obtained from young untrained, aged untrained and aged lifelong exercise trained men, and from whole-body PGC-1α knockout mice and their littermate controls that were either lifelong exercise trained or sedentary young and aged. Lifelong exercise training prevented the aging-induced reduction in PGC-1α, p53 and p21 mRNA as well as the increase in LC3II and BNIP3 protein in mouse skeletal muscle, while aging decreased the BAX/Bcl-2 ratio, LC3I and BAX protein in mouse skeletal muscle without effects of lifelong exercise training. In humans, aging was associated with reduced PGC-1α mRNA as well as decreased p62 and p21 protein in skeletal muscle, while lifelong exercise training increased BNIP3 protein and decreased p53 mRNA. In conclusion, there was a divergent regulation of autophagy and apoptosis in mouse muscle with aging and lifelong exercise training, whereas healthy aged human skeletal muscle seemed rather robust to changes in apoptosis, autophagy and mitophagy markers compared with mouse muscle at the investigated age.

Clonal yeast biofilms can reap competitive advantages through cell differentiation without being obligatorily multicellular

How differentiation between cell types evolved is a fundamental question in biology, but few studies have explored single-gene phenotypes that mediate first steps towards division of labour with selective advantage for groups of cells. Here, we show that differential expression of the FLO11 gene produces stable fractions of Flo11⁺ and Flo11^- cells in clonal Saccharomyces cerevisiae biofilm colonies on medium with intermediate viscosity. Differentiated Flo11^+/- colonies, consisting of adhesive and non-adhesive cells, obtain a fourfold growth advantage over undifferentiated colonies by overgrowing glucose resources before depleting them, rather than depleting them while they grow as undifferentiated Flo11^- colonies do. Flo11^+/- colonies maintain their structure and differentiated state by switching non-adhesive cells to adhesive cells with predictable probability. Mixtures of Flo11⁺ and Flo11^- cells from mutant strains that are unable to use this epigenetic switch mechanism produced neither integrated colonies nor growth advantages, so the condition-dependent selective advantages of differentiated FLO11 expression can only be reaped by clone-mate cells. Our results show that selection for cell differentiation in clonal eukaryotes can evolve before the establishment of obligate undifferentiated multicellularity, and without necessarily leading to more advanced organizational complexity.

Genome-wide Purification of Extrachromosomal Circular DNA from Eukaryotic Cells

Extrachromosomal circular DNAs (eccDNAs) are common genetic elements in Saccharomyces cerevisiae and are reported in other eukaryotes as well. EccDNAs contribute to genetic variation among somatic cells in multicellular organisms and to evolution of unicellular eukaryotes. Sensitive methods for detecting eccDNA are needed to clarify how these elements affect genome stability and how environmental and biological factors induce their formation in eukaryotic cells. This video presents a sensitive eccDNA-purification method called Circle-Seq. The method encompasses column purification of circular DNA, removal of remaining linear chromosomal DNA, rolling-circle amplification of eccDNA, deep sequencing, and mapping. Extensive exonuclease treatment was required for sufficient linear chromosomal DNA degradation. The rolling-circle amplification step by φ29 polymerase enriched for circular DNA over linear DNA. Validation of the Circle-Seq method on three S. cerevisiae CEN.PK populations of 10(10) cells detected hundreds of eccDNA profiles in sizes larger than 1 kilobase. Repeated findings of ASP3-1, COS111, CUP1, RSC30, HXT6, HXT7 genes on circular DNA in both S288c and CEN.PK suggests that DNA circularization is conserved between strains at these loci. In sum, the Circle-Seq method has broad applicability for genome-scale screening for eccDNA in eukaryotes as well as for detecting specific eccDNA types.

Saccharomyces cerevisiae biofilm tolerance towards systemic antifungals depends on growth phase

Background

Biofilm-forming Candida species cause infections that can be difficult to eradicate, possibly because of antifungal drug tolerance mechanisms specific to biofilms. In spite of decades of research, the connection between biofilm and drug tolerance is not fully understood.

Results

We used Saccharomyces cerevisiae as a model for drug susceptibility of yeast biofilms. Confocal laser scanning microscopy showed that S. cerevisiae and C. glabrata form similarly structured biofilms and that the viable cell numbers were significantly reduced by treatment of mature biofilms with amphotericin B but not voriconazole, flucytosine, or caspofungin. We showed that metabolic activity in yeast biofilm cells decreased with time, as visualized by FUN-1 staining, and mature, 48-hour biofilms contained cells with slow metabolism and limited growth. Time-kill studies showed that in exponentially growing planktonic cells, voriconazole had limited antifungal activity, flucytosine was fungistatic, caspofungin and amphotericin B were fungicidal. In growth-arrested cells, only amphotericin B had antifungal activity. Confocal microscopy and colony count viability assays revealed that the response of growing biofilms to antifungal drugs was similar to the response of exponentially growing planktonic cells. The response in mature biofilm was similar to that of non-growing planktonic cells. These results confirmed the importance of growth phase on drug efficacy.

Conclusions

We showed that in vitro susceptibility to antifungal drugs was independent of biofilm or planktonic growth mode. Instead, drug tolerance was a consequence of growth arrest achievable by both planktonic and biofilm populations. Our results suggest that efficient strategies for treatment of yeast biofilm might be developed by targeting of non-dividing cells.

Genetic basis for Saccharomyces cerevisiae biofilm in liquid medium

Biofilm-forming microorganisms switch between two forms: free-living planktonic and sessile multicellular. Sessile communities of yeast biofilms in liquid medium provide a primitive example of multicellularity and are clinically important because biofilms tend to have other growth characteristics than free-living cells. We investigated the genetic basis for yeast, Saccharomyces cerevisiae, biofilm on solid surfaces in liquid medium by screening a comprehensive deletion mutant collection in the Σ1278b background and found 71 genes that were essential for biofilm development. Quantitative northern blots further revealed that AIM1, ASG1, AVT1, DRN1, ELP4, FLO8, FMP10, HMT1, KAR5, MIT1, MRPL32, MSS11, NCP1, NPR1, PEP5, PEX25, RIM8, RIM101, RGT1, SNF8, SPC2, STB6, STP22, TEC1, VID24, VPS20, VTC3, YBL029W, YBL029C-A, YFL054C, YGR161W-C, YIL014C-A, YIR024C, YKL151C, YNL200C, YOR034C-A, and YOR223W controlled biofilm through FLO11 induction. Almost all deletion mutants that were unable to form biofilms in liquid medium also lost the ability to form surface-spreading biofilm colonies (mats) on agar and 69% also lost the ability to grow invasively. The protein kinase A isoform Tpk3p functioned specifically in biofilm and mat formation. In a tpk3 mutant, transcription of FLO11 was induced three-fold compared with wild-type, but biofilm development and cell-cell adhesion was absent, suggesting that Tpk3p regulates FLO11 positive posttranscriptionally and negative transcriptionally.The study provides a resource of biofilm-influencing genes for additional research on biofilm development and suggests that the regulation of FLO11 is more complex than previously anticipated.

The synthetic amphipathic peptidomimetic LTX109 is a potent fungicide that disturbs plasma membrane integrity in a sphingolipid dependent manner

The peptidomimetic LTX109 (arginine-tertbutyl tryptophan-arginine-phenylethan) was previously shown to have antibacterial properties. Here, we investigated the activity of this novel antimicrobial peptidomimetic on the yeast Saccharomyces cerevisiae. We found that LTX109 was an efficient fungicide that killed all viable cells in an exponentially growing population as well as a large proportion of cells in biofilm formed on an abiotic surface. LTX109 had similar killing kinetics to the membrane-permeabilizing fungicide amphotericin B, which led us to investigate the ability of LTX109 to disrupt plasma membrane integrity. S. cerevisiae cells exposed to a high concentration of LTX109 showed rapid release of potassium and amino acids, suggesting that LTX109 acted by destabilizing the plasma membrane. This was supported by the finding that cells were permeable to the fluorescent nucleic acid stain SYTOX Green after a few minutes of LTX109 treatment. We screened a haploid S. cerevisiae gene deletion library for mutants resistant to LTX109 to uncover potential molecular targets. Eight genes conferred LTX109 resistance when deleted and six were involved in the sphingolipid biosynthetic pathway (SUR1, SUR2, SKN1, IPT1, FEN1 and ORM2). The involvement of all of these genes in the biosynthetic pathway for the fungal-specific lipids mannosylinositol phosphorylceramide (MIPC) and mannosyl di-(inositol phosphoryl) ceramide (M(IP)2C) suggested that these lipids were essential for LTX109 sensitivity. Our observations are consistent with a model in which LTX109 kills S. cerevisiae by nonspecific destabilization of the plasma membrane through direct or indirect interaction with the sphingolipids.

A model for generating several adaptive phenotypes from a single genetic event: Saccharomyces cerevisiae GAP1 as a potential bet-hedging switch

Microbial populations adapt to environmental fluctuations through random switching of fitness-related traits in individual cells. This increases the likelihood that a subpopulation will be adaptive in a future milieu. However, populations are particularly challenged when several environment factors change simultaneously. We suggest that a population can rapidly adapt to multiple environmental changes if individual members stochastically flip a hub-switch that controls a set of adaptive phenotypes in a single event. This mechanism of coupling phenotypic outcomes via a hub-switch can protect a population against large fluctuations in size. Here we report that the general amino acid transporter Gap1 is a potential hub-switch. The GAP1 gene is flanked by two direct repeats that can lead to GAP1 deletions (∆gap1) and a self-replicating GAP1 circle. Thus, an isogenic GAP1 population can differentiate into two variant, reversible genotypes, ∆gap1 or GAP1 (circle). These subpopulations have different phenotypic advantages. A ∆gap1 population has a selective advantage on allantoin or ammonium as a nitrogen source and high stress tolerance. Advantages of the GAP1 population include amino acid uptake, fast energy recruitment by trehalose mobilization, and in some cases, adherent biofilm growth. Our proposed model of a hub-switch locus enhances the bet-hedging model of population dynamics.

Amino acid transporter genes are essential for FLO11-dependent and FLO11-independent biofilm formation and invasive growth in Saccharomyces cerevisiae

Amino acids can induce yeast cell adhesion but how amino acids are sensed and signal the modulation of the FLO adhesion genes is not clear. We discovered that the budding yeast Saccharomyces cerevisiae CEN.PK evolved invasive growth ability under prolonged nitrogen limitation. Such invasive mutants were used to identify amino acid transporters as regulators of FLO11 and invasive growth. One invasive mutant had elevated levels of FLO11 mRNA and a Q320STOP mutation in the SFL1 gene that encodes a protein kinase A pathway regulated repressor of FLO11. Glutamine-transporter genes DIP5 and GNP1 were essential for FLO11 expression, invasive growth and biofilm formation in this mutant. Invasive growth relied on known regulators of FLO11 and the Ssy1-Ptr3-Ssy5 complex that controls DIP5 and GNP1, suggesting that Dip5 and Gnp1 operates downstream of the Ssy1-Ptr3-Ssy5 complex for regulation of FLO11 expression in a protein kinase A dependent manner. The role of Dip5 and Gnp1 appears to be conserved in the S. cerevisiae strain ∑1278b since the dip5 gnp1 ∑1278b mutant showed no invasive phenotype. Secondly, the amino acid transporter gene GAP1 was found to influence invasive growth through FLO11 as well as other FLO genes. Cells carrying a dominant loss-of-function PTR3(647::CWNKNPLSSIN) allele had increased transcription of the adhesion genes FLO1, 5, 9, 10, 11 and the amino acid transporter gene GAP1. Deletion of GAP1 caused loss of FLO11 expression and invasive growth. However, deletions of FLO11 and genes encoding components of the mitogen-activated protein kinase pathway or the protein kinase A pathway were not sufficient to abolish invasive growth, suggesting involvement of other FLO genes and alternative pathways. Increased intracellular amino acid pools in the PTR3(647::CWNKNPLSSIN)-containing strain opens the possibility that Gap1 regulates the FLO genes through alteration of the amino acid pool sizes.

Pseudomonas aeruginosa and Saccharomyces cerevisiae biofilm in flow cells

Many microbial cells have the ability to form sessile microbial communities defined as biofilms that have altered physiological and pathological properties compared to free living microorganisms. Biofilms in nature are often difficult to investigate and reside under poorly defined conditions¹. Using a transparent substratum it is possible to device a system where simple biofilms can be examined in a non-destructive way in real-time: here we demonstrate the assembly and operation of a flow cell model system, for in vitro 3D studies of microbial biofilms generating high reproducibility under well-defined conditions^2,3.

The system consists of a flow cell that serves as growth chamber for the biofilm. The flow cell is supplied with nutrients and oxygen from a medium flask via a peristaltic pump and spent medium is collected in a waste container. This construction of the flow system allows a continuous supply of nutrients and administration of e.g. antibiotics with minimal disturbance of the cells grown in the flow chamber. Moreover, the flow conditions within the flow cell allow studies of biofilm exposed to shear stress. A bubble trapping device confines air bubbles from the tubing which otherwise could disrupt the biofilm structure in the flow cell.

The flow cell system is compatible with Confocal Laser Scanning Microscopy (CLSM) and can thereby provide highly detailed 3D information about developing microbial biofilms. Cells in the biofilm can be labeled with fluorescent probes or proteins compatible with CLSM analysis. This enables online visualization and allows investigation of niches in the developing biofilm. Microbial interrelationship, investigation of antimicrobial agents or the expression of specific genes, are of the many experimental setups that can be investigated in the flow cell system.

The roles of galactitol, galactose-1-phosphate, and phosphoglucomutase in galactose-induced toxicity in Saccharomyces cerevisiae

The uptake and catabolism of galactose by the yeast Saccharomyces cerevisiae is much lower than for glucose and fructose, and in applications of this yeast for utilization of complex substrates that contain galactose, for example, lignocellulose and raffinose, this causes prolonged fermentations. Galactose is metabolized via the Leloir pathway, and besides the industrial interest in improving the flux through this pathway it is also of medical relevance to study the Leloir pathway. Thus, genetic disorders in the genes encoding galactose-1-phosphate uridylyltransferase or galactokinase result in galactose toxicity both in patients with galactosemia and in yeast. In order to elucidate galactose related toxicity, which may explain the low uptake and catabolic rates of S. cerevisiae, we have studied the physiological characteristics and intracellular metabolite profiles of recombinant S. cerevisiae strains with improved or impaired growth on galactose. Aerobic batch cultivations on galactose of strains with different combinations of overexpression of the genes GAL1, GAL2, GAL7, and GAL10, which encode proteins that together convert extracellular galactose into glucose-1-phosphate, revealed a decrease in the maximum specific growth rate when compared to the reference strain. The hypothesized toxic intermediate galactose-1-phosphate cannot be the sole cause of galactose related toxicity, but indications were found that galactose-1-phosphate might cause a negative effect through inhibition of phosphoglucomutase. Furthermore, we show that galactitol is formed in S. cerevisiae, and that the combination of elevated intracellular galactitol concentration, and the ratio between galactose-1-phosphate concentration and phosphoglucomutase activity seems to be important for galactose related toxicity causing decreased growth rates.

Growth-rate regulated genes have profound impact on interpretation of transcriptome profiling in Saccharomyces cerevisiae

BACKGROUND

Growth rate is central to the development of cells in all organisms. However, little is known about the impact of changing growth rates. We used continuous cultures to control growth rate and studied the transcriptional program of the model eukaryote Saccharomyces cerevisiae, with generation times varying between 2 and 35 hours.

RESULTS

A total of 5930 transcripts were identified at the different growth rates studied. Consensus clustering of these revealed that half of all yeast genes are affected by the specific growth rate, and that the changes are similar to those found when cells are exposed to different types of stress (>80% overlap). Genes with decreased transcript levels in response to faster growth are largely of unknown function (>50%) whereas genes with increased transcript levels are involved in macromolecular biosynthesis such as those that encode ribosomal proteins. This group also covers most targets of the transcriptional activator RAP1, which is also known to be involved in replication. A positive correlation between the location of replication origins and the location of growth-regulated genes suggests a role for replication in growth rate regulation.

CONCLUSION

Our data show that the cellular growth rate has great influence on transcriptional regulation. This, in turn, implies that one should be cautious when comparing mutants with different growth rates. Our findings also indicate that much of the regulation is coordinated via the chromosomal location of the affected genes, which may be valuable information for the control of heterologous gene expression in metabolic engineering.

Global transcriptional and physiological responses of Saccharomyces cerevisiae to ammonium, L-alanine, or L-glutamine limitation

The yeast Saccharomyces cerevisiae encounters a range of nitrogen sources at various concentrations in its environment. The impact of these two parameters on transcription and metabolism was studied by growing S. cerevisiae in chemostat cultures with l-glutamine, l-alanine, or l-ammonium in limitation and by growing cells in an excess of ammonium. Cells grown in l-alanine-limited cultures had higher biomass yield per nitrogen mole (19%) than those from ammonium-limited cultures. Whole-genome transcript profiles were analyzed with a genome-scale metabolic model that suggested increased anabolic activity in l-alanine-limited cells. The changes in these cells were found to be focused around pyruvate, acetyl coenzyme A, glyoxylate, and alpha-ketoglutarate via increased levels of ALT1, DAL7, PYC1, GDH2, and ADH5 and decreased levels of GDH3, CIT2, and ACS1 transcripts. The transcript profiles were then clustered. Approximately 1,400 transcripts showed altered levels when amino acid-grown cells were compared to those from ammonium. Another 400 genes had low transcript levels when ammonium was in excess. Overrepresentation of the GATAAG element in their promoters suggests that nitrogen catabolite repression (NCR) may be responsible for this regulation. Ninety-one genes had transcript levels on both l-glutamine and ammonium that were decreased compared to those on l-alanine, independent of the concentration. The GATAAG element in these genes suggests two groups of NCR-responsive genes, those that respond to high levels of nitrogen and those that respond to levels below 30 muM. In conclusion, our results reveal that the nitrogen source has substantial influence on the transcriptome of yeasts and that transcriptional changes may be correlated to physiology via a metabolic model.

In silico aided metabolic engineering of Saccharomyces cerevisiae for improved bioethanol production

In silico genome-scale cell models are promising tools for accelerating the design of cells with improved and desired properties. We demonstrated this by using a genome-scale reconstructed metabolic network of Saccharomyces cerevisiae to score a number of strategies for metabolic engineering of the redox metabolism that will lead to decreased glycerol and increased ethanol yields on glucose under anaerobic conditions. The best-scored strategies were predicted to completely eliminate formation of glycerol and increase ethanol yield with 10%. We successfully pursued one of the best strategies by expressing a non-phosphorylating, NADP(+)-dependent glyceraldehyde-3-phosphate dehydrogenase in S. cerevisiae. The resulting strain had a 40% lower glycerol yield on glucose while the ethanol yield increased with 3% without affecting the maximum specific growth rate. Similarly, expression of GAPN in a strain harbouring xylose reductase and xylitol dehydrogenase led to an improvement in ethanol yield by up to 25% on xylose/glucose mixtures.

Deletion of RTS1, encoding a regulatory subunit of protein phosphatase 2A, results in constitutive amino acid signaling via increased Stp1p processing

In Saccharomyces cerevisiae, extracellular amino acids are sensed at the plasma membrane by the SPS sensor, consisting of the transporter homologue Ssy1p, Ptr3p, and the endoprotease Ssy5p. Amino acid sensing results in proteolytic truncation of the transcription factors Stp1p and Stp2p, followed by their relocation from the cytoplasm to the nucleus, where they activate transcription of amino acid permease genes. We screened a transposon mutant library for constitutively signaling mutants, with the aim of identifying down-regulating components of the SPS-mediated pathway. Three isolated mutants were carrying a transposon in the RTS1 gene, which encodes a regulatory subunit of protein phosphatase 2A. We investigated the basal activity of the AGP1 and BAP2 promoters in rts1delta cells and found increased transcription from these promoters, as well as increased Stp1p processing, even in the absence of amino acids. Based on our findings we propose that the phosphatase complex containing Rts1p keeps the SPS-mediated pathway down-regulated in the absence of extracellular amino acids by dephosphorylating a component of the pathway.

Improvement of galactose uptake in Saccharomyces cerevisiae through overexpression of phosphoglucomutase: example of transcript analysis as a tool in inverse metabolic engineering

Through genome-wide transcript analysis of a reference strain and two recombinant Saccharomyces cerevisiae strains with different rates of galactose uptake, we obtained information about the global transcriptional response to metabolic engineering of the GAL gene regulatory network. One of the recombinant strains overexpressed the gene encoding the transcriptional activator Gal4, and in the other strain the genes encoding Gal80, Gal6, and Mig1, which are negative regulators of the GAL system, were deleted. Even though the galactose uptake rates were significantly different in the three strains, we surprisingly did not find any significant changes in the expression of the genes encoding the enzymes catalyzing the first steps of the pathway (i.e., the genes encoding Gal2, Gal1, Gal7, and Gal10). We did, however, find that PGM2, encoding the major isoenzyme of phosphoglucomutase, was slightly up-regulated in the two recombinant strains with higher galactose uptake rates. This indicated that PGM2 is a target for overexpression in terms of increasing the flux through the Leloir pathway, and through overexpression of PGM2 the galactose uptake rate could be increased by 70% compared to that of the reference strain. Based on our findings, we concluded that phosphoglucomutase plays a key role in controlling the flux through the Leloir pathway, probably due to increased conversion of glucose-1-phosphate to glucose-6-phosphate. This conclusion was supported by measurements of sugar phosphates, which showed that there were increased concentrations of glucose-6-phosphate, galactose-6-phosphate, and fructose-6-phosphate in the strain construct overexpressing PGM2.

Robust multi-scale clustering of large DNA microarray datasets with the consensus algorithm

MOTIVATION

Hierarchical and relocation clustering (e.g. K-means and self-organizing maps) have been successful tools in the display and analysis of whole genome DNA microarray expression data. However, the results of hierarchical clustering are sensitive to outliers, and most relocation methods give results which are dependent on the initialization of the algorithm. Therefore, it is difficult to assess the significance of the results. We have developed a consensus clustering algorithm, where the final result is averaged over multiple clustering runs, giving a robust and reproducible clustering, capable of capturing small signal variations. The algorithm preserves valuable properties of hierarchical clustering, which is useful for visualization and interpretation of the results.

RESULTS

We show for the first time that one can take advantage of multiple clustering runs in DNA microarray analysis by collecting re-occurring clustering patterns in a co-occurrence matrix. The results show that consensus clustering obtained from clustering multiple times with Variational Bayes Mixtures of Gaussians or K-means significantly reduces the classification error rate for a simulated dataset. The method is flexible and it is possible to find consensus clusters from different clustering algorithms. Thus, the algorithm can be used as a framework to test in a quantitative manner the homogeneity of different clustering algorithms. We compare the method with a number of state-of-the-art clustering methods. It is shown that the method is robust and gives low classification error rates for a realistic, simulated dataset. The algorithm is also demonstrated for real datasets. It is shown that more biological meaningful transcriptional patterns can be found without conservative statistical or fold-change exclusion of data.

AVAILABILITY

Matlab source code for the clustering algorithm ClusterLustre, and the simulated dataset for testing are available upon request from T.G. and O.W.

Grr1p is required for transcriptional induction of amino acid permease genes and proper transcriptional regulation of genes in carbon metabolism of Saccharomyces cerevisiae

The F-box protein Grr1p is involved in cell cycle regulation, glucose repression and transcriptional induction of the amino acid permease (AAP) gene AGP1. We investigated the role of Grr1p in amino acid-mediated induction of AAP genes by performing batch cultivations with a wild-type strain and a grr1Delta strain and adding citrulline in the exponential phase. Whole-genome transcription analyses were performed on samples from each cultivation, both immediately before and 30 min after citrulline addition. Transcriptional induction of the AAP genes AGP1, BAP2, BAP3, DIP5, GNP1 and TAT1 is fully dependent on Grr1p. Comparison of the grr1Delta strain with the reference strain in the absence of citrulline revealed that GRR1 disruption leads to increased transcription of numerous genes. These encode enzymes in the tricarboxylic acid cycle, the pentose-phosphate pathway and both glucose and starch metabolism. Promoter analysis showed that many of the genes with increased transcription display Mig1p- and/or Msn2p/Msn4p-binding sites. Increased expression of glucose-repressed genes in the grr1Delta strain may be explained by the reduced expression of the hexose transporter genes HXT1, HXT2, HXT3 and HXT4 and a subsequent lowering of the glucose uptake; and the effect of GRR1 deletion on general carbon metabolism may therefore be indirect. Finally, none of the genes known to be primarily involved in cell cycle regulation displayed different expression levels in the grr1Delta cells as compared with the reference strain, suggesting that the role of Grr1p in cell cycle regulation does not include any transcriptional component.

Use of laminar flow patterning for miniaturised biochemical assays

Laminar flow in microfluidic chambers was used to construct low (one dimensional) density arrays suitable for miniaturized biochemical assays. By varying the ratio of flows of two guiding streams flanking a sample stream, precise focusing and positioning of the latter was achieved, and reactive species carried in the sample stream were deposited on functionalized chip surfaces as discrete 50 microm wide lanes. Using different model systems we have confirmed the method's suitability for qualitative screening and quantification tasks in receptor-ligand assays, recording biotin-streptavidin interactions, DNA-hybridization and DNA-triplex formation. The system is simple, fast, reproducible, flexible, and has small sample requirements.

Transcriptional profiling of extracellular amino acid sensing in Saccharomyces cerevisiae and the role of Stp1p and Stp2p

S. cerevisiae responds to the presence of amino acids in the environment through the membrane-bound complex SPS, by altering transcription of several genes. Global transcription analysis shows that 46 genes are induced by L-citrulline. Under the given conditions there appears to be only one pathway for induction with L-citrulline, and this pathway is completely dependent on the SPS component, Ssy1p, and either of the transcription factors, Stp1p and Stp2p. Besides the effects on amino acid permease genes, an ssy1 and an stp1 stp2 mutant exhibit a number of other transcriptional phenotypes, such as increased expression of genes subject to nitrogen catabolite repression and genes involved in stress response. A group of genes involved in the upper part of the glycolysis, including those encoding hexose transporters Hxt4p, Hxt5p, Hxt6p, Hxt7p, hexokinase Hxk1p, glyceraldehyde 3-phosphate dehydrogenase Tdh1p and glucokinase (Glk1p), shows increased transcription levels in either or both of the mutants. Also, most of the structural genes involved in trehalose and glycogen synthesis and a few genes in the glyoxylate cycle and the pentose phosphate pathway are derepressed in the ssy1 and stp1 stp2 strains.

Genome-wide transcriptional response of a Saccharomyces cerevisiae strain with an altered redox metabolism

The genome-wide transcriptional response of a Saccharomyces cerevisiae strain deleted in GDH1 that encodes a NADP(+)-dependent glutamate dehydrogenase was compared to a wild-type strain under anaerobic steady-state conditions. The GDH1-deleted strain has a significantly reduced NADPH requirement, and therefore, an altered redox metabolism. Identification of genes with significantly changed expression using a t-test and a Bonferroni correction yielded only 16 transcripts when accepting two false-positives, and 7 of these were Open Reading Frames (ORFs) with unknown function. Among the 16 transcripts the only one with a direct link to redox metabolism was GND1, encoding phosphogluconate dehydrogenase. To extract additional information we analyzed the transcription data for a gene subset consisting of all known genes encoding metabolic enzymes that use NAD(+) or NADP(+). The subset was analyzed for genes with significantly changed expression again with a t-test and correction for multiple testing. This approach was found to enrich the analysis since GND1, ZWF1 and ALD6, encoding the most important enzymes for regeneration of NADPH under anaerobic conditions, were down-regulated together with eight other genes encoding NADP(H)-dependent enzymes. This indicates a possible common redox-dependent regulation of these genes. Furthermore, we showed that it might be necessary to analyze the expression of a subset of genes to extract all available information from global transcription analysis.

Transcriptional, proteomic, and metabolic responses to lithium in galactose-grown yeast cells

Lithium is highly toxic to yeast when grown in galactose medium mainly because phosphoglucomutase, a key enzyme of galactose metabolism, is inhibited. We studied the global protein and gene expression profiles of Saccharomyces cerevisiae grown in galactose in different time intervals after addition of lithium. These results were related to physiological studies where both secreted and intracellular metabolites were determined. Microarray analysis showed that 664 open reading frames were down-regulated and 725 up-regulated in response to addition of lithium. Genes involved in transcription, translation, and nucleotide metabolism were down-regulated at the transcriptional level, whereas genes responsive to different stresses as well as genes from energy reserve metabolism and monosaccharide metabolism were up-regulated. Compared with the proteomic data, 26% of the down-regulated and 48% of the up-regulated proteins were also identified as being changed on the mRNA level. Functional clusters obtained from proteome data were coincident with transcriptional clusters. Physiological studies showed that acetate, glycerol, and glycogen accumulate in response to lithium, as reflected in expression data, whereas a change from respiro-fermentative to respiratory growth could not be predicted from the expression analyses.

Reproducibility of oligonucleotide microarray transcriptome analyses. An interlaboratory comparison using chemostat cultures of Saccharomyces cerevisiae

Assessment of reproducibility of DNA-microarray analysis from published data sets is complicated by the use of different microbial strains, cultivation techniques, and analytical procedures. Because intra- and interlaboratory reproducibility is highly relevant for application of DNA-microarray analysis in functional genomics and metabolic engineering, we designed a set of experiments to specifically address this issue. Saccharomyces cerevisiae CEN.PK113-7D was grown under defined conditions in glucose-limited chemostats, followed by transcriptome analysis with Affymetrix GeneChip arrays. In each of the laboratories, three independent replicate cultures were grown aerobically as well as anaerobically. Although variations introduced by in vitro handling steps were small and unbiased, greater variation from replicate cultures underscored that, to obtain reliable information, experimental replication is essential. Under aerobic conditions, 86% of the most highly expressed yeast genes showed an average intralaboratory coefficient of variation of 0.23. This is significantly lower than previously reported for shake-flask-culture transcriptome analyses and probably reflects the strict control of growth conditions in chemostats. Using the triplicate data sets and appropriate statistical analysis, the change calls from anaerobic versus aerobic comparisons yielded an over 95% agreement between the laboratories for transcripts that changed by over 2-fold, leaving only a small fraction of genes that exhibited laboratory bias.

Amino acid residues important for substrate specificity of the amino acid permeases Can1p and Gnp1p in Saccharomyces cerevisiae

Deletion of the general amino acid permease gene GAP1 abolishes uptake of L-citrulline in Saccharomyces cerevisiae, resulting in the inability to grow on L-citrulline as sole nitrogen source. Selection for suppressor mutants that restored growth on L-citrulline led to isolation of 21 mutations in the arginine permease gene CAN1. One similar mutation was found in the glutamine-asparagine permease gene GNP1. L-[(14)C]citrulline uptake measurements confirmed that suppressor mutations in CAN1 conferred uptake of this amino acid, while none of the mutant permeases had lost the ability to transport L-[(14)C]arginine. Substrate specificity seemed to remain narrow in most cases, and broad substrate specificity was only observed in the cases where mutations affect two proline residues (P148 and P313) that are both conserved in the amino acid-polyamine-choline (APC) transporter superfamily. We found mutations affecting six predicted domains (helices III and X, and loops 1, 2, 6 and 7) of the permeases. Helix III and loop 7 are candidates for domains in direct contact with thetransported amino acid. Helix III was affected in both CAN1 (Y173H, Y173D) and GNP1 (W239C) mutants and has previously been found to be important for substrate preference in other members of the family. Furthermore, the mutations affecting loop 7 (residue T354, S355, Y356) are close to a glutamate side chain (E367) potentially interacting with the positively charged substrate, a notion supported by conservation of the side chain in permeases for cationic substrates.

GAP1, a novel selection and counter-selection marker for multiple gene disruptions in Saccharomyces cerevisiae

We report on the use of a new homologous marker for use in multiple gene deletions in S. cerevisiae, the general amino acid permease gene (GAP1). A GAP1 strain can utilize L-citrulline as the sole nitrogen source but cannot grow in the presence of the toxic amino acid D-histidine. L-citrulline as well as D-histidine uptake is mediated solely by the general amino acid permease, and a gap1 strain is therefore able to grow in the presence of D-histidine but cannot utilize L-citrulline. Gene disruption is effected by transforming a gap1 strain with a gene cassette generated by PCR, containing GAP1 flanked by short (60 bp) stretches of the gene in question. Through homologous recombination, the cassette will integrate into the target gene, which is thus replaced by GAP1, and mutants are selected for on minimal L-citrulline medium. When propagated under non-selective conditions, some cells will lose the GAP1 gene. This is caused by recombination between two Ashbya gossypii AgLEU2 [corrected] direct repeats embracing GAP1, and will result in a sub-population of gap1 cells. Such cells are selected on a medium containing D-histidine, and may subsequently be used for a second gene disruption. Hence, multiple gene disruptions can be made fast, cheaply and easily in a gap1 strain, with two positive selection steps for each disruption.

Substrate specificity and gene expression of the amino-acid permeases in Saccharomyces cerevisiae

All known amino-acid permeases (AAPs) in Saccharomyces cerevisiae belong to a single family of homologous proteins. Genes of 15 AAPs were overexpressed in different strains, and the ability to take up one or more of the 20 common L-alpha-amino acids was studied in order to obtain a complete picture of the substrate specificity for these permeases. Radiolabelled amino-acid uptake measurements showed that Agp1p is a general permease for most uncharged amino acids (Ala, Gly, Ser, Thr, Cys, Met, Phe, Tyr, Ile, Leu, Val, Gln and Asn). Gnp1p, which is closely related to Agp1p, has a somewhat less-broad specificity, transporting Leu, Ser, Thr, Cys, Met, Gln and Asn, while Bap2p and Bap3p, which are also closely related to Agp1p, are able to transport Ile, Leu, Val, Cys, Met, Phe, Tyr and Trp. All four permeases are transcriptionally induced by an extracellular amino acid, but differ in expression with respect to the nitrogen source. On a non-repressive nitrogen source, AGP1 is induced, while GLN1, BAP2 and BAP3 are not. Except for Dip5p, which is a transporter for Glu, Asp, Gln, Asn, Ser, Ala and Gly, the rest of the permeases exhibit narrow specificity. Tat2p can take up Phe, Trp and Tyr; Put4p can transport Ala, Gly and Pro; while Can1p, Lyp1p and the previously uncharacterized Alp1p are specific for the cationic amino acids. These findings modify the prevalent view that S. cerevisiae only contains one general amino-acid permease, Gap1p, and a number of permeases that are specific for a single or a few amino acids.

Cysteine uptake by Saccharomyces cerevisiae is accomplished by multiple permeases

Uptake by Saccharomyces cerevisiae of the sulphur-containing amino acid L-cysteine was found to be non-saturable under various conditions, and uptake kinetics suggested the existence of two or more transport systems in addition to the general amino-acid permease, Gap1p. Overexpression studies identified BAP2, BAP3, AGP1 and GNP1 as genes encoding transporters of cysteine. Uptake studies with disruption mutants confirmed this, and identified two additional genes for transporters of cysteine, TAT1 and TAT2, both very homologous to BAP2, BAP3, AGP1 and GNP1. While Gap1p and Agp1p appear to be the main cysteine transporters on the non-repressing nitrogen source proline, Bap2p, Bap3p, Tat1p, Tat2p, Agp1p and Gnp1p are all important for cysteine uptake on ammonium-based medium. Furthermore, whereas Bap2p, Bap3p, Tat1p and Tat2p seem most important under amino acid-rich conditions, Agp1p contributes significantly when only ammonium is present, and Gnp1p only contributes under the latter condition.