Whole-Exome Sequencing (WES) for Illumina Short Read Sequencers Using Solution-Based Capture.

Next-generation sequencing (NGS) is transforming clinical research and diagnostics, vastly enhancing our ability to identify novel disease-causing genetic mutations and perform comprehensive diagnostic testing in the clinic. Whole-exome sequencing (WES) is a commonly used method which captures the majority of coding regions of the genome for sequencing, as these regions contain the majority of disease-causing mutations. The clinical applications of WES are not limited to diagnosis; the technique can be employed to help determine an optimal therapeutic strategy for a patient considering their mutation profile. WES may also be used to predict a patient's risk of developing a disease, e.g., type 2 diabetes (T2D), and can therefore be used to tailor advice for the patient about lifestyle choices that could mitigate those risks. Thus, genome sequencing strategies, such as WES, underpin the emerging field of personalized medicine. Initiatives also exist for sharing WES data in public repositories, e.g., the Exome Aggregation Consortium (ExAC) database. In time, by mining these valuable data resources, we will acquire a better understanding of the roles of both single rare mutations and specific combinations of common mutations (mutation signatures) in the pathology of complex diseases such as diabetes.Herein, we describe a protocol for performing WES on genomic DNA extracted from blood or saliva. Starting with gDNA extraction, we document preparation of a library for sequencing on Illumina instruments and the enrichment of the protein-coding regions from the library using the Roche NimbleGen SeqCap EZ Exome v3 kit; a solution-based capture method. We include details of how to efficiently purify the products of each step using the AMPure XP System and describe how to use qPCR to test the efficiency of capture, and thus determine finished library quality.

Gene expression elucidates functional impact of polygenic risk for schizophrenia.

Over 100 genetic loci harbor schizophrenia-associated variants, yet how these variants confer liability is uncertain. The CommonMind Consortium sequenced RNA from dorsolateral prefrontal cortex of people with schizophrenia (N = 258) and control subjects (N = 279), creating a resource of gene expression and its genetic regulation. Using this resource, ∼20% of schizophrenia loci have variants that could contribute to altered gene expression and liability. In five loci, only a single gene was involved: FURIN, TSNARE1, CNTN4, CLCN3 or SNAP91. Altering expression of FURIN, TSNARE1 or CNTN4 changed neurodevelopment in zebrafish; knockdown of FURIN in human neural progenitor cells yielded abnormal migration. Of 693 genes showing significant case-versus-control differential expression, their fold changes were ≤ 1.33, and an independent cohort yielded similar results. Gene co-expression implicates a network relevant for schizophrenia. Our findings show that schizophrenia is polygenic and highlight the utility of this resource for mechanistic interpretations of genetic liability for brain diseases.

Informed decision-making among students analyzing their personal genomes on a whole genome sequencing course: a longitudinal cohort study.

BACKGROUND: Multiple laboratories now offer clinical whole genome sequencing (WGS). We anticipate WGS becoming routinely used in research and clinical practice. Many institutions are exploring how best to educate geneticists and other professionals about WGS. Providing students in WGS courses with the option to analyze their own genome sequence is one strategy that might enhance students' engagement and motivation to learn about personal genomics. However, if this option is presented to students, it is vital they make informed decisions, do not feel pressured into analyzing their own genomes by their course directors or peers, and feel free to analyze a third-party genome if they prefer. We therefore developed a 26-hour introductory genomics course in part to help students make informed decisions about whether to receive personal WGS data in a subsequent advanced genomics course. In the advanced course, they had the option to receive their own personal genome data, or an anonymous genome, at no financial cost to them. Our primary aims were to examine whether students made informed decisions regarding analyzing their personal genomes, and whether there was evidence that the introductory course enabled the students to make a more informed decision. METHODS: This was a longitudinal cohort study in which students (N = 19) completed questionnaires assessing their intentions, informed decision-making, attitudes and knowledge before (T1) and after (T2) the introductory course, and before the advanced course (T3). Informed decision-making was assessed using the Decisional Conflict Scale. RESULTS: At the start of the introductory course (T1), most (17/19) students intended to receive their personal WGS data in the subsequent course, but many expressed conflict around this decision. Decisional conflict decreased after the introductory course (T2) indicating there was an increase in informed decision-making, and did not change before the advanced course (T3). This suggests that it was the introductory course content rather than simply time passing that had the effect. In the advanced course, all (19/19) students opted to receive their personal WGS data. No changes in technical knowledge of genomics were observed. Overall attitudes towards WGS were broadly positive. CONCLUSIONS: Providing students with intensive introductory education about WGS may help them make informed decisions about whether or not to work with their personal WGS data in an educational setting.

Genome-wide mapping of methylated adenine residues in pathogenic Escherichia coli using single-molecule real-time sequencing.

Single-molecule real-time (SMRT) DNA sequencing allows the systematic detection of chemical modifications such as methylation but has not previously been applied on a genome-wide scale. We used this approach to detect 49,311 putative 6-methyladenine (m6A) residues and 1,407 putative 5-methylcytosine (m5C) residues in the genome of a pathogenic Escherichia coli strain. We obtained strand-specific information for methylation sites and a quantitative assessment of the frequency of methylation at each modified position. We deduced the sequence motifs recognized by the methyltransferase enzymes present in this strain without prior knowledge of their specificity. Furthermore, we found that deletion of a phage-encoded methyltransferase-endonuclease (restriction-modification; RM) system induced global transcriptional changes and led to gene amplification, suggesting that the role of RM systems extends beyond protecting host genomes from foreign DNA.

Genetic diagnosis of neuroacanthocytosis disorders using exome sequencing.

Neuroacanthocytoses are neurodegenerative disorders marked by phenotypic and genetic heterogeneity. There are several associated genetic loci, and many defects, including gene deletions and insertions, and missense, nonsense, and splicing mutations, have been found spread over hundreds of kilobases of genomic DNA. In some cases, specific diagnosis is unclear, particularly in the early stages of disease or when there is an atypical presentation. Determination of the precise genetic defect allows assignment of the diagnosis and permits carrier detection and genetic counseling. The objective of this report was to utilize exome sequencing for genetic diagnosis in the neuroacanthocytosis syndromes. Genomic DNA from 2 patients with clinical features of chorea-acanthocytosis was subjected to targeted exon capture. Captured DNA was subjected to ultrahigh throughput next-generation sequencing. Sequencing data were assembled, filtered against known human variant genetic databases, and results were analyzed. Both patients were compound heterozygotes for mutations in the VPS13A gene, the gene associated with chorea-acanthocytosis. Patient 1 had a 4-bp deletion that removes the 5' donor splice site of exon 58 and a nucleotide substitution that disrupts the 5' donor splice site of exon 70. Patient 2 had a dinucleotide deletion in exon 16 and a dinucleotide insertion in exon 33. No mutations were identified in the XK, PANK2, or JPH3 gene loci. Exome sequencing is a valuable diagnostic tool in the neuroacanthocytosis syndromes. These studies may provide a better understanding of the function of the associated proteins and provide insight into the pathogenesis of these disorders.

A large gene network in immature erythroid cells is controlled by the myeloid and B cell transcriptional regulator PU.1.

PU.1 is a hematopoietic transcription factor that is required for the development of myeloid and B cells. PU.1 is also expressed in erythroid progenitors, where it blocks erythroid differentiation by binding to and inhibiting the main erythroid promoting factor, GATA-1. However, other mechanisms by which PU.1 affects the fate of erythroid progenitors have not been thoroughly explored. Here, we used ChIP-Seq analysis for PU.1 and gene expression profiling in erythroid cells to show that PU.1 regulates an extensive network of genes that constitute major pathways for controlling growth and survival of immature erythroid cells. By analyzing fetal liver erythroid progenitors from mice with low PU.1 expression, we also show that the earliest erythroid committed cells are dramatically reduced in vivo. Furthermore, we find that PU.1 also regulates many of the same genes and pathways in other blood cells, leading us to propose that PU.1 is a multifaceted factor with overlapping, as well as distinct, functions in several hematopoietic lineages.

Parallel declines in cognition, motivation, and locomotion in aging mice: association with immune gene upregulation in the medial prefrontal cortex.

Aging in humans is associated with parallel changes in cognition, motivation, and motoric performance. Based on the human aging literature, we hypothesized that this constellation of age-related changes is mediated by the medial prefrontal cortex and that it would be observed in aging mice. Toward this end, we performed detailed assessments of cognition, motivation, and motoric behavior in aging mice. We assessed behavioral and cognitive performance in C57Bl/6 mice aged 6, 18, and 24 months, and followed this with microarray analysis of tissue from the medial prefrontal cortex and analysis of serum cytokine levels. Multivariate modeling of these data suggested that the age-related changes in cognition, motivation, motor performance, and prefrontal immune gene expression were highly correlated. Peripheral cytokine levels were also correlated with these variables, but less strongly than measures of prefrontal immune gene upregulation. To determine whether the observed immune gene expression changes were due to prefrontal microglial cells, we isolated CD11b-positive cells from the prefrontal cortex and subject them to next-generation RNA sequencing. Many of the immune changes present in whole medial prefrontal cortex were enriched in this cell population. These data suggest that, as in humans, cognition, motivation, and motoric performance in the mouse change together with age and are strongly associated with CNS immune gene upregulation.

A multiprotein complex necessary for both transcription and DNA replication at the ß-globin locus.

DNA replication, repair, transcription and chromatin structure are intricately associated nuclear processes, but the molecular links between these events are often obscure. In this study, we have surveyed the protein complexes that bind at ß-globin locus control region, and purified and characterized the function of one such multiprotein complex from human erythroleukemic K562 cells. We further validated the existence of this complex in human CD34+ cell-derived normal erythroid cells. This complex contains ILF2/ILF3 transcription factors, p300 acetyltransferase and proteins associated with DNA replication, transcription and repair. RNAi knockdown of ILF2, a DNA-binding component of this complex, abrogates the recruitment of the complex to its cognate DNA sequence and inhibits transcription, histone acetylation and usage of the origin of DNA replication at the ß-globin locus. These results imply a direct link between mammalian DNA replication, transcription and histone acetylation mediated by a single multiprotein complex.

X chromosome-wide analyses of genomic DNA methylation states and gene expression in male and female neutrophils.

The DNA methylation status of human X chromosomes from male and female neutrophils was identified by high-throughput sequencing of HpaII and MspI digested fragments. In the intergenic and intragenic regions on the X chromosome, the sites outside CpG islands were heavily hypermethylated to the same degree in both genders. Nearly half of X chromosome promoters were either hypomethylated or hypermethylated in both females and males. Nearly one third of X chromosome promoters were a mixture of hypomethylated and heterogeneously methylated sites in females and were hypomethylated in males. Thus, a large fraction of genes that are silenced on the inactive X chromosome are hypomethylated in their promoter regions. These genes frequently belong to the evolutionarily younger strata of the X chromosome. The promoters that were hypomethylated at more than two sites contained most of the genes that escaped silencing on the inactive X chromosome. The overall levels of expression of X-linked genes were indistinguishable in females and males, regardless of the methylation state of the inactive X chromosome. Thus, in addition to DNA methylation, other factors are involved in the fine tuning of gene dosage compensation in neutrophils.

Chromatin architecture and transcription factor binding regulate expression of erythrocyte membrane protein genes.

Erythrocyte membrane protein genes serve as excellent models of complex gene locus structure and function, but their study has been complicated by both their large size and their complexity. To begin to understand the intricate interplay of transcription, dynamic chromatin architecture, transcription factor binding, and genomic organization in regulation of erythrocyte membrane protein genes, we performed chromatin immunoprecipitation (ChIP) coupled with microarray analysis and ChIP coupled with massively parallel DNA sequencing in both erythroid and nonerythroid cells. Unexpectedly, most regions of GATA-1 and NF-E2 binding were remote from gene promoters and transcriptional start sites, located primarily in introns. Cooccupancy with FOG-1, SCL, and MTA-2 was found at all regions of GATA-1 binding, with cooccupancy of SCL and MTA-2 also found at regions of NF-E2 binding. Cooccupancy of GATA-1 and NF-E2 was found frequently. A common signature of histone H3 trimethylation at lysine 4, GATA-1, NF-E2, FOG-1, SCL, and MTA-2 binding and consensus GATA-1-E-box binding motifs located 34 to 90 bp away from NF-E2 binding motifs was found frequently in erythroid cell-expressed genes. These results provide insights into our understanding of membrane protein gene regulation in erythropoiesis and the regulation of complex genetic loci in erythroid and nonerythroid cells and identify numerous candidate regions for mutations associated with membrane-linked hemolytic anemia.

Dynamics of alpha-globin locus chromatin structure and gene expression during erythroid differentiation of human CD34(+) cells in culture.

OBJECTIVE: The aim of the present study has been to establish serum-free culture conditions for ex vivo expansion and differentiation of human CD34(+) cells into erythroid lineage and to study the chromatin structure, gene expression, and transcription factor recruitment at the alpha-globin locus in the developing erythron. MATERIALS AND METHODS: A basal Iscove's modified Dulbecco's medium cell culture medium with 1% bovine serum albumin as a serum replacement and a combination of cytokines and growth factors was used for expansion and differentiation of the CD34(+) cells. Expression patterns of the alpha- and beta-like genes at various stages of erythropoiesis was studied by reverse transcriptase quantitative polymerase chain reaction analysis, profile of key erythroid transcription factors was investigated by Western blotting, and the chromatin structure and transcription factor recruitment at the alpha-globin locus was investigated by chromatin immunoprecipitation quantitative polymerase chain reaction analysis. RESULTS: Human CD34(+) cells in the serum-free medium undergo near synchronous erythroid differentiation to yield large amount of cells at different differentiation stages. We observe distinct patterns of the histone modifications and transcription factor binding at the alpha-globin locus during erythroid differentiation of CD34(+) cells. Nuclear factor erythroid-derived 2 (NF-E2) was present at upstream activator sites even before addition of erythropoietin (EPO), while bound GATA-1 was only detectable after EPO treatment. After 7 days of EPO treatment, H3K4Me2 modification uniformly increases throughout the alpha-globin locus. Acetylation at H3K9 and binding of Pol II, NF-E2, and GATA-1 were restricted to certain hypersensitive sites of the enhancer and theta gene, and were conspicuously low at the alpha-like globin promoters. Rearrangement of the insulator binding factor CTCF took place at and around the alpha-globin locus as CD34(+) cells differentiated into erythroid pathway. CONCLUSION: Our results indicate that remodeling of the upstream elements may be the primary event in activation of alpha-globin gene expression. Activation of alpha-globin genes upon EPO treatment involves initial binding of Pol II, downregulation of pre-existing factors like NF-E2, removal of CTCF from the locus, then rebinding of CTCF in an altered pattern, and concurrent or subsequent binding of transcription factors like GATA-1.

A genomic analysis of RNA polymerase II modification and chromatin architecture related to 3' end RNA polyadenylation.

Genomic analyses have been applied extensively to analyze the process of transcription initiation in mammalian cells, but less to transcript 3' end formation and transcription termination. We used a novel approach to prepare 3' end fragments from polyadenylated RNA, and mapped the position of the poly(A) addition site using oligonucleotide arrays tiling 1% of the human genome. This approach revealed more 3' ends than had been annotated. The distribution of these ends relative to RNA polymerase II (PolII) and di- and trimethylated lysine 4 and lysine 36 of histone H3 was compared. A substantial fraction of unannotated 3' ends of RNA are intronic and antisense to the embedding gene. Poly(A) ends of annotated messages lie on average 2 kb upstream of the end of PolII binding (termination). Near the termination sites, and in some internal sites, unphosphorylated and C-terminal domain (CTD) serine 2 phosphorylated PolII (POLR2A) accumulate, suggesting pausing of the polymerase and perhaps dephosphorylation prior to release. Lysine 36 trimethylation occurs across transcribed genes, sometimes alternating with stretches of DNA in which lysine 36 dimethylation is more prominent. Lysine 36 methylation decreases at or near the site of polyadenylation, sometimes disappearing before disappearance of phosphorylated RNA PolII or release of PolII from DNA. Our results suggest that transcription termination loss of histone 3 lysine 36 methylation and later release of RNA polymerase. The latter is often associated with polymerase pausing. Overall, our study reveals extensive sites of poly(A) addition and provides insights into the events that occur during 3' end formation.

Control of beta globin genes.

The developmental changes in expression of the beta like genes from embryonic to adult stages of human life are controlled at least partially at the level of the promoter sequences of these genes and their binding factors, and competition for promoter specific interactions with the locus control region (LCR). In recent years, the control of beta globin genes has also been investigated at the level of chromatin structure involving the chemical modification of histones and their remodelling by DNA dependent ATPases (SMARCA) containing protein complexes. The role of intergenic RNA is also being investigated with renewed interest. Although a wealth of information on the structure/function relationship of the LCR and globin promoters has been gathered over more than two decades, the fundamental nature of the control of these genes at the molecular level is still not completely understood. In the following pages, we intend to briefly describe the progress made in the field and discuss future directions.

A limited number of genes are involved in the differentiation of germinal center B cells.

Mature B cells, upon activation, progressively differentiate through centroblasts into centrocytes and finally to plasmacytes that express large amounts of selected immunoglobulins. A significant part of this maturation is thought to involve induction of the unfolded protein response (UPR). We have compared gene expression in normal germinal center centroblasts, centrocytes, lymphoblastoid cells undergoing induced UPR, and the CCL155 plasmacytoma cell line. In the centroblast to centrocyte transition there is a change in the expression of a relatively small number of genes. These include a limited subset of the genes upregulated by a fully activated UPR as well as a small number of other transcription factors, some disulphide isomerases, and other genes. This is consistent with a model in which this transition is mediated by changes in the levels of expression of transcription factor B-lymphocyte-induced maturation protein 1 (Blimp1) (PRDM1), BACH2, X-box binding protein 1 (XBP1), interferon regulatory factor 4 (IRF4), and possibly vitamin D receptor (VDR) expression, together with post-transcriptional changes and a limited induction of aspects of the UPR.

Multi-protein complexes at the beta-globin locus.

Biochemical analysis of the beta-globin gene function has led to the identification of several multi-protein complexes at the locus control region (LCR), insulator and promoters. This review briefly summarizes these multi-protein complexes and discusses their contribution towards the regulation of the beta-globin gene expression.

Heterogeneous nuclear ribonucleoprotein C1/C2, MeCP1, and SWI/SNF form a chromatin remodeling complex at the beta-globin locus control region.

Locus control regions (LCRs) are regulatory DNA sequences that are situated many kilobases away from their cognate promoters. LCRs protect transgenes from position effect variegation and heterochromatinization and also promote copy-number dependence of the levels of transgene expression. In this work, we describe the biochemical purification of a previously undescribed LCR-associated remodeling complex (LARC) that consists of heterogeneous nuclear ribonucleoprotein C1/C2, nucleosome remodeling SWI/SNF, and nucleosome remodeling and deacetylating (NuRD)/MeCP1 as a single homogeneous complex. LARC binds to the hypersensitive 2 (HS2)-Maf recognition element (MARE) DNA in a sequence-specific manner and remodels nucleosomes. Heterogeneous nuclear ribonucleoprotein C1/C2, previously known as a general RNA binding protein, provides a sequence-specific DNA recognition element for LARC, and the LARC DNA-recognition sequence is essential for the enhancement of transcription by HS2. Independently of the initiation of transcription, LARC becomes associated with beta-like globin promoters.

GATA-1 binding sites mapped in the beta-globin locus by using mammalian chIp-chip analysis.

The expression of the beta-like globin genes is intricately regulated by a series of both general and tissue-restricted transcription factors. The hemapoietic lineage-specific transcription factor GATA-1 is important for erythroid differentiation and has been implicated in regulating the expression of the erythroid-specific genes including the genes of the beta-globin locus. In the human erythroleukemic K562 cell line, only one DNA region has been identified previously as a putative site of GATA-1 interaction by in vivo footprinting studies. We mapped GATA-1 binding throughout the beta-globin locus by using chIp-chip analysis of K562 cells. We found that GATA-1 binds in a region encompassing the HS2 core element, as was previously identified, and an additional region of GATA-1 binding upstream of the gammaG gene. This approach will be of general utility for mapping transcription factor binding sites within the beta-globin locus and throughout the genome.

DNA-dependent adenosine triphosphatase (helicaselike transcription factor) activates beta-globin transcription in K562 cells.

Correct developmental regulation of beta-like globin gene expression is achieved by preferential transcription of a gene at a given developmental stage, silencing of other beta-like gene promoters, and competition among these promoters for interaction with the locus control region (LCR). Several evolutionarily conserved DNA elements in the promoters of the beta-like genes and LCR have been studied in detail, and the role of their binding factors has been investigated. However, the beta-globin promoter includes additional evolutionarily conserved sequences of unknown function. The present study examined the properties of a 21-base pair (bp) promoter-conserved sequence (PCS) located at positions -115 to -136 bp relative to the transcription start site of the beta-globin gene. A helicaselike transcription factor (HLTF) belonging to the SWI2/SNF2 family of proteins binds to the PCS and a partly homologous sequence in the enhancer region of the LCR hypersensitive site 2 (HS2). Elevation of the level of HLTF in K562 erythroleukemic cells increases beta-promoter activity in transient transfection experiments, and mutations in the PCS that remove HLTF-binding regions abolish this effect, suggesting that HLTF is an activator of beta-globin transcription. Overexpression of HLTF in K562 cells does not affect the endogenous levels of gamma- and epsilon-globin message, but it markedly activates beta-globin transcription. In conclusion, this study reports a transcription factor belonging to the SWI2/SNF2 family, which preferentially activates chromosomal beta-globin gene transcription and which has not previously been implicated in globin gene regulation.