Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6.

The cytochrome P450-2D6 (CYP2D6) enzyme metabolizes ∼25% of common medications, yet homologous pseudogenes and copy number variants (CNVs) make interrogating the polymorphic CYP2D6 gene with short-read sequencing challenging. Therefore, we developed a novel long-read, full gene CYP2D6 single molecule real-time (SMRT) sequencing method using the Pacific Biosciences platform. Long-range PCR and CYP2D6 SMRT sequencing of 10 previously genotyped controls identified expected star (*) alleles, but also enabled suballele resolution, diplotype refinement, and discovery of novel alleles. Coupled with an optimized variant-calling pipeline, CYP2D6 SMRT sequencing was highly reproducible as triplicate intra- and inter-run nonreference genotype results were completely concordant. Importantly, targeted SMRT sequencing of upstream and downstream CYP2D6 gene copies characterized the duplicated allele in 15 control samples with CYP2D6 CNVs. The utility of CYP2D6 SMRT sequencing was further underscored by identifying the diplotypes of 14 samples with discordant or unclear CYP2D6 configurations from previous targeted genotyping, which again included suballele resolution, duplicated allele characterization, and discovery of a novel allele and tandem arrangement. Taken together, long-read CYP2D6 SMRT sequencing is an innovative, reproducible, and validated method for full-gene characterization, duplication allele-specific analysis, and novel allele discovery, which will likely improve CYP2D6 metabolizer phenotype prediction for both research and clinical testing applications.

Interrogation of a context-specific transcription factor network identifies novel regulators of pluripotency.

The predominant view of pluripotency regulation proposes a stable ground state with coordinated expression of key transcription factors (TFs) that prohibit differentiation. Another perspective suggests a more complexly regulated state involving competition between multiple lineage-specifying TFs that define pluripotency. These contrasting views were developed from extensive analyses of TFs in pluripotent cells in vitro. An experimentally validated, genome-wide repertoire of the regulatory interactions that control pluripotency within the in vivo cellular contexts is yet to be developed. To address this limitation, we assembled a TF interactome of adult human male germ cell tumors (GCTs) using the Algorithm for the Accurate Reconstruction of Cellular Pathways (ARACNe) to analyze gene expression profiles of 141 tumors comprising pluripotent and differentiated subsets. The network (GCT(Net)) comprised 1,305 TFs, and its ingenuity pathway analysis identified pluripotency and embryonal development as the top functional pathways. We experimentally validated GCT(Net) by functional (silencing) and biochemical (ChIP-seq) analysis of the core pluripotency regulatory TFs POU5F1, NANOG, and SOX2 in relation to their targets predicted by ARACNe. To define the extent of the in vivo pluripotency network in this system, we ranked all TFs in the GCT(Net) according to sharing of ARACNe-predicted targets with those of POU5F1 and NANOG using an odds-ratio analysis method. To validate this network, we silenced the top 10 TFs in the network in H9 embryonic stem cells. Silencing of each led to downregulation of pluripotency and induction of lineage; 7 of the 10 TFs were identified as pluripotency regulators for the first time.

Direct ChIP-Seq significance analysis improves target prediction.

BACKGROUND: Chromatin immunoprecipitation followed by sequencing of protein-bound DNA fragments (ChIP-Seq) is an effective high-throughput methodology for the identification of context specific DNA fragments that are bound by specific proteins in vivo. Despite significant progress in the bioinformatics analysis of this genome-scale data, a number of challenges remain as technology-dependent biases, including variable target accessibility and mappability, sequence-dependent variability, and non-specific binding affinity must be accounted for. RESULTS AND DISCUSSION: We introduce a nonparametric method for scoring consensus regions of aligned immunoprecipitated DNA fragments when appropriate control experiments are available. Our method uses local models for null binding; these are necessary because binding prediction scores based on global models alone fail to properly account for specialized features of genomic regions and chance pull downs of specific DNA fragments, thus disproportionally rewarding some genomic regions and decreasing prediction accuracy. We make no assumptions about the structure or amplitude of bound peaks, yet we show that our method outperforms leading methods developed using either global or local null hypothesis models for random binding. We test prediction performance by comparing analyses of ChIP-seq, ChIP-chip, motif-based binding-site prediction, and shRNA assays, showing high reproducibility, binding-site enrichment in predicted target regions, and functional regulation of predicted targets. CONCLUSIONS: Given appropriate controls, a direct nonparametric method for identifying transcription-factor targets from ChIP-Seq assays may lead to both higher sensitivity and higher specificity, and should be preferred or used in conjunction with methods that use parametric models for null binding.

Deletion of ERF and CIC causes abnormal skull morphology and global developmental delay.

The ETS2 repressor factor (ERF) is a transcription factor in the RAS-MEK-ERK signal transduction cascade that regulates cell proliferation and differentiation, and pathogenic sequence variants in the ERF gene cause variable craniosynostosis inherited in an autosomal dominant pattern. The reported ERF variants are largely loss-of-function, implying haploinsufficiency as a primary disease mechanism; however, ERF gene deletions have not been reported previously. Here we describe three probands with macrocephaly, craniofacial dysmorphology, and global developmental delay. Clinical genetic testing for fragile X and other relevant sequencing panels were negative; however, chromosomal microarray identified heterozygous deletions (63.7-583.2 kb) on Chromosome 19q13.2 in each proband that together included five genes associated with Mendelian diseases (ATP1A3, ERF, CIC, MEGF8, and LIPE). Parental testing indicated that the aberrations were apparently de novo in two of the probands and were inherited in the one proband with the smallest deletion. Deletion of ERF is consistent with the reported loss-of-function ERF variants, prompting clinical copy-number-variant classifications of likely pathogenic. Moreover, the recent characterization of heterozygous loss-of-function CIC sequence variants as a cause of intellectual disability and neurodevelopmental disorders inherited in an autosomal dominant pattern is also consistent with the developmental delays and intellectual disabilities identified among the two probands with CIC deletions. Taken together, this case series adds to the previously reported patients with ERF and/or CIC sequence variants and supports haploinsufficiency of both genes as a mechanism for a variable syndromic cranial phenotype with developmental delays and intellectual disability inherited in an autosomal dominant pattern.

Cooperative binding of the yeast Spt10p activator to the histone upstream activating sequences is mediated through an N-terminal dimerization domain.

The yeast Spt10p activator is a putative histone acetyltransferase (HAT) possessing a sequence-specific DNA-binding domain (DBD) which binds to the upstream activation sequences (UAS elements) in the histone gene promoters. Spt10p binds to a pair of histone UAS elements with extreme positive cooperativity. The molecular basis of this cooperativity was addressed. Spt10p (640 residues) is an elongated dimer, but the isolated DBD (residues 283-396) is a monomer and binds non-cooperatively to DNA. A Spt10p fragment comprising the N-terminal domain (NTD), HAT domain and DBD (residues 1-396) binds cooperatively and is a dimer, whereas an overlapping Spt10p fragment comprising the DBD and C-terminal domains (residues 283-640) binds non-cooperatively and is a monomer. These observations imply that cooperative binding requires dimerization. The isolated NTD (residues 1-98) is a dimer and is responsible for dimerization. We propose that cooperativity involves a conformational change in the Spt10p dimer which facilitates the simultaneous recognition of two UAS elements. In vivo, deletion of the NTD results in poor growth, but does not prevent the binding at the HTA1 promoter, suggesting that dimerization is biologically important. Residues 1-396 are sufficient for normal growth, indicating that the critical functions of Spt10p reside in the N-terminal domains.

Integrated CYP2D6 interrogation for multiethnic copy number and tandem allele detection.

AIM: To comprehensively interrogate CYP2D6 by integrating genotyping, copy number analysis and novel strategies to identify CYP2D6*36 and characterize CYP2D6 duplications. METHODS: Genotyping of 16 CYP2D6 alleles, multiplex ligation-dependent probe amplification (MLPA) and CYP2D6*36 and duplication allele-specific genotyping were performed on 427 African-American, Asian, Caucasian, Hispanic, and Ashkenazi Jewish individuals. RESULTS: A novel PCR strategy determined that almost half of all CYP2D6*10 (100C>T) alleles are actually *36 (isolated or in tandem with *10) and all identified duplication alleles were characterized. Integrated results from all testing platforms enabled the refinement of genotype frequencies across all studied populations. CONCLUSION: The polymorphic CYP2D6 gene requires comprehensive interrogation to characterize allelic variation across ethnicities, which was enabled in this study by integrating multiplexed genotyping, MLPA copy number analysis, novel PCR strategies and duplication allele-specific genotyping.

Global regulation by the yeast Spt10 protein is mediated through chromatin structure and the histone upstream activating sequence elements.

The yeast SPT10 gene encodes a putative histone acetyltransferase (HAT) implicated as a global transcription regulator acting through basal promoters. Here we address the mechanism of this global regulation. Although microarray analysis confirmed that Spt10p is a global regulator, Spt10p was not detected at any of the most strongly affected genes in vivo. In contrast, the presence of Spt10p at the core histone gene promoters in vivo was confirmed. Since Spt10p activates the core histone genes, a shortage of histones could occur in spt10Delta cells, resulting in defective chromatin structure and a consequent activation of basal promoters. Consistent with this hypothesis, the spt10Delta phenotype can be rescued by extra copies of the histone genes and chromatin is poorly assembled in spt10Delta cells, as shown by irregular nucleosome spacing and reduced negative supercoiling of the endogenous 2mum plasmid. Furthermore, Spt10p binds specifically and highly cooperatively to pairs of upstream activating sequence elements in the core histone promoters [consensus sequence, (G/A)TTCCN(6)TTCNC], consistent with a direct role in histone gene regulation. No other high-affinity sites are predicted in the yeast genome. Thus, Spt10p is a sequence-specific activator of the histone genes, possessing a DNA-binding domain fused to a likely HAT domain.

Dual Diagnosis of Ellis-van Creveld Syndrome and Hearing Loss in a Consanguineous Family.

Multilocus analysis of rare or genetically heterogeneous diseases is a distinct advantage of next-generation sequencing (NGS) over conventional single-gene investigations. Recent studies have begun to uncover an under-recognized prevalence of dual molecular diagnoses in patients with a "blended" phenotype that is the result of 2 clinical diagnoses involving 2 separate genetic loci. This blended phenotype could be mistakenly interpreted as a novel clinical extension of a single-gene disorder. In this study, we ascertained a proband from a large consanguineous Iranian family who manifests postlingual, progressive, moderate hearing loss in addition to suspected Ellis-van Creveld syndrome phenotype. NGS with a customized skeletal dysplasia panel containing over 370 genes and subsequent bioinformatics analysis disclosed 2 homozygous mutations in EVC2 (c.2653C>T; p.Arg885*) and COL11A2 (c.966dup; p.Thr323Hisfs*19), respectively. This study highlights a dual molecular diagnosis in a patient with a blending of 2 distinct phenotypes and illustrates the advantage and importance of this staple technology to facilitate rapid and comprehensive genetic dissection of a heterogeneous phenotype. The differentiation between phenotypic expansion of a genetic disorder and a blended phenotype that is due to more than one distinct genetic aberration is essential in order to reduce the diagnostic odyssey endured by patients.

Cross-platform assessment of genomic imbalance confirms the clinical relevance of genomic complexity and reveals loci with potential pathogenic roles in diffuse large B-cell lymphoma.

Genomic copy number alterations (CNAs) in diffuse large B-cell lymphoma (DLBCL) have roles in disease pathogenesis, but overall clinical relevance remains unclear. Herein, an unbiased algorithm was uniformly applied across three genome profiling datasets comprising 392 newly-diagnosed DLBCL specimens that defined 32 overlapping CNAs, involving 36 minimal common regions (MCRs). Scoring criteria were established for 50 aberrations within the MCRs while considering peak gains/losses. Application of these criteria to independent datasets revealed novel candidate genes with coordinated expression, such as CNOT2, potentially with pathogenic roles. No one single aberration significantly associated with patient outcome across datasets, but genomic complexity, defined by imbalance in more than one MCR, significantly portended adverse outcome in two of three independent datasets. Thus, the standardized scoring of CNAs currently developed can be uniformly applied across platforms, affording robust validation of genomic imbalance and complexity in DLBCL and overall clinical utility as biomarkers of patient outcome.

Mechanism and Role of SOX2 Repression in Seminoma: Relevance to Human Germline Specification.

Human male germ cell tumors (GCTs) are derived from primordial germ cells (PGCs). The master pluripotency regulator and neuroectodermal lineage effector transcription factor SOX2 is repressed in PGCs and the seminoma (SEM) subset of GCTs. The mechanism of SOX2 repression and its significance to GC and GCT development currently are not understood. Here, we show that SOX2 repression in SEM-derived TCam-2 cells is mediated by the Polycomb repressive complex (PcG) and the repressive H3K27me3 chromatin mark that are enriched at its promoter. Furthermore, SOX2 repression in TCam-2 cells can be abrogated by recruitment of the constitutively expressed H3K27 demethylase UTX to the SOX2 promoter through retinoid signaling, leading to expression of neuronal and other lineage genes. SOX17 has been shown to initiate human PGC specification, with its target PRDM1 suppressing mesendodermal genes. Our results are consistent with a role for SOX2 repression in normal germline development by suppressing neuroectodermal genes.

Genomic imbalance defines three prognostic groups for risk stratification of patients with chronic lymphocytic leukemia.

Array comparative genomic hybridization (aCGH) has yet to be fully leveraged in a prognostic setting in chronic lymphocytic leukemia (CLL). Genomic imbalance was assessed in 288 CLL specimens using a targeted array. Based on 20 aberrations in a hierarchical manner, all 228 treatment-naive specimens were classified into a group with poor outcome (20.6%) exhibiting at least one aberration that was univariately associated with adverse outcome (gain: 2p, 3q, 8q, 17q, loss: 7q, 8p, 11q, 17p, 18p), good outcome (32.5%) showing 13q14 loss without any of the other 10 aberrations (gain: 1p, 7p, 12, 18p, 18q, 19, loss: 4p, 5p, 6q, 7p) or intermediate outcome (remainder). The three groups were significantly separated with respect to time to first treatment and overall survival (p < 0.001), and validation of the stratification scheme was performed in two independent datasets. Gain of 3q and 8q, and 17p loss were determined to be independent unfavorable prognostic biomarkers. TP53, NOTCH1 and SF3B1 mutations correlated with the presence of one poor outcome aCGH marker, at a considerably higher frequency than when only considering poor risk aberrations routinely detected by fluorescence in situ hybridization (FISH). These data support genomic imbalance evaluation in CLL by aCGH to assist in risk stratification.

The DNA-binding domain of the yeast Spt10p activator includes a zinc finger that is homologous to foamy virus integrase.

The yeast SPT10 gene encodes a putative histone acetyltransferase that binds specifically to pairs of upstream activating sequence (UAS) elements found only in the histone gene promoters. Here, we demonstrate that the DNA-binding domain of Spt10p is located between residues 283 and 396 and includes a His(2)-Cys(2) zinc finger. The binding of Spt10p to the histone UAS is zinc-dependent and is disabled by a zinc finger mutation (C388S). The isolated DNA-binding domain binds to single histone UAS elements with high affinity. In contrast, full-length Spt10p binds with high affinity only to pairs of UAS elements with very strong positive cooperativity and is unable to bind to a single UAS element. This implies the presence of a "blocking" domain in full-length Spt10p, which forces it to search for a pair of UAS elements. Chromatin immunoprecipitation experiments indicate that, unlike wild-type Spt10p, the C388S protein does not bind to the promoter of the gene encoding histone H2A (HTA1) in vivo. The C388S mutant has a phenotype similar to that of the spt10Delta mutant: poor growth and global aberrations in gene expression. Thus, the C388S mutation disables the DNA-binding function of Spt10p in vitro and in vivo. The zinc finger of Spt10p is homologous to that of foamy virus integrase, perhaps suggesting that this integrase is also a sequence-specific DNA-binding protein.