Search

Base editing is an efficient and precise gene-editing technique, by which a single base can be changed without introducing double-strand breaks, and it is currently widely used in studies of various species. In this study, we used hA3A-BE3-Y130F to simultaneously introduce premature stop codons (TAG, TGA, and TAA) into three tumor suppressor genes, TP53, PTEN, and APC, in large white porcine fetal fibroblasts (PFFs). Among the isolated 290 single-cell colonies, 232 (80%) had premature stop codons in all the three genes. C−to−T conversion was found in 98.6%, 92.8%, and 87.2% of these cell colonies for TP53, PTEN, and APC, respectively. High frequencies of bystander C−to−T edits were observed within the editing window (positions 3–8), and there were nine (3.01%) clones with the designed simultaneous three-gene C−to−T conversion without bystander conversion. C−to−T conversion outside the editing window was found in 9.0%, 14.1%, and 26.2% of the 290 cell colonies for TP53, PTEN, and APC, respectively. Low-frequency C−to−G or C−to−A transversion occurred in APC. The mRNA levels of the three genes showed significant declines in triple-gene-mutant (Tri-Mut) cells as expected. No PTEN and a significantly lower (p < 0.05) APC protein expression were detected in Tri-Mut cells. Interestingly, the premature stop codon introduced into the TP53 gene did not eliminate the expression of its full-length protein in the Tri-Mut cells, suggesting that stop codon read-through occurred. Tri-Mut cells showed a significantly higher (p < 0.05) proliferation rate than WT cells. Furthermore, we identified 1418 differentially expressed genes (DEGs) between the Tri-Mut and WT groups, which were mainly involved in functions such as tumor progression, cell cycle, and DNA repair. This study indicates that hA3A-BE3-Y130F can be a powerful tool to create diverse knockout cell models without double-strand breaks (DSBs), with further possibilities to produce porcine models with various purposes.

Many pregnancy disorders, including early-onset preeclampsia (EOPE), are associated with defects in placental trophoblast cell invasion and differentiation during early placental development. Bone morphogenetic protein 2 (BMP2) belongs to the TGF-β superfamily and controls various physiological and developmental processes. However, the expression of BMP2 in the placenta and underlying molecular mechanisms of how BMP2 regulates trophoblast function remain unclear. In this study, we analyzed several publicly available microarray and RNA-seq datasets and revealed differences in expression of TGF-β superfamily members between gestational age-matched non-preeclamptic control and EOPE placentas. Importantly, BMP2 levels were significantly reduced in EOPE placentas compared with controls, and RNAscope in situ hybridization further demonstrated BMP2 expression was disrupted in EOPE placental villi. To explore the molecular mechanisms of BMP2-regulated early trophoblast differentiation, we examined BMP2 expression in first-trimester human placenta and found it to be localized to all subtypes of trophoblasts and the decidua. RNA-seq analysis on control and BMP2-treated primary human trophoblast cells identified 431 differentially expressed genes, including several canonical TGF-β/BMP signaling targets (BAMBI, ID1, INHBA, IGFBP3). Gene ontology annotations revealed that differentially expressed genes were involved in cell adhesion and extracellular matrix organization. Furthermore, we identified adhesion molecule with IgG-like domain 2 (AMIGO2) as a novel target for BMP2 that contributed to BMP2-induced trophoblast invasion and endothelial-like tube formation. Overall, our findings provide insight into the molecular processes controlled by BMP2 during early placental development that may contribute to the pathogenesis of EOPE.

Betaine is trimethylglycine and a universal methyl donor which could provide methyl and glycine for cells and animals. As a new star in epigenetics, N6-Methyladenosine has been reported to regulate multiple biological activities, but the regulatory mechanism of betaine on N6-Methyladenosine as well as myogenesis was little studied. In this study, we treated chicken primary myoblast cells with different concentrations of betaine (0, 10, 25, and 50 mmol/L) and found that myoblast cell proliferation was inhibited, although the cell cycle was promoted in the S phase by betaine, where the myotube area was increased as well as the differentiation marker genes MyoD, MyoG, MyHC, Myomarker, and Ckm. RNA sequencing obtained a total of 61 differentially expressed genes (DEGs); DEGs caused by 50 mmol/L betaine were mainly enriched in the regulation of skeletal muscle tissue regeneration and some amino acid metabolic processes. The gene expression pattern trends of all DEGs were mainly clustered into 2 profiles, with the increase in betaine concentration, the gene expression pattern either increased or decreased continuously. Overall, a low concentration betaine can increase the N6-Methyladenosine modification level and myotube area but depresses myoblast cell proliferation in vitro.

Hirschmanniella mucronata is a plant‐parasitic nematode that is widespread in rice production areas and causes 10–25% yield losses a year on average. Here, we investigated the mechanism of resistance to this nematode by comparing the transcriptomes of roots from resistant (Jiabali) and susceptible (Bawangbian) varieties of rice. Of 39 233 unigenes, 2243. exhibited altered total expression levels between control and infected resistant and susceptible varieties. Significant differences were observed in the expression levels of genes related to stress, peptidase regulation or inhibition, oxidoreductase activity, peroxidase activity and antioxidant activity. The up‐regulated genes related to plant secondary metabolites, such as phenylpropanoid, lignin, cellulose or hemicellulose, may result in an increase in the degree of resistance of Jiabali to the H. mucronata infection compared with that of Bawangbian by affecting cell wall organization or biogenesis. Of the genes that responded similarly to H. mucronata infection, ~252 (~76.59%) showed greater changes (whether induced or suppressed) in RN155 (susceptible varieties infected by rice root nematode) than in RN51 (resistance varieties infected by rice root nematode). Nineteen pathogenesis‐related genes belonging to nine pathogenesis‐related gene families were significantly induced by H. mucronata in the infected roots of Jiabali and Bawangbian, and 13 differentially expressed genes showed changes in their abundance only in the susceptible Bawangbian variety. This study may help enhance our understanding of the mechanisms underlying plant resistance to nematodes.

Gluconacetobacter diazotrophicus has been the focus of several studies aiming to understand the mechanisms behind this endophytic diazotrophic bacterium. The present study is the first global analysis of the early transcriptional response of exponentially growing G. diazotrophicus to iron, an essential cofactor for many enzymes involved in various metabolic pathways. RNA-seq, targeted gene mutagenesis and computational motif discovery tools were used to define the G. diazotrophicusfur regulon. The data analysis showed that genes encoding functions related to iron homeostasis were significantly upregulated in response to iron limitations. Certain genes involved in secondary metabolism were overexpressed under iron-limited conditions. In contrast, it was observed that the expression of genes involved in Fe-S cluster biosynthesis, flagellar biosynthesis and type IV secretion systems were downregulated in an iron-depleted culture medium. Our results support a model that controls transcription in G. diazotrophicus by fur function. The G. diazotrophicusfur protein was able to complement an E. colifur mutant. These results provide new insights into the effects of iron on the metabolism of G. diazotrophicus, as well as demonstrate the essentiality of this micronutrient for the main characteristics of plant growth promotion by G. diazotrophicus.

Parkinson’s disease (PD) is the second most common age-related neurodegenerative disease. It is presently only accurately diagnosed at an advanced stage by a series of motor deficits, which are predated by a litany of non-motor symptoms manifesting over years or decades. Aberrant epigenetic modifications exist across a range of diseases and are non-invasively detectable in blood as potential markers of disease. We performed comparative analyses of the methylome and transcriptome in blood from PD patients and matched controls. Our aim was to characterize DNA methylation and gene expression patterns in whole blood from PD patients as a foundational step toward the future goal of identifying molecular markers that could predict, accurately diagnose, or track the progression of PD. We found that differentially expressed genes (DEGs) were involved in the processes of transcription and mitochondrial function and that PD methylation profiles were readily distinguishable from healthy controls, even in whole-blood DNA samples. Differentially methylated regions (DMRs) were functionally varied, including near transcription factor nuclear transcription factor Y subunit alpha (NFYA), receptor tyrosine kinase DDR1, RING finger ubiquitin ligase (RNF5), acetyltransferase AGPAT1, and vault RNA VTRNA2-1. Expression quantitative trait methylation sites were found at long non-coding RNA PAX8-AS1 and transcription regulator ZFP57 among others. Functional epigenetic modules were highlighted by IL18R1, PTPRC, and ITGB2. We identified patterns of altered disease-specific DNA methylation and associated gene expression in whole blood. Our combined analyses extended what we learned from the DEG or DMR results alone. These studies provide a foundation to support the characterization of larger sample cohorts, with the goal of building a thorough, accurate, and non-invasive molecular PD biomarker.

Messenger RNA-seq data of tomato lines affected in the LTP2 gene. RNA isolation To isolate RNA, plant tissues were flash frozen in liquid nitrogen immediately after harvesting and stored at −80°C prior to RNA extraction. Frozen pellets were grounded by using a mortar and pestle before immersion in QIAzol Lysis Reagent (Qiagen). RNA was isolated and purified into two separate fractions (>200 nt and 200 nt RNA an on-column treatment was included using the RNase-free DNase set (Qiagen). mRNA sequencing protocol A poly-A enrichment was performed on the >200 nt RNA fraction using the NEBNext Poly(A) mRNA Magnetic Isolation Module (New England BioLabs). RNA-Seq libraries were generated according to the manufacturers’ protocols using the NEBNext Ultra II Directional RNA Library Prep Kit for Illumina and NEBNext Multiplex Oligos for Illumina (Unique Dual Index Primer Pairs) (New England BioLabs). The size distribution of the libraries with indexed adapters was assessed using a 2200 TapeStation System with Agilent D1000 ScreenTapes (Agilent Technologies). The libraries were quantified on a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific) using the NEBNext Library Quant Kit for Illumina (New England BioLabs) according to the instructions of the manufacturer. The libraries were clustered and sequenced (75 bp) on a NextSeq 550 Sequencing System (Illumina) using a NextSeq 500/550 High Output Kit v2.5 (75 Cycles) (Illumina). The sample list can be found in list_of_samples.csv.

Tra catfish (Pangasianodon hypophthalmus), also known as striped catfish, is a facultative air-breather that uses its swim bladder as an air-breathing organ (ABO). A related species in the same order (Siluriformes), channel catfish (Ictalurus punctatus), does not possess an ABO and thus cannot breathe in the air. Tra and channel catfish serve as great comparative models for investigating possible genetic underpinnings of aquatic to land transitions, as well as for understanding genes that are crucial for the development of the swim bladder and the function of air-breathing in tra catfish. In this study, hypoxia challenge and microtomy experiments collectively revealed critical time points for the development of the air-breathing function and swim bladder in tra catfish. Seven developmental stages in tra catfish were selected for RNA-seq analysis based on their transition to a stage that could live at 0 ppm oxygen. More than 587 million sequencing clean reads were generated, and a total of 21,448 unique genes were detected. A comparative genomic analysis between channel catfish and tra catfish revealed 76 genes that were present in tra catfish, but absent from channel catfish. In order to further narrow down the list of these candidate genes, gene expression analysis was performed for these tra catfish-specific genes. Fourteen genes were inferred to be important for air-breathing. Of these,HRG,GRP, andCX3CL1were identified to be the most likely genes related to air-breathing ability in tra catfish. This study provides a foundational data resource for functional genomic studies in air-breathing function in tra catfish and sheds light on the adaptation of aquatic organisms to the terrestrial environment.

The occurrence of adventitious roots is necessary for the survival of cuttings. In this study, comparative transcriptome analysis between two ramie (Boehmeria nivea L.) varieties with different adventitious root (AR) patterns was performed by mRNA-Seq before rooting (control, CK) and 10 days water-induced adventitious rooting (treatment, T) to reveal the regulatory mechanism of rooting. Characterization of the two ramie cultivars, Zhongzhu No 2 (Z2) and Huazhu No 4 (H4), indicated that Z2 had a high adventitious rooting rate but H4 had a low rooting rate. Twelve cDNA libraries of the two varieties were constructed, and a total of 26,723 genes were expressed. In the non-water culture condition, the number of the distinctive genes in H4 was 2.7 times of that in Z2, while in the water culture condition, the number of the distinctive genes in Z2 was nearly 2 times of that in H4. A total of 4411 and 5195 differentially expressed genes (DEGs) were identified in the comparison of H4CK vs. H4T and Z2CK vs. Z2T, respectively. After the water culture, more DEGs were upregulated in Z2, but more DEGs were downregulated in H4. Gene ontology (GO) functional analysis of the DEGs indicated that the polysaccharide metabolic process, carbohydrate metabolic process, cellular carbohydrate metabolic process, cell wall macromolecule metabolic process, and photosystem GO terms were distinctively significantly enriched in H4. Simultaneously, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis, photosynthesis antenna proteins, and starch and sucrose metabolism pathways were distinctively significantly enriched in H4. Moreover, KEGG analysis showed that jasmonic acid (JA) could interact with ethylene to regulate the occurrence and number of AR in Z2. This study reveals the transcriptomic divergence of two ramie varieties with high and low adventitious rooting rates, and provides insights into the molecular regulatory mechanism of AR formation in ramie.

Background Chinese Yunnan Province, located in the Yunnan–Guizhou Plateau, is a famous tourist paradise where acute high-altitude illness common occurs among lowland people visitors due to non-acclimatization to the acute hypobaric hypoxia (AHH) conditions. Traditional Chinese medicine, such as Qi-Long-Tian (QLT) formula, has shown effectiveness and safety in the treatment of acute high-altitude diseases. The aim of this study was to clarify the therapeutic mechanisms of this traditional formula using a rat model in a simulated plateau environment. Methods Following testing, lung tissue samples were evaluated by hematoxylin–eosin staining and for biochemical characteristics. mRNA-Seq was used to compare differentially expressed genes in control rats, and in rats exposed to AHH and AHH with QLT treatment. Results Inflammation-related effectors induced following QLT treatment for AHH included MMP9 and TIMP1, and involved several phosphorylation signaling pathways implicated in AHH pathogenesis such as PI3K/AKT and MAPK signaling. Conclusion This study provides insights into the major signaling pathways induced by AHH and in the protective mechanisms involved in QLT formula activity.

Background: Whether retinal microglia can maintain or restore immune homeostasis during and after inflammation is unclear. We performed single-eye mRNA-sequencing on microglia at different timepoints following a single inflammatory stimulus to characterise their transcriptome during and after resolution of endotoxin-induced uveitis (EIU).Experimental Approach: Cx3cr1CreER:R26-tdTomato (C57BL/6) male heterozygotes were administered tamoxifen via different regimes at 4–5 weeks of age. Four weeks post-tamoxifen, mice were injected intravitreally with 10 ng lipopolysaccharide (endotoxin induced uveitis, EIU). Six-hundred retinal microglia were obtained by FACS from individual naïve retinas and at 4 h, 18 h, and 2 weeks following EIU induction. Samples were sequenced to a depth of up to 16.7 million reads using the SMART-Seq v4 Ultra Low Input RNA kit. The data was analysed using Partek software and Ingenuity Pathway Analysis. Genes were considered differentially-expressed (DEG) if the FDR step-up p-value was ≤0.05 and the fold-change was ≥±2.Results: Flow cytometric analysis indicates that the Cx3cr1CreER:R26-tdTomato strain is both sensitive (>95% tagging) and specific (>95% specificity) for microglia when tamoxifen is administered topically to the eye for 3 days. During “early” activation, 613 DEGs were identified. In contrast, 537 DEGs were observed during peak cellular infiltrate and none at 2 weeks, compared to baseline controls (1,069 total unique DEGs). Key marker changes were validated by qPCR, flow cytometry, and fluorescence microscopy. C5AR1 was identified and validated as a robust marker of differentiating microglial subsets during an LPS response.Conclusion: Using EIU to provide a single defined inflammatory stimulus, mRNA-Seq identified acute transcriptional changes in retinal microglia which returned to their original transcriptome after 2 weeks. Yolk-sac derived microglia are capable of restoring their homeostatic state after acute inflammation.

Background Biological adhesion (bioadhesion), enables organisms to attach to surfaces as well as to a range of other targets. Bioadhesion evolved numerous times independently and is ubiquitous throughout the kingdoms of life. To date, investigations have focussed on various taxa of animals, plants and bacteria, but the fundamental processes underlying bioadhesion and the degree of conservation in different biological systems remain poorly understood. This study had two aims: 1) To characterise tissue-specific gene regulation in the pedal disc of the model cnidarian Exaiptasia pallida, and 2) to elucidate putative genes involved in pedal disc adhesion. Results Five hundred and forty-seven genes were differentially expressed in the pedal disc compared to the rest of the animal. Four hundred and twenty-seven genes were significantly upregulated and 120 genes were significantly downregulated. Forty-one condensed gene ontology terms and 19 protein superfamily classifications were enriched in the pedal disc. Eight condensed gene ontology terms and 11 protein superfamily classifications were depleted. Enriched superfamilies were consistent with classifications identified previously as important for the bioadhesion of unrelated marine invertebrates. A host of genes involved in regulation of extracellular matrix generation and degradation were identified, as well as others related to development and immunity. Ab initio prediction identified 173 upregulated genes that putatively code for extracellularly secreted proteins. Conclusion The analytical workflow facilitated identification of genes putatively involved in adhesion, immunity, defence and development of the E. pallida pedal disc. When defence, immunity and development-related genes were identified, those remaining corresponded most closely to formation of the extracellular matrix (ECM), implicating ECM in the adhesion of anemones to surfaces. This study therefore provides a valuable high-throughput resource for the bioadhesion community and lays a foundation for further targeted research to elucidate bioadhesion in the Cnidaria. Electronic supplementary material The online version of this article (10.1186/s12864-019-5917-5) contains supplementary material, which is available to authorized users.

Background Feline injection-site sarcoma (FISS), an aggressive iatrogenic subcutaneous malignancy, is challenging to manage clinically and little is known about the molecular basis of its pathogenesis. Tumor transcriptome profiling has proved valuable for gaining insights into the molecular basis of cancers and for identifying new therapeutic targets. Here, we report the first study of the FISS transcriptome and the first cross-species comparison of the FISS transcriptome with those of anatomically similar soft-tissue sarcomas in dogs and humans. Methods Using high-throughput short-read paired-end sequencing, we comparatively profiled FISS tumors vs. normal tissue samples as well as cultured FISS-derived cell lines vs. skin-derived fibroblasts. We analyzed the mRNA-seq data to compare cancer/normal gene expression level, identify biological processes and molecular pathways that are associated with the pathogenesis of FISS, and identify multimegabase genomic regions of potential somatic copy number alteration (SCNA) in FISS. We additionally conducted cross-species analyses to compare the transcriptome of FISS to those of soft-tissue sarcomas in dogs and humans, at the level of cancer/normal gene expression ratios. Results We found: (1) substantial differential expression biases in feline orthologs of human oncogenes and tumor suppressor genes suggesting conserved functions in FISS; (2) a genomic region with recurrent SCNA in human sarcomas that is syntenic to a feline genomic region of probable SCNA in FISS; and (3) significant overlap of the pattern of transcriptional alterations in FISS with the patterns of transcriptional alterations in soft-tissue sarcomas in humans and in dogs. We demonstrated that a protein, BarH-like homeobox 1 (BARX1), has increased expression in FISS cells at the protein level. We identified 11 drugs and four target proteins as potential new therapies for FISS, and validated that one of them (GSK-1059615) inhibits growth of FISS-derived cells in vitro. Conclusions (1) Window-based analysis of mRNA-seq data can uncover SCNAs. (2) The transcriptome of FISS-derived cells is highly consistent with that of FISS tumors. (3) FISS is highly similar to soft-tissue sarcomas in dogs and humans, at the level of gene expression. This work underscores the potential utility of comparative oncology in improving understanding and treatment of FISS. Electronic supplementary material The online version of this article (10.1186/s12885-019-5501-z) contains supplementary material, which is available to authorized users.

The development of techniques for measuring gene expression globally has greatly expanded our understanding of gene regulatory mechanisms in depth and scale. We can now quantify every intermediate and transition in the canonical pathway of gene expression-from DNA to mRNA to protein-genome-wide. Employing such measurements in parallel can produce rich datasets, but extracting the most information requires careful experimental design and analysis. Here, we argue for the value of genome-wide studies that measure multiple outputs of gene expression over many timepoints during the course of a natural developmental process. We discuss our findings from a highly parallel gene expression dataset of meiotic differentiation, and those of others, to illustrate how leveraging these features can provide new and surprising insight into fundamental mechanisms of gene regulation.

The emergence and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in significant global morbidity, mortality, and societal disruption. A better understanding of virus-host interactions may potentiate therapeutic insights toward limiting this infection. Here we investigated the dynamics of the systemic response to SARS-CoV-2 in hamsters by histological analysis and transcriptional profiling. Infection resulted in consistently high levels of virus in the upper and lower respiratory tracts and sporadic occurrence in other distal tissues. A longitudinal cohort revealed a wave of inflammation, including a type I interferon (IFN-I) response, that was evident in all tissues regardless of viral presence but was insufficient to prevent disease progression. Bolstering the antiviral response with intranasal administration of recombinant IFN-I reduced viral disease, prevented transmission, and lowered inflammation in vivo. This study defines the systemic host response to SARS-CoV-2 infection and supports use of intranasal IFN-I as an effective means of early treatment., Graphical abstract, The host response to SARS-CoV-2 results in significant inflammation. To understand this biology, Hoagland et al. utilize infected hamsters to elucidate transcriptional footprints across tissues longitudinally, showing an inflammatory response beyond the site of acute replication. Local administration of IFN-I reduces virus load and improves immune infiltrate.

Naringenin, a natural flavonoid widely found in citrus fruits, has been reported to possess anti-oxidant, anti-inflammatory, and hepatoprotective properties as a natural dietary supplement. However, the regulatory mechanism of naringenin in human liver remains unclear. In the present study, messenger RNA sequencing (mRNA-seq), microRNA sequencing (miRNA-seq), and real-time qPCR were used to distinguish the expression differences between control and naringenin-treated HepaRG cells. We obtained 1037 differentially expressed mRNAs and 234 miRNAs. According to the target prediction and integration analysis in silico, we found 20 potential miRNA-mRNA pairs involved in liver metabolism. This study is the first to provide a perspective of miRNA–mRNA interactions in the regulation of naringenin via an integrated analysis of mRNA-seq and miRNA-seq in HepaRG cells, which further characterizes the nutraceutical value of naringenin as a food additive.

Acquired drug resistance is the major reason why patients fail to respond to cancer therapies. It is a challenging task to determine the tipping point of endocrine resistance and detect the associated molecules. Derived from new systems biology theory, the dynamic network biomarker (DNB) method is designed to quantitatively identify the tipping point of a drastic system transition and can theoretically identify DNB genes that play key roles in acquiring drug resistance. We analyzed time-course mRNA sequence data generated from the tamoxifen-treated estrogen receptor (ER)-positive MCF-7 cell line, and identified the tipping point of endocrine resistance with its leading molecules. The results show that there is interplay between gene mutations and DNB genes, in which the accumulated mutations eventually affect the DNB genes that subsequently cause the change of transcriptional landscape, enabling full-blown drug resistance. Survival analyses based on clinical datasets validated that the DNB genes were associated with the poor survival of breast cancer patients. The results provided the detection for the pre-resistance state or early signs of endocrine resistance. Our predictive method may greatly benefit the scheduling of treatments for complex diseases in which patients are exposed to considerably different drugs and may become drug resistant.

Female adults of many mosquito species possess distinct physiological features adapting to blood feeding for successful reproduction. The disease pathogens that are transmitted by mosquitoes have evolved to take advantages of the indispensable blood feedings to complete their transmission cycles and to survive attacks from the mosquito's innate immune system. Normal egg development and mosquito immunity are tightly controlled by tissue- and stage-specific gene expression and coordinated by many signal molecules in the mosquito. Understanding gene regulation affecting mosquito reproduction and malaria parasites infection is of paramount importance for developing novel malaria control strategies. A growing body of evidence indicates that microRNAs (miRNAs) are involved in egg maturation and immune reactions against invading pathogens in mosquitoes. However, the molecular mechanisms by which specific miRNAs selectively modulate reproduction and the survival of pathogens are largely unknown. The miRNA-induced gene-silencing pathway in mosquitoes was mostly extrapolated from the studies of flies. To explore the dynamics of miRNAs in reproduction, I used small RNAs sequencing to monitor miRNAs expression and their association with Argonaute 1 (Ago1) and Argonaute 2 (Ago2) in the malaria mosquito Anopheles gambiae (An. gambiae) during the 72-h period immediately after blood feeding. I found the abundance and Ago loading of most of the mature miRNAs were relatively stable after blood ingestion. However, miRNAs of the miR-309/286/2944 cluster were considerably upregulated after blood feeding. I confirmed that miR-309 is essential for normal egg development by depletion of endogenous miR-309 with a specific antagomir. In addition, my results showed that the Ago association of some miRNAs was not proportional to their cellular abundance implying additional regulation at miRNA integration. To investigate the functional roles of miRNAs and define context-dependent miRNA-mRNA interactions during the reproductive process, I have applied an innovative experimental approach to study miRNA-mRNA interactome. CLEAR (covalent ligation of endogenous Argonaute-bound RNAs)-CLIP can generate miRNA-mRNA chimeras from UV-irradiation stabilized Ago-miRNA-mRNA complex. My results have defined tens of thousands of miRNA-mRNA interactions in mosquitoes, including novel targets for mosquito-specific miRNAs. Verification of the predicted interactions using mRNA-seq, ribosome-profiling, and luciferase reporter assay revealed a reliable miRNA-mRNA interaction network. Based on the detected interactions, I refined the paring rules for mosquito miRNAs and illustrated the dynamic pairing between different regions of miRNAs with their targets in vivo. The miRNA-mRNA interactions were compared using this approach at multiple time points before and after blood feeding. Importantly, this study showed that the interactions were dynamic and enriched in genes that are involved in metabolisms, supporting the proposed functions of miRNAs in coordinating the gene regulation in mosquito reproduction. Plasmodium falciparum (P. falciparum) is a major human malaria parasite. To understand the functions of miRNAs in the mosquito resistance to Plasmodium infection, we analyzed the miRNA-mRNA interactions after female mosquitoes taking a P. falciparum-infected blood meal or an uninfected blood meal. Comparison of the interactions revealed enhanced miRNA-mRNA interactions after P. falciparum infection involving a group of immunity-related genes. In summary, this study has provided a systematic view and significantly advanced our understanding of the miRNA functions in mosquito reproduction and P. falciparum infection. PHD

The development of techniques for measuring gene expression globally has greatly expanded our understanding of gene regulatory mechanisms in depth and scale. We can now quantify every intermediate and transition in the canonical pathway of gene expression-from DNA to mRNA to protein-genome-wide. Employing such measurements in parallel can produce rich datasets, but extracting the most information requires careful experimental design and analysis. Here, we argue for the value of genome-wide studies that measure multiple outputs of gene expression over many timepoints during the course of a natural developmental process. We discuss our findings from a highly parallel gene expression dataset of meiotic differentiation, and those of others, to illustrate how leveraging these features can provide new and surprising insight into fundamental mechanisms of gene regulation.

The environment has profound effects on the expression of many traits and reaction norms describe the expression dynamics of a trait across a broad range of environmental conditions. Here, we analyze gene expression in Drosophila melanogaster across four different developmental temperatures (13–29 °C). Gene expression is highly plastic with 83.3% of the genes being differentially expressed. We distinguished three components of plasticity: 1) Dynamics of gene expression intensity (sum of change), 2) direction of change, and 3) curvature of the reaction norm (linear vs. quadratic). Studying their regulatory architecture we found that all three plasticity components were most strongly affected by the number of different transcription factors (TFs) binding to the target gene. More TFs were found in genes with less expression changes across temperatures. Although the effect of microRNAs was weaker, we consistently noted a trend in the opposite direction. The most plastic genes were regulated by fewer TFs and more microRNAs than less plastic genes. Different patterns of plasticity were also reflected by their functional characterization based on gene ontology. Our results suggest that reaction norms provide an important key to understand the functional requirements of natural populations exposed to variable environmental conditions.

Background Oleaginous microalgae are promising production hosts for the sustainable generation of lipid-based bioproducts and as bioenergy carriers such as biodiesel. Transcriptomics of the lipid accumulation phase, triggered efficiently by nitrogen starvation, is a valuable approach for the identification of gene targets for metabolic engineering. Results An explorative analysis of the detailed transcriptional response to different stages of nitrogen availability was performed in the oleaginous green alga Monoraphidium neglectum. Transcript data were correlated with metabolic data for cellular contents of starch and of different lipid fractions. A pronounced transcriptional down-regulation of photosynthesis became apparent in response to nitrogen starvation, whereas glucose catabolism was found to be up-regulated. An in-depth reconstruction and analysis of the pathways for glycerolipid, central carbon, and starch metabolism revealed that distinct transcriptional changes were generally found only for specific steps within a metabolic pathway. In addition to pathway analyses, the transcript data were also used to refine the current genome annotation. The transcriptome data were integrated into a database and complemented with data for other microalgae which were also subjected to nitrogen starvation. It is available at https://tdbmn.cebitec.uni-bielefeld.de. Conclusions Based on the transcriptional responses to different stages of nitrogen availability, a model for triacylglycerol and lipid hyperaccumulation is proposed, which involves transcriptional induction of thioesterases, differential regulation of lipases, and a re-routing of the central carbon metabolism. Over-expression of distinct thioesterases was identified to be a potential strategy to increase the oleaginous phenotype of M. neglectum, and furthermore specific lipases were identified as potential targets for future metabolic engineering approaches. Electronic supplementary material The online version of this article (doi:10.1186/s13068-017-0882-1) contains supplementary material, which is available to authorized users.

Cultivated tomato (Solanum lycopersicum L.) is susceptible to abiotic stresses, including drought and chilling stress, while its wild relative (Solanum habrochaites) exhibits tolerance to many abiotic stresses. Chilling roots to 6°C induces rapid-onset water stress by impeding water movement from roots to shoots. Wild S. habrochaites responds to root chilling by closing stomata and maintaining shoot turgor, while cultivated tomato fails to close stomata and wilts. This phenotypic response (shoot turgor maintenance under root chilling) is controlled by a major QTL stm9 on chromosome 9 from S. habrochaites that was previously high-resolution mapped to a 0.32 cM region, but its effects on transcriptional regulation were unknown. Here we used paired near isogenic lines (NILs) differing only for the presence or absence of the S. habrochaites introgression containing stm9 in an otherwise S. lycopersicum background to investigate global transcriptional regulation in response to rapid-onset water stress induced by root chilling. NIL175 contains the S. habrochaites introgression and exhibits tolerance to root chilling stress, while NIL163 does not contain the introgression and is susceptible. RNA from roots of the two NILs was obtained at five time points during exposure to root chilling and mRNA-Seq performed. Differential expression analysis and hierarchical clustering of transcript levels were used to determine patterns of and changes in mRNA levels. Our results show that the transcriptional response of roots exposed to chilling stress is complex, with both overlapping and unique responses in tolerant and susceptible lines. In general, susceptible NIL 163 had a more complex transcriptional response to root chilling, while NIL175 exhibited a more targeted response to the imposed stress. Our evidence suggests that both the tolerant and susceptible NILs may be primed for response to root-chilling, with many of these response genes located on chromosome 9. Furthermore, serine/threonine kinase activity likely has an important role in the root chilling response of tolerant NIL175.

CD24 is an anchored cell surface marker that is highly expressed in cancer cells (Lee et al., 2009) and its expression is associated with poorer outcome of cancer patients (Kristiansen et al., 2003). Phenotype comparison between two subpopulations derived from the Mvt1 cell line, CD24(-) cells (with no CD24 cell surface expression) and the CD24(+) cells, identified high tumorigenic capacity for the CD24(+) cells. In order to reveal the transcripts that support the CD24(+) aggressive and invasive phenotype we compared the gene profiles of these two subpopulations. mRNA profiles of CD24(-) and CD24(+) cells were generated by deep sequencing, in triplicate, using an Illumina HiSeq 2500. Here we provide a detailed description of the mRNA-seq analysis from our recent study (Rostoker et al., 2015). The mRNA-seq data have been deposited in the NCBI GEO database (accession number GSE68746).

Despite its eradication over 30 years ago, smallpox (as well as other orthopox viruses) remains a pathogen of interest both in terms of biodefense and for its use as a vector for vaccines and immunotherapies. Here we describe the application of mRNA-Seq transcriptome profiling to understanding immune responses in smallpox vaccine recipients. Contrary to other studies examining gene expression in virally infected cell lines, we utilized a mixed population of peripheral blood mononuclear cells in order to capture the essential intercellular interactions that occur in vivo, and would otherwise be lost, using single cell lines or isolated primary cell subsets. In this mixed cell population we were able to detect expression of all annotated vaccinia genes. On the host side, a number of genes encoding cytokines, chemokines, complement factors and intracellular signaling molecules were downregulated upon viral infection, whereas genes encoding histone proteins and the interferon response were upregulated. We also identified a small number of genes that exhibited significantly different expression profiles in subjects with robust humoral immunity compared with those with weaker humoral responses. Our results provide evidence that differential gene regulation patterns may be at work in individuals with robust humoral immunity compared with those with weaker humoral immune responses.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide-spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based Integrated Fluidic Circuit (IFC) that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to a variety of stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

Using mRNA-Seq and clinical data for 469 clear cell Renal Cell Carcinoma (ccRCC) samples from The Cancer Genome Atlas (TCGA), we develop a protocol to identify patients likely to have early recurrence of their disease. We first split the data into two sets, with 380 samples in the training set and 89 samples in the test set. Using the training set, we identify genes whose outlier status (high or low mRNA expression) is predictive of recurrence, based on Kaplan-Meier recurrence free survival log-rank p-value. We find a significant overlap among genes identified as predictive biomarkers in Reads per Kilobase Million (RPKM) normalized data and Raw Reads mRNA-Seq data. Using 80 consensus genes predictive in both RPKM and Raw Reads data, we define an outlier-based risk score R to stratify patients into two groups, a high-risk (early recurrence) group (R 2). The KM recurrence curve using this stratification shows excellent separation in training and test sets. Restricting the analysis to patients who had recurrence within two years (109 cases) and those who had no recurrence in five years (107 cases) we find that the risk predictor achieves ca. 80 percent sensitivity and specificity. The 80 genes identified by the outlier analysis were used to develop a more intuitive classifier based on Generalized Matrix Learning Vector Quantization (GMLVQ). This method stratifies samples into risk classes based on defining prototypes in feature space and an appropriate distance metric. GMLVQ identified a subset of 12 genes that have high accuracy in predicting recurrence, which suggests that an assay with a small number of genes might be able to predict recurrence in ccRCC.

Background Fish reared under intensive conditions are repeatedly exposed to stress, which negatively impacts growth. Although most fish follow a conserved pattern of stress response, with increased concentrations of cortisol, each species presents specificities in the cell response and stress tolerance. Therefore, culturing new species requires a detailed knowledge of these specific responses. The red cusk-eel (Genypterus chilensis) is a new economically important marine species for the Chilean aquaculture industry. However, there is no information on the stress- and cortisol-induced mechanisms that decrease skeletal muscle growth in this teleost. Results Using Illumina RNA-seq technology, skeletal muscle sequence reads for G. chilensis were generated under control and handling stress conditions. Reads were mapped onto a reference transcriptome, resulting in the in silico identification of 785 up-regulated and 167 down-regulated transcripts. Gene ontology enrichment analysis revealed a significant up-regulation of catabolic genes associated with skeletal muscle atrophy. These results were validated by RT-qPCR analysis for ten candidates genes involved in ubiquitin-mediated proteolysis, autophagy and skeletal muscle growth. Additionally, using a primary culture of fish skeletal muscle cells, the effect of cortisol was evaluated in relation to red cusk-eel skeletal muscle atrophy. Conclusions The present data demonstrated that handling stress promotes skeletal muscle atrophy in the marine teleost G. chilensis through the expression of components of the ubiquitin-proteasome and autophagy-lysosome systems. Furthermore, cortisol was a powerful inductor of skeletal muscle atrophy in fish myotubes. This study is an important step towards understanding the atrophy system in non-model teleost species and provides novel insights on the cellular and molecular mechanisms that control skeletal muscle growth in early vertebrates. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2232-7) contains supplementary material, which is available to authorized users.

Notch signaling as a conserved cell fate regulator is involved in the regulation of cell quiescence, proliferation, differentiation and postnatal tissue regeneration. However, how Notch signaling regulates porcine satellite cells (PSCs) has not been elucidated. We stably transfected Notch1 intracellular domain (N1ICD) into PSCs to analyze the gene expression profile and miRNA-seq. The analysis of the gene expression profile identified 295 differentially-expressed genes (DEGs) in proliferating-N1ICD PSCs (P-N1ICD) and nine DEGs on differentiating-N1ICD PSCs (D-N1ICD), compared with that in control groups (P-Control and D-Control, respectively). Analyzing the underlying function of DEGs showed that most of the upregulated DEGs enriched in P-N1ICD PSCs are related to the cell cycle. Forty-four and 12 known differentially-expressed miRNAs (DEMs) were identified in the P-N1ICD PSCs and D-N1ICD PSCs group, respectively. Furthermore, we constructed the gene-miRNA network of the DEGs and DEMs. In P-N1ICD PSCs, miR-125a, miR-125b, miR-10a-5p, ssc-miR-214, miR-423 and miR-149 are downregulated hub miRNAs, whose corresponding hub genes are marker of proliferation Ki-67 (MKI67) and nuclear receptor binding SET domain protein 2 (WHSC1). By contrast, miR-27a, miR-146a-5p and miR-221-3p are upregulated hub miRNAs, whose hub genes are RUNX1 translocation partner 1 (RUNX1T1) and fibroblast growth factor 2 (FGF2). All the hub miRNAs and genes are associated with cell proliferation. Quantitative RT-PCR results are consistent with the gene expression profile and miRNA-seq results. The results of our study provide valuable information for understanding the molecular mechanisms underlying Notch signaling in PSCs and skeletal muscle development.

Background Invasive aspergillosis is started after germination of Aspergillus fumigatus conidia that are inhaled by susceptible individuals. Fungal hyphae can grow in the lung through the epithelial tissue and disseminate hematogenously to invade into other organs. Low fungaemia indicates that fungal elements do not reside in the bloodstream for long. Results We analyzed whether blood represents a hostile environment to which the physiology of A. fumigatus has to adapt. An in vitro model of A. fumigatus infection was established by incubating mycelium in blood. Our model allowed to discern the changes of the gene expression profile of A. fumigatus at various stages of the infection. The majority of described virulence factors that are connected to pulmonary infections appeared not to be activated during the blood phase. Three active processes were identified that presumably help the fungus to survive the blood environment in an advanced phase of the infection: iron homeostasis, secondary metabolism, and the formation of detoxifying enzymes. Conclusions We propose that A. fumigatus is hardly able to propagate in blood. After an early stage of sensing the environment, virtually all uptake mechanisms and energy-consuming metabolic pathways are shut-down. The fungus appears to adapt by trans-differentiation into a resting mycelial stage. This might reflect the harsh conditions in blood where A. fumigatus cannot take up sufficient nutrients to establish self-defense mechanisms combined with significant growth. Electronic supplementary material The online version of this article (doi10.1186/s12864-015-1853-1) contains supplementary material, which is available to authorized users.

Background Invasive aspergillosis is started after germination of Aspergillus fumigatus conidia that are inhaled by susceptible individuals. Fungal hyphae can grow in the lung through the epithelial tissue and disseminate hematogenously to invade into other organs. Low fungaemia indicates that fungal elements do not reside in the bloodstream for long. Results We analyzed whether blood represents a hostile environment to which the physiology of A. fumigatus has to adapt. An in vitro model of A. fumigatus infection was established by incubating mycelium in blood. Our model allowed to discern the changes of the gene expression profile of A. fumigatus at various stages of the infection. The majority of described virulence factors that are connected to pulmonary infections appeared not to be activated during the blood phase. Three active processes were identified that presumably help the fungus to survive the blood environment in an advanced phase of the infection: iron homeostasis, secondary metabolism, and the formation of detoxifying enzymes. Conclusions We propose that A. fumigatus is hardly able to propagate in blood. After an early stage of sensing the environment, virtually all uptake mechanisms and energy-consuming metabolic pathways are shut-down. The fungus appears to adapt by trans-differentiation into a resting mycelial stage. This might reflect the harsh conditions in blood where A. fumigatus cannot take up sufficient nutrients to establish self-defense mechanisms combined with significant growth. peerReviewed

Accepted: 2013-12-31, 資料番号: SA1140056000

Background Neurodegenerative diseases are progressive and irreversible and they can be initiated by mutations in specific genes. Spalt-like genes (Sall) encode transcription factors expressed in the central nervous system. In humans, SALL mutations are associated with hereditary syndromes characterized by mental retardation, sensorineural deafness and motoneuron problems, among others. Drosophila sall mutants exhibit severe neurodegeneration of the central nervous system at embryonic stage 16, which surprisingly reverts later in development at embryonic stage 17, suggesting a potential to recover from neurodegeneration. We hypothesize that this recovery is mediated by a reorganization of the transcriptome counteracting SALL lost. To identify genes associated to neurodegeneration and neuroprotection, we used mRNA-Seq to compare the transcriptome of Drosophila sall mutant and wild type embryos from neurodegeneration and reversal stages. Results Neurodegeneration stage is associated with transcriptional changes in 220 genes, of which only 5% were already described as relevant for neurodegeneration. Genes related to the groups of Redox, Lifespan/Aging and Mitochondrial diseases are significantly represented at this stage. By contrast, neurodegeneration reversal stage is associated with significant changes in 480 genes, including 424 not previously associated with neuroprotection. Immune response and Salt stress are the most represented groups at this stage. Conclusions We identify new genes associated to neurodegeneration and neuroprotection by using an mRNA-Seq approach. The strong homology between Drosophila and human genes raises the possibility to unveil novel genes involved in neurodegeneration and neuroprotection also in humans.

Background The filamentous fungus Aspergillus fumigatus has become the most important airborne fungal pathogen causing life-threatening infections in immuno-compromised patients. Recently developed high-throughput transcriptome and proteome technologies, such as microarrays, RNA deep-sequencing, and LC-MS/MS of peptide mixtures, are of enormous value for systematically investigating pathogenic organisms. In the field of infection biology, one of the priorities is to collect and standardise data, in order to generate datasets that can be used to investigate and compare pathways and gene responses involved in pathogenicity. The “omics” era provides a multitude of inputs that need to be integrated and assessed. We therefore evaluated the potential of paired-end mRNA-Seq for investigating the regulatory role of the central mitogen activated protein kinase (MpkA). This kinase is involved in the cell wall integrity signalling pathway of A. fumigatus and essential for maintaining an intact cell wall in response to stress. Results The comparison of the transcriptome and proteome of an A. fumigatus wild-type strain with an mpkA null mutant strain revealed that 70.4% of the genome was found to be expressed and that MpkA plays a significant role in the regulation of many genes involved in cell wall remodelling, oxidative stress and iron starvation response, and secondary metabolite biosynthesis. Moreover, absence of the mpkA gene also strongly affects the expression of genes involved in primary metabolism. The data were further processed to evaluate the potential of the mRNA-Seq technique. We comprehensively matched up our data to published transcriptome studies and were able to show an improved data comparability of mRNA-Seq experiments independently of the technique used. Analysis of transcriptome and proteome data revealed only a weak correlation between mRNA and protein abundance. Conclusions High-throughput analysis of MpkA-dependent gene expression confirmed many previous findings that this kinase is important for regulating many genes involved in metabolic pathways. Our analysis showed more than 2000 differentially regulated genes. RNA deep-sequencing is less error-prone than established microarray-based technologies. It also provides additional information in A. fumigatus studies and as a result is more suitable for the creation of extensive datasets.

mRNA-Seq is a precise and highly reproducible technique for measurement of transcripts levels and yields sequence information of a transcriptome at a single nucleotide base-level thus enabling us to determine splice junctions and alternative splicing events with high confidence. Often analysis of mRNA-Seq data does not attempt to quantify the expressions at isoform level. In this paper our objective would be use the mRNA-Seq data to infer expression at isoform level, where splicing patterns of a gene is assumed to be known. A Bayesian latent variable based modeling framework is proposed here, where the parameterization enables us to infer at various levels. For example, expression variability of an isoform across different conditions; the model parameterization also allows us to carry out two-sample comparisons, e.g., using a Bayesian t-test, in addition simple presence or absence of an isoform can also be estimated by the use of the latent variables present in the model. In this paper we would carry out inference on isoform expression under different normalization techniques, since it has been recently shown that one of the most prominent sources of variation in differential call using mRNA-Seq data is the normalization method used. The statistical framework is developed for multiple isoforms and easily extends to reads mapping to multiple genes. This could be achieved by slight conceptual modifications in definitions of what we consider as a gene and what as an exon. Additionally proposed framework can be extended by appropriate modeling of the design matrix to infer about yet unknown novel transcripts. However such attempts should be made judiciously since the input date used in the proposed model does not use reads from splice junctions.

Background Transcriptome analysis is increasingly being used to study the evolutionary origins and ecology of non-model plants. One issue for both transcriptome assembly and differential gene expression analyses is the common occurrence in plants of hybridisation and whole genome duplication (WGD) and hybridization resulting in allopolyploidy. The divergence of duplicated genes following WGD creates near identical homeologues that can be problematic for de novo assembly and also reference based assembly protocols that use short reads (35 - 100 bp). Results Here we report a successful strategy for the assembly of two transcriptomes made using 75 bp Illumina reads from Pachycladon fastigiatum and Pachycladon cheesemanii. Both are allopolyploid plant species (2n = 20) that originated in the New Zealand Alps about 0.8 million years ago. In a systematic analysis of 19 different coverage cutoffs and 20 different k-mer sizes we showed that i) none of the genes could be assembled across all of the parameter space ii) assembly of each gene required an optimal set of parameter values and iii) these parameter values could be explained in part by different gene expression levels and different degrees of similarity between genes. Conclusions To obtain optimal transcriptome assemblies for allopolyploid plants, k-mer size and k-mer coverage need to be considered simultaneously across a broad parameter space. This is important for assembling a maximum number of full length ESTs and for avoiding chimeric assemblies of homeologous and paralogous gene copies.

To compare the two RNA-sequencing protocols, ribo-minus RNA-sequencing (rmRNA-seq) and polyA-selected RNA-sequencing (mRNA-seq), we acquired transcriptomic data–52 and 32 million alignable reads of 35 bases in length–from the mouse cerebrum, respectively. We found that a higher proportion, 44% and 25%, of the uniquely alignable rmRNA-seq reads, is in intergenic and intronic regions, respectively, as compared to 23% and 15% from the mRNA-seq dataset. Further analysis made an additional discovery of transcripts of protein-coding genes (such as Histone, Heg1, and Dux), ncRNAs, snoRNAs, snRNAs, and novel ncRNAs as well as repeat elements in rmRNA-seq dataset. This result suggests that rmRNA-seq method should detect more polyA- or bimorphic transcripts. Finally, through comparative analyses of gene expression profiles among multiple datasets, we demonstrated that different RNA sample preparations may result in significant variations in gene expression profiles.

The accumulating evidence of the beneficial effects of cinnamon (Cinnamomum burmanni) in type-2 diabetes, a chronic age-associated disease, has prompted the commercialisation of various supplemental forms of the spice. One such supplement, Cinnulin PF®, represents the water soluble fraction containing relatively high levels of the double-linked procyanidin type-A polymers of flavanoids. The overall aim of this study was to utilize genome-wide mRNA-Seq analysis to characterise the changes in gene expression caused by Cinnulin PF in immortalised human keratinocytes and microvascular endothelial cells, which are relevant with respect to diabetic complications. In summary, our findings provide insights into the mechanisms of action of Cinnulin PF in diabetes and diabetic complications. More generally, we identify relevant candidate genes which could provide the basis for further investigation. Keywords: cinnulin PF; diabetes; diabetic complications; diabetic ulcer; mRNA-Seq (Published: 23 February 2012) Citation: Pathobiology of Aging & Age-related Diseases 2012, 2 : 11905 - DOI: 10.3402/pba.v2i0.11905 To access the supplementary material to this article 'Supplementary tables 1-3' please see the Supplementary files in the column to the right (under Reading Tools)