Genome-wide analysis of nucleotide-level variation in commonly used Saccharomyces cerevisiae strains

Ten years have passed since the genome of Saccharomyces cerevisiae-more precisely, the S288c strain-was completely sequenced. However, experimental work in yeast is commonly performed using strains that are of unknown genetic relationship to S288c. Here, we characterized the nucleotide-level similarity between S288c and seven commonly used lab strains (A364A, W303, FL100, CEN.PK, summation 1278b, SK1 and BY4716) using 25mer oligonucleotide microarrays that provide complete and redundant coverage of the approximately 12 Mb Saccharomyces cerevisiae genome. Using these data, we assessed the frequency and distribution of nucleotide variation in comparison to the sequenced reference genome. These data allow us to infer the relationships between experimentally important strains of yeast and provide insight for experimental designs that are sensitive to sequence variation. We propose a rational approach for near complete sequencing of strains related to the reference using these data and directed re-sequencing. These data and new visualization tools are accessible online in a new resource: the Yeast SNPs Browser (YSB; http://gbrowse.princeton.edu/cgi-bin/gbrowse/yeast_strains_snps) that is available to all researchers.     

Structure of a conjugative element in Streptococcus pneumoniae.

We have cloned and mapped a 69-kilobase (kb) region of the chromosome of Streptococcus pneumoniae DP1322, which carries the conjugative omega (cat-tet) insertion from S. pneumoniae BM6001. This element proved to be 65.5 kb in size. Location of the junctions was facilitated by cloning a preferred target region from the wild-type strain Rx1 recipient genome. This target site was preferred by both the BM6001 element and the cat-erm-tet element from Streptococcus agalactiae B109. Within the BM6001 element cat and tet were separated by 30 kb, and cat was flanked by two copies of a sequence that was also present in the recipient strain Rx1 DNA. Another sequence at least 2.4 kb in size was found inside the BM6001 element and at two places in the Rx1 genome. Its role is unknown. The ends of the BM6001 element appear to be the same as those of the B109 element, both as seen after transfer to S. pneumoniae and as mapped by others in pDP5 after transposition in Streptococcus faecalis. We see no homology between the ends of the BM6001 element and find no evidence suggesting that it ever circularizes.     

Evolution of plant mitochondrial genomes via substoichiometric intermediates

Comparison of the modern fertile maize mitochondrial genome (N) with an ancestral maize mitochondrial genome (RU) reveals a 12 kb duplication (containing the atpA gene) in the modern genome that is absent from the ancestor. Cloning, mapping, and sequencing of the relevant portions of the ancestral genome shows that this duplication probably arose via a three-stage recombination process involving substoichiometric intermediates. Comparison with analogous observations on yeast mitochondrial genomes suggests that this three-stage model of genome reorganization can be generally applied to plant mitochondrial genomes to explain both deletions and the creation of novel repeats, common features of plant mitochondrial genome evolution.     

Paired-End Sequencing of Long-Range DNA Fragments for De Novo Assembly of Large,Complex Mammalian Genomes by Direct Intra-Molecule Ligation

Background

The relatively short read lengths from next generation sequencing (NGS) technologies still pose a challenge for de novo assembly of complex mammal genomes. One important solution is to use paired-end (PE) sequence information experimentally obtained from long-range DNA fragments (>1 kb). Here, we characterize and extend a long-range PE library construction method based on direct intra-molecule ligation (or molecular linker-free circularization) for NGS.

Results

We found that the method performs stably for PE sequencing of 2- to 5- kb DNA fragments, and can be extended to 10–20 kb (and even in extremes, up to ∼35 kb). We also characterized the impact of low quality input DNA on the method, and develop a whole-genome amplification (WGA) based protocol using limited input DNA (<1 µg). Using this PE dataset, we accurately assembled the YanHuang (YH) genome, the first sequenced Asian genome, into a scaffold N50 size of >2 Mb, which is over100-times greater than the initial size produced with only small insert PE reads(17 kb). In addition, we mapped two 7- to 8- kb insertions in the YH genome using the larger insert sizes of the long-range PE data.

Conclusions

In conclusion, we demonstrate here the effectiveness of this long-range PE sequencing method and its use for the de novo assembly of a large, complex genome using NGS short reads.     

金针菇基因组测序与组装策略分析

采用二代和三代测序技术分别对金针菇单核体菌株“6-3”进行测序,应用4种组装策略进行基因组的de novo组装,对比组装效果。基因组组装的参数方面,仅使用二代测序组装的效果最差,长度大于10kb的Contig全长只有24.6Mb,Contig N50只有23kb,组装率只有59.27%。采用三代组装二代校正的组装策略效果最好,长度大于10kb的Contig全长为38.3Mb,Contig N50为2.8Mb,组装率高达92.16%。保守单拷贝基因拼接效果方面,4种组装策略获得基因组序列与BUSCO数据库里的担子菌的保守单拷贝基因比对,基因完整性均大于94%。在组装准确性方面,经过PCR扩增、Sanger测序验证,三代组装二代校正的基因组序列完整并且连续,同时序列上碱基的SNP、InDel数量最少。综上所述,三代组装二代校正得到的基因组序列具有Contig N50值大、组装率高、碱基准确性高的特点,是食用菌基因组测序较为理想的方案。     

Direct sequencing and characterization of a clinical isolate of Epstein-Barr virus from nasopharyngeal carcinoma tissue by using next-generation sequencing technology

Epstein-Barr virus (EBV)-encoded molecules have been detected in the tumor tissues of several cancers, including nasopharyngeal carcinoma (NPC), suggesting that EBV plays an important role in tumorigenesis. However, the nature of EBV with respect to genome width in vivo and whether EBV undergoes clonal expansion in the tumor tissues are still poorly understood. In this study, next-generation sequencing (NGS) was used to sequence DNA extracted directly from the tumor tissue of a patient with NPC. Apart from the human sequences, a clinically isolated EBV genome 164.7 kb in size was successfully assembled and named GD2 (GenBank accession number HQ020558). Sequence and phylogenetic analyses showed that GD2 was closely related to GD1, a previously assembled variant derived from a patient with NPC. GD2 contains the most prevalent EBV variants reported in Cantonese patients with NPC, suggesting that it might be the prevalent strain in this population. Furthermore, GD2 could be grouped into a single subtype according to common classification criteria and contains only 6 heterozygous point mutations, suggesting the monoclonal expansion of GD2 in NPC. This study represents the first genome-wide analysis of a clinical isolate of EBV directly extracted from NPC tissue. Our study reveals that NGS allows the characterization of genome-wide variations of EBV in clinical tumors and provides evidence of monoclonal expansion of EBV in vivo. The pipeline could also be applied to the study of other pathogen-related malignancies. With additional NGS studies of NPC, it might be possible to uncover the potential causative EBV variant involved in NPC.     

Binding of the glucose-dependent Mig1p repressor to the GAL1 and GAL4 promoters in vivo: regulationby glucose and chromatin structure.

Binding of the MIG1 repressor to the glucose-repressible GAL1 and GAL4 promoters was analyzed in vivo by cyclobutane dimer footprinting in two yeast strains that show different glucose repression responses. Mig1p binding to the two promoters in both strains was glucose-induced. In cells subject to rapid and stringent glucose repression (S288c), long-term Mig1p binding in glucose-grown cells was inhibited by the formation of a competing chromatin structure. Under conditions where glucose repression was only partially effective (gal80 - or low glucose), the chromatin structure did not form and long-term Mig1p binding was observed The same long-term binding of Mig1p was seen in cells of a different strain (W303A) that shows only partial glucose repression of the GAL1 promoter. We conclude from these experiments that Mig1p binding to glucose-repressed promoters is glucose-dependent but transient. We suggest that Mig1p functions at an early step in repression, but is not required to maintain the repressed state.     

Comparative genome analysis of a Saccharomyces cerevisiae wine strain

Many industrial strains of Saccharomyces cerevisiae have been selected primarily for their ability to convert sugars into ethanol efficiently despite exposure to a variety of stresses. To begin investigation of the genetic basis of phenotypic variation in industrial strains of S. cerevisiae, we have sequenced the genome of a wine yeast, AWRI1631, and have compared this sequence with both the laboratory strain S288c and the human pathogenic isolate YJM789. AWRI1631 was found to be substantially different from S288c and YJM789, especially at the level of single-nucleotide polymorphisms, which were present, on average, every 150 bp between all three strains. In addition, there were major differences in the arrangement and number of Ty elements between the strains, as well as several regions of DNA that were specific to AWRI1631 and that were predicted to encode proteins that are unique to this industrial strain.     

Crosstalk of prolyl isomerases, Pin1/Ess1, and cyclophilin A.

Previous studies have indicated that Ess1/Pin1, a gene in the parvulin family of peptidyl-prolyl isomerases (PPIases), plays an important role in regulating the G(2)/M transition of the cell cycle by binding cell-cycle-regulating proteins in eukaryotic cells. Although the ess1 gene has been considered to be essential in yeast, we have isolated viable ess1 deletion mutants and demonstrated, via analysis of yeast gene expression profiles using microarray techniques, a novel regulatory role for ESS1 in the G(1) phase. Although the overall expression profiles in the tested strains (C110-1, W303, S288c, and RAY-3AD) were similar, marked changes were detected for a number of genes involved in the molecular action of ESS1. Among these, the expression levels of a cyclophilin A gene, also a member of the PPIase family, increased in the ess1 null mutant derived from C110-1. Subsequent treatment with cyclosporin A significantly retarded growth, which suggests that ESS1 and cyclophilin A are functionally linked in yeast cells and play important roles at the G(1) phase of the cell cycle.     

Comparative genotyping of the Saccharomyces cerevisiae laboratory strains S288C and CEN.PK113-7D using oligonucleotide microarrays

To analyse the reliability and accuracy of genotype analysis with high-density oligonucleotide microarrays, this method and other experimental approaches were used to analyse genomic DNA of two popular Saccharomyces cerevisiae laboratory strains. S288C was used for systematic sequencing of 'the' S. cerevisiae genome; CEN.PK113-7D is a popular strain for physiological studies and functional genomics. Random amplified polymorphic DNA, electrophoretic karyotyping and microarray analysis all indicated a high level of sequence similarity between the two strains. In the microarray analysis, as few as 288 (4.5%) of the ca. 6300 represented yeast genes were identified that yielded significantly different hybridisation intensities between the two strains. These could be classified as amplified, absent, or with sequence polymorphism in CEN.PK113-7D compared to S288C. A detailed analysis focused on the subset of 25 genes called absent in CEN.PK113-7D. Among these absent genes, 17 were clustered together on five chromosomes, mainly in subtelomeric regions. Thorough analysis of these regions by polymerase chain reaction (PCR) and restriction fragment length polymorphism confirmed the absence of these genes in CEN.PK113-7D. Surprisingly, three of these regions were not smaller in CEN.PK113-7D chromosomes, indicating that they may harbour unidentified and potentially new sequences. In addition, eight genes called absent by the microarrays were scattered over the chromosomes. Using diagnostic PCR most of these genes were actually found to be present in CEN.PK113-7D, but after sequencing were found to differ significantly at the DNA level from S288C, explaining the poor hybridisation to the arrays. Our results indicate that DNA microarrays are a powerful tool for determining genotypic similarity between different yeast strains. However, to obtain meaningful information at the individual gene level, this method should be backed up by additional techniques.     

采用诱导环化酶制备Illumina Mate-paired DNA文库

新一代测序技术(NGS)的文库制备方法在基因组的拼装中起着重要作用。但是NGS技术制备的普通DNA文库片段只有500 bp左右,难以满足复杂基因组的从头(de novo)拼装要求。三代测序技术的读长可以达到20 kb,但是其高错误率及测序成本过高使得其又不易推广。因此二代测序的Mate-paired文库制备技术一直在基因组的de novo拼装中扮演着非常重要的角色。目前主流的NGS平台Illumina制备的Mate-paired文库的片段范围只有2~5 kb,为了得到更长的可用于Illumina平台测序的Mate-paired文库,本研究首次整合并优化了Illumina和Roche/454两种测序平台的Mate-paired文库制备技术,采用诱导环化酶来提高基因组长片段DNA的环化效率,成功建立了20 kb Mate-paired文库制备技术,并已将该技术应用于人类基因组20 kb Mate-paired文库制备。该技术为Illumina平台制备长片段Mate-paired库提供了方法指导。     

Presence of 16S rRNA genes in multiple replicons in Azospirillum brasilense

Rhizosphere and endophytic Azospirillum brasilense isolates recovered from sugarcane plants and the reference strains Sp7 and Cd were analyzed for plasmid occurrence. All of the 26 A. brasilense isolates analyzed harbored from five to eight replicons. Several strains contained small plasmids from 45 to 70 kb, but all of the isolates harbored other plasmids ranging from 100 to 290 kb and two megareplicons of approximately 1700 and over 1800 kb. Most of the strains contained a replicon with a size of either 570 or 630 kb, and another large 910- or 980-kb replicon. The 1700-kb megareplicon and some others around 600 kb strongly hybridized to 16S rDNA genes, while the 910- or 980-kb replicons hybridized only slightly. This suggests that the A. brasilense genome is composed of multiple minichromosomes instead of a single circular chromosome. The apparent genome complexity of A. brasilense deserves to eventually be resolved by complete genome sequencing.     

Synthetic genome engineering forging new frontiers for wine yeast

Over the past 15 years, the seismic shifts caused by the convergence of biomolecular, chemical, physical, mathematical, and computational sciences alongside cutting-edge developments in information technology and engineering have erupted into a new field of scientific endeavor dubbed Synthetic Biology. Recent rapid advances in high-throughput DNA sequencing and DNA synthesis techniques are enabling the design and construction of new biological parts (genes), devices (gene networks) and modules (biosynthetic pathways), and the redesign of biological systems (cells and organisms) for useful purposes. In 2014, the budding yeast Saccharomyces cerevisiae became the first eukaryotic cell to be equipped with a fully functional synthetic chromosome. This was achieved following the synthesis of the first viral (poliovirus in 2002 and bacteriophage Phi-X174 in 2003) and bacterial (Mycoplasma genitalium in 2008 and Mycoplasma mycoides in 2010) genomes, and less than two decades after revealing the full genome sequence of a laboratory (S288c in 1996) and wine (AWRI1631 in 2008) yeast strain. A large international project – the Synthetic Yeast Genome (Sc2.0) Project – is now underway to synthesize all 16 chromosomes (～12?Mb carrying ～6000 genes) of the sequenced S288c laboratory strain by 2018. If successful, S. cerevisiae will become the first eukaryote to cross the horizon of in silico design of complex cells through de novo synthesis, reshuffling, and editing of genomes. In the meantime, yeasts are being used as cell factories for the semi-synthetic production of high-value compounds, such as the potent antimalarial artemisinin, and food ingredients, such as resveratrol, vanillin, stevia, nootkatone, and saffron. As a continuum of previously genetically engineered industrially important yeast strains, precision genome engineering is bound to also impact the study and development of wine yeast strains supercharged with synthetic DNA. The first taste of what the future holds is the de novo production of the raspberry ketone aroma compound, 4-[4-hydroxyphenyl]butan-2-one, in a wine yeast strain (AWRI1631), which was recently achieved via metabolic pathway engineering and synthetic enzyme fusion. A peek over the horizon is revealing that the future of “Wine Yeast 2.0” is already here. Therefore, this article seeks to help prepare the wine industry – an industry rich in history and tradition on the one hand, and innovation on the other – for the inevitable intersection of the ancient art practiced by winemakers and the inventive science of pioneering “synthetic genomicists”. It would be prudent to proactively engage all stakeholders – researchers, industry practitioners, policymakers, regulators, commentators, and consumers – in a meaningful dialog about the potential challenges and opportunities emanating from Synthetic Biology. To capitalize on the new vistas of synthetic yeast genomics, this paper presents wine yeast research in a fresh context, raises important questions and proposes new directions.     

An evolutionary approach uncovers a diverse response of tRNA 2-thiolation to elevated temperatures in yeast

Chemical modifications of transfer RNA (tRNA) molecules are evolutionarily well conserved and critical for translation and tRNA structure. Little is known how these nucleoside modifications respond to physiological stress. Using mass spectrometry and complementary methods, we defined tRNA modification levels in six yeast species in response to elevated temperatures. We show that 2-thiolation of uridine at position 34 (s²U₃₄) is impaired at temperatures exceeding 30°C in the commonly used Saccharomyces cerevisiae laboratory strains S288C and W303, and in Saccharomyces bayanus. Upon stress relief, thiolation levels recover and we find no evidence that modified tRNA or s²U₃₄ nucleosides are actively removed. Our results suggest that loss of 2-thiolation follows accumulation of newly synthesized tRNA that lack s²U₃₄ modification due to temperature sensitivity of the URM1 pathway in S. cerevisiae and S. bayanus. Furthermore, our analysis of the tRNA modification pattern in selected yeast species revealed two alternative phenotypes. Most strains moderately increase their tRNA modification levels in response to heat, possibly constituting a common adaptation to high temperatures. However, an overall reduction of nucleoside modifications was observed exclusively in S288C. This surprising finding emphasizes the importance of studies that utilize the power of evolutionary biology, and highlights the need for future systematic studies on tRNA modifications in additional model organisms.     

Both yeast W double-stranded RNA and its single-stranded form 20S RNA are linear.

Most yeast strains carry a cytoplasmic double-stranded RNA (dsRNA) molecule called W, of 2.5 kb in size. We have cloned and sequenced most of W genome (1), and we proposed that W (+) strands were identical to 20S RNA, a single-stranded RNA (ssRNA) species, whose copy number is highly induced under stress conditions. Recently it was proposed that 20S RNA was circular (2). In this paper, however, we demonstrate that both W dsRNA and 20S RNA are linear. Linearity of W dsRNA is shown by the stoichiometric labelling of both strands of W with 32P-pCp and T4 RNA ligase. The last 3' end nucleotide of both strands is about 70 to 80% C and 20 to 30% A. Linearity of 20S RNA is directly demonstrated by a site-specific cleavage of 20S RNA with RNase H, using an oligodeoxynucleotide complementary to an internal site of 20S RNA. The cleavage produced not one but two RNA fragments expected from the linearity of 20S RNA.     

Walking, cloning, and mapping with yeast artificial chromosomes: a contig encompassing D21S13 and D21S16.

Chromosome 21 has often been used as a model system for the development of genome mapping and cloning strategies in humans. In this report methods for systematic chromosome walking, cloning, and mapping are exemplified in the construction of a 1.5-Mb yeast artificial chromosome (YAC) contig encompassing and extending 400 kb beyond each of the genetic loci D21S13 and D21S16. Isolation of insert-terminal sequences from YACs in this contig provides a set of closely spaced physical markers. These have been used to generate a long-range genomic restriction map.     

Differential regulation and substrate preferences in two peptide transporters of Saccharomyces cerevisiae

Dal5p has been shown previously to act as an allantoate/ureidosuccinate permease and to play a role in the utilization of certain dipeptides as a nitrogen source in Saccharomyces cerevisiae. Here, we provide direct evidence that dipeptides are transported by Dal5p, although the affinity of Dal5p for allantoate and ureidosuccinate is higher than that for dipeptides. Allantoate, ureidosuccinate, and to a lesser extent allantoin competed with dipeptide transport by reducing the toxicity of the peptide Ala-Eth and decreasing the accumulation of [(14)C]Gly-Leu. In contrast to the well-studied di/tripeptide transporter Ptr2p, whose substrate specificity is very broad, Dal5p preferred to transport non-N-end rule dipeptides. S. cerevisiae W303 was sensitive to the toxic peptide Ala-Eth (non-N-end rule peptide) but not Leu-Eth (N-end rule peptide). Non-N-end rule dipeptides showed better competition with the uptake of [(14)C]Gly-Leu than N-end rule dipeptides. Similar to the regulation of PTR2, DAL5 expression was influenced by the addition of Leu and by the CUP9 gene. However, DAL5 expression was downregulated in the presence of leucine and the absence of CUP9, whereas PTR2 was upregulated. Toxic dipeptide and uptake assays indicated that either Ptr2p or Dal5p was predominantly used for dipeptide transport in the common laboratory strains S288c and W303, respectively. These studies highlight the complementary activities of two dipeptide transport systems under different regulatory controls in common laboratory yeast strains, suggesting that dipeptide transport pathways evolved to respond to different environmental conditions.     

Whole-genome sequencing of the efficient industrial fuel-ethanol fermentative Saccharomyces cerevisiae strain CAT-1

The Saccharomyces cerevisiae strains widely used for industrial fuel-ethanol production have been developed by selection, but their underlying beneficial genetic polymorphisms remain unknown. Here, we report the draft whole-genome sequence of the S. cerevisiae strain CAT-1, which is a dominant fuel-ethanol fermentative strain from the sugarcane industry in Brazil. Our results indicate that strain CAT-1 is a highly heterozygous diploid yeast strain, and the ~12-Mb genome of CAT-1, when compared with the reference S228c genome, contains ~36,000 homozygous and ~30,000 heterozygous single nucleotide polymorphisms, exhibiting an uneven distribution among chromosomes due to large genomic regions of loss of heterozygosity (LOH). In total, 58 % of the 6,652 predicted protein-coding genes of the CAT-1 genome constitute different alleles when compared with the genes present in the reference S288c genome. The CAT-1 genome contains a reduced number of transposable elements, as well as several gene deletions and duplications, especially at telomeric regions, some correlated with several of the physiological characteristics of this industrial fuel-ethanol strain. Phylogenetic analyses revealed that some genes were likely associated with traits important for bioethanol production. Identifying and characterizing the allelic variations controlling traits relevant to industrial fermentation should provide the basis for a forward genetics approach for developing better fermenting yeast strains.