RefSeq: an update on mammalian reference sequences

Identification of X-chromosomal genes that drive sex differences in embryonic stem cells through a hierarchical CRISPR screening approach

BACKGROUND

X-chromosomal genes contribute to sex differences, in particular during early development, when both X chromosomes are active in females. Double X-dosage shifts female pluripotent cells towards the naive stem cell state by increasing pluripotency factor expression, inhibiting the differentiation-promoting MAP kinase (MAPK) signaling pathway, and delaying differentiation.

RESULTS

To identify the genetic basis of these sex differences, we use a two-step CRISPR screening approach to comprehensively identify X-linked genes that cause the female pluripotency phenotype in murine embryonic stem cells. A primary chromosome-wide CRISPR knockout screen and three secondary screens assaying for different aspects of the female pluripotency phenotype allow us to uncover multiple genes that act in concert and to disentangle their relative roles. Among them, we identify Dusp9 and Klhl13 as two central players. While Dusp9 mainly affects MAPK pathway intermediates, Klhl13 promotes pluripotency factor expression and delays differentiation, with both factors jointly repressing MAPK target gene expression.

CONCLUSIONS

Here, we elucidate the mechanisms that drive sex-induced differences in pluripotent cells and our approach serves as a blueprint to discover the genetic basis of the phenotypic consequences of other chromosomal effects.

Unearthing Neanderthal population history using nuclear and mitochondrial DNA from cave sediments

Bones and teeth are important sources of Pleistocene hominin DNA, but are rarely recovered at archaeological sites. Mitochondrial DNA (mtDNA) has been retrieved from cave sediments but provides limited value for studying population relationships. We therefore developed methods for the enrichment and analysis of nuclear DNA from sediments and applied them to cave deposits in western Europe and southern Siberia dated to between 200,000 and 50,000 years ago. We detected a population replacement in northern Spain about 100,000 years ago, which was accompanied by a turnover of mtDNA. We also identified two radiation events in Neanderthal history during the early part of the Late Pleistocene. Our work lays the ground for studying the population history of ancient hominins from trace amounts of nuclear DNA in sediments.

Clinical trait-connected network analysis reveals transcriptional markers of active psoriasis treatment with Liangxue-Jiedu decoction

ETHNOPHARMACOLOGICAL RELEVANCE

Psoriasis is a complex recurrent inflammatory skin disease with different pathological changes in different stages. Psoriasis in its active stage, which is comparable to the blood-heat type in traditional Chinese medicine (TCM), has been treated by Liangxue Jiedu Decoction (LJD) in TCM for decades, with proven efficacy. According to TCM theories, LJD has the function of removing heat and pathogenic factors from the blood.

AIM OF THE STUDY

We aimed to investigate the molecular features associated with the active stage psoriasis and identify genes responding to LJD treatment accompanied by lesion remission.

MATERIALS AND METHODS

Healthy volunteers and psoriasis patients who met specific diagnostic criteria were recruited. Twenty-six transcriptomes were profiled from the peripheral blood mononuclear cells (PBMCs) of 10 psoriasis patients (pre- and post-treatment) and 6 healthy volunteers. RNA sequencing data were analyzed using an integrated approach combining differential gene expression analysis (DGEA) and weighted gene co-expression network analysis (WGCNA), by which gene expression was linked to multiple clinical traits, including psoriasis area and severity index (PASI), as well as the improvement rate of skin lesions (ΔPASI). The actions of LJD were then verified using an in vitro cell assay coupled to flow cytometric analysis and RT-PCR.

RESULTS

We identified four network modules with statistical significance (P < 0.05), two of which connected to the PASI score, while the other two connected to 8-week treatment and ΔPASI, respectively. In psoriasis patients, activated inflammatory pathways and inhibited G-protein signaling genes (GTPase IMAP family member and G protein-coupled receptor) co-occurred, with high expression of CD83 and CD69, and low expression of CD160 and CD180, compared with the health. Accompanying LJD treatment and lesion remission, the expression of CD69 and cell cycle-related genes, including CCNA2, CCNB2, CDK1, and TOP2A, was down-regulated. The inhibitory role of LJD on CD69 expression was confirmed by the decline of activating naïve CD4⁺ T lymphocytes.

CONCLUSION

Our study suggests that active psoriasis is characterized by unbalanced immune status with dendrite cell and lymphocyte-associated inflammatory activation as well as NK cell- and B cell-associated defense response aberrance. LJD played an inhibitory role in T cell activation, a process located downstream pathological cascade of psoriasis.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

A chromosome-level genome of Astyanax mexicanus surface fish for comparing population-specific genetic differences contributing to trait evolution

Identifying the genetic factors that underlie complex traits is central to understanding the mechanistic underpinnings of evolution. Cave-dwelling Astyanax mexicanus populations are well adapted to subterranean life and many populations appear to have evolved troglomorphic traits independently, while the surface-dwelling populations can be used as a proxy for the ancestral form. Here we present a high-resolution, chromosome-level surface fish genome, enabling the first genome-wide comparison between surface fish and cavefish populations. Using this resource, we performed quantitative trait locus (QTL) mapping analyses and found new candidate genes for eye loss such as dusp26. We used CRISPR gene editing in A. mexicanus to confirm the essential role of a gene within an eye size QTL, rx3, in eye formation. We also generated the first genome-wide evaluation of deletion variability across cavefish populations to gain insight into this potential source of cave adaptation. The surface fish genome reference now provides a more complete resource for comparative, functional and genetic studies of drastic trait differences within a species.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

REV-ERBα alters circadian rhythms by modulating mTOR signaling

REV-ERBα is a nuclear receptor that inhibits Bmal1 transcription as part of the circadian clock molecular mechanism. Mammalian target of rapamycin (mTOR) complex 1 (mTORC1) is a master regulator of cell and whole-body energy homeostasis, that serves as an important link between metabolism and circadian clock, in part, by regulating BMAL1 activity. While the connection of REV-ERBα to the circadian clock molecular mechanism is well characterized, the interaction between mTORC1, REV-ERBα and the circadian clock machinery is not very clear. We used leucine and rapamycin to modulate mTORC1 activation and evaluate this effect on circadian rhythms. In the liver, mTORC1 was inhibited by leucine. REV-ERBα overexpression activated the mTORC1 signaling pathway via transcription inhibition of mTORC1 inhibitor, Tsc1, antagonizing the effect of leucine, while its silencing downregulated mTORC1 signaling. Activation of mTORC1 led to increased BMAL1 phosphorylation. Activation as well as inhibition of mTORC1 led to altered circadian rhythms in mouse muscle. Inhibition of liver mTORC1 by leucine or rapamycin led to low-amplitude circadian rhythms. In summary, our study shows that leucine inhibits liver mTORC1 pathway leading to dampened circadian rhythms. REV-ERBα activates the mTORC1 pathway, leading to phosphorylation of the clock protein BMAL1.

MicroExonator enables systematic discovery and quantification of microexons across mouse embryonic development

BACKGROUND

Microexons, exons that are ≤ 30 nucleotides, are a highly conserved and dynamically regulated set of cassette exons. They have key roles in nervous system development and function, as evidenced by recent results demonstrating the impact of microexons on behaviour and cognition. However, microexons are often overlooked due to the difficulty of detecting them using standard RNA-seq aligners.

RESULTS

Here, we present MicroExonator, a novel pipeline for reproducible de novo discovery and quantification of microexons. We process 289 RNA-seq datasets from eighteen mouse tissues corresponding to nine embryonic and postnatal stages, providing the most comprehensive survey of microexons available for mice. We detect 2984 microexons, 332 of which are differentially spliced throughout mouse embryonic brain development, including 29 that are not present in mouse transcript annotation databases. Unsupervised clustering of microexons based on their inclusion patterns segregates brain tissues by developmental time, and further analysis suggests a key function for microexons in axon growth and synapse formation. Finally, we analyse single-cell RNA-seq data from the mouse visual cortex, and for the first time, we report differential inclusion between neuronal subpopulations, suggesting that some microexons could be cell type-specific.

CONCLUSIONS

MicroExonator facilitates the investigation of microexons in transcriptome studies, particularly when analysing large volumes of data. As a proof of principle, we use MicroExonator to analyse a large collection of both mouse bulk and single-cell RNA-seq datasets. The analyses enabled the discovery of previously uncharacterized microexons, and our study provides a comprehensive microexon inclusion catalogue during mouse development.

A comprehensive analysis of the phylogenetic signal in ramp sequences in 211 vertebrates

Ramp sequences increase translational speed and accuracy when rare, slowly-translated codons are found at the beginnings of genes. Here, the results of the first analysis of ramp sequences in a phylogenetic construct are presented. Ramp sequences were compared from 247 vertebrates (114 Mammalian and 133 non-mammalian), where the presence and absence of ramp sequences was analyzed as a binary character in a parsimony and maximum likelihood framework. Additionally, ramp sequences were mapped to the Open Tree of Life synthetic tree to determine the number of parallelisms and reversals that occurred, and those results were compared to random permutations. Parsimony and maximum likelihood analyses of the presence and absence of ramp sequences recovered phylogenies that are highly congruent with established phylogenies. Additionally, 81% of vertebrate mammalian ramps and 81.2% of other vertebrate ramps had less parallelisms and reversals than the mean from 1000 randomly permuted trees. A chi-square analysis of completely orthologous ramp sequences resulted in a p-value < 0.001 as compared to random chance. Ramp sequences recover comparable phylogenies as other phylogenomic methods. Although not all ramp sequences appear to have a phylogenetic signal, more ramp sequences track speciation than expected by random chance. Therefore, ramp sequences may be used in conjunction with other phylogenomic approaches if many orthologs are taken into account. However, phylogenomic methods utilizing few orthologs should be cautious in incorporating ramp sequences because individual ramp sequences may provide conflicting signals.

HGFDB: a collective database of helmeted guinea fowl genomics

As a vigorous and hardy and an almost disease-free game bird, the domestic helmeted guinea fowl (Numida meleagris, hereafter HGF) has attracted considerable attention in a large number of genetic study projects. However, none of the current/recent avian databases are related to this agriculturally and commercially important poultry species. To address this data gap, we developed Helmeted Guinea Fowl Database (HGFDB), which manages and shares HGF genomic and genetic data. By processing the data of genome assembly, sequencing reads and genetic variations, we organized them into eight modules, which correspond to ‘Home’, ‘Genome’, ‘Re-sequence’, ‘Gene’, ‘Variation’, ‘Download’, ‘Tools’ and ‘Help’, HGFDB provides the most comprehensive view of the HGF genome to date and will be relevant for future studies on HGF structural and functional genomics and genetic improvement.

Database URL:http://hgfdb.ynau.edu.cn/

OMA orthology in 2021: website overhaul, conserved isoforms, ancestral gene order and more

OMA is an established resource to elucidate evolutionary relationships among genes from currently 2326 genomes covering all domains of life. OMA provides pairwise and groupwise orthologs, functional annotations, local and global gene order conservation (synteny) information, among many other functions. This update paper describes the reorganisation of the database into gene-, group- and genome-centric pages. Other new and improved features are detailed, such as reporting of the evolutionarily best conserved isoforms of alternatively spliced genes, the inferred local order of ancestral genes, phylogenetic profiling, better cross-references, fast genome mapping, semantic data sharing via RDF, as well as a special coronavirus OMA with 119 viruses from the Nidovirales order, including SARS-CoV-2, the agent of the COVID-19 pandemic. We conclude with improvements to the documentation of the resource through primers, tutorials and short videos. OMA is accessible at https://omabrowser.org.

Software Options for the Analysis of MS-Proteomic Data

Mass spectrometry (MS)-based proteomics is currently the most successful approach to measure and compare peptides and proteins in a large variety of biological samples. Modern mass spectrometers, equipped with high-resolution analyzers, provide large amounts of data output. This is the case of shotgun/bottom-up proteomics, which consists in the enzymatic digestion of protein into peptides that are then measured by MS-instruments through a data dependent acquisition (DDA) mode. Dedicated bioinformatic tools and platforms have been developed to face the increasing size and complexity of raw MS data that need to be processed and interpreted for large-scale protein identification and quantification. This chapter illustrates the most popular bioinformatics solution for the analysis of shotgun MS-proteomics data. A general description will be provided on the data preprocessing options and the different search engines available, including practical suggestions on how to optimize the parameters for peptide search, based on hands-on experience.

A long-term study of AAV gene therapy in dogs with hemophilia A identifies clonal expansions of transduced liver cells

Nine dogs with hemophilia A were treated with adeno-associated viral (AAV) gene therapy and followed for up to 10 years. Administration of AAV8 or AAV9 vectors expressing canine factor VIII (AAV-cFVIII) corrected the FVIII deficiency to 1.9-11.3% of normal FVIII levels. In two of nine dogs, levels of FVIII activity increased gradually starting about 4 years after treatment. None of the dogs showed evidence of tumors or altered liver function. Analysis of integration sites in liver samples from six treated dogs identified 1,741 unique AAV integration events in genomic DNA and expanded cell clones in five dogs, with 44% of the integrations near genes involved in cell growth. All recovered integrated vectors were partially deleted and/or rearranged. Our data suggest that the increase in FVIII protein expression in two dogs may have been due to clonal expansion of cells harboring integrated vectors. These results support the clinical development of liver-directed AAV gene therapy for hemophilia A, while emphasizing the importance of long-term monitoring for potential genotoxicity.

Phenotypic differences of mutation-negative cases in Gitelman syndrome clinically diagnosed in adulthood

Gitelman syndrome (GS), an autosomal recessive kidney disorder, is characterized by hypokalemia, hypomagnesemia, hypocalciuria, and metabolic alkalosis. Generally, diagnosis is made in school-aged children but multiple cases have been diagnosed in adulthood. This study examines the phenotypic differences between genetically confirmed cases and mutation-negative cases in adults. A comprehensive screening of 168 genes, including GS-related genes, was performed for 84 independent individuals who were referred to our institute with a clinical diagnosis of GS. The cases of pseudo-Bartter syndrome (BS)/GS because of diuretic abuse or other causes, which was determined based on patients' medical records, were excluded during registration. Of these 70 eligible cases for analysis, 27 (38.6%) had genetic confirmation of GS, while 37 (52.8%) had no known variants associated with GS and were considered to be unsolved cases. Note that unsolved cases comprised older, mostly female, individuals with decreased kidney function and multiple basic features of GS. The phenotype of unsolved cases is similar to that of pseudo BS/GS cases, although these cases were excluded in advance. However, the genetic and autoimmune profiles of these unsolved cases have not yet been investigated to date. Therefore, these cases may be categorized into new disease groups.

Predicting the Disease Risk of Protein Mutation Sequences With Pre-training Model

Accurately identifying the missense mutations is of great help to alleviate the loss of protein function and structural changes, which might greatly reduce the risk of disease for tumor suppressor genes (e.g., BRCA1 and PTEN). In this paper, we propose a hybrid framework, called BertVS, that predicts the disease risk for the missense mutation of proteins. Our framework is able to learn sequence representations from the protein domain through pre-training BERT models, and also integrates with the hydrophilic properties of amino acids to obtain the sequence representations of biochemical characteristics. The concatenation of two learned representations are then sent to the classifier to predict the missense mutations of protein sequences. Specifically, we use the protein family database (Pfam) as a corpus to train the BERT model to learn the contextual information of protein sequences, and our pre-training BERT model achieves a value of 0.984 on accuracy in the masked language model prediction task. We conduct extensive experiments on BRCA1 and PTEN datasets. With comparison to the baselines, results show that BertVS achieves higher performance of 0.920 on AUROC and 0.915 on AUPR in the functionally critical domain of the BRCA1 gene. Additionally, the extended experiment on the ClinVar dataset can illustrate that gene variants with known clinical significance can also be efficiently classified by our method. Therefore, BertVS can learn the functional information of the protein sequences and effectively predict the disease risk of variants with an uncertain clinical significance.

Annotation of Full-Length Long Noncoding RNAs with Capture Long-Read Sequencing (CLS)

Metazoan genomes produce thousands of long-noncoding RNAs (lncRNAs), of which just a small fraction have been well characterized. Understanding their biological functions requires accurate annotations, or maps of the precise location and structure of genes and transcripts in the genome. Current lncRNA annotations are limited by compromises between quality and size, with many gene models being fragmentary or uncatalogued. To overcome this, the GENCODE consortium has developed RNA capture long-read sequencing (CLS), an approach combining targeted RNA capture with third-generation long-read sequencing. CLS provides accurate annotations at high-throughput rates. It eliminates the need for noisy transcriptome assembly from short reads, and requires minimal manual curation. The full-length transcript models produced are of quality comparable to present-day manually curated annotations. Here we describe a detailed CLS protocol, from probe design through long-read sequencing to creation of final annotations.

Mitochondria Autoimmunity and MNRR1 in Breast Carcinogenesis: A Review

We review here the evidence for participation of mitochondrial autoimmunity in BC inception and progression and propose a new paradigm that may challenge the prevailing thinking in oncogenesis by suggesting that mitochondrial autoimmunity is a major contributor to breast carcinogenesis and probably to the inception and progression of other solid tumors. It has been shown that MNRR1 mediated mitochondrial-nuclear function promotes BC cell growth and migration and the development of metastasis and constitutes a proof of concept supporting the participation of mitochondrial autoimmunity in breast carcinogenesis. The resemblance of the autoantibody profile in BC detected by IFA with that in the rheumatic autoimmune diseases suggested that studies on the autoantibody response to tumor associated antigens and the characterization of the mtDNA- and nDNA-encoded antigens may provide functional data on breast carcinogenesis. We also review the studies supporting the view that a panel of autoreactive nDNA-encoded mitochondrial antigens in addition to MNRR1 may be involved in breast carcinogenesis. These include GAPDH, PKM2, GSTP1, SPATA5, MFF, ncRNA PINK1-AS/DDOST as probably contributing to BC progression and metastases and the evidence suggesting that DDX21 orchestrates a complex signaling network with participation of JUND and ATF3 driving chronic inflammation and breast tumorigenesis. We suggest that the widespread autoreactivity of mtDNA- and nDNA-encoded mitochondrial proteins found in BC sera may be the reflection of autoimmunity triggered by mitochondrial and non-mitochondrial tumor associated antigens involved in multiple tumorigenic pathways. Furthermore, we suggest that mitochondrial proteins may contribute to mitochondrial dysfunction in BC even if mitochondrial respiration is found to be within normal limits. However, although the studies show that mitochondrial autoimmunity is a major factor in breast cancer inception and progression, it is not the only factor since there is a multiplex autoantibody profile targeting centrosome and stem cell antigens as well as anti-idiotypic antibodies, revealing the complex signaling network involved in breast carcinogenesis. In summary, the studies reviewed here open new, unexpected therapeutic avenues for cancer prevention and treatment of patients with cancer derived from an entirely new perspective of breast carcinogenesis.

Frameshift Variant in Novel Adenosine-A1-Receptor Homolog Associated With Bovine Spastic Syndrome/Late-Onset Bovine Spastic Paresis in Holstein Sires

Since their first description almost 100 years ago, bovine spastic paresis (BSP) and bovine spastic syndrome (BSS) are assumed to be inherited neuronal-progressive diseases in cattle. Affected animals are characterized by (frequent) spasms primarily located in the hind limbs, accompanied by severe pain symptoms and reduced vigor, thus initiating premature slaughter or euthanasia. Due to the late onset of BSP and BSS and the massively decreased lifespan of modern cattle, the importance of these diseases is underestimated. In the present study, BSP/BSS-affected German Holstein breeding sires from artificial insemination centers were collected and pedigree analysis, genome-wide association studies, whole genome resequencing, protein-protein interaction network analysis, and protein-homology modeling were performed to elucidate the genetic background. The analysis of 46 affected and 213 control cattle revealed four significantly associated positions on chromosome 15 (BTA15), i.e., AC_000172.1:g.83465449A>G (-log₁₀P = 19.17), AC_000172.1:g.81871849C>T (-log₁₀P = 8.31), AC_000172.1:g.81872621A>T (-log₁₀P = 6.81), and AC_000172.1:g.81872661G>C (-log₁₀P = 6.42). Two additional loci were significantly associated located on BTA8 and BTA19, i.e., AC_000165.1:g.71177788T>C and AC_000176.1:g.30140977T>G, respectively. Whole genome resequencing of five affected individuals and six unaffected relatives (two fathers, two mothers, a half sibling, and a full sibling) belonging to three different not directly related families was performed. After filtering, a homozygous loss of function variant was identified in the affected cattle, causing a frameshift in the so far unknown gene locus LOC100848076 encoding an adenosine-A1-receptor homolog. An allele frequency of the variant of 0.74 was determined in 3,093 samples of the 1000 Bull Genomes Project.

Newly established gastrointestinal cancer cell lines retain the genomic and immunophenotypic landscape of their parental cancers

Human cancer cell lines are frequently used as model systems to study molecular mechanisms and genetic changes in cancer. However, the model is repeatedly criticized for its lack of proximity to original patient tumors. Therefore, understanding to what extent cell lines cultured under artificial conditions reflect the phenotypic and genomic profiles of their corresponding parental tumors is crucial when analyzing their biological properties. To directly compare molecular alterations between patient tumors and derived cell lines, we have established new cancer cell lines from four patients with gastrointestinal tumors. Tumor entities comprised esophageal cancer, colon cancer, rectal cancer and pancreatic cancer. Phenotype and genotype of both patient tumors and derived low-passage cell lines were characterized by immunohistochemistry (22 different antibodies), array-based comparative genomic hybridization and targeted next generation sequencing (48-gene panel). The immunophenotype was highly consistent between patient tumors and derived cell lines; the expression of most markers in cell lines was concordant with the respective parental tumor and characteristic for the respective tumor entities in general. The chromosomal aberration patterns of the parental tumors were largely maintained in the cell lines and the distribution of gains and losses was typical for the respective cancer entity, despite a few distinct differences. Cancer gene mutations (e.g., KRAS, TP53) and microsatellite status were also preserved in the respective cell line derivates. In conclusion, the four examined newly established cell lines exhibited a phenotype and genotype closely recapitulating their parental tumor. Hence, newly established cancer cell lines may be useful models for further pharmacogenomic studies.

Non-invasive epigenomic molecular phenotyping of the human brain via liquid biopsy of cerebrospinal fluid and next generation sequencing

Our goal was to undertake a genome-wide epigenomic liquid biopsy of cerebrospinal fluid (CSF) for the comprehensive analysis of cell-free DNA (cfDNA) methylation signatures in the human central nervous system (CNS). Solution-phase hybridization and massively parallel sequencing of bisulfite converted human DNA was employed to compare methylation signatures of cfDNA obtained from CSF with plasma. Recovery of cfDNA from CSF was relatively low (68-840 pg/mL) compared to plasma (2720-8390 pg/mL) and cfDNA fragments from CSF were approximately 20 bp shorter than their plasma-derived counterparts. Distributions of CpG methylation signatures were significantly altered between CSF and plasma, both globally and at the level of functional elements including exons, introns, CpG islands, and shores. Sliding window analysis was used to identify differentially methylated regions. We found numerous gene/locus-specific differences in CpG methylation between cfDNA from CSF and plasma. These loci were more frequently hypomethylated in CSF compared to plasma. Differentially methylated CpGs in CSF were identified in genes related to branching of neurites and neuronal development. Using the GTEx RNA expression database, we found clear association between tissue-specific gene expression in the CNS and cfDNA methylation patterns in CSF. Ingenuity pathway analysis of differentially methylated regions identified an enrichment of functional pathways related to neurobiology. In conclusion, we present a genome-wide analysis of DNA methylation in human CSF. Our methods and the resulting data demonstrate the potential of epigenomic liquid biopsy of the human CNS for molecular phenotyping of brain-derived DNA methylation signatures.

Splice-Junction-Based Mapping of Alternative Isoforms in the Human Proteome

The protein-level translational status and function of many alternative splicing events remain poorly understood. We use an RNA sequencing (RNA-seq)-guided proteomics method to identify protein alternative splicing isoforms in the human proteome by constructing tissue-specific protein databases that prioritize transcript splice junction pairs with high translational potential. Using the custom databases to reanalyze ~80 million mass spectra in public proteomics datasets, we identify more than 1,500 noncanonical protein isoforms across 12 human tissues, including ~400 sequences undocumented on TrEMBL and RefSeq databases. We apply the method to original quantitative mass spectrometry experiments and observe widespread isoform regulation during human induced pluripotent stem cell cardiomyocyte differentiation. On a proteome scale, alternative isoform regions overlap frequently with disordered sequences and post-translational modification sites, suggesting that alternative splicing may regulate protein function through modulating intrinsically disordered regions. The described approach may help elucidate functional consequences of alternative splicing and expand the scope of proteomics investigations in various systems.

The translation and function of many alternative splicing events await confirmation at the protein level. Lau et al. use an integrated proteotranscriptomics approach to identify non-canonical and undocumented isoforms from 12 organs in the human proteome. Alternative isoforms interfere with functional sequence features and are differentially regulated during iPSC cardiomyocyte differentiation.

Protein Databases Related to Liquid-Liquid Phase Separation

Liquid-liquid phase separation (LLPS) of biomolecules, which underlies the formation of membraneless organelles (MLOs) or biomolecular condensates, has been investigated intensively in recent years. It contributes to the regulation of various physiological processes and related disease development. A rapidly increasing number of studies have recently focused on the biological functions, driving, and regulating mechanisms of LLPS in cells. Based on the mounting data generated in the investigations, six databases (LLPSDB, PhaSePro, PhaSepDB, DrLLPS, RNAgranuleDB, HUMAN CELL MAP) have been developed, which are designed directly based on LLPS studies or the component identification of MLOs. These resources are invaluable for a deeper understanding of the cellular function of biomolecular phase separation, as well as the development of phase-separating protein prediction and design. In this review, we compare the data contents, annotations, and organization of these databases, highlight their unique features, overlaps, and fundamental differences, and discuss their suitable applications.

PCOSKBR2: a database of genes, diseases, pathways, and networks associated with polycystic ovary syndrome

PolyCystic Ovary Syndrome KnowledgeBase (PCOSKBR2) is a manually curated database with information on 533 genes, 145 SNPs, 29 miRNAs, 1,150 pathways, and 1,237 diseases associated with PCOS. This data has been retrieved based on evidence gleaned by critically reviewing literature and related records available for PCOS in databases such as KEGG, DisGeNET, OMIM, GO, Reactome, STRING, and dbSNP. Since PCOS is associated with multiple genes and comorbidities, data mining algorithms for comorbidity prediction and identification of enriched pathways and hub genes are integrated in PCOSKBR2, making it an ideal research platform for PCOS. PCOSKBR2 is freely accessible at http://www.pcoskb.bicnirrh.res.in/ .

Influenza A virus utilizes noncanonical cap-snatching to diversify its mRNA/ncRNA

Influenza A virus (IAV) utilizes cap-snatching to obtain host capped small RNAs for priming viral mRNA synthesis, generating capped hybrid mRNAs for translation. Previous studies have been focusing on canonical cap-snatching, which occurs at the very 5' end of viral mRNAs. Here we discovered noncanonical cap-snatching, which generates capped hybrid mRNAs/noncoding RNAs mapped to the region ∼300 nucleotides (nt) upstream of each mRNA 3' end, and to the 5' region, primarily starting at the second nt, of each virion RNAs (vRNA). Like canonical cap-snatching, noncanonical cap-snatching utilizes a base-pairing between the last nt G of host capped RNAs and a nt C of template RNAs to prime RNA synthesis. However, the nt upstream of this template C is usually A/U rather than just U; prime-realignment occurs less frequently. We also demonstrate that IAV can snatch capped IAV RNAs in addition to host RNAs. Noncanonical cap-snatching likely generates novel mRNAs with start AUG encoded in viral or host RNAs. These findings expand our understanding of cap-snatching mechanisms and suggest that IAV may utilize noncanonical cap-snatching to diversify its mRNAs/ncRNAs.

Genome-Wide Transcriptional Analysis Reveals Alternative Splicing Event Profiles in Hepatocellular Carcinoma and Their Prognostic Significance

Accumulating evidence indicates an unexpected role of aberrant splicing in hepatocellular carcinoma (HCC) that has been seriously neglected in previous studies. There is a need for a detailed analysis of alternative splicing (AS) and its underlying biological and clinical relevance in HCC. In this study, clinical information and corresponding RNA sequencing data of HCC patients were obtained from The Cancer Genome Atlas. Percent spliced in (PSI) values and transcriptional splicing patterns of genes were determined from the original RNA sequencing data using SpliceSeq. Then, based on the PSI values of AS events in different patients, a series of bioinformatics methods was used to identify differentially expressed AS events (DEAS), determine potential regulatory relationships, and investigate the correlation between DEAS and the patients' clinicopathological features. Finally, 25,934 AS events originating from 8,795 genes were screened with high reliability; 263 of these AS events were identified as DEAS. The parent genes of these DEAS formed an intricate network with roles in the regulation of cancer-related pathway and liver metabolism. In HCC, 36 splicing factors were involved in the dysregulation of part DEAS, 100 DEAS events were correlated with overall survival, and 71 DEAS events were correlated with disease-free survival. Stratifying HCC patients according to DEAS resulted in four clusters with different survival patterns. Significant variations in AS occurred during HCC initiation and maintenance; these are likely to be vital both for biological processes and in prognosis. The HCC-related AS events identified here and the splicing networks constructed will be valuable in deciphering the underlying role of AS in HCC.

Exploring the functional impact of alternative splicing on human protein isoforms using available annotation sources

In recent years, the emphasis of scientific inquiry has shifted from whole-genome analyses to an understanding of cellular responses specific to tissue, developmental stage or environmental conditions. One of the central mechanisms underlying the diversity and adaptability of the contextual responses is alternative splicing (AS). It enables a single gene to encode multiple isoforms with distinct biological functions. However, to date, the functions of the vast majority of differentially spliced protein isoforms are not known. Integration of genomic, proteomic, functional, phenotypic and contextual information is essential for supporting isoform-based modeling and analysis. Such integrative proteogenomics approaches promise to provide insights into the functions of the alternatively spliced protein isoforms and provide high-confidence hypotheses to be validated experimentally. This manuscript provides a survey of the public databases supporting isoform-based biology. It also presents an overview of the potential global impact of AS on the human canonical gene functions, molecular interactions and cellular pathways.

A practical guide for DNase-seq data analysis: from data management to common applications

Deoxyribonuclease I (DNase I)-hypersensitive site sequencing (DNase-seq) has been widely used to determine chromatin accessibility and its underlying regulatory lexicon. However, exploring DNase-seq data requires sophisticated downstream bioinformatics analyses. In this study, we first review computational methods for all of the major steps in DNase-seq data analysis, including experimental design, quality control, read alignment, peak calling, annotation of cis-regulatory elements, genomic footprinting and visualization. The challenges associated with each step are highlighted. Next, we provide a practical guideline and a computational pipeline for DNase-seq data analysis by integrating some of these tools. We also discuss the competing techniques and the potential applications of this pipeline for the analysis of analogous experimental data. Finally, we discuss the integration of DNase-seq with other functional genomics techniques.

RNA profiling of human testicular cells identifies syntenic lncRNAs associated with spermatogenesis

STUDY QUESTION

Is the noncoding transcriptional landscape during spermatogenesis conserved between human and rodents?

SUMMARY ANSWER

We identified a core group of 113 long noncoding RNAs (lncRNAs) and 20 novel genes dynamically and syntenically transcribed during spermatogenesis.

WHAT IS KNOWN ALREADY

Spermatogenesis is a complex differentiation process driven by a tightly regulated and highly specific gene expression program. Recently, several studies in various species have established that a large proportion of known lncRNAs are preferentially expressed during meiosis and spermiogenesis in a testis-specific manner.

STUDY DESIGN, SIZE, DURATION

To further investigate lncRNA expression in human spermatogenesis, we carried out a cross-species RNA profiling study using isolated testicular cells.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Human testes were obtained from post-mortem donors (N = 8, 51 years old on average) or from prostate cancer patients with no hormonal treatment (N = 9, 80 years old on average) and only patients with full spermatogenesis were used to prepare enriched populations of spermatocytes, spermatids, Leydig cells, peritubular cells and Sertoli cells. To minimize potential biases linked to inter-patient variations, RNAs from two or three donors were pooled prior to RNA-sequencing (paired-end, strand-specific). Resulting reads were mapped to the human genome, allowing for assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

Our RNA-sequencing analysis of pools of isolated human testicular cells enabled us to reconstruct over 25 000 transcripts. Among them we identified thousands of lncRNAs, as well as many previously unidentified genes (novel unannotated transcripts) that share many properties of lncRNAs. Of note is that although noncoding genes showed much lower synteny than protein-coding ones, a significant fraction of syntenic lncRNAs displayed conserved expression during spermatogenesis.

LARGE SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI Gene Expression Omnibus under accession number GSE74896.

LIMITATIONS, REASONS FOR CAUTION

Isolation procedures may alter the physiological state of testicular cells, especially for somatic cells, leading to substantial changes at the transcriptome level. We therefore cross-validated our findings with three previously published transcriptomic analyses of human spermatogenesis. Despite the use of stringent filtration criteria, i.e. expression cut-off of at least three fragments per kilobase of exon model per million reads mapped, fold-change of at least three and false discovery rate adjusted P-values of less than <1%, the possibility of assembly artifacts and false-positive transcripts cannot be fully ruled out.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of a large number of conserved germline-associated lncRNAs that are potentially important for spermatogenesis and sexual reproduction. In addition to further substantiating the basis of the human testicular physiology, our study provides new candidate genes for male infertility of genetic origin. This is likely to be relevant for identifying interesting diagnostic and prognostic biomarkers and also potential novel therapeutic targets for male contraception.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by l'Institut national de la santé et de la recherche médicale (Inserm); l'Université de Rennes 1; l'Ecole des hautes études en santé publique (EHESP); INERIS-STORM to B.J. [N 10028NN]; Rennes Métropole 'Défis scientifiques émergents' to F.C (2011) and A.D.R (2013). The authors have no competing financial interests.

Diversification of CpG-Island Promoters Revealed by Comparative Analysis Between Human and Rhesus Monkey Genomes

While CpG dinucleotides are significantly reduced compared to other dinucleotides in mammalian genomes, they can congregate and form CpG islands, which localize around the 5' regions of genes, where they function as promoters. CpG-island promoters are generally unmethylated and are often found in housekeeping genes. However, their nucleotide sequences and existence per se are not conserved between humans and mice, which may be due to evolutionary gain and loss of the regulatory regions. In this study, human and rhesus monkey genomes, with moderately conserved sequences, were compared at base resolution. Using transcription start site data, we first validated our methods' ability to identify orthologous promoters and indicated a limitation using the 5' end of curated gene models, such as NCBI RefSeq, as their transcription start sites. We found that, in addition to deamination mutations, insertions and deletions of bases, repeats, and long fragments contributed to the mutations of CpG dinucleotides. We also observed that the G + C contents tended to change in CpG-poor environments, while CpG content was altered in G + C-rich environments. While loss of CpG islands can be caused by gradual decreases in CpG sites, gain of these islands appear to require two distinct nucleotide altering steps. Taken together, our findings provide novel insights into the process of acquisition and diversification of CpG-island promoters in vertebrates.

MeCP2 regulates gene expression through recognition of H3K27me3

MeCP2 plays a multifaceted role in gene expression regulation and chromatin organization. Interaction between MeCP2 and methylated DNA in the regulation of gene expression is well established. However, the widespread distribution of MeCP2 suggests it has additional interactions with chromatin. Here we demonstrate, by both biochemical and genomic analyses, that MeCP2 directly interacts with nucleosomes and its genomic distribution correlates with that of H3K27me3. In particular, the methyl-CpG-binding domain of MeCP2 shows preferential interactions with H3K27me3. We further observe that the impact of MeCP2 on transcriptional changes correlates with histone post-translational modification patterns. Our findings indicate that MeCP2 interacts with genomic loci via binding to DNA as well as histones, and that interaction between MeCP2 and histone proteins plays a key role in gene expression regulation.

Nursing Genetic Research: New Insights Linking Breast Cancer Genetics and Bone Density

Nursing research is expected to provide options for the primary prevention of disease and health promotion, regardless of pathology or disease. Nurses have the skills to develop and lead research that addresses the relationship between genetic factors and health. Increasing genetic knowledge and research capacity through interdisciplinary cooperation as well as the development of research resources, will accelerate the rate at which nurses contribute to the knowledge about genetics and health. There are currently different fields in which knowledge can be expanded by research developed from the nursing field. Here, we present an emerging field of research in which it is hypothesized that genetics may affect bone metabolism. Better insight of genetic factors that are contributing to metabolic bone diseases would allow for focused nursing care and preventive interventions.

Type II Alexander disease caused by splicing errors and aberrant overexpression of an uncharacterized GFAP isoform

Alexander disease results from gain-of-function mutations in the gene encoding glial fibrillary acidic protein (GFAP). At least eight GFAP isoforms have been described, however, the predominant alpha isoform accounts for ∼90% of GFAP protein. We describe exonic variants identified in three unrelated families with Type II Alexander disease that alter the splicing of GFAP pre-messenger RNA (mRNA) and result in the upregulation of a previously uncharacterized GFAP lambda isoform (NM_001363846.1). Affected members of Family 1 and Family 2 shared the same missense variant, NM_001363846.1:c.1289G>A;p.(Arg430His) while in Family 3 we identified a synonymous variant in the adjacent nucleotide, NM_001363846.1:c.1290C>A;p.(Arg430Arg). Using RNA and protein analysis of brain autopsy samples, and a mini-gene splicing reporter assay, we demonstrate both variants result in the upregulation of the lambda isoform. Our approach demonstrates the importance of characterizing the effect of GFAP variants on mRNA splicing to inform future pathophysiologic and therapeutic study for Alexander disease.

An Independent Locus Upstream of ASIP Controls Variation in the Shade of the Bay Coat Colour in Horses

Novel coat colour phenotypes often emerge during domestication, and there is strong evidence of genetic selection for the two main genes that control base coat colour in horses-ASIP and MC1R. These genes direct the type of pigment produced, red pheomelanin (MC1R) or black eumelanin (ASIP), as well as the relative concentration and the temporal-spatial distribution of melanin pigment deposits in the skin and hair coat. Here, we describe a genome-wide association study (GWAS) to identify novel genic regions involved in the determination of the shade of bay. In total, 126 horses from five different breeds were ranked according to the extent of the distribution of eumelanin: spanning variation in phenotype from black colour restricted only to the extremities to the presence of some black pigment across nearly all the body surface. We identified a single region associated with the shade of bay ranking spanning approximately 0.5 MB on ECA22, just upstream of the ASIP gene (p = 9.76 × 10-15). This candidate region encompasses the distal 5' end of the ASIP transcript (as predicted from other species) as well as the RALY gene. Both loci are viable candidates based on the presence of similar alleles in other species. These results contribute to the growing understanding of coat colour genetics in the horse and to the mapping of genetic determinants of pigmentation on a molecular level. Given pleiotropic phenotypes in behaviour and obesity for ASIP alleles, especially those in the 5' regulatory region, improved understanding of this new Shade allele may have implications for health management in the horse.

Genome-Wide Association Meta-Analysis of Single-Nucleotide Polymorphisms and Symptomatic Venous Thromboembolism during Therapy for Acute Lymphoblastic Leukemia and Lymphoma in Caucasian Children

Symptomatic venous thromboembolism (VTE) occurs in five percent of children treated for acute lymphoblastic leukemia (ALL), but whether a genetic predisposition exists across different ALL treatment regimens has not been well studied.

METHODS

We undertook a genome-wide association study (GWAS) meta-analysis for VTE in consecutively treated children in the Nordic/Baltic acute lymphoblastic leukemia 2008 (ALL2008) cohort and the Australian Evaluation of Risk of ALL Treatment-Related Side-Effects (ERASE) cohort. A total of 92 cases and 1481 controls of European ancestry were included.

RESULTS

No SNPs reached genome-wide significance (p < 5 × 10-8) in either cohort. Among the top 34 single-nucleotide polymorphisms (SNPs) (p < 1 × 10-6), two loci had concordant effects in both cohorts: ALOX15B (rs1804772) (MAF: 1%; p = 3.95 × 10-7) that influences arachidonic acid metabolism and thus platelet aggregation, and KALRN (rs570684) (MAF: 1%; p = 4.34 × 10-7) that has been previously associated with risk of ischemic stroke, atherosclerosis, and early-onset coronary artery disease.

CONCLUSION

This represents the largest GWAS meta-analysis conducted to date associating SNPs to VTE in children and adolescents treated on childhood ALL protocols. Validation of these findings is needed and may then lead to patient stratification for VTE preventive interventions. As VTE hemostasis involves multiple pathways, a more powerful GWAS is needed to detect combination of variants associated with VTE.

Codon Pairs are Phylogenetically Conserved: A comprehensive analysis of codon pairing conservation across the Tree of Life

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques.

COXPRESdb v7: a gene coexpression database for 11 animal species supported by 23 coexpression platforms for technical evaluation and evolutionary inference

The advent of RNA-sequencing and microarray technologies has led to rapid growth of transcriptome data generated for a wide range of organisms, under various cellular, organ and individual conditions. Since the number of possible combinations of intercellular and extracellular conditions is almost unlimited, cataloging all transcriptome conditions would be an immeasurable challenge. Gene coexpression refers to the similarity of gene expression patterns under various conditions, such as disease states, tissue types, and developmental stages. Since the quality of gene coexpression data depends on the quality and quantity of transcriptome data, timely usage of the growing data is key to promoting individual research in molecular biology. COXPRESdb (http://coxpresdb.jp) is a database providing coexpression information for 11 animal species. One characteristic feature of COXPRESdb is its ability to compare multiple coexpression data derived from different transcriptomics technologies and different species, which strongly reduces false positive relationships in individual gene coexpression data. Here, we summarized the current version of this database, including 23 coexpression platforms with the highest-level quality till date. Using various functionalities in COXPRESdb, the new coexpression data would support a broader area of research from molecular biology to medical sciences.

REV-ERBα activates the mTOR signalling pathway and promotes myotubes differentiation

BACKGROUND INFORMATION

Mammalian target of rapamycin (mTOR) complex 1 (mTORC1) is a master regulator of cell and whole-body energy homoeostasis. REV-ERBα is a nuclear receptor that plays an important role in metabolism. While mTORC1 activation is necessary for muscle differentiation, the role of REV-ERBα is less clear.

RESULTS

We studied the effect of REV-ERBα overexpression and silencing as well as mTORC1 activation and inhibition on the differentiation of C2C12 myoblasts to myotubes. mTOR, myogenin and REV-ERBα were induced during differentiation of myoblasts into myotubes. REV-ERBα was found to activate mTORC1 during the differentiation process even in the absence of the differentiation medium. This activation was presumably through the downregulation of the expression of TSC1, an mTORC1 inhibitor.

CONCLUSION

Herein we show that REV-ERBα promotes myoblasts differentiation via the activation of the mTORC1 signalling pathway.

SIGNIFICANCE

REV-ERBα modulation can activate mTORC1 signalling and promote myoblasts differentiation.

DolphinNext: a distributed data processing platform for high throughput genomics

BACKGROUND

The emergence of high throughput technologies that produce vast amounts of genomic data, such as next-generation sequencing (NGS) is transforming biological research. The dramatic increase in the volume of data, the variety and continuous change of data processing tools, algorithms and databases make analysis the main bottleneck for scientific discovery. The processing of high throughput datasets typically involves many different computational programs, each of which performs a specific step in a pipeline. Given the wide range of applications and organizational infrastructures, there is a great need for highly parallel, flexible, portable, and reproducible data processing frameworks. Several platforms currently exist for the design and execution of complex pipelines. Unfortunately, current platforms lack the necessary combination of parallelism, portability, flexibility and/or reproducibility that are required by the current research environment. To address these shortcomings, workflow frameworks that provide a platform to develop and share portable pipelines have recently arisen. We complement these new platforms by providing a graphical user interface to create, maintain, and execute complex pipelines. Such a platform will simplify robust and reproducible workflow creation for non-technical users as well as provide a robust platform to maintain pipelines for large organizations.

RESULTS

To simplify development, maintenance, and execution of complex pipelines we created DolphinNext. DolphinNext facilitates building and deployment of complex pipelines using a modular approach implemented in a graphical interface that relies on the powerful Nextflow workflow framework by providing 1. A drag and drop user interface that visualizes pipelines and allows users to create pipelines without familiarity in underlying programming languages. 2. Modules to execute and monitor pipelines in distributed computing environments such as high-performance clusters and/or cloud 3. Reproducible pipelines with version tracking and stand-alone versions that can be run independently. 4. Modular process design with process revisioning support to increase reusability and pipeline development efficiency. 5. Pipeline sharing with GitHub and automated testing 6. Extensive reports with R-markdown and shiny support for interactive data visualization and analysis.

CONCLUSION

DolphinNext is a flexible, intuitive, web-based data processing and analysis platform that enables creating, deploying, sharing, and executing complex Nextflow pipelines with extensive revisioning and interactive reporting to enhance reproducible results.

FilTar: using RNA-Seq data to improve microRNA target prediction accuracy in animals

MOTIVATION

MicroRNA (miRNA) target prediction algorithms do not generally consider biological context and therefore generic target prediction based on seed binding can lead to a high level of false-positive predictions. Here, we present FilTar, a method that incorporates RNA-Seq data to make miRNA target prediction specific to a given cell type or tissue of interest.

RESULTS

We demonstrate that FilTar can be used to: (i) provide sample specific 3'-UTR reannotation; extending or truncating default annotations based on RNA-Seq read evidence and (ii) filter putative miRNA target predictions by transcript expression level, thus removing putative interactions where the target transcript is not expressed in the tissue or cell line of interest. We test the method on a variety of miRNA transfection datasets and demonstrate increased accuracy versus generic miRNA target prediction methods.

AVAILABILITY AND IMPLEMENTATION

FilTar is freely available and can be downloaded from https://github.com/TBradley27/FilTar. The tool is implemented using the Python and R programming languages, and is supported on GNU/Linux operating systems.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transient vitamin B5 starving improves mammalian cell homeostasis and protein production

Maintaining a metabolic steady state is essential for an organism's fitness and survival when confronted with environmental stress, and metabolic imbalance can be reversed by exposing the organism to fasting. Here, we attempted to apply this physiological principle to mammalian cell cultures to improve cellular fitness and consequently their ability to express recombinant proteins. We showed that transient vitamin B5 deprivation, an essential cofactor of central cellular metabolism, can quickly and irreversibly affect mammalian cell growth and division. A selection method was designed that relies on mammalian cell dependence on vitamin B5 for energy production, using the co-expression of the B5 transporter SLC5A6 and a gene of interest. We demonstrated that vitamin B5 selection persistently activates peroxisome proliferator-activated receptors (PPAR), a family of transcription factors involved in energy homeostasis, thereby altering lipid metabolism, improving cell fitness and therapeutic protein production. Thus, stable PPAR activation may constitute a cellular memory of past deprivation state, providing increased resistance to further potential fasting events. In other words, our results imply that cultured cells, once exposed to metabolic starvation, may display an improved metabolic fitness as compared to non-exposed cells, allowing increased resistance to cellular stress.

HIPPIE2: a method for fine-scale identification of physically interacting chromatin regions

Most regulatory chromatin interactions are mediated by various transcription factors (TFs) and involve physically interacting elements such as enhancers, insulators or promoters. To map these elements and interactions at a fine scale, we developed HIPPIE2 that analyzes raw reads from high-throughput chromosome conformation (Hi-C) experiments to identify precise loci of DNA physically interacting regions (PIRs). Unlike standard genome binning approaches (e.g. 10-kb to 1-Mb bins), HIPPIE2 dynamically infers the physical locations of PIRs using the distribution of restriction sites to increase analysis precision and resolution. We applied HIPPIE2 to in situ Hi-C datasets across six human cell lines (GM12878, IMR90, K562, HMEC, HUVEC, NHEK) with matched ENCODE/Roadmap functional genomic data. HIPPIE2 detected 1042 738 distinct PIRs, with high resolution (average PIR length of 1006 bp) and high reproducibility (92.3% in GM12878). PIRs are enriched for epigenetic marks (H3K27ac, H3K4me1) and open chromatin, suggesting active regulatory roles. HIPPIE2 identified 2.8 million significant PIR–PIR interactions, 27.2% of which were enriched for TF binding sites. 50 608 interactions were enhancer–promoter interactions and were enriched for 33 TFs, including known DNA looping/long-range mediators. These findings demonstrate that the novel dynamic approach of HIPPIE2 (https://bitbucket.com/wanglab-upenn/HIPPIE2) enables the characterization of chromatin and regulatory interactions with high resolution and reproducibility.

The genome-wide landscape of C:G > T:A polymorphism at the CpG contexts in the human population

BACKGROUND

The C:G > T:A substitution at the CpG dinucleotide contexts is the most frequent substitution type in genome evolution. The mutational process is obviously ongoing in the human germline; however, its impact on common and rare genomic polymorphisms has not been comprehensively investigated yet. Here we observed the landscape and dynamics of C:G > T:A substitutions from population-scale human genome sequencing datasets including ~ 4300 whole-genomes from the 1000 Genomes and the pan-cancer analysis of whole genomes (PCAWG) Project and ~ 60,000 whole-exomes from the Exome Aggregation Consortium (ExAC) database.

RESULTS

Of the 28,084,558 CpG sites in the human reference genome, 26.0% show C:G > T:A substitution in the dataset. Remarkably, CpGs in CpG islands (CGIs) have a much lower frequency of such mutations (5.6%). Interestingly, the mutation frequency of CGIs is not uniform with a significantly higher C:G > T:A substitution rate for intragenic CGIs compared to other types. For non-CGI CpGs, the mutation rate was positively correlated with the distance from the nearest CGI up to 2 kb. Finally, we found the impact of negative selection for coding CpG mutations resulting in amino acid change.

CONCLUSIONS

This study provides the first unbiased rate of C:G > T:A substitution at the CpG dinucleotide contexts, using population-scale human genome sequencing data. Our findings provide insights into the dynamics of the mutation acquisition in the human genome.

Formation of human long intergenic non-coding RNA genes, pseudogenes, and protein genes: Ancestral sequences are key players

Pathways leading to formation of non-coding RNA and protein genes are varied and complex. We report finding a conserved repeat sequence present in human and chimpanzee genomes that appears to have originated from a common primate ancestor. This sequence is repeatedly copied in human chromosome 22 (chr22) low copy repeats (LCR22) or segmental duplications and forms twenty-one different genes, which include the human long intergenic non-coding RNA (lincRNA) family FAM230, a newly discovered lincRNA gene family termed conserved long intergenic non-coding RNAs (clincRNA), pseudogene families, as well as the gamma-glutamyltransferase (GGT) protein gene family and the RNA pseudogenes that originate from GGT sequences. Of particular interest are the GGT5 and USP18 protein genes that appear to have formed from an homologous repeat sequence that also forms the clincRNA gene family. The data point to ancestral DNA sequences, conserved through evolution and duplicated in humans by chromosomal repeat sequences that may serve as functional genomic elements in the development of diverse genes.

Targeted sequencing of the LRRTM gene family in suicide attempters with bipolar disorder

Glutamatergic signaling is the primary excitatory neurotransmission pathway in the brain, and its relationship to neuropsychiatric disorders is of considerable interest. Our previous attempted suicide genome-wide association study, and numerous studies investigating gene expression, genetic variation, and DNA methylation have implicated aberrant glutamatergic signaling in suicide risk. The glutamatergic pathway gene LRRTM4 was an associated gene identified in our attempted suicide genome-wide association study, with association support seen primarily in females. Recent evidence has also shown that glutamatergic signaling is partly regulated by sex-related hormones. The LRRTM gene family encodes neuronal leucine-rich transmembrane proteins that localize to and promote glutamatergic synapse development. In this study, we sequenced the coding and regulatory regions of all four LRRTM gene members plus a large intronic region of LRRTM4 in 476 bipolar disorder suicide attempters and 473 bipolar disorder nonattempters. We identified two male-specific variants, one female- and five male-specific haplotypes significantly associated with attempted suicide in LRRTM4. Furthermore, variants within significant haplotypes may be brain expression quantitative trait loci for LRRTM4 and some of these variants overlap with predicted hormone response elements. Overall, these results provide supporting evidence for a sex-specific association of genetic variation in LRRTM4 with attempted suicide.

Twin enzymes, divergent control: The cholesterogenic enzymes DHCR14 and LBR are differentially regulated transcriptionally and post-translationally

Cholesterol synthesis is a tightly regulated process, both transcriptionally and post-translationally. Transcriptional control of cholesterol synthesis is relatively well-understood. However, of the ∼20 enzymes in cholesterol biosynthesis, post-translational regulation has only been examined for a small number. Three of the four sterol reductases in cholesterol production, 7-dehydrocholesterol reductase (DHCR7), 14-dehydrocholesterol reductase (DHCR14), and lamin-B receptor (LBR), share evolutionary ties with a high level of sequence homology and predicted structural homology. DHCR14 and LBR uniquely share the same Δ-14 reductase activity in cholesterol biosynthesis, yet little is known about their post-translational regulation. We have previously identified specific modes of post-translational control of DHCR7, but it is unknown whether these regulatory mechanisms are shared by DHCR14 and LBR. Using CHO-7 cells stably expressing epitope-tagged DHCR14 or LBR, we investigated the post-translational regulation of these enzymes. We found that DHCR14 and LBR undergo differential post-translational regulation, with DHCR14 being rapidly turned over, triggered by cholesterol and other sterol intermediates, whereas LBR remained stable. DHCR14 is degraded via the ubiquitin-proteasome system, and we identified several DHCR14 and DHCR7 putative interaction partners, including a number of E3 ligases that modulate DHCR14 levels. Interestingly, we found that gene expression across an array of human tissues showed a negative relationship between the C14-sterol reductases; one enzyme or the other tends to be predominantly expressed in each tissue. Overall, our findings indicate that whereas LBR tends to be the constitutively active C14-sterol reductase, DHCR14 levels are tunable, responding to the local cellular demands for cholesterol.

Integrated omics profiling reveals novel patterns of epigenetic programming in cancer-associated myofibroblasts

There is increasing evidence that stromal myofibroblasts play a key role in the tumour development however, the mechanisms by which they become reprogrammed to assist in cancer progression remain unclear. As cultured cancer-associated myofibroblasts (CAMs) retain an ability to enhance the proliferation and migration of cancer cells in vitro, it is possible that epigenetic reprogramming of CAMs within the tumour microenvironment may confer long-term pro-tumourigenic changes in gene expression. This study reports the first comparative multi-omics analysis of cancer-related changes in gene expression and DNA methylation in primary myofibroblasts derived from gastric and oesophageal tumours. In addition, we identify novel CAM-specific DNA methylation signatures, which are not observed in patient-matched adjacent tissue-derived myofibroblasts, or corresponding normal tissue-derived myofibroblasts. Analysis of correlated changes in DNA methylation and gene expression shows that different patterns of gene-specific DNA methylation have the potential to confer pro-tumourigenic changes in metabolism, cell signalling and differential responses to hypoxia. These molecular signatures provide new insights into potential mechanisms of stromal reprogramming in gastric and oesophageal cancer, while also providing a new resource to facilitate biomarker identification and future hypothesis-driven studies into mechanisms of stromal reprogramming and tumour progression in solid tumours.

The Clinical Genome and Ancestry Report: An interactive web application for prioritizing clinically implicated variants from genome sequencing data with ancestry composition

Genome sequencing is positioned as a routine clinical work-up for diverse clinical conditions. A commonly used approach to highlight candidate variants with potential clinical implication is to search over locus- and gene-centric knowledge databases. Most web-based applications allow a federated query across diverse databases for a single variant; however, sifting through a large number of genomic variants with combination of filtering criteria is a substantial challenge. Here we describe the Clinical Genome and Ancestry Report (CGAR), an interactive web application developed to follow clinical interpretation workflows by organizing variants into seven categories: (1) reported disease-associated variants, (2) rare- and high-impact variants in putative disease-associated genes, (3) secondary findings that the American College of Medical Genetics and Genomics recommends reporting back to patients, (4) actionable pharmacogenomic variants, (5) focused reports for candidate genes, (6) de novo variant candidates for trio analysis, and (7) germline and somatic variants implicated in cancer risk, diagnosis, treatment and prognosis. For each variant, a comprehensive list of external links to variant-centric and phenotype databases are provided. Furthermore, genotype-derived ancestral composition is used to highlight allele frequencies from a matched population since some disease-associated variants show a wide variation between populations. CGAR is an open-source software and is available at https://tom.tch.harvard.edu/apps/cgar/.