Ensembl regulation resources

Bag1 protein loss sensitizes mouse embryonic fibroblasts to glutathione depletion

Bcl2-associated athanogene-1 protein (Bag1) acts as a co-chaperone of heat shock protein 70 and heat shock cognate 70 and regulates multiple cellular processes, including cell proliferation, apoptosis, environmental stress response, and drug resistance. Since Bag1 knockout mice exhibited fetal lethality, the in vivo function of Bag1 remains unclear. In this study, we established a mouse line expressing Bag1 gene missing exon 5, which corresponds to an encoding region for the interface of heat shock protein 70/heat shock cognate 70. Despite mice carrying homoalleles of the Bag1 mutant (Bag1Δex5) expressing undetectable levels of Bag1, Bag1Δex5 homozygous mice developed without abnormalities. Bag1Δex5 protein was found to be highly unstable in cells and in vitro. We found that the growth of mouse embryonic fibroblasts derived from Bag1Δex5-homo mice was attenuated by doxorubicin and a glutathione (GSH) synthesis inhibitor, buthionine sulfoximine. In response to buthionine sulfoximine, Bag1Δex5-mouse embryonic fibroblasts exhibited a higher dropping rate of GSH relative to the oxidized glutathione level. In addition, Bag1 might mitigate cellular hydrogen peroxide levels. Taken together, our results demonstrate that the loss of Bag1 did not affect mouse development and that Bag1 is involved in intracellular GSH homeostasis, namely redox homeostasis.

Genetic Factors Contribute to the Phenotypic Variability in GJB2-Related Hearing Impairment

Recessive variants in GJB2 are the most important genetic cause of sensorineural hearing impairment (SNHI) worldwide. Phenotypes vary significantly in GJB2-related SNHI, even in patients with identical variants. For instance, patients homozygous for the GJB2 p.V37I variant, which is highly prevalent in the Asian populations, usually present with mild-to-moderate SNHI; yet severe-to-profound SNHI is occasionally observed in approximately 10% of p.V37I homozygotes. To investigate the genomic underpinnings of the phenotypic variability, we performed next-generation sequencing of GJB2 and other deafness genes in 63 p.V37I homozygotes with extreme phenotypic severities. Additional pathogenic variants of other deafness genes were identified in five of the 35 patients with severe-to-profound SNHI. Furthermore, case-control association analyses were conducted for 30 unrelated p.V37I homozygotes with severe-to-profound SNHI against 28 p.V37I homozygotes with mild-to-moderate SNHI, and 120 population controls from the Taiwan Biobank. The severe-to-profound group exhibited a higher frequency of the crystallin lambda 1 (CRYL1) variant (rs14236), located upstream of GJB2, than the mild-to-moderate and Taiwan Biobank groups. Our results demonstrated that pathogenic variants in other deafness genes and a possible modifier, the CRYL1 rs14236 variant, may contribute to phenotypic variability in GJB2-realted SNHI, highlighting the importance of comprehensive genomic surveys to delineate the genotype-phenotype correlations.

Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits

BACKGROUND

Insights into the genetic basis of complex traits and disease in both human and livestock species have been achieved over the past decade through detection of genetic variants in genome-wide association studies (GWAS). A majority of such variants were found located in noncoding genomic regions, and though the involvement of numerous regulatory elements (REs) has been predicted across multiple tissues in domesticated animals, their evolutionary conservation and effects on complex traits have not been fully elucidated, particularly in ruminants. Here, we systematically analyzed 137 epigenomic and transcriptomic datasets of six mammals, including cattle, sheep, goats, pigs, mice, and humans, and then integrated them with large-scale GWAS of complex traits.

RESULTS

Using 40 ChIP-seq datasets of H3K4me3 and H3K27ac, we detected 68,479, 58,562, 63,273, 97,244, 111,881, and 87,049 REs in the liver of cattle, sheep, goats, pigs, humans and mice, respectively. We then systematically characterized the dynamic functional landscapes of these REs by integrating multi-omics datasets, including gene expression, chromatin accessibility, and DNA methylation. We identified a core set (n = 6359) of ruminant-specific REs that are involved in liver development, metabolism, and immune processes. Genes with more complex cis-REs exhibited higher gene expression levels and stronger conservation across species. Furthermore, we integrated expression quantitative trait loci (eQTLs) and GWAS from 44 and 52 complex traits/diseases in cattle and humans, respectively. These results demonstrated that REs with different degrees of evolutionary conservation across species exhibited distinct enrichments for GWAS signals of complex traits.

CONCLUSIONS

We systematically annotated genome-wide functional REs in liver across six mammals and demonstrated the evolution of REs and their associations with transcriptional output and conservation. Detecting lineage-specific REs allows us to decipher the evolutionary and genetic basis of complex phenotypes in livestock and humans, which may benefit the discovery of potential biomedical models for functional variants and genes of specific human diseases.

RefSeq Functional Elements as experimentally assayed nongenic reference standards and functional interactions in human and mouse

Eukaryotic genomes contain many nongenic elements that function in gene regulation, chromosome organization, recombination, repair, or replication, and mutation of those elements can affect genome function and cause disease. Although numerous epigenomic studies provide high coverage of gene regulatory regions, those data are not usually exposed in traditional genome annotation and can be difficult to access and interpret without field-specific expertise. The National Center for Biotechnology Information (NCBI) therefore provides RefSeq Functional Elements (RefSeqFEs), which represent experimentally validated human and mouse nongenic elements derived from the literature. The curated data set is comprised of richly annotated sequence records, descriptive records in the NCBI Gene database, reference genome feature annotation, and activity-based interactions between nongenic regions, target genes, and each other. The data set provides succinct functional details and transparent experimental evidence, leverages data from multiple experimental sources, is readily accessible and adaptable, and uses a flexible data model. The data have multiple uses for basic functional discovery, bioinformatics studies, genetic variant interpretation; as known positive controls for epigenomic data evaluation; and as reference standards for functional interactions. Comparisons to other gene regulatory data sets show that the RefSeqFE data set includes a wider range of feature types representing more areas of biology, but it is comparatively smaller and subject to data selection biases. RefSeqFEs thus provide an alternative and complementary resource for experimentally assayed functional elements, with future data set growth expected.

Sequence ontology terminology for gene regulation

The Sequence Ontology (SO) is a structured, controlled vocabulary that provides terms and definitions for genomic annotation. The Gene Regulation Ensemble Effort for the Knowledge Commons (GREEKC) initiative has gathered input from many groups of researchers, including the SO, the Gene Ontology (GO), and gene regulation experts, with the goal of curating information about how gene expression is regulated at the molecular level. Here we discuss recent updates to the SO reflecting current knowledge. We have developed more accurate human-readable terms (also known as classes), including new definitions, and relationships related to the expression of genes. New findings continue to give us insight into the biology of gene regulation, including the order of events, and participants in those events. These updates to the SO support logical reasoning with the current understanding of gene expression regulation at the molecular level.

A zebrafish model of granulin deficiency reveals essential roles in myeloid cell differentiation

Granulin is a pleiotropic protein involved in inflammation, wound healing, neurodegenerative disease, and tumorigenesis. These roles in human health have prompted research efforts to use granulin to treat rheumatoid arthritis and frontotemporal dementia and to enhance wound healing. But how granulin contributes to each of these diverse biological functions remains largely unknown. Here, we have uncovered a new role for granulin during myeloid cell differentiation. We have taken advantage of the tissue-specific segregation of the zebrafish granulin paralogues to assess the functional role of granulin in hematopoiesis without perturbing other tissues. By using our zebrafish model of granulin deficiency, we revealed that during normal and emergency myelopoiesis, myeloid progenitors are unable to terminally differentiate into neutrophils and macrophages in the absence of granulin a (grna), failing to express the myeloid-specific genes cebpa, rgs2, lyz, mpx, mpeg1, mfap4, and apoeb. Functionally, macrophages fail to recruit to the wound, resulting in abnormal healing. Our CUT&RUN experiments identify Pu.1, which together with Irf8, positively regulates grna expression. In vivo imaging and RNA sequencing experiments show that grna inhibits the expression of gata1, leading to the repression of the erythroid program. Importantly, we demonstrated functional conservation between the mammalian granulin and the zebrafish ortholog grna. Our findings uncover a previously unrecognized role for granulin during myeloid cell differentiation, which opens a new field of study that can potentially have an impact on different aspects of human health and expand the therapeutic options for treating myeloid disorders such as neutropenia or myeloid leukemia.

LINE retrotransposons characterize mammalian tissue-specific and evolutionarily dynamic regulatory regions

BACKGROUND

To investigate the mechanisms driving regulatory evolution across tissues, we experimentally mapped promoters, enhancers, and gene expression in the liver, brain, muscle, and testis from ten diverse mammals.

RESULTS

The regulatory landscape around genes included both tissue-shared and tissue-specific regulatory regions, where tissue-specific promoters and enhancers evolved most rapidly. Genomic regions switching between promoters and enhancers were more common across species, and less common across tissues within a single species. Long Interspersed Nuclear Elements (LINEs) played recurrent evolutionary roles: LINE L1s were associated with tissue-specific regulatory regions, whereas more ancient LINE L2s were associated with tissue-shared regulatory regions and with those switching between promoter and enhancer signatures across species.

CONCLUSIONS

Our analyses of the tissue-specificity and evolutionary stability among promoters and enhancers reveal how specific LINE families have helped shape the dynamic mammalian regulome.

Machine learning reveals bilateral distribution of somatic L1 insertions in human neurons and glia

Retrotransposons can cause somatic genome variation in the human nervous system, which is hypothesized to have relevance to brain development and neuropsychiatric disease. However, the detection of individual somatic mobile element insertions presents a difficult signal-to-noise problem. Using a machine-learning method (RetroSom) and deep whole-genome sequencing, we analyzed L1 and Alu retrotransposition in sorted neurons and glia from human brains. We characterized two brain-specific L1 insertions in neurons and glia from a donor with schizophrenia. There was anatomical distribution of the L1 insertions in neurons and glia across both hemispheres, indicating retrotransposition occurred during early embryogenesis. Both insertions were within the introns of genes (CNNM2 and FRMD4A) inside genomic loci associated with neuropsychiatric disorders. Proof-of-principle experiments revealed these L1 insertions significantly reduced gene expression. These results demonstrate that RetroSom has broad applications for studies of brain development and may provide insight into the possible pathological effects of somatic retrotransposition.

Ensembl 2021

The Ensembl project (https://www.ensembl.org) annotates genomes and disseminates genomic data for vertebrate species. We create detailed and comprehensive annotation of gene structures, regulatory elements and variants, and enable comparative genomics by inferring the evolutionary history of genes and genomes. Our integrated genomic data are made available in a variety of ways, including genome browsers, search interfaces, specialist tools such as the Ensembl Variant Effect Predictor, download files and programmatic interfaces. Here, we present recent Ensembl developments including two new website portals. Ensembl Rapid Release (http://rapid.ensembl.org) is designed to provide core tools and services for genomes as soon as possible and has been deployed to support large biodiversity sequencing projects. Our SARS-CoV-2 genome browser (https://covid-19.ensembl.org) integrates our own annotation with publicly available genomic data from numerous sources to facilitate the use of genomics in the international scientific response to the COVID-19 pandemic. We also report on other updates to our annotation resources, tools and services. All Ensembl data and software are freely available without restriction.

A comprehensive microsatellite landscape of human Y-DNA at kilobase resolution

Background

Though interest in human simple sequence repeats (SSRs) is increasing, little is known about the exact distributional features of numerous SSRs in human Y-DNA at chromosomal level. Herein, totally 540 maps were established, which could clearly display SSR landscape in every bin of 1 k base pairs (Kbp) along the sequenced part of human reference Y-DNA (NC_000024.10), by our developed differential method for improving the existing method to reveal SSR distributional characteristics in large genomic sequences.

Results

The maps show that SSRs accumulate significantly with forming density peaks in at least 2040 bins of 1 Kbp, which involve different coding, noncoding and intergenic regions of the Y-DNA, and 10 especially high density peaks were reported to associate with biological significances, suggesting that the other hundreds of especially high density peaks might also be biologically significant and worth further analyzing. In contrast, the maps also show that SSRs are extremely sparse in at least 207 bins of 1 Kbp, including many noncoding and intergenic regions of the Y-DNA, which is inconsistent with the widely accepted view that SSRs are mostly rich in these regions, and these sparse distributions are possibly due to powerfully regional selection. Additionally, many regions harbor SSR clusters with same or similar motif in the Y-DNA.

Conclusions

These 540 maps may provide the important information of clearly position-related SSR distributional features along the human reference Y-DNA for better understanding the genome structures of the Y-DNA. This study may contribute to further exploring the biological significance and distribution law of the huge numbers of SSRs in human Y-DNA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07389-5.

Interpretation of psychiatric genome-wide association studies with multispecies heterogeneous functional genomic data integration

Genome-wide association studies and other discovery genetics methods provide a means to identify previously unknown biological mechanisms underlying behavioral disorders that may point to new therapeutic avenues, augment diagnostic tools, and yield a deeper understanding of the biology of psychiatric conditions. Recent advances in psychiatric genetics have been made possible through large-scale collaborative efforts. These studies have begun to unearth many novel genetic variants associated with psychiatric disorders and behavioral traits in human populations. Significant challenges remain in characterizing the resulting disease-associated genetic variants and prioritizing functional follow-up to make them useful for mechanistic understanding and development of therapeutics. Model organism research has generated extensive genomic data that can provide insight into the neurobiological mechanisms of variant action, but a cohesive effort must be made to establish which aspects of the biological modulation of behavioral traits are evolutionarily conserved across species. Scalable computing, new data integration strategies, and advanced analysis methods outlined in this review provide a framework to efficiently harness model organism data in support of clinically relevant psychiatric phenotypes.

VINYL: Variant prIoritizatioN by survivaL analysis

MOTIVATION

Clinical applications of genome re-sequencing technologies typically generate large amounts of data that need to be carefully annotated and interpreted to identify genetic variants potentially associated with pathological conditions. In this context, accurate and reproducible methods for the functional annotation and prioritization of genetic variants are of fundamental importance.

RESULTS

In this paper, we present VINYL, a flexible and fully automated system for the functional annotation and prioritization of genetic variants. Extensive analyses of both real and simulated datasets suggest that VINYL can identify clinically relevant genetic variants in a more accurate manner compared to equivalent state of the art methods, allowing a more rapid and effective prioritization of genetic variants in different experimental settings. As such we believe that VINYL can establish itself as a valuable tool to assist healthcare operators and researchers in clinical genomics investigations.

AVAILABILITY

VINYL is available at http://beaconlab.it/VINYL and https://github.com/matteo14c/VINYL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ensembl 2019

The Ensembl project (https://www.ensembl.org) makes key genomic data sets available to the entire scientific community without restrictions. Ensembl seeks to be a fundamental resource driving scientific progress by creating, maintaining and updating reference genome annotation and comparative genomics resources. This year we describe our new and expanded gene, variant and comparative annotation capabilities, which led to a 50% increase in the number of vertebrate genomes we support. We have also doubled the number of available human variants and added regulatory regions for many mouse cell types and developmental stages. Our data sets and tools are available via the Ensembl website as well as a through a RESTful webservice, Perl application programming interface and as data files for download.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14+ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14+ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14+ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

Ensembl 2018

The Ensembl project has been aggregating, processing, integrating and redistributing genomic datasets since the initial releases of the draft human genome, with the aim of accelerating genomics research through rapid open distribution of public data. Large amounts of raw data are thus transformed into knowledge, which is made available via a multitude of channels, in particular our browser (http://www.ensembl.org). Over time, we have expanded in multiple directions. First, our resources describe multiple fields of genomics, in particular gene annotation, comparative genomics, genetics and epigenomics. Second, we cover a growing number of genome assemblies; Ensembl Release 90 contains exactly 100. Third, our databases feed simultaneously into an array of services designed around different use cases, ranging from quick browsing to genome-wide bioinformatic analysis. We present here the latest developments of the Ensembl project, with a focus on managing an increasing number of assemblies, supporting efforts in genome interpretation and improving our browser.

Significantly altered peripheral blood cell DNA methylation profile as a result of immediate effect of metformin use in healthy individuals

BACKGROUND

Metformin is a widely prescribed antihyperglycemic agent that has been also associated with multiple therapeutic effects in various diseases, including several types of malignancies. There is growing evidence regarding the contribution of the epigenetic mechanisms in reaching metformin's therapeutic goals; however, the effect of metformin on human cells in vivo is not comprehensively studied. The aim of our study was to examine metformin-induced alterations of DNA methylation profiles in white blood cells of healthy volunteers, employing a longitudinal study design.

RESULTS

Twelve healthy metformin-naïve individuals where enrolled in the study. Genome-wide DNA methylation pattern was estimated at baseline, 10 h and 7 days after the start of metformin administration. The whole-genome DNA methylation analysis in total revealed 125 differentially methylated CpGs, of which 11 CpGs and their associated genes with the most consistent changes in the DNA methylation profile were selected: POFUT2, CAMKK1, EML3, KIAA1614, UPF1, MUC4, LOC727982, SIX3, ADAM8, SNORD12B, VPS8, and several differentially methylated regions as novel potential epigenetic targets of metformin. The main functions of the majority of top-ranked differentially methylated loci and their representative cell signaling pathways were linked to the well-known metformin therapy targets: regulatory processes of energy homeostasis, inflammatory responses, tumorigenesis, and neurodegenerative diseases.

CONCLUSIONS

Here we demonstrate for the first time the immediate effect of short-term metformin administration at therapeutic doses on epigenetic regulation in human white blood cells. These findings suggest the DNA methylation process as one of the mechanisms involved in the action of metformin, thereby revealing novel targets and directions of the molecular mechanisms underlying the various beneficial effects of metformin.

TRIAL REGISTRATION

EU Clinical Trials Register, 2016-001092-74. Registered 23 March 2017, https://www.clinicaltrialsregister.eu/ctr-search/trial/2016-001092-74/LV .

Epigenetic influences on aging: a longitudinal genome-wide methylation study in old Swedish twins

Age-related changes in DNA methylation were observed in cross-sectional studies, but longitudinal evidence is still limited. Here, we aimed to characterize longitudinal age-related methylation patterns using 1011 blood samples collected from 385 Swedish twins (age at entry: mean 69 and standard deviation 9.7, 73 monozygotic and 96 dizygotic pairs) up to five times (mean 2.6) over 20 years (mean 8.7). We identified 1316 age-associated methylation sites (P<1.3×10⁻⁷) using a longitudinal epigenome-wide association study design. We measured how estimated cellular compositions changed with age and how much they confounded the age effect. We validated the results in two independent longitudinal cohorts, where 118 CpGs were replicated in Prospective Investigation of the Vasculature in Uppsala Seniors (PIVUS, 390 samples) (P<3.9×10⁻⁵), 594 in Lothian Birth Cohort (LBC, 3018 samples) (P<5.1×10⁻⁵) and 63 in both. Functional annotation of age-associated CpGs showed enrichment in CCCTC-binding factor (CTCF) and other transcription factor binding sites. We further investigated genetic influences on methylation and found no interaction between age and genetic effects in the 1316 age-associated CpGs. Moreover, in the same CpGs, methylation differences within twin pairs increased with 6.4% over 10 years, where monozygotic twins had smaller intra-pair differences than dizygotic twins. In conclusion, we show that age-related methylation changes persist in a longitudinal perspective, and are fairly stable across cohorts. The changes are under genetic influence, although this effect is independent of age. Moreover, methylation variability increase over time, especially in age-associated CpGs, indicating the increase of environmental contributions on DNA methylation with age.

The chromosome-level genome assemblies of two rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Background

Calamus simplicifolius and Daemonorops jenkinsiana are two representative rattans, the most significant material sources for the rattan industry. However, the lack of reference genome sequences is a major obstacle for basic and applied biology on rattan.

Findings

We produced two chromosome-level genome assemblies of C. simplicifolius and D. jenkinsiana using Illumina, Pacific Biosciences, and Hi-C sequencing data. A total of ∼730 Gb and ∼682 Gb of raw data covered the predicted genome lengths (∼1.98 Gb of C. simplicifolius and ∼1.61 Gb of D. jenkinsiana) to ∼372 × and ∼426 × read depths, respectively. The two de novo genome assemblies, ∼1.94 Gb and ∼1.58 Gb, were generated with scaffold N50s of ∼160 Mb and ∼119 Mb in C. simplicifolius and D. jenkinsiana, respectively. The C. simplicifolius and D. jenkinsiana genomes were predicted to harbor  51,235 and  53,342 intact protein-coding gene models, respectively. Benchmarking Universal Single-Copy Orthologs evaluation demonstrated that genome completeness reached 96.4% and 91.3% in the C. simplicifolius and D. jenkinsiana genomes, respectively. Genome evolution showed that four Arecaceae plants clustered together, and the divergence time between the two rattans was ∼19.3 million years ago. Additionally, we identified 193 and 172 genes involved in the lignin biosynthesis pathway in the C. simplicifolius and D. jenkinsiana genomes, respectively.

Conclusions

We present the first de novo assemblies of two rattan genomes (C. simplicifolius and D. jenkinsiana). These data will not only provide a fundamental resource for functional genomics, particularly in promoting germplasm utilization for breeding, but also serve as reference genomes for comparative studies between and among different species.

Local Delivery of miR-21 Stabilizes Fibrous Caps in Vulnerable Atherosclerotic Lesions

miRNAs are potential regulators of carotid artery stenosis and concordant vulnerable atherosclerotic plaques. Hence, we analyzed miRNA expression in laser captured micro-dissected fibrous caps of either ruptured or stable plaques (n = 10 each), discovering that miR-21 was significantly downregulated in unstable lesions. To functionally evaluate miR-21 in plaque vulnerability, miR-21 and miR-21/apolipoprotein-E double-deficient mice (Apoe-/-miR-21-/-) were assessed. miR-21-/- mice lacked sufficient smooth muscle cell proliferation in response to carotid ligation injury. When exposing Apoe-/-miR-21-/- mice to an inducible plaque rupture model, they presented with more atherothrombotic events (93%) compared with miR-21+/+Apoe-/- mice (57%). We discovered that smooth muscle cell fate in experimentally induced advanced lesions is steered via a REST-miR-21-REST feedback signaling pathway. Furthermore, Apoe-/-miR-21-/- mice presented with more pronounced atherosclerotic lesions, greater foam cell formation, and substantially higher levels of arterial macrophage infiltration. Local delivery of a miR-21 mimic using ultrasound-targeted microbubbles into carotid plaques rescued the vulnerable plaque rupture phenotype. In the present study, we identify miR-21 as a key modulator of pathologic processes in advanced atherosclerosis. Targeted, lesion site-specific overexpression of miR-21 can stabilize vulnerable plaques.

Technological Developments in lncRNA Biology

It is estimated that more than 90% of the mammalian genome is transcribed as non-coding RNAs. Recent evidences have established that these non-coding transcripts are not junk or just transcriptional noise, but they do serve important biological purpose. One of the rapidly expanding fields of this class of transcripts is the regulatory lncRNAs, which had been a major challenge in terms of their molecular functions and mechanisms of action. The emergence of high-throughput technologies and the development in various conventional approaches have led to the expansion of the lncRNA world. The combination of multidisciplinary approaches has proven to be essential to unravel the complexity of their regulatory networks and helped establish the importance of their existence. Here, we review the current methodologies available for discovering and investigating functions of long non-coding RNAs (lncRNAs) and focus on the powerful technological advancement available to specifically address their functional importance.

Ensembl variation resources

The major goal of sequencing humans and many other species is to understand the link between genomic variation, phenotype and disease. There are numerous valuable and well-established variation resources, but collating and making sense of non-homogeneous, often large-scale data sets from disparate sources remains a challenge. Without a systematic catalogue of these data and appropriate query and annotation tools, understanding the genome sequence of an individual and assessing their disease risk is impossible. In Ensembl, we substantially solve this problem: we develop methods to facilitate data integration and broad access; aggregate information in a consistent manner and make it available a variety of standard formats, both visually and programmatically; build analysis pipelines to compare variants to comprehensive genomic annotation sets; and make all tools and data publicly available.

The Ensembl Genome Browser: Strategies for Accessing Eukaryotic Genome Data

The Ensembl Genome Browser provides a wealth of freely available genomic data that can be accessed for many purposes by genetics, genomics, and molecular biology researchers. Herein we present two protocols for exploring different aspects of these data: a phenotype and its associated variants and genes, and a promoter and the epigenetic marks and protein-binding activity associated with it. These workflows illustrate a subset of the data types available through the Ensembl Browser, and can be considered a springboard for further exploration.

GeneHancer: genome-wide integration of enhancers and target genes in GeneCards

A major challenge in understanding gene regulation is the unequivocal identification of enhancer elements and uncovering their connections to genes. We present GeneHancer, a novel database of human enhancers and their inferred target genes, in the framework of GeneCards. First, we integrated a total of 434 000 reported enhancers from four different genome-wide databases: the Encyclopedia of DNA Elements (ENCODE), the Ensembl regulatory build, the functional annotation of the mammalian genome (FANTOM) project and the VISTA Enhancer Browser. Employing an integration algorithm that aims to remove redundancy, GeneHancer portrays 285 000 integrated candidate enhancers (covering 12.4% of the genome), 94 000 of which are derived from more than one source, and each assigned an annotation-derived confidence score. GeneHancer subsequently links enhancers to genes, using: tissue co-expression correlation between genes and enhancer RNAs, as well as enhancer-targeted transcription factor genes; expression quantitative trait loci for variants within enhancers; and capture Hi-C, a promoter-specific genome conformation assay. The individual scores based on each of these four methods, along with gene–enhancer genomic distances, form the basis for GeneHancer’s combinatorial likelihood-based scores for enhancer–gene pairing. Finally, we define ‘elite’ enhancer–gene relations reflecting both a high-likelihood enhancer definition and a strong enhancer–gene association.

GeneHancer predictions are fully integrated in the widely used GeneCards Suite, whereby candidate enhancers and their annotations are displayed on every relevant GeneCard. This assists in the mapping of non-coding variants to enhancers, and via the linked genes, forms a basis for variant–phenotype interpretation of whole-genome sequences in health and disease.

Database URL:http://www.genecards.org/

Fibrinogen and clot-related phenotypes determined by fibrinogen polymorphisms: Independent and IL-6-interactive associations

Interleukin-6 (IL-6) induces the expression of fibrinogen, and polymorphic variation within the fibrinogen genes is believed to alter the magnitude of this expression. The identification of the functional relevance of individual fibrinogen single nucleotide polymorphisms (SNPs) has been hindered by the high linkage disequilibrium (LD) reported in the European fibrinogen gene locus. This study investigated two novel and 12 known fibrinogen SNPs of potential functional relevance, in 2010 Tswana individuals known to have low LD. We aimed to identify functional polymorphisms that contribute to clot-related phenotypes and total and γ’ fibrinogen concentrations independently and through their interaction with IL-6, by taking advantage of the high fibrinogen and IL-6 concentrations and the low LD reported in black South Africans. Fibrinogen was significantly associated with IL-6, thereby mediating associations of IL-6 with clot formation and structure, although attenuating the association of IL-6 with clot lysis time. None of the common European fibrinogen haplotypes was present in this study population. Putative functional fibrinogen SNPs FGB–rs7439150, rs1800789 (–1420_G/A) and rs1800787 (–148_C/T) were significantly associated with fibrinogen concentration and altered clot properties, with several associations significantly influenced by IL-6 concentrations. The impact of harbouring several minor fibrinogen SNP alleles on the association of IL-6 and fibrinogen concentration was cumulative, with possession of each additional minor allele showing a stronger relationship of IL-6 with fibrinogen. This was also reflected in differences in clot properties, suggesting potential clinical relevance. Therefore, when investigating the effect of fibrinogen genetics on fibrinogen concentrations and CVD outcome, the possible interactions with modulating factors and the fact that SNP effects seem to be additive should be taken into account.

Interferon-related genetic markers of necroinflammatory activity in chronic hepatitis C

INTRODUCTION

Chronic hepatitis C (CHC) is a major cause of liver disease worldwide which often leads to progressive liver inflammation, fibrosis, cirrhosis and hepatocellular carcinoma (HCC). CHC displays heterogeneous progression depending on a broad set of factors, some of them intrinsic to each individual such as the patient's genetic profile. This study aims to evaluate the contribution of certain genetic variants of crucial interferon alpha and lambda signaling pathways to the hepatic necroinflammatory activity (NIA) grade of CHC patients.

METHODS

NIA was evaluated in 119 CHC patients by METAVIR scale and classified as low (NIA = 0-2, n = 80) or high grade (NIA = 3, n = 39). In a candidate gene approach, 64 SNPs located in 30 different genes related to interferon pathways (IL-28B, IFNAR1-2, JAK-STAT and OAS1-3, among others) were genotyped using the Illumina GoldenGate® Genotyping Assay. Statistical association was determined by logistic regression and expressed as OR and 95% CI. Those SNPs significantly associated were further adjusted by other covariates.

RESULTS

Seven SNPs located in IL-28B (rs12979860), JAK1 (rs11576173 and rs1497056), TYK2 (rs280519), OAS1 (rs2057778), SOCS1 (rs33932899) and RNASEL (rs3738579) genes were significantly related to severe NIA grade (p<0.05). Regarding to clinical variables, elevated NIA was notably associated with aspartate aminotransferase (AST) serum levels >40 IU/L (p<0.05) but not with other clinical factors. Multivariate logistic regression analysis of these factors reflected that AST (>40 IU/L), TYK2 rs280519 (G allele) and RNASEL rs3738579 (G allele) were factors independently associated with elevated NIA (p<0.05). AST concentration showed a moderate AUC value (AUC = 0.63), similar to TYK2 (rs280519) and RNASEL (rs3738579) SNPs (AUC = 0.61, both) in the ROC_AUC analysis. Interestingly, the model including all significant variables reached a considerable predictive value (AUC = 0.74).

CONCLUSION

The identified genetic variants in interferon signaling pathways may constitute useful prognostic markers of CHC progression. Further validation in larger cohorts of patients is needed.

Hypermethylated CpG sites in the MTR gene promoter in preterm placenta

AIM

Altered maternal one-carbon metabolism influences placental DNA methylation patterns and 'programs' the fetus for noncommunicable diseases in adult life.

EXPERIMENTAL PROCEDURES

Levels of plasma folate, vitamin B₁₂, homocysteine, mRNA and protein levels of MTHFR and MTR enzymes in placenta were compared among women delivering preterm (n = 83) and term (n = 75). MTR promoter CpG methylation was undertaken.

RESULTS

MTHFR and MTR mRNA levels were higher while protein levels were lower, and MTR CpG sites were hypermethylated in the preterm group, as compared with the term group. Methylated CpG sites were negatively associated with maternal plasma vitamin B₁₂ levels.

CONCLUSION

Study suggests a dysregulation of enzyme genes in remethylation arm of the one-carbon metabolism in placenta of women delivering preterm.

Perspectives on Gene Regulatory Network Evolution

Animal development proceeds through the activity of genes and their cis-regulatory modules (CRMs) working together in sets of gene regulatory networks (GRNs). The emergence of species-specific traits and novel structures results from evolutionary changes in GRNs. Recent work in a wide variety of animal models, and particularly in insects, has started to reveal the modes and mechanisms of GRN evolution. I discuss here various aspects of GRN evolution and argue that developmental system drift (DSD), in which conserved phenotype is nevertheless a result of changed genetic interactions, should regularly be viewed from the perspective of GRN evolution. Advances in methods to discover related CRMs in diverse insect species, a critical requirement for detailed GRN characterization, are also described.

GeneHancer: genome-wide integration of enhancers and target genes in GeneCards

A major challenge in understanding gene regulation is the unequivocal identification of enhancer elements and uncovering their connections to genes. We present GeneHancer, a novel database of human enhancers and their inferred target genes, in the framework of GeneCards. First, we integrated a total of 434 000 reported enhancers from four different genome-wide databases: the Encyclopedia of DNA Elements (ENCODE), the Ensembl regulatory build, the functional annotation of the mammalian genome (FANTOM) project and the VISTA Enhancer Browser. Employing an integration algorithm that aims to remove redundancy, GeneHancer portrays 285 000 integrated candidate enhancers (covering 12.4% of the genome), 94 000 of which are derived from more than one source, and each assigned an annotation-derived confidence score. GeneHancer subsequently links enhancers to genes, using: tissue co-expression correlation between genes and enhancer RNAs, as well as enhancer-targeted transcription factor genes; expression quantitative trait loci for variants within enhancers; and capture Hi-C, a promoter-specific genome conformation assay. The individual scores based on each of these four methods, along with gene–enhancer genomic distances, form the basis for GeneHancer’s combinatorial likelihood-based scores for enhancer–gene pairing. Finally, we define ‘elite’ enhancer–gene relations reflecting both a high-likelihood enhancer definition and a strong enhancer–gene association. GeneHancer predictions are fully integrated in the widely used GeneCards Suite, whereby candidate enhancers and their annotations are displayed on every relevant GeneCard. This assists in the mapping of non-coding variants to enhancers, and via the linked genes, forms a basis for variant–phenotype interpretation of whole-genome sequences in health and disease.

DATABASE URL

http://www.genecards.org/.

Ensembl core software resources: storage and programmatic access for DNA sequence and genome annotation

The Ensembl software resources are a stable infrastructure to store, access and manipulate genome assemblies and their functional annotations. The Ensembl 'Core' database and Application Programming Interface (API) was our first major piece of software infrastructure and remains at the centre of all of our genome resources. Since its initial design more than fifteen years ago, the number of publicly available genomic, transcriptomic and proteomic datasets has grown enormously, accelerated by continuous advances in DNA-sequencing technology. Initially intended to provide annotation for the reference human genome, we have extended our framework to support the genomes of all species as well as richer assembly models. Cross-referenced links to other informatics resources facilitate searching our database with a variety of popular identifiers such as UniProt and RefSeq. Our comprehensive and robust framework storing a large diversity of genome annotations in one location serves as a platform for other groups to generate and maintain their own tailored annotation. We welcome reuse and contributions: our databases and APIs are publicly available, all of our source code is released with a permissive Apache v2.0 licence at http://github.com/Ensembl and we have an active developer mailing list ( http://www.ensembl.org/info/about/contact/index.html ).

Database URL

http://www.ensembl.org.

Ensembl 2017

Ensembl (www.ensembl.org) is a database and genome browser for enabling research on vertebrate genomes. We import, analyse, curate and integrate a diverse collection of large-scale reference data to create a more comprehensive view of genome biology than would be possible from any individual dataset. Our extensive data resources include evidence-based gene and regulatory region annotation, genome variation and gene trees. An accompanying suite of tools, infrastructure and programmatic access methods ensure uniform data analysis and distribution for all supported species. Together, these provide a comprehensive solution for large-scale and targeted genomics applications alike. Among many other developments over the past year, we have improved our resources for gene regulation and comparative genomics, and added CRISPR/Cas9 target sites. We released new browser functionality and tools, including improved filtering and prioritization of genome variation, Manhattan plot visualization for linkage disequilibrium and eQTL data, and an ontology search for phenotypes, traits and disease. We have also enhanced data discovery and access with a track hub registry and a selection of new REST end points. All Ensembl data are freely released to the scientific community and our source code is available via the open source Apache 2.0 license.

Lineage-Specific Genome Architecture Links Enhancers and Non-coding Disease Variants to Target Gene Promoters

Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases.

Computational Analysis of G-Quadruplex Forming Sequences across Chromosomes Reveals High Density Patterns Near the Terminal Ends

G-quadruplex structures (G4) are found throughout the human genome and are known to play a regulatory role in a variety of molecular processes. Structurally, they have many configurations and can form from one or more DNA strands. At the gene level, they regulate gene expression and protein synthesis. In this paper, chromosomal-level patterns of distribution are analyzed on the human genome to identify high-level distribution patterns potentially related to global functional processes. Here we show unique high density banding patterns on individual chromosomes that are highly correlated, appearing in a mirror pattern, across forward and reverse DNA strands. The highest density of G4 sequences occurs within four megabases of one end of most chromosomes and contains G4 motifs that bind with zinc finger proteins. These findings suggest that G4 may play a role in global chromosomal processes such as those found in meiosis.

The state of play in higher eukaryote gene annotation

A genome sequence is worthless if it cannot be deciphered; therefore, efforts to describe - or 'annotate' - genes began as soon as DNA sequences became available. Whereas early work focused on individual protein-coding genes, the modern genomic ocean is a complex maelstrom of alternative splicing, non-coding transcription and pseudogenes. Scientists - from clinicians to evolutionary biologists - need to navigate these waters, and this has led to the design of high-throughput, computationally driven annotation projects. The catalogues that are being produced are key resources for genome exploration, especially as they become integrated with expression, epigenomic and variation data sets. Their creation, however, remains challenging.

The Ensembl gene annotation system

The Ensembl gene annotation system has been used to annotate over 70 different vertebrate species across a wide range of genome projects. Furthermore, it generates the automatic alignment-based annotation for the human and mouse GENCODE gene sets. The system is based on the alignment of biological sequences, including cDNAs, proteins and RNA-seq reads, to the target genome in order to construct candidate transcript models. Careful assessment and filtering of these candidate transcripts ultimately leads to the final gene set, which is made available on the Ensembl website. Here, we describe the annotation process in detail.Database URL: http://www.ensembl.org/index.html.