Genome-wide amplifications caused by chromosomal rearrangements play a major role in the adaptive evolution of natural yeast

The relative importance of gross chromosomal rearrangements to adaptive evolution has not been precisely defined. The Saccharomyces cerevisiae flor yeast strains offer significant advantages for the study of molecular evolution since they have recently evolved to a high degree of specialization in a very restrictive environment. Using DNA microarray technology, we have compared the genomes of two prominent variants of S. cerevisiae flor yeast strains. The strains differ from one another in the DNA copy number of 116 genomic regions that comprise 38% of the genome. In most cases, these regions are amplicons flanked by repeated sequences or other recombination hotspots previously described as regions where double-strand breaks occur. The presence of genes that confer specific characteristics to the flor yeast within the amplicons supports the role of chromosomal rearrangements as a major mechanism of adaptive evolution in S. cerevisiae. We propose that nonallelic interactions are enhanced by ethanol- and acetaldehyde-induced double-strand breaks in the chromosomal DNA, which are repaired by pathways that yield gross chromosomal rearrangements. This mechanism of chromosomal evolution could also account for the sexual isolation shown among the flor yeast.     

Clustering Deviation Index (CDI): a robust and accurate internal measure for evaluating scRNA-seq data clustering

Single-cell high-throughput chromatin conformation capture methodologies (scHi-C) enable profiling of long-range genomic interactions. However, data from these technologies are prone to technical noise and biases that hinder downstream analysis. We develop a normalization approach, BandNorm, and a deep generative modeling framework, scVI-3D, to account for scHi-C specific biases. In benchmarking experiments, BandNorm yields leading performances in a time and memory efficient manner for cell-type separation, identification of interacting loci, and recovery of cell-type relationships, while scVI-3D exhibits advantages for rare cell types and under high sparsity scenarios. Application of BandNorm coupled with gene-associating domain analysis reveals scRNA-seq validated sub-cell type identification.     

Systematic bias in high-throughput sequencing data and its correction by BEADS

Genomic sequences obtained through high-throughput sequencing are not uniformly distributed across the genome. For example, sequencing data of total genomic DNA show significant, yet unexpected enrichments on promoters and exons. This systematic bias is a particular problem for techniques such as chromatin immunoprecipitation, where the signal for a target factor is plotted across genomic features. We have focused on data obtained from Illumina's Genome Analyser platform, where at least three factors contribute to sequence bias: GC content, mappability of sequencing reads, and regional biases that might be generated by local structure. We show that relying on input control as a normalizer is not generally appropriate due to sample to sample variation in bias. To correct sequence bias, we present BEADS (bias elimination algorithm for deep sequencing), a simple three-step normalization scheme that successfully unmasks real binding patterns in ChIP-seq data. We suggest that this procedure be done routinely prior to data interpretation and downstream analyses.     

Genome-wide prediction of genetic interactions in a metazoan

Genetic interactions provide information about genes and processes with overlapping functions in biological systems. For Saccharomyces cerevisiae, computational integration of multiple types of functional genomic data is used to generate genome-wide predictions of genetic interactions. However, this methodology cannot be applied to the vastly more complex genome of metazoans, and only recently has the first metazoan genome-wide prediction of genetic interactions been reported. The prediction for Caenorhabditis elegans was generated by computationally integrating functional genomic data from S. cerevisiae, C. elegans and Drosophila melanogaster. This achievement is an important step toward system-level understanding of biological systems and human diseases.     

Group normalization for genomic data

Data normalization is a crucial preliminary step in analyzing genomic datasets. The goal of normalization is to remove global variation to make readings across different experiments comparable. In addition, most genomic loci have non-uniform sensitivity to any given assay because of variation in local sequence properties. In microarray experiments, this non-uniform sensitivity is due to different DNA hybridization and cross-hybridization efficiencies, known as the probe effect. In this paper we introduce a new scheme, called Group Normalization (GN), to remove both global and local biases in one integrated step, whereby we determine the normalized probe signal by finding a set of reference probes with similar responses. Compared to conventional normalization methods such as Quantile normalization and physically motivated probe effect models, our proposed method is general in the sense that it does not require the assumption that the underlying signal distribution be identical for the treatment and control, and is flexible enough to correct for nonlinear and higher order probe effects. The Group Normalization algorithm is computationally efficient and easy to implement. We also describe a variant of the Group Normalization algorithm, called Cross Normalization, which efficiently amplifies biologically relevant differences between any two genomic datasets.     

Selective isolation of genomic loci from complex genomes by transformation-associated recombination cloning in the yeast Saccharomyces cerevisiae

Here, we describe a protocol for the selective isolation of any genomic fragment or gene of interest up to 250 kb in size from complex genomes as a circular yeast artificial chromosome (YAC). The method is based on transformation-associated recombination (TAR) in the yeast Saccharomyces cerevisiae between genomic DNA and a linearized TAR cloning vector containing targeting sequences homologous to a region of interest. Recombination between the vector and homologous sequences in the co-transformed mammalian DNA results in the establishment of a YAC that is able to propagate, segregate and be selected for in yeast. Yield of gene-positive clones varies from 1% to 5%. The entire procedure takes 2 weeks to complete once the TAR vector is constructed and genomic DNA is prepared. The TAR cloning method has a broad application in functional and comparative genomics, long-range haplotyping and characterization of chromosomal rearrangements, including copy number variations.     

Identifying multi-locus chromatin contacts in human cells using tethered multiple 3C

Background

Several recently developed experimental methods, each an extension of the chromatin conformation capture (3C) assay, have enabled the genome-wide profiling of chromatin contacts between pairs of genomic loci in 3D. Especially in complex eukaryotes, data generated by these methods, coupled with other genome-wide datasets, demonstrated that non-random chromatin folding correlates strongly with cellular processes such as gene expression and DNA replication.

Results

We describe a genome architecture assay, tethered multiple 3C (TM3C), that maps genome-wide chromatin contacts via a simple protocol of restriction enzyme digestion and religation of fragments upon agarose gel beads followed by paired-end sequencing. In addition to identifying contacts between pairs of loci, TM3C enables identification of contacts among more than two loci simultaneously. We use TM3C to assay the genome architectures of two human cell lines: KBM7, a near-haploid chronic leukemia cell line, and NHEK, a normal diploid human epidermal keratinocyte cell line. We confirm that the contact frequency maps produced by TM3C exhibit features characteristic of existing genome architecture datasets, including the expected scaling of contact probabilities with genomic distance, megabase scale chromosomal compartments and sub-megabase scale topological domains. We also confirm that TM3C captures several known cell type-specific contacts, ploidy shifts and translocations, such as Philadelphia chromosome formation (Ph+) in KBM7. We confirm a subset of the triple contacts involving the IGF2-H19 imprinting control region (ICR) using PCR analysis for KBM7 cells. Our genome-wide analysis of pairwise and triple contacts demonstrates their preference for linking open chromatin regions to each other and for linking regions with higher numbers of DNase hypersensitive sites (DHSs) to each other. For near-haploid KBM7 cells, we infer whole genome 3D models that exhibit clustering of small chromosomes with each other and large chromosomes with each other, consistent with previous studies of the genome architectures of other human cell lines.

Conclusion

TM3C is a simple protocol for ascertaining genome architecture and can be used to identify simultaneous contacts among three or four loci. Application of TM3C to a near-haploid human cell line revealed large-scale features of chromosomal organization and multi-way chromatin contacts that preferentially link regions of open chromatin.

Electronic supplementary material

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

Building and analyzing protein interactome networks by cross-species comparisons

Background  

A genomic catalogue of protein-protein interactions is a rich source of information, particularly for exploring the relationships between proteins. Numerous systems-wide and small-scale experiments have been conducted to identify interactions; however, our knowledge of all interactions for any one species is incomplete, and alternative means to expand these network maps is needed. We therefore took a comparative biology approach to predict protein-protein interactions across five species (human, mouse, fly, worm, and yeast) and developed InterologFinder for research biologists to easily navigate this data. We also developed a confidence score for interactions based on available experimental evidence and conservation across species.     

酿酒酵母ADH3基因的敲除

设计含有与酿酒酵母(Saccharomyces cerevisiae)编码乙醇脱氢酶Ⅲ的ADH3基因ORF两侧序列同源的长引物,以质粒pUG6为模板进行PCR构建带有Cre/loxP系统的敲除组件。转化酿酒酵母YS3(Saccharomyces cerevisiae),并将质粒pSH65转入阳性克隆子。半乳糖诱导表达Cre酶切除Kanr基因,在YPD培养基中连续传代培养丢失pSH65质粒,在原ORF处留下一个loxP位点,获得ADH3单倍体缺陷型菌株。利用同样的方法再次敲除双倍体的另一个等位基因。最终获得ADH3双倍体基因缺陷型突变株YS3-ADH3。     

SynteBase/SynteView: a tool to visualize gene order conservation in prokaryotic genomes

Background

Most microarray studies are made using labelling with one or two dyes which allows the hybridization of one or two samples on the same slide. In such experiments, the most frequently used dyes areCy3 andCy5. Recent improvements in the technology (dye-labelling, scanner and, image analysis) allow hybridization up to four samples simultaneously. The two additional dyes areAlexa488 andAlexa494. The triple-target or four-target technology is very promising, since it allows more flexibility in the design of experiments, an increase in the statistical power when comparing gene expressions induced by different conditions and a scaled down number of slides. However, there have been few methods proposed for statistical analysis of such data. Moreover the lowess correction of the global dye effect is available for only two-color experiments, and even if its application can be derived, it does not allow simultaneous correction of the raw data.

Results

We propose a two-step normalization procedure for triple-target experiments. First the dye bleeding is evaluated and corrected if necessary. Then the signal in each channel is normalized using a generalized lowess procedure to correct a global dye bias. The normalization procedure is validated using triple-self experiments and by comparing the results of triple-target and two-color experiments. Although the focus is on triple-target microarrays, the proposed method can be used to normalizepdifferently labelled targets co-hybridized on a same array, for any value ofpgreater than 2.

Conclusion

The proposed normalization procedure is effective: the technical biases are reduced, the number of false positives is under control in the analysis of differentially expressed genes, and the triple-target experiments are more powerful than the corresponding two-color experiments. There is room for improving the microarray experiments by simultaneously hybridizing more than two samples.     

CODEX: a normalization and copy number variation detection method for whole exome sequencing

High-throughput sequencing of DNA coding regions has become a common way of assaying genomic variation in the study of human diseases. Copy number variation (CNV) is an important type of genomic variation, but detecting and characterizing CNV from exome sequencing is challenging due to the high level of biases and artifacts. We propose CODEX, a normalization and CNV calling procedure for whole exome sequencing data. The Poisson latent factor model in CODEX includes terms that specifically remove biases due to GC content, exon capture and amplification efficiency, and latent systemic artifacts. CODEX also includes a Poisson likelihood-based recursive segmentation procedure that explicitly models the count-based exome sequencing data. CODEX is compared to existing methods on a population analysis of HapMap samples from the 1000 Genomes Project, and shown to be more accurate on three microarray-based validation data sets. We further evaluate performance on 222 neuroblastoma samples with matched normals and focus on a well-studied rare somatic CNV within the ATRX gene. We show that the cross-sample normalization procedure of CODEX removes more noise than normalizing the tumor against the matched normal and that the segmentation procedure performs well in detecting CNVs with nested structures.     

Polymorphic extracellular glucoamylase genes and their evolutionary origin in the yeast Saccharomyces diastaticus.

DNA coding for extracellular glucoamylase genes STA1 and STA3 was isolated from DNA libraries of two Saccharomyces diastaticus strains, each carrying STA1 or STA3. Cells transformed with a plasmid carrying either the STA1 or STA3 gene secreted glucoamylases having the same enzymatic and immunological properties and the same electrophoretic mobilities in acrylamide gel electrophoresis as those of authentic glucoamylases. Southern blot analysis of genomic DNA from S. diastaticus and a glucoamylase-non-secreting yeast, Saccharomyces cerevisiae, revealed that the STA1 and STA3 loci of S. diastaticus showed a high degree of homology, and that both yeast species (S. diastaticus and S. cerevisiae) contained DNA segments highly homologous to those of the extracellular glucoamylase genes. Restriction maps of the homologous DNA segments suggested that the extracellular glucoamylase genes of S. diastaticus may have arisen from recombination among the resident DNA segments in S. cerevisiae.     

The genome sequence of the wine yeast VIN7 reveals an allotriploid hybrid genome with Saccharomyces cerevisiae and Saccharomyces kudriavzevii origins

The vast majority of wine fermentations are performed principally by Saccharomyces cerevisiae. However, there are a growing number of instances in which other species of Saccharomyces play a predominant role. Interestingly, the presence of these other yeast species generally occurs via the formation of interspecific hybrids that contain genomic contributions from both S.?cerevisiae and non-S.?cerevisiae species. However, despite the large number of wine strains that are characterized at the genomic level, there remains limited information regarding the detailed genomic structure of hybrids used in winemaking. To address this, we describe the genome sequence of the thiol-releasing commercial wine yeast hybrid VIN7. VIN7 is shown to be an almost complete allotriploid interspecific hybrid that is comprised of a heterozygous diploid complement of S.?cerevisiae chromosomes and a haploid Saccharomyces kudriavzevii genomic contribution. Both parental strains appear to be of European origin, with the S.?cerevisiae parent being closely related to, but distinct from, the commercial wine yeasts QA23 and EC1118. In addition, several instances of chromosomal rearrangement between S.?cerevisiae and S.?kudriavzevii sequences were observed that may mark the early stages of hybrid genome consolidation.     

ORC- and Cdc6-dependent complexes at active and inactive chromosomal replication origins in Saccharomyces cerevisiae.

We have developed a genomic footprinting protocol which allows us to examine protein-DNA interactions at single copy chromosomal origins of DNA replication in the budding yeast Saccharomyces cerevisiae. We show that active replication origins oscillate between two chromatin states during the cell cycle: an origin recognition complex (ORC)-dependent post-replicative state and a Cdc6p-dependent pre-replicative state. Furthermore, we show that both post- and pre-replicative complexes can form efficiently on closely apposed replicators. Surprisingly, ARS301 which is active as an origin on plasmids but not in its normal chromosomal location, forms ORC- and Cdc6p-dependent complexes in both its active and inactive contexts. Thus, although ORC and Cdc6p are essential for initiation, their binding is not sufficient to dictate origin use.     

Inhibition of HIV derived lentiviral production by TAR RNA binding domain of TAT protein

Background

Retroviruses have a diploid genome and recombine at high frequency. Recombinant proviruses can be generated when two genetically different RNA genomes are packaged into the same retroviral particle. It was shown in several studies that recombinant proviruses could be generated in each round of HIV-1 replication, whereas the recombination rates of SNV and Mo-MuLV are 5 to 10-fold lower. The reason for these differences is not clear. One possibility is that these retroviruses may differ in their ability to copackage genomic RNAs produced at different chromosomal loci.

Results

To investigate whether there is a difference in the efficiency of heterodimer formation when two proviruses have the same or different chromosomal locations, we introduced two different Mo-MuLV-based retroviral vectors into the packaging cell line using either the cotransfection or sequential transfection procedure. The comparative study has shown that the frequency of recombination increased about four-fold when the cotransfection procedure was used. This difference was not associated with possible recombination of retroviral vectors during or after cotransfection and the ratios of retroviral virion RNAs were the same for two variants of transfection.

Conclusions

The results of this study indicate that a mechanism exists to enable the preferential copackaging of Mo-MuLV genomic RNA molecules that are transcribed on the same DNA template. The properties of Mo-MuLV genomic RNAs transport, processing or dimerization might be responsible for this preference. The data presented in this report can be useful when designing methods to study different aspects of replication and recombination of a diploid retroviral genome.     

Functional dissection of in vivo interchromosome association in Saccharomyces cerevisiae

Homologue pairing mediates both recombination and segregation of chromosomes at meiosis I. The recognition of nucleic-acid-sequence homology within the somatic nucleus has an impact on DNA repair and epigenetic control of gene expression. Here we investigate interchromosomal interactions using a non-invasive technique that allows tagging and visualization of DNA sequences in vegetative and meiotic live yeast cells. In non-meiotic cells, chromosomes are ordered in the nucleus, but preferential pairing between homologues is not observed. Association of tagged chromosomal domains occurs irrespective of their genomic location, with some preference for similar chromosomal positions. Here we describe a new phenomenon that promotes associations between sequence-identical ectopic tags with a tandem-repeat structure. These associations, termed interchromosome trans-associations, may underlie epigenetic phenomena.