Recombination dynamics of a human Y-chromosomal palindrome:rapid GC-biased gene conversion, multi-kilobase conversion tracts, and rare inversions

The male-specific region of the human Y chromosome (MSY) includes eight large inverted repeats (palindromes) in which arm-to-arm similarity exceeds 99.9%, due to gene conversion activity. Here, we studied one of these palindromes, P6, in order to illuminate the dynamics of the gene conversion process. We genotyped ten paralogous sequence variants (PSVs) within the arms of P6 in 378 Y chromosomes whose evolutionary relationships within the SNP-defined Y phylogeny are known. This allowed the identification of 146 historical gene conversion events involving individual PSVs, occurring at a rate of 2.9-8.4×10(-4) events per generation. A consideration of the nature of nucleotide change and the ancestral state of each PSV showed that the conversion process was significantly biased towards the fixation of G or C nucleotides (GC-biased), and also towards the ancestral state. Determination of haplotypes by long-PCR allowed likely co-conversion of PSVs to be identified, and suggested that conversion tract lengths are large, with a mean of 2068 bp, and a maximum in excess of 9 kb. Despite the frequent formation of recombination intermediates implied by the rapid observed gene conversion activity, resolution via crossover is rare: only three inversions within P6 were detected in the sample. An analysis of chimpanzee and gorilla P6 orthologs showed that the ancestral state bias has existed in all three species, and comparison of human and chimpanzee sequences with the gorilla outgroup confirmed that GC bias of the conversion process has apparently been active in both the human and chimpanzee lineages

A Model-Based Analysis of GC-Biased Gene Conversion in the Human and Chimpanzee Genomes

GC-biased gene conversion (gBGC) is a recombination-associated process that favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is hypothesized to contribute to variation in GC content, rapidly evolving sequences, and the fixation of deleterious mutations, but its prevalence and general functional consequences remain poorly understood. gBGC is difficult to incorporate into models of molecular evolution and so far has primarily been studied using summary statistics from genomic comparisons. Here, we introduce a new probabilistic model that captures the joint effects of natural selection and gBGC on nucleotide substitution patterns, while allowing for correlations along the genome in these effects. We implemented our model in a computer program, called phastBias, that can accurately detect gBGC tracts about 1 kilobase or longer in simulated sequence alignments. When applied to real primate genome sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide differences. These tracts fall in clusters, particularly in subtelomeric regions; they are enriched for recombination hotspots and fast-evolving sequences; and they display an ongoing fixation preference for G and C alleles. They are also significantly enriched for disease-associated polymorphisms, suggesting that they contribute to the fixation of deleterious alleles. The gBGC tracts provide a unique window into historical recombination processes along the human and chimpanzee lineages. They supply additional evidence of long-term conservation of megabase-scale recombination rates accompanied by rapid turnover of hotspots. Together, these findings shed new light on the evolutionary, functional, and disease implications of gBGC. The phastBias program and our predicted tracts are freely available. © 2013 Capra et al

ReCombine: A Suite of Programs for Detection and Analysis of Meiotic Recombination in Whole-Genome Datasets

In meiosis, the exchange of DNA between chromosomes by homologous recombination is a critical step that ensures proper chromosome segregation and increases genetic diversity. Products of recombination include reciprocal exchanges, known as crossovers, and non-reciprocal gene conversions or non-crossovers. The mechanisms underlying meiotic recombination remain elusive, largely because of the difficulty of analyzing large numbers of recombination events by traditional genetic methods. These traditional methods are increasingly being superseded by high-throughput techniques capable of surveying meiotic recombination on a genome-wide basis. Next-generation sequencing or microarray hybridization is used to genotype thousands of polymorphic markers in the progeny of hybrid yeast strains. New computational tools are needed to perform this genotyping and to find and analyze recombination events. We have developed a suite of programs, ReCombine, for using short sequence reads from next-generation sequencing experiments to genotype yeast meiotic progeny. Upon genotyping, the program CrossOver, a component of ReCombine, then detects recombination products and classifies them into categories based on the features found at each location and their distribution among the various chromatids. CrossOver is also capable of analyzing segregation data from microarray experiments or other sources. This package of programs is designed to allow even researchers without computational expertise to use high-throughput, whole-genome methods to study the molecular mechanisms of meiotic recombination

Budding yeast ATM/ATR control meiotic double-strand break (DSB) levels by down-regulating Rec114, an essential component of the DSB-machinery

An essential feature of meiosis is Spo11 catalysis of programmed DNA double strand breaks (DSBs). Evidence suggests that the number of DSBs generated per meiosis is genetically determined and that this ability to maintain a pre-determined DSB level, or "DSB homeostasis", might be a property of the meiotic program. Here, we present direct evidence that Rec114, an evolutionarily conserved essential component of the meiotic DSB-machinery, interacts with DSB hotspot DNA, and that Tel1 and Mec1, the budding yeast ATM and ATR, respectively, down-regulate Rec114 upon meiotic DSB formation through phosphorylation. Mimicking constitutive phosphorylation reduces the interaction between Rec114 and DSB hotspot DNA, resulting in a reduction and/or delay in DSB formation. Conversely, a non-phosphorylatable rec114 allele confers a genome-wide increase in both DSB levels and in the interaction between Rec114 and the DSB hotspot DNA. These observations strongly suggest that Tel1 and/or Mec1 phosphorylation of Rec114 following Spo11 catalysis down-regulates DSB formation by limiting the interaction between Rec114 and DSB hotspots. We also present evidence that Ndt80, a meiosis specific transcription factor, contributes to Rec114 degradation, consistent with its requirement for complete cessation of DSB formation. Loss of Rec114 foci from chromatin is associated with homolog synapsis but independent of Ndt80 or Tel1/Mec1 phosphorylation. Taken together, we present evidence for three independent ways of regulating Rec114 activity, which likely contribute to meiotic DSBs-homeostasis in maintaining genetically determined levels of breaks

The effect of tightly-bound water molecules on scaffold diversity in computer-aided de novo ligand design of CDK2 inhibitors

We have determined the effects that tightly bound water molecules have on the de novo design of cyclin-dependent kinase-2 (CDK2) ligands. In particular, we have analyzed the impact of a specific structural water molecule on the chemical diversity and binding mode of ligands generated through a de novo structure-based ligand generation method in the binding site of CDK2. The tightly bound water molecule modifies the size and shape of the binding site and we have found that it also imposed constraints on the observed binding modes of the generated ligands. This in turn had the indirect effect of reducing the chemical diversity of the underlying molecular scaffolds that were able to bind to the enzyme satisfactorily

Morpholino Gene Knockdown in Adult Fundulus heteroclitus: Role of SGK1 in Seawater Acclimation

The Atlantic killifish (Fundulus heteroclitus) is an environmental sentinel organism used extensively for studies on environmental toxicants and salt (NaCl) homeostasis. Previous research in our laboratory has shown that rapid acclimation of killifish to seawater is mediated by trafficking of CFTR chloride channels from intracellular vesicles to the plasma membrane in the opercular membrane within the first hour in seawater, which enhances chloride secretion into seawater, thereby contributing to salt homeostasis. Acute transition to seawater is also marked by an increase in both mRNA and protein levels of serum glucocorticoid kinase 1 (SGK1) within 15 minutes of transfer. Although the rise in SGK1 in gill and its functional analog, the opercular membrane, after seawater transfer precedes the increase in membrane CFTR, a direct role of SGK1 in elevating membrane CFTR has not been established in vivo. To test the hypothesis that SGK1 mediates the increase in plasma membrane CFTR we designed two functionally different vivo-morpholinos to knock down SGK1 in gill, and developed and validated a vivo-morpholino knock down technique for adult killifish. Injection (intraperitoneal, IP) of the splice blocking SGK1 vivo-morpholino reduced SGK1 mRNA in the gill after transition from fresh to seawater by 66%. The IP injection of the translational blocking and splice blocking vivo-morpholinos reduced gill SGK1 protein abundance in fish transferred from fresh to seawater by 64% and 53%, respectively. Moreover, knock down of SGK1 completely eliminated the seawater induced rise in plasma membrane CFTR, demonstrating that the increase in SGK1 protein is required for the trafficking of CFTR from intracellular vesicles in mitochondrion rich cells to the plasma membrane in the gill during acclimation to seawater. This is the first report of the use of vivo-morpholinos in adult killifish and demonstrates that vivo-morpholinos are a valuable genetic tool for this environmentally relevant model organism

The surprising negative correlation of gene length and optimal codon use - disentangling translational selection from GC-biased gene conversion in yeast

<p>Abstract</p> <p>Background</p> <p>Surprisingly, in several multi-cellular eukaryotes optimal codon use correlates negatively with gene length. This contrasts with the expectation under selection for translational accuracy. While suggested explanations focus on variation in strength and efficiency of translational selection, it has rarely been noticed that the negative correlation is reported only in organisms whose optimal codons are biased towards codons that end with G or C (-GC). This raises the question whether forces that affect base composition - such as GC-biased gene conversion - contribute to the negative correlation between optimal codon use and gene length.</p> <p>Results</p> <p>Yeast is a good organism to study this as equal numbers of optimal codons end in -GC and -AT and one may hence compare frequencies of optimal GC- with optimal AT-ending codons to disentangle the forces. Results of this study demonstrate in yeast frequencies of GC-ending (optimal AND non-optimal) codons decrease with gene length and increase with recombination. A decrease of GC-ending codons along genes contributes to the negative correlation with gene length. Correlations with recombination and gene expression differentiate between GC-ending and optimal codons, and also substitution patterns support effects of GC-biased gene conversion.</p> <p>Conclusion</p> <p>While the general effect of GC-biased gene conversion is well known, the negative correlation of optimal codon use with gene length has not been considered in this context before. Initiation of gene conversion events in promoter regions and the presence of a gene conversion gradient most likely explain the observed decrease of GC-ending codons with gene length and gene position.</p

High-throughput 454 resequencing for allele discovery and recombination mapping in Plasmodium falciparum

<p>Abstract</p> <p>Background</p> <p>Knowledge of the origins, distribution, and inheritance of variation in the malaria parasite (<it>Plasmodium falciparum</it>) genome is crucial for understanding its evolution; however the 81% (A+T) genome poses challenges to high-throughput sequencing technologies. We explore the viability of the Roche 454 Genome Sequencer FLX (GS FLX) high throughput sequencing technology for both whole genome sequencing and fine-resolution characterization of genetic exchange in malaria parasites.</p> <p>Results</p> <p>We present a scheme to survey recombination in the haploid stage genomes of two sibling parasite clones, using whole genome pyrosequencing that includes a sliding window approach to predict recombination breakpoints. Whole genome shotgun (WGS) sequencing generated approximately 2 million reads, with an average read length of approximately 300 bp. <it>De novo </it>assembly using a combination of WGS and 3 kb paired end libraries resulted in contigs ≤ 34 kb. More than 8,000 of the 24,599 SNP markers identified between parents were genotyped in the progeny, resulting in a marker density of approximately 1 marker/3.3 kb and allowing for the detection of previously unrecognized crossovers (COs) and many non crossover (NCO) gene conversions throughout the genome.</p> <p>Conclusions</p> <p>By sequencing the 23 Mb genomes of two haploid progeny clones derived from a genetic cross at more than 30× coverage, we captured high resolution information on COs, NCOs and genetic variation within the progeny genomes. This study is the first to resequence progeny clones to examine fine structure of COs and NCOs in malaria parasites.</p

Mimotopes selected with a neutralizing antibody against urease B from Helicobacter pylori induce enzyme inhibitory antibodies in mice upon vaccination

<p>Abstract</p> <p>Background</p> <p>Urease B is an important virulence factor that is required for <it>Helicobacter pylori </it>to colonise the gastric mucosa. Mouse monoclonal antibodies (mAbs) that inhibit urease B enzymatic activity will be useful as vaccines for the prevention and treatment of <it>H. pylori </it>infection. Here, we produced murine mAbs against urease B that neutralize the enzyme's activity. We mapped their epitopes by phage display libraries and investigated the immunogenicity of the selected mimotopes <it>in vivo</it>.</p> <p>Results</p> <p>The urease B gene was obtained (GenBank accession No. <ext-link ext-link-id="DQ141576" ext-link-type="gen">DQ141576</ext-link>) and the recombinant pGEX-4T-1/UreaseB protein was expressed in <it>Escherichia coli </it>as a 92-kDa recombinant fusion protein with glutathione-S-transferase (GST). Five mAbs U001-U005 were produced by a hybridoma-based technique with urease B-GST as an immunogen. Only U001 could inhibit urease B enzymatic activity. Immunoscreening via phage display libraries revealed two different mimotopes of urease B protein; EXXXHDM from ph.D.12-library and EXXXHSM from ph.D.C7C that matched the urease B proteins at 347-353 aa. The antiserum induced by selected phage clones clearly recognised the urease B protein and inhibited its enzymatic activity, which indicated that the phagotope-induced immune responses were antigen specific.</p> <p>Conclusions</p> <p>The present work demonstrated that phage-displayed mimotopes were accessible to the mouse immune system and triggered a humoral response. The urease B mimotope could provide a novel and promising approach for the development of a vaccine for the diagnosis and treatment of <it>H. pylori </it>infection.</p

Mimotopes selected with neutralizing antibodies against multiple subtypes of influenza A

<p>Abstract</p> <p>Background</p> <p>The mimotopes of viruses are considered as the good targets for vaccine design. We prepared mimotopes against multiple subtypes of influenza A and evaluate their immune responses in flu virus challenged Balb/c mice.</p> <p>Methods</p> <p>The mimotopes of influenza A including pandemic H1N1, H3N2, H2N2 and H1N1 swine-origin influenza virus were screened by peptide phage display libraries, respectively. These mimotopes were engineered in one protein as multi- epitopes in Escherichia coli (E. coli) and purified. Balb/c mice were immunized using the multi-mimotopes protein and specific antibody responses were analyzed using hemagglutination inhibition (HI) assay and enzyme-linked immunosorbent assay (ELISA). The lung inflammation level was evaluated by hematoxylin and eosin (HE).</p> <p>Results</p> <p>Linear heptopeptide and dodecapeptide mimotopes were obtained for these influenza virus. The recombinant multi-mimotopes protein was a 73 kDa fusion protein. Comparing immunized infected groups with unimmunized infected subsets, significant differences were observed in the body weight loss and survival rate. The antiserum contained higher HI Ab titer against H1N1 virus and the lung inflammation level were significantly decreased in immunized infected groups.</p> <p>Conclusions</p> <p>Phage-displayed mimotopes against multiple subtypes of influenza A were accessible to the mouse immune system and triggered a humoral response to above virus.</p