Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences

Resolving the three-dimensional interactome of Human Accelerated Regions during human and chimpanzee neurodevelopment

Human Accelerated Regions (HARs) are highly conserved across species but exhibit a significant excess of human-specific sequence changes, suggesting they may have gained novel functions in human evolution. HARs include transcriptional enhancers with human-specific activity and have been implicated in the evolution of the human brain. However, our understanding of how HARs contributed to uniquely human features of the brain is hindered by a lack of insight into the genes and pathways that HARs regulate. It is unclear whether HARs acted by altering the expression of gene targets conserved between HARs and their chimpanzee orthologs or by gaining new gene targets in human, a mechanism termed enhancer hijacking. We generated a high-resolution map of chromatin interactions for 1,590 HARs and their orthologs in human and chimpanzee neural stem cells (NSCs) to comprehensively identify gene targets in both species. HARs and their chimpanzee orthologs targeted a conserved set of 2,963 genes enriched for neurodevelopmental processes including neurogenesis and synaptic transmission. Changes in HAR enhancer activity were correlated with changes in conserved gene target expression. Conserved targets were enriched among genes differentially expressed between human and chimpanzee NSCs or between human and non-human primate developing and adult brain. Species-specific HAR gene targets did not converge on known biological functions and were not significantly enriched among differentially expressed genes, suggesting that HARs did not alter gene expression via enhancer hijacking. HAR gene targets, including differentially expressed targets, also showed cell type-specific expression patterns in the developing human brain, including outer radial glia, which are hypothesized to contribute to human cortical expansion. Our findings support that HARs influenced human brain evolution by altering the expression of conserved gene targets and provide the means to functionally link HARs with novel human brain features.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Phenotypic screening models for rapid diagnosis of genetic variants and discovery of personalized therapeutics

Precision medicine strives for highly individualized treatments for disease under the notion that each individual's unique genetic makeup and environmental exposures imprints upon them not only a disposition to illness, but also an optimal therapeutic approach. In the realm of rare disorders, genetic predisposition is often the predominant mechanism driving disease presentation. For such, mostly, monogenic disorders, a causal gene to phenotype association is likely. As a result, it becomes important to query the patient's genome for the presence of pathogenic variations that are likely to cause the disease. Determining whether a variant is pathogenic or not is critical to these analyses and can be challenging, as many disease-causing variants are novel and, ergo, have no available functional data to help categorize them. This problem is exacerbated by the need for rapid evaluation of pathogenicity, since many genetic diseases present in young children who will experience increased morbidity and mortality without rapid diagnosis and therapeutics. Here, we discuss the utility of animal models, with a focus mainly on C. elegans, as a contrast to tissue culture and in silico approaches, with emphasis on how these systems are used in determining pathogenicity of variants with uncertain significance and then used to screen for novel therapeutics.

Alignment-Free Sequence Comparison With Multiple k Values

Alignment-free sequence comparison approaches have become increasingly popular in computational biology, because alignment-based approaches are inefficient to process large-scale datasets. Still, there is no way to determine the optimal value of the critical parameter k for alignment-free approaches in general. In this article, we tried to solve the problem by involving multiple k values simultaneously. The method counts the occurrence of each k-mer with different k values in a sequence. Two weighting schemes, based on maximizing deviation method and genetic algorithm, are then used on these counts. We applied the method to enhance the three common alignment-free approaches D₂, D₂^S, and D₂^*, and evaluated its performance on similarity search and functionally related regulatory sequences recognition. The enhanced approaches achieve better performance than the original approaches in all cases, and much better performance than some other common measures, such as Pcc, Eu, Ma, Ch, Kld, and Cos.

Accurate prediction of cis-regulatory modules reveals a prevalent regulatory genome of humans

cis-regulatory modules(CRMs) formed by clusters of transcription factor (TF) binding sites (TFBSs) are as important as coding sequences in specifying phenotypes of humans. It is essential to categorize all CRMs and constituent TFBSs in the genome. In contrast to most existing methods that predict CRMs in specific cell types using epigenetic marks, we predict a largely cell type agonistic but more comprehensive map of CRMs and constituent TFBSs in the gnome by integrating all available TF ChIP-seq datasets. Our method is able to partition 77.47% of genome regions covered by available 6092 datasets into a CRM candidate (CRMC) set (56.84%) and a non-CRMC set (43.16%). Intriguingly, the predicted CRMCs are under strong evolutionary constraints, while the non-CRMCs are largely selectively neutral, strongly suggesting that the CRMCs are likely cis-regulatory, while the non-CRMCs are not. Our predicted CRMs are under stronger evolutionary constraints than three state-of-the-art predictions (GeneHancer, EnhancerAtlas and ENCODE phase 3) and substantially outperform them for recalling VISTA enhancers and non-coding ClinVar variants. We estimated that the human genome might encode about 1.47M CRMs and 68M TFBSs, comprising about 55% and 22% of the genome, respectively; for both of which, we predicted 80%. Therefore, the cis-regulatory genome appears to be more prevalent than originally thought.

The Impact of the FKBP5 Gene Polymorphisms on the Relationship between Traumatic Life Events and Psychotic-Like Experiences in Non-Clinical Adults

Common variations of the FKBP5 gene are implicated in psychotic disorders, by modulating the hypothalamic–pituitary–adrenal axis reactivity to stress. It has been demonstrated that some of them might moderate the effects of childhood trauma on psychosis proneness. However, these associations have not been investigated with respect to traumatic life events (TLEs). Therefore, we aimed to explore whether the FKBP5 polymorphisms moderate the effects of TLEs on the level of psychotic-like experiences (PLEs). A total of 535 non-clinical adults were approached for participation, and genotyping of six FKBP5 polymorphisms (rs3800373, rs9470080, rs4713902, rs737054, rs1360780 and rs9296158) was performed. The Prodromal Questionnaire-16 (PQ-16) and the Traumatic Events Checklist (TEC) were administered to assess PLEs and TLEs, respectively. Among the rs1360780 CC homozygotes, a history of physical abuse was associated with significantly higher PQ-16 scores. This difference was not significant in the rs1360780 T allele carriers. Similarly, a history of physical abuse was associated with significantly higher PQ-16 scores in the rs9296158 GG homozygotes but not in the rs9296158 A allele carriers. Finally, emotional neglect was related to significantly higher PQ-16 scores in the rs737054 T allele carriers but not in the rs737054 CC homozygotes. The present study indicates that variation in the FKBP5 gene might moderate the effects of lifetime traumatic events on psychosis proneness.

MicroRNA-298 reduces levels of human amyloid-β precursor protein (APP), β-site APP-converting enzyme 1 (BACE1) and specific tau protein moieties

Alzheimer's disease (AD) is the most common age-related form of dementia, associated with deposition of intracellular neuronal tangles consisting primarily of hyperphosphorylated microtubule-associated protein tau (p-tau) and extracellular plaques primarily comprising amyloid- β (Aβ) peptide. The p-tau tangle unit is a posttranslational modification of normal tau protein. Aβ is a neurotoxic peptide excised from the amyloid-β precursor protein (APP) by β-site APP-cleaving enzyme 1 (BACE1) and the γ-secretase complex. MicroRNAs (miRNAs) are short, single-stranded RNAs that modulate protein expression as part of the RNA-induced silencing complex (RISC). We identified miR-298 as a repressor of APP, BACE1, and the two primary forms of Aβ (Aβ40 and Aβ42) in a primary human cell culture model. Further, we discovered a novel effect of miR-298 on posttranslational levels of two specific tau moieties. Notably, miR-298 significantly reduced levels of ~55 and 50 kDa forms of the tau protein without significant alterations of total tau or other forms. In vivo overexpression of human miR-298 resulted in nonsignificant reduction of APP, BACE1, and tau in mice. Moreover, we identified two miR-298 SNPs associated with higher cerebrospinal fluid (CSF) p-tau and lower CSF Aβ42 levels in a cohort of human AD patients. Finally, levels of miR-298 varied in postmortem human temporal lobe between AD patients and age-matched non-AD controls. Our results suggest that miR-298 may be a suitable target for AD therapy.

A Unique Epigenomic Landscape Defines Human Erythropoiesis

Mammalian erythropoiesis yields a highly specialized cell type, the mature erythrocyte, evolved to meet the organismal needs of increased oxygen-carrying capacity. To better understand the regulation of erythropoiesis, we performed genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomics using a recently developed strategy to obtain highly purified populations of primary human erythroid cells. The integration of gene expression, DNA methylation, and chromatin state dynamics reveals that stage-specific gene regulation during erythropoiesis is a stepwise and hierarchical process involving many cis-regulatory elements. Erythroid-specific, nonpromoter sites of chromatin accessibility are linked to erythroid cell phenotypic variation and inherited disease. Comparative analyses of stage-specific chromatin accessibility indicate that there is limited early chromatin priming of erythroid genes during hematopoiesis. The epigenome of terminally differentiating erythroid cells defines a distinct subset of highly specialized cells that are vastly dissimilar from other hematopoietic and nonhematopoietic cell types. These epigenomic and transcriptome data are powerful tools to study human erythropoiesis.

Schulz et al. use genome-wide studies of chromatin accessibility, DNA methylation, and transcriptomes in primary human erythroid cells to reveal important characteristics of erythropoiesis. Chromatin accessibility of terminal erythroid differentiation is markedly dissimilar from other hematopoietic cell types. Epigenomic changes are linked to erythroid cell traits and disease genes.

Creating New β-Globin-Expressing Lentiviral Vectors by High-Resolution Mapping of Locus Control Region Enhancer Sequences

Hematopoietic stem cell gene therapy is a promising approach for treating disorders of the hematopoietic system. Identifying combinations of cis-regulatory elements that do not impede packaging or transduction efficiency when included in lentiviral vectors has proven challenging. In this study, we deploy LV-MPRA (lentiviral vector-based, massively parallel reporter assay), an approach that simultaneously analyzes thousands of synthetic DNA fragments in parallel to identify sequence-intrinsic and lineage-specific enhancer function at near-base-pair resolution. We demonstrate the power of LV-MPRA in elucidating the boundaries of previously unknown intrinsic enhancer sequences of the human β-globin locus control region. Our approach facilitated the rapid assembly of novel therapeutic βAS3-globin lentiviral vectors harboring strong lineage-specific recombinant control elements capable of correcting a mouse model of sickle cell disease. LV-MPRA can be used to map any genomic locus for enhancer activity and facilitates the rapid development of therapeutic vectors for treating disorders of the hematopoietic system or other specific tissues and cell types.

KTN1 variants and risk for attention deficit hyperactivity disorder

Individuals with attention deficit hyperactivity disorder (ADHD) show gray matter volume (GMV) reduction in the putamen. KTN1 variants may regulate kinectin 1 expression in the putamen and influence putamen structure and function. We aim to test the hypothesis that the KTN1 variants may represent a genetic risk factor of ADHD. Two independent family-based Caucasian samples were analyzed, including 922 parent-child trios (a total of 2,757 subjects with 924 ADHD children) and 735 parent-child trios (a total of 1,383 subjects with 613 ADHD children). The association between ADHD and a total of 143 KTN1 SNPs was analyzed in the first sample, and the nominally-significant (p < .05) risk SNPs were classified into independent haplotype blocks. All SNPs, including imputed SNPs within these blocks, and haplotypes across each block, were explored for replication of associations in both samples. The potential biological functions of all risk SNPs were predicted using a series of bioinformatics analyses, their regulatory effects on the putamen volumes were tested, and the KTN1 mRNA expression was examined in three independent human putamen tissue samples. We found that fifteen SNPs were nominally associated with ADHD (p < .05) in the first sample, and three of them remained significant even after correction for multiple testing (1.3 × 10^-10  ≤ p ≤ 1.2 × 10^-4 ; α = 2.5 × 10^-3 ). These 15 risk SNPs were located in five haplotype blocks, and 13 SNPs within four of these blocks were associated with ADHD in the second sample. Six haplotypes within these blocks were also significantly (1.2 × 10^-7  ≤ p ≤ .009) associated with ADHD in these samples. These risk variants were located in disease-related transposons and/or transcription-related functional regions. Major alleles of these risk variants significantly increased putamen volumes. Finally, KTN1 mRNA was significantly expressed in putamen across three independent cohorts. We concluded that the KTN1 variants were significantly associated with ADHD. KTN1 may play a functional role in the development of ADHD.

Population genetic models of GERP scores suggest pervasive turnover of constrained sites across mammalian evolution

Comparative genomic approaches have been used to identify sites where mutations are under purifying selection and of functional consequence by searching for sequences that are conserved across distantly related species. However, the performance of these approaches has not been rigorously evaluated under population genetic models. Further, short-lived functional elements may not leave a footprint of sequence conservation across many species. We use simulations to study how one measure of conservation, the Genomic Evolutionary Rate Profiling (GERP) score, relates to the strength of selection (N_es). We show that the GERP score is related to the strength of purifying selection. However, changes in selection coefficients or functional elements over time (i.e. functional turnover) can strongly affect the GERP distribution, leading to unexpected relationships between GERP and N_es. Further, we show that for functional elements that have a high turnover rate, adding more species to the analysis does not necessarily increase statistical power. Finally, we use the distribution of GERP scores across the human genome to compare models with and without turnover of sites where mutations are under purifying selection. We show that mutations in 4.51% of the noncoding human genome are under purifying selection and that most of this sequence has likely experienced changes in selection coefficients throughout mammalian evolution. Our work reveals limitations to using comparative genomic approaches to identify deleterious mutations. Commonly used GERP score thresholds miss over half of the noncoding sites in the human genome where mutations are under purifying selection.

One of the most significant and challenging tasks in modern genomics is to assess the functional consequences of a particular nucleotide change in a genome. A common approach to address this challenge prioritizes sequences that share similar nucleotides across distantly related species, with the rationale that mutations at such positions were deleterious and removed from the population by purifying natural selection. Our manuscript shows that one popular measure of sequence conservation, the GERP score, performs well at identifying selected mutations if mutations at a site were under selection across all of mammalian evolution. Changes in selection at a given site dramatically reduces the power of GERP to detect selected mutations in humans. We also combine population genetic models with the distribution of GERP scores at noncoding sites across the human genome to show that the degree of selection at individual sites has changed throughout mammalian evolution. Importantly, we demonstrate that at least 80 Mb of noncoding sequence under purifying selection in humans will not have extreme GERP scores and will likely be missed by modern comparative genomic approaches. Our work argues that new approaches, potentially based on genetic variation within species, will be required to identify deleterious mutations.

A likelihood approach to testing hypotheses on the co-evolution of epigenome and genome

Central questions to epigenome evolution include whether interspecies changes of histone modifications are independent of evolutionary changes of DNA, and if there is dependence whether they depend on any specific types of DNA sequence changes. Here, we present a likelihood approach for testing hypotheses on the co-evolution of genome and histone modifications. The gist of this approach is to convert evolutionary biology hypotheses into probabilistic forms, by explicitly expressing the joint probability of multispecies DNA sequences and histone modifications, which we refer to as a class of Joint Evolutionary Model for the Genome and the Epigenome (JEMGE). JEMGE can be summarized as a mixture model of four components representing four evolutionary hypotheses, namely dependence and independence of interspecies epigenomic variations to underlying sequence substitutions and to underlying sequence insertions and deletions (indels). We implemented a maximum likelihood method to fit the models to the data. Based on comparison of likelihoods, we inferred whether interspecies epigenomic variations depended on substitution or indels in local genomic sequences based on DNase hypersensitivity and spermatid H3K4me3 ChIP-seq data from human and rhesus macaque. Approximately 5.5% of homologous regions in the genomes exhibited H3K4me3 modification in either species, among which approximately 67% homologous regions exhibited local-sequence-dependent interspecies H3K4me3 variations. Substitutions accounted for less local-sequence-dependent H3K4me3 variations than indels. Among transposon-mediated indels, ERV1 insertions and L1 insertions were most strongly associated with H3K4me3 gains and losses, respectively. By initiating probabilistic formulation on the co-evolution of genomes and epigenomes, JEMGE helps to bring evolutionary biology principles to comparative epigenomic studies.

Epigenetic modifications play a significant role in gene regulations and thus heavily influence phenotypic outcomes. Whereas cross-species epigenomic comparisons have been fruitful in revealing the function of epigenetic modifications, it still remains unclear how the epigenome changes across species. A central question in epigenome evolution studies is whether interspecies epigenomic variations rely on genomic changes in cis and, if partially yes, whether different genomic changes have distinct impacts. To tackle this question, we initiated a likelihood-based approach, in which different hypotheses related to the co-evolution of the genome and the epigenome could be converted into probabilistic models. By fitting the models to actual data, each model yielded a likelihood, and the hypothesis corresponded to the largest likelihood was selected as most supported by observed data. In this work, we focused on the influence of two types of underlying sequence changes: substitutions, and insertions and deletions (indels). We quantitatively assessed the dependence of H3K4me3 variations on substitutions and indels between human and rhesus, and separated their relative impacts within each genomic region with H3K4me3. The methodology presented here provides a framework for modeling the epigenome together with the genome and a quantitative approach to test different evolutionary hypotheses.

Impact of genetic variations in ADORA2A gene on depression and symptoms: a cross-sectional population-based study

Genetic variants involved in adenosine metabolism and its receptors were associated with increased risk for psychiatric disorders, including anxiety, depression, and schizophrenia. Here, we examined an association between a single nucleotide polymorphism in A_2A receptor gene (ADORA2A, rs2298383 SNP) with current depressive episode and symptom profile. A total of 1253 individuals from a cross-sectional population-based study were analyzed by the Mini International Neuropsychiatric Interview 5.0. Our data showed that the TT genotype of ADORA2A rs2298383 SNP was associated with reduced risk for depression when compared to the CC/CT genotypes (p = 0.020). This association remained significant after adjusting for confounding variables such as smoking, gender, socioeconomic class, and ethnicity (OR = 0.631 (95% CI 0.425-0.937); p = 0.022). Regarding the symptoms associated with depression, we evaluated the impact of the ADORA2A SNP in the occurrence of sad/discouraged mood, anhedonia, appetite changes, sleep disturbances, motion changes, energy loss, feelings of worthless or guilty, difficulty in concentrating, and presence of bad thoughts. Notably, the TT genotype was independently associated with reduced sleep disturbances (OR = 0.438 (95% CI 0.258-0.743); p = 0.002) and less difficulty in concentrating (OR = 0.534 (95% CI 0.316-0.901; p = 0.019). The cross-sectional design cannot evaluate the cause-effect relationship and did not evaluate the functional consequences of this polymorphism. Our data support an important role for ADORA2A rs2298383 SNP in clinical heterogeneity associated with depression. The presence of the TT genotype was associated with decrease risk for current depression and disturbances in sleep and attention, two of the most common symptoms associated with this disorder.

Single-Base Resolution Map of Evolutionary Constraints and Annotation of Conserved Elements across Major Grass Genomes

Conserved noncoding sequences (CNSs) are evolutionarily conserved DNA sequences that do not encode proteins but may have potential regulatory roles in gene expression. CNS in crop genomes could be linked to many important agronomic traits and ecological adaptations. Compared with the relatively mature exon annotation protocols, efficient methods are lacking to predict the location of noncoding sequences in the plant genomes. We implemented a computational pipeline that is tailored to the comparisons of plant genomes, yielding a large number of conserved sequences using rice genome as the reference. In this study, we used 17 published grass genomes, along with five monocot genomes as well as the basal angiosperm genome of Amborella trichopoda. Genome alignments among these genomes suggest that at least 12.05% of the rice genome appears to be evolving under constraints in the Poaceae lineage, with close to half of the evolutionarily constrained sequences located outside protein-coding regions. We found evidence for purifying selection acting on the conserved sequences by analyzing segregating SNPs within the rice population. Furthermore, we found that known functional motifs were significantly enriched within CNS, with many motifs associated with the preferred binding of ubiquitous transcription factors. The conserved elements that we have curated are accessible through our public database and the JBrowse server. In-depth functional annotations and evolutionary dynamics of the identified conserved sequences provide a solid foundation for studying gene regulation, genome evolution, as well as to inform gene isolation for cereal biologists.

Integrative genomic analysis of methylphenidate response in attention-deficit/hyperactivity disorder

Methylphenidate (MPH) is the most frequently used pharmacological treatment in children with attention-deficit/hyperactivity disorder (ADHD). However, a considerable interindividual variability exists in clinical outcome. Thus, we performed a genome-wide association study of MPH efficacy in 173 ADHD paediatric patients. Although no variant reached genome-wide significance, the set of genes containing single-nucleotide polymorphisms (SNPs) nominally associated with MPH response (P < 0.05) was significantly enriched for candidates previously studied in ADHD or treatment outcome. We prioritised the nominally significant SNPs by functional annotation and expression quantitative trait loci (eQTL) analysis in human brain, and we identified 33 SNPs tagging cis-eQTL in 32 different loci (referred to as eSNPs and eGenes, respectively). Pathway enrichment analyses revealed an over-representation of genes involved in nervous system development and function among the eGenes. Categories related to neurological diseases, psychological disorders and behaviour were also significantly enriched. We subsequently meta-analysed the association with clinical outcome for the 33 eSNPs across the discovery sample and an independent cohort of 189 ADHD adult patients (target sample) and we detected 15 suggestive signals. Following this comprehensive strategy, our results provide a better understanding of the molecular mechanisms implicated in MPH treatment effects and suggest promising candidates that may encourage future studies.

Transcriptome analysis of developing lens reveals abundance of novel transcripts and extensive splicing alterations

Lens development involves a complex and highly orchestrated regulatory program. Here, we investigate the transcriptomic alterations and splicing events during mouse lens formation using RNA-seq data from multiple developmental stages, and construct a molecular portrait of known and novel transcripts. We show that the extent of novelty of expressed transcripts decreases significantly in post-natal lens compared to embryonic stages. Characterization of novel transcripts into partially novel transcripts (PNTs) and completely novel transcripts (CNTs) (novelty score ≥ 70%) revealed that the PNTs are both highly conserved across vertebrates and highly expressed across multiple stages. Functional analysis of PNTs revealed their widespread role in lens developmental processes while hundreds of CNTs were found to be widely expressed and predicted to encode for proteins. We verified the expression of four CNTs across stages. Examination of splice isoforms revealed skipped exon and retained intron to be the most abundant alternative splicing events during lens development. We validated by RT-PCR and Sanger sequencing, the predicted splice isoforms of several genes Banf1, Cdk4, Cryaa, Eif4g2, Pax6, and Rbm5. Finally, we present a splicing browser Eye Splicer (http://www.iupui.edu/~sysbio/eye-splicer/), to facilitate exploration of developmentally altered splicing events and to improve understanding of post-transcriptional regulatory networks during mouse lens development.

Associations of rare nicotinic cholinergic receptor gene variants to nicotine and alcohol dependence

Nicotine's rewarding effects are mediated through distinct subunits of nAChRs, encoded by different nicotinic cholinergic receptor (CHRN) genes and expressed in discrete regions in the brain. In the present study, we aimed to test the associations between rare variants at CHRN genes and nicotine dependence (ND), and alcohol dependence (AD). A total of 26,498 subjects with nine different neuropsychiatric disorders in 15 independent cohorts, which were genotyped on Illumina, Affymetrix, or PERLEGEN microarray platforms, were analyzed. Associations between rare variants (minor allele frequency (MAF) <0.05) at CHRN genes and nicotine dependence, and alcohol dependence were tested. The mRNA expression of all Chrn genes in whole mouse brain and 10 specific brain areas was investigated. All CHRN genes except the muscle-type CHRNB1, including eight genomic regions containing 11 neuronal CHRN genes and three genomic regions containing four muscle-type CHRN genes, were significantly associated with ND, and/or AD. All of these genes were expressed in the mouse brain. We conclude that CHRNs are associated with ND (mainly) and AD, supporting the hypothesis that the full catalog of ND/AD risk genes may contain most neuronal nAChRs-encoding genes. © 2016 Wiley Periodicals, Inc.

Replicated Risk Nicotinic Cholinergic Receptor Genes for Nicotine Dependence

It has been hypothesized that the nicotinic acetylcholine receptors (nAChRs) play important roles in nicotine dependence (ND) and influence the number of cigarettes smoked per day (CPD) in smokers. We compiled the associations between nicotinic cholinergic receptor genes (CHRNs) and ND/CPD that were replicated across different studies, reviewed the expression of these risk genes in human/mouse brains, and verified their expression using independent samples of both human and mouse brains. The potential functions of the replicated risk variants were examined using cis-eQTL analysis or predicted using a series of bioinformatics analyses. We found replicated and significant associations for ND/CPD at 19 SNPs in six genes in three genomic regions (CHRNB3-A6, CHRNA5-A3-B4 and CHRNA4). These six risk genes are expressed in at least 18 distinct areas of the human/mouse brain, with verification in our independent human and mouse brain samples. The risk variants might influence the transcription, expression and splicing of the risk genes, alter RNA secondary or protein structure. We conclude that the replicated associations between CHRNB3-A6, CHRNA5-A3-B4,CHRNA4 and ND/CPD are very robust. More research is needed to examine how these genetic variants contribute to the risk for ND/CPD.

GAPDH-mediated posttranscriptional regulations of sodium channel Scn1a and Scn3a genes under seizure and ketogenic diet conditions

Abnormal expressions of sodium channel SCN1A and SCN3A genes alter neural excitability that are believed to contribute to the pathogenesis of epilepsy, a long-term risk of recurrent seizures. Ketogenic diet (KD), a high-fat and low-carbohydrate treatment for difficult-to-control (refractory) epilepsy in children, has been suggested to reverse gene expression patterns. Here, we reveal a novel role of GAPDH on the posttranscriptional regulation of mouse Scn1a and Scn3a expressions under seizure and KD conditions. We show that GAPDH binds to a conserved region in the 3' UTRs of human and mouse SCN1A and SCN3A genes, which decreases and increases genes' expressions by affecting mRNA stability through SCN1A 3' UTR and SCN3A 3' UTR, respectively. In seizure mice, the upregulation and phosphorylation of GAPDH enhance its binding to the 3' UTR, which lead to downregulation of Scn1a and upregulation of Scn3a. Furthermore, administration of KD generates β-hydroxybutyric acid which rescues the abnormal expressions of Scn1a and Scn3a by weakening the GAPDH's binding to the element. Taken together, these data suggest that GAPDH-mediated expression regulation of sodium channel genes may be associated with epilepsy and the anticonvulsant action of KD.

Association study of polymorphisms in FOXO3, AKT1 and IGF-2R genes with human longevity in a Han Chinese population

FOXO3, AKT1 and IGF-2R are critical members of the insulin/IGF-1 signaling pathway. Previous studies showed that polymorphisms (SNPs) in FOXO3, AKT1 and IGF-2R were associated with human longevity in Caucasian population. However, the association of these SNPs in different ethnic groups is often inconsistent. Here, we investigated the association of genetic variants in three genes with human longevity in Han Chinese population. Twelve SNPs from FOXO3, AKT1 and IGF-2R were selected and genotyped in 1202 long-lived individuals (nonagenarians and centenarians) and younger individuals. Rs9486902 of FOXO3 was found to be associated with human longevity in both genders combined in this study (allelic P = 0.002, corrected P = 0.024). The other eleven SNPs were not significantly associated with human longevity in Han Chinese population. The haplotypes TTCTT, CCTTC and CTCCT of FOXO3 as well as GGTCGG and GGTCAG of AKT1 were shown to have a significant difference between case and control (P =0.006, 2.78×10-5, 4.68×10-6, 0.003,0.005, respectively). The estimated prevalence of diabetes and prediabetes in long-lived individuals was significantly lower than in common adult populations (P = 0.001, 2.3×10-26) .Therefore, the search for longevity-associated genes provides the identification of new potential targets beneficial for the treatment of diabetes.

Influence of ADORA2A gene polymorphism on leukoencephalopathy risk in MTX-treated pediatric patients affected by hematological malignancies

BACKGROUND

Methotrexate (MTX) can lead to neurotoxicity and asymptomatic leukoencephalopathy. However, the mechanism of MTX-related leukoencephalopathy is obscure. MTX and its metabolites inhibit 5-aminoimidazole-4-carboxamide ribonucleotide formiltransferase (ATIC) and promote adenosine release. Recently, it has been reported that adenosine and its receptor are related to certain central nervous system diseases. We investigated whether adenosine pathway gene polymorphisms and clinical factors were related to MTX-related leukoencephalopathy in pediatric patients affected by hematological malignancies.

PROCEDURE

Fifty-six Japanese childhood acute lymphoblastic leukemia or lymphoma patients were investigated. Patients were evaluated by magnetic resonance imaging of the brain before maintenance therapy or stem cell transplantation. Gene polymorphisms within the adenosine pathway (ATIC, adenosine A2A receptor [ADORA2A]) and the MTX pathway (methylenetetrahydrofolate reductase [MTHFR] and ABCB1) were genotyped using TaqMan assays. Clinical data were collected by accessing the medical records. MTX-related leukoencephalopathy was evaluated by a pediatric neurologist.

RESULTS

Twenty-one (37%) of 56 patients developed MTX-related leukoencephalopathy. Four of 21 patients developed clinical neurotoxicity. The minor allele CC genotype of rs2298383 (ADORA2A) was associated with MTX-related leukoencephalopathy (P = 0.010, odds ratio = 5.81, 95% confidence interval 1.50-22.50). High cumulative dose of systemic MTX was associated with MTX-related leukoencephalopathy after adjusting for sex, ADORA2A polymorphism, and prolonged high MTX concentration (P = 0.042, odds ratio = 1.18, 95% confidence interval 1.01-1.37).

CONCLUSIONS

ADORA2A rs2298383 and high cumulative dose of systemic MTX administration were significantly associated with MTX-related leukoencephalopathy. Our results indicate that pharmacological intervention within the adenosine pathway may be both a treatment and preventative option for MTX-related leukoencephalopathy.

Integrative Tissue-Specific Functional Annotations in the Human Genome Provide Novel Insights on Many Complex Traits and Improve Signal Prioritization in Genome Wide Association Studies

Extensive efforts have been made to understand genomic function through both experimental and computational approaches, yet proper annotation still remains challenging, especially in non-coding regions. In this manuscript, we introduce GenoSkyline, an unsupervised learning framework to predict tissue-specific functional regions through integrating high-throughput epigenetic annotations. GenoSkyline successfully identified a variety of non-coding regulatory machinery including enhancers, regulatory miRNA, and hypomethylated transposable elements in extensive case studies. Integrative analysis of GenoSkyline annotations and results from genome-wide association studies (GWAS) led to novel biological insights on the etiologies of a number of human complex traits. We also explored using tissue-specific functional annotations to prioritize GWAS signals and predict relevant tissue types for each risk locus. Brain and blood-specific annotations led to better prioritization performance for schizophrenia than standard GWAS p-values and non-tissue-specific annotations. As for coronary artery disease, heart-specific functional regions was highly enriched of GWAS signals, but previously identified risk loci were found to be most functional in other tissues, suggesting a substantial proportion of still undetected heart-related loci. In summary, GenoSkyline annotations can guide genetic studies at multiple resolutions and provide valuable insights in understanding complex diseases. GenoSkyline is available at http://genocanyon.med.yale.edu/GenoSkyline.

IDIOPATHIC PULMONARY FIBROSIS IS A COMPLEX GENETIC DISORDER

Idiopathic pulmonary fibrosis (IPF) is a complex, heterogeneous genetic disorder that is associated with rare and common sequence variants in many genes (MUC5B, SFTPC, SFTPA2, RTEL1, TERT, and hTR), 11 novel loci, and multiple emerging epigenetic and transcriptional profiles. In the past 5 years, we have found that: 1) genetic risk variants play major and similar roles in the development of both familial and sporadic fibrotic idiopathic interstitial pneumonia, accounting for up to 35% of the risk of idiopathic interstitial pneumonia (a disease that was previously thought to be idiopathic); 2) a promoter variant in MUC5B rs35705950 is the strongest risk factor for the development of IIP and IPF; however, rs35705950 has a low penetrance; and 3) IPF is a complex genetic disease with 11 independent loci contributing to the development of this disease, pronounced changes in DNA methylation, and transcriptional subtypes. In aggregate, these findings suggest that IPF is a heterogeneous disease and that genetic and molecular subtypes of IPF will provide essential clues to disease pathogenesis, prognosis, treatment, and survival, all of which remain major problems in understanding and treating patients with IPF. Although the basic biological mechanisms involved in IPF are emerging, the disease is heterogeneous pathologically and the final common pathways of fibrogenesis are not well understood. These observations lead us to postulate that the etiology and severity/extent of this complex condition will best be understood through an integrated approach that accounts for inherited factors, epigenetic marks, and dynamic changes in the transcriptome.

Validation of Next-Generation Sequencing of Entire Mitochondrial Genomes and the Diversity of Mitochondrial DNA Mutations in Oral Squamous Cell Carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases.

METHODS

We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base).

RESULTS

We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases.

CONCLUSIONS

We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.

A statistical framework to predict functional non-coding regions in the human genome through integrated analysis of annotation data

Identifying functional regions in the human genome is a major goal in human genetics. Great efforts have been made to functionally annotate the human genome either through computational predictions, such as genomic conservation, or high-throughput experiments, such as the ENCODE project. These efforts have resulted in a rich collection of functional annotation data of diverse types that need to be jointly analyzed for integrated interpretation and annotation. Here we present GenoCanyon, a whole-genome annotation method that performs unsupervised statistical learning using 22 computational and experimental annotations thereby inferring the functional potential of each position in the human genome. With GenoCanyon, we are able to predict many of the known functional regions. The ability of predicting functional regions as well as its generalizable statistical framework makes GenoCanyon a unique and powerful tool for whole-genome annotation. The GenoCanyon web server is available at http://genocanyon.med.yale.edu

A novel promoter mutation (HBB: c.-75G>T) was identified as a cause of β(+)-thalassemia

We report a novel β-globin gene promoter mutation in a Chinese family identified using fluorescence resolution melting curve analysis and gene sequencing. The proband, who showed the phenotype of β-thalassemia intermedia (β-TI), was found to be a compound heterozygote for the novel mutation -25 (G>T) (HBB: c.-75G>T) and a codon 17 (HBB: c.52A>T) mutation. Moreover, conservation analysis using phyloP and phastCons indicated that the mutated base in the proband was conserved. This novel point mutation on the β-globin gene is in close proximity to the conserved ATAA sequence located at position -25 relative to the mRNA Cap site. We performed a further comparative analysis of the clinical phenotypes and hematological parameters in this pedigree and found that the father was a carrier of the novel point mutation and showed low levels of hemoglobin (Hb), mean corpuscular volume (MCV) and mean corpuscular Hb (MCH). Thus, the available evidence suggests that this novel mutation, -25, results in β(+)-thalassemia (β(+)-thal).

Body mass index associated to rs2021966 ENPP1 polymorphism increases the risk for gestational diabetes mellitus

Gestational diabetes mellitus (GDM) is a condition of impaired glucose tolerance occurring in 1-14% of all pregnancies. This wide range reflects pathological involvement of single nucleotide polymorphisms (SNPs) and maternal weight as risk factors. This study evaluated the association of genetic component and maternal factors to identify women with higher risk of developing GDM. About 240 pregnant women characterized by negative Oral Glucose Tolerance Test (-OGTT) and 38 with positive OGGT (+OGTT) were enrolled. SNPs for ENPP1, NRF1, VEGFA, CEBPA, and PIK3R1 were analyzed by SNP genotyping. An association study was performed and differences in genotype and allele frequencies between cases and controls were analyzed by χ(2) test. +OGTT was associated to high values of pre-gestational body mass index (BMI) and age. SNP for ENPP1 gene was associated to +OGTT, while genetic variants for other genes did not correlate to GDM. ENPP1 homozygous for A allele and heterozygous showed altered frequencies in +OGTT when compared with -OGTT. Association of both pre-gestational BMI and age with AA homozygous genotype increased significantly the risk to +OGTT. Our results demonstrate that correlation of age and pre-gestational BMI with homozygous for A allele increased significantly the risk of impaired glucose tolerance and GDM.

RRHP: a tag-based approach for 5-hydroxymethylcytosine mapping at single-site resolution

Current methods for genomic mapping of 5-hydroxymethylcytosine (5hmC) have been limited by either costly sequencing depth, high DNA input, or lack of single-base resolution. We present an approach called Reduced Representation 5-Hydroxymethylcytosine Profiling (RRHP) to map 5hmC sites at single-base resolution by exploiting the use of beta-glucosyltransferase to inhibit enzymatic digestion at the junction where adapters are ligated to a genomic library. Therefore, only library fragments presenting glucosylated 5hmC residues at the junction are sequenced. RRHP can detect sites with low 5hmC abundance, and when combined with RRBS data, 5-methylcytosine and 5-hydroxymethylcytosine can be compared at a specific site.

A simple repeat polymorphism in the MITF-M promoter is a key regulator of white spotting in dogs

The white spotting locus (S) in dogs is colocalized with the MITF (microphtalmia-associated transcription factor) gene. The phenotypic effects of the four S alleles range from solid colour (S) to extreme white spotting (s^w). We have investigated four candidate mutations associated with the s^w allele, a SINE insertion, a SNP at a conserved site and a simple repeat polymorphism all associated with the MITF-M promoter as well as a 12 base pair deletion in exon 1B. The variants associated with white spotting at all four loci were also found among wolves and we conclude that none of these could be a sole causal mutation, at least not for extreme white spotting. We propose that the three canine white spotting alleles are not caused by three independent mutations but represent haplotype effects due to different combinations of causal polymorphisms. The simple repeat polymorphism showed extensive diversity both in dogs and wolves, and allele-sharing was common between wolves and white spotted dogs but was non-existent between solid and spotted dogs as well as between wolves and solid dogs. This finding was unexpected as Solid is assumed to be the wild-type allele. The data indicate that the simple repeat polymorphism has been a target for selection during dog domestication and breed formation. We also evaluated the significance of the three MITF-M associated polymorphisms with a Luciferase assay, and found conclusive evidence that the simple repeat polymorphism affects promoter activity. Three alleles associated with white spotting gave consistently lower promoter activity compared with the allele associated with solid colour. We propose that the simple repeat polymorphism affects cooperativity between transcription factors binding on either flanking sides of the repeat. Thus, both genetic and functional evidence show that the simple repeat polymorphism is a key regulator of white spotting in dogs.

Deep sequencing reveals a novel class of bidirectional promoters associated with neuronal genes

Background

Comprehensive annotation of transcripts expressed in a given tissue is a critical step towards the understanding of regulatory and functional pathways that shape the transcriptome.

Results

Here, we reconstructed a cumulative transcriptome of the human prefrontal cortex (PFC) based on approximately 300 million strand-specific RNA sequence (RNA-seq) reads collected at different stages of postnatal development. We find that more than 50% of reconstructed transcripts represent novel transcriptome elements, including 8,343 novel exons and exon extensions of annotated coding genes, 11,217 novel antisense transcripts and 29,541 novel intergenic transcripts or their fragments showing canonical features of long non-coding RNAs (lncRNAs). Our analysis further led to a surprising discovery of a novel class of bidirectional promoters (NBiPs) driving divergent transcription of mRNA and novel lncRNA pairs and displaying a distinct set of sequence and epigenetic features. In contrast to known bidirectional and unidirectional promoters, NBiPs are strongly associated with genes involved in neuronal functions and regulated by neuron-associated transcription factors.

Conclusions

Taken together, our results demonstrate that large portions of the human transcriptome remain uncharacterized. The distinct sequence and epigenetic features of NBiPs, as well as their specific association with neuronal genes, further suggest existence of regulatory pathways specific to the human brain.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-457) contains supplementary material, which is available to authorized users.

Determining causality and consequence of expression quantitative trait loci

Expression quantitative trait loci (eQTLs) are currently the most abundant and systematically-surveyed class of functional consequence for genetic variation. Recent genetic studies of gene expression have identified thousands of eQTLs in diverse tissue types for the majority of human genes. Application of this large eQTL catalog provides an important resource for understanding the molecular basis of common genetic diseases. However, only now has both the availability of individuals with full genomes and corresponding advances in functional genomics provided the opportunity to dissect eQTLs to identify causal regulatory variants. Resolving the properties of such causal regulatory variants is improving understanding of the molecular mechanisms that influence traits and guiding the development of new genome-scale approaches to variant interpretation. In this review, we provide an overview of current computational and experimental methods for identifying causal regulatory variants and predicting their phenotypic consequences.

Defining functional DNA elements in the human genome

With the completion of the human genome sequence, attention turned to identifying and annotating its functional DNA elements. As a complement to genetic and comparative genomics approaches, the Encyclopedia of DNA Elements Project was launched to contribute maps of RNA transcripts, transcriptional regulator binding sites, and chromatin states in many cell types. The resulting genome-wide data reveal sites of biochemical activity with high positional resolution and cell type specificity that facilitate studies of gene regulation and interpretation of noncoding variants associated with human disease. However, the biochemically active regions cover a much larger fraction of the genome than do evolutionarily conserved regions, raising the question of whether nonconserved but biochemically active regions are truly functional. Here, we review the strengths and limitations of biochemical, evolutionary, and genetic approaches for defining functional DNA segments, potential sources for the observed differences in estimated genomic coverage, and the biological implications of these discrepancies. We also analyze the relationship between signal intensity, genomic coverage, and evolutionary conservation. Our results reinforce the principle that each approach provides complementary information and that we need to use combinations of all three to elucidate genome function in human biology and disease.

Systems biology approach to integrative comparative genomics

As more and more systems biology approaches are used to investigate the different types of biological macromolecules, increasing numbers of whole genomic studies are now available for a large array of organisms. Whether it is genomics, transcriptomics, proteomics, interactomics or metabolomics, the full complement of genomic information on all different levels can be juxtaposed between different organisms to reveal similarities or differences, and even to provide consensus models. At the intersection of comparative genomics and systems biology lies great possibility for discovery, analysis and prediction. This paper explores this nexus and the relationship from four general levels: DNA, RNA, protein and extragenomic. For each level, we provide an overview of the methods, discuss the potential challenges and survey the current research. Finally, we suggest some organizing principles and make proposals for new areas that will be important for future research.

Dense genotyping of immune-related disease regions identifies 14 new susceptibility loci for juvenile idiopathic arthritis

We used the Immunochip array to analyze 2,816 individuals with juvenile idiopathic arthritis (JIA), comprising the most common subtypes (oligoarticular and rheumatoid factor-negative polyarticular JIA), and 13,056 controls. We confirmed association of 3 known JIA risk loci (the human leukocyte antigen (HLA) region, PTPN22 and PTPN2) and identified 14 loci reaching genome-wide significance (P < 5 × 10(-8)) for the first time. Eleven additional new regions showed suggestive evidence of association with JIA (P < 1 × 10(-6)). Dense mapping of loci along with bioinformatics analysis refined the associations to one gene in each of eight regions, highlighting crucial pathways, including the interleukin (IL)-2 pathway, in JIA disease pathogenesis. The entire Immunochip content, the HLA region and the top 27 loci (P < 1 × 10(-6)) explain an estimated 18, 13 and 6% of the risk of JIA, respectively. In summary, this is the largest collection of JIA cases investigated so far and provides new insight into the genetic basis of this childhood autoimmune disease.

Machine learning and genome annotation: a match meant to be?

By its very nature, genomics produces large, high-dimensional datasets that are well suited to analysis by machine learning approaches. Here, we explain some key aspects of machine learning that make it useful for genome annotation, with illustrative examples from ENCODE.

Systematic dissection of regulatory motifs in 2000 predicted human enhancers using a massively parallel reporter assay

Genome-wide chromatin annotations have permitted the mapping of putative regulatory elements across multiple human cell types. However, their experimental dissection by directed regulatory motif disruption has remained unfeasible at the genome scale. Here, we use a massively parallel reporter assay (MPRA) to measure the transcriptional levels induced by 145-bp DNA segments centered on evolutionarily conserved regulatory motif instances within enhancer chromatin states. We select five predicted activators (HNF1, HNF4, FOXA, GATA, NFE2L2) and two predicted repressors (GFI1, ZFP161) and measure reporter expression in erythroleukemia (K562) and liver carcinoma (HepG2) cell lines. We test 2104 wild-type sequences and 3314 engineered enhancer variants containing targeted motif disruptions, each using 10 barcode tags and two replicates. The resulting data strongly confirm the enhancer activity and cell-type specificity of enhancer chromatin states, the ability of 145-bp segments to recapitulate both, the necessary role of regulatory motifs in enhancer function, and the complementary roles of activator and repressor motifs. We find statistically robust evidence that (1) disrupting the predicted activator motifs abolishes enhancer function, while silent or motif-improving changes maintain enhancer activity; (2) evolutionary conservation, nucleosome exclusion, binding of other factors, and strength of the motif match are predictive of enhancer activity; (3) scrambling repressor motifs leads to aberrant reporter expression in cell lines where the enhancers are usually inactive. Our results suggest a general strategy for deciphering cis-regulatory elements by systematic large-scale manipulation and provide quantitative enhancer activity measurements across thousands of constructs that can be mined to develop predictive models of gene expression.

Genetic variants of MARCO are associated with susceptibility to pulmonary tuberculosis in a Gambian population

BACKGROUND

The two major class A scavenger receptors are scavenger receptor A (SRA), which is constitutively expressed on most macrophage populations, and macrophage receptor with collagenous structure (MARCO), which is constitutively expressed on a more restricted subset of macrophages, (e.g. alveolar macrophages) but whose expression increases on most macrophages during the course of infection. Although the primary role of SRA appears to be clearance of modified host proteins and lipids, mice defective in expression of either MARCO or SRA are immunocompromised in multiple models of infection and in vitro assays, the scavenger receptors have been demonstrated to bind bacteria and to enhance pro-inflammatory signalling to many bacterial lung pathogens; however their importance in Mycobacterium tuberculosis infection, is less clear.

METHODS

To determine whether polymorphisms in either SRA or MARCO were associated with tuberculosis, a case-control study of was performed. DNA samples from newly-detected, smear-positive, pulmonary tuberculosis cases were collected from The Gambia. Controls for this study consisted of DNA from cord bloods obtained from routine births at local Gambian health clinics. Informed written consent was obtained from patients or their parents or guardians. Ethical approval was provided by the joint The Gambian Government/MRC Joint Ethics Committee.

RESULTS

We studied the frequencies of 25 polymorphisms of MSR1 (SRA) and 22 in MARCO in individuals with tuberculosis (n=1284) and matched controls (n=1349). No SNPs within the gene encoding or within 1 kb of the promoter sequence of MSR1 were associated with either susceptibility or resistance to tuberculosis. Three SNPs in MARCO (rs4491733, Mantel-Haenszel 2x2 χ2 = 6.5, p = 0.001, rs12998782, Mantel-Haenszel 2x2 χ2 = 6.59, p = 0.001, rs13389814 Mantel-Haenszel 2x2 χ2 = 6.9, p = 0.0009) were associated with susceptibility to tuberculosis and one (rs7559955, Mantel-Haenszel 2x2 χ2 = 6.9, p = 0.0009) was associated with resistance to tuberculosis.

CONCLUSIONS

These findings identify MARCO as a potentially important receptor in the host response to tuberculosis.

ADORA2A polymorphism predisposes children to encephalopathy with febrile status epilepticus

OBJECTIVE

Acute encephalopathy with biphasic seizures and late reduced diffusion (AESD) is a childhood encephalopathy following severe febrile seizures, leaving neurologic sequelae in many patients. However, its pathogenesis remains unclear. In this study, we clarified that genetic variation in the adenosine A2A receptor (ADORA2A), whose activation is involved in excitotoxicity, may be a predisposing factor of AESD.

METHODS

We analyzed 4 ADORA2A single nucleotide polymorphisms in 85 patients with AESD. The mRNA expression in brain samples, mRNA and protein expression in lymphoblasts, as well as the production of cyclic adenosine monophosphate (cAMP) by lymphoblasts in response to adenosine were compared among ADORA2A diplotypes.

RESULTS

Four single nucleotide polymorphisms were completely linked, which resulted in 2 haplotypes, A and B. Haplotype A (C at rs2298383, T at rs5751876, deletion at rs35320474, and C at rs4822492) frequency in patients was significantly higher than in controls (p = 0.005). Homozygous haplotype A (AA diplotype) had a higher risk of developing AESD (odds ratio 2.32, 95% confidence interval 1.32-4.08; p = 0.003) via a recessive model. mRNA expression was significantly higher in AA than AB and BB diplotypes, both in the brain (p = 0.003 and 0.002, respectively) and lymphoblasts (p = 0.035 and 0.003, respectively). In lymphoblasts, ADORA2A protein expression (p = 0.024), as well as cellular cAMP production (p = 0.0006), was significantly higher in AA than BB diplotype.

CONCLUSIONS

AA diplotype of ADORA2A is associated with AESD and may alter the intracellular adenosine/cAMP cascade, thereby promoting seizures and excitotoxic brain damage in patients.

Identification of biologically relevant enhancers in human erythroid cells

Identification of cell type-specific enhancers is important for understanding the regulation of programs controlling cellular development and differentiation. Enhancers are typically marked by the co-transcriptional activator protein p300 or by groups of cell-expressed transcription factors. We hypothesized that a unique set of enhancers regulates gene expression in human erythroid cells, a highly specialized cell type evolved to provide adequate amounts of oxygen throughout the body. Using chromatin immunoprecipitation followed by massively parallel sequencing, genome-wide maps of candidate enhancers were constructed for p300 and four transcription factors, GATA1, NF-E2, KLF1, and SCL, using primary human erythroid cells. These data were combined with gene expression analyses, and candidate enhancers were identified. Consistent with their predicted function as candidate enhancers, there was statistically significant enrichment of p300 and combinations of co-localizing erythroid transcription factors within 1-50 kb of the transcriptional start site (TSS) of genes highly expressed in erythroid cells. Candidate enhancers were also enriched near genes with known erythroid cell function or phenotype. Candidate enhancers exhibited moderate conservation with mouse and minimal conservation with nonplacental vertebrates. Candidate enhancers were mapped to a set of erythroid-associated, biologically relevant, SNPs from the genome-wide association studies (GWAS) catalogue of NHGRI, National Institutes of Health. Fourteen candidate enhancers, representing 10 genetic loci, mapped to sites associated with biologically relevant erythroid traits. Fragments from these loci directed statistically significant expression in reporter gene assays. Identification of enhancers in human erythroid cells will allow a better understanding of erythroid cell development, differentiation, structure, and function and provide insights into inherited and acquired hematologic disease.

Integrative annotation of chromatin elements from ENCODE data

The ENCODE Project has generated a wealth of experimental information mapping diverse chromatin properties in several human cell lines. Although each such data track is independently informative toward the annotation of regulatory elements, their interrelations contain much richer information for the systematic annotation of regulatory elements. To uncover these interrelations and to generate an interpretable summary of the massive datasets of the ENCODE Project, we apply unsupervised learning methodologies, converting dozens of chromatin datasets into discrete annotation maps of regulatory regions and other chromatin elements across the human genome. These methods rediscover and summarize diverse aspects of chromatin architecture, elucidate the interplay between chromatin activity and RNA transcription, and reveal that a large proportion of the genome lies in a quiescent state, even across multiple cell types. The resulting annotation of non-coding regulatory elements correlate strongly with mammalian evolutionary constraint, and provide an unbiased approach for evaluating metrics of evolutionary constraint in human. Lastly, we use the regulatory annotations to revisit previously uncharacterized disease-associated loci, resulting in focused, testable hypotheses through the lens of the chromatin landscape.

Evolution of hemoglobin and its genes

Insights into the evolution of hemoglobins and their genes are an abundant source of ideas regarding hemoglobin function and regulation of globin gene expression. This article presents the multiple genes and gene families encoding human globins, summarizes major events in the evolution of the hemoglobin gene clusters, and discusses how these studies provide insights into regulation of globin genes. Although the genes in and around the α-like globin gene complex are relatively stable, the β-like globin gene clusters are more dynamic, showing evidence of transposition to a new locus and frequent lineage-specific expansions and deletions. The cis-regulatory modules controlling levels and timing of gene expression are a mix of conserved and lineage-specific DNA, perhaps reflecting evolutionary constraint on core regulatory functions shared broadly in mammals and adaptive fine-tuning in different orders of mammals.

Nuclear RNA sequencing of the mouse erythroid cell transcriptome

In addition to protein coding genes a substantial proportion of mammalian genomes are transcribed. However, most transcriptome studies investigate steady-state mRNA levels, ignoring a considerable fraction of the transcribed genome. In addition, steady-state mRNA levels are influenced by both transcriptional and posttranscriptional mechanisms, and thus do not provide a clear picture of transcriptional output. Here, using deep sequencing of nuclear RNAs (nucRNA-Seq) in parallel with chromatin immunoprecipitation sequencing (ChIP-Seq) of active RNA polymerase II, we compared the nuclear transcriptome of mouse anemic spleen erythroid cells with polymerase occupancy on a genome-wide scale. We demonstrate that unspliced transcripts quantified by nucRNA-seq correlate with primary transcript frequencies measured by RNA FISH, but differ from steady-state mRNA levels measured by poly(A)-enriched RNA-seq. Highly expressed protein coding genes showed good correlation between RNAPII occupancy and transcriptional output; however, genome-wide we observed a poor correlation between transcriptional output and RNAPII association. This poor correlation is due to intergenic regions associated with RNAPII which correspond with transcription factor bound regulatory regions and a group of stable, nuclear-retained long non-coding transcripts. In conclusion, sequencing the nuclear transcriptome provides an opportunity to investigate the transcriptional landscape in a given cell type through quantification of unspliced primary transcripts and the identification of nuclear-retained long non-coding RNAs.

A misplaced lncRNA causes brachydactyly in humans

Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis- and in trans-acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.

PPARG binding landscapes in macrophages suggest a genome-wide contribution of PU.1 to divergent PPARG binding in human and mouse

Background

Genome-wide comparisons of transcription factor binding sites in different species can be used to evaluate evolutionary constraints that shape gene regulatory circuits and to understand how the interaction between transcription factors shapes their binding landscapes over evolution.

Results

We have compared the PPARG binding landscapes in macrophages to investigate the evolutionary impact on PPARG binding diversity in mouse and humans for this important nuclear receptor. Of note, only 5% of the PPARG binding sites were shared between the two species. In contrast, at the gene level, PPARG target genes conserved between both species constitute more than 30% of the target genes regulated by PPARG ligand in human macrophages. Moreover, the majority of all PPARG binding sites (55–60%) in macrophages show co-occupancy of the lineage-specification factor PU.1 in both species. Exploring the evolutionary dynamics of PPARG binding sites, we observed that PU.1 co-binding to PPARG sites appears to be important for possible PPARG ancestral functions such as lipid metabolism. Thus we speculate that PU.1 may have guided utilization of these species-specific PPARG conserved binding sites in macrophages during evolution.

Conclusions

We propose a model in which PU.1 sites may have served as “anchor” loci for the formation of new and functionally relevant PPARG binding sites throughout evolution. As PU.1 is an essential factor in macrophage biology, such an evolutionary mechanism would allow for the establishment of relevant PPARG regulatory modules in a PU.1-dependent manner and yet permit for nuanced regulatory changes in individual species.

Small noncoding differentially methylated copy-number variants, including lncRNA genes, cause a lethal lung developmental disorder

An unanticipated and tremendous amount of the noncoding sequence of the human genome is transcribed. Long noncoding RNAs (lncRNAs) constitute a significant fraction of non-protein-coding transcripts; however, their functions remain enigmatic. We demonstrate that deletions of a small noncoding differentially methylated region at 16q24.1, including lncRNA genes, cause a lethal lung developmental disorder, alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV), with parent-of-origin effects. We identify overlapping deletions 250 kb upstream of FOXF1 in nine patients with ACD/MPV that arose de novo specifically on the maternally inherited chromosome and delete lung-specific lncRNA genes. These deletions define a distant cis-regulatory region that harbors, besides lncRNA genes, also a differentially methylated CpG island, binds GLI2 depending on the methylation status of this CpG island, and physically interacts with and up-regulates the FOXF1 promoter. We suggest that lung-transcribed 16q24.1 lncRNAs may contribute to long-range regulation of FOXF1 by GLI2 and other transcription factors. Perturbation of lncRNA-mediated chromatin interactions may, in general, be responsible for position effect phenomena and potentially cause many disorders of human development.

Uncovering transcription factor modules using one- and three-dimensional analyses

Transcriptional regulation is a critical mediator of many normal cellular processes, as well as disease progression. Transcription factors (TFs) often co-localize at cis-regulatory elements on the DNA, form protein complexes, and collaboratively regulate gene expression. Machine learning and Bayesian approaches have been used to identify TF modules in a one-dimensional context. However, recent studies using high throughput technologies have shown that TF interactions should also be considered in three-dimensional nuclear space. Here, we describe methods for identifying TF modules and discuss how moving from a one-dimensional to a three-dimensional paradigm, along with integrated experimental and computational approaches, can lead to a better understanding of TF association networks.

A novel sarcoidosis risk locus for Europeans on chromosome 11q13.1

RATIONALE

Sarcoidosis is a complex inflammatory disease with a heterogeneous clinical picture. Among others, an acute and chronic clinical course can be distinguished, for which specific genetic risk factors are known.

OBJECTIVES

To identify additional risk loci for sarcoidosis and its acute and chronic subforms, we analyzed imputed data from a genome-wide association scan for these phenotypes.

METHODS

After quality control, the genome-wide association scan comprised nearly 1.3 million imputed single-nucleotide polymorphisms based on an Affymetrix 6.0 Gene Chip dataset of 564 German sarcoidosis cases, including 176 acute and 354 chronic cases and 1,575 control subjects.

MEASUREMENTS AND MAIN RESULTS

We identified chromosome 11q13.1 (rs479777) as a novel locus influencing susceptibility to sarcoidosis with genome-wide significance. The marker was significantly associated in three distinct German case-control populations and in an additional German family sample with odds ratios ranging from 0.67 to 0.77. This finding was further replicated in two independent European case-control populations from the Czech Republic (odds ratio, 0.75) and from Sweden (odds ratio, 0.79). In a meta-analysis of the included European case-control samples the marker yielded a P value of 2.68 × 10(-18). The locus was previously reported to be associated with Crohn disease, psoriasis, alopecia areata, and leprosy. For sarcoidosis, fine-mapping and expression analysis suggest KCNK4, PRDX5, PCLB3, and most promising CCDC88B as candidates for the underlying risk gene in the associated region.

CONCLUSIONS

This study provides striking evidence for association of chromosome 11q13.1 with sarcoidosis in Europeans, and thus identified a further genetic risk locus shared by sarcoidosis, Crohn disease and psoriasis.

A guide for functional analysis of BRCA1 variants of uncertain significance

Germline mutations in the tumor suppressor gene BRCA1 confer an estimated lifetime risk of 56-80% for breast cancer and 15-60% for ovarian cancer. Since the mid 1990s when BRCA1 was identified, genetic testing has revealed over 1,500 unique germline variants. However, for a significant number of these variants, the effect on protein function is unknown making it difficult to infer the consequences on risks of breast and ovarian cancers. Thus, many individuals undergoing genetic testing for BRCA1 mutations receive test results reporting a variant of uncertain clinical significance (VUS), leading to issues in risk assessment, counseling, and preventive care. Here, we describe functional assays for BRCA1 to directly or indirectly assess the impact of a variant on protein conformation or function and how these results can be used to complement genetic data to classify a VUS as to its clinical significance. Importantly, these methods may provide a framework for genome-wide pathogenicity assignment.

Characterization of enhancer function from genome-wide analyses

There has been a recent surge in the use of genome-wide methodologies to identify and annotate the transcriptional regulatory elements in the human genome. Here we review some of these methodologies and the conceptual insights about transcription regulation that have been gained from the use of genome-wide studies. It has become clear that the binding of transcription factors is itself a highly regulated process, and binding does not always appear to have functional consequences. Numerous properties have now been associated with regulatory elements that may be useful in their identification. Several aspects of enhancer function have been shown to be more widespread than was previously appreciated, including the highly combinatorial nature of transcription factor binding, the postinitiation regulation of many target genes, and the binding of enhancers at early stages to maintain their competence during development. Going forward, the integration of multiple genome-wide data sets should become a standard approach to elucidate higher-order regulatory interactions.