Mycoparasitism related targets of Tmk1 indicate stimulating regulatory functions of this MAP kinase in Trichoderma atroviride.

Mycoparasitism is a key feature of Trichoderma (Hypocreales, Ascomycota) biocontrol agents. Recent studies of intracellular signal transduction pathways of the potent mycoparasite Trichoderma atroviride revealed the involvement of Tmk1, a mitogen-activated protein kinase (MAPK), in triggering the mycoparasitic response. We previously showed that mutants missing Tmk1 exhibit reduced mycoparasitic activity against several plant pathogenic fungi. In this study, we identified the most robustly regulated targets that were governed by Tmk1 during mycoparasitism using transcriptome and proteome profiling. Tmk1 mainly exerts a stimulating function for T. atroviride during its mycoparasitic interaction with the fungal plant pathogen Rhizoctonia solani, as reflected by 89% of strongly differently responding genes in the ∆tmk1 mutant compared to the wild type. Specifically, 54% of these genes showed strong downregulation in the response with a deletion of the tmk1 gene, whereas in the wild type the same genes were strongly upregulated during the interaction with the fungal host. These included the gene encoding the mycoparasitism-related proteinase Prb1; genes involved in signal transduction pathways such as a candidate coding for a conserved 14-3-3 protein, and a gene coding for Tmk2, the T. atroviride cell-wall integrity MAP kinase; genes encoding a specific siderophore synthetase, and multiple FAD-dependent oxidoreductases and aminotransferases. Due to the phosphorylating activity of Tmk1, different (phospho-)proteomics approaches were applied and identified proteins associated with cellular metabolism, energy production, protein synthesis and fate, and cell organization. Members of FAD- and NAD/NADP-binding-domain proteins, vesicular trafficking of molecules between cellular organelles, fungal translational, as well as protein folding apparatus were among others found to be phosphorylated by Tmk1 during mycoparasitism. Outstanding downregulation in the response of the ∆tmk1 mutant to the fungal host compared to the wild type at both the transcriptome and the proteome levels was observed for nitrilase, indicating that its defense and detoxification functions might be greatly dependent on Tmk1 during T. atroviride mycoparasitism. An intersection network analysis between the identified transcripts and proteins revealed a strong involvement of Tmk1 in molecular functions with GTPase and oxidoreductase activity. These data suggest that during T. atroviride mycoparasitism this MAPK mainly governs processes regulating cell responses to extracellular signals and those involved in reactive oxygen stress.

ViLoN-a multi-layer network approach to data integration demonstrated for patient stratification.

With more and more data being collected, modern network representations exploit the complementary nature of different data sources as well as similarities across patients. We here introduce the Variation of information fused Layers of Networks algorithm (ViLoN), a novel network-based approach for the integration of multiple molecular profiles. As a key innovation, it directly incorporates prior functional knowledge (KEGG, GO). In the constructed network of patients, patients are represented by networks of pathways, comprising genes that are linked by common functions and joint regulation in the disease. Patient stratification remains a key challenge both in the clinic and for research on disease mechanisms and treatments. We thus validated ViLoN for patient stratification on multiple data type combinations (gene expression, methylation, copy number), showing substantial improvements and consistently competitive performance for all. Notably, the incorporation of prior functional knowledge was critical for good results in the smaller cohorts (rectum adenocarcinoma: 90, esophageal carcinoma: 180), where alternative methods failed.

Silicone Breast Implant Surface Texture Impacts Gene Expression in Periprosthetic Fibrous Capsules.

BACKGROUND: Silicone breast implants with smooth outer shells are associated with higher rates of capsular contracture, whereas textured implants have been linked to the development of breast implant-associated anaplastic large cell lymphoma. By assessing the gene expression profile of fibrous capsules formed in response to smooth and textured implants, insight into the development of breast implant-associated abnormalities can be gained. METHODS: Miniature smooth or textured silicone implants were surgically inserted into female rats ( n = 10) and harvested for the surrounding capsules at postoperative week 6. RNA sequencing and quantitative polymerase chain reaction were performed to identify genes differentially expressed between smooth and textured capsules. For clinical correlation, the expression of candidate genes was assayed in implant capsules harvested from human patients with and without capsular contracture. RESULTS: Of 18,555 differentially expressed transcripts identified, three candidate genes were selected: matrix metalloproteinase-3 ( MMP3 ), troponin-T3 ( TNNT3 ), and neuregulin-1 ( NRG1 ). In textured capsules, relative gene expression and immunostaining of MMP3 and TNNT3 was up-regulated, whereas NRG1 was down-regulated compared to smooth capsules [mean relative fold change, 8.79 ( P = 0.0059), 4.81 ( P = 0.0056), and 0.40 ( P < 0.0001), respectively]. Immunostaining of human specimens with capsular contracture revealed similar gene expression patterns to those of animal-derived smooth capsules. CONCLUSIONS: An expression pattern of low MMP3 /low TNNT3 /high NRG1 is specifically associated with smooth implant capsules and human implant capsules with capsular contracture. The authors' clinically relevant breast implant rat model provides a strong foundation to further explore the molecular genetics of implant texture and its effect on breast implant-associated abnormalities. CLINICAL RELEVANCE STATEMENT: The authors have demonstrated that there are distinct gene expression profiles in response to smooth versus textured breast implants. Since surface texture may be linked to implant-related pathology, further molecular analysis of periprosthetic capsules may yield strategies to mitigate implant-related complications.

Expression Analyses in the Rachis Hint towards Major Cell Wall Modifications in Grape Clusters Showing Berry Shrivel Symptoms.

Berry shrivel (BS) is one of the prominent and still unresolved ripening physiological disorders in grapevine. The causes of BS are unclear, and previous studies focused on the berry metabolism or histological studies, including cell viability staining in the rachis and berries of BS clusters. Herein, we studied the transcriptional modulation induced by BS in the rachis of pre-symptomatic and symptomatic clusters with a custom-made microarray qPCR in relation to a previous RNASeq study of BS berries. Gene set analysis of transcript expression in symptomatic rachis tissue determined suppression of cell wall biosynthesis, which could also be confirmed already in pre-symptomatic BS rachis by CESA8 qPCR analyses, while in BS berries, a high number of SWITCH genes were suppressed at veraison. Additionally, genes associated with the cell wall were differently affected by BS in berries. A high percentage of hydrolytic enzymes were induced in BS grapes in rachis and berries, while other groups such as, e.g., xyloglucan endotransglucosylase/hydrolase, were suppressed in BS rachis. In conclusion, we propose that modulated cell wall biosynthesis and cell wall assembly in pre-symptomatic BS rachis have potential consequences for cell wall strength and lead to a forced degradation of cell walls in symptomatic grape clusters. The similarity to sugar starvation transcriptional profiles provides a link to BS berries, which are low in sugar accumulation. However, further studies remain necessary to investigate the temporal and spatial coordination in both tissues.

Molecular Mechanisms of Fetal Tendon Regeneration Versus Adult Fibrous Repair.

Tendinopathies are painful, disabling conditions that afflict 25% of the adult human population. Filling an unmet need for realistic large-animal models, we here present an ovine model of tendon injury for the comparative study of adult scarring repair and fetal regeneration. Complete regeneration of the fetal tendon within 28 days is demonstrated, while adult tendon defects remained macroscopically and histologically evident five months post-injury. In addition to a comprehensive histological assessment, proteome analyses of secretomes were performed. Confirming histological data, a specific and pronounced inflammation accompanied by activation of neutrophils in adult tendon defects was observed, corroborated by the significant up-regulation of pro-inflammatory factors, neutrophil attracting chemokines, the release of potentially tissue-damaging antimicrobial and extracellular matrix-degrading enzymes, and a response to oxidative stress. In contrast, secreted proteins of injured fetal tendons included proteins initiating the resolution of inflammation or promoting functional extracellular matrix production. These results demonstrate the power and relevance of our novel ovine fetal tendon regeneration model, which thus promises to accelerate research in the field. First insights from the model already support our molecular understanding of successful fetal tendon healing processes and may guide improved therapeutic strategies.

A comprehensive rat transcriptome built from large scale RNA-seq-based annotation.

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.

Sex-Specific Human Cardiomyocyte Gene Regulation in Left Ventricular Pressure Overload.

OBJECTIVE: To assess gene expression in cardiomyocytes isolated from patients with aortic stenosis, hypothesizing that maladaptive remodeling and inflammation-related genes are higher in male vs female patients. PATIENTS AND METHODS: In this study, 34 patients with aortic stenosis undergoing aortic valve replacement from March 20, 2016, through May 24, 2017, at the German Heart Centre in Berlin, Germany, were included. Isolated cardiomyocytes from interventricular septum samples were used for gene expression analysis. Clinical and echocardiographic data were collected preoperatively. RESULTS: Age, body mass index, systolic and diastolic blood pressure, comorbidities, and medication were similar between the 17 male and 17 female patients. The mean ± SD left ventricular end-diastolic diameter (52±9 vs 45±4 mm; P=.007) and posterior wall thickness (14.2±2.5 vs 12.1±1.6 mm; P=.03) were higher in male vs female patients, while ejection fraction was lower in male patients (49%±14% vs 59%±5%; P=.01). Focusing on structural genes involved in the development of cardiac hypertrophy and remodeling, we found that most were expressed higher in male vs female patients. Our modeling analysis revealed that 2 inflammation-related genes, CCN2 and NFKB1, were negatively related to ejection fraction, with this effect being male specific (P=.03 and P=.02, respectively). CONCLUSION: These findings provide novel insight into cardiomyocyte-specific molecular changes related to sex differences in pressure overload and a significant male-specific association between cardiac function and inflammation-related genes. Considering these sex differences may contribute toward a more accurate design of research and the development of more appropriate therapeutic approaches for both male and female patients.

Loss of the ribosomal RNA methyltransferase NSUN5 impairs global protein synthesis and normal growth.

Modifications of ribosomal RNA expand the nucleotide repertoire and thereby contribute to ribosome heterogeneity and translational regulation of gene expression. One particular m5C modification of 25S ribosomal RNA, which is introduced by Rcm1p, was previously shown to modulate stress responses and lifespan in yeast and other small organisms. Here, we report that NSUN5 is the functional orthologue of Rcm1p, introducing m5C3782 into human and m5C3438 into mouse 28S ribosomal RNA. Haploinsufficiency of the NSUN5 gene in fibroblasts from William Beuren syndrome patients causes partial loss of this modification. The N-terminal domain of NSUN5 is required for targeting to nucleoli, while two evolutionary highly conserved cysteines mediate catalysis. Phenotypic consequences of NSUN5 deficiency in mammalian cells include decreased proliferation and size, which can be attributed to a reduction in total protein synthesis by altered ribosomes. Strikingly, Nsun5 knockout in mice causes decreased body weight and lean mass without alterations in food intake, as well as a trend towards reduced protein synthesis in several tissues. Together, our findings emphasize the importance of single RNA modifications for ribosome function and normal cellular and organismal physiology.

QuaPra: Efficient transcript assembly and quantification using quadratic programming with Apriori algorithm.

RNA sequencing (RNA-seq) has greatly facilitated the exploring of transcriptome landscape for diverse organisms. However, transcriptome reconstruction is still challenging due to various limitations of current tools and sequencing technologies. Here, we introduce an efficient tool, QuaPra (Quadratic Programming combined with Apriori), for accurate transcriptome assembly and quantification. QuaPra could detect at least 26.5% more low abundance (0.1-1 FPKM) transcripts with over 2.1% increase of sensitivity and precision on simulated data compared to other currently popular tools. Moreover, around one-quarter more known transcripts were correctly assembled by QuaPra than other assemblers on real sequencing data. QuaPra is freely available at https://doi.org/www.megabionet.org/QuaPra/ .

The Gpr1-regulated Sur7 family protein Sfp2 is required for hyphal growth and cell wall stability in the mycoparasite Trichoderma atroviride.

Mycoparasites, e.g. fungi feeding on other fungi, are prominent within the genus Trichoderma and represent a promising alternative to chemical fungicides for plant disease control. We previously showed that the seven-transmembrane receptor Gpr1 regulates mycelial growth and asexual development and governs mycoparasitism-related processes in Trichoderma atroviride. We now describe the identification of genes being targeted by Gpr1 under mycoparasitic conditions. The identified gene set includes a candidate, sfp2, encoding a protein of the fungal-specific Sur7 superfamily, whose upregulation in T. atroviride upon interaction with a fungal prey is dependent on Gpr1. Sur7 family proteins are typical residents of membrane microdomains such as the membrane compartment of Can1 (MCC)/eisosome in yeast. We found that GFP-labeled Gpr1 and Sfp2 proteins show partly overlapping localization patterns in T. atroviride hyphae, which may point to shared functions and potential interaction during signal perception and endocytosis. Deletion of sfp2 caused heavily altered colony morphology, defects in polarized growth, cell wall integrity and endocytosis, and significantly reduced mycoparasitic activity, whereas sfp2 overexpression enhanced full overgrowth and killing of the prey. Transcriptional activation of a chitinase specific for hyphal growth and network formation and strong downregulation of chitin synthase-encoding genes were observed in Δsfp2. Taken together, these findings imply crucial functions of Sfp2 in hyphal morphogenesis of T. atroviride and its interaction with prey fungi.

Fetal articular cartilage regeneration versus adult fibrocartilaginous repair: secretome proteomics unravels molecular mechanisms in an ovine model.

Osteoarthritis (OA), a degenerative joint disease characterized by progressive cartilage degeneration, is one of the leading causes of disability worldwide owing to the limited regenerative capacity of adult articular cartilage. Currently, there are no disease-modifying pharmacological or surgical therapies for OA. Fetal mammals, in contrast to adults, are capable of regenerating injured cartilage in the first two trimesters of gestation. A deeper understanding of the properties intrinsic to the response of fetal tissue to injury would allow us to modulate the way in which adult tissue responds to injury. In this study, we employed secretome proteomics to compare fetal and adult protein regulation in response to cartilage injury using an ovine cartilage defect model. The most relevant events comprised proteins associated with the immune response and inflammation, proteins specific for cartilage tissue and cartilage development, and proteins involved in cell growth and proliferation. Alarmins S100A8, S100A9 and S100A12 and coiled-coil domain containing 88A (CCDC88A), which are associated with inflammatory processes, were found to be significantly upregulated following injury in adult, but not in fetal animals. By contrast, cartilage-specific proteins like proteoglycan 4 were upregulated in response to injury only in fetal sheep postinjury. Our results demonstrate the power and relevance of the ovine fetal cartilage regeneration model presented here for the first time. The identification of previously unrecognized modulatory proteins that plausibly affect the healing process holds great promise for potential therapeutic interventions.

Transcriptomic changes in CHO cells after adaptation to suspension growth in protein-free medium analysed by a species-specific microarray.

Chinese Hamster Ovary (CHO) cells are the preferred cell line for production of biopharmaceuticals. These cells are capable to grow without serum supplementation, but drastic changes in their phenotype occur during adaptation to protein-free growth, which typically include the change to a suspension phenotype with reduced growth rate. A possible approach to understand this transformation, with the intention to counteract the reduction in growth by targeted supplementation of protein-free media, is gene expression profiling. The increasing availability of genome-scale data for CHO now facilitates quests for a better understanding of metabolic pathways and gene networks. So far, systematic large-scale expression profiling in CHO cells by microarray was limited due to lack of publicly available array designs and limitations of alternative approaches. Based on the recent release of CHO and Chinese Hamster genome sequences, including an annotated RefSeq genome, we have constructed a publicly available microarray design for effective genome-scale expression profiling. The design employed microarray probes optimized for uniformity, sensitivity, and specificity, with probe properties computed using the latest thermodynamic models. We validated the platform in an analysis of gene expression changes in response to serum-free adaptation. The observed effects on the lipid metabolism as well as on nucleotide synthesis were used to successfully select media supplements that were able to increase growth rate.

Sensitivity, specificity, and reproducibility of RNA-Seq differential expression calls.

BACKGROUND: The MAQC/SEQC consortium has recently compiled a key benchmark that can serve for testing the latest developments in analysis tools for microarray and RNA-seq expression profiling. Such objective benchmarks are required for basic and applied research, and can be critical for clinical and regulatory outcomes. Going beyond the first comparisons presented in the original SEQC study, we here present extended benchmarks including effect strengths typical of common experiments. RESULTS: With artefacts removed by factor analysis and additional filters, for genome scale surveys, the reproducibility of differential expression calls typically exceed 80% for all tool combinations examined. This directly reflects the robustness of results and reproducibility across different studies. Similar improvements are observed for the top ranked candidates with the strongest relative expression change, although here some tools clearly perform better than others, with typical reproducibility ranging from 60 to 93%. CONCLUSIONS: In our benchmark of alternative tools for RNA-seq data analysis we demonstrated the benefits that can be gained by analysing results in the context of other experiments employing a reference standard sample. This allowed the computational identification and removal of hidden confounders, for instance, by factor analysis. In itself, this already substantially improved the empirical False Discovery Rate (eFDR) without changing the overall landscape of sensitivity. Further filtering of false positives, however, is required to obtain acceptable eFDR levels. Appropriate filters noticeably improved agreement of differentially expressed genes both across sites and between alternative differential expression analysis pipelines. REVIEWERS: An extended abstract of this research paper was selected for the CAMDA Satellite Meeting to ISMB 2015 by the CAMDA Programme Committee. The full research paper then underwent one round of Open Peer Review under a responsible CAMDA Programme Committee member, Lan Hu, PhD (Bio-Rad Laboratories, Digital Biology Center-Cambridge). Open Peer Review was provided by Charlotte Soneson, PhD (University of Zürich) and Michal Okoniewski, PhD (ETH Zürich). The Reviewer Comments section shows the full reviews and author responses.

MALDI TOF/TOF-Based Approach for the Identification of d- Amino Acids in Biologically Active Peptides and Proteins.

Several biologically active peptides contain a d- amino acid in a well-defined position, which is position 2 in all peptide epimers isolated to date from vertebrates and also some from invertebrates. The detection of such D- residues by standard analytical techniques is challenging. In tandem mass spectrometric (MS) analysis, although fragment masses are the same for all stereoisomers, peak intensities are known to depend on chirality. Here, we observe that the effect of a d- amino acid in the second N-terminal position on the fragmentation pattern in matrix assisted laser desorption time-of-flight spectrometry (MALDI-TOF/TOF MS) strongly depends on the peptide sequence. Stereosensitive fragmentation (SF) is correlated to a neighborhood effect, but the d- residue also exerts an overall effect influencing distant bonds. In a fingerprint analysis, multiple peaks can thus serve to identify the chirality of a sample in short time and potentially high throughput. Problematic variations between individual spots could be successfully suppressed by cospotting deuterated analogues of the epimers. By identifying the [d-Leu2] isomer of the predicted peptide GH-2 (gene derived bombininH) in skin secretions of the toad Bombina orientalis, we demonstrated the analytical power of SF-MALDI-TOF/TOF measurements. In conclusion, SF-MALDI-TOF/TOF MS combines high sensitivity, versatility, and the ability to complement other methods.

Experiences with workflows for automating data-intensive bioinformatics.

High-throughput technologies, such as next-generation sequencing, have turned molecular biology into a data-intensive discipline, requiring bioinformaticians to use high-performance computing resources and carry out data management and analysis tasks on large scale. Workflow systems can be useful to simplify construction of analysis pipelines that automate tasks, support reproducibility and provide measures for fault-tolerance. However, workflow systems can incur significant development and administration overhead so bioinformatics pipelines are often still built without them. We present the experiences with workflows and workflow systems within the bioinformatics community participating in a series of hackathons and workshops of the EU COST action SeqAhead. The organizations are working on similar problems, but we have addressed them with different strategies and solutions. This fragmentation of efforts is inefficient and leads to redundant and incompatible solutions. Based on our experiences we define a set of recommendations for future systems to enable efficient yet simple bioinformatics workflow construction and execution.

Polarization-Sensitive Optical Coherence Tomography and Conventional Retinal Imaging Strategies in Assessing Foveal Integrity in Geographic Atrophy.

PURPOSE: To compare current imaging methods with respect to their ability to detect the condition of the fovea in patients with geographic atrophy (GA). METHODS: The retinas of 176 eyes with GA were imaged using two spectral-domain optical coherence tomography (SD-OCT) systems, Cirrus HD-OCT and Spectralis HRA+OCT, and fundus autofluorescence (FAF) and infrared imaging (IR) was used in the scanning laser ophthalmoscope (SLO) mode. Polarization-sensitive OCT (PS-OCT), which selectively visualizes the RPE in addition to SD-OCT features, was used to image 95 eyes. Geographic atrophy lesions were categorized as fovea spared, involved, or not quantifiable (grades 0, 1, and 2). Morphologic gradings were subsequently correlated with best-corrected visual acuity (BCVA) measurements to independently identify the corresponding functional condition of the fovea. Cohen's κ statistics with a bootstrap method was applied to compare retinal imaging methods. RESULTS: In PS-OCT, 84% of eyes with BCVA greater than or equal to 20/40 were detected, whereas in conventional retinal imaging the rate ranged from 27% in FAF to 45% in the SD-OCT segment. Cohen's κ statistics revealed significant differences between the gradings of PS-OCT and conventional imaging with κ = 0.488 and a global Hotelling's T2 statistic of 17.9 with a P value of P = 0.003. Statistical tests revealed no statistically significant differences between the conventional retinal imaging modalities. CONCLUSIONS: Polarization-sensitive OCT can better allow correct grading of the fovea in relation to BCVA and identify foveal sparing than other imaging modalities. The differences in imaging precision should be considered in diagnostic and therapeutic evaluations.

Microarray profiling of preselected CHO host cell subclones identifies gene expression patterns associated with increased production capacity.

Over the last three decades, product yields from CHO cells have increased dramatically, yet specific productivity (qP) remains a limiting factor. In a previous study, using repeated cell-sorting, we have established different host cell subclones that show superior transient qP over their respective parental cell lines (CHO-K1, CHO-S). The transcriptome of the resulting six cell lines in different biological states (untransfected, mock transfected, plasmid transfected) was first explored by hierarchical clustering and indicated that gene activity associated with increased qP did not stem from a certain cellular state but seemed to be inherent for a high qP host line. We then performed a novel gene regression analysis identifying drivers for an increase in qP. Genes significantly implicated were first systematically tested for enrichment of GO terms using a Bayesian approach incorporating the hierarchical structure of the GO term tree. Results indicated that specific cellular components such as nucleus, ER, and Golgi are relevant for cellular productivity. This was complemented by targeted GSA that tested functionally homogeneous, manually curated subsets of KEGG pathways known to be involved in transcription, translation, and protein processing. Significantly implicated pathways included mRNA surveillance, proteasome, protein processing in the ER and SNARE interactions in vesicular transport.