High sensitivity EndoV mutation scanning through real-time ligase proofreading

The ability to associate mutations in cancer genes with the disease and its subtypes is critical for understanding oncogenesis and identifying biomarkers for clinical diagnosis. A two-step mutation scanning method that sequentially used endonuclease V (EndoV) to nick at mismatches and DNA ligase to reseal incorrectly or nonspecifically nicked sites was previously developed in our laboratory. Herein we report an optimized single-step assay that enables ligase to proofread EndoV cleavage in real-time under a compromise between buffer conditions. Real-time proofreading results in a dramatic reduction of background cleavage. A universal PCR strategy that employs both unlabeled gene-specific primers and labeled universal primers, allows for multiplexed gene amplification and precludes amplification of primer dimers. Internally labeled PCR primers eliminate EndoV cleavage at the 5′ terminus, enabling high-throughput capillary electrophoresis readout. Furthermore, signal intensity is increased and artifacts are reduced by generating heteroduplexes containing only one of the two possible mismatches (e.g. either A/C or G/T). The single-step assay improves sensitivity to 1:50 and 1:100 (mutant:wild type) for unknown mutations in the p53 and K-ras genes, respectively, opening prospects as an early detection tool

EndoV/DNA ligase mutation scanning assay using microchip capillary electrophoresis and dual-color laser-induced fluorescence detection

We report the ability to detect with high sensitivity sporadic mutations using a mutation scanning assay, which employs thermostable endonuclease V (EndoV) and DNA ligase. The products of the mutation scanning assay were separated using microchip capillary electrophoresis (mu CE) and detected with a dual-color laser-induced fluorescence (LIF) detector. PCR products from mutant and wild-type DNA of p53 exon 8 were generated using Cy3-labeled forward and Cy5-labeled reverse primers to allow LIF detection with mCE. EndoV recognizes and primarily cleaves heteroduplexed DNA one base 30 to a mismatch and can nick matched sites at low levels as well. DNA ligase is used to reseal nicks generated at matched sites, which creates a highly sensitive and specific assay for analyzing sporadic mutations in genomic DNA. Heteroduplexed DNA samples were treated with EndoV alone and with both EndoV and DNA ligase and separated using a 4% (w/v) linear polyacrylamide gel constituted in 1x TTE buffer, 7 M urea, and 0.05% (w/v) methyl hydroxyethyl cellulose, which was used to suppress the EOF in the microchip. Sizing of the bands appearing in the electropherogram revealed the approximate position of the mutation. In this study, mutations present in p53 exon 8 generated Cy3-labeled cleavage products of 158 nt and Cy5-labeled cleavage products of 195 nt. The DNA fragments were simultaneously monitored at their respective color using a dual-color LIF system with the 158 and 195 nt fragments detected along with heteroduplexed fragments of 350 nt. The microchip separation was completed within 7 min, almost tenfold shorter time compared to conventional capillary gel electrophoresis.close7

Protocol for high-throughput DNA methylation profiling in rat tissues using automated reduced representation bisulfite sequencing

Summary: Although reduced representation bisulfite sequencing (RRBS) measures DNA methylation (DNAme) across CpG-rich genomic regions with high sensitivity, the assay can be time-consuming and prone to batch effects. Here, we present a high-throughput, automated RRBS protocol starting with DNA extraction from frozen rat tissues. We describe steps for RRBS library preparation, library quality control, and sequencing. We also detail an optimized pipeline for sequencing data processing. This protocol has been applied successfully to DNAme profiling across multiple rat tissues.For complete details on the use and execution of this protocol, please refer to Nair et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics

Modeling and high-throughput experimental data uncover the mechanisms underlying Fshb gene sensitivity to gonadotropin-releasing hormone pulse frequency

Neuroendocrine control of reproduction by brain-secreted pulses of gonadotropin-releasing hormone (GnRH) represents a longstanding puzzle about extracellular signal decoding mechanisms. GnRH regulates the pituitary gonadotropin's follicle-stimulating hormone (FSH) and luteinizing hormone (LH), both of which are heterodimers specified by unique β subunits (FSHβ/LHβ). Contrary to Lhb, Fshb gene induction has a preference for low-frequency GnRH pulses. To clarify the underlying regulatory mechanisms, we developed three biologically anchored mathematical models: 1) parallel activation of Fshb inhibitory factors (e.g. inhibin α and VGF nerve growth factor-inducible), 2) activation of a signaling component with a refractory period (e.g. G protein), and 3) inactivation of a factor needed for Fshb induction (e.g. growth differentiation factor 9). Simulations with all three models recapitulated the Fshb expression levels obtained in pituitary gonadotrope cells perifused with varying GnRH pulse frequencies. Notably, simulations altering average concentration, pulse duration, and pulse frequency revealed that the apparent frequency-dependent pattern of Fshb expression in model 1 actually resulted from variations in average GnRH concentration. In contrast, models 2 and 3 showed "true" pulse frequency sensing. To resolve which components of this GnRH signal induce Fshb, we developed a high-throughput parallel experimental system. We analyzed over 4,000 samples in experiments with varying near-physiological GnRH concentrations and pulse patterns. Whereas Egr1 and Fos genes responded only to variations in average GnRH concentration, Fshb levels were sensitive to both average concentration and true pulse frequency. These results provide a foundation for understanding the role of multiple regulatory factors in modulating Fshb gene activity

Regulatory Architecture of the LβT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone

The LβT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LβT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LβT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LβT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at −40 to −90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study

Association of survival and disease progression with chromosomal instability: A genomic exploration of colorectal cancer

During disease progression the cells that comprise solid malignancies undergo significant changes in gene copy number and chromosome structure. Colorectal cancer provides an excellent model to study this process. To indentify and characterize chromosomal abnormalities in colorectal cancer, we performed a statistical analysis of 299 expression and 130 SNP arrays profiled at different stages of the disease, including normal tissue, adenoma, stages 1–4 adenocarcinoma, and metastasis. We identified broad (> 1/2 chromosomal arm) and focal (< 1/2 chromosomal arm) events. Broad amplifications were noted on chromosomes 7, 8q, 13q, 20, and X and broad deletions on chromosomes 4, 8p, 14q, 15q, 17p, 18, 20p, and 22q. Focal events (gains or losses) were identified in regions containing known cancer pathway genes, such as VEGFA, MYC, MET, FGF6, FGF23, LYN, MMP9, MYBL2, AURKA, UBE2C, and PTEN. Other focal events encompassed potential new candidate tumor suppressors (losses) and oncogenes (gains), including CCDC68, CSMD1, POLR1D, and PMEPA1. From the expression data, we identified genes whose expression levels reflected their copy number changes and used this relationship to impute copy number changes to samples without accompanying SNP data. This analysis provided the statistical power to show that deletions of 8p, 4p, and 15q are associated with survival and disease progression, and that samples with simultaneous deletions in 18q, 8p, 4p, and 15q have a particularly poor prognosis. Annotation analysis reveals that the oxidative phosphorylation pathway shows a strong tendency for decreased expression in the samples characterized by poor prognosis.Gilbert Family FoundationHilton-Ludwig Cancer Metastasis InitiativeRidgefield FoundationNational Cancer Institut