CGAT: computational genomics analysis toolkit

The rate of W chromosome degeneration across multiple avian neo-sex chromosomes

When sex chromosomes evolve recombination suppression, the sex-limited chromosome (Y/W) commonly degenerate by losing functional genes. The rate of Y/W degeneration is believed to slow down over time as the most essential genes are maintained by purifying selection, but supporting data are scarce especially for ZW systems. Here, we study W degeneration in Sylvioidea songbirds where multiple autosomal translocations to the sex chromosomes, and multiple recombination suppression events causing separate evolutionary strata, have occurred during the last ~ 28.1-4.5 million years (Myr). We show that the translocated regions have maintained 68.3-97.7% of their original gene content, compared to only 4.2% on the much older ancestral W chromosome. By mapping W gene losses onto a dated phylogeny, we estimate an average gene loss rate of 1.0% per Myr, with only moderate variation between four independent lineages. Consistent with previous studies, evolutionarily constrained and haploinsufficient genes were preferentially maintained on W. However, the gene loss rate did not show any consistent association with strata age or with the number of W genes at strata formation. Our study provides a unique account on the pace of W gene loss and reinforces the significance of purifying selection in maintaining essential genes on sex chromosomes.

In silico genome wide identification of long non-coding RNAs differentially expressed during Candida auris host pathogenesis

Candida auris is an invasive fungal pathogen of high concern due to acquired drug tolerance against antifungals used in clinics. The prolonged persistence on biotic and abiotic surfaces can result in onset of hospital outbreaks causing serious health threat. An in depth understanding of pathology of C. auris is highly desirable for development of efficient therapeutics. Non-coding RNAs play crucial role in fungal pathology. However, the information about ncRNAs is scanty to be utilized. Herein our aim is to identify long noncoding RNAs with potent role in pathobiology of C. auris. Thereby, we analyzed the transcriptomics data of C. auris infection in blood for identification of potential lncRNAs with regulatory role in determining invasion, survival or drug tolerance under infection conditions. Interestingly, we found 275 lncRNAs, out of which 253 matched with lncRNAs reported in Candidamine, corroborating for our accurate data analysis pipeline. Nevertheless, we obtained 23 novel lncRNAs not reported earlier. Three lncRNAs were found to be under expressed throughout the course of infection, in the transcriptomics data. 16 of potent lncRNAs were found to be coexpressed with coding genes, emphasizing for their functional role. Noteworthy, these ncRNAs are expressed from intergenic regions of the genes associated with transporters, metabolism, cell wall biogenesis. This study recommends for possible association between lncRNA expression and C. auris pathogenesis.

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

Benchmarking single-cell RNA-seq (scRNA-seq) and single-cell Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) computational tools demands simulators to generate realistic sequencing reads. However, none of the few read simulators aim to mimic real data. To fill this gap, we introduce scReadSim, a single-cell RNA-seq and ATAC-seq read simulator that allows user-specified ground truths and generates synthetic sequencing reads (in a FASTQ or BAM file) by mimicking real data. At both read-sequence and read-count levels, scReadSim mimics real scRNA-seq and scATAC-seq data. Moreover, scReadSim provides ground truths, including unique molecular identifier (UMI) counts for scRNA-seq and open chromatin regions for scATAC-seq. In particular, scReadSim allows users to design cell-type-specific ground-truth open chromatin regions for scATAC-seq data generation. In benchmark applications of scReadSim, we show that UMI-tools achieves the top accuracy in scRNA-seq UMI deduplication, and HMMRATAC and MACS3 achieve the top performance in scATAC-seq peak calling.

MYC is a clinically significant driver of mTOR inhibitor resistance in breast cancer

Targeting the PI3K-AKT-mTOR pathway is a promising therapeutic strategy for breast cancer treatment. However, low response rates and development of resistance to PI3K-AKT-mTOR inhibitors remain major clinical challenges. Here, we show that MYC activation drives resistance to mTOR inhibitors (mTORi) in breast cancer. Multiomic profiling of mouse invasive lobular carcinoma (ILC) tumors revealed recurrent Myc amplifications in tumors that acquired resistance to the mTORi AZD8055. MYC activation was associated with biological processes linked to mTORi response and counteracted mTORi-induced translation inhibition by promoting translation of ribosomal proteins. In vitro and in vivo induction of MYC conferred mTORi resistance in mouse and human breast cancer models. Conversely, AZD8055-resistant ILC cells depended on MYC, as demonstrated by the synergistic effects of mTORi and MYCi combination treatment. Notably, MYC status was significantly associated with poor response to everolimus therapy in metastatic breast cancer patients. Thus, MYC is a clinically relevant driver of mTORi resistance that may stratify breast cancer patients for mTOR-targeted therapies.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

Pervasive downstream RNA hairpins dynamically dictate start-codon selection

Translational reprogramming allows organisms to adapt to changing conditions. Upstream start codons (uAUGs), which are prevalently present in mRNAs, have crucial roles in regulating translation by providing alternative translation start sites1-4. However, what determines this selective initiation of translation between conditions remains unclear. Here, by integrating transcriptome-wide translational and structural analyses during pattern-triggered immunity in Arabidopsis, we found that transcripts with immune-induced translation are enriched with upstream open reading frames (uORFs). Without infection, these uORFs are selectively translated owing to hairpins immediately downstream of uAUGs, presumably by slowing and engaging the scanning preinitiation complex. Modelling using deep learning provides unbiased support for these recognizable double-stranded RNA structures downstream of uAUGs (which we term uAUG-ds) being responsible for the selective translation of uAUGs, and allows the prediction and rational design of translating uAUG-ds. We found that uAUG-ds-mediated regulation can be generalized to human cells. Moreover, uAUG-ds-mediated start-codon selection is dynamically regulated. After immune challenge in plants, induced RNA helicases that are homologous to Ded1p in yeast and DDX3X in humans resolve these structures, allowing ribosomes to bypass uAUGs to translate downstream defence proteins. This study shows that mRNA structures dynamically regulate start-codon selection. The prevalence of this RNA structural feature and the conservation of RNA helicases across kingdoms suggest that mRNA structural remodelling is a general feature of translational reprogramming.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Tissue-dependent transcriptional and bacterial associations in primary sclerosing cholangitis-associated inflammatory bowel disease

Background: Patients with primary sclerosing cholangitis (PSC) frequently have co-ocurring ulcerative colitis (UC) and develop colorectal cancer. Colorectal cancer risk in patients with PSC-associated ulcerative colitis (PSC/UC) is elevated relative to patients with ulcerative colitis (UC) alone, reasons for which remain obscure. Understanding the molecular and microbial basis for differences between these two patient groups and how these vary across intestinal sites is important for the development of therapies to prevent colorectal cancer development in at-risk individuals.   Methods: We employed ribonucleic acid sequencing (RNA-seq) analysis of biopsy samples across three intestinal tissue locations (ileum, caecum and rectum) in patients with PSC/UC (ileum n = 7, caecum n = 7, rectum n = 7), UC (ileum n = 9, caecum n = 10, rectum n = 10) and healthy controls (ileum n = 11, caecum n = 9, rectum n = 12) to determine tissue-dependent transcriptional alterations in PSC/UC. We also performed 16S ribosomal RNA (rRNA) amplicon sequencing to determine bacterial associations with PSC/UC. Results: Tissue-defining transcriptional signatures revealed that the ileum was enriched for genes involved in lipid and drug metabolism, the caecum for activated immune cells and the rectum for enteric neurogenesis. Transcriptional alterations relative to healthy control samples were largely shared between patients with PSC/UC or UC although were distinct across tissue locations. Nevertheless, we observed reduced expression of gamma-glutamyl transferase 1 ( GGT1) specifically in the ileum and caecum of patients with PSC/UC. Analysis of the bacterial component of the microbiome revealed high inter-individual variability of microbiome composition and little evidence for tissue-dependency. We observed a reduction in Parabacteroides relative abundance in the rectum of patients with PSC/UC. Conclusions: The role of gamma-glutamyl transferase in maintaining the redox environment through the glutathione salvage pathway makes our observed alterations a potential pathway to PSC-associated colorectal cancer.

A mechanism for oxidative damage repair at gene regulatory elements

Oxidative genome damage is an unavoidable consequence of cellular metabolism. It arises at gene regulatory elements by epigenetic demethylation during transcriptional activation^1,2. Here we show that promoters are protected from oxidative damage via a process mediated by the nuclear mitotic apparatus protein NuMA (also known as NUMA1). NuMA exhibits genomic occupancy approximately 100 bp around transcription start sites. It binds the initiating form of RNA polymerase II, pause-release factors and single-strand break repair (SSBR) components such as TDP1. The binding is increased on chromatin following oxidative damage, and TDP1 enrichment at damaged chromatin is facilitated by NuMA. Depletion of NuMA increases oxidative damage at promoters. NuMA promotes transcription by limiting the polyADP-ribosylation of RNA polymerase II, increasing its availability and release from pausing at promoters. Metabolic labelling of nascent RNA identifies genes that depend on NuMA for transcription including immediate-early response genes. Complementation of NuMA-deficient cells with a mutant that mediates binding to SSBR, or a mitotic separation-of-function mutant, restores SSBR defects. These findings underscore the importance of oxidative DNA damage repair at gene regulatory elements and describe a process that fulfils this function.

Assessing and assuring interoperability of a genomics file format

Motivation

Bioinformatics software tools operate largely through the use of specialized genomics file formats. Often these formats lack formal specification, making it difficult or impossible for the creators of these tools to robustly test them for correct handling of input and output. This causes problems in interoperability between different tools that, at best, wastes time and frustrates users. At worst, interoperability issues could lead to undetected errors in scientific results.

Results

We developed a new verification system, Acidbio, which tests for correct behavior in bioinformatics software packages. We crafted tests to unify correct behavior when tools encounter various edge cases—potentially unexpected inputs that exemplify the limits of the format. To analyze the performance of existing software, we tested the input validation of 80 Bioconda packages that parsed the Browser Extensible Data (BED) format. We also used a fuzzing approach to automatically perform additional testing. Of 80 software packages examined, 75 achieved less than 70% correctness on our test suite. We categorized multiple root causes for the poor performance of different types of software. Fuzzing detected other errors that the manually designed test suite could not. We also created a badge system that developers can use to indicate more precisely which BED variants their software accepts and to advertise the software’s performance on the test suite.

Availability and implementation

Acidbio is available at https://github.com/hoffmangroup/acidbio.

Supplementary information

Supplementary data are available at Bioinformatics online.

Chromosomal level genome of Ilex asprella and insight into antiviral triterpenoid pathway

Ilex asprella is a widely used herbs in Traditional Chinese Medicine for treating viral infection and relieving inflammation. Due to the earlier fruiting period of I. asprella, it is the major food source for frugivores in summer. Despite its pharmacological and ecological importance, a reference genome for I. asprella is lacking. By using Illumina, stLFR and Omni-C sequencing data, we present the first chromosomal-level assembly for I. asprella. The genome assembly size is 804 Mbp, with Benchmarking Universal Single-Copy Orthologs (BUSCO) score 94.4% for eudicotyledon single copy genes. Transcriptomes of leaves, stems, flowers, premature fruits and roots were analyzed, providing 39,215 gene models. The complete set of genes involved in the triterpenoids production is disclosed for the first time. We have also found the oxidosqualene cyclases (OSCs), CYP716s and UDP-glycosyltransferases (UGTs), which are responsible for the modification of triterpenoid backbones, resulting in the high variety of triterpenoid saponins.

The long non-coding RNA landscape of Candida yeast pathogens

Long non-coding RNAs (lncRNAs) constitute a poorly studied class of transcripts with emerging roles in key cellular processes. Despite efforts to characterize lncRNAs across a wide range of species, these molecules remain largely unexplored in most eukaryotic microbes, including yeast pathogens of the Candida clade. Here, we analyze thousands of publicly available sequencing datasets to infer and characterize the lncRNA repertoires of five major Candida pathogens: Candida albicans, Candida tropicalis, Candida parapsilosis, Candida auris and Candida glabrata. Our results indicate that genomes of these species encode hundreds of lncRNAs that show levels of evolutionary constraint intermediate between those of intergenic genomic regions and protein-coding genes. Despite their low sequence conservation across the studied species, some lncRNAs are syntenic and are enriched in shared sequence motifs. We find co-expression of lncRNAs with certain protein-coding transcripts, hinting at potential functional associations. Finally, we identify lncRNAs that are differentially expressed during infection of human epithelial cells for four of the studied species. Our comprehensive bioinformatic analyses of Candida lncRNAs pave the way for future functional characterization of these transcripts.

Long non-coding RNAs (lncRNAs) play roles in key cellular processes, but remain largely unexplored in fungal pathogens such as Candida. Here, Hovhannisyan and Gabaldón analyze thousands of sequencing datasets to infer and characterize the lncRNA repertoires of five Candida species, paving the way for their future functional characterization.

Loss of mutual protection between human osteoclasts and chondrocytes in damaged joints initiates osteoclast-mediated cartilage degradation by MMPs

Osteoclasts are multinucleated, bone-resorbing cells. However, they also digest cartilage during skeletal maintenance, development and in degradative conditions including osteoarthritis, rheumatoid arthritis and primary bone sarcoma. This study explores the mechanisms behind the osteoclast–cartilage interaction. Human osteoclasts differentiated on acellular human cartilage expressed osteoclast marker genes (e.g. CTSK, MMP9) and proteins (TRAP, VNR), visibly damaged the cartilage surface and released glycosaminoglycan in a contact-dependent manner. Direct co-culture with chondrocytes during differentiation increased large osteoclast formation (p < 0.0001) except when co-cultured on dentine, when osteoclast formation was inhibited (p = 0.0002). Osteoclasts cultured on dentine inhibited basal cartilage degradation (p = 0.012). RNA-seq identified MMP8 overexpression in osteoclasts differentiated on cartilage versus dentine (8.89-fold, p = 0.0133), while MMP9 was the most highly expressed MMP. Both MMP8 and MMP9 were produced by osteoclasts in osteosarcoma tissue. This study suggests that bone-resident osteoclasts and chondrocytes exert mutually protective effects on their ‘native’ tissue. However, when osteoclasts contact non-native cartilage they cause degradation via MMPs. Understanding the role of osteoclasts in cartilage maintenance and degradation might identify new therapeutic approaches for pathologies characterized by cartilage degeneration.

Productive visualization of high-throughput sequencing data using the SeqCode open portable platform

Large-scale sequencing techniques to chart genomes are entirely consolidated. Stable computational methods to perform primary tasks such as quality control, read mapping, peak calling, and counting are likewise available. However, there is a lack of uniform standards for graphical data mining, which is also of central importance. To fill this gap, we developed SeqCode, an open suite of applications that analyzes sequencing data in an elegant but efficient manner. Our software is a portable resource written in ANSI C that can be expected to work for almost all genomes in any computational configuration. Furthermore, we offer a user-friendly front-end web server that integrates SeqCode functions with other graphical analysis tools. Our analysis and visualization toolkit represents a significant improvement in terms of performance and usability as compare to other existing programs. Thus, SeqCode has the potential to become a key multipurpose instrument for high-throughput professional analysis; further, it provides an extremely useful open educational platform for the world-wide scientific community. SeqCode website is hosted at http://ldicrocelab.crg.eu, and the source code is freely distributed at https://github.com/eblancoga/seqcode.

Tissue-dependent transcriptional and bacterial associations in primary sclerosing cholangitis-associated inflammatory bowel disease

Background: Primary sclerosing cholangitis (PSC) is a disease of the bile duct and liver. However, patients frequently have co-morbidities including inflammatory bowel disease (IBD) and colorectal cancer. Colorectal cancer risk in patients with PSC-associated ulcerative colitis (PSC/UC) is elevated relative to patients with ulcerative colitis (UC) alone, reasons for which remain obscure. Further, clinical and immunological features, and involved intestinal sites differ between PSC/UC and UC. Understanding the molecular and microbial basis for differences in cancer risk between these two patient groups and how these differ across intestinal sites is important for the development of therapies to prevent colorectal cancer development in at-risk individuals.   Methods: We employed ribonucleic acid sequencing (RNA-seq) analysis of biopsy samples across three intestinal tissue locations (ileum, caecum and rectum) in patients with PSC/UC (n = 8), UC (n = 10) and healthy controls (n = 12) to determine tissue-dependent transcriptional alterations in PSC/UC. We also performed 16S ribosomal RNA (rRNA) amplicon sequencing to determine bacterial associations with PSC/UC and host-microbiome associations. Results: Tissue-defining transcriptional signatures revealed that the ileum was enriched for genes involved in lipid and drug metabolism, the caecum for activated immune cells and the rectum for enteric neurogenesis. Transcriptional alterations relative to healthy control samples were largely shared between patients with PSC/UC or UC although were distinct across tissue locations. Nevertheless, we observed reduced expression of gamma-glutamyl transferase 1 ( GGT1) specifically in the ileum and caecum of patients with PSC/UC. Analysis of the bacterial component of the microbiome revealed high inter-individual variability of microbiome composition and little evidence for tissue-dependency. We observed a reduction in Parabacteroides relative abundance in the rectum of patients with PSC/UC. Conclusions: The role of gamma-glutamyl transferase in maintaining the redox environment through the glutathione salvage pathway makes our observed alterations a potential pathway to PSC-associated colorectal cancer.

m6 A modification of HSATIII lncRNAs regulates temperature-dependent splicing

Nuclear stress bodies (nSBs) are nuclear membraneless organelles formed around stress-inducible HSATIII architectural long noncoding RNAs (lncRNAs). nSBs repress splicing of hundreds of introns during thermal stress recovery, which are partly regulated by CLK1 kinase phosphorylation of temperature-dependent Ser/Arg-rich splicing factors (SRSFs). Here, we report a distinct mechanism for this splicing repression through protein sequestration by nSBs. Comprehensive identification of RNA-binding proteins revealed HSATIII association with proteins related to N⁶ -methyladenosine (m⁶ A) RNA modification. 11% of the first adenosine in the repetitive HSATIII sequence were m⁶ A-modified. nSBs sequester the m⁶ A writer complex to methylate HSATIII, leading to subsequent sequestration of the nuclear m⁶ A reader, YTHDC1. Sequestration of these factors from the nucleoplasm represses m⁶ A modification of pre-mRNAs, leading to repression of m⁶ A-dependent splicing during stress recovery phase. Thus, nSBs serve as a common platform for regulation of temperature-dependent splicing through dual mechanisms employing two distinct ribonucleoprotein modules with partially m⁶ A-modified architectural lncRNAs.

The Prolyl-tRNA Synthetase Inhibitor Halofuginone Inhibits SARS-CoV-2 Infection

We identify the prolyl-tRNA synthetase (PRS) inhibitor halofuginone ¹ , a compound in clinical trials for anti-fibrotic and anti-inflammatory applications ² , as a potent inhibitor of SARS-CoV-2 infection and replication. The interaction of SARS-CoV-2 spike protein with cell surface heparan sulfate (HS) promotes viral entry ³ . We find that halofuginone reduces HS biosynthesis, thereby reducing spike protein binding, SARS-CoV-2 pseudotyped virus, and authentic SARS-CoV-2 infection. Halofuginone also potently suppresses SARS-CoV-2 replication post-entry and is 1,000-fold more potent than Remdesivir ⁴ . Inhibition of HS biosynthesis and SARS-CoV-2 infection depends on specific inhibition of PRS, possibly due to translational suppression of proline-rich proteins. We find that pp1a and pp1ab polyproteins of SARS-CoV-2, as well as several HS proteoglycans, are proline-rich, which may make them particularly vulnerable to halofuginone's translational suppression. Halofuginone is orally bioavailable, has been evaluated in a phase I clinical trial in humans and distributes to SARS-CoV-2 target organs, including the lung, making it a near-term clinical trial candidate for the treatment of COVID-19.

Dissecting the Role of BET Bromodomain Proteins BRD2 and BRD4 in Human NK Cell Function

Natural killer (NK) cells are innate lymphocytes that play a pivotal role in the immune surveillance and elimination of transformed or virally infected cells. Using a chemo-genetic approach, we identify BET bromodomain containing proteins BRD2 and BRD4 as central regulators of NK cell functions, including direct cytokine secretion, NK cell contact-dependent inflammatory cytokine secretion from monocytes as well as NK cell cytolytic functions. We show that both BRD2 and BRD4 control inflammatory cytokine production in NK cells isolated from healthy volunteers and from rheumatoid arthritis patients. In contrast, knockdown of BRD4 but not of BRD2 impairs NK cell cytolytic responses, suggesting BRD4 as critical regulator of NK cell mediated tumor cell elimination. This is supported by pharmacological targeting where the first-generation pan-BET bromodomain inhibitor JQ1(+) displays anti-inflammatory effects and inhibit tumor cell eradication, while the novel bivalent BET bromodomain inhibitor AZD5153, which shows differential activity towards BET family members, does not. Given the important role of both cytokine-mediated inflammatory microenvironment and cytolytic NK cell activities in immune-oncology therapies, our findings present a compelling argument for further clinical investigation.

Gene expression in diapausing rotifer eggs in response to divergent environmental predictability regimes

In unpredictable environments in which reliable cues for predicting environmental variation are lacking, a diversifying bet-hedging strategy for diapause exit is expected to evolve, whereby only a portion of diapausing forms will resume development at the first occurrence of suitable conditions. This study focused on diapause termination in the rotifer Brachionus plicatilis s.s., addressing the transcriptional profile of diapausing eggs from environments differing in the level of predictability and the relationship of such profiles with hatching patterns. RNA-Seq analyses revealed significant differences in gene expression between diapausing eggs produced in the laboratory under combinations of two contrasting selective regimes of environmental fluctuation (predictable vs unpredictable) and two different diapause conditions (passing or not passing through forced diapause). The results showed that the selective regime was more important than the diapause condition in driving differences in the transcriptome profile. Most of the differentially expressed genes were upregulated in the predictable regime and mostly associated with molecular functions involved in embryo morphological development and hatching readiness. This was in concordance with observations of earlier, higher, and more synchronous hatching in diapausing eggs produced under the predictable regime.

Mountain stoneflies may tolerate warming streams: Evidence from organismal physiology and gene expression

Rapid glacier recession is altering the physical conditions of headwater streams. Stream temperatures are predicted to rise and become increasingly variable, putting entire meltwater-associated biological communities at risk of extinction. Thus, there is a pressing need to understand how thermal stress affects mountain stream insects, particularly where glaciers are likely to vanish on contemporary timescales. In this study, we measured the critical thermal maximum (CT_MAX ) of stonefly nymphs representing multiple species and a range of thermal regimes in the high Rocky Mountains, USA. We then collected RNA-sequencing data to assess how organismal thermal stress translated to the cellular level. Our focal species included the meltwater stonefly, Lednia tumana, which was recently listed under the U.S. Endangered Species Act due to climate-induced habitat loss. For all study species, critical thermal maxima (CT_MAX  > 20°C) far exceeded the stream temperatures mountain stoneflies experience (<10°C). Moreover, while evidence for a cellular stress response was present, we also observed constitutive expression of genes encoding proteins known to underlie thermal stress (i.e., heat shock proteins) even at low temperatures that reflected natural conditions. We show that high-elevation aquatic insects may not be physiologically threatened by short-term exposure to warm temperatures and that longer-term physiological responses or biotic factors (e.g., competition) may better explain their extreme distributions.

Correction of amyotrophic lateral sclerosis related phenotypes in induced pluripotent stem cell-derived motor neurons carrying a hexanucleotide expansion mutation in C9orf72 by CRISPR/Cas9 genome editing using homology-directed repair

The G4C2 hexanucleotide repeat expansion (HRE) in C9orf72 is the commonest cause of familial amyotrophic lateral sclerosis (ALS). A number of different methods have been used to generate isogenic control lines using clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 and non-homologous end-joining by deleting the repeat region, with the risk of creating indels and genomic instability. In this study, we demonstrate complete correction of an induced pluripotent stem cell (iPSC) line derived from a C9orf72-HRE positive ALS/frontotemporal dementia patient using CRISPR/Cas9 genome editing and homology-directed repair (HDR), resulting in replacement of the excised region with a donor template carrying the wild-type repeat size to maintain the genetic architecture of the locus. The isogenic correction of the C9orf72 HRE restored normal gene expression and methylation at the C9orf72 locus, reduced intron retention in the edited lines and abolished pathological phenotypes associated with the C9orf72 HRE expansion in iPSC-derived motor neurons (iPSMNs). RNA sequencing of the mutant line identified 2220 differentially expressed genes compared with its isogenic control. Enrichment analysis demonstrated an over-representation of ALS relevant pathways, including calcium ion dependent exocytosis, synaptic transport and the Kyoto Encyclopedia of Genes and Genomes ALS pathway, as well as new targets of potential relevance to ALS pathophysiology. Complete correction of the C9orf72 HRE in iPSMNs by CRISPR/Cas9-mediated HDR provides an ideal model to study the earliest effects of the hexanucleotide expansion on cellular homeostasis and the key pathways implicated in ALS pathophysiology.

Genome wide distribution of G-quadruplexes and their impact on gene expression in malaria parasites

Mechanisms of transcriptional control in malaria parasites are still not fully understood. The positioning patterns of G-quadruplex (G4) DNA motifs in the parasite's AT-rich genome, especially within the var gene family which encodes virulence factors, and in the vicinity of recombination hotspots, points towards a possible regulatory role of G4 in gene expression and genome stability. Here, we carried out the most comprehensive genome-wide survey, to date, of G4s in the Plasmodium falciparum genome using G4Hunter, which identifies G4 forming sequences (G4FS) considering their G-richness and G-skewness. We show an enrichment of G4FS in nucleosome-depleted regions and in the first exon of var genes, a pattern that is conserved within the closely related Laverania Plasmodium parasites. Under G4-stabilizing conditions, i.e., following treatment with pyridostatin (a high affinity G4 ligand), we show that a bona fide G4 found in the non-coding strand of var promoters modulates reporter gene expression. Furthermore, transcriptional profiling of pyridostatin-treated parasites, shows large scale perturbations, with deregulation affecting for instance the ApiAP2 family of transcription factors and genes involved in ribosome biogenesis. Overall, our study highlights G4s as important DNA secondary structures with a role in Plasmodium gene expression regulation, sub-telomeric recombination and var gene biology.

Histone H3K27me3 demethylases regulate human Th17 cell development and effector functions by impacting on metabolism

T cells control many immune functions, with Th17 cells critical in regulating inflammation. Following activation, T cells undergo metabolic reprogramming and utilize glycolysis to increase the ATP availability. Epigenetic mechanisms controlling metabolic functions in T cells are currently not well-defined. Here, we establish an epigenetic link between the histone H3K27me3 demethylases KDM6A/B and the coordination of a metabolic response. Inhibition of KDM6A/B leads to global increases in the repressive H3K27me3 histone mark, resulting in down-regulation of key transcription factors, followed by metabolic reprogramming and anergy. This work suggests a critical role of H3K27 demethylase enzymes in maintaining Th17 functions by controlling metabolic switches. Short-term treatment with KDM6 enzyme inhibitors may be useful in the therapy of chronic inflammatory diseases.

T helper (Th) cells are CD4⁺ effector T cells that play a critical role in immunity by shaping the inflammatory cytokine environment in a variety of physiological and pathological situations. Using a combined chemico-genetic approach, we identify histone H3K27 demethylases KDM6A and KDM6B as central regulators of human Th subsets. The prototypic KDM6 inhibitor GSK-J4 increases genome-wide levels of the repressive H3K27me3 chromatin mark and leads to suppression of the key transcription factor RORγt during Th17 differentiation. In mature Th17 cells, GSK-J4 induces an altered transcriptional program with a profound metabolic reprogramming and concomitant suppression of IL-17 cytokine levels and reduced proliferation. Single-cell analysis reveals a specific shift from highly inflammatory cell subsets toward a resting state upon demethylase inhibition. The root cause of the observed antiinflammatory phenotype in stimulated Th17 cells is reduced expression of key metabolic transcription factors, such as PPRC1. Overall, this leads to reduced mitochondrial biogenesis, resulting in a metabolic switch with concomitant antiinflammatory effects. These data are consistent with an effect of GSK-J4 on Th17 T cell differentiation pathways directly related to proliferation and include regulation of effector cytokine profiles. This suggests that inhibiting KDM6 demethylases may be an effective, even in the short term, therapeutic target for autoimmune diseases, including ankylosing spondylitis.

LncRNA-dependent nuclear stress bodies promote intron retention through SR protein phosphorylation

A number of long noncoding RNAs (lncRNAs) are induced in response to specific stresses to construct membrane-less nuclear bodies; however, their function remains poorly understood. Here, we report the role of nuclear stress bodies (nSBs) formed on highly repetitive satellite III (HSATIII) lncRNAs derived from primate-specific satellite III repeats upon thermal stress exposure. A transcriptomic analysis revealed that depletion of HSATIII lncRNAs, resulting in elimination of nSBs, promoted splicing of 533 retained introns during thermal stress recovery. A HSATIII-Comprehensive identification of RNA-binding proteins by mass spectrometry (ChIRP-MS) analysis identified multiple splicing factors in nSBs, including serine and arginine-rich pre-mRNA splicing factors (SRSFs), the phosphorylation states of which affect splicing patterns. SRSFs are rapidly de-phosphorylated upon thermal stress exposure. During stress recovery, CDC like kinase 1 (CLK1) was recruited to nSBs and accelerated the re-phosphorylation of SRSF9, thereby promoting target intron retention. Our findings suggest that HSATIII-dependent nSBs serve as a conditional platform for phosphorylation of SRSFs by CLK1 to promote the rapid adaptation of gene expression through intron retention following thermal stress exposure.

The MITF-SOX10 regulated long non-coding RNA DIRC3 is a melanoma tumour suppressor

The MITF and SOX10 transcription factors regulate the expression of genes important for melanoma proliferation, invasion and metastasis. Despite growing evidence of the contribution of long noncoding RNAs (lncRNAs) in cancer, including melanoma, their functions within MITF-SOX10 transcriptional programmes remain poorly investigated. Here we identify 245 candidate melanoma associated lncRNAs whose loci are co-occupied by MITF-SOX10 and that are enriched at active enhancer-like regions. Our work suggests that one of these, Disrupted In Renal Carcinoma 3 (DIRC3), may be a clinically important MITF-SOX10 regulated tumour suppressor. DIRC3 depletion in human melanoma cells leads to increased anchorage-independent growth, a hallmark of malignant transformation, whilst melanoma patients classified by low DIRC3 expression have decreased survival. DIRC3 is a nuclear lncRNA that activates expression of its neighbouring IGFBP5 tumour suppressor through modulating chromatin structure and suppressing SOX10 binding to putative regulatory elements within the DIRC3 locus. In turn, DIRC3 dependent regulation of IGFBP5 impacts the expression of genes involved in cancer associated processes and is needed for DIRC3 control of anchorage-independent growth. Our work indicates that lncRNA components of MITF-SOX10 networks are an important new class of melanoma regulators and candidate therapeutic targets that can act not only as downstream mediators of MITF-SOX10 function but as feedback regulators of MITF-SOX10 activity.

Transcriptome-wide dynamics of extensive m6A mRNA methylation during Plasmodium falciparum blood-stage development

Malaria pathogenesis results from the asexual replication of Plasmodium falciparum within human red blood cells, which relies on a precisely timed cascade of gene expression over a 48-h life cycle. Although substantial post-transcriptional regulation of this hardwired program has been observed, it remains unclear how these processes are mediated on a transcriptome-wide level. To this end, we identified mRNA modifications in the P. falciparum transcriptome and performed a comprehensive characterization of N⁶-methyladenosine (m⁶A) over the course of blood-stage development. Using mass spectrometry and m⁶A RNA sequencing, we demonstrate that m⁶A is highly developmentally regulated, exceeding m⁶A levels known in any other eukaryote. We characterize a distinct m⁶A writer complex and show that knockdown of the putative m⁶A methyltransferase, PfMT-A70, by CRISPR interference leads to increased levels of transcripts that normally contain m⁶A. In accordance, we find an inverse correlation between m⁶A methylation and mRNA stability or translational efficiency. We further identify two putative m⁶A-binding YTH proteins that are likely to be involved in the regulation of these processes across the parasite's life cycle. Our data demonstrate unique features of an extensive m⁶A mRNA methylation programme in malaria parasites and reveal its crucial role in dynamically fine-tuning the transcriptional cascade of a unicellular eukaryote.

Bromodomain inhibition of the coactivators CBP/EP300 facilitate cellular reprogramming

Silencing of the somatic cell type-specific genes is a critical yet poorly understood step in reprogramming. To uncover pathways that maintain cell identity, we performed a reprogramming screen using inhibitors of chromatin factors. Here, we identify acetyl-lysine competitive inhibitors targeting the bromodomains of coactivators CREB (cyclic-AMP response element binding protein) binding protein (CBP) and E1A binding protein of 300 kDa (EP300) as potent enhancers of reprogramming. These inhibitors accelerate reprogramming, are critical during its early stages and, when combined with DOT1L inhibition, enable efficient derivation of human induced pluripotent stem cells (iPSCs) with OCT4 and SOX2. In contrast, catalytic inhibition of CBP/EP300 prevents iPSC formation, suggesting distinct functions for different coactivator domains in reprogramming. CBP/EP300 bromodomain inhibition decreases somatic-specific gene expression, histone H3 lysine 27 acetylation (H3K27Ac) and chromatin accessibility at target promoters and enhancers. The master mesenchymal transcription factor PRRX1 is one such functionally important target of CBP/EP300 bromodomain inhibition. Collectively, these results show that CBP/EP300 bromodomains sustain cell-type-specific gene expression and maintain cell identity.

An essential role for the Zn2+ transporter ZIP7 in B cell development

Despite the known importance of zinc for human immunity, molecular insights into its roles have remained limited. Here we report a novel autosomal recessive disease characterized by absent B cells, agammaglobulinemia and early onset infections in five unrelated families. The immunodeficiency results from hypomorphic mutations of SLC39A7, which encodes the endoplasmic reticulum-to-cytoplasm zinc transporter ZIP7. Using CRISPR-Cas9 mutagenesis we have precisely modeled ZIP7 deficiency in mice. Homozygosity for a null allele caused embryonic death, but hypomorphic alleles reproduced the block in B cell development seen in patients. B cells from mutant mice exhibited a diminished concentration of cytoplasmic free zinc, increased phosphatase activity and decreased phosphorylation of signaling molecules downstream of the pre-B cell and B cell receptors. Our findings highlight a specific role for cytosolic Zn2+ in modulating B cell receptor signal strength and positive selection.

Comprehensive identification of RNA-protein interactions in any organism using orthogonal organic phase separation (OOPS)

Existing high-throughput methods to identify RNA-binding proteins (RBPs) are based on capture of polyadenylated RNAs and cannot recover proteins that interact with nonadenylated RNAs, including long noncoding RNA, pre-mRNAs and bacterial RNAs. We present orthogonal organic phase separation (OOPS), which does not require molecular tagging or capture of polyadenylated RNA, and apply it to recover cross-linked protein-RNA and free protein, or protein-bound RNA and free RNA, in an unbiased way. We validated OOPS in HEK293, U2OS and MCF10A human cell lines, and show that 96% of proteins recovered were bound to RNA. We show that all long RNAs can be cross-linked to proteins, and recovered 1,838 RBPs, including 926 putative novel RBPs. OOPS is approximately 100-fold more efficient than existing methods and can enable analyses of dynamic RNA-protein interactions. We also characterize dynamic changes in RNA-protein interactions in mammalian cells following nocodazole arrest, and present a bacterial RNA-interactome for Escherichia coli. OOPS is compatible with downstream proteomics and RNA sequencing, and can be applied in any organism.

The m6A-methylase complex recruits TREX and regulates mRNA export

N⁶-methyladenosine (m⁶A) is the most abundant internal modification of eukaryotic mRNA. This modification has previously been shown to alter the export kinetics for mRNAs though the molecular details surrounding this phenomenon remain poorly understood. Recruitment of the TREX mRNA export complex to mRNA is driven by transcription, 5' capping and pre-mRNA splicing. Here we identify a fourth mechanism in human cells driving the association of TREX with mRNA involving the m⁶A methylase complex. We show that the m⁶A complex recruits TREX to m⁶A modified mRNAs and this process is essential for their efficient export. TREX also stimulates recruitment of the m⁶A reader protein YTHDC1 to the mRNA and the m⁶A complex influences the interaction of TREX with YTHDC1. Together our studies reveal a key role for TREX in the export of m⁶A modified mRNAs.

Foxp3+ T reg cells control psoriasiform inflammation by restraining an IFN-I-driven CD8+ T cell response

Psoriasis is a complex inflammatory skin disease affecting ∼3% of the population worldwide. Although type I interferons (IFN-I) are thought to be involved in its pathogenesis, the details of this relationship remain elusive. Here we show that in a murine model of imiquimod-driven psoriatic skin inflammation, Foxp3+ regulatory T cells (T reg cells) control inflammation severity by restraining IFN-I. Depletion of T reg cells induces IFN-I and IFN-stimulated gene expression, and leads to accumulation of CD8+ T cells in lesional skin. Mononuclear phagocytes (MNPs) were the source of IFN-I, and their depletion reversed the effect of T reg cell depletion. Blockade of IFN-I signaling abolished CD8+ T cell infiltration and excess inflammation in the skin of T reg cell-depleted mice. Depletion of CD8+ T cells attenuated pathology, confirming their role as critical effector cells downstream of IFN-I. Our results describe an unexpected role for T reg cells in restraint of an MNP-IFN-I-driven CD8+ T cell response during psoriasiform skin inflammation. These findings highlight a pathway with potential relevance for the treatment of early-stage disease.

CpG island composition differences are a source of gene expression noise indicative of promoter responsiveness

BACKGROUND

Population phenotypic variation can arise from genetic differences between individuals, or from cellular heterogeneity in an isogenic group of cells or organisms. The emergence of gene expression differences between genetically identical cells is referred to as gene expression noise, the sources of which are not well understood.

RESULTS

In this work, by studying gene expression noise between multiple cell lineages and mammalian species, we find consistent evidence of a role for CpG islands as sources of gene expression noise. Variation in noise among CpG island promoters can be partially attributed to differences in island size, in which short islands have noisier gene expression. Building on these findings, we investigate the potential for short CpG islands to act as fast response elements to environmental stimuli. Specifically, we find that these islands are enriched amongst primary response genes in SWI/SNF-independent stimuli, suggesting that expression noise is an indicator of promoter responsiveness.

CONCLUSIONS

Thus, through the integration of single-cell RNA expression profiling, chromatin landscape and temporal gene expression dynamics, we have uncovered a role for short CpG island promoters as fast response elements.

The Light Chain IgLV3-21 Defines a New Poor Prognostic Subgroup in Chronic Lymphocytic Leukemia: Results of a Multicenter Study

Purpose: Unmutated (UM) immunoglobulin heavy chain variable region (IgHV) status or IgHV3-21 gene usage is associated with poor prognosis in chronic lymphocytic leukemia (CLL) patients. Interestingly, IgHV3-21 is often co-expressed with light chain IgLV3-21, which is potentially able to trigger cell-autonomous BCR-mediated signaling. However, this light chain has never been characterized independently of the heavy chain IgHV3-21.Experimental Design: We performed total RNA sequencing in 32 patients and investigated IgLV3-21 prognostic impact in terms of treatment-free survival (TFS) and overall survival (OS) in 3 other independent cohorts for a total of 813 patients. IgLV3-21 presence was tested by real-time PCR and confirmed by Sanger sequencing.Results: Using total RNA sequencing to characterize 32 patients with high-risk CLL, we found a high frequency (28%) of IgLV3-21 rearrangements. Gene set enrichment analysis revealed that these patients express higher levels of genes responsible for ribosome biogenesis and translation initiation (P < 0.0001) as well as MYC target genes (P = 0.0003). Patients with IgLV3-21 rearrangements displayed a significantly shorter TFS and OS (P < 0.05), particularly those with IgHV mutation. In each of the three independent validation cohorts, we showed that IgLV3-21 rearrangements-similar to UM IgHV status-conferred poor prognosis compared with mutated IgHV (P < 0.0001). Importantly, we confirmed by multivariate analysis that this was independent of IgHV mutational status or subset #2 stereotyped receptor (P < 0.0001).Conclusions: We have demonstrated for the first time that a light chain can affect CLL prognosis and that IgLV3-21 light chain usage defines a new subgroup of CLL patients with poor prognosis. Clin Cancer Res; 24(20); 5048-57. ©2018 AACR.

Investigation into the role of the germline epigenome in the transmission of glucocorticoid-programmed effects across generations

Background

Early life exposure to adverse environments affects cardiovascular and metabolic systems in the offspring. These programmed effects are transmissible to a second generation through both male and female lines, suggesting germline transmission. We have previously shown that prenatal overexposure to the synthetic glucocorticoid dexamethasone (Dex) in rats reduces birth weight in the first generation (F1), a phenotype which is transmitted to a second generation (F2), particularly through the male line. We hypothesize that Dex exposure affects developing germ cells, resulting in transmissible alterations in DNA methylation, histone marks and/or small RNA in the male germline.

Results

We profile epigenetic marks in sperm from F1 Sprague Dawley rats expressing a germ cell-specific GFP transgene following Dex or vehicle treatment of the mothers, using methylated DNA immunoprecipitation sequencing, small RNA sequencing and chromatin immunoprecipitation sequencing for H3K4me3, H3K4me1, H3K27me3 and H3K9me3. Although effects on birth weight are transmitted to the F2 generation through the male line, no differences in DNA methylation, histone modifications or small RNA were detected between germ cells and sperm from Dex-exposed animals and controls.

Conclusions

Although the phenotype is transmitted to a second generation, we are unable to detect specific changes in DNA methylation, common histone modifications or small RNA profiles in sperm. Dex exposure is associated with more variable 5mC levels, particularly at non-promoter loci. Although this could be one mechanism contributing to the observed phenotype, other germline epigenetic modifications or non-epigenetic mechanisms may be responsible for the transmission of programmed effects across generations in this model.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1422-4) contains supplementary material, which is available to authorized users.

Punctuated evolution of canonical genomic aberrations in uveal melanoma

Cancer is thought to arise through the accumulation of genomic aberrations evolving under Darwinian selection. However, it remains unclear when the aberrations associated with metastasis emerge during tumor evolution. Uveal melanoma (UM) is the most common primary eye cancer and frequently leads to metastatic death, which is strongly linked to BAP1 mutations. Accordingly, UM is ideally suited for studying the clonal evolution of metastatic competence. Here we analyze sequencing data from 151 primary UM samples using a customized bioinformatic pipeline, to improve detection of BAP1 mutations and infer the clonal relationships among genomic aberrations. Strikingly, we find BAP1 mutations and other canonical genomic aberrations usually arise in an early punctuated burst, followed by neutral evolution extending to the time of clinical detection. This implies that the metastatic proclivity of UM is "set in stone" early in tumor evolution and may explain why advances in primary treatment have not improved survival.

Inhibition of histone H3K27 demethylases selectively modulates inflammatory phenotypes of natural killer cells

Natural killer (NK) cells are innate lymphocytes, important in immune surveillance and elimination of stressed, transformed, or virus-infected cells. They critically shape the inflammatory cytokine environment to orchestrate interactions of cells of the innate and adaptive immune systems. Some studies have reported that NK cell activation and cytokine secretion are controlled epigenetically but have yielded only limited insight into the mechanisms. Using chemical screening with small-molecule inhibitors of chromatin methylation and acetylation, further validated by knockdown approaches, we here identified Jumonji-type histone H3K27 demethylases as key regulators of cytokine production in human NK cell subsets. The prototypic JMJD3/UTX (Jumonji domain–containing protein 3) H3K27 demethylase inhibitor GSK-J4 increased global levels of the repressive H3K27me3 mark around transcription start sites of effector cytokine genes. Moreover, GSK-J4 reduced IFN-γ, TNFα, granulocyte–macrophage colony-stimulating factor (GM-CSF), and interleukin-10 levels in cytokine-stimulated NK cells while sparing their cytotoxic killing activity against cancer cells. The anti-inflammatory effect of GSK-J4 in NK cell subsets, isolated from peripheral blood or tissue from individuals with rheumatoid arthritis (RA), coupled with an inhibitory effect on formation of bone-resorbing osteoclasts, suggested that histone demethylase inhibition has broad utility for modulating immune and inflammatory responses. Overall, our results indicate that H3K27me3 is a dynamic and important epigenetic modification during NK cell activation and that JMJD3/UTX-driven H3K27 demethylation is critical for NK cell function.

N6-methyladenosine demethylase FTO targets pre-mRNAs and regulates alternative splicing and 3'-end processing

N6-methyladenosine (m⁶A) is the most abundant base modification found in messenger RNAs (mRNAs). The discovery of FTO as the first m⁶A mRNA demethylase established the concept of reversible RNA modification. Here, we present a comprehensive transcriptome-wide analysis of RNA demethylation and uncover FTO as a potent regulator of nuclear mRNA processing events such as alternative splicing and 3΄ end mRNA processing. We show that FTO binds preferentially to pre-mRNAs in intronic regions, in the proximity of alternatively spliced (AS) exons and poly(A) sites. FTO knockout (KO) results in substantial changes in pre-mRNA splicing with prevalence of exon skipping events. The alternative splicing effects of FTO KO anti-correlate with METTL3 knockdown suggesting the involvement of m⁶A. Besides, deletion of intronic region that contains m⁶A-linked DRACH motifs partially rescues the FTO KO phenotype in a reporter system. All together, we demonstrate that the splicing effects of FTO are dependent on the catalytic activity in vivo and are mediated by m⁶A. Our results reveal for the first time the dynamic connection between FTO RNA binding and demethylation activity that influences several mRNA processing events.

U1 snRNP telescripting regulates a size-function-stratified human genome

U1 snRNP (U1) functions in splicing introns and telescripting, which suppresses premature cleavage and polyadenylation (PCPA). Using U1 inhibition in human cells, we show that U1 telescripting is selectively required for sustaining long-distance transcription elongation in introns of large genes (median 39 kb). Evidence of widespread PCPA in the same locations in normal tissues reveals that large genes incur natural transcription attrition. Underscoring the importance of U1 telescripting as a gene-size-based mRNA-regulation mechanism, small genes were not sensitive to PCPA, and the spliced-mRNA productivity of ∼1,000 small genes (median 6.8 kb) increased upon U1 inhibition. Notably, these small, upregulated genes were enriched in functions related to acute stimuli and cell-survival response, whereas genes subject to PCPA were enriched in cell-cycle progression and developmental functions. This gene size-function polarization increased in metazoan evolution by enormous intron expansion. We propose that telescripting adds an overarching layer of regulation to size-function-stratified genomes, leveraged by selective intron expansion to rapidly shift gene expression priorities.

High concordance between Illumina HiSeq2500 and NextSeq500 for reduced representation bisulfite sequencing (RRBS)

Reduced representation bisulfite sequencing (RRBS) provides an efficient method for measuring DNA methylation at single base resolution in regions of high CpG density. This technique has been extensively tested on the HiSeq2500, which uses a 4-colour detection method, however it is unclear if the method will also work on the NextSeq500 platform, which employs a 2-colour detection system. We created an RRBS library and sequenced it on both the HiSeq2500 and NextSeq500, and found no significant difference in the base composition of reads derived from either machine. Moreover, the methylation calls made from the data of each instrument were highly concordant, with methylation patterns across the genome appearing as expected. Therefore, RRBS can be sequenced on the Nextseq500 with comparable quality to that of the HiSeq2500. All sequencing data are deposited in the GEO database under accession number GSE87097.

Tuning the Transcriptional Response to Hypoxia by Inhibiting Hypoxia-inducible Factor (HIF) Prolyl and Asparaginyl Hydroxylases

The hypoxia-inducible factor (HIF) system orchestrates cellular responses to hypoxia in animals. HIF is an α/β-heterodimeric transcription factor that regulates the expression of hundreds of genes in a tissue context-dependent manner. The major hypoxia-sensing component of the HIF system involves oxygen-dependent catalysis by the HIF hydroxylases; in humans there are three HIF prolyl hydroxylases (PHD1-3) and an asparaginyl hydroxylase (factor-inhibiting HIF (FIH)). PHD catalysis regulates HIFα levels, and FIH catalysis regulates HIF activity. How differences in HIFα hydroxylation status relate to variations in the induction of specific HIF target gene transcription is unknown. We report studies using small molecule HIF hydroxylase inhibitors that investigate the extent to which HIF target gene expression is induced by PHD or FIH inhibition. The results reveal substantial differences in the role of prolyl and asparaginyl hydroxylation in regulating hypoxia-responsive genes in cells. PHD inhibitors with different structural scaffolds behave similarly. Under the tested conditions, a broad-spectrum 2-oxoglutarate dioxygenase inhibitor is a better mimic of the overall transcriptional response to hypoxia than the selective PHD inhibitors, consistent with an important role for FIH in the hypoxic transcriptional response. Indeed, combined application of selective PHD and FIH inhibitors resulted in the transcriptional induction of a subset of genes not fully responsive to PHD inhibition alone. Thus, for the therapeutic regulation of HIF target genes, it is important to consider both PHD and FIH activity, and in the case of some sets of target genes, simultaneous inhibition of the PHDs and FIH catalysis may be preferable.

High-fat diet disrupts metabolism in two generations of rats in a parent-of-origin specific manner

Experimental and epidemiological evidence demonstrate that ancestral diet might contribute towards offspring health. This suggests that nutrition may be able to modify genetic or epigenetic information carried by germ cells (GCs). To examine if a parental high fat diet (HFD) influences metabolic health in two generations of offspring, GC-eGFP Sprague Dawley rats were weaned onto HFD (45% fat) or Control Diet (CD; 10% fat). At 19 weeks, founders (F0) were bred with controls, establishing the F1 generation. HFD resulted in 9.7% and 14.7% increased weight gain in male and female F0 respectively. F1 offspring of HFD mothers and F1 daughters of HFD-fed fathers had increased weight gain compared to controls. F1 rats were bred with controls at 19 weeks to generate F2 offspring. F2 male offspring derived from HFD-fed maternal grandfathers exhibited increased adiposity, plasma leptin and luteinising hormone to testosterone ratio. Despite transmission via the founding male germline, we did not find significant changes in the F0 intra-testicular GC transcriptome. Thus, HFD consumption by maternal grandfathers results in a disrupted metabolic and reproductive hormone phenotype in grandsons in the absence of detectable changes in the intra-testicular GC transcriptome.

Defining the microbial transcriptional response to colitis through integrated host and microbiome profiling

The gut microbiome is significantly altered in inflammatory bowel diseases, but the basis of these changes is not well understood. We have combined metagenomic and metatranscriptomic profiling of the gut microbiome to assess modifications to both bacterial community structure and transcriptional activity in a mouse model of colitis. By using transcriptomic analysis of colonic tissue and luminal RNA derived from the host, we have also characterised how host transcription relates to the microbial transcriptional response in inflammation. In colitis, increased abundance and transcription of diverse microbial gene families involved in responses to nutrient deprivation, antimicrobial peptide production and oxidative stress support an adaptation of multiple commensal genera to withstand a diverse set of environmental stressors in the inflammatory environment. These data are supported by a transcriptional signature of activated macrophages and granulocytes in the gut lumen during colitis, a signature that includes the transcription of the key antimicrobial genes S100a8 and S100a9 (calprotectin). Genes involved in microbial resistance to oxidative stress, including Dps/ferritin, Fe-dependent peroxidase and glutathione S-transferase were identified as changing to a greater extent at the level of transcription than would be predicted by DNA abundance changes, implicating a role for increased oxygen tension and/or host-derived reactive oxygen species in driving transcriptional changes in commensal microbes.

Temporal transcriptomics suggest that twin-peaking genes reset the clock

The mammalian suprachiasmatic nucleus (SCN) drives daily rhythmic behavior and physiology, yet a detailed understanding of its coordinated transcriptional programmes is lacking. To reveal the finer details of circadian variation in the mammalian SCN transcriptome we combined laser-capture microdissection (LCM) and RNA-seq over a 24 hr light / dark cycle. We show that 7-times more genes exhibited a classic sinusoidal expression signature than previously observed in the SCN. Another group of 766 genes unexpectedly peaked twice, near both the start and end of the dark phase; this twin-peaking group is significantly enriched for synaptic transmission genes that are crucial for light-induced phase shifting of the circadian clock. 341 intergenic non-coding RNAs, together with novel exons of annotated protein-coding genes, including Cry1, also show specific circadian expression variation. Overall, our data provide an important chronobiological resource (www.wgpembroke.com/shiny/SCNseq/) and allow us to propose that transcriptional timing in the SCN is gating clock resetting mechanisms.

DOI:http://dx.doi.org/10.7554/eLife.10518.001

The daily cycles of life in mammals are driven by a small region of the brain called the suprachiasmatic nucleus (or SCN). The SCN receives signals from sunlight and other environmental factors to help coordinate most aspects of daily biological activity and behaviour. To work correctly, it is essential that the SCN switches certain genes on and off at exactly the right time. However, many questions remain over the identity of these genes and how their levels of activity change during a 24-hour period.

When a gene is active (or “being expressed”), it is used as a template to build the molecules of RNA that are needed to make proteins and to help to control how cells work. Pembroke et al. have now sequenced the RNA molecules made in the SCN of mice (which plays the same role as the equivalent human brain region) over a 24-hour period. The mice spent half of each day in the light, and half in the dark. This revealed that the expression levels of over a quarter of all the genes that are found in the SCN fluctuate over a 24-hour period. One particular group of genes peak in activity twice a day; Pembroke et al. suggest that these genes are important for controlling how an animal can adjust its body clock to light.

Further research is now needed to find out which of the newly discovered fluctuating genes play the most important roles in daily activity rhythms, and which might play a part in disease.

DOI:http://dx.doi.org/10.7554/eLife.10518.002

CGAT: a model for immersive personalized training in computational genomics

How should the next generation of genomics scientists be trained while simultaneously pursuing high quality and diverse research? CGAT, the Computational Genomics Analysis and Training programme, was set up in 2010 by the UK Medical Research Council to complement its investment in next-generation sequencing capacity. CGAT was conceived around the twin goals of training future leaders in genome biology and medicine, and providing much needed capacity to UK science for analysing genome scale data sets. Here we outline the training programme employed by CGAT and describe how it dovetails with collaborative research projects to launch scientists on the road towards independent research careers in genomics.

Population and single-cell genomics reveal the Aire dependency, relief from Polycomb silencing, and distribution of self-antigen expression in thymic epithelia

Promiscuous gene expression (PGE) by thymic epithelial cells (TEC) is essential for generating a diverse T cell antigen receptor repertoire tolerant to self-antigens, and thus for avoiding autoimmunity. Nevertheless, the extent and nature of this unusual expression program within TEC populations and single cells are unknown. Using deep transcriptome sequencing of carefully identified mouse TEC subpopulations, we discovered a program of PGE that is common between medullary (m) and cortical TEC, further elaborated in mTEC, and completed in mature mTEC expressing the autoimmune regulator gene (Aire). TEC populations are capable of expressing up to 19,293 protein-coding genes, the highest number of genes known to be expressed in any cell type. Remarkably, in mouse mTEC, Aire expression alone positively regulates 3980 tissue-restricted genes. Notably, the tissue specificities of these genes include known targets of autoimmunity in human AIRE deficiency. Led by the observation that genes induced by Aire expression are generally characterized by a repressive chromatin state in somatic tissues, we found these genes to be strongly associated with H3K27me3 marks in mTEC. Our findings are consistent with AIRE targeting and inducing the promiscuous expression of genes previously epigenetically silenced by Polycomb group proteins. Comparison of the transcriptomes of 174 single mTEC indicates that genes induced by Aire expression are transcribed stochastically at low cell frequency. Furthermore, when present, Aire expression-dependent transcript levels were 16-fold higher, on average, in individual TEC than in the mTEC population.

Personalized medicine in hematology - A landmark from bench to bed

Personalized medicine is the cornerstone of medical practice. It tailors treatments for specific conditions of an affected individual. The borders of personalized medicine are defined by limitations in technology and our understanding of biology, physiology and pathology of various conditions. Current advances in technology have provided physicians with the tools to investigate the molecular makeup of the disease. Translating these molecular make-ups to actionable targets has led to the development of small molecular inhibitors. Also, detailed understanding of genetic makeup has allowed us to develop prognostic markers, better known as companion diagnostics. Current attempts in the development of drug delivery systems offer the opportunity of delivering specific inhibitors to affected cells in an attempt to reduce the unwanted side effects of drugs.