Identification and utilization of arbitrary correlations in models of recombination signal sequences

BACKGROUND: A significant challenge in bioinformatics is to develop methods for detecting and modeling patterns in variable DNA sequence sites, such as protein-binding sites in regulatory DNA. Current approaches sometimes perform poorly when positions in the site do not independently affect protein binding. We developed a statistical technique for modeling the correlation structure in variable DNA sequence sites. The method places no restrictions on the number of correlated positions or on their spatial relationship within the site. No prior empirical evidence for the correlation structure is necessary. RESULTS: We applied our method to the recombination signal sequences (RSS) that direct assembly of B-cell and T-cell antigen-receptor genes via V(D)J recombination. The technique is based on model selection by cross-validation and produces models that allow computation of an information score for any signal-length sequence. We also modeled RSS using order zero and order one Markov chains. The scores from all models are highly correlated with measured recombination efficiencies, but the models arising from our technique are better than the Markov models at discriminating RSS from non-RSS. CONCLUSIONS: Our model-development procedure produces models that estimate well the recombinogenic potential of RSS and are better at RSS recognition than the order zero and order one Markov models. Our models are, therefore, valuable for studying the regulation of both physiologic and aberrant V(D)J recombination. The approach could be equally powerful for the study of promoter and enhancer elements, splice sites, and other DNA regulatory sites that are highly variable at the level of individual nucleotide positions

Prospective Estimation of Recombination Signal Efficiency and Identification of Functional Cryptic Signals in the Genome by Statistical Modeling

The recombination signals (RS) that guide V(D)J recombination are phylogenetically conserved but retain a surprising degree of sequence variability, especially in the nonamer and spacer. To characterize RS variability, we computed the position-wise information, a measure correlated with sequence conservation, for each nucleotide position in an RS alignment and demonstrate that most position-wise information is present in the RS heptamers and nonamers. We have previously demonstrated significant correlations between RS positions and here show that statistical models of the correlation structure that underlies RS variability efficiently identify physiologic and cryptic RS and accurately predict the recombination efficiencies of natural and synthetic RS. In scans of mouse and human genomes, these models identify a highly conserved family of repetitive DNA as an unexpected source of frequent, cryptic RS that rearrange both in extrachromosomal substrates and in their genomic context

Multiple, conserved cryptic recombination signals in VH gene segments: detection of cleavage products only in pro–B cells

Receptor editing is believed to play the major role in purging newly formed B cell compartments of autoreactivity by the induction of secondary V(D)J rearrangements. In the process of immunoglobulin heavy (H) chain editing, these secondary rearrangements are mediated by direct VH-to-JH joining or cryptic recombination signals (cRSs) within VH gene segments. Using a statistical model of RS, we have identified potential cRSs within VH gene segments at conserved sites flanking complementarity-determining regions 1 and 2. These cRSs are active in extrachromosomal recombination assays and cleaved during normal B cell development. Cleavage of multiple VH cRSs was observed in the bone marrow of C57BL/6 and RAG2:GFP and μMT congenic animals, and we determined that cRS cleavage efficiencies are 30–50-fold lower than a physiological RS. cRS signal ends are abundant in pro–B cells, including those recovered from μMT mice, but undetectable in pre– or immature B cells. Thus, VH cRS cleavage regularly occurs before the generation of functional preBCR and BCR. Conservation of cRSs distal from the 3′ end of VH gene segments suggests a function for these cryptic signals other than VH gene replacement

Structure-Guided Molecular Grafting Of A Complex Broadly Neutralizing Viral Epitope

Antigenic variation and viral evolution have thwarted traditional influenza vaccination strategies. The broad protection afforded by a “universal” influenza vaccine may come from immunogens that elicit humoral immune responses targeting conserved epitopes on the viral hemagglutinin (HA), such as the receptor-binding site (RBS). Here, we engineered candidate immunogens that use noncirculating, avian influenza HAs as molecular scaffolds to present the broadly neutralizing RBS epitope from historical, circulating H1 influenzas. These “resurfaced” HAs (rsHAs) remove epitopes potentially targeted by strain-specific responses in immune-experienced individuals. Through structure-guided optimization, we improved two antigenically different scaffolds to bind a diverse panel of pan-H1 and H1/H3 cross-reactive bnAbs with high affinity. Subsequent serological and single germinal center B cell analyses from murine prime-boost immunizations show that the rsHAs are both immunogenic and can augment the quality of elicited RBS-directed antibodies. Our structure-guided, RBS grafting approach provides candidate immunogens for selectively presenting a conserved viral epitope

A functional variant in the serotonin receptor 7 gene (HTR7), rs7905446, is associated with good response to SSRIs in bipolar and unipolar depression.

Predicting antidepressant response has been a clinical challenge for mood disorder. Although several genome-wide association studies have suggested a number of genetic variants to be associated with antidepressant response, the sample sizes are small and the results are difficult to replicate. Previous animal studies have shown that knockout of the serotonin receptor 7 gene (HTR7) resulted in an antidepressant-like phenotype, suggesting it was important to antidepressant action. In this report, in the first stage, we used a cost-effective pooled-sequencing strategy to sequence the entire HTR7 gene and its regulatory regions to investigate the association of common variants in HTR7 and clinical response to four selective serotonin reuptake inhibitors (SSRIs: citalopram, paroxetine, fluoxetine and sertraline) in a retrospective cohort mainly consisting of subjects with bipolar disorder (n = 359). We found 80 single-nucleotide polymorphisms (SNPs) with false discovery rate < 0.05 associated with response to paroxetine. Among the significant SNPs, rs7905446 (T/G), which is located at the promoter region, also showed nominal significance (P < 0.05) in fluoxetine group. GG/TG genotypes for rs7905446 and female gender were associated with better response to two SSRIs (paroxetine and fluoxetine). In the second stage, we replicated this association in two independent prospective samples of SSRI-treated patients with major depressive disorder: the MARS (n = 253, P = 0.0169) and GENDEP studies (n = 432, P = 0.008). The GG/TG genotypes were consistently associated with response in all three samples. Functional study of rs7905446 showed greater activity of the G allele in regulating expression of HTR7. The G allele displayed higher luciferase activity in two neuronal-related cell lines, and estrogen treatment decreased the activity of only the G allele. Electrophoretic mobility shift assay suggested that the G allele interacted with CCAAT/enhancer-binding protein beta transcription factor (TF), while the T allele did not show any interaction with any TFs. Our results provided novel pharmacogenomic evidence to support the role of HTR7 in association with antidepressant response