The FANTOM5 collection, a data series underpinning mammalian transcriptome atlases in diverse cell types

Disrupting the RNA polymerase II transcription cycle through CDK7 inhibition ameliorates inflammatory arthritis

Macrophages are key drivers of inflammation and tissue damage in autoimmune diseases including rheumatoid arthritis. The rate-limiting step for transcription of more than 70% of inducible genes in macrophages is RNA polymerase II (Pol II) promoter-proximal pause release; however, the specific role of Pol II early elongation control in inflammation, and whether it can be modulated therapeutically, is unknown. Genetic ablation of a pause-stabilizing negative elongation factor (NELF) in macrophages did not affect baseline Pol II occupancy but enhanced the transcriptional response of paused anti-inflammatory genes to lipopolysaccharide followed by secondary attenuation of inflammatory signaling in vitro and in the K/BxN serum transfer mouse model of arthritis. To pharmacologically disrupt the Pol II transcription cycle, we used two covalent inhibitors of the transcription factor II H-associated cyclin-dependent kinase 7 (CDK7), THZ1 and YKL-5-124. Both reduced Pol II pausing in murine and human macrophages, broadly suppressed induction of pro- but not anti-inflammatory genes, and rapidly reversed preestablished inflammatory macrophage polarization. In mice, CDK7 inhibition ameliorated both acute and chronic progressive inflammatory arthritis. Lastly, CDK7 inhibition down-regulated a pathogenic gene expression signature in synovial explants from patients with rheumatoid arthritis. We propose that interfering with Pol II early elongation by targeting CDK7 represents a therapeutic opportunity for rheumatoid arthritis and other inflammatory diseases.

Latent space arithmetic on data embeddings from healthy multi-tissue human RNA-seq decodes disease modules

Computational analyses of transcriptomic data have dramatically improved our understanding of complex diseases. However, such approaches are limited by small sample sets of disease-affected material. We asked if a variational autoencoder trained on large groups of healthy human RNA sequencing (RNA-seq) data can capture the fundamental gene regulation system and generalize to unseen disease changes. Importantly, we found this model to successfully compress unseen transcriptomic changes from 25 independent disease datasets. We decoded disease-specific signals from the latent space and found them to contain more disease-specific genes than the corresponding differential expression analysis in 20 of 25 cases. Finally, we matched these disease signals with known drug targets and extracted sets of known and potential pharmaceutical candidates. In summary, our study demonstrates how data-driven representation learning enables the arithmetic deconstruction of the latent space, facilitating the dissection of disease mechanisms and drug targets.

Mechanisms of epigenomic and functional convergence between glucocorticoid- and IL4-driven macrophage programming

Macrophages adopt distinct phenotypes in response to environmental cues, with type-2 cytokine interleukin-4 promoting a tissue-repair homeostatic state (M2IL4). Glucocorticoids (GC), widely used anti-inflammatory therapeutics, reportedly impart a similar phenotype (M2GC), but how such disparate pathways may functionally converge is unknown. We show using integrative functional genomics that M2IL4 and M2GC transcriptomes share a striking overlap mirrored by a shift in chromatin landscape in both common and signal-specific gene subsets. This core homeostatic program is enacted by transcriptional effectors KLF4 and the glucocorticoid receptor, whose genome-wide occupancy and actions are integrated in a stimulus-specific manner by the nuclear receptor cofactor GRIP1. Indeed, many of the M2IL4:M2GC-shared transcriptomic changes were GRIP1-dependent. Consistently, GRIP1 loss attenuated phagocytic activity of both populations in vitro and macrophage tissue-repair properties in the murine colitis model in vivo. These findings provide a mechanistic framework for homeostatic macrophage programming by distinct signals, to better inform anti-inflammatory drug design.

Perspectives in the investigation of Cockayne syndrome group B neurological disease: the utility of patient-derived brain organoid models

Cockayne syndrome (CS) group B (CSB), which results from mutations in the excision repair cross-complementation group 6 (ERCC6) genes, which produce CSB protein, is an autosomal recessive disease characterized by multiple progressive disorders including growth failure, microcephaly, skin photosensitivity, and premature aging. Clinical data show that brain atrophy, demyelination, and calcification are the main neurological manifestations of CS, which progress with time. Neuronal loss and calcification occur in various brain areas, particularly the cerebellum and basal ganglia, resulting in dyskinesia, ataxia, and limb tremors in CSB patients. However, the understanding of neurodevelopmental defects in CS has been constrained by the lack of significant neurodevelopmental and functional abnormalities observed in CSB-deficient mice. In this review, we focus on elucidating the protein structure and distribution of CSB and delve into the impact of CSB mutations on the development and function of the nervous system. In addition, we provide an overview of research models that have been instrumental in exploring CS disorders, with a forward-looking perspective on the substantial contributions that brain organoids are poised to further advance this field.

Decoding bioactive signals of the RNA secretome: the cell-free messenger RNA catalogue

Despite gene-expression profiling being one of the most common methods to evaluate molecular dysregulation in tissues, the utilization of cell-free messenger RNA (cf-mRNA) as a blood-based non-invasive biomarker analyte has been limited compared to other RNA classes. Recent advancements in low-input RNA-sequencing and normalization techniques, however, have enabled characterization as well as accurate quantification of cf-mRNAs allowing direct pathological insights. The molecular profile of the cell-free transcriptome in multiple diseases has subsequently been characterized including, prenatal diseases, neurological disorders, liver diseases and cancers suggesting this biological compartment may serve as a disease agnostic platform. With mRNAs packaged in a myriad of extracellular vesicles and particles, these signals may be used to develop clinically actionable, non-invasive disease biomarkers. Here, we summarize the recent scientific developments of extracellular mRNA, biology of extracellular mRNA carriers, clinical utility of cf-mRNA as disease biomarkers, as well as proposed functions in cell and tissue pathophysiology.

Mechanisms of Epigenomic and Functional Convergence Between Glucocorticoid- and IL4-Driven Macrophage Programming

Macrophages adopt distinct phenotypes in response to environmental cues, with type-2 cytokine interleukin-4 promoting a tissue-repair homeostatic state (M2IL4). Glucocorticoids, widely used anti-inflammatory therapeutics, reportedly impart a similar phenotype (M2GC), but how such disparate pathways may functionally converge is unknown. We show using integrative functional genomics that M2IL4 and M2GC transcriptomes share a striking overlap mirrored by a shift in chromatin landscape in both common and signal-specific gene subsets. This core homeostatic program is enacted by transcriptional effectors KLF4 and the GC receptor, whose genome-wide occupancy and actions are integrated in a stimulus-specific manner by the nuclear receptor cofactor GRIP1. Indeed, many of the M2IL4:M2GC-shared transcriptomic changes were GRIP1-dependent. Consistently, GRIP1 loss attenuated phagocytic activity of both populations in vitro and macrophage tissue-repair properties in the murine colitis model in vivo. These findings provide a mechanistic framework for homeostatic macrophage programming by distinct signals, to better inform anti-inflammatory drug design.

Major urinary protein (Mup) gene family deletion drives sex-specific alterations in the house-mouse gut microbiota

The gut microbiota is shaped by host metabolism. In house mice (Mus musculus), major urinary protein (MUP) pheromone production represents a considerable energy investment, particularly in sexually mature males. Deletion of the Mup gene family shifts mouse metabolism toward an anabolic state, marked by lipogenesis, lipid accumulation, and body mass increases. Given the metabolic implications of MUPs, they may also influence the gut microbiota. Here, we investigated the effect of a deletion of the Mup gene family on the gut microbiota of sexually mature mice. Shotgun metagenomics revealed distinct taxonomic and functional profiles between wild-type and knockout males but not females. Deletion of the Mup gene cluster significantly reduced diversity in microbial families and functions in male mice. Additionally, a species of Ruminococcaceae and several microbial functions, such as transporters involved in vitamin B5 acquisition, were significantly depleted in the microbiota of Mup knockout males. Altogether, these results show that MUPs significantly affect the gut microbiota of house mouse in a sex-specific manner.IMPORTANCEThe community of microorganisms that inhabits the gastrointestinal tract can have profound effects on host phenotypes. The gut microbiota is in turn shaped by host genes, including those involved with host metabolism. In adult male house mice, expression of the major urinary protein (Mup) gene cluster represents a substantial energy investment, and deletion of the Mup gene family leads to fat accumulation and weight gain in males. We show that deleting Mup genes also alters the gut microbiota of male, but not female, mice in terms of both taxonomic and functional compositions. Male mice without Mup genes harbored fewer gut bacterial families and reduced abundance of a species of Ruminococcaceae, a family that has been previously shown to reduce obesity risk. Studying the impact of the Mup gene family on the gut microbiota has the potential to reveal the ways in which these genes affect host phenotypes.

Pig-eRNAdb: a comprehensive enhancer and eRNA dataset of pigs

Enhancers and the enhancer RNAs (eRNAs) have been strongly implicated in regulations of transcriptions. Based the multi-omics data (ATAC-seq, ChIP-seq and RNA-seq) from public databases, Pig-eRNAdb is a dataset that comprehensively integrates enhancers and eRNAs for pigs using the machine learning strategy, which incorporates 82,399 enhancers and 37,803 eRNAs from 607 samples across 15 tissues of pigs. This user-friendly dataset covers a comprehensive depth of enhancers and eRNAs annotation for pigs. The coordinates of enhancers and the expression patterns of eRNAs are downloadable. Besides, thousands of regulators on eRNAs, the target genes of eRNAs, the tissue-specific eRNAs, and the housekeeping eRNAs are also accessible as well as the sequence similarity of eRNAs with humans. Moreover, the tissue-specific eRNA-trait associations encompass 652 traits are also provided. It will crucially facilitate investigations on enhancers and eRNAs with Pig-eRNAdb as a reference dataset in pigs.

MethMarkerDB: a comprehensive cancer DNA methylation biomarker database

DNA methylation plays a crucial role in tumorigenesis and tumor progression, sparking substantial interest in the clinical applications of cancer DNA methylation biomarkers. Cancer-related whole-genome bisulfite sequencing (WGBS) data offers a promising approach to precisely identify these biomarkers with differentially methylated regions (DMRs). However, currently there is no dedicated resource for cancer DNA methylation biomarkers with WGBS data. Here, we developed a comprehensive cancer DNA methylation biomarker database (MethMarkerDB, https://methmarkerdb.hzau.edu.cn/), which integrated 658 WGBS datasets, incorporating 724 curated DNA methylation biomarker genes from 1425 PubMed published articles. Based on WGBS data, we documented 5.4 million DMRs from 13 common types of cancer as candidate DNA methylation biomarkers. We provided search and annotation functions for these DMRs with different resources, such as enhancers and SNPs, and developed diagnostic and prognostic models for further biomarker evaluation. With the database, we not only identified known DNA methylation biomarkers, but also identified 781 hypermethylated and 5245 hypomethylated pan-cancer DMRs, corresponding to 693 and 2172 genes, respectively. These novel potential pan-cancer DNA methylation biomarkers hold significant clinical translational value. We hope that MethMarkerDB will help identify novel cancer DNA methylation biomarkers and propel the clinical application of these biomarkers.

Intra-promoter switch of transcription initiation sites in proliferation signaling-dependent RNA metabolism

Global changes in transcriptional regulation and RNA metabolism are crucial features of cancer development. However, little is known about the role of the core promoter in defining transcript identity and post-transcriptional fates, a potentially crucial layer of transcriptional regulation in cancer. In this study, we use CAGE-seq analysis to uncover widespread use of dual-initiation promoters in which non-canonical, first-base-cytosine (C) transcription initiation occurs alongside first-base-purine initiation across 59 human cancers and healthy tissues. C-initiation is often followed by a 5' terminal oligopyrimidine (5'TOP) sequence, dramatically increasing the range of genes potentially subjected to 5'TOP-associated post-transcriptional regulation. We show selective, dynamic switching between purine and C-initiation site usage, indicating transcription initiation-level regulation in cancers. We additionally detail global metabolic changes in C-initiation transcripts that mark differentiation status, proliferative capacity, radiosensitivity, and response to irradiation and to PI3K-Akt-mTOR and DNA damage pathway-targeted radiosensitization therapies in colorectal cancer organoids and cancer cell lines and tissues.

The mycosis fungoides cutaneous microenvironment shapes dysfunctional cell trafficking, antitumor immunity, matrix interactions, and angiogenesis

Malignant T lymphocyte proliferation in mycosis fungoides (MF) is largely restricted to the skin, implying that malignant cells are dependent on their specific cutaneous tumor microenvironment (TME), including interactions with non-malignant immune and stromal cells, cytokines, and other immunomodulatory factors. To explore these interactions, we performed a comprehensive transcriptome analysis of the TME in advanced-stage MF skin tumors by single-cell RNA sequencing. Our analysis identified cell-type compositions, cellular functions, and cell-to-cell interactions in the MF TME that were distinct from those from healthy skin and benign dermatoses. While patterns of gene expression were common among patient samples, high transcriptional diversity was also observed in immune and stromal cells, with dynamic interactions and crosstalk between these cells and malignant T lymphocytes. This heterogeneity mapped to processes such as cell trafficking, matrix interactions, angiogenesis, immune functions, and metabolism that affect cancer cell growth, migration, and invasion, as well as antitumor immunity. By comprehensively characterizing the transcriptomes of immune and stromal cells within the cutaneous microenvironment of individual MF tumors, we have identified patterns of dysfunction common to all tumors that represent a resource for identifying candidates with therapeutic potential as well as patient-specific heterogeneity that has important implications for personalized disease management.

Broad H3K4me3 domains: Maintaining cellular identity and their implication in super-enhancer hijacking

The human and mouse genomes are complex from a genomic standpoint. Each cell has the same genomic sequence, yet a wide array of cell types exists due to the presence of a plethora of regulatory elements in the non-coding genome. Recent advances in epigenomic profiling have uncovered non-coding gene proximal promoters and distal enhancers of transcription genome-wide. Extension of promoter-associated H3K4me3 histone mark across the gene body, known as a broad H3K4me3 domain (H3K4me3-BD), is a signature of constitutive expression of cell-type-specific regulation and of tumour suppressor genes in healthy cells. Recently, it has been discovered that the presence of H3K4me3-BDs over oncogenes is a cancer-specific feature associated with their dysregulated gene expression and tumourigenesis. Moreover, it has been shown that the hijacking of clusters of enhancers, known as super-enhancers (SE), by proto-oncogenes results in the presence of H3K4me3-BDs over the gene body. Therefore, H3K4me3-BDs and SE crosstalk in healthy and cancer cells therefore represents an important mechanism to identify future treatments for patients with SE driven cancers.

The status of the human gene catalogue

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

Clorfl86/RHEX Is a Negative Regulator of SCF/KIT Signaling in Human Skin Mast Cells

Mast cells (MCs) are key effector cells in allergic and inflammatory diseases, and the SCF/KIT axis regulates most aspects of the cells' biology. Using terminally differentiated skin MCs, we recently reported on proteome-wide phosphorylation changes initiated by KIT dimerization. C1orf186/RHEX was revealed as one of the proteins to become heavily phosphorylated. Its function in MCs is undefined and only some information is available for erythroblasts. Using public databases and our own data, we now report that RHEX exhibits highly restricted expression with a clear dominance in MCs. While expression is most pronounced in mature MCs, RHEX is also abundant in immature/transformed MC cell lines (HMC-1, LAD2), suggesting early expression with further increase during differentiation. Using RHEX-selective RNA interference, we reveal that RHEX unexpectedly acts as a negative regulator of SCF-supported skin MC survival. This finding is substantiated by RHEX's interference with KIT signal transduction, whereby ERK1/2 and p38 both were more strongly activated when RHEX was attenuated. Comparing RHEX and capicua (a recently identified repressor) revealed that each protein preferentially suppresses other signaling modules elicited by KIT. Induction of immediate-early genes strictly requires ERK1/2 in SCF-triggered MCs; we now demonstrate that RHEX diminution translates to this downstream event, and thereby enhances NR4A2, JUNB, and EGR1 induction. Collectively, our study reveals RHEX as a repressor of KIT signaling and function in MCs. As an abundant and selective lineage marker, RHEX may have various roles in the lineage, and the provided framework will enable future work on its involvement in other crucial processes.

The status of the human gene catalogue

Scientists have been trying to identify all of the genes in the human genome since the initial draft of the genome was published in 2001. Over the intervening years, much progress has been made in identifying protein-coding genes, and the estimated number has shrunk to fewer than 20,000, although the number of distinct protein-coding isoforms has expanded dramatically. The invention of high-throughput RNA sequencing and other technological breakthroughs have led to an explosion in the number of reported non-coding RNA genes, although most of them do not yet have any known function. A combination of recent advances offers a path forward to identifying these functions and towards eventually completing the human gene catalogue. However, much work remains to be done before we have a universal annotation standard that includes all medically significant genes, maintains their relationships with different reference genomes, and describes clinically relevant genetic variants.

Lentiviral gene therapy for X-linked chronic granulomatous disease recapitulates endogenous CYBB regulation and expression

X-linked chronic granulomatous disease (X-CGD) is a primary immunodeficiency caused by mutations in the CYBB gene, resulting in the inability of phagocytic cells to eliminate infections. To design a lentiviral vector (LV) capable of recapitulating the endogenous regulation and expression of CYBB, a bioinformatics-guided approach was used to elucidate the cognate enhancer elements regulating the native CYBB gene. Using this approach, we analyzed a 600-kilobase topologically associated domain of the CYBB gene and identified endogenous enhancer elements to supplement the CYBB promoter to develop MyeloVec, a physiologically regulated LV for the treatment of X-CGD. When compared with an LV currently in clinical trials for X-CGD, MyeloVec showed improved expression, superior gene transfer to hematopoietic stem and progenitor cells (HSPCs), corrected an X-CGD mouse model leading to complete protection against Burkholderia cepacia infection, and restored healthy donor levels of antimicrobial oxidase activity in neutrophils derived from HSPCs from patients with X-CGD. Our findings validate the bioinformatics-guided design approach and have yielded a novel LV with clinical promise for the treatment of X-CGD.

Exploring the potential mechanisms of impairment on genitourinary system associated with coronavirus disease 2019 infection: Bioinformatics and molecular simulation analyses

Objective

The novel coronavirus (severe acute respiratory syndrome coronavirus 2) has been spreading worldwide since December 2019, posing a serious danger to human health and socioeconomic development. A large number of clinical trials have revealed that coronavirus disease 2019 (COVID-19) results in multi-organ damage including the urogenital system. This study aimed to explore the potential mechanisms of genitourinary damage associated with COVID-19 infection through bioinformatics and molecular simulation analysis.

Methods

We used multiple publicly available databases to explore the expression patterns of ACE2, TMPRSS2, and CD147 (Basigin [BSG]) in major organs in the healthy and disease-specific populations, particularly the genitourinary organs. Single-cell RNA sequencing was used to analyze the cell-specific expression patterns of ACE2, TMPRSS2, CD147, cytokine receptors, and cytokine interacting proteins in genitourinary organs, such as the bladder, kidney, prostate, and testis. Additionally, gene set enrichment analysis was used to investigate the relationship between testosterone levels and COVID-19 vulnerability in patients with prostate cancer.

Results

The results revealed that ACE2, TMPRSS2, and CD147 were highly expressed in normal urogenital organs. Then, they were also highly expressed in multiple tumors and chronic kidney diseases. Additionally, ACE2, TMPRSS2, and CD147 were significantly expressed in a range of cells in urogenital organs according to single-cell RNA sequencing. Cytokine receptors and cytokine interacting proteins, especially CCL2, JUN, and TIMP1, were commonly highly expressed in urogenital organs. Finally, gene set enrichment analysis results showed that high testosterone levels in prostate cancer patients were significantly related to the JAK/STAT signaling pathway and the Toll-like receptor signaling pathway which were associated with COVID-19.

Conclusion

Our study provides new insights into the potential mechanisms of severe acute respiratory syndrome coronavirus 2 damage to urogenital organs from multiple perspectives, which may draw the attention of urologists to COVID-19 and contribute to the development of targeted drugs.

ChromLoops: a comprehensive database for specific protein-mediated chromatin loops in diverse organisms

Chromatin loops (or chromatin interactions) are important elements of chromatin structures. Disruption of chromatin loops is associated with many diseases, such as cancer and polydactyly. A few methods, including ChIA-PET, HiChIP and PLAC-Seq, have been proposed to detect high-resolution, specific protein-mediated chromatin loops. With rapid progress in 3D genomic research, ChIA-PET, HiChIP and PLAC-Seq datasets continue to accumulate, and effective collection and processing for these datasets are urgently needed. Here, we developed a comprehensive, multispecies and specific protein-mediated chromatin loop database (ChromLoops, https://3dgenomics.hzau.edu.cn/chromloops), which integrated 1030 ChIA-PET, HiChIP and PLAC-Seq datasets from 13 species, and documented 1 491 416 813 high-quality chromatin loops. We annotated genes and regions overlapping with chromatin loop anchors with rich functional annotations, such as regulatory elements (enhancers, super-enhancers and silencers), variations (common SNPs, somatic SNPs and eQTLs), and transcription factor binding sites. Moreover, we identified genes with high-frequency chromatin interactions in the collected species. In particular, we identified genes with high-frequency interactions in cancer samples. We hope that ChromLoops will provide a new platform for studying chromatin interaction regulation in relation to biological processes and disease.

Differential nuclear import sets the timing of protein access to the embryonic genome

The development of a fertilized egg to an embryo requires the proper temporal control of gene expression. During cell differentiation, timing is often controlled via cascades of transcription factors (TFs). However, in early development, transcription is often inactive, and many TF levels stay constant, suggesting that alternative mechanisms govern the observed rapid and ordered onset of gene expression. Here, we find that in early embryonic development access of maternally deposited nuclear proteins to the genome is temporally ordered via importin affinities, thereby timing the expression of downstream targets. We quantify changes in the nuclear proteome during early development and find that nuclear proteins, such as TFs and RNA polymerases, enter the nucleus sequentially. Moreover, we find that the timing of nuclear proteins' access to the genome corresponds to the timing of downstream gene activation. We show that the affinity of proteins to importin is a major determinant in the timing of protein entry into embryonic nuclei. Thus, we propose a mechanism by which embryos encode the timing of gene expression in early development via biochemical affinities. This process could be critical for embryos to organize themselves before deploying the regulatory cascades that control cell identities.

A systematic evaluation of the computational tools for ligand-receptor-based cell-cell interaction inference

Cell-cell interactions (CCIs) are essential for multicellular organisms to coordinate biological processes and functions. One classical type of CCI interaction is between secreted ligands and cell surface receptors, i.e. ligand-receptor (LR) interactions. With the recent development of single-cell technologies, a large amount of single-cell ribonucleic acid (RNA) sequencing (scRNA-Seq) data has become widely available. This data availability motivated the single-cell-resolution study of CCIs, particularly LR-based CCIs. Dozens of computational methods and tools have been developed to predict CCIs by identifying LR-based CCIs. Many of these tools have been theoretically reviewed. However, there is little study on current LR-based CCI prediction tools regarding their performance and running results on public scRNA-Seq datasets. In this work, to fill this gap, we tested and compared nine of the most recent computational tools for LR-based CCI prediction. We used 15 well-studied scRNA-Seq samples that correspond to approximately 100K single cells under different experimental conditions for testing and comparison. Besides briefing the methodology used in these nine tools, we summarized the similarities and differences of these tools in terms of both LR prediction and CCI inference between cell types. We provided insight into using these tools to make meaningful discoveries in understanding cell communications.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

Bioinformatics network analyses of growth differentiation factor 11

Growth differentiation factor 11 (GDF11) has been implicated in rejuvenating functions in age-related diseases. The molecular mechanisms connecting GDF11 with these anti-aging phenomena, including reverse age-related cardiac hypertrophy and vascular and neurogenic rejuvenation, remain unclear. In this study, we sought to uncover the molecular functions of GDF11 using bioinformatics and network-driven analyses at the human gene and transcription levels using the gene co-expression network analysis, the protein–protein interaction network analysis, and the transcription factor network analysis. Our findings suggested that GDF11 is involved in a variety of functions, such as apoptosis, DNA repair, telomere maintenance, and interaction with key transcription factors, such as MYC proto-oncogene, specificity protein 1, and ETS proto-oncogene 2. The human skin fibroblast premature senescence model was established by UVB. The treatment with 10 ng/mL GDF11 in this cell model could reduce cell damage, reduce the apoptosis rate and the expression of caspase-3, and increase the length of telomeres. Therefore, our findings shed light on the functions of GDF11 and provide insights into the roles of GDF11 in aging.

Genome-wide annotation of protein-coding genes in pig

BACKGROUND

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

RESULTS

An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented.

CONCLUSIONS

Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Transcriptome screening followed by integrated physicochemical and structural analyses for investigating RNA-mediated berberine activity

Non-coding RNAs are regarded as promising targets for the discovery of innovative drugs due to their abundance in the genome and their involvement in many biological processes. Phytochemicals (PCs) are the primary source of ligand-based drugs due to their broad spectrum of biological activities. Since many PCs are heterocyclic and have chemical groups potentially involved in the interaction with nucleic acids, detailed interaction analysis between PCs and RNA is crucial to explore the effect of PCs on RNA functions. In this study, an integrated approach for investigating interactions between PCs and RNAs were demonstrated to verify the RNA-mediated PCs functions by using berberine (BRB) as a model PC. RNA screening of a transcriptome library followed by sequence refinement found minimal RNA motif consisting of a cytosine bulge with U-A and G-U neighbouring base pairs for interaction with BRB. NMR-based structure determination and physicochemical analyses using chemical analogues of BRB demonstrated the importance of electrostatic and stacking interactions for sequence selective interaction and RNA stabilization. The selective interaction with a relatively small RNA motif based on a chemical structure of a planer heterocyclic highlights the biological activities of various PCs mediated by the interactions with particular functional RNAs. In addition, the systematic and quantitative investigations demonstrated in this study could be useful for the development of therapeutic chemicals targeting functional RNAs, based on the PCs, in the future.

Positive Selection in Gene Regulatory Factors Suggests Adaptive Pleiotropic Changes During Human Evolution

Gene regulatory factors (GRFs), such as transcription factors, co-factors and histone-modifying enzymes, play many important roles in modifying gene expression in biological processes. They have also been proposed to underlie speciation and adaptation. To investigate potential contributions of GRFs to primate evolution, we analyzed GRF genes in 27 publicly available primate genomes. Genes coding for zinc finger (ZNF) proteins, especially ZNFs with a Krüppel-associated box (KRAB) domain were the most abundant TFs in all genomes. Gene numbers per TF family differed between all species. To detect signs of positive selection in GRF genes we investigated more than 3,000 human GRFs with their more than 70,000 orthologs in 26 non-human primates. We implemented two independent tests for positive selection, the branch-site-model of the PAML suite and aBSREL of the HyPhy suite, focusing on the human and great ape branch. Our workflow included rigorous procedures to reduce the number of false positives: excluding distantly similar orthologs, manual corrections of alignments, and considering only genes and sites detected by both tests for positive selection. Furthermore, we verified the candidate sites for selection by investigating their variation within human and non-human great ape population data. In order to approximately assign a date to positively selected sites in the human lineage, we analyzed archaic human genomes. Our work revealed with high confidence five GRFs that have been positively selected on the human lineage and one GRF that has been positively selected on the great ape lineage. These GRFs are scattered on different chromosomes and have been previously linked to diverse functions. For some of them a role in speciation and/or adaptation can be proposed based on the expression pattern or association with human diseases, but it seems that they all contributed independently to human evolution. Four of the positively selected GRFs are KRAB-ZNF proteins, that induce changes in target genes co-expression and/or through arms race with transposable elements. Since each positively selected GRF contains several sites with evidence for positive selection, we suggest that these GRFs participated pleiotropically to phenotypic adaptations in humans.

LanCLs add glutathione to dehydroamino acids generated at phosphorylated sites in the proteome

Enzyme-mediated damage repair or mitigation, while common for nucleic acids, is rare for proteins. Examples of protein damage are elimination of phosphorylated Ser/Thr to dehydroalanine/dehydrobutyrine (Dha/Dhb) in pathogenesis and aging. Bacterial LanC enzymes use Dha/Dhb to form carbon-sulfur linkages in antimicrobial peptides, but the functions of eukaryotic LanC-like (LanCL) counterparts are unknown. We show that LanCLs catalyze the addition of glutathione to Dha/Dhb in proteins, driving irreversible C-glutathionylation. Chemo-enzymatic methods were developed to site-selectively incorporate Dha/Dhb at phospho-regulated sites in kinases. In human MAPK-MEK1, such "elimination damage" generated aberrantly activated kinases, which were deactivated by LanCL-mediated C-glutathionylation. Surveys of endogenous proteins bearing damage from elimination (the eliminylome) also suggest it is a source of electrophilic reactivity. LanCLs thus remove these reactive electrophiles and their potentially dysregulatory effects from the proteome. As knockout of LanCL in mice can result in premature death, repair of this kind of protein damage appears important physiologically.

Fish-Ing for Enhancers in the Heart

Precise control of gene expression is crucial to ensure proper development and biological functioning of an organism. Enhancers are non-coding DNA elements which play an essential role in regulating gene expression. They contain specific sequence motifs serving as binding sites for transcription factors which interact with the basal transcription machinery at their target genes. Heart development is regulated by intricate gene regulatory network ensuring precise spatiotemporal gene expression program. Mutations affecting enhancers have been shown to result in devastating forms of congenital heart defect. Therefore, identifying enhancers implicated in heart biology and understanding their mechanism is key to improve diagnosis and therapeutic options. Despite their crucial role, enhancers are poorly studied, mainly due to a lack of reliable way to identify them and determine their function. Nevertheless, recent technological advances have allowed rapid progress in enhancer discovery. Model organisms such as the zebrafish have contributed significant insights into the genetics of heart development through enabling functional analyses of genes and their regulatory elements in vivo. Here, we summarize the current state of knowledge on heart enhancers gained through studies in model organisms, discuss various approaches to discover and study their function, and finally suggest methods that could further advance research in this field.

Genetic basis of mitochondrial diseases

Mitochondrial disorders are monogenic disorders characterized by a defect in oxidative phosphorylation and caused by pathogenic variants in one of over 340 different genes. The implementation of whole-exome sequencing has led to a revolution in their diagnosis, duplicated the number of associated disease genes, and significantly increased the diagnosed fraction. However, the genetic etiology of a substantial fraction of patients exhibiting mitochondrial disorders remains unknown, highlighting limitations in variant detection and interpretation, which calls for improved computational and DNA sequencing methods, as well as the addition of OMICS tools. More intriguingly, this also suggests that some pathogenic variants lie outside of the protein-coding genes and that the mechanisms beyond the Mendelian inheritance and the mtDNA are of relevance. This review covers the current status of the genetic basis of mitochondrial diseases, discusses current challenges and perspectives, and explores the contribution of factors beyond the protein-coding regions and monogenic inheritance in the expansion of the genetic spectrum of disease.

CTD Anatomy: analyzing chemical-induced phenotypes and exposures from an anatomical perspective, with implications for environmental health studies

The Comparative Toxicogenomics Database (CTD) is a freely available public resource that curates and interrelates chemical, gene/protein, phenotype, disease, organism, and exposure data. CTD can be used to address toxicological mechanisms for environmental chemicals and facilitate the generation of testable hypotheses about how exposures affect human health. At CTD, manually curated interactions for chemical-induced phenotypes are enhanced with anatomy terms (tissues, fluids, and cell types) to describe the physiological system of the reported event. These same anatomy terms are used to annotate the human media (e.g., urine, hair, nail, blood, etc.) in which an environmental chemical was assayed for exposure. Currently, CTD uses more than 880 unique anatomy terms to contextualize over 255,000 chemical-phenotype interactions and 167,000 exposure statements. These annotations allow chemical-phenotype interactions and exposure data to be explored from a novel, anatomical perspective. Here, we describe CTD's anatomy curation process (including the construction of a controlled, interoperable vocabulary) and new anatomy webpages (that coalesce and organize the curated chemical-phenotype and exposure data sets). We also provide examples that demonstrate how this feature can be used to identify system- and cell-specific chemical-induced toxicities, help inform exposure data, prioritize phenotypes for environmental diseases, survey tissue and pregnancy exposomes, and facilitate data connections with external resources. Anatomy annotations advance understanding of environmental health by providing new ways to explore and survey chemical-induced events and exposure studies in the CTD framework.

STARRPeaker: uniform processing and accurate identification of STARR-seq active regions

STARR-seq technology has employed progressively more complex genomic libraries and increased sequencing depths. An issue with the increased complexity and depth is that the coverage in STARR-seq experiments is non-uniform, overdispersed, and often confounded by sequencing biases, such as GC content. Furthermore, STARR-seq readout is confounded by RNA secondary structure and thermodynamic stability. To address these potential confounders, we developed a negative binomial regression framework for uniformly processing STARR-seq data, called STARRPeaker. Moreover, to aid our effort, we generated whole-genome STARR-seq data from the HepG2 and K562 human cell lines and applied STARRPeaker to comprehensively and unbiasedly call enhancers in them.

MFSD12 mediates the import of cysteine into melanosomes and lysosomes

Dozens of genes contribute to the wide variation in human pigmentation. Many of these genes encode proteins that localize to the melanosome-the organelle, related to the lysosome, that synthesizes pigment-but have unclear functions1,2. Here we describe MelanoIP, a method for rapidly isolating melanosomes and profiling their labile metabolite contents. We use this method to study MFSD12, a transmembrane protein of unknown molecular function that, when suppressed, causes darker pigmentation in mice and humans3,4. We find that MFSD12 is required to maintain normal levels of cystine-the oxidized dimer of cysteine-in melanosomes, and to produce cysteinyldopas, the precursors of pheomelanin synthesis made in melanosomes via cysteine oxidation5,6. Tracing and biochemical analyses show that MFSD12 is necessary for the import of cysteine into melanosomes and, in non-pigmented cells, lysosomes. Indeed, loss of MFSD12 reduced the accumulation of cystine in lysosomes of fibroblasts from patients with cystinosis, a lysosomal-storage disease caused by inactivation of the lysosomal cystine exporter cystinosin7-9. Thus, MFSD12 is an essential component of the cysteine importer for melanosomes and lysosomes.

Multiple Expression Assessments of ACE2 and TMPRSS2 SARS-CoV-2 Entry Molecules in the Urinary Tract and Their Associations with Clinical Manifestations of COVID-19

Background

Since December 2019, the novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first spread quickly in Wuhan, China, then globally. Based on previously published evidence, ACE2 and TMPRSS2 are both pivotal entry molecules that enable cellular infection by SARS-CoV-2. Also, increased expression of pro-inflammatory cytokines, or a “cytokine storm,” is associated with multiple organ dysfunction syndrome often observed in critically ill patients.

Methods

We investigated the expression pattern of ACE2 and TMPRSS2 in major organs in the human body, especially in specific disease conditions. Multiple sequence alignment of ACE2 in different species was used to explain animal susceptibility. Moreover, the cell-specific expression patterns of ACE2 and cytokine receptors in the urinary tract were assessed using single-cell RNA sequencing (scRNA-seq). Additional biological relevance was determined through Gene Set Enrichment Analysis (GSEA) using an ACE2-specific signature.

Results

Our results revealed that ACE2 and TMPRSS2 were highly expressed in genitourinary organs. ACE2 was highly and significantly expressed in the kidney among individuals with chronic kidney diseases or diabetic nephropathy. In single cells, ACE2 was primarily enriched in gametocytes in the testis and renal proximal tubules. The receptors for pro-inflammatory cytokines, especially IL6ST, were notably concentrated in endothelial cells, macrophages, spermatogonial stem cells in the testis, and renal endothelial cells, which suggested the occurrence of alternative damaging autoimmune mechanisms.

Conclusion

This study provided new insights into the pathogenic mechanisms of SARS-CoV-2 that underlie the clinical manifestations observed in the human testis and kidney. These observations might substantially facilitate the development of effective treatments for this rapidly spreading disease.

A systematic evaluation of the design and context dependencies of massively parallel reporter assays

Massively parallel reporter assays (MPRAs) functionally screen thousands of sequences for regulatory activity in parallel. To date, there are limited studies that systematically compare differences in MPRA design. Here, we screen a library of 2,440 candidate liver enhancers and controls for regulatory activity in HepG2 cells using nine different MPRA designs. We identify subtle but significant differences that correlate with epigenetic and sequence-level features, as well as differences in dynamic range and reproducibility. We also validate that enhancer activity is largely independent of orientation, at least for our library and designs. Finally, we assemble and test the same enhancers as 192-mers, 354-mers and 678-mers and observe sizable differences. This work provides a framework for the experimental design of high-throughput reporter assays, suggesting that the extended sequence context of tested elements and to a lesser degree the precise assay, influence MPRA results.

Differential Characteristics of HMGB2 Versus HMGB1 and their Perspectives in Ovary and Prostate Cancer

We have summarized common and differential functions of HMGB1 and HMGB2 proteins with reference to pathological processes, with a special focus on cancer. Currently, several "omic" approaches help us compare the relative expression of these 2 proteins in healthy and cancerous human specimens, as well as in a wide range of cancer-derived cell lines, or in fetal versus adult cells. Molecules that interfere with HMGB1 functions, though through different mechanisms, have been extensively tested as therapeutic agents in animal models in recent years, and their effects are summarized. The review concludes with a discussion on the perspectives of HMGB molecules as targets in prostate and ovarian cancers.

Whole-exome sequencing in multiplex preeclampsia families identifies novel candidate susceptibility genes

OBJECTIVE

Preeclampsia is a common and serious heritable disorder of human pregnancy. Although there have been notable successes in identification of maternal susceptibility genes a large proportion of the heritability of preeclampsia remains unaccounted for. It is has been postulated that rare variation may account for some of this missing heritability. In this study, we performed whole-exome sequencing (WES) in multiplex families to identify rare exonic risk variants.

METHODS

We conducted WES in 244 individuals from 34 Australian/New Zealand multiplex preeclampsia families. Variants were tested for association with preeclampsia using a threshold model and logistic regression.

RESULTS

We found significant association for two moderately rare missense variants, rs145743393 (Padj = 0.0032, minor allele frequency = 0.016) in the chromosome 1 open reading frame 35 (C1orf35) gene, and rs34270076 (Padj = 0.0128, minor allele frequency = 0.024) in the pyroglutamylated RFamide peptide receptor (QRFPR) gene. To replicate these associations we performed imputation in our Australian genome wide association scan for preeclampsia and found no significant exonic variants in either C1orf35 or QRFPR. However, 11 variants demonstrating nominal significance (P < 0.05) in the genomic region between QRFPR and annexin A5 (ANXA5) were identified. We further leveraged publicly available genome-wide available summary data from the UK Biobank to investigate association of these two variants with the underlying clinical phenotypes of preeclampsia and detected nominal association of the QRFPR variant (rs34270076, P = 0.03) with protein levels in females.

CONCLUSION

The study represents the first to use WES in multiplex families for preeclampsia and identifies two novel genes (QRFPR and C1orf35) not previously associated with preeclampsia and find nominal association of rs34270076 with protein levels, a key clinical feature of preeclampsia. We find further support for ANXA5 previously associated with pregnancy complications, including preeclampsia.

Effect of oral urea supplementation on the endometrial transcriptome of mares

An intravenous large dose of protein led to an increased blood urea nitrogen (BUN), resulting in a lesser uterine pH and altered uterine gene expression in mares. The objective of the present study was to evaluate effects of a more physiological methodology to increase BUN on the endometrium of mares. Mares were fed hay and a treatment or control diet (n = 11 mares/treatment) in a crossover design starting at time of ovulation detection (D0) and continuing until D7. Mares of the treated group were fed urea (0.4 g/kg BW) with sweet feed and molasses, and those of the control group were fed sweet feed and molasses. Blood samples were collected daily, 1 hour after feeding, for BUN determination. Uterine and vaginal pH were determined after the last feeding on D7, and endometrial biopsies were performed. The RNA sequencing of the endometrium of a subset of mares (n = 6/treatment) was conducted. Differentially expressed genes (DEGs) between treatments were calculated (FDR-adjusted P-value<0.1). Urea-treated mares had greater BUN (P < 0.05), with no differences in uterine and vaginal pH compared to control mares. A total of 60 DEGs were characterized, those with largest fold change were SIK1, ATF3, SPINK7, NR4A1 and EGR3. Processes related to necrosis and cellular movement were predicted with the DEGs. Dietary administration of urea resulted in transcriptomic changes in the endometrium of mares related to necrosis, tissue remodeling and concentration of lipids. The observed changes in gene expression after an increased BUN might result in a disruption to the endometrium.

DASHR 2.0: integrated database of human small non-coding RNA genes and mature products

Motivation

Small non-coding RNAs (sncRNAs, <100 nts) are highly abundant RNAs that regulate diverse and often tissue-specific cellular processes by associating with transcription factor complexes or binding to mRNAs. While thousands of sncRNA genes exist in the human genome, no single resource provides searchable, unified annotation, expression and processing information for full sncRNA transcripts and mature RNA products derived from these larger RNAs.

Results

Our goal is to establish a complete catalog of annotation, expression, processing, conservation, tissue-specificity and other biological features for all human sncRNA genes and mature products derived from all major RNA classes. DASHR (Database of small human non-coding RNAs) v2.0 database is the first that integrates human sncRNA gene and mature products profiles obtained from multiple RNA-seq protocols. Altogether, 185 tissues/cell types and sncRNA annotations and >800 curated experiments from ENCODE and GEO/SRA across multiple RNA-seq protocols for both GRCh38/hg38 and GRCh37/hg19 assemblies are integrated in DASHR. Moreover, DASHR is the first to contain both known and novel, previously un-annotated sncRNA loci identified by unsupervised segmentation (13 times more loci with 1 678 800 total). Additionally, DASHR v2.0 adds >3 200 000 annotations for non-small RNA genes and other genomic features (long-noncoding RNAs, mRNAs, promoters, repeats). Furthermore, DASHR v2.0 introduces an enhanced user interface, interactive experiment-by-locus table view, sncRNA locus sorting and filtering by biological features. All annotation and expression information directly downloadable and accessible as UCSC genome browser tracks.

Availability and implementation

DASHR v2.0 is freely available at https://lisanwanglab.org/DASHRv2.

Supplementary information

Supplementary data are available at Bioinformatics online.

A genome-wide transcriptomic analysis of protein-coding genes in human blood cells

Blood is the predominant source for molecular analyses in humans, both in clinical and research settings. It is the target for many therapeutic strategies, emphasizing the need for comprehensive molecular maps of the cells constituting human blood. In this study, we performed a genome-wide transcriptomic analysis of protein-coding genes in sorted blood immune cell populations to characterize the expression levels of each individual gene across the blood cell types. All data are presented in an interactive, open-access Blood Atlas as part of the Human Protein Atlas and are integrated with expression profiles across all major tissues to provide spatial classification of all protein-coding genes. This allows for a genome-wide exploration of the expression profiles across human immune cell populations and all major human tissues and organs.

Functional annotation of the cattle genome through systematic discovery and characterization of chromatin states and butyrate-induced variations

BACKGROUND

The functional annotation of genomes, including chromatin accessibility and modifications, is important for understanding and effectively utilizing the increased amount of genome sequences reported. However, while such annotation has been well explored in a diverse set of tissues and cell types in human and model organisms, relatively little data are available for livestock genomes, hindering our understanding of complex trait variation, domestication, and adaptive evolution. Here, we present the first complete global landscape of regulatory elements in cattle and explore the dynamics of chromatin states in rumen epithelial cells induced by the rumen developmental regulator-butyrate.

RESULTS

We established the first global map of regulatory elements (15 chromatin states) and defined their coordinated activities in cattle, through genome-wide profiling for six histone modifications, RNA polymerase II, CTCF-binding sites, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial primary cells (REPC), rumen tissues, and Madin-Darby bovine kidney epithelial cells (MDBK). We demonstrated that each chromatin state exhibited specific enrichment for sequence ontology, transcription, methylation, trait-associated variants, gene expression-associated variants, selection signatures, and evolutionarily conserved elements, implying distinct biological functions. After butyrate treatments, we observed that the weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states, occurred concomitantly with significant alterations in gene expression and DNA methylation, which was significantly associated with heifer conception rate and stature economic traits.

CONCLUSION

Our results demonstrate the crucial role of functional genome annotation for understanding genome regulation, complex trait variation, and adaptive evolution in livestock. Using butyrate to induce the dynamics of the epigenomic landscape, we were able to establish the correlation among nutritional elements, chromatin states, gene activities, and phenotypic outcomes.

The landscape of transcription initiation across latent and lytic KSHV genomes

Precise promoter annotation is required for understanding the mechanistic basis of transcription initiation. In the context of complex genomes, such as herpesviruses where there is extensive genic overlap, identification of transcription start sites (TSSs) is particularly problematic and cannot be comprehensively accessed by standard RNA sequencing approaches. Kaposi's sarcoma-associated herpesvirus (KSHV) is an oncogenic gammaherpesvirus and the etiological agent of Kaposi's sarcoma and the B cell lymphoma primary effusion lymphoma (PEL). Here, we leverage RNA annotation and mapping of promoters for analysis of gene expression (RAMPAGE) and define KSHV TSSs transcriptome-wide and at nucleotide resolution in two widely used models of KSHV infection, namely iSLK.219 cells and the PEL cell line TREx-BCBL1-RTA. By mapping TSSs over a 96 h time course of reactivation we confirm 48 of 50 previously identified TSSs. Moreover, we identify over 100 novel transcription start site clusters (TSCs) in each cell line. Our analyses identified cell-type specific differences in TSC positions as well as promoter strength, and defined motifs within viral core promoters. Collectively, by defining TSSs at high resolution we have greatly expanded the transcriptional landscape of the KSHV genome and identified transcriptional control mechanisms at play during KSHV lytic reactivation.

Dynamics of Dual Specificity Phosphatases and Their Interplay with Protein Kinases in Immune Signaling

Dual specificity phosphatases (DUSPs) have a well-known role as regulators of the immune response through the modulation of mitogen-activated protein kinases (MAPKs). Yet the precise interplay between the various members of the DUSP family with protein kinases is not well understood. Recent multi-omics studies characterizing the transcriptomes and proteomes of immune cells have provided snapshots of molecular mechanisms underlying innate immune response in unprecedented detail. In this study, we focus on deciphering the interplay between members of the DUSP family with protein kinases in immune cells using publicly available omics datasets. Our analysis resulted in the identification of potential DUSP-mediated hub proteins including MAPK7, MAPK8, AURKA, and IGF1R. Furthermore, we analyzed the association of DUSP expression with TLR4 signaling and identified VEGF, FGFR, and SCF-KIT pathway modules to be regulated by the activation of TLR4 signaling. Finally, we identified several important kinases including LRRK2, MAPK8, and cyclin-dependent kinases as potential DUSP-mediated hubs in TLR4 signaling. The findings from this study have the potential to aid in the understanding of DUSP signaling in the context of innate immunity. Further, this will promote the development of therapeutic modalities for disorders with aberrant DUSP signaling.

Identification of rheumatoid arthritis causal genes using functional genomics

Over the past decade, genome-wide association studies have contributed a wealth of knowledge to our understanding of polygenic disorders such as rheumatoid arthritis. As the size of sample cohorts has improved so too have the computational and experimental methods used to robustly define variants associated with disease susceptibility. The challenge now remains to translate these findings into improved understanding of disease aetiology and patient care. Whilst much of the focus of translating the findings of genome-wide association studies has been on global analysis of all variants identified, careful functional study of individual disease susceptibility loci will be required in order to refine our understanding of how individual variants contribute to disease risk. Here, we present the argument behind such an approach and describe some of the novel tools being used to investigate risk loci. This includes the use of chromosomal conformation capture techniques and modifications of the CRISPR-Cas9 system, with several examples of their implementation being described.

Functional analysis of mammalian phospholipase D enzymes

Phosphatidylcholine (PC)-specific phospholipase D (PLD) hydrolyzes the phosphodiester bond of the PC to generate phosphatidic acid (PA) and regulates several subcellular functions. Mammalian genomes contain two genes encoding distinct isoforms of PLD in contrast with invertebrate genomes that include a single PLD gene. However, the significance of two genes within a genome encoding the same biochemical activity remains unclear. Recently, loss of function in the only PLD gene in Drosophila was reported to result in reduced PA levels and a PA-dependent collapse of the photoreceptor plasma membrane due to defects in vesicular transport. Phylogenetic analysis reveals that human PLD1 (hPLD1) is evolutionarily closer to dPLD than human PLD2 (hPLD2). In the present study, we expressed hPLD1 and hPLD2 in Drosophila and found that while reconstitution of hPLD1 is able to completely rescue retinal degeneration in a loss of function dPLD mutant, hPLD2 was less effective in its ability to mediate a rescue. Using a newly developed analytical method, we determined the acyl chain composition of PA species produced by each enzyme. While dPLD was able to restore the levels of most PA species in dPLD^3.1 cells, hPLD1 and hPLD2 each were unable to restore the levels of a subset of unique species of PA. Finally, we found that in contrast with hPLD2, dPLD and hPLD1 are uniquely distributed to the subplasma membrane region in photoreceptors. In summary, hPLD1 likely represents the ancestral PLD in mammalian genomes while hPLD2 represents neofunctionalization to generate PA at distinct subcellular membranes.

SLIC-CAGE: high-resolution transcription start site mapping using nanogram-levels of total RNA

Cap analysis of gene expression (CAGE) is a methodology for genome-wide quantitative mapping of mRNA 5′ ends to precisely capture transcription start sites at a single nucleotide resolution. In combination with high-throughput sequencing, CAGE has revolutionized our understanding of the rules of transcription initiation, led to discovery of new core promoter sequence features, and discovered transcription initiation at enhancers genome-wide. The biggest limitation of CAGE is that even the most recently improved version (nAnT-iCAGE) still requires large amounts of total cellular RNA (5 µg), preventing its application to scarce biological samples such as those from early embryonic development or rare cell types. Here, we present SLIC-CAGE, a Super-Low Input Carrier-CAGE approach to capture 5′ ends of RNA polymerase II transcripts from as little as 5–10 ng of total RNA. This dramatic increase in sensitivity is achieved by specially designed, selectively degradable carrier RNA. We demonstrate the ability of SLIC-CAGE to generate data for genome-wide promoterome with 1000-fold less material than required by existing CAGE methods, by generating a complex, high-quality library from mouse embryonic day 11.5 primordial germ cells.

High-throughput mouse phenomics for characterizing mammalian gene function

We are entering a new era of mouse phenomics, driven by large-scale and economical generation of mouse mutants coupled with increasingly sophisticated and comprehensive phenotyping. These studies are generating large, multidimensional gene-phenotype data sets, which are shedding new light on the mammalian genome landscape and revealing many hitherto unknown features of mammalian gene function. Moreover, these phenome resources provide a wealth of disease models and can be integrated with human genomics data as a powerful approach for the interpretation of human genetic variation and its relationship to disease. In the future, the development of novel phenotyping platforms allied to improved computational approaches, including machine learning, for the analysis of phenotype data will continue to enhance our ability to develop a comprehensive and powerful model of mammalian gene-phenotype space.

Whole-exome sequencing for variant discovery in blepharospasm

BACKGROUND

Blepharospasm (BSP) is a type of focal dystonia characterized by involuntary orbicularis oculi spasms that are usually bilateral, synchronous, and symmetrical. Despite strong evidence for genetic contributions to BSP, progress in the field has been constrained by small cohorts, incomplete penetrance, and late age of onset. Although several genetic etiologies for dystonia have been identified through whole-exome sequencing (WES), none of these are characteristically associated with BSP as a singular or predominant manifestation.

METHODS

We performed WES on 31 subjects from 21 independent pedigrees with BSP. The strongest candidate sequence variants derived from in silico analyses were confirmed with bidirectional Sanger sequencing and subjected to cosegregation analysis.

RESULTS

Cosegregating deleterious variants (GRCH37/hg19) in CACNA1A (NM_001127222.1: c.7261_7262delinsGT, p.Pro2421Val), REEP4 (NM_025232.3: c.109C>T, p.Arg37Trp), TOR2A (NM_130459.3: c.568C>T, p.Arg190Cys), and ATP2A3 (NM_005173.3: c.1966C>T, p.Arg656Cys) were identified in four independent multigenerational pedigrees. Deleterious variants in HS1BP3 (NM_022460.3: c.94C>A, p.Gly32Cys) and GNA14 (NM_004297.3: c.989_990del, p.Thr330ArgfsTer67) were identified in a father and son with segmental cranio-cervical dystonia first manifest as BSP. Deleterious variants in DNAH17, TRPV4, CAPN11, VPS13C, UNC13B, SPTBN4, MYOD1, and MRPL15 were found in two or more independent pedigrees. To our knowledge, none of these genes have previously been associated with isolated BSP, although other CACNA1A mutations have been associated with both positive and negative motor disorders including ataxia, episodic ataxia, hemiplegic migraine, and dystonia.

CONCLUSIONS

Our WES datasets provide a platform for future studies of BSP genetics which will demand careful consideration of incomplete penetrance, pleiotropy, population stratification, and oligogenic inheritance patterns.

Non-coding RNAs: long non-coding RNAs and microRNAs in endocrine-related cancers

The human genome is 'pervasively transcribed' leading to a complex array of non-coding RNAs (ncRNAs) that far outnumber coding mRNAs. ncRNAs have regulatory roles in transcription and post-transcriptional processes as well numerous cellular functions that remain to be fully described. Best characterized of the 'expanding universe' of ncRNAs are the ~22 nucleotide microRNAs (miRNAs) that base-pair to target mRNA's 3' untranslated region within the RNA-induced silencing complex (RISC) and block translation and may stimulate mRNA transcript degradation. Long non-coding RNAs (lncRNAs) are classified as >200 nucleotides in length, but range up to several kb and are heterogeneous in genomic origin and function. lncRNAs fold into structures that interact with DNA, RNA and proteins to regulate chromatin dynamics, protein complex assembly, transcription, telomere biology and splicing. Some lncRNAs act as sponges for miRNAs and decoys for proteins. Nuclear-encoded lncRNAs can be taken up by mitochondria and lncRNAs are transcribed from mtDNA. Both miRNAs and lncRNAs are dysregulated in endocrine cancers. This review provides an overview on the current understanding of the regulation and function of selected lncRNAs and miRNAs, and their interaction, in endocrine-related cancers: breast, prostate, endometrial and thyroid.