Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

Gene expression plasticity as a mechanism of coral adaptation to a variable environment

Local adaptation is ubiquitous, but the molecular mechanisms that give rise to this ecological phenomenon remain largely unknown. A year-long reciprocal transplant of mustard hill coral (Porites astreoides) between a highly environmentally variable inshore habitat and a more stable offshore habitat demonstrated that populations exhibit phenotypic signatures that are consistent with local adaptation. We characterized the genomic basis of this adaptation in both coral hosts and their intracellular symbionts (Symbiodinium sp.) using genome-wide gene expression profiling. Populations differed primarily in their capacity for plasticity: following transplantation to a novel environment, inshore-origin coral expression profiles became significantly more similar to the local population's profiles than those in offshore-origin corals. Furthermore, elevated plasticity of the environmental stress response expression was correlated with lower susceptibility to a natural summer bleaching event, suggesting that plasticity is adaptive in the inshore environment. Our results reveal a novel genomic mechanism of resilience to a variable environment, demonstrating that corals are capable of a more diverse molecular response to stress than previously thought.

Piphillin: Improved Prediction of Metagenomic Content by Direct Inference from Human Microbiomes

Functional analysis of a clinical microbiome facilitates the elucidation of mechanisms by which microbiome perturbation can cause a phenotypic change in the patient. The direct approach for the analysis of the functional capacity of the microbiome is via shotgun metagenomics. An inexpensive method to estimate the functional capacity of a microbial community is through collecting 16S rRNA gene profiles then indirectly inferring the abundance of functional genes. This inference approach has been implemented in the PICRUSt and Tax4Fun software tools. However, those tools have important limitations since they rely on outdated functional databases and uncertain phylogenetic trees and require very specific data pre-processing protocols. Here we introduce Piphillin, a straightforward algorithm independent of any proposed phylogenetic tree, leveraging contemporary functional databases and not obliged to any singular data pre-processing protocol. When all three inference tools were evaluated against actual shotgun metagenomics, Piphillin was superior in predicting gene composition in human clinical samples compared to both PICRUSt and Tax4Fun (p<0.01 and p<0.001, respectively) and Piphillin’s ability to predict disease associations with specific gene orthologs exhibited a 15% increase in balanced accuracy compared to PICRUSt. From laboratory animal samples, no performance advantage was observed for any one of the tools over the others and for environmental samples all produced unsatisfactory predictions. Our results demonstrate that functional inference using the direct method implemented in Piphillin is preferable for clinical biospecimens. Piphillin is publicly available for academic use at http://secondgenome.com/Piphillin.

De novo assembly and analysis of changes in the protein-coding transcriptome of the freshwater shrimp Paratya australiensis (Decapoda: Atyidae) in response to acid sulfate drainage water

Background

The atyid shrimp Paratya australiensis occurs in surface freshwater habitats throughout eastern Australia and has been used to study the ecotoxicology of contaminants such as pesticides and metals. The acidification of surface water that can occur after acid sulfate material in soils and sediments is oxidised and subsequently re-wetted is a serious environmental issue in coastal regions and inland riverine floodplains worldwide. Solubilisation of soil-associated minerals can result in high waterborne concentrations of mineral salts and dissolved metals, which together with low pH represent a potential threat to aquatic ecosystems in affected regions. The aims of the present study were to gain insight into stress responses induced by exposure to acid drainage water (ADW) in P. australiensis by determining changes in the abundance of protein-coding transcripts and to generate a comprehensive transcriptomic resource to facilitate further research into gene regulation or protein structure and function in this species. Adult P. australiensis were exposed for 24 h to undiluted ADW, 50 % ADW diluted in river water, or to river water as control, and high-throughput mRNA sequencing (RNA-Seq) conducted on whole-body tissues. A reference transcriptome was generated using de novo assembly and putative protein-coding regions were identified and annotated. Changes in transcript abundance in response to ADW exposure were determined by aligning reads to the reference transcriptome and quantifying coverage.

Results

A high proportion of arthropod benchmarking universal single-copy orthologues were present in the reference transcriptome. Functions associated with cuticle biosynthesis and oxidative stress were significantly enriched in the lists of transcripts exhibiting differential abundance in either direction after exposure to 50 % or 100 % ADW. Transcripts involved in osmoregulation exhibited decreased abundance following exposure to ADW. The transcriptome contained full-length coding sequences for numerous proteins known to be involved in environmental response pathways, including two putative metallothioneins, four glutathione peroxidases and 19 nuclear receptors.

Conclusions

The results of the present study provide insight into stress response pathways induced in crustaceans by short-term exposure to multiple stressors present in ADW such as low pH, high salinity and dissolved metals, and represent a resource for future toxicogenomics and protein functional studies in P. australiensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3208-y) contains supplementary material, which is available to authorized users.

Recessive mutations in the kinase ZAK cause a congenital myopathy with fibre type disproportion

Congenital myopathies define a heterogeneous group of neuromuscular diseases with neonatal or childhood hypotonia and muscle weakness. The genetic cause is still unknown in many patients, precluding genetic counselling and better understanding of the physiopathology. To identify novel genetic causes of congenital myopathies, exome sequencing was performed in three consanguineous families. We identified two homozygous frameshift mutations and a homozygous nonsense mutation in the mitogen-activated protein triple kinase ZAK. In total, six affected patients carry these mutations. Reverse transcription polymerase chain reaction and transcriptome analyses suggested nonsense mRNA decay as a main impact of mutations. The patients demonstrated a generalized slowly progressive muscle weakness accompanied by decreased vital capacities. A combination of proximal contractures with distal joint hyperlaxity is a distinct feature in one family. The low endurance and compound muscle action potential amplitude were strongly ameliorated on treatment with anticholinesterase inhibitor in another patient. Common histopathological features encompassed fibre size variation, predominance of type 1 fibre and centralized nuclei. A peculiar subsarcolemmal accumulation of mitochondria pointing towards the centre of the fibre was a novel histological hallmark in one family. These findings will improve the molecular diagnosis of congenital myopathies and implicate the mitogen-activated protein kinase (MAPK) signalling as a novel pathway altered in these rare myopathies.

Conditional estimation of local pooled dispersion parameter in small-sample RNA-Seq data improves differential expression test

High throughput sequencing technology in transcriptomics studies contribute to the understanding of gene regulation mechanism and its cellular function, but also increases a need for accurate statistical methods to assess quantitative differences between experiments. Many methods have been developed to account for the specifics of count data: non-normality, a dependence of the variance on the mean, and small sample size. Among them, the small number of samples in typical experiments is still a challenge. Here we present a method for differential analysis of count data, using conditional estimation of local pooled dispersion parameters. A comprehensive evaluation of our proposed method in the aspect of differential gene expression analysis using both simulated and real data sets shows that the proposed method is more powerful than other existing methods while controlling the false discovery rates. By introducing conditional estimation of local pooled dispersion parameters, we successfully overcome the limitation of small power and enable a powerful quantitative analysis focused on differential expression test with the small number of samples.

The founder-cell transcriptome in the Arabidopsis apetala1 cauliflower inflorescence meristem

BACKGROUND

Although the pattern of lateral organ formation from apical meristems establishes species-specific plant architecture, the positional information that confers cell fate to cells as they transit to the meristem flanks where they differentiate, remains largely unknown. We have combined fluorescence-activated cell sorting and RNA-seq to characterise the cell-type-specific transcriptome at the earliest developmental time-point of lateral organ formation using DORNRÖSCHEN-LIKE::GFP to mark founder-cell populations at the periphery of the inflorescence meristem (IM) in apetala1 cauliflower double mutants, which overproliferate IMs.

RESULTS

Within the lateral organ founder-cell population at the inflorescence meristem, floral primordium identity genes are upregulated and stem-cell identity markers are downregulated. Additional differentially expressed transcripts are involved in polarity generation and boundary formation, and in epigenetic and post-translational changes. However, only subtle transcriptional reprogramming within the global auxin network was observed.

CONCLUSIONS

The transcriptional network of differentially expressed genes supports the hypothesis that lateral organ founder-cell specification involves the creation of polarity from the centre to the periphery of the IM and the establishment of a boundary from surrounding cells, consistent with bract initiation. However, contrary to the established paradigm that sites of auxin response maxima pre-pattern lateral organ initiation in the IM, auxin response might play a minor role in the earliest stages of lateral floral initiation.

RNA Sequencing of Contaminated Seeds Reveals the State of the Seed Permissive for Pre-Harvest Aflatoxin Contamination and Points to a Potential Susceptibility Factor

Pre-harvest aflatoxin contamination (PAC) is a major problem facing peanut production worldwide. Produced by the ubiquitous soil fungus, Aspergillus flavus, aflatoxin is the most naturally occurring known carcinogen. The interaction between fungus and host resulting in PAC is complex, and breeding for PAC resistance has been slow. It has been shown that aflatoxin production can be induced by applying drought stress as peanut seeds mature. We have implemented an automated rainout shelter that controls temperature and moisture in the root and peg zone to induce aflatoxin production. Using polymerase chain reaction (PCR) and high performance liquid chromatography (HPLC), seeds meeting the following conditions were selected: infected with Aspergillus flavus and contaminated with aflatoxin; and not contaminated with aflatoxin. RNA sequencing analysis revealed groups of genes that describe the transcriptional state of contaminated vs. uncontaminated seed. These data suggest that fatty acid biosynthesis and abscisic acid (ABA) signaling are altered in contaminated seeds and point to a potential susceptibility factor, ABR1, as a repressor of ABA signaling that may play a role in permitting PAC.

Insecticide resistance mediated by an exon skipping event

Many genes increase coding capacity by alternate exon usage. The gene encoding the insect nicotinic acetylcholine receptor (nAChR) α6 subunit, target of the bio‐insecticide spinosad, is one example of this and expands protein diversity via alternative splicing of mutually exclusive exons. Here, we show that spinosad resistance in the tomato leaf miner, Tuta absoluta is associated with aberrant regulation of splicing of Taα6 resulting in a novel form of insecticide resistance mediated by exon skipping. Sequencing of the α6 subunit cDNA from spinosad selected and unselected strains of T. absoluta revealed all Taα6 transcripts of the selected strain were devoid of exon 3, with comparison of genomic DNA and mRNA revealing this is a result of exon skipping. Exon skipping cosegregated with spinosad resistance in survival bioassays, and functional characterization of this alteration using modified human nAChR α7, a model of insect α6, demonstrated that exon 3 is essential for receptor function and hence spinosad sensitivity. DNA and RNA sequencing analyses suggested that exon skipping did not result from genetic alterations in intronic or exonic cis‐regulatory elements, but rather was associated with a single epigenetic modification downstream of exon 3a, and quantitative changes in the expression of trans‐acting proteins that have known roles in the regulation of alternative splicing. Our results demonstrate that the intrinsic capacity of the α6 gene to generate transcript diversity via alternative splicing can be readily exploited during the evolution of resistance and identifies exon skipping as a molecular alteration conferring insecticide resistance.

De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection

Background

Piscine reovirus (PRV) has been associated with the serious disease known as Heart and Skeletal Muscle Inflammation (HSMI) in cultured Atlantic salmon Salmo salar in Norway. PRV is also prevalent in wild and farmed salmon without overt disease manifestations, suggesting multifactorial triggers or PRV variant-specific factors are required to initiate disease. In this study, we explore the head kidney transcriptome of Sockeye salmon Oncorhynchus nerka during early PRV infection to identify host responses in the absence of disease in hopes of elucidating mechanisms by which PRV may directly alter host functions and contribute to the development of a disease state. We further investigate the role of PRV as a coinfecting agent following superinfection with infectious hematopoietic necrosis virus (IHNV) – a highly pathogenic rhabdovirus endemic to the west coast of North America.

Results

Challenge of Sockeye salmon with PRV resulted in high quantities of viral transcripts to become present in the blood and kidney of infected fish without manifestations of disease. De novo transcriptome assembly of over 2.3 billion paired RNA-seq reads from the head kidneys of 36 fish identified more than 320,000 putative unigenes, of which less than 20 were suggested to be differentially expressed in response to PRV at either 2 or 3 weeks post challenge by DESeq2 and edgeR analysis. Of these, only one, Ependymin, was confirmed to be differentially expressed by qPCR in an expanded sample set. In contrast, IHNV induced substantial transcriptional changes (differential expression of > 20,000 unigenes) which included transcripts involved in antiviral and inflammatory response pathways. Prior infection with PRV had no significant effect on host responses to superinfecting IHNV, nor did host responses initiated by IHNV exposure influence increasing PRV loads.

Conclusions

PRV does not substantially alter the head kidney transcriptome of Sockeye salmon during early (2 to 3 week) infection and dissemination in a period of significant increasing viral load, nor does the presence of PRV change the host transcriptional response to an IHNV superinfection. Further, concurrent infections of PRV and IHNV do not appear to significantly influence the infectivity or severity of IHNV associated disease, or conversely, PRV load.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3196-y) contains supplementary material, which is available to authorized users.

Two-miRNA classifiers differentiate mutation-negative follicular thyroid carcinomas and follicular thyroid adenomas in fine needle aspirations with high specificity

Diagnosis of thyroid by fine needle aspiration is challenging for the "indeterminate" category and can be supported by molecular testing. We set out to identify miRNA markers that could be used in a diagnostic setting to improve the discrimination of mutation-negative indeterminate fine needle aspirations. miRNA high-throughput sequencing was performed for freshly frozen tissue samples of 19 RAS and PAX8/PPARG mutation-negative follicular thyroid carcinomas, and 23 RAS and PAX8/PPARG mutation-negative follicular adenomas. Differentially expressed miRNAs were validated by quantitative polymerase chain reaction in a set of 44 fine needle aspiration samples representing 24 follicular thyroid carcinomas and 20 follicular adenomas. Twenty-six miRNAs characterized by a significant differential expression between follicular thyroid carcinomas and follicular adenomas were identified. Nevertheless, since no single miRNA had satisfactory predictive power, classifiers comprising two differentially expressed miRNAs were designed with the aim to improve the classification. Six two-miRNA classifiers were established and quantitative polymerase chain reaction validated in fine needle aspiration samples. Four out of six classifiers were characterized by a high specificity (≥94 %). The best two-miRNA classifier (miR-484/miR-148b-3p) identified thyroid malignancy with a sensitivity of 89 % and a specificity of 87 %. The high-throughput sequencing allowed the identification of subtle differences in the miRNA expression profiles of follicular thyroid carcinomas and follicular adenomas. While none of the differentially expressed miRNAs could be used as a stand-alone malignancy marker, the validation results for two-miRNA classifiers in an independent set of fine needle aspirations are very promising. The ultimate evaluation of these classifiers for their capability of discriminating mutation-negative indeterminate fine needle aspirations will require the evaluation of a sufficiently large number of fine needle aspirations with histological confirmation.

Neural Stem Cells Derived from Human Parthenogenetic Stem Cells Engraft and Promote Recovery in a Nonhuman Primate Model of Parkinson's Disease

Cell therapy has attracted considerable interest as a promising therapeutic alternative for patients with Parkinson's disease (PD). Clinical studies have shown that grafted fetal neural tissue can achieve considerable biochemical and clinical improvements in PD. However, the source of fetal tissue grafts is limited and ethically controversial. Human parthenogenetic stem cells offer a good alternative because they are derived from unfertilized oocytes without destroying potentially viable human embryos and can be used to generate an unlimited supply of neural cells for transplantation. We have previously reported that human parthenogenetic stem cell-derived neural stem cells (hpNSCs) successfully engraft, survive long term, and increase brain dopamine (DA) levels in rodent and nonhuman primate models of PD. Here we report the results of a 12-month transplantation study of hpNSCs in 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-lesioned African green monkeys with moderate to severe clinical parkinsonian symptoms. The hpNSCs manufactured under current good manufacturing practice (cGMP) conditions were injected bilaterally into the striatum and substantia nigra of immunosuppressed monkeys. Transplantation of hpNSCs was safe and well tolerated by the animals with no dyskinesia, tumors, ectopic tissue formation, or other test article-related serious adverse events. We observed that hpNSCs promoted behavioral recovery; increased striatal DA concentration, fiber innervation, and number of dopaminergic neurons; and induced the expression of genes and pathways downregulated in PD compared to vehicle control animals. These results provide further evidence for the clinical translation of hpNSCs and support the approval of the world's first pluripotent stem cell-based phase I/IIa study for the treatment of PD (Clinical Trial Identifier NCT02452723).

Increased irritability, anxiety, and immune reactivity in transgenic Huntington's disease monkeys

Although the most notable clinical symptoms of Huntington's disease (HD) are motor disturbances and brain atrophy, other symptoms include cognitive dysfunction, emotional and hormonal dysregulation. Emotional dysregulation (irritability, anger/aggression, and anxiety) and increased inflammation are early emerging symptoms which can be detected decades before the onset of motor symptoms in HD patients. Despite the advances in understanding the genetic causes of HD there is still no cure or preventative treatment. Thus, to better understand the pathogenesis of HD and develop effective treatments, a holistic understanding of HD is needed, as well as animal models that replicate the full spectrum of HD symptoms. The current study examined the emotional, hormonal, and gene expression responses to an acute stressor of adult male transgenic HD rhesus monkeys (n=2) as compared to wild-type controls (n=2). Results revealed that HD monkeys expressed increased anxiety and irritability/aggression as compared to controls. Reactive cortisol response to the stressor was similar between groups. However, HD monkeys exhibited increased pro-inflammatory cytokines and higher induction of immune pathway genes as compared to controls. Overall, results reveal that HD monkeys exhibit these early emerging symptoms of HD and may be an effective animal model to facilitate the development of new therapeutics for HD patients.

An essential role for the IL-2 receptor in Treg cell function

Regulatory T cells (Treg cells), which have abundant expression of the interleukin 2 receptor (IL-2R), are reliant on IL-2 produced by activated T cells. This feature indicates a key role for a simple network based on the consumption of IL-2 by Treg cells in their suppressor function. However, congenital deficiency in IL-2R results in reduced expression of the Treg cell lineage-specification factor Foxp3, which has confounded experimental efforts to understand the role of IL-2R expression and signaling in the suppressor function of Treg cells. Using genetic gain- and loss-of-function approaches, we found that capture of IL-2 was dispensable for the control of CD4+ T cells but was important for limiting the activation of CD8+ T cells, and that IL-2R-dependent activation of the transcription factor STAT5 had an essential role in the suppressor function of Treg cells separable from signaling via the T cell antigen receptor.

A Two-Component Regulatory System Impacts Extracellular Membrane-Derived Vesicle Production in Group A Streptococcus

Export of macromolecules via extracellular membrane-derived vesicles (MVs) plays an important role in the biology of Gram-negative bacteria. Gram-positive bacteria have also recently been reported to produce MVs; however, the composition and mechanisms governing vesiculogenesis in Gram-positive bacteria remain undefined. Here, we describe MV production in the Gram-positive human pathogen group A streptococcus (GAS), the etiological agent of necrotizing fasciitis and streptococcal toxic shock syndrome. M1 serotype GAS isolates in culture exhibit MV structures both on the cell wall surface and in the near vicinity of bacterial cells. A comprehensive analysis of MV proteins identified both virulence-associated protein substrates of the general secretory pathway in addition to "anchorless surface proteins." Characteristic differences in the contents, distributions, and fatty acid compositions of specific lipids between MVs and GAS cell membrane were also observed. Furthermore, deep RNA sequencing of vesicular RNAs revealed that GAS MVs contained differentially abundant RNA species relative to bacterial cellular RNA. MV production by GAS strains varied in a manner dependent on an intact two-component system, CovRS, with MV production negatively regulated by the system. Modulation of MV production through CovRS was found to be independent of both GAS cysteine protease SpeB and capsule biosynthesis. Our data provide an explanation for GAS secretion of macromolecules, including RNAs, lipids, and proteins, and illustrate a regulatory mechanism coordinating this secretory response.

IMPORTANCE

Group A streptococcus (GAS) is a Gram-positive bacterial pathogen responsible for more than 500,000 deaths annually. Establishment of GAS infection is dependent on a suite of proteins exported via the general secretory pathway. Here, we show that GAS naturally produces extracellular vesicles with a unique lipid composition that are laden with proteins and RNAs. Interestingly, both virulence-associated proteins and RNA species were found to be differentially abundant in vesicles relative to the bacteria. Furthermore, we show that genetic disruption of the virulence-associated two-component regulator CovRS leads to an increase in vesicle production. This study comprehensively describes the protein, RNA, and lipid composition of GAS-secreted MVs and alludes to a regulatory system impacting this process.

The transcription factor, Nuclear factor, erythroid 2 (Nfe2), is a regulator of the oxidative stress response during Danio rerio development

Development is a complex and well-defined process characterized by rapid cell proliferation and apoptosis. At this stage in life, a developmentally young organism is more sensitive to toxicants as compared to an adult. In response to pro-oxidant exposure, members of the Cap'n'Collar (CNC) basic leucine zipper (b-ZIP) transcription factor family (including Nfe2 and Nfe2-related factors, Nrfs) activate the expression of genes whose protein products contribute to reduced toxicity. Here, we studied the role of the CNC protein, Nfe2, in the developmental response to pro-oxidant exposure in the zebrafish (Danio rerio). Following acute waterborne exposures to diquat or tert-buytlhydroperoxide (tBOOH) at one of three developmental stages, wildtype (WT) and nfe2 knockout (KO) embryos and larvae were morphologically scored and their transcriptomes sequenced. Early in development, KO animals suffered from hypochromia that was made more severe through exposure to pro-oxidants; this phenotype in the KO may be linked to decreased expression of alas2, a gene involved in heme synthesis. WT and KO eleutheroembryos and larvae were phenotypically equally affected by exposure to pro-oxidants, where tBOOH caused more pronounced phenotypes as compared to diquat. Comparing diquat and tBOOH exposed embryos relative to the WT untreated control, a greater number of genes were up-regulated in the tBOOH condition as compared to diquat (tBOOH: 304 vs diquat: 148), including those commonly found to be differentially regulated in the vertebrate oxidative stress response (OSR) (e.g. hsp70.2, txn1, and gsr). When comparing WT and KO across all treatments and times, there were 1170 genes that were differentially expressed, of which 33 are known targets of the Nrf proteins Nrf1 and Nrf2. More specifically, in animals exposed to pro-oxidants a total of 968 genes were differentially expressed between WT and KO across developmental time, representing pathways involved in coagulation, embryonic organ development, body fluid level regulation, erythrocyte differentiation, and oxidation-reduction, amongst others. The greatest number of genes that changed in expression between WT and KO occurred in animals exposed to diquat at 2h post fertilization (hpf). Across time and treatment, there were six genes (dhx40, cfap70, dnajb9b, slc35f4, spi-c, and gpr19) that were significantly up-regulated in KO compared to WT and four genes (fhad1, cyp4v7, nlrp12, and slc16a6a) that were significantly down-regulated. None of these genes have been previously identified as targets of Nfe2 or the Nrf family. These results demonstrate that the zebrafish Nfe2 may be a regulator of both primitive erythropoiesis and the OSR during development.

The drought response of potato reference cultivars with contrasting tolerance

Systems responses to drought stress of four potato reference cultivars with differential drought tolerance (Solanum tuberosum L.) were investigated by metabolome profiling and RNA sequencing. Systems analysis was based on independent field and greenhouse trials. Robust differential drought responses across all cultivars under both conditions comprised changes of proline, raffinose, galactinol, arabitol, arabinonic acid, chlorogenic acid and 102 transcript levels. The encoded genes contained a high proportion of heat shock proteins and proteins with signalling or regulatory functions, for example, a homolog of abscisic acid receptor PYL4. Constitutive differences of the tolerant compared with the sensitive cultivars included arbutin, octopamine, ribitol and 248 transcripts. The gene products of many of these transcripts were pathogen response related, such as receptor kinases, or regulatory proteins, for example, a homolog of the Arabidopsis FOUR LIPS MYB-regulator of stomatal cell proliferation. Functional enrichment analyses imply heat stress as a major acclimation component of potato leaves to long-term drought stress. Enhanced heat stress during drought can be caused by loss of transpiration cooling. This effect and CO₂ limitation are the main consequences of drought-induced or abscisic acid-induced stomatal closure. Constitutive differences in metabolite and transcript levels between tolerant and sensitive cultivars indicate interactions of drought tolerance and pathogen resistance in potato.

Chemokine (C-C motif) receptor 2-positive monocytes aggravate the early phase of acetaminophen-induced acute liver injury

Acetaminophen (APAP, paracetamol) poisoning is a leading cause of acute liver failure (ALF) in humans and induces hepatocyte necrosis, followed by activation of the innate immune system, further aggravating liver injury. The role of infiltrating monocytes during the early phase of ALF is still ambiguous. Upon experimental APAP overdose in mice, monocyte-derived macrophages (MoMFs) massively accumulated in injured liver within 12-24 hours, whereas the number of tissue-resident macrophages (Kupffer cells) decreased. Influx of MoMFs is dependent on the chemokine receptor, chemokine (C-C motif) receptor 2 (CCR2), given that Ccr2^-/- mice display reduced infiltration of monocytes and attenuated liver injury post-APAP overdose at early time points. As evidenced by intravital multiphoton microscopy of Ccr2 reporter mice, CCR2⁺ monocytes infiltrate liver as early as 8-12 hours post-APAP overdose and form dense cellular clusters around necrotic areas. CCR2⁺ MoMFs express a distinct pattern of inflammatory, but also repair-associated, genes in injured livers. Adoptive transfer experiments revealed that MoMFs primarily exert proinflammatory functions early post-APAP, thereby aggravating liver injury. Consequently, early pharmacological inhibition of either chemokine (C-C motif) ligand (CCL2; by the inhibitor, mNOX-E36) or CCR2 (by the orally available dual CCR2/CCR5 inhibitor, cenicriviroc) reduces monocyte infiltration and APAP-induced liver injury (AILI) in mice. Importantly, neither the early nor continuous inhibition of CCR2 hinder repair processes during resolution from injury. In line with this, human livers of ALF patients requiring liver transplantation reveal increased CD68⁺ hepatic macrophage numbers with massive infiltrates of periportal CCR2⁺ macrophages that display a proinflammatory polarization.

CONCLUSION

Infiltrating monocyte-derived macrophages aggravate APAP hepatotoxicity, and the pharmacological inhibition of either CCL2 or CCR2 might bear therapeutic potential by reducing the inflammatory reaction during the early phase of AILI. (Hepatology 2016;64:1667-1682).

RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells

RUNX1 is a transcription factor functioning both as an oncogene and a tumor suppressor in breast cancer. RUNX1 alters chromatin structure in cooperation with chromatin modifier and remodeling enzymes. In this study, we examined the relationship between RUNX1-mediated transcription and genome organization. We characterized genome-wide RUNX1 localization and performed RNA-seq and Hi-C in RUNX1-depleted and control MCF-7 breast cancer cells. RNA-seq analysis showed that RUNX1 depletion led to up-regulation of genes associated with chromatin structure and down-regulation of genes related to extracellular matrix biology, as well as NEAT1 and MALAT1 lncRNAs. Our ChIP-Seq analysis supports a prominent role for RUNX1 in transcriptional activation. About 30% of all RUNX1 binding sites were intergenic, indicating diverse roles in promoter and enhancer regulation and suggesting additional functions for RUNX1. Hi-C analysis of RUNX1-depleted cells demonstrated that overall three-dimensional genome organization is largely intact, but indicated enhanced association of RUNX1 near Topologically Associating Domain (TAD) boundaries and alterations in long-range interactions. These results suggest an architectural role for RUNX1 in fine-tuning local interactions rather than in global organization. Our results provide novel insight into RUNX1-mediated perturbations of higher-order genome organization that are functionally linked with RUNX1-dependent compromised gene expression in breast cancer cells.

Zinc finger protein 407 overexpression upregulates PPAR target gene expression and improves glucose homeostasis in mice

The peroxisome proliferator-activated receptor (PPAR) family of nuclear receptors is central to the pathophysiology and treatment of metabolic disease through the receptors' ability to regulate the expression of genes involved in glucose homeostasis, adipogenesis, and lipid metabolism. However, the mechanism by which PPAR is regulated remains incompletely understood. We generated a transgenic mouse strain (ZFP-TG) that overexpressed Zfp407 primarily in muscle and heart. Transcriptome analysis by RNA-Seq identified 1,300 differentially expressed genes in the muscle of ZFP-TG mice, among which PPAR target genes were significantly enriched. Among the physiologically important PPARγ target genes, Glucose transporter (Glut)-4 mRNA and protein levels were increased in heart and muscle. The increase in Glut4 and other transcriptional effects of Zfp407 overexpression together decreased body weight and lowered plasma glucose, insulin, and HOMA-IR scores relative to control littermates. When placed on high-fat diet, ZFP-TG mice remained more glucose tolerant than their wild-type counterparts. Cell-based assays demonstrated that Zfp407 synergistically increased the transcriptional activity of all PPAR subtypes, PPARα, PPARγ, and PPARδ. The increased PPAR activity was not associated with increased PPAR mRNA or protein levels, suggesting that Zfp407 posttranslationally regulates PPAR activity. Collectively, these results demonstrate that Zfp407 overexpression improved glucose homeostasis. Thus, Zfp407 represents a new drug target for treating metabolic disease.

The genomic landscape of schwannoma

Schwannomas are common peripheral nerve sheath tumors that can cause debilitating morbidities. We performed an integrative analysis to determine genomic aberrations common to sporadic schwannomas. Exome sequence analysis with validation by targeted DNA sequencing of 125 samples uncovered, in addition to expected NF2 disruption, recurrent mutations in ARID1A, ARID1B and DDR1. RNA sequencing identified a recurrent in-frame SH3PXD2A-HTRA1 fusion in 12/125 (10%) cases, and genomic analysis demonstrated the mechanism as resulting from a balanced 19-Mb chromosomal inversion on chromosome 10q. The fusion was associated with male gender predominance, occurring in one out of every six men with schwannoma. Methylation profiling identified distinct molecular subgroups of schwannomas that were associated with anatomical location. Expression of the SH3PXD2A-HTRA1 fusion resulted in elevated phosphorylated ERK, increased proliferation, increased invasion and in vivo tumorigenesis. Targeting of the MEK-ERK pathway was effective in fusion-positive Schwann cells, suggesting a possible therapeutic approach for this subset of tumors.

Exosomes from bulk and stem cells from human prostate cancer have a differential microRNA content that contributes cooperatively over local and pre-metastatic niche

The different prostate cancer (PCa) cell populations (bulk and cancer stem cells, CSCs) release exosomes that contain miRNAs that could modify the local or premetastatic niche. The analysis of the differential expression of miRNAs in exosomes allows evaluating the differential biological effect of both populations on the niche, and the identification of potential biomarkers and therapeutic targets. Five PCa primary cell cultures were established to originate bulk and CSCs cultures. From them, exosomes were purified by precipitation for miRNAs extraction to perform a comparative profile of miRNAs by next generation sequencing in an Illumina platform. 1839 miRNAs were identified in the exosomes. Of these 990 were known miRNAs, from which only 19 were significantly differentially expressed: 6 were overexpressed in CSCs and 13 in bulk cells exosomes. miR-100-5p and miR-21-5p were the most abundant miRNAs. Bioinformatics analysis indicated that differentially expressed miRNAs are highly related with PCa carcinogenesis, fibroblast proliferation, differentiation and migration, and angiogenesis. Besides, miRNAs from bulk cells affects osteoblast differentiation. Later, their effect was evaluated in normal prostate fibroblasts (WPMY-1) where transfection with miR-100-5p, miR-21-5p and miR-139-5p increased the expression of metalloproteinases (MMPs) -2, -9 and -13 and RANKL and fibroblast migration. The higher effect was achieved with miR21 transfection. As conclusion, miRNAs have a differential pattern between PCa bulk and CSCs exosomes that act collaboratively in PCa progression and metastasis. The most abundant miRNAs in PCa exosomes are interesting potential biomarkers and therapeutic targets.

Retrograde signalling caused by heritable mitochondrial dysfunction is partially mediated by ANAC017 and improves plant performance

Mitochondria are crucial for plant viability and are able to communicate information on their functional status to the cellular nucleus via retrograde signalling, thereby affecting gene expression. It is currently unclear if retrograde signalling in response to constitutive mitochondrial biogenesis defects is mediated by the same pathways as those triggered during acute mitochondrial dysfunction. Furthermore, it is unknown if retrograde signalling can effectively improve plant performance when mitochondrial function is constitutively impaired. Here we show that retrograde signalling in mutants defective in mitochondrial proteins RNA polymerase rpotmp or prohibitin atphb3 can be suppressed by knocking out the transcription factor ANAC017. Genome-wide RNA-seq expression analysis revealed that ANAC017 is almost solely responsible for the most dramatic transcriptional changes common to rpotmp and atphb3 mutants, regulating classical marker genes such as alternative oxidase 1a (AOX1a) and also previously-uncharacterised DUF295 genes that appear to be new retrograde markers. In contrast, ANAC017 does not regulate intra-mitochondrial gene expression or transcriptional changes unique to either rpotmp or atphb3 genotype, suggesting the existence of currently unknown signalling cascades. The data show that ANAC017 function extends beyond common retrograde transcriptional responses and affects downstream protein abundance and enzyme activity of alternative oxidase, as well as steady-state energy metabolism in atphb3 plants. Furthermore, detailed growth analysis revealed that ANAC017-dependent retrograde signalling provides benefits for growth and productivity in plants with mitochondrial defects. In conclusion, ANAC017 plays a key role in both biogenic and operational mitochondrial retrograde signalling, and improves plant performance when mitochondrial function is constitutively impaired.

Efficient generation of pancreatic β-like cells from the mouse gallbladder

Direct reprogramming is a promising approach for the replacement of β cells in diabetes. Reprogramming of cells originating from the endodermal lineage, such as acinar cells in the pancreas, liver cells and gallbladder cells has been of particular interest because of their developmental proximity to β cells. Our previous work showed that mouse gallbladder epithelium can be partially reprogrammed in vitro to generate islet-like cells (rGBC1). Here, the reprogramming protocol was substantially improved, yielding cells (rGBC2) closer to functional β cells than the 1st generation method with higher conversion efficiency and insulin expression. In addition to insulin synthesis and processing, rGBC2 presented many hallmark features of β cells, including insulin secretion in response to high glucose stimulation. Gene expression analysis indicated that rGBC2 clustered closer with β cells and had a metabolic gene expression profile resembling neonatal β cells. When transplanted into immune-deficient animals, rGBC2 were stable for at least 5months and further matured in vivo. Taken together, this approach provides further understanding of endodermal lineage conversion and potential for development of cell replacement therapy for type 1 diabetes patients.

DNA Methylation Influences the Expression of DICER-LIKE4 Isoforms, Which Encode Proteins of Alternative Localization and Function

Plant RNA silencing operates via RNA-directed DNA-methylation (RdDM) to repress transcription or by targeting mRNAs via posttranscriptional gene silencing (PTGS). These pathways rely on distinct Dicer-like (DCL) proteins that process double-stranded RNA (dsRNA) into small-interfering RNAs (siRNAs). Here, we explored the expression and subcellular localization of Arabidopsis thaliana DCL4. DCL4 expression predominates as a transcription start site isoform encoding a cytoplasmic protein, which also represents the ancestral form in plants. A longer DCL4 transcript isoform encoding a nuclear localization signal, DCL4^NLS, is present in Arabidopsis, but DNA methylation normally suppresses its expression. Hypomethylation caused by mutation, developmental reprogramming, and biotic stress correlates with enhanced DCL4^NLS expression, while hypermethylation of a DCL4 transgene causes a reduction in DCL4^NLS expression. DCL4^NLS functions in a noncanonical siRNA pathway, producing a unique set of 21-nucleotide-long "disiRNAs," for DCL4^NLS isoform-dependent siRNAs, through the nuclear RdDM dsRNA synthesis pathway. disiRNAs originate mostly from transposable elements (TEs) and TE-overlapping/proximal genes, load into the PTGS effector ARGONAUTE1 (AGO1), and display a subtle effect on transcript accumulation together with overlapping 24-nucleotide siRNAs. We propose that, via PTGS, disiRNAs could help to tighten the expression of epigenetically activated TEs and genes using the methylation-state-responsive DCL4^NLS.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.

Polycationic Nanofibers for Nucleic Acid Scavenging

Dying cells release nucleic acids (NA) and NA-containing complexes that activate inflammatory pathways of immune cells. Sustained activation of these pathways contributes to chronic inflammation frequently encountered in autoimmune and inflammatory diseases. In this study, grafting of cationic polymers onto a nanofibrous mesh enabled local scavenging of negatively charged pro-inflammatory molecules in the extracellular space. Nucleic acid scavenging nanofibers (NASFs) formed from poly(styrene-alt-maleic anhydride) conjugated with 1.8 kDa bPEI resulted in nanofibers of diameters 486 ± 9 nm. NASFs inhibited the NF-κB response stimulated by the negatively charged agonists, CpG and poly(I:C), in Ramos-blue cells but not Pam3CSK4, a nonanionic agonist. Moreover, NASFs significantly impeded NF-κB activation in cells stimulated with damage-associated molecular pattern molecules (DAMPs) released from doxorubicin killed cancer cells. In vivo application of NASFs to open wounds demonstrated nucleic acid scavenging in wounds of diabetic mice infected with Pseudomonas aeruginosa, suggesting the in vivo efficacy of NASFs. This simple technique of generating NASF results in effective localized anti-inflammation in vitro and local nucleic acid scavenging in vivo.

Genome-scale deletion screening of human long non-coding RNAs using a paired-guide RNA CRISPR-Cas9 library

CRISPR-Cas9 screens have been widely adopted to analyze coding-gene functions, but high-throughput screening of non-coding elements using this method is more challenging because indels caused by a single cut in non-coding regions are unlikely to produce a functional knockout. A high-throughput method to produce deletions of non-coding DNA is needed. We report a high-throughput genomic deletion strategy to screen for functional long non-coding RNAs (lncRNAs) that is based on a lentiviral paired-guide RNA (pgRNA) library. Applying our screening method, we identified 51 lncRNAs that can positively or negatively regulate human cancer cell growth. We validated 9 of 51 lncRNA hits using CRISPR-Cas9-mediated genomic deletion, functional rescue, CRISPR activation or inhibition and gene-expression profiling. Our high-throughput pgRNA genome deletion method will enable rapid identification of functional mammalian non-coding elements.

A mucosal imprint left by prior Escherichia coli bladder infection sensitizes to recurrent disease

Recurrent bacterial infections are a significant burden worldwide, and prior history of infection is often a significant risk factor for developing new infections. For urinary tract infection (UTI), a history of two or more episodes is an independent risk factor for acute infection. However, mechanistic knowledge of UTI pathogenesis has come almost exclusively from studies in naive mice. Here we show that, in mice, an initial Escherichia coli UTI, whether chronic or self-limiting, leaves a long-lasting molecular imprint on the bladder tissue that alters the pathophysiology of subsequent infections, affecting host susceptibility and disease outcome. In bladders of previously infected versus non-infected, antibiotic-treated mice, we found (1) an altered transcriptome and defects in cell maturation, (2) a remodelled epithelium that confers resistance to intracellular bacterial colonization, and (3) changes to cyclooxygenase-2-dependent inflammation. Furthermore, in mice with a history of chronic UTI, cyclooxygenase-2-dependent inflammation allowed a variety of clinical E. coli isolates to circumvent intracellular colonization resistance and cause severe recurrent UTI, which could be prevented by cyclooxygenase-2 inhibition or vaccination. This work provides mechanistic insight into how a history of infection can impact the risk for developing recurrent infection and has implications for the development of therapeutics for recurrent UTI.

A Sample-to-Sequence Protocol for Genus Targeted Transcriptomic Profiling: Application to Marine Synechococcus

Recent studies using whole community metagenomic and metatranscriptomic approaches are revealing important new insights into the functional potential and activity of natural marine microbial communities. Here, we complement these approaches by describing a complete ocean sample-to-sequence protocol, specifically designed to target a single bacterial genus for purposes of both DNA and RNA profiling using fluorescence activated cell sorting (FACS). The importance of defining and understanding the effects of a sampling protocol are critical if we are to gain meaningful data from environmental surveys. Rigorous pipeline trials with a cultured isolate, Synechococcus sp. BL107 demonstrate that water filtration has a well-defined but limited impact on the transcriptomic profile of this organism, whilst sample storage and multiple rounds of cell sorting have almost no effect on the resulting RNA sequence profile. Attractively, the means to replicate the sampling strategy is within the budget and expertise of most researchers.

Dietary Gluten-Induced Gut Dysbiosis Is Accompanied by Selective Upregulation of microRNAs with Intestinal Tight Junction and Bacteria-Binding Motifs in Rhesus Macaque Model of Celiac Disease

The composition of the gut microbiome reflects the overall health status of the host. In this study, stool samples representing the gut microbiomes from 6 gluten-sensitive (GS) captive juvenile rhesus macaques were compared with those from 6 healthy, age- and diet-matched peers. A total of 48 samples representing both groups were studied using V4 16S rRNA gene DNA analysis. Samples from GS macaques were further characterized based on type of diet administered: conventional monkey chow, i.e., wheat gluten-containing diet (GD), gluten-free diet (GFD), barley gluten-derived diet (BOMI) and reduced gluten barley-derived diet (RGB). It was hypothesized that the GD diet would lower the gut microbial diversity in GS macaques. This is the first report illustrating the reduction of gut microbial alpha-diversity (p < 0.05) following the consumption of dietary gluten in GS macaques. Selected bacterial families (e.g., Streptococcaceae and Lactobacillaceae) were enriched in GS macaques while Coriobacteriaceae was enriched in healthy animals. Within several weeks after the replacement of the GD by the GFD diet, the composition (beta-diversity) of gut microbiome in GS macaques started to change (p = 0.011) towards that of a normal macaque. Significance for alpha-diversity however, was not reached by the day 70 when the feeding experiment ended. Several inflammation-associated microRNAs (miR-203, -204, -23a, -23b and -29b) were upregulated (p < 0.05) in jejunum of 4 biopsied GS macaques fed GD with predicted binding sites on 16S ribosomal RNA of Lactobacillus reuteri (accession number: NR_025911), Prevotella stercorea (NR_041364) and Streptococcus luteciae (AJ297218) that were overrepresented in feces. Additionally, claudin-1, a validated tight junction protein target of miR-29b was significantly downregulated in jejunal epithelium of GS macaques. Taken together, we predict that with the introduction of effective treatments in future studies the diversity of gut microbiomes in GS macaques will approach those of healthy individuals. Further studies are needed to elucidate the regulatory pathways of inflammatory miRNAs in intestinal mucosa of GS macaques and to correlate their expression with gut dysbiosis.

RNA-Seq Transcriptomic Responses of Full-Thickness Dermal Excision Wounds to Pseudomonas aeruginosa Acute and Biofilm Infection

Pseudomonas aeruginosa infections of wounds in clinical settings are major complications whose outcomes are influenced by host responses that are not completely understood. Herein we evaluated transcriptomic changes of wounds as they counter P. aeruginosa infection—first active infection, and then chronic biofilm infection. We used the dermal full-thickness, rabbit ear excisional wound model. We studied the wound response: towards acute infection at 2, 6, and 24 hrs after inoculating 10⁶ bacteria into day-3 wounds; and, towards more chronic biofilm infection of wounds similarly infected for 24 hrs but then treated with topical antibiotic to coerce biofilm growth and evaluated at day 5 and 9 post-infection. The wounds were analyzed for bacterial counts, expression of P. aeruginosa virulence and biofilm-synthesis genes, biofilm morphology, infiltrating immune cells, re-epithelialization, and genome-wide gene expression (RNA-Seq transcriptome). This analysis revealed that 2 hrs after bacterial inoculation into day-3 wounds, the down-regulated genes (infected vs. non-infected) of the wound edge were nearly all non-coding RNAs (ncRNAs), comprised of snoRNA, miRNA, and RNU6 pseudogenes, and their down-regulation preceded a general down-regulation of skin-enriched coding gene expression. As the active infection intensified, ncRNAs remained overrepresented among down-regulated genes; however, at 6 and 24 hrs they changed to a different set, which overlapped between these times, and excluded RNU6 pseudogenes but included snRNA components of the major and minor spliceosomes. Additionally, the raw counts of multiple types of differentially-expressed ncRNAs increased on post-wounding day 3 in control wounds, but infection suppressed this increase. After 5 and 9 days, these ncRNA counts in control wounds decreased, whereas they increased in the infected, healing-impaired wounds. These data suggest a sequential and coordinated change in the levels of transcripts of multiple major classes of ncRNAs in wound cells transitioning from inflammation to the proliferation phase of healing.

Long-Term Warming Alters Carbohydrate Degradation Potential in Temperate Forest Soils

As Earth's climate warms, soil carbon pools and the microbial communities that process them may change, altering the way in which carbon is recycled in soil. In this study, we used a combination of metagenomics and bacterial cultivation to evaluate the hypothesis that experimentally raising soil temperatures by 5°C for 5, 8, or 20 years increased the potential for temperate forest soil microbial communities to degrade carbohydrates. Warming decreased the proportion of carbohydrate-degrading genes in the organic horizon derived from eukaryotes and increased the fraction of genes in the mineral soil associated with Actinobacteria in all studies. Genes associated with carbohydrate degradation increased in the organic horizon after 5 years of warming but had decreased in the organic horizon after warming the soil continuously for 20 years. However, a greater proportion of the 295 bacteria from 6 phyla (10 classes, 14 orders, and 34 families) isolated from heated plots in the 20-year experiment were able to depolymerize cellulose and xylan than bacterial isolates from control soils. Together, these findings indicate that the enrichment of bacteria capable of degrading carbohydrates could be important for accelerated carbon cycling in a warmer world.

IMPORTANCE The massive carbon stocks currently held in soils have been built up over millennia, and while numerous lines of evidence indicate that climate change will accelerate the processing of this carbon, it is unclear whether the genetic repertoire of the microbes responsible for this elevated activity will also change. In this study, we showed that bacteria isolated from plots subject to 20 years of 5°C of warming were more likely to depolymerize the plant polymers xylan and cellulose, but that carbohydrate degradation capacity is not uniformly enriched by warming treatment in the metagenomes of soil microbial communities. This study illustrates the utility of combining culture-dependent and culture-independent surveys of microbial communities to improve our understanding of the role changing microbial communities may play in soil carbon cycling under climate change.

TET-mediated DNA demethylation controls gastrulation by regulating Lefty-Nodal signalling

Mammalian genomes undergo epigenetic modifications, including cytosine methylation by DNA methyltransferases (DNMTs). Oxidation of 5-methylcytosine by the Ten-eleven translocation (TET) family of dioxygenases can lead to demethylation. Although cytosine methylation has key roles in several processes such as genomic imprinting and X-chromosome inactivation, the functional significance of cytosine methylation and demethylation in mouse embryogenesis remains to be fully determined. Here we show that inactivation of all three Tet genes in mice leads to gastrulation phenotypes, including primitive streak patterning defects in association with impaired maturation of axial mesoderm and failed specification of paraxial mesoderm, mimicking phenotypes in embryos with gain-of-function Nodal signalling. Introduction of a single mutant allele of Nodal in the Tet mutant background partially restored patterning, suggesting that hyperactive Nodal signalling contributes to the gastrulation failure of Tet mutants. Increased Nodal signalling is probably due to diminished expression of the Lefty1 and Lefty2 genes, which encode inhibitors of Nodal signalling. Moreover, reduction in Lefty gene expression is linked to elevated DNA methylation, as both Lefty-Nodal signalling and normal morphogenesis are largely restored in Tet-deficient embryos when the Dnmt3a and Dnmt3b genes are disrupted. Additionally, a point mutation in Tet that specifically abolishes the dioxygenase activity causes similar morphological and molecular abnormalities as the null mutation. Taken together, our results show that TET-mediated oxidation of 5-methylcytosine modulates Lefty-Nodal signalling by promoting demethylation in opposition to methylation by DNMT3A and DNMT3B. These findings reveal a fundamental epigenetic mechanism featuring dynamic DNA methylation and demethylation crucial to regulation of key signalling pathways in early body plan formation.

The epigenetic landscape of T cell exhaustion

Exhausted T cells in cancer and chronic viral infection express distinctive patterns of genes, including sustained expression of programmed cell death protein 1 (PD-1). However, the regulation of gene expression in exhausted T cells is poorly understood. Here, we define the accessible chromatin landscape in exhausted CD8⁺ T cells and show that it is distinct from functional memory CD8⁺ T cells. Exhausted CD8⁺ T cells in humans and a mouse model of chronic viral infection acquire a state-specific epigenetic landscape organized into functional modules of enhancers. Genome editing shows that PD-1 expression is regulated in part by an exhaustion-specific enhancer that contains essential RAR, T-bet, and Sox3 motifs. Functional enhancer maps may offer targets for genome editing that alter gene expression preferentially in exhausted CD8⁺ T cells.

A Viral Immunity Chromosome in the Marine Picoeukaryote, Ostreococcus tauri

Micro-algae of the genus Ostreococcus and related species of the order Mamiellales are globally distributed in the photic zone of world's oceans where they contribute to fixation of atmospheric carbon and production of oxygen, besides providing a primary source of nutrition in the food web. Their tiny size, simple cells, ease of culture, compact genomes and susceptibility to the most abundant large DNA viruses in the sea render them attractive as models for integrative marine biology. In culture, spontaneous resistance to viruses occurs frequently. Here, we show that virus-producing resistant cell lines arise in many independent cell lines during lytic infections, but over two years, more and more of these lines stop producing viruses. We observed sweeping over-expression of all genes in more than half of chromosome 19 in resistant lines, and karyotypic analyses showed physical rearrangements of this chromosome. Chromosome 19 has an unusual genetic structure whose equivalent is found in all of the sequenced genomes in this ecologically important group of green algae.

We propose that chromosome 19 of O. tauri is specialized in defence against viral attack, a constant threat for all planktonic life, and that the most likely cause of resistance is the over-expression of numerous predicted glycosyltransferase genes. O. tauri thus provides an amenable model for molecular analysis of genome evolution under environmental stress and for investigating glycan-mediated host-virus interactions, such as those seen in herpes, influenza, HIV, PBCV and mimivirus.

Transcriptomic responses of a simplified soil microcosm to a plant pathogen and its biocontrol agent reveal a complex reaction to harsh habitat

BACKGROUND

Soil microorganisms are key determinants of soil fertility and plant health. Soil phytopathogenic fungi are one of the most important causes of crop losses worldwide. Microbial biocontrol agents have been extensively studied as alternatives for controlling phytopathogenic soil microorganisms, but molecular interactions between them have mainly been characterised in dual cultures, without taking into account the soil microbial community. We used an RNA sequencing approach to elucidate the molecular interplay of a soil microbial community in response to a plant pathogen and its biocontrol agent, in order to examine the molecular patterns activated by the microorganisms.

RESULTS

A simplified soil microcosm containing 11 soil microorganisms was incubated with a plant root pathogen (Armillaria mellea) and its biocontrol agent (Trichoderma atroviride) for 24 h under controlled conditions. More than 46 million paired-end reads were obtained for each replicate and 28,309 differentially expressed genes were identified in total. Pathway analysis revealed complex adaptations of soil microorganisms to the harsh conditions of the soil matrix and to reciprocal microbial competition/cooperation relationships. Both the phytopathogen and its biocontrol agent were specifically recognised by the simplified soil microcosm: defence reaction mechanisms and neutral adaptation processes were activated in response to competitive (T. atroviride) or non-competitive (A. mellea) microorganisms, respectively. Moreover, activation of resistance mechanisms dominated in the simplified soil microcosm in the presence of both A. mellea and T. atroviride. Biocontrol processes of T. atroviride were already activated during incubation in the simplified soil microcosm, possibly to occupy niches in a competitive ecosystem, and they were not further enhanced by the introduction of A. mellea.

CONCLUSIONS

This work represents an additional step towards understanding molecular interactions between plant pathogens and biocontrol agents within a soil ecosystem. Global transcriptional analysis of the simplified soil microcosm revealed complex metabolic adaptation in the soil environment and specific responses to antagonistic or neutral intruders.

Epithelial Coculture and l-Lactate Promote Growth of Helicobacter cinaedi under H2-Free Aerobic Conditions

Helicobacter cinaedi is an emerging opportunistic pathogen associated with infections of diverse anatomic sites. Nevertheless, the species demonstrates fastidious axenic growth; it has been described as requiring a microaerobic atmosphere, along with a strong preference for supplemental H₂ gas. In this context, we examined the hypothesis that in vitro growth of H. cinaedi could be enhanced by coculture with human epithelial cells. When inoculated (in Ham's F12 medium) over Caco-2 monolayers, the type strain (ATCC BAA-847) gained the ability to proliferate under H₂-free aerobic conditions. Identical results were observed during coculture with several other monolayer types (LS-174T, AGS, and HeLa). Under chemically defined conditions, 40 amino acids and carboxylates were screened for their effect on the organism's atmospheric requirements. Several molecules promoted H₂-free aerobic proliferation, although it occurred most prominently with millimolar concentrations of l-lactate. The growth response of H. cinaedi to Caco-2 cells and l-lactate was confirmed with a collection of 12 human-derived clinical strains. mRNA sequencing was next performed on the type strain under various growth conditions. In addition to providing a whole-transcriptome profile of H. cinaedi, this analysis demonstrated strong constitutive expression of the l-lactate utilization locus, as well as differential transcription of terminal respiratory proteins as a function of Caco-2 coculture and l-lactate supplementation. Overall, these findings challenge traditional views of H. cinaedi as an obligate microaerophile.

IMPORTANCE

H. cinaedi is an increasingly recognized pathogen in people with compromised immune systems. Atypical among other members of its bacterial class, H. cinaedi has been associated with infections of diverse anatomic sites. Growing H. cineadi in the laboratory is quite difficult, due in large part to the need for a specialized atmosphere. The suboptimal growth of H. cinaedi is an obstacle to clinical diagnosis, and it also limits investigation into the organism's biology. The current work shows that H. cinaedi has more flexible atmospheric requirements in the presence of host cells and a common host-derived molecule. This nutritional interplay raises new questions about how the organism behaves during human infections and provides insights for how to optimize its laboratory cultivation.

RNA-Seq Count Data Modelling by Grey Relational Analysis and Nonparametric Gaussian Process

This paper introduces an approach to classification of RNA-seq read counts using grey relational analysis (GRA) and Bayesian Gaussian process (GP) models. Read counts are transformed to microarray-like data to facilitate normal-based statistical methods. GRA is designed to select differentially expressed genes by integrating outcomes of five individual feature selection methods including two-sample t-test, entropy test, Bhattacharyya distance, Wilcoxon test and receiver operating characteristic curve. GRA performs as an aggregate filter method through combining advantages of the individual methods to produce significant feature subsets that are then fed into a nonparametric GP model for classification. The proposed approach is verified by using two benchmark real datasets and the five-fold cross-validation method. Experimental results show the performance dominance of the GRA-based feature selection method as well as GP classifier against their competing methods. Moreover, the results demonstrate that GRA-GP considerably dominates the sparse Poisson linear discriminant analysis classifiers, which were introduced specifically for read counts, on different number of features. The proposed approach therefore can be implemented effectively in real practice for read count data analysis, which is useful in many applications including understanding disease pathogenesis, diagnosis and treatment monitoring at the molecular level.

Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia

Sexual dimorphism is an interesting biological phenomenon. Previous studies showed that DNA methylation might play a role in sexual dimorphism. However, the overall picture of the genome-wide methylation landscape in sexually dimorphic species remains unclear. We analyzed the DNA methylation landscape and transcriptome in hybrid tilapia (Oreochromis spp.) using whole genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq). We found 4,757 sexually dimorphic differentially methylated regions (DMRs), with significant clusters of DMRs located on chromosomal regions associated with sex determination. CpG methylation in promoter regions was negatively correlated with the gene expression level. MAPK/ERK pathway was upregulated in male tilapia. We also inferred active cis-regulatory regions (ACRs) in skeletal muscle tissues from WGBS datasets, revealing sexually dimorphic cis-regulatory regions. These results suggest that DNA methylation contribute to sex-specific phenotypes and serve as resources for further investigation to analyze the functions of these regions and their contributions towards sexual dimorphisms.

Analyses of RNA-Seq and sRNA-Seq data reveal a complex network of anti-viral defense in TCV-infected Arabidopsis thaliana

In order to identify specific plant anti-viral genes related to the miRNA regulatory pathway, RNA-Seq and sRNA-Seq were performed using Arabidopsis WT and dcl1-9 mutant line. A total of 5,204 DEGs were identified in TCV-infected WT plants. In contrast, only 595 DEGs were obtained in the infected dcl1-9 mutant plants. GO enrichment analysis of the shared DEGs and dcl1-9 unique DEGs showed that a wide range of biological processes were affected in the infected WT plants. In addition, miRNAs displayed different patterns between mock and infected WT plants. This is the first global view of dcl1-9 transcriptome which provides TCV responsive miRNAs data. In conclusion, our results indicated the significance of DCL1 and suggested that PPR genes may play an important role in plant anti-viral defense.

A two-step integrated approach to detect differentially expressed genes in RNA-Seq data

One of the primary objectives of ribonucleic acid (RNA) sequencing or RNA-Seq experiment is to identify differentially expressed (DE) genes in two or more treatment conditions. It is a common practice to assume that all read counts from RNA-Seq data follow overdispersed (OD) Poisson or negative binomial (NB) distribution, which is sometimes misleading because within each condition, some genes may have unvarying transcription levels with no overdispersion. In such a case, it is more appropriate and logical to consider two sets of genes: OD and non-overdispersed (NOD). We propose a new two-step integrated approach to distinguish DE genes in RNA-Seq data using standard Poisson and NB models for NOD and OD genes, respectively. This is an integrated approach because this method can be merged with any other NB-based methods for detecting DE genes. We design a simulation study and analyze two real RNA-Seq data to evaluate the proposed strategy. We compare the performance of this new method combined with the three [Formula: see text]-software packages namely edgeR, DESeq2, and DSS with their default settings. For both the simulated and real data sets, integrated approaches perform better or at least equally well compared to the regular methods embedded in these [Formula: see text]-packages.

Pan-cancer transcriptomic analysis associates long non-coding RNAs with key mutational driver events

Thousands of long non-coding RNAs (lncRNAs) lie interspersed with coding genes across the genome, and a small subset has been implicated as downstream effectors in oncogenic pathways. Here we make use of transcriptome and exome sequencing data from thousands of tumours across 19 cancer types, to identify lncRNAs that are induced or repressed in relation to somatic mutations in key oncogenic driver genes. Our screen confirms known coding and non-coding effectors and also associates many new lncRNAs to relevant pathways. The associations are often highly reproducible across cancer types, and while many lncRNAs are co-expressed with their protein-coding hosts or neighbours, some are intergenic and independent. We highlight lncRNAs with possible functions downstream of the tumour suppressor TP53 and the master antioxidant transcription factor NFE2L2. Our study provides a comprehensive overview of lncRNA transcriptional alterations in relation to key driver mutational events in human cancers.

Long non-coding RNAs are implicated in multiple aspects of tumourigenesis. Here, the authors generate a landscape of these macromolecules in a wide array of cancer types and examine which RNAs are transcriptionally altered in relation to somatic driver mutations in established coding cancer genes.

Detection of generic differential RNA processing events from RNA-seq data

RNA-seq data analysis has revealed abundant alternative splicing in eukaryotic mRNAs. However, splicing is only one of many processing events that transcripts may undergo during their lifetime. We present here RNAprof (RNA profile analysis), a program for the detection of differential processing events from the comparison of RNA-seq experiments. RNAprof implements a specific gene-level normalization procedure and compares RNA-seq coverage profiles at nucleotide resolution to detect regions of significant coverage differences, independently of splice sites or other gene features. We used RNAprof to analyze the effect of alternative-splicing regulators NSRa and NSRb on the Arabidopsis thaliana transcriptome. A number of intron retention events and alternative transcript structures were specifically detected by RNAprof and confirmed by qRT-PCR. Further tests using a public Mus musculus RNA-seq dataset and comparisons with other RNA isoform predictors showed that RNAprof uniquely identified sets of highly significant processing events as well as other relevant library-specific differences in RNA-seq profiles. This highlights an important layer of variation that remains undetected by current protocols for RNA-seq analysis.

Relationships between PROMPT and gene expression

Most mammalian protein-coding gene promoters are divergent, yielding promoter upstream transcripts (PROMPTs) in the reverse direction from their conventionally produced mRNAs. PROMPTs are rapidly degraded by the RNA exosome rendering a general function of these molecules elusive. Yet, levels of certain PROMPTs are altered in stress conditions, like the DNA damage response (DDR), suggesting a possible regulatory role for at least a subset of these molecules. Here we manipulate PROMPT levels by either exosome depletion or UV treatment and analyze possible effects on their neighboring genes. For the CTSZ and DAP genes we find that TFIIB and TBP promoter binding decrease when PROMPTs accumulate. Moreover, DNA methylation increases concomitant with the recruitment of the DNA methyltransferase DNMT3B. Thus, although a correlation between increased PROMPT levels and decreased gene activity is generally absent, some promoters may have co-opted their divergent transcript production for regulatory purposes.

Potent and selective bivalent inhibitors of BET bromodomains

Proteins of the bromodomain and extraterminal (BET) family, in particular bromodomain-containing protein 4 (BRD4), are of great interest as biological targets. BET proteins contain two separate bromodomains, and existing inhibitors bind to them monovalently. Here we describe the discovery and characterization of probe compound biBET, capable of engaging both bromodomains simultaneously in a bivalent, in cis binding mode. The evidence provided here was obtained in a variety of biophysical and cellular experiments. The bivalent binding results in very high cellular potency for BRD4 binding and pharmacological responses such as disruption of BRD4-mediator complex subunit 1 foci with an EC₅₀ of 100 pM. These compounds will be of considerable utility as BET/BRD4 chemical probes. This work illustrates a novel concept in ligand design-simultaneous targeting of two separate domains with a drug-like small molecule-providing precedent for a potentially more effective paradigm for developing ligands for other multi-domain proteins.

Sabotaging of the oxidative stress response by an oncogenic noncoding RNA

Overexpression of the multiple myeloma set domain (MMSET) Wolf-Hirschhorn syndrome candidate 1 gene, which contains an orphan box H/ACA class small nucleolar RNA, ACA11, in an intron, is associated with several cancer types, including multiple myeloma (MM). ACA11 and MMSET are overexpressed cotranscriptionally as a result of the t(4;14) chromosomal translocation in a subset of patients with MM. RNA sequencing of CD138⁺ tumor cells from t(4;14)-positive and -negative MM patient bone marrow samples revealed an enhanced oxidative phosphorylation mRNA signature. Supporting these data, ACA11 overexpression in a t(4;14)-negative MM cell line, MM1.S, demonstrated enhanced reactive oxygen species (ROS) levels. In addition, an enhancement of cell proliferation, increased soft agar colony size, and elevated ERK1/2 phosphorylation were observed. This ACA11-driven hyperproliferative phenotype depended on increased ROS levels as exogenously added antioxidants attenuate the increased proliferation. A major transcriptional regulator of the cellular antioxidant response, nuclear factor (erythroid-derived 2)-like 2 (NRF2), shuttled to the nucleus, as expected, in response to ACA11-driven increases in ROS; however, transcriptional up-regulation of some of NRF2's antioxidant target genes was abrogated in the presence of ACA11 overexpression. These data show for the first time that ACA11 promotes proliferation through inhibition of NRF2 function resulting in sustained ROS levels driving cancer cell proliferation.-Mahajan, N., Wu, H.-J., Bennett, R. L., Troche, C., Licht, J. D., Weber, J. D., Maggi, L. B., Jr., Tomasson, M. H. Sabotaging of the oxidative stress response by an oncogenic noncoding RNA.

GPR133 (ADGRD1), an adhesion G-protein-coupled receptor, is necessary for glioblastoma growth

Glioblastoma (GBM) is a deadly primary brain malignancy with extensive intratumoral hypoxia. Hypoxic regions of GBM contain stem-like cells and are associated with tumor growth and angiogenesis. The molecular mechanisms that regulate tumor growth in hypoxic conditions are incompletely understood. Here, we use primary human tumor biospecimens and cultures to identify GPR133 (ADGRD1), an orphan member of the adhesion family of G-protein-coupled receptors, as a critical regulator of the response to hypoxia and tumor growth in GBM. GPR133 is selectively expressed in CD133+ GBM stem cells (GSCs) and within the hypoxic areas of PPN in human biospecimens. GPR133 mRNA is transcriptionally upregulated by hypoxia in hypoxia-inducible factor 1α (Hif1α)-dependent manner. Genetic inhibition of GPR133 with short hairpin RNA reduces the prevalence of CD133+ GSCs, tumor cell proliferation and tumorsphere formation in vitro. Forskolin rescues the GPR133 knockdown phenotype, suggesting that GPR133 signaling is mediated by cAMP. Implantation of GBM cells with short hairpin RNA-mediated knockdown of GPR133 in the mouse brain markedly reduces tumor xenograft formation and increases host survival. Analysis of the TCGA data shows that GPR133 expression levels are inversely correlated with patient survival. These findings indicate that GPR133 is an important mediator of the hypoxic response in GBM and has significant protumorigenic functions. We propose that GPR133 represents a novel molecular target in GBM and possibly other malignancies where hypoxia is fundamental to pathogenesis.

Short- and long-term effects of 56Fe irradiation on cognition and hippocampal DNA methylation and gene expression

Background

Astronauts are exposed to ⁵⁶Fe ions that may pose a significant health hazard during and following prolonged missions in deep space. We showed previously that object recognition requiring the hippocampus, a structure critical for cognitive function, is affected in 2-month-old mice irradiated with ⁵⁶Fe ions. Here we examined object recognition in 6-month-old mice irradiated with ⁵⁶Fe ions, a biological age more relevant to the typical ages of astronauts. Moreover, because the mechanisms mediating the detrimental effects of ⁵⁶Fe ions on hippocampal function are unclear, we examined changes in hippocampal networks involved in synaptic plasticity and memory, gene expression, and epigenetic changes in cytosine methylation (5mC) and hydroxymethylation (5hmC) that could accompany changes in gene expression. We assessed the effects of whole body ⁵⁶Fe ion irradiation at early (2 weeks) and late (20 weeks) time points on hippocampus-dependent memory and hippocampal network stability, and whether these effects are associated with epigenetic changes in hippocampal DNA methylation (both 5mC and 5hmC) and gene expression.

Results

At the two-week time point, object recognition and network stability were impaired following irradiation at the 0.1 and 0.4 Gy dose, but not following irradiation at the 0.2 Gy dose. No impairments in object recognition or network stability were seen at the 20-week time point at any irradiation dose used. Consistent with this pattern, the significance of pathways for gene categories for 5hmC was lower, though not eliminated, at the 20-week time point compared to the 2-week time point. Similarly, significant changes were observed for 5mC gene pathways at the 2-week time point, but no significant gene categories were observed at the 20-week time point. Only the 5hmC changes tracked with gene expression changes.

Conclusions

Dose- and time-dependent epigenomic remodeling in the hippocampus following ⁵⁶Fe ion exposure correlates with behavioral changes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3110-7) contains supplementary material, which is available to authorized users.