Genetic effects on molecular network states explain complex traits.

The complexity of many cellular and organismal traits results from the integration of genetic and environmental factors via molecular networks. Network structure and effect propagation are best understood at the level of functional modules, but so far, no concept has been established to include the global network state. Here, we show when and how genetic perturbations lead to molecular changes that are confined to small parts of a network versus when they lead to modulation of network states. Integrating multi-omics profiling of genetically heterogeneous budding and fission yeast strains with an array of cellular traits identified a central state transition of the yeast molecular network that is related to PKA and TOR (PT) signaling. Genetic variants affecting this PT state globally shifted the molecular network along a single-dimensional axis, thereby modulating processes including energy and amino acid metabolism, transcription, translation, cell cycle control, and cellular stress response. We propose that genetic effects can propagate through large parts of molecular networks because of the functional requirement to centrally coordinate the activity of fundamental cellular processes.

Does rapid sequence divergence preclude RNA structure conservation in vertebrates?

Accelerated evolution of any portion of the genome is of significant interest, potentially signaling positive selection of phenotypic traits and adaptation. Accelerated evolution remains understudied for structured RNAs, despite the fact that an RNA's structure is often key to its function. RNA structures are typically characterized by compensatory (structure-preserving) basepair changes that are unexpected given the underlying sequence variation, i.e., they have evolved through negative selection on structure. We address the question of how fast the primary sequence of an RNA can change through evolution while conserving its structure. Specifically, we consider predicted and known structures in vertebrate genomes. After careful control of false discovery rates, we obtain 13 de novo structures (and three known Rfam structures) that we predict to have rapidly evolving sequences-defined as structures where the primary sequences of human and mouse have diverged at least twice as fast (1.5 times for Rfam) as nearby neutrally evolving sequences. Two of the three known structures function in translation inhibition related to infection and immune response. We conclude that rapid sequence divergence does not preclude RNA structure conservation in vertebrates, although these events are relatively rare.

Systematic inference of indirect transcriptional regulation by protein kinases and phosphatases.

Gene expression is controlled by pathways of regulatory factors often involving the activity of protein kinases on transcription factor proteins. Despite this well established mechanism, the number of well described pathways that include the regulatory role of protein kinases on transcription factors is surprisingly scarce in eukaryotes. To address this, PhosTF was developed to infer functional regulatory interactions and pathways in both simulated and real biological networks, based on linear cyclic causal models with latent variables. GeneNetWeaverPhos, an extension of GeneNetWeaver, was developed to allow the simulation of perturbations in known networks that included the activity of protein kinases and phosphatases on gene regulation. Over 2000 genome-wide gene expression profiles, where the loss or gain of regulatory genes could be observed to perturb gene regulation, were then used to infer the existence of regulatory interactions, and their mode of regulation in the budding yeast Saccharomyces cerevisiae. Despite the additional complexity, our inference performed comparably to the best methods that inferred transcription factor regulation assessed in the DREAM4 challenge on similar simulated networks. Inference on integrated genome-scale data sets for yeast identified ∼ 8800 protein kinase/phosphatase-transcription factor interactions and ∼ 6500 interactions among protein kinases and/or phosphatases. Both types of regulatory predictions captured statistically significant numbers of known interactions of their type. Surprisingly, kinases and phosphatases regulated transcription factors by a negative mode or regulation (deactivation) in over 70% of the predictions.

Emergence of Phenotypically Distinct Subpopulations Is a Factor in Adaptation of Recombinant Saccharomyces cerevisiae under Glucose-Limited Conditions.

Cells cultured in a nutrient-limited environment can undergo adaptation, which confers improved fitness under long-term energy limitation. We have shown previously how a recombinant Saccharomyces cerevisiae strain, producing a heterologous insulin product, under glucose-limited conditions adapts over time at the average population level. Here, we investigated this adaptation at the single-cell level by application of fluorescence-activated cell sorting (FACS) and showed that the following three apparent phenotypes underlie the adaptive response observed at the bulk level: (i) cells that drastically reduced insulin production (23%), (ii) cells with reduced enzymatic capacity in central carbon metabolism (46%), and (iii) cells that exhibited pseudohyphal growth (31%). We speculate that the phenotypic heterogeneity is a result of different mechanisms to increase fitness. Cells with reduced insulin productivity have increased fitness by reducing the burden of the heterologous insulin production, and the populations with reduced enzymatic capacity of the central carbon metabolism and pseudohyphal growth have increased fitness toward the glucose-limited conditions. The results highlight the importance of considering population heterogeneity when studying adaptation and evolution. IMPORTANCE The yeast Saccharomyces cerevisiae is an attractive microbial host for industrial production and is used widely for manufacturing, e.g., pharmaceuticals. Chemostat cultivation mode is an efficient cultivation strategy for industrial production processes as it ensures a constant, well-controlled cultivation environment. Nevertheless, both the production of a heterologous product and the constant cultivation environment in the chemostat impose a selective pressure on the production organism, which may result in adaptation and loss of productivity. The exact mechanisms behind the observed adaptation and loss of performance are often unidentified. We used a recombinant S. cerevisiae strain producing heterologous insulin and investigated the adaptation occurring during chemostat growth at the single-cell level. We showed that three apparent phenotypes underlie the adaptive response observed at the bulk level in the chemostat. These findings highlight the importance of considering population heterogeneity when studying adaptation in industrial bioprocesses.

Innate IL-23/Type 17 immune responses mediate the effect of the 17q21 locus on childhood asthma.

BACKGROUND: Several childhood asthma risk loci that relate to immune function have been identified by genome-wide association studies (GWAS), but the underlying mechanisms remain unknown. OBJECTIVE: Here, we examined whether perturbed innate immune responses mediate the association between known genetic risk variants and development of childhood asthma. METHODS: Peripheral blood mononuclear cells from 336 six-month-old infants from the Copenhagen Prospective Studies on Asthma in Childhood (COPSAC2000 ) cohort were stimulated in vitro with six different innate ligands (LPS, CpG, poly(I:C), R848, HDMAPP and aluminium hydroxide together with low levels of LPS) followed by quantification of 18 released cytokines and chemokines 40 h after the stimulations. The innate immune response profiles were decomposed by principal component (PC) analysis, and PC1-5 were used in mediation analyses of the effect of 25 known genetic risk variants on childhood asthma until age 7. RESULTS: The effects of two variants from the 17q21 locus (rs7216389, rs2305480) on asthma and exacerbation risk were significantly mediated by immune parameters induced in response to ligands mimicking intracellular colonization; bacterial DNA (CpG) and double-stranded viral RNA (poly(I:C)). The Th17 and innate lymphoid cell type 3-amplifying cytokine IL-23 was the most prominent cytokine involved. CONCLUSION: The 17q21 effect on childhood asthma and exacerbations was partly mediated by deregulation of IL-23 in response to intracellular microbial ligands, which may suggest ineffective clearance of intracellular pathogens in the lungs.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response.

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

The identification and functional annotation of RNA structures conserved in vertebrates.

Structured elements of RNA molecules are essential in, e.g., RNA stabilization, localization, and protein interaction, and their conservation across species suggests a common functional role. We computationally screened vertebrate genomes for conserved RNA structures (CRSs), leveraging structure-based, rather than sequence-based, alignments. After careful correction for sequence identity and GC content, we predict ∼516,000 human genomic regions containing CRSs. We find that a substantial fraction of human-mouse CRS regions (1) colocalize consistently with binding sites of the same RNA binding proteins (RBPs) or (2) are transcribed in corresponding tissues. Additionally, a CaptureSeq experiment revealed expression of many of our CRS regions in human fetal brain, including 662 novel ones. For selected human and mouse candidate pairs, qRT-PCR and in vitro RNA structure probing supported both shared expression and shared structure despite low abundance and low sequence identity. About 30,000 CRS regions are located near coding or long noncoding RNA genes or within enhancers. Structured (CRS overlapping) enhancer RNAs and extended 3' ends have significantly increased expression levels over their nonstructured counterparts. Our findings of transcribed uncharacterized regulatory regions that contain CRSs support their RNA-mediated functionality.

A scored human protein-protein interaction network to catalyze genomic interpretation.

Genome-scale human protein-protein interaction networks are critical to understanding cell biology and interpreting genomic data, but challenging to produce experimentally. Through data integration and quality control, we provide a scored human protein-protein interaction network (InWeb_InBioMap, or InWeb_IM) with severalfold more interactions (>500,000) and better functional biological relevance than comparable resources. We illustrate that InWeb_InBioMap enables functional interpretation of >4,700 cancer genomes and genes involved in autism.

Characterization of glutathione proteome in CHO cells and its relationship with productivity and cholesterol synthesis.

Glutathione (GSH) plays a central role in the redox balance maintenance in mammalian cells. Previous studies of industrial Chinese hamster ovary cell lines have demonstrated a relationship between GSH metabolism and clone productivity. However, a thorough investigation is required to understand this relationship and potentially highlight new targets for cell engineering. In this study, we have modulated the GSH intracellular content of an industrial cell line under bioprocess conditions to further elucidate the role of the GSH synthesis pathway. Two strategies were used: the variation of cystine supply and the direct inhibition of the GSH synthesis using buthionine sulfoximine (BSO). Over time of the bioprocess, a correlation between intracellular GSH and product titer has been observed. Analysis of metabolites uptake/secretion rates and proteome comparison between BSO-treated cells and nontreated cells has highlighted a slowdown of the tricarboxylic acid cycle leading to a secretion of lactate and alanine in the extracellular environment. Moreover, an adaptation of the GSH-related proteome has been observed with an upregulation of the regulatory subunit of glutamate-cysteine ligase and a downregulation of a specific GSH transferase subgroup, the Mu family. Surprisingly, the main impact of BSO treatment was observed on a global downregulation of the cholesterol synthesis pathways. As cholesterol is required for protein secretion, it could be the missing piece of the puzzle to finally elucidate the link between GSH synthesis and productivity.

Fluctuations in glucose availability prevent global proteome changes and physiological transition during prolonged chemostat cultivations of Saccharomyces cerevisiae.

Chemostat cultivation mode imposes selective pressure on the cells, which may result in slow adaptation in the physiological state over time. We applied a two-compartment scale-down chemostat system imposing feast-famine conditions to characterize the long-term (100 s of hours) response of Saccharomyces cerevisiae to fluctuating glucose availability. A wild-type strain and a recombinant strain, expressing an insulin precursor, were cultured in the scale-down system, and analyzed at the physiological and proteomic level. Phenotypes of both strains were compared with those observed in a well-mixed chemostat. Our results show that S. cerevisiae subjected to long-term chemostat conditions undergoes a global reproducible shift in its cellular state and that this transition occurs faster and is larger in magnitude for the recombinant strain including a significant decrease in the expression of the insulin product. We find that the transition can be completely avoided in the presence of fluctuations in glucose availability as the strains subjected to feast-famine conditions under otherwise constant culture conditions exhibited constant levels of the measured proteome for over 250 hr. We hypothesize possible mechanisms responsible for the observed phenotypes and suggest experiments that could be used to test these mechanisms.

Familial co-occurrence of congenital heart defects follows distinct patterns.

Aims: Congenital heart defects (CHD) affect almost 1% of all live born children and the number of adults with CHD is increasing. In families where CHD has occurred previously, estimates of recurrence risk, and the type of recurring malformation are important for counselling and clinical decision-making, but the recurrence patterns in families are poorly understood. We aimed to determine recurrence patterns, by investigating the co-occurrences of CHD in 1163 families with known malformations, comprising 3080 individuals with clinically confirmed diagnosis. Methods and results: We calculated rates of concordance and discordance for 41 specific types of malformations, observing a high variability in the rates of concordance and discordance. By calculating odds ratios for each of 1640 pairs of discordant lesions observed between affected family members, we were able to identify 178 pairs of malformations that co-occurred significantly more or less often than expected in families. The data show that distinct groups of cardiac malformations co-occur in families, suggesting influence from underlying developmental mechanisms. Analysis of human and mouse susceptibility genes showed that they were shared in 19% and 20% of pairs of co-occurring discordant malformations, respectively, but none of malformations that rarely co-occur, suggesting that a significant proportion of co-occurring lesions in families is caused by overlapping susceptibility genes. Conclusion: Familial CHD follow specific patterns of recurrence, suggesting a strong influence from genetically regulated developmental mechanisms. Co-occurrence of malformations in families is caused by shared susceptibility genes.

Correction to: Industrial antifoam agents impair ethanol fermentation and induce stress responses in yeast cells.

The article "Industrial antifoam agents impair ethanol fermentation and induce stress responses in yeast cells" was originally published Online First without open access. After publication in volume 101, issue 22, page 8237-8248, the author decided to opt for Open Choice and to make the article an open access publication.

Major differences between human atopic dermatitis and murine models, as determined by using global transcriptomic profiling.

BACKGROUND: Atopic dermatitis (AD) is caused by a complex interplay between immune and barrier abnormalities. Murine models of AD are essential for preclinical assessments of new treatments. Although many models have been used to simulate AD, their transcriptomic profiles are not fully understood, and a comparison of these models with the human AD transcriptomic fingerprint is lacking. OBJECTIVE: We sought to evaluate the transcriptomic profiles of 6 common murine models and determine how they relate to human AD skin. METHODS: Transcriptomic profiling was performed by using microarrays and quantitative RT-PCR on biopsy specimens from NC/Nga, flaky tail, Flg-mutated, ovalbumin-challenged, oxazolone-challenged, and IL-23-injected mice. Gene expression data of patients with AD, psoriasis, and contact dermatitis were obtained from previous patient cohorts. Criteria of a fold change of 2 or greater and a false discovery rate of 0.05 or less were used for gene arrays. RESULTS: IL-23-injected, NC/Nga, and oxazolone-challenged mice show the largest homology with our human meta-analysis-derived AD transcriptome (37%, 18%, 17%, respectively). Similar to human AD, robust TH1, TH2, and also TH17 activation are seen in IL-23-injected and NC/Nga mice, with similar but weaker inflammation in ovalbumin-challenged mice. Oxazolone-challenged mice show a TH1-centered reaction, and flaky tail mice demonstrate a strong TH17 polarization. Flg-mutated mice display filaggrin downregulation without significant inflammation. CONCLUSION: No single murine model fully captures all aspects of the AD profile; instead, each model reflects different immune or barrier disease aspects. Overall, among the 6 murine models, IL-23-injected mice best simulate human AD; still, the translational focus of the investigation should determine which model is most applicable.

miR-18b overexpression identifies mantle cell lymphoma patients with poor outcome and improves the MIPI-B prognosticator.

Recent studies show that mantle cell lymphoma (MCL) express aberrant microRNA (miRNA) profiles; however, the clinical effect of miRNA expression has not previously been examined and validated in large prospective homogenously treated cohorts. We performed genome-wide miRNA microarray profiling of 74 diagnostic MCL samples from the Nordic MCL2 trial (screening cohort). Prognostic miRNAs were validated in diagnostic MCL samples from 94 patients of the independent Nordic MCL3 trial (validation cohort). Three miRNAs (miR-18b, miR-92a, and miR-378d) were significantly differentially expressed in patients who died of MCL in both cohorts. MiR-18b was superior to miR-92a and miR-378d in predicting high risk. Thus, we generated a new biological MCL International Prognostic Index (MIPI-B)-miR prognosticator, combining expression levels of miR-18b with MIPI-B data. Compared to the MIPI-B, this prognosticator improved identification of high-risk patients with regard to cause-specific, overall, and progression-free survival. Transfection of 2 MCL cell lines with miR-18b decreased their proliferation rate without inducing apoptosis, suggesting that miR-18b may render MCL cells resistant to chemotherapy by decelerating cell proliferation. We conclude that overexpression of miR-18b identifies patients with poor prognosis in 2 large prospective MCL cohorts and adds prognostic information to the MIPI-B. MiR-18b may reduce the proliferation rate of MCL cells as a mechanism of chemoresistance.

Industrial antifoam agents impair ethanol fermentation and induce stress responses in yeast cells.

The Brazilian sugarcane industry constitutes one of the biggest and most efficient ethanol production processes in the world. Brazilian ethanol production utilizes a unique process, which includes cell recycling, acid wash, and non-aseptic conditions. Process characteristics, such as extensive CO2 generation, poor quality of raw materials, and frequent contaminations, all lead to excessive foam formation during fermentations, which is treated with antifoam agents (AFA). In this study, we have investigated the impact of industrial AFA treatments on the physiology and transcriptome of the industrial ethanol strain Saccharomyces cerevisiae CAT-1. The investigated AFA included industrially used AFA acquired from Brazilian ethanol plants and commercially available AFA commonly used in the fermentation literature. In batch fermentations, it was shown that industrial AFA compromised growth rates and glucose uptake rates, while commercial AFA had no effect in concentrations relevant for defoaming purposes. Industrial AFA were further tested in laboratory scale simulations of the Brazilian ethanol production process and proved to decrease cell viability compared to the control, and the effects were intensified with increasing AFA concentrations and exposure time. Transcriptome analysis showed that AFA treatments induced additional stress responses in yeast cells compared to the control, shown by an up-regulation of stress-specific genes and a down-regulation of lipid biosynthesis, especially ergosterol. By documenting the detrimental effects associated with chemical AFA, we highlight the importance of developing innocuous systems for foam control in industrial fermentation processes.

Lysine deacetylase inhibition prevents diabetes by chromatin-independent immunoregulation and ß-cell protection.

Type 1 diabetes is due to destruction of pancreatic ß-cells. Lysine deacetylase inhibitors (KDACi) protect ß-cells from inflammatory destruction in vitro and are promising immunomodulators. Here we demonstrate that the clinically well-tolerated KDACi vorinostat and givinostat revert diabetes in the nonobese diabetic (NOD) mouse model of type 1 diabetes and counteract inflammatory target cell damage by a mechanism of action consistent with transcription factor--rather than global chromatin--hyperacetylation. Weaning NOD mice received low doses of vorinostat and givinostat in their drinking water until 100-120 d of age. Diabetes incidence was reduced by 38% and 45%, respectively, there was a 15% increase in the percentage of islets without infiltration, and pancreatic insulin content increased by 200%. Vorinostat treatment increased the frequency of functional regulatory T-cell subsets and their transcription factors Gata3 and FoxP3 in parallel to a decrease in inflammatory dendritic cell subsets and their cytokines IL-6, IL-12, and TNF-α. KDACi also inhibited LPS-induced Cox-2 expression in peritoneal macrophages from C57BL/6 and NOD mice. In insulin-producing ß-cells, givinostat did not upregulate expression of the anti-inflammatory genes Socs1-3 or sirtuin-1 but reduced levels of IL-1ß + IFN-γ-induced proinflammatory Il1a, Il1b, Tnfα, Fas, Cxcl2, and reduced cytokine-induced ERK phosphorylation. Further, NF-κB genomic iNos promoter binding was reduced by 50%, and NF-κB-dependent mRNA expression was blocked. These effects were associated with NF-κB subunit p65 hyperacetylation. Taken together, these data provide a rationale for clinical trials of safety and efficacy of KDACi in patients with autoimmune disease such as type 1 diabetes.

15-year follow-up of the Second Nordic Mantle Cell Lymphoma trial (MCL2): prolonged remissions without survival plateau.

In recent decades, the prognosis of Mantle Cell Lymphoma (MCL) has been significantly improved by intensified first-line regimens containing cytarabine, rituximab and consolidation with high-dose-therapy and autologous stem cell transplantation. One such strategy is the Nordic MCL2 regimen, developed by the Nordic Lymphoma Group. We here present the 15-year updated results of the Nordic MCL2 study after a median follow-up of 11·4 years: For all patients on an intent-to-treat basis, the median overall and progression-free survival was 12·7 and 8·5 years, respectively. The MCL International Prognostic Index (MIPI), biological MIPI, including Ki67 expression (MIPI-B) and the MIPI-B including mIR-18b expression (MIPI-B-miR), in particular, significantly divided patients into distinct risk groups. Despite very long response durations of the low and intermediate risk groups, we observed a continuous pattern of relapse and the survival curves never reached a plateau. In conclusion, despite half of the patients being still alive and 40% in first remission after more than 12 years, we still see an excess disease-related mortality, even among patients experiencing long remissions. Even though we consider the Nordic regimen as a very good choice of regimen, we recommend inclusion in prospective studies to explore the benefit of novel agents in the frontline treatment of MCL.

Quantification of oxidative stress phenotypes based on high-throughput growth profiling of protein kinase and phosphatase knockouts.

Cellular responses to oxidative stress are important for restoring redox balance and ensuring cell survival. Genetic defects in response factors can lead to impaired response to oxidative damage and contribute to disease and aging. In single cell organisms, such as yeasts, the integrity of the oxidative stress response can be observed through its influences on growth characteristics. In this study, we investigated the time-dependent batch growth effects as a function of oxidative stress levels in protein kinase and phosphatase deletion backgrounds of Saccharomyces cerevisiae. In total, 41 different protein kinases and phosphatase mutants were selected for their known activities in oxidative stress or other stress response pathways and were investigated for their dosage-dependent response to hydrogen peroxide. Detailed growth profiles were analyzed after the induction of stress for growth rate, lag time duration and growth efficiency, and by a novel method to identify stress-induced diauxic shift delay. This approach extracts more phenotypic information than traditional plate-based methods due to the assessment of time dynamics in the time scale of minutes. With this approach, we were able to identify surprisingly diverse sensitivity and resistance patterns as a function of gene knockout.

A DNA-binding-site landscape and regulatory network analysis for NAC transcription factors in Arabidopsis thaliana.

Target gene identification for transcription factors is a prerequisite for the systems wide understanding of organismal behaviour. NAM-ATAF1/2-CUC2 (NAC) transcription factors are amongst the largest transcription factor families in plants, yet limited data exist from unbiased approaches to resolve the DNA-binding preferences of individual members. Here, we present a TF-target gene identification workflow based on the integration of novel protein binding microarray data with gene expression and multi-species promoter sequence conservation to identify the DNA-binding specificities and the gene regulatory networks of 12 NAC transcription factors. Our data offer specific single-base resolution fingerprints for most TFs studied and indicate that NAC DNA-binding specificities might be predicted from their DNA-binding domain's sequence. The developed methodology, including the application of complementary functional genomics filters, makes it possible to translate, for each TF, protein binding microarray data into a set of high-quality target genes. With this approach, we confirm NAC target genes reported from independent in vivo analyses. We emphasize that candidate target gene sets together with the workflow associated with functional modules offer a strong resource to unravel the regulatory potential of NAC genes and that this workflow could be used to study other families of transcription factors.

Natural genetic variation impacts expression levels of coding, non-coding, and antisense transcripts in fission yeast.

Our current understanding of how natural genetic variation affects gene expression beyond well-annotated coding genes is still limited. The use of deep sequencing technologies for the study of expression quantitative trait loci (eQTLs) has the potential to close this gap. Here, we generated the first recombinant strain library for fission yeast and conducted an RNA-seq-based QTL study of the coding, non-coding, and antisense transcriptomes. We show that the frequency of distal effects (trans-eQTLs) greatly exceeds the number of local effects (cis-eQTLs) and that non-coding RNAs are as likely to be affected by eQTLs as protein-coding RNAs. We identified a genetic variation of swc5 that modifies the levels of 871 RNAs, with effects on both sense and antisense transcription, and show that this effect most likely goes through a compromised deposition of the histone variant H2A.Z. The strains, methods, and datasets generated here provide a rich resource for future studies.