Sequence Features Contributing to Chromosomal Rearrangements in Neisseria gonorrhoeae

Through whole genome sequence alignments, breakpoints in chromosomal synteny can be identified and the sequence features associated with these determined. Alignments of the genome sequences of Neisseria gonorrhoeae strain FA1090, N. gonorrhoeae strain NCCP11945, and N. gonorrhoeae strain TCDC-NG08107 reveal chromosomal rearrangements that have occurred. Based on these alignments and dot plot pair-wise comparisons, the overall chromosomal arrangement of strain NCCP11945 and TCDC-NG08107 are very similar, with no large inversions or translocations. The insertion of the Gonococcal Genetic Island in strain NCCP11945 is the most prominent distinguishing feature differentiating these strains. When strain NCCP11945 is compared to strain FA1090, however, 14 breakpoints in chromosomal synteny are identified between these gonococcal strains. The majority of these, 11 of 14, are associated with a prophage, IS elements, or IS-like repeat enclosed elements which appear to have played a role in the rearrangements observed. Additional rearrangements of small regions of the genome are associated with pilin genes. Evidence presented here suggests that the rearrangements of blocks of sequence are mediated by activation of prophage and associated IS elements and reintegration elsewhere in the genome or by homologous recombination between IS-like elements that have generated inversions.     

Reactivation of the ATCase domain of the URA2 gene complex: a positive selection method for Ty insertions and chromosomal rearrangements in Saccharomyces cerevisiae

Genetic rearrangements such as deletions or duplications of DNA sequences are rarely detected in the yeast Saccharomyces cerevisiae. We have developed a screening system using the URA2 gene coding for the bi-functional CPSase-ATCase (carbamyl phosphate synthetase — aspartate transcarbamylase) to select positively for these kinds of events. Nonsense mutations in the CPSase region cause a complete loss of the ATCase activity because of their strong polar effect. Thirty-seven ATCase⁺ revertants were isolated from a strain containing three nonsense mutations in the proximal CPSase region. Genetic and structural analysis of the URA2 locus in these strains allowed us to characterize two major classes of revertants. In the first, an entire copy of a Ty transposon was found to be inserted in the CPSase coding domain. This event, which represents a new form of Ty-mediated gene activation was further analysed by mapping the Ty integration site in 26 strains. In a second class of revertants, we observed chromosomal rearrangements and, in particular, duplication of the ATCase region and its integration in a new chromosomal environment in which this sequence becomes active.     

Defining the Sequence Elements and Candidate Genes for the Coloboma Mutation

The chicken coloboma mutation exhibits features similar to human congenital developmental malformations such as ocular coloboma, cleft-palate, dwarfism, and polydactyly. The coloboma-associated region and encoded genes were investigated using advanced genomic, genetic, and gene expression technologies. Initially, the mutation was linked to a 990 kb region encoding 11 genes; the application of the genetic and genomic tools led to a reduction of the linked region to 176 kb and the elimination of 7 genes. Furthermore, bioinformatics analyses of capture array-next generation sequence data identified genetic elements including SNPs, insertions, deletions, gaps, chromosomal rearrangements, and miRNA binding sites within the introgressed causative region relative to the reference genome sequence. Coloboma-specific variants within exons, UTRs, and splice sites were studied for their contribution to the mutant phenotype. Our compiled results suggest three genes for future studies. The three candidate genes, SLC30A5 (a zinc transporter), CENPH (a centromere protein), and CDK7 (a cyclin-dependent kinase), are differentially expressed (compared to normal embryos) at stages and in tissues affected by the coloboma mutation. Of these genes, two (SLC30A5 and CENPH) are considered high-priority candidate based upon studies in other vertebrate model systems.     

Conditional genomic rearrangement by designed meiotic recombination using VDE (PI-SceI) in yeast

Meiotic recombination plays critical roles in the acquisition of genetic diversity and has been utilized for conventional breeding of livestock and crops. The frequency of meiotic recombination is normally low, and is extremely low in regions called “recombination cold domains”. Here, we describe a new and highly efficient method to modulate yeast meiotic gene rearrangements using VDE (PI-SceI), an intein-encoded endonuclease that causes an efficient unidirectional meiotic gene conversion at its recognition sequence (VRS). We designed universal targeting vectors, by use of which the strain that inserts the VRS at a desired site is acquired. Meiotic induction of the strains provided unidirectional gene conversions and frequent genetic rearrangements of flanking genes with little impact on cell viability. This system thus opens the way for the designed modulation of meiotic gene rearrangements, regardless of recombinational activity of chromosomal domains. Finally, the VDE–VRS system enabled us to conduct meiosis-specific conditional knockout of genes where VDE-initiated gene conversion disrupts the target gene during meiosis, serving as a novel approach to examine the functions of genes during germination of resultant spores.     

Learning a Prior on Regulatory Potential from eQTL Data

Genome-wide RNA expression data provide a detailed view of an organism's biological state; hence, a dataset measuring expression variation between genetically diverse individuals (eQTL data) may provide important insights into the genetics of complex traits. However, with data from a relatively small number of individuals, it is difficult to distinguish true causal polymorphisms from the large number of possibilities. The problem is particularly challenging in populations with significant linkage disequilibrium, where traits are often linked to large chromosomal regions containing many genes. Here, we present a novel method, Lirnet, that automatically learns a regulatory potential for each sequence polymorphism, estimating how likely it is to have a significant effect on gene expression. This regulatory potential is defined in terms of “regulatory features”—including the function of the gene and the conservation, type, and position of genetic polymorphisms—that are available for any organism. The extent to which the different features influence the regulatory potential is learned automatically, making Lirnet readily applicable to different datasets, organisms, and feature sets. We apply Lirnet both to the human HapMap eQTL dataset and to a yeast eQTL dataset and provide statistical and biological results demonstrating that Lirnet produces significantly better regulatory programs than other recent approaches. We demonstrate in the yeast data that Lirnet can correctly suggest a specific causal sequence variation within a large, linked chromosomal region. In one example, Lirnet uncovered a novel, experimentally validated connection between Puf3—a sequence-specific RNA binding protein—and P-bodies—cytoplasmic structures that regulate translation and RNA stability—as well as the particular causative polymorphism, a SNP in Mkt1, that induces the variation in the pathway.     

Identification of a novel type of WRKY transcription factor binding site in elicitor-responsive cis-sequences from Arabidopsis thaliana

Using a combination of bioinformatics and synthetic promoters, novel elicitor-responsive cis-sequences were discovered in promoters of pathogen-upregulated genes from Arabidopsis thaliana. One group of functional sequences contains the conserved core sequence GACTTTT. This core sequence and adjacent nucleotides are essential for elicitor-responsive gene expression in a parsley protoplast system. By yeast one-hybrid screening, WRKY70 was selected with a cis-sequence harbouring the core sequence GACTTTT but no known WRKY binding site (W-box). Transactivation experiments, mutation analyses, and electrophoretic mobility shift assays demonstrate that the sequence CGACTTTT is the binding site for WRKY70 in the investigated cis-sequence and is required for WRKY70-activated gene expression. Using several cis-sequences in transactivation experiments and binding studies, the CGACTTTT sequence can be extended to propose YGACTTTT as WRKY70 binding site. This binding site, designated WT-box, is enriched in promoters of genes upregulated in a WRKY70 overexpressing line. Interestingly, functional WRKY70 binding sites are present in the promoter of WRKY30, supporting recent evidence that both factors play a role in the same regulatory network.     

Rad51p and Rad54p,but not Rad52p,elevate gene repair in Saccharomyces cerevisiae directed by modified single-stranded oligonucleotide vectors

Synthetic single-stranded DNA vectors have been used to correct point and frameshift mutations in episomal or chromosomal targets in the yeast Saccharomyces cerevisiae. Certain parameters, such as the length of the vector and the genetic background of the organism, have a significant impact on the process of targeted gene repair, and point mutations are corrected at a higher frequency than frameshift mutations. Genetic analyses reveal that expression levels of the recombination/repair genes RAD51, RAD52 and RAD54 can affect the frequency of gene repair. Overexpression of RAD51 enhances the frequency 4-fold for correction of an episomal target and 5-fold for correction of a chromosomal target; overexpression of RAD54 is also effective in stimulating gene repair, to the same extent as RAD51 in the chromosomal target. In sharp contrast, RAD52 gene expression serves to reduce gene repair activity in rescue experiments and in experiments where RAD52 is overexpressed in a wild-type strain. This may suggest an antagonist role for Rad52p. Consistent with this notion, the highest level of targeted repair occurs when the RAD51 gene is overexpressed in a strain of yeast deficient in RAD52 gene function.     

Inactivation of Mre11 does not affect VSG gene duplication mediated by homologous recombination in Trypanosoma brucei

We demonstrate, by gene deletion analysis, that Mre11 has a critical role in maintaining genomic integrity in Trypanosoma brucei. mre11(-/-) null mutant strains exhibited retarded growth but no delay or disruption of cell cycle progression. They showed also a weak hyporecombination phenotype and the accumulation of gross chromosomal rearrangements, which did not involve sequence translocation, telomere loss, or formation of new telomeres. The trypanosome mre11(-/-) strains were hypersensitive to phleomycin, a mutagen causing DNA double strand breaks (DSBs) but, in contrast to mre11(-/-) null mutants in other organisms and T. brucei rad51(-/-) null mutants, displayed no hypersensitivity to methyl methanesulfonate, which causes point mutations and DSBs. Mre11 therefore is important for the repair of chromosomal damage and DSBs in trypanosomes, although in this organism the intersection of repair pathways appears to differ from that in other organisms. Mre11 inactivation appears not to affect VSG gene switching during antigenic variation of a laboratory strain, which is perhaps surprising given the importance of homologous recombination during this process.     

Rapid Evolution of Recombinant Saccharomyces cerevisiae for Xylose Fermentation through Formation of Extra-chromosomal Circular DNA

Circular DNA elements are involved in genome plasticity, particularly of tandem repeats. However, amplifications of DNA segments in Saccharomyces cerevisiae reported so far involve pre-existing repetitive sequences such as ribosomal DNA, Ty elements and Long Terminal Repeats (LTRs). Here, we report the generation of an eccDNA, (extrachromosomal circular DNA element) in a region without any repetitive sequences during an adaptive evolution experiment. We performed whole genome sequence comparison between an efficient D-xylose fermenting yeast strain developed by metabolic and evolutionary engineering, and its parent industrial strain. We found that the heterologous gene XylA that had been inserted close to an ARS sequence in the parent strain has been amplified about 9 fold in both alleles of the chromosomal locus of the evolved strain compared to its parent. Analysis of the amplification process during the adaptive evolution revealed formation of a XylA-carrying eccDNA, pXI2-6, followed by chromosomal integration in tandem arrays over the course of the evolutionary adaptation. Formation of the eccDNA occurred in the absence of any repetitive DNA elements, probably using a micro-homology sequence of 8 nucleotides flanking the amplified sequence. We isolated the pXI2-6 eccDNA from an intermediate strain of the evolutionary adaptation process, sequenced it completely and showed that it confers high xylose fermentation capacity when it is transferred to a new strain. In this way, we have provided clear evidence that gene amplification can occur through generation of eccDNA without the presence of flanking repetitive sequences and can serve as a rapid means of adaptation to selection pressure.     

Expansion and contraction of the DUP240 multigene family in Saccharomyces cerevisiae populations

The influence of duplicated sequences on chromosomal stability is poorly understood. To characterize chromosomal rearrangements involving duplicated sequences, we compared the organization of tandem repeats of the DUP240 gene family in 15 Saccharomyces cerevisiae strains of various origins. The DUP240 gene family consists of 10 members of unknown function in the reference strain S288C. Five DUP240 paralogs on chromosome I and two on chromosome VII are arranged as tandem repeats that are highly polymorphic in copy number and sequence. We characterized DNA sequences that are likely involved in homologous or nonhomologous recombination events and are responsible for intra- and interchromosomal rearrangements that cause the creation and disappearance of DUP240 paralogs. The tandemly repeated DUP240 genes seem to be privileged sites of gene birth and death.     

Genome Sequence of the Lager Brewing Yeast, an Interspecies Hybrid

This work presents the genome sequencing of the lager brewing yeast (Saccharomyces pastorianus) Weihenstephan 34/70, a strain widely used in lager beer brewing. The 25 Mb genome comprises two nuclear sub-genomes originating from Saccharomyces cerevisiae and Saccharomyces bayanus and one circular mitochondrial genome originating from S. bayanus. Thirty-six different types of chromosomes were found including eight chromosomes with translocations between the two sub-genomes, whose breakpoints are within the orthologous open reading frames. Several gene loci responsible for typical lager brewing yeast characteristics such as maltotriose uptake and sulfite production have been increased in number by chromosomal rearrangements. Despite an overall high degree of conservation of the synteny with S. cerevisiae and S. bayanus, the syntenies were not well conserved in the sub-telomeric regions that contain lager brewing yeast characteristic and specific genes. Deletion of larger chromosomal regions, a massive unilateral decrease of the ribosomal DNA cluster and bilateral truncations of over 60 genes reflect a post-hybridization evolution process. Truncations and deletions of less efficient maltose and maltotriose uptake genes may indicate the result of adaptation to brewing. The genome sequence of this interspecies hybrid yeast provides a new tool for better understanding of lager brewing yeast behavior in industrial beer production.Key words: Saccharomyces pastorianus, beer, genome, interspecies hybrid, larger yeast     

Convergent adaptation of Saccharomyces uvarum to sulfite,an antimicrobial preservative widely used in human-driven fermentations

Different species can find convergent solutions to adapt their genome to the same evolutionary constraints, although functional convergence promoted by chromosomal rearrangements in different species has not previously been found. In this work, we discovered that two domesticated yeast species, Saccharomyces cerevisiae, and Saccharomyces uvarum, acquired chromosomal rearrangements to convergently adapt to the presence of sulfite in fermentation environments. We found two new heterologous chromosomal translocations in fermentative strains of S. uvarum at the SSU1 locus, involved in sulfite resistance, an antimicrobial additive widely used in food production. These are convergent events that share similarities with other SSU1 locus chromosomal translocations previously described in domesticated S. cerevisiae strains. In S. uvarum, the newly described VII^XVI and XI^XVI chromosomal translocations generate an overexpression of the SSU1 gene and confer increased sulfite resistance. This study highlights the relevance of chromosomal rearrangements to promote the adaptation of yeast to anthropic environments.     

Functional Knowledge Transfer for High-accuracy Prediction of Under-studied Biological Processes

A key challenge in genetics is identifying the functional roles of genes in pathways. Numerous functional genomics techniques (e.g. machine learning) that predict protein function have been developed to address this question. These methods generally build from existing annotations of genes to pathways and thus are often unable to identify additional genes participating in processes that are not already well studied. Many of these processes are well studied in some organism, but not necessarily in an investigator''s organism of interest. Sequence-based search methods (e.g. BLAST) have been used to transfer such annotation information between organisms. We demonstrate that functional genomics can complement traditional sequence similarity to improve the transfer of gene annotations between organisms. Our method transfers annotations only when functionally appropriate as determined by genomic data and can be used with any prediction algorithm to combine transferred gene function knowledge with organism-specific high-throughput data to enable accurate function prediction.We show that diverse state-of-art machine learning algorithms leveraging functional knowledge transfer (FKT) dramatically improve their accuracy in predicting gene-pathway membership, particularly for processes with little experimental knowledge in an organism. We also show that our method compares favorably to annotation transfer by sequence similarity. Next, we deploy FKT with state-of-the-art SVM classifier to predict novel genes to 11,000 biological processes across six diverse organisms and expand the coverage of accurate function predictions to processes that are often ignored because of a dearth of annotated genes in an organism. Finally, we perform in vivo experimental investigation in Danio rerio and confirm the regulatory role of our top predicted novel gene, wnt5b, in leftward cell migration during heart development. FKT is immediately applicable to many bioinformatics techniques and will help biologists systematically integrate prior knowledge from diverse systems to direct targeted experiments in their organism of study.     

RAPYD--rapid annotation platform for yeast data

Lower eukaryotes of the kingdom Fungi include a variety of biotechnologically important yeast species that are in the focus of genome research for more than a decade. Due to the rapid progress in ultra-fast sequencing technologies, the amount of available yeast genome data increases steadily. Thus, an efficient bioinformatics platform is required that covers genome assembly, eukaryotic gene prediction, genome annotation, comparative yeast genomics, and metabolic pathway reconstruction. Here, we present a bioinformatics platform for yeast genomics named RAPYD addressing the key requirements of extensive yeast sequence data analysis. The first step is a comprehensive regional and functional annotation of a yeast genome. A region prediction pipeline was implemented to obtain reliable and high-quality predictions of coding sequences and further genome features. Functions of coding sequences are automatically determined using a configurable prediction pipeline. Based on the resulting functional annotations, a metabolic pathway reconstruction module can be utilized to rapidly generate an overview of organism-specific features and metabolic blueprints. In a final analysis step shared and divergent features of closely related yeast strains can be explored using the comparative genomics module. An in-depth application example of the yeast Meyerozyma guilliermondii illustrates the functionality of RAPYD. A user-friendly web interface is available at https://rapyd.cebitec.uni-bielefeld.de.     

Mechanisms for gene conversion and homologous recombination: The double-strand break repair model and the successive half crossing-over model

Two mechanisms for gene conversion and homologous recombination were discussed. (1) The double-strand break repair model. A double-strand break is expanded to a gap, which is then repaired by copying a homologous sequence. The gene conversion is often accompanied by crossing-over of the flanking sequences. We obtained evidence for this model in Red pathway of bacteriophage lambda and RecE pathway of E. coli. (2) The successive half crossing-over model. Half crossing-over leaves one recombinant duplex and one or two end(s) out of two parental duplexes. The resulting ends are, in turn, recombinogenic. Successive rounds of the half crossing-over mechanism explains why apparent plasmid gene conversion in RecF pathway of E. coli is not accompanied by crossing-over. This model can explain chromosomal gene conversion if we assume that the donor is first replicated. Gene conversion during mating-type switching in yeast, antigenic variation in unicellular microorganisms, and chromosomal gene conversion in mammalian somatic cells are explained by this model. Distinguishing between these two mechanisms is important in understanding recombination in yeast and mammalian cells and also in its application to gene targeting.     

Reactivation of the ATCase domain of the URA2 gene complex: a positive selection method for Ty insertions and chromosomal rearrangements in Saccharomyces cerevisiae

Genetic rearrangements such as deletions or duplications of DNA sequences are rarely detected in the yeast Saccharomyces cerevisiae. We have developed a screening system using the URA2 gene coding for the bi-functional CPSase-ATCase (carbamyl phosphate synthetase — aspartate transcarbamylase) to select positively for these kinds of events. Nonsense mutations in the CPSase region cause a complete loss of the ATCase activity because of their strong polar effect. Thirty-seven ATCase⁺ revertants were isolated from a strain containing three nonsense mutations in the proximal CPSase region. Genetic and structural analysis of the URA2 locus in these strains allowed us to characterize two major classes of revertants. In the first, an entire copy of a Ty transposon was found to be inserted in the CPSase coding domain. This event, which represents a new form of Ty-mediated gene activation was further analysed by mapping the Ty integration site in 26 strains. In a second class of revertants, we observed chromosomal rearrangements and, in particular, duplication of the ATCase region and its integration in a new chromosomal environment in which this sequence becomes active.     

Mismatch Repair Proteins Regulate Heteroduplex Formation during Mitotic Recombination in Yeast

Mismatch repair (MMR) proteins actively inhibit recombination between diverged sequences in both prokaryotes and eukaryotes. Although the molecular basis of the antirecombination activity exerted by MMR proteins is unclear, it presumably involves the recognition of mismatches present in heteroduplex recombination intermediates. This recognition could be exerted during the initial stage of strand exchange, during the extension of heteroduplex DNA, or during the resolution of recombination intermediates. We previously used an assay system based on 350-bp inverted-repeat substrates to demonstrate that MMR proteins strongly inhibit mitotic recombination between diverged sequences in Saccharomyces cerevisiae. The assay system detects only those events that reverse the orientation of the region between the recombination substrates, which can occur as a result of either intrachromatid crossover or sister chromatid conversion. In the present study we sequenced the products of mitotic recombination between 94%-identical substrates in order to map gene conversion tracts in wild-type versus MMR-defective yeast strains. The sequence data indicate that (i) most recombination occurs via sister chromatid conversion and (ii) gene conversion tracts in an MMR-defective strain are significantly longer than those in an isogenic wild-type strain. The shortening of conversion tracts observed in a wild-type strain relative to an MMR-defective strain suggests that at least part of the antirecombination activity of MMR proteins derives from the blockage of heteroduplex extension in the presence of mismatches.     

Segments missing from the draft human genome sequence can be isolated by transformation-associated recombination cloning in yeast

The reported draft human genome sequence includes many contigs that are separated by gaps of unknown sequence. These gaps may be due to chromosomal regions that are not present in the Escherichia coli libraries used for DNA sequencing because they cannot be cloned efficiently, if at all, in bacteria. Using a yeast artificial chromosome (YAC)/ bacterial artificial chromosome (BAC) library generated in yeast, we found that approximately 6% of human DNA sequences tested transformed E. coli cells less efficiently than yeast cells, and were less stable in E. coli than in yeast. When the ends of several YAC/BAC isolates cloned in yeast were sequenced and compared with the reported draft sequence, major inconsistencies were found with the sequences of those YAC/BAC isolates that transformed E. coli cells inefficiently. Two human genomic fragments were re-isolated from human DNA by transformation-associated recombination (TAR) cloning. Re-sequencing of these regions showed that the errors in the draft are the results of both missassembly and loss of specific DNA sequences during cloning in E. coli. These results show that TAR cloning might be a valuable method that could be widely used during the final stages of the Human Genome Project.