Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources

Transcriptomic Analysis of Cyclamen persicum to Identify Invovled Genes in Triterpene Secondary Metabolites Pathway

Saponins are considered as a diverse group of natural active compounds, which are widely found in crops. Mevalonate pathway (MVA) is regarded as the main pathway for synthesis of saponins in crops. This study aims to compare transcriptome of the leaf with tuber of crop including tubers and roots. First, more than 166 million reads were generated. The existence of 36,678 unigenes in the two samples out of 48,936 assembled ones showed a significant difference in expression. Finally, 310 and 290 highly up-regulated genes in leaf and tuber were selected for the next analysis. In addition, the expression profiles of 13 key genes in the MVA pathway were compared in RNA sequencing and reverse transcription-quantitative polymerase chain reaction analysis. The results indicated that cyclamen tuber has a higher level of expression of MVA pathway genes. The topological analysis for gene co-expression network involved in triterpenoid synthesis represented that the genes at the beginning of such pathway play a critical role so that the reduction of their expression challenges triterpenoid synthesis severely. The tuber of the cyclamen appears to be the major site of triterpene synthesis, and transfer of excess Isopentenyl pyrophosphate (IPP) from tuber to leaf activates downstream genes in leaf of crop.

Transpulmonary Expression of Exosomal microRNAs in Idiopathic and Congenital Heart Disease-Related Pulmonary Arterial Hypertension

BACKGROUND

Pulmonary artery hypertension (PAH) is a fatal disease characterized by a complex pathogenesis. Exosomes containing microRNAs (miRs) have emerged as a novel biomarker. Transpulmonary exosomal miRs offer valuable insights into pulmonary circulation microenvironments. Hereby, we aimed to explore the potentials of transpulmonary exosomal miRs as differentiating factors between idiopathic PAH and congenital heart disease (CHD)-related PAH.

METHODS AND RESULTS

During right heart catheterization, we collected exosomes at pulmonary arteries in 25 patients diagnosed with idiopathic PAH and 20 patients with CHD-related PAH. Next-generation sequencing identified several candidate exosomal miRs. Using quantitative polymerase chain reaction, we validated the expressions of these miRs and revealed significantly elevated expressions of miR-21, miR-139-5p, miR-155-5p, let-7f-5p, miR-328-3p, miR-330-3p, and miR-103a-3p in patients with CHD-related PAH, in contrast to patients with idiopathic PAH. Among these miRs, miR-21 exhibited the highest expression in patients with CHD-related PAH. These findings were further corroborated in an external cohort comprising 10 patients with idiopathic PAH and 8 patients with CHD-related PAH. Using an in vitro flow model simulating the shear stress experienced by pulmonary endothelial cells, we observed a significant upregulation of miR-21. Suppressing miR-21 rescued the shear stress-induced downregulation of the RAS/phosphatidylinositol 3-kinase/protein kinase B pathway, leading to a mitigation of apoptosis.

CONCLUSIONS

Our study identified a pronounced expression of transpulmonary exosomal miR-21, particularly in patients with CHD-related PAH, through next-generation sequencing analysis. Further investigation is warranted to elucidate the regulatory mechanisms involving miR-21 in the pathophysiology of PAH.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

A high-throughput sequencing approach identifies immunotherapeutic targets for bacterial meningitis in neonates

BACKGROUND

Worldwide, Escherichia coli is the leading cause of neonatal Gram-negative bacterial meningitis, but full understanding of the pathogenesis of this disease is not yet achieved. Moreover, to date, no vaccine is available against bacterial neonatal meningitis.

METHODS

Here, we used Transposon Sequencing of saturated banks of mutants (TnSeq) to evaluate E. coli K1 genetic fitness in murine neonatal meningitis. We identified E. coli K1 genes encoding for factors important for systemic dissemination and brain infection, and focused on products with a likely outer-membrane or extra-cellular localization, as these are potential vaccine candidates. We used in vitro and in vivo models to study the efficacy of active and passive immunization.

RESULTS

We selected for further study the conserved surface polysaccharide Poly-β-(1-6)-N-Acetyl Glucosamine (PNAG), as a strong candidate for vaccine development. We found that PNAG was a virulence factor in our animal model. We showed that both passive and active immunization successfully prevented and/or treated meningitis caused by E. coli K1 in neonatal mice. We found an excellent opsonophagocytic killing activity of the antibodies to PNAG and in vitro these antibodies were also able to decrease binding, invasion and crossing of E. coli K1 through two blood brain barrier cell lines. Finally, to reinforce the potential of PNAG as a vaccine candidate in bacterial neonatal meningitis, we demonstrated that Group B Streptococcus, the main cause of neonatal meningitis in developed countries, also produced PNAG and that antibodies to PNAG could protect in vitro and in vivo against this major neonatal pathogen.

INTERPRETATION

Altogether, these results indicate the utility of a high-throughput DNA sequencing method to identify potential immunotherapy targets for a pathogen, including in this study a potential broad-spectrum target for prevention of neonatal bacterial infections.

FUNDINGS

ANR Seq-N-Vaq, Charles Hood Foundation, Hearst Foundation, and Groupe Pasteur Mutualité.

Comparative de novo transcriptome analysis of flower and root of Oliveria decumbens Vent. to identify putative genes in terpenes biosynthesis pathway

The Oliveria decumbens Vent. is a wild, rare, annual medicinal plant and endemic plant of Iran that has metabolites (mostly terpenes) which make it a precious plant in Persian Traditional Medicine and also a potential chemotherapeutic agent. The lack of genetic resources has slowed the discovery of genes involved in the terpenes biosynthesis pathway. It is a wild relative of Daucus carota. In this research, we performed the transcriptomic differences between two samples, flower and root of Oliveria decumbens, and also analyze the expression value of the genes involved in terpenoid biosynthesis by RNA-seq and its essential oil’s phytochemicals analyzed by GC/MS. In total, 136,031,188 reads from two samples of flower and root have been produced. The result shows that the MEP pathway is mostly active in the flower and the MVA in the root. Three genes of GPP, FPPS, and GGPP that are the precursors in the synthesis of mono, di, and triterpenes are upregulated in root and 23 key genes were identified that are involved in the biosynthesis of terpenes. Three genes had the highest upregulation in the root including, and on the other hand, another three genes had the expression only in the flower. Meanwhile, 191 and 185 upregulated genes in the flower and root of the plant, respectively, were selected for the gene ontology analysis and reconstruction of co-expression networks. The current research is the first of its kind on Oliveria decumbens transcriptome and discussed 67 genes that have been deposited into the NCBI database. Collectively, the information obtained in this study unveils the new insights into characterizing the genetic blueprint of Oliveria decumbens Vent. which paved the way for medical/plant biotechnology and the pharmaceutical industry in the future.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Transcriptome Analysis Reveals Potential Mechanisms for Ethylene-Inducible Pedicel–Fruit Abscission Zone Activation in Non-Climacteric Sweet Cherry (Prunus avium L.)

The harvesting of sweet cherry (Prunus avium L.) fruit is a labor-intensive process. The mechanical harvesting of sweet cherry fruit is feasible; however, it is dependent on the formation of an abscission zone at the fruit–pedicel junction. The natural propensity for pedicel—fruit abscission zone (PFAZ) activation varies by cultivar, and the general molecular basis for PFAZ activation is not well characterized. In this study, ethylene-inducible change in pedicel fruit retention force (PFRF) was recorded in a developmental time-course with a concomitant analysis of the PFAZ transcriptome from three sweet cherry cultivars. In ‘Skeena’, mean PFRF for both control and treatment fruit dropped below the 0.40 kg-force (3.92 N) threshold for mechanical harvesting, indicating the activation of a discrete PFAZ. In ‘Bing’, mean PFRF for both control and treatment groups decreased over time. However, a mean PFRF conducive to mechanical harvesting was achieved only in the ethylene-treated fruit. While in ‘Chelan’ the mean PFRF of the control and treatment groups did not meet the threshold required for efficient mechanical harvesting. Transcriptome analysis of the PFAZ region followed by the functional annotation, differential expression analysis, and gene ontology (GO) enrichment analyses of the data facilitated the identification of phytohormone-responsive and abscission-related transcripts, as well as processes that exhibited differential expression and enrichment in a cultivar-dependent manner over the developmental time-course. Additionally, read alignment-based variant calling revealed several short variants in differentially expressed genes, associated with enriched gene ontologies and associated metabolic processes, lending potential insight into the genetic basis for different abscission responses between the cultivars. These results provide genetic targets for the induction or inhibition of PFAZ activation, depending on the desire to harvest the fruit with or without the stem attached. Understanding the genetic mechanisms underlying the development of the PFAZ will inform future cultivar development while laying a foundation for mechanized sweet cherry harvest.

Trehalose-Induced Remodelling of the Human Microbiota Affects Clostridioides difficile Infection Outcome in an In Vitro Colonic Model: A Pilot Study

Within the human intestinal tract, dietary, microbial- and host-derived compounds are used as signals by many pathogenic organisms, including Clostridioides difficile. Trehalose has been reported to enhance virulence of certain C. difficile ribotypes; however, such variants are widespread and not correlated with clinical outcomes for patients suffering from C. difficile infection (CDI). Here, we make preliminary observations on how trehalose supplementation affects the microbiota in an in vitro model and show that trehalose-induced changes can reduce the outgrowth of C. difficile, preventing simulated CDI. Three clinically reflective human gut models simulated the effects of sugar (trehalose or glucose) or saline ingestion on the microbiota. Models were instilled with sugar or saline and further exposed to C. difficile spores. The recovery of the microbiota following antibiotic treatment and CDI induction was monitored in each model. The human microbiota remodelled to utilise the bioavailable trehalose. Clindamycin induction caused simulated CDI in models supplemented with either glucose or saline; however, trehalose supplementation did not result in CDI, although limited spore germination did occur. The absence of CDI in trehalose model was associated with enhanced abundances of Finegoldia, Faecalibacterium and Oscillospira, and reduced abundances of Klebsiella and Clostridium spp., compared with the other models. Functional analysis of the microbiota in the trehalose model revealed differences in the metabolic pathways, such as amino acid metabolism, which could be attributed to prevention of CDI. Our data show that trehalose supplementation remodelled the microbiota, which prevented simulated CDI, potentially due to enhanced recovery of nutritionally competitive microbiota against C. difficile.

Transcriptional and Metabolic Response of Wine-Related Lactiplantibacillus plantarum to Different Conditions of Aeration and Nitrogen Availability

Lactic acid bacteria (LAB) perform the process of malolactic fermentation (MLF) in wine. Availability of oxygen and nitrogen nutrients could influence LAB growth, malolactic activity, and other metabolic pathways, impacting the subsequent wine quality. The impact of these two factors has received limited investigation within LAB, especially on a transcriptome level. The aim of this study was to evaluate metabolic changes in the strain Lactiplantibacillus plantarum IWBT B063, growing in synthetic grape juice medium (GJM) under different oxygen exposure conditions, and with low availability of nitrogen-based nutrients. Next-generation sequencing was used to analyze expression across the transcriptome (RNA-seq), in combination with conventional microbiological and chemical analysis. L. plantarum consumed the malic acid present in all the conditions evaluated, with a slight delay and impaired growth for nitrogen limitation and for anaerobiosis. Comparison of L. plantarum transcriptome during growth in GJM with and without O2 revealed differential expression of 148 functionally annotated genes, which were mostly involved in carbohydrate metabolism, genetic information processing, and signaling and cellular processes. In particular, genes with a protective role against oxidative stress and genes related to amino acid metabolism were differentially expressed. This study confirms the suitability of L. plantarum IWBT B063 to carry out MLF in different environmental conditions due to its potential adaption to the stress conditions tested and provides a better understanding of the genetic background of an industrially relevant strain.

Unveiling the shade nature of cyanic leaves: A view from the "blue absorbing side" of anthocyanins

Anthocyanins have long been suggested as having great potential in offering photoprotection to plants facing high light irradiance. Nonetheless, their effective ability in protecting the photosynthetic apparatus from supernumerary photons has been questioned by some authors, based upon the inexact belief that anthocyanins almost exclusively absorb green photons, which are poorly absorbed by chlorophylls. Here we focus on the blue light absorbing features of anthocyanins, a neglected issue in anthocyanin research. Anthocyanins effectively absorb blue photons: the absorbance of blue relative to green photons increases from tri- to mono-hydroxy B-ring substituted structures, reaching up to 50% of green photons absorption. We offer a comprehensive picture of the molecular events activated by low blue-light availability, extending our previous analysis in purple and green basil, which we suggest to be responsible for the "shade syndrome" displayed by cyanic leaves. While purple leaves display overexpression of genes promoting chlorophyll biosynthesis and light harvesting, in green leaves it is the genes involved in the stability/repair of photosystems that are largely overexpressed. As a corollary, this adds further support to the view of an effective photoprotective role of anthocyanins. We discuss the profound morpho-anatomical adjustments imposed by the epidermal anthocyanin shield, which reflect adjustments in light harvesting capacity under imposed shade and make complex the analysis of the photosynthetic performance of cyanic versus acyanic leaves.

The androgen receptor is a tumor suppressor in estrogen receptor-positive breast cancer

The role of the androgen receptor (AR) in estrogen receptor (ER)-α-positive breast cancer is controversial, constraining implementation of AR-directed therapies. Using a diverse, clinically relevant panel of cell-line and patient-derived models, we demonstrate that AR activation, not suppression, exerts potent antitumor activity in multiple disease contexts, including resistance to standard-of-care ER and CDK4/6 inhibitors. Notably, AR agonists combined with standard-of-care agents enhanced therapeutic responses. Mechanistically, agonist activation of AR altered the genomic distribution of ER and essential co-activators (p300, SRC-3), resulting in repression of ER-regulated cell cycle genes and upregulation of AR target genes, including known tumor suppressors. A gene signature of AR activity positively predicted disease survival in multiple clinical ER-positive breast cancer cohorts. These findings provide unambiguous evidence that AR has a tumor suppressor role in ER-positive breast cancer and support AR agonism as the optimal AR-directed treatment strategy, revealing a rational therapeutic opportunity.

A piggybacking mechanism enables peroxisomal localization of the glyoxylate cycle enzyme Mdh2 in yeast

Eukaryotic cells have evolved organelles that allow the compartmentalization and regulation of metabolic processes. Knowledge of molecular mechanisms that allow temporal and spatial organization of enzymes within organelles is therefore crucial for understanding eukaryotic metabolism. Here, we show that the yeast malate dehydrogenase 2 (Mdh2) is dually localized to the cytosol and to peroxisomes and is targeted to peroxisomes via association with Mdh3 and a Pex5-dependent piggybacking mechanism. This dual localization of Mdh2 contributes to our understanding of the glyoxylate cycle and provides a new perspective on compartmentalization of cellular metabolism, which is critical for the perception of metabolic disorders and aging.

Overexpression of AtAHL20 causes delayed flowering in Arabidopsis via repression of FT expression

BACKGROUND

The 29-member Arabidopsis AHL gene family is classified into three main classes based on nucleotide and protein sequence evolutionary differences. These differences include the presence or absence of introns, type and/or number of conserved AT-hook and PPC domains. AHL gene family members are divided into two phylogenetic clades, Clade-A and Clade-B. A majority of the 29 members remain functionally uncharacterized. Furthermore, the biological significance of the DNA and peptide sequence diversity, observed in the conserved motifs and domains found in the different AHL types, is a subject area that remains largely unexplored.

RESULTS

Transgenic plants overexpressing AtAHL20 flowered later than the wild type under both short and long days. Transcript accumulation analyses showed that 35S:AtAHL20 plants contained reduced FT, TSF, AGL8 and SPL3 mRNA levels. Similarly, overexpression of AtAHL20's orthologue in Camelina sativa, Arabidopsis' closely related Brassicaceae family member species, conferred a late-flowering phenotype via suppression of CsFT expression. However, overexpression of an aberrant AtAHL20 gene harboring a missense mutation in the AT-hook domain's highly conserved R-G-R core motif abolished the late-flowering phenotype. Data from targeted yeast-two-hybrid assays showed that AtAHL20 interacted with itself and several other Clade-A Type-I AHLs which have been previously implicated in flowering-time regulation: AtAHL19, AtAHL22 and AtAHL29.

CONCLUSION

We showed via gain-of-function analysis that AtAHL20 is a negative regulator of FT expression, as well as other downstream flowering time regulating genes. A similar outcome in Camelina sativa transgenic plants overexpressing CsAHL20 suggest that this is a conserved function. Our results demonstrate that AtAHL20 acts as a photoperiod-independent negative regulator of transition to flowering.

Characterization of the leaf rust responsive ARF genes in wheat (Triticum aestivum L.)

Genome-wide identification, classification, functional characterization and expression analysis of Auxin Responsive Factor (ARF) gene family in wheat reveal their attributes and role during leaf rust infection. Auxins are important plant growth regulators that also impact plant-pathogen interaction. Auxin responsive factors (ARF) are plant specific transcription factors that control responses to auxins. Whole genome investigation of ARF gene family is limited in allohexaploid wheat (Triticum aestivum L.). Comprehensive study of this gene family was carried out by employing the currently available reference genome sequence of wheat. In total, 27 ARF genes were identified and located on the wheat genome as well as were positioned on wheat chromosome arms. Additionally, examination of the predicted genes unveiled a decent degree of relatedness within and among the phylogenetic clades. Leaf rust, caused by the obligate biotrophic fungal pathogen Puccinia triticina, is responsible for drastic loss of wheat crop worldwide reducing grain yield by 10-90%. Expression profiling of ARF genes in retort to leaf rust infection indicated their differential regulation during this plant-pathogen interaction. Highest expression of ARF genes were observed at 12 hpi that was maintained up to 72 hpi during incompatible interaction, whereas the high expression levels receded at 48 hpi during compatible interactions. Few of the identified ARF genes were likely to be post-transcriptionally regulated by microRNAs. Many light and stress responsive elements were detected in the promoter regions of ARF genes. Microsynteny analysis showed the conservation of ARF genes within the members of the Poaceae family. This study provides fundamental details for understanding the different types of ARF genes in wheat and there putative roles during leaf rust-wheat interaction.

Circulating MicroRNAs from Serum Exosomes May Serve as a Putative Biomarker in the Diagnosis and Treatment of Patients with Focal Cortical Dysplasia

Focal cortical dysplasia (FCD) is a congenital malformation of cortical development where the cortical neurons located in the brain area fail to migrate in the proper formation. Epilepsy, particularly medically refractory epilepsy, is the most common clinical presentation for all types of FCD. This study aimed to explore the expression change of circulating miRNAs in patients with FCD from serum exosomes. A total of nine patients with FCD and four healthy volunteers were enrolled in this study. The serum exosomes were isolated from the peripheral blood of the subjects. Transmission electron microscopy (TEM) was used to identify the exosomes. Both exosomal markers and neuronal markers were detected by Western blotting analysis to prove that we could obtain central nervous system-derived exosomes from the circulation. The expression profiles of circulating exosomal miRNAs were assessed using next-generation sequencing analysis (NGS). We obtained a total of 107 miRNAs with dominant fold change (>2-fold) from both the annotated 5p-arm and 3p-arm of 2780 mature miRNAs. Based on the integrated platform of HMDD v3.2, miRway DB and DIANA-miRPath v3.0 online tools, and confirmed by MiRBase analysis, four potentially predicted miRNAs from serum exosomes in patients with FCD were identified, including miR194-2-5p, miR15a-5p, miR-132-3p, and miR-145-5p. All four miRNAs presented upregulated expression in patients with FCD compared with controls. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and pathway category of four target miRNAs, we found eight possible signaling pathways that may be related to FCD. Among them, we suggest that the mTOR signaling pathway, PI3K-Akt signaling pathway, p53 signaling pathway, and cell cycle regulation and TGF-beta signaling pathway are high-risk pathways that play a crucial role in the pathogenesis of FCD and refractory epilepsy. Our results suggest that the circulating miRNAs from exosomes may provide a potential biomarker for diagnostic, prognostic, and therapeutic adjuncts in patients with FCD and refractory epilepsy.

Comparative transcriptome profiling reveals cold stress responsiveness in two contrasting Chinese jujube cultivars

BACKGROUND

Low temperature is a major factor influencing the growth and development of Chinese jujube (Ziziphus jujuba Mill.) in cold winter and spring. Little is known about the molecular mechanisms enabling jujube to cope with different freezing stress conditions. To elucidate the freezing-related molecular mechanism, we conducted comparative transcriptome analysis between 'Dongzao' (low freezing tolerance cultivar) and 'Jinsixiaozao' (high freezing tolerance cultivar) using RNA-Seq.

RESULTS

More than 20,000 genes were detected at chilling (4 °C) and freezing (- 10 °C, - 20 °C, - 30 °C and - 40 °C) stress between the two cultivars. The numbers of differentially expressed genes (DEGs) between the two cultivars were 1831, 2030, 1993, 1845 and 2137 under the five treatments. Functional enrichment analysis suggested that the metabolic pathway, response to stimulus and catalytic activity were significantly enriched under stronger freezing stress. Among the DEGs, nine participated in the Ca²⁺ signal pathway, thirty-two were identified to participate in sucrose metabolism, and others were identified to participate in the regulation of ROS, plant hormones and antifreeze proteins. In addition, important transcription factors (WRKY, AP2/ERF, NAC and bZIP) participating in freezing stress were activated under different degrees of freezing stress.

CONCLUSIONS

Our research first provides a more comprehensive understanding of DEGs involved in freezing stress at the transcriptome level in two Z. jujuba cultivars with different freezing tolerances. These results may help to elucidate the molecular mechanism of freezing tolerance in jujube and also provides new insights and candidate genes for genetically enhancing freezing stress tolerance.

The role of myoglobin in epithelial cancers: Insights from transcriptomics

The muscle-associated respiratory protein myoglobin (MB) is expressed in multiple types of cancer, including breast and prostate tumors. In Kaplan-Meier analyses of the two tumor types, MB positivity is associated with favorable prognoses. Despite its well-characterized function in myocytes, the role of MB in cancer remains unclear. To study the impact of endogenous MB expression, small interfering RNA MB-knockdown cells were engineered using breast, prostate and colon cancer cell lines (MDA-MB468, LNCaP, DLD-1), and their transcriptomes were investigated using RNA-Seq at different oxygen levels. In MB-positive cells, increased expression of glycolytic genes was observed, which was possibly mediated by a higher activity of hypoxia-inducible factor 1α. In addition, the results of the gene set enrichment analysis suggested that MB contributed to fatty acid transport and turnover. MB-positive, wild-type-p53 LNCaP cells also exhibited increased expression of p53 target genes involved in cell cycle checkpoint control and prevention of cell migration. MB-positive cells expressing mutant p53 exhibited upregulation of genes associated with prolonged cancer cell viability and motility. Therefore, it was hypothesized that these transcriptomic differences may result from MB-mediated generation of nitric oxide or reactive oxygen species, thus employing established enzymatic activities of the globin. In summary, the transcriptome comparisons identified potential molecular functions of MB in carcinogenesis by highlighting the interaction of MB with key metabolic and regulatory processes.

RNA-seq Analysis of Salt-Stressed Versus Non Salt-Stressed Transcriptomes of Chenopodium quinoa Landrace R49

Quinoa (Chenopodium quinoa Willd.), a model halophytic crop species, was used to shed light on salt tolerance mechanisms at the transcriptomic level. An RNA-sequencing analysis of genotype R49 at an early vegetative stage was performed by Illumina paired-ends method comparing high salinity and control conditions in a time-course pot experiment. Genome-wide transcriptional salt-induced changes and expression profiling of relevant salt-responsive genes in plants treated or not with 300 mM NaCl were analyzed after 1 h and 5 days. We obtained up to 49 million pairs of short reads with an average length of 101 bp, identifying a total of 2416 differentially expressed genes (DEGs) based on the treatment and time of sampling. In salt-treated vs. control plants, the total number of up-regulated and down-regulated genes was 945 and 1471, respectively. The number of DEGs was higher at 5 days than at 1 h after salt treatment, as reflected in the number of transcription factors, which increased with time. We report a strong transcriptional reprogramming of genes involved in biological processes like oxidation-reduction, response to stress and response to abscisic acid (ABA), and cell wall organization. Transcript analyses by real-time RT- qPCR supported the RNA-seq results and shed light on the contribution of roots and shoots to the overall transcriptional response. In addition, it revealed a time-dependent response in the expression of the analyzed DEGs, including a quick (within 1 h) response for some genes, suggesting a "stress-anticipatory preparedness" in this highly salt-tolerant genotype.

Transcriptomic Responses to Water Deficit and Nematode Infection in Mycorrhizal Tomato Roots

Climate changes include the intensification of drought in many parts of the world, increasing its frequency, severity and duration. However, under natural conditions, environmental stresses do not occur alone, and, in addition, more stressed plants may become more susceptible to attacks by pests and pathogens. Studies on the impact of the arbuscular mycorrhizal (AM) symbiosis on tomato response to water deficit showed that several drought-responsive genes are differentially regulated in AM-colonized tomato plants (roots and leaves) during water deficit. To date, global changes in mycorrhizal tomato root transcripts under water stress conditions have not been yet investigated. Here, changes in root transcriptome in the presence of an AM fungus, with or without water stress (WS) application, have been evaluated in a commercial tomato cultivar already investigated for the water stress response during AM symbiosis. Since root-knot nematodes (RKNs, Meloidogyne incognita) are obligate endoparasites and cause severe yield losses in tomato, the impact of the AM fungal colonization on RKN infection at 7 days post-inoculation was also evaluated. Results offer new information about the response to AM symbiosis, highlighting a functional redundancy for several tomato gene families, as well as on the tomato and fungal genes involved in WS response during symbiosis, underlying the role of the AM fungus. Changes in the expression of tomato genes related to nematode infection during AM symbiosis highlight a role of AM colonization in triggering defense responses against RKN in tomato. Overall, new datasets on the tomato response to an abiotic and biotic stress during AM symbiosis have been obtained, providing useful data for further researches.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.

Dissecting Adaptation Mechanisms to Contrasting Solar Irradiance in the Mediterranean Shrub Cistus incanus

Molecular mechanisms that are the base of the strategies adopted by Mediterranean plants to cope with the challenges imposed by limited or excessive solar radiation during the summer season have received limited attention. In our study, conducted on C. incanus plants growing in the shade or in full sunlight, we performed measurements of relevant physiological traits, such as leaf water potential, gas exchange and PSII photochemistry, RNA-Seq with de-novo assembly, and the analysis of differentially expressed genes. We also identified and quantified photosynthetic pigments, abscisic acid, and flavonoids. Here, we show major mechanisms regulating light perception and signaling which, in turn, sustain the shade avoidance syndrome displayed by the 'sun loving' C. incanus. We offer clear evidence of the detrimental effects of excessive light on both the assembly and the stability of PSII, and the activation of a suite of both repair and effective antioxidant mechanisms in sun-adapted leaves. For instance, our study supports the view of major antioxidant functions of zeaxanthin in sunny plants concomitantly challenged by severe drought stress. Finally, our study confirms the multiple functions served by flavonoids, both flavonols and flavanols, in the adaptive mechanisms of plants to the environmental pressures associated to Mediterranean climate.

Changes in the transcriptome profile of breast cancer cells grown as spheroids

BACKGROUND

Multicellular tumor spheroids mimic the functional organization of tumors in vivo, providing biological readouts that predict the behavior of cancer cells more accurately. The current study aimed to evaluate the transcriptome (mRNAs and long non-coding RNAs) of multicellular tumor spheroids from breast cancer cells.

METHODS

MCF-7 cell spheroids were used; the transcriptome was analyzed using RNAseq and RNA microarrays; the secretion of macrophage migration inhibitor (MIF), a cytokine exported by the cholesterol efflux regulatory protein, was measured by ELISA. Linc00052 was inhibited using short-hairpin RNAs (shRNAs).

RESULTS

We found several differentially regulated mRNAs and lncRNAs in MCF-7 cell spheroids. We also found significant enrichment of the Wnt/B-catenin death receptor and the cholesterol metabolic processes. Interestingly, we also found an increased concentration of MIF. Further, at 12 and 20 days of 3D culture we found 221 and 1146 dysregulated lncRNAs, respectively; including linc00052 (long intergenic non-protein coding RNA 52), which has been involved in breast cancer. Linc00052 knock-down experiments suggest that it could be a key regulator of cholesterol pathways in breast cancer.

CONCLUSIONS

Our data shows that tumor spheroids can induce changes in the transcriptome of the cultured cells, including both mRNAs and ncRNA. One of the major changes included the deregulation of cholesterol pathways, of which linc00052 is apparently a key regulator.

Transcriptome profiling of induced susceptibility effects on soybean-soybean aphid (Hemiptera: Aphididae) interaction

OBJECTIVES

Soybean aphid (Aphis glycines Matsumura; SBA) is the most economically damaging insect of soybean (Glycine max) in the United States. One previous study demonstrated that avirulent (biotype 1) and virulent (biotype 2) biotypes could co-occur and interact on resistant (i.e., Rag1) and susceptible soybean resulting in induced susceptibility after 11 days of feeding. The main objective of this research was to employ RNA sequencing (RNA-seq) technique to compare the induced susceptibility effect of biotype 2 on susceptible and resistant soybean at day 1 and day 11 (i.e., both susceptible and resistant soybean were initially challenged by biotype 2 and the effect was monitored through biotype 1 populations).

DATA DESCRIPTION

We investigated susceptible and Rag1 transcriptome response to SBA feeding in soybean plants colonized by biotype 1 in the presence or absence of an inducer population (i.e., biotype 2). Ten RNA datasets are reported with 266,535,654 sequence reads (55.2 GB) obtained from pooled samples derived from the leaves collected at day 1 and day 11 post SBA infestation. A comprehensive understanding of these transcriptome data will enhance our understanding of interactions among soybean and two different biotypes of soybean aphids at the molecular level.

Gibberellin causes wide transcriptional modifications in the early stage of grape cluster development

Yaghooti grape of Sistan is seedless, early ripening but has compact clusters. To study gibberellin effect on cluster compactness of Yaghooti grape, it has been studied transcriptomic changes in three developmental stages (cluster formation, berry formation and final size of cluster). We found out that 5409 of 22,756 genes in cluster tissue showed significant changes under gibberellin. Finally, it was showed that 2855, 2862 and 497 genes have critically important role on above developmental stages, respectively. GO enrichment analysis showed that gibberellin enhances biochemical pathways activity. Moreover, genes involved in ribosomal structure and photosynthesis rate in cluster tissue were up- and down- regulated, respectively. In addition, we observed location of metabolomic activities was transferred from nucleus to cytoplasm and from cytoplasm to cell wall and intercellular spaces during cluster development; but there is not such situation in gibberellin treated samples.

OMICfpp: a fuzzy approach for paired RNA-Seq counts

Background

RNA sequencing is a widely used technology for differential expression analysis. However, the RNA-Seq do not provide accurate absolute measurements and the results can be different for each pipeline used. The major problem in statistical analysis of RNA-Seq and in the omics data in general, is the small sample size with respect to the large number of variables. In addition, experimental design must be taken into account and few tools consider it.

Results

We propose OMICfpp, a method for the statistical analysis of RNA-Seq paired design data. First, we obtain a p-value for each case-control pair using a binomial test. These p-values are aggregated using an ordered weighted average (OWA) with a given orness previously chosen. The aggregated p-value from the original data is compared with the aggregated p-value obtained using the same method applied to random pairs. These new pairs are generated using between-pairs and complete randomization distributions. This randomization p-value is used as a raw p-value to test the differential expression of each gene. The OMICfpp method is evaluated using public data sets of 68 sample pairs from patients with colorectal cancer. We validate our results through bibliographic search of the reported genes and using simulated data set. Furthermore, we compared our results with those obtained by the methods edgeR and DESeq2 for paired samples. Finally, we propose new target genes to validate these as gene expression signatures in colorectal cancer. OMICfpp is available at http://www.uv.es/ayala/software/OMICfpp_0.2.tar.gz.

Conclusions

Our study shows that OMICfpp is an accurate method for differential expression analysis in RNA-Seq data with paired design. In addition, we propose the use of randomized p-values pattern graphic as a powerful and robust method to select the target genes for experimental validation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5496-5) contains supplementary material, which is available to authorized users.

Transcriptome datasets from leaves and fruits of the sweet cherry cultivars ‘Bing’, ‘Lapins’ and ‘Rainier’

Sweet cherry fruits from different cultivars have different pre- and post-harvest qualities. Here we present the transcriptome profile datasets of leaves and mature fruits of three sweet cherry cultivars (‘Bing’, ‘Lapin’ and ‘Rainier’). Using 454 GS-FLX technology (454 Life Sciences, Roche), transcriptomes of leaves and mature fruits were obtained from these cultivars. These transcriptome data sets are reported here.

Glycerol as a substrate for Saccharomyces cerevisiae based bioprocesses - Knowledge gaps regarding the central carbon catabolism of this 'non-fermentable' carbon source

Glycerol is an interesting alternative carbon source in industrial bioprocesses due to its higher degree of reduction per carbon atom compared to sugars. During the last few years, significant progress has been made in improving the well-known industrial platform organism Saccharomyces cerevisiae with regard to its glycerol utilization capability, particularly in synthetic medium. This provided a basis for future metabolic engineering focusing on the production of valuable chemicals from glycerol. However, profound knowledge about the central carbon catabolism in synthetic glycerol medium is a prerequisite for such incentives. As a matter of fact, the current assumptions about the actual in vivo fluxes active on glycerol as the sole carbon source have mainly been based on omics data collected in complex media or were even deduced from studies with other non-fermentable carbon sources, such as ethanol or acetate. A number of uncertainties have been identified which particularly regard the role of the glyoxylate cycle, the subcellular localization of the respective enzymes, the contributions of mitochondrial transporters and the active anaplerotic reactions under these conditions. The review scrutinizes the current knowledge, highlights the necessity to collect novel experimental data using cells growing in synthetic glycerol medium and summarizes the current state of the art with regard to the production of valuable fermentation products from a carbon source that has been considered so far as 'non-fermentable' for the yeast S. cerevisiae.

Database on spermatozoa transcriptogram of catagorised Frieswal crossbred (Holstein Friesian X Sahiwal) bulls

Bull spermatozoa contain different functional genes and many of them plays important roles in different stages of spermatogenesis, spermatozoa kinetics, fertilization as well as embryonic development. RNA deep sequencing is one of the preferred tools for absolute quantification of messenger RNA. The intention of the current study was to investigate the abundance of spermatozoal transcripts in categorized Frieswal (Holstein-Friesian X Sahiwal) crossbred bull semen through RNA deep sequencing. A total 1546561 and 1019308 numbers of reads were identified among good and poor quality bull spermatozoa based on their conception rate. Post mapping with Bos taurus reference genome identified 1,321,236 and 842,022 number of transcripts among good and poor quality RNA libraries, respectively. However, a total number of 3510 and 6759 functional transcripts were identified among good and poor quality bull spermatozoa, respectively. Most of the identified transcripts were related to spermatozoa functions, embryonic development and other functional aspects of fertilization. Wet laboratory validation of the top five selected transcripts (AKAP4, PRM1, ATP2B4, TRIM71 and SLC9B2) illustrated the significant (p < 0.01) level of expression in the good quality crossbred bull semen than the poor quality counterparts. The present study with comprehensive profiling of spermatozoal transcripts provides a useful non-invasive tool to understand the causes of as well as an effective way to predict male infertility in crossbred bulls.

MiR-203a is differentially expressed during branching morphogenesis and EMT in breast progenitor cells and is a repressor of peroxidasin

MicroRNAs regulate developmental events such as branching morphogenesis, epithelial to mesenchymal transition (EMT) and its reverse process mesenchymal to epithelial transition (MET). In this study, we performed small RNA sequencing of a breast epithelial progenitor cell line (D492), and its mesenchymal derivative (D492M) cultured in three-dimensional microenvironment. Among the most downregulated miRNAs in D492M was miR-203a, a miRNA that plays an important role in epithelial differentiation. Increased expression of miR-203a was seen in D492, concomitant with increased complexity of branching. When miR-203a was overexpressed in D492M, a partial reversion towards epithelial phenotype was seen. Gene expression analysis of D492M and D492MmiR-203a revealed peroxidasin, a collagen IV cross-linker, as the most significantly downregulated gene in D492MmiR-203a. Collectively, we demonstrate that miR-203a expression temporally correlates with branching morphogenesis and is suppressed in D492M. Overexpression of miR-203a in D492M induces a partial MET and reduces the expression of peroxidasin. Furthermore, we demonstrate that miR-203a is a novel repressor of peroxidasin. MiR-203-peroxidasin axis may be an important regulator in branching morphogenesis, EMT/MET and basement membrane remodeling.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation

Background

Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes “walks” to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life.

Results

Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1–2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish.

Conclusions

This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5355-9) contains supplementary material, which is available to authorized users.

Symbiotic microbes affect the expression of male reproductive genes in Glossina m. morsitans

Background

Tsetse flies (Diptera, Glossinidae) display unique reproductive biology traits. Females reproduce through adenotrophic viviparity, nourishing the growing larva into their modified uterus until parturition. Males transfer their sperm and seminal fluid, produced by both testes and male accessory glands, in a spermatophore capsule transiently formed within the female reproductive tract upon mating. Both sexes are obligate blood feeders and have evolved tight relationships with endosymbionts, already shown to provide essential nutrients lacking in their diet. However, the partnership between tsetse and its symbionts has so far been investigated, at the molecular, genomic and metabolomics level, only in females, whereas the roles of microbiota in male reproduction are still unexplored.

Results

Here we begin unravelling the impact of microbiota on Glossina m. morsitans (G. morsitans) male reproductive biology by generating transcriptomes from the reproductive tissues of males deprived of their endosymbionts (aposymbiotic) via maternal antibiotic treatment and dietary supplementation. We then compared the transcriptional profiles of genes expressed in the male reproductive tract of normal and these aposymbiotic flies. We showed that microbiota removal impacts several male reproductive genes by depressing the activity of genes in the male accessory glands (MAGs), including sequences encoding seminal fluid proteins, and increasing expression of genes in the testes. In the MAGs, in particular, the expression of genes related to mating, immunity and seminal fluid components’ synthesis is reduced. In the testes, the absence of symbionts activates genes involved in the metabolic apparatus at the basis of male reproduction, including sperm production, motility and function.

Conclusions

Our findings mirrored the complementary roles male accessory glands and testes play in supporting male reproduction and open new avenues for disentangling the interplay between male insects and endosymbionts. From an applied perspective, unravelling the metabolic and functional relationships between tsetse symbionts and male reproductive physiology will provide fundamental information useful to understanding the biology underlying improved male reproductive success in tsetse. This information is of particular importance in the context of tsetse population control via Sterile Insect Technique (SIT) and its impact on trypanosomiasis transmission.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1289-2) contains supplementary material, which is available to authorized users.

Characterization of Dendritic Cell-Derived Extracellular Vesicles During Dengue Virus Infection

The dengue virus (DENV), transmitted by Aedes spp. mosquitoes, is one of the most important arboviral infections in the world. Dengue begins as a febrile condition, and in certain patients, it can evolve severe clinical outcomes, such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). The reasons why certain patients develop DHF or DSS have not been thoroughly elucidated to date, and both patient and viral factors have been implicated. Previous work has shown that a severe immune dysfunction involving dendritic cells and T cells plays a key role in increasing the disease severity, especially in secondary heterologous infections. Extracellular vesicles (EVs) are membranous particles that are secreted by several cell types involved in homeostatic and pathological processes. Secretion of EVs by infected cells can enhance immune responses or favor viral evasion. In this study, we compare the molecular content of EVs that are secreted by human primary dendritic cells under different conditions: uninfected or infected with DENV3 strains isolated from patients with different infection phenotypes (a severe case involving DSS and a mild case). Human monocyte-derived dendritic cells (mdDCs) were infected with the dengue virus strains DENV3 5532 (severe) or DENV3 290 (mild), and the EVs were isolated. The presence of cup-shaped EVs was confirmed by electron microscopy and immunostaining with CD9, CD81, and CD83. The RNA content from the mdDC-infected cells contained several mRNAs and miRNAs related to immune responses compared to the EVs from mock-infected mdDCs. A number of these RNAs were detected exclusively during infection with DENV3 290 or DENV3 5532. This result suggests that the differential immune modulation of mdDCs by dengue strains can be achieved through the EV pathway. Additionally, we observed an association of EVs with DENV-infectious particles that seem to be protected from antibodies targeting the DENV envelope protein. We also showed that EVs derived from cells treated with IFN alpha have a protective effect against DENV infection in other cells. These results suggested that during DENV infection, the EV pathway could be exploited to favor viral viability, although immune mechanisms to counteract viral infection can also involve DC-derived EVs.

Unraveling the transcriptional complexity of compactness in sistan grape cluster

Yaghooti grape of Sistan is the earliest ripening grape in Iran, harvested every May annually. It is adapted to dry conditions in Sistan region and its water requirement is less than the other grape cultivars. The transcriptional complexity of this grape was studied in three stages of cluster development. Totally, 24121 genes were expressed in different cluster development steps (step 1: cluster formation, step 2: berry formation, step 3: final size of cluster) of which 3040 genes in the first stage, 2381 genes in the second stage and 2400 genes in the third stage showed a significant increase in expression. GO analysis showed that when the clusters are ripening, the activity of the nucleus, cytoplasmic, cytosol, membrane and chloroplast genes in the cluster architecture cells decreases. In contrast, the activity of the endoplasmic reticulum, vacuole and extracellular region genes enhances. When Yaghooti grape is growing and developing, some of metabolic pathways were activated in the response to biotic and abiotic stresses. Gene co-expression network reconstruction showed that AGAMOUS is a key gene in compactness of Sistan grape cluster, because it influences on expression of GA gene which leads to increase cluster length and berries size.

Basal polarization of the immune responses to Streptococcus agalactiae susceptible and resistant tilapia (Oreochromis niloticus)

One of the highest priority areas for improvement is the development of effective strategies for decreasing disease mortality levels in aquaculture production, a better understanding of the components of the fish immune system and their functions in the context of pathogen invasion is needed. Tilapia is the most common fish in South China, and Streptococcus agalactiae has become the most serious disease problem for tilapia industry in China. Here, we profiled gene expression differences between tilapia differing in their susceptibility to S. agalactiae both basally (before infection) and at three early timepoints post-infection (5 h, 50 h, and 7 d). Between group comparisons revealed 5756 unique genes differentially expressed greater than 2-fold at one or more timepoints. And the resistant fish showed much more strong ability in pathogen recognition, antigen presentation, immune activation, while the susceptible fish showed fast activation of apoptosis. Taken together, the immune profiles expand our knowledge for molecular mechanisms for disease resistance, as well as provide solid molecular resources for further identification of the candidate markers for disease-resistant selection and evaluation of disease prevention and treatment options for tilapia industry.

TRAIN (Transcription of Repeats Activates INterferon) in response to chromatin destabilization induced by small molecules in mammalian cells

Cellular responses to the loss of genomic stability are well-established, while how mammalian cells respond to chromatin destabilization is largely unknown. We previously found that DNA demethylation on p53-deficient background leads to transcription of repetitive heterochromatin elements, followed by an interferon response, a phenomenon we named TRAIN (Transcription of Repeats Activates INterferon). Here, we report that curaxin, an anticancer small molecule, destabilizing nucleosomes via disruption of histone/DNA interactions, also induces TRAIN. Furthermore, curaxin inhibits oncogene-induced transformation and tumor growth in mice in an interferon-dependent manner, suggesting that anticancer activity of curaxin, previously attributed to p53-activation and NF-kappaB-inhibition, may also involve induction of interferon response to epigenetic derepression of the cellular 'repeatome'. Moreover, we observed that another type of drugs decondensing chromatin, HDAC inhibitor, also induces TRAIN. Thus, we proposed that TRAIN may be one of the mechanisms ensuring epigenetic integrity of mammalian cells via elimination of cells with desilenced chromatin.

Tsetse fly (Glossina pallidipes) midgut responses to Trypanosoma brucei challenge

BACKGROUND

Tsetse flies (Glossina spp.) are the prominent vector of African trypanosome parasites (Trypanosoma spp.) in sub-Saharan Africa, and Glossina pallidipes is the most widely distributed species in Kenya. This species displays strong resistance to infection by parasites, which are typically eliminated in the midgut shortly after acquisition from the mammalian host. Although extensive molecular information on immunity for the related species Glossina morsitans morsitans exists, similar information is scarce for G. pallidipes.

METHODS

To determine temporal transcriptional responses of G. pallidipes to Trypanosoma brucei brucei challenge, we conducted Illumina based RNA-seq on midgut organ and carcass from teneral females G. pallidipes at 24 and 48 h post-challenge (hpc) with T. b. brucei relative to their respective controls that received normal blood meals (without the parasite). We used a suite of bioinformatics tools to determine differentially expressed and enriched transcripts between and among tissues, and to identify expanded transcripts in G. pallidipes relative to their orthologs G. m. morsitans.

RESULTS

Midgut transcripts induced at 24 hpc encoded proteins were associated with lipid remodelling, proteolysis, collagen metabolism, apoptosis, and cell growth. Midgut transcripts induced at 48 hpc encoded proteins linked to embryonic growth and development, serine endopeptidases and proteosomal degradation of the target protein, mRNA translation and neuronal development. Temporal expression of immune responsive transcripts at 48 relative to 24 hpc was pronounced, indicative of a gradual induction of host immune responses the following challenge. We also searched for G. m. morsitans orthologous groups that may have experienced expansions in the G. pallidipes genome. We identified ten expanded groups in G. pallidipes with putative immunity-related functions, which may play a role in the higher refractoriness exhibited by this species.

CONCLUSIONS

There appears to be a lack of strong immune responses elicited by gut epithelia of teneral adults. This in combination with a compromised peritrophic matrix at this stage during the initial phase of T. b. brucei challenge may facilitate the increased parasite infection establishment noted in teneral flies relative to older adults. Although teneral flies are more susceptible than older adults, the majority of tenerals are still able to eliminate parasite infections. Hence, robust responses elicited at a later time point, such as 72 hpc, may clear parasite infections from the majority of flies. The expanded G. m. morsitans orthologous groups in G. pallidipes may also be functionally associated with the enhanced refractoriness to trypanosome infections reported in G. pallidipes relative to G. m. morsitans.

Modifying SAMseq to account for asymmetry in the distribution of effect sizes when identifying differentially expressed genes

RNA-Seq is a developing technology for generating gene expression data by directly sequencing mRNA molecules in a sample. RNA-Seq data consist of counts of reads recorded to a particular gene that are often used to identify differentially expressed (DE) genes. A common statistical method used to analyze RNA-Seq data is Significance Analysis of Microarray with emphasis on RNA-Seq data (SAMseq). SAMseq is a nonparametric method that uses a resampling technique to account for differences in sequencing depths when identifying DE genes. We propose a modification of this method that takes into account asymmetry in the distribution of the effect sizes by taking into account the sign of the test statistics. Through simulation studies, we showthat the proposed method, comparedwith the traditional SAMseqmethod and other existing methods provides better power for identifying truly DE genes or more sufficiently controls FDR in most settings where asymmetry is present. We illustrate the use of the proposed method by analyzing an RNA-Seq data set containing C57BL/6J (B6) and DBA/2J (D2) mouse strains samples.

Identification of tomato miRNAs responsive to root colonization by endophytic Pochonia chlamydosporia

The molecular mechanisms active during the endophytic phase of the fungus Pochonia chlamydosporia are still poorly understood. In particular, few data are available on the links between the endophyte and the root response, as modulated by noncoding small RNAs. In this study, we describe the microRNAs (miRNAs) that are differentially expressed (DE) in the roots of tomato, colonized by P. chlamydosporia. A genome-wide NGS expression profiling of small RNAs in roots, either colonized or not by the fungus, showed 26 miRNAs upregulated in inoculated roots. Their predicted target genes are involved in the plant information processing system, which recognizes, percepts, and transmits signals, with higher representations in processes such as apoptosis and plant defense regulation. RNAseq data showed that predicted miRNA target genes were downregulated in tomato roots after 4, 7, 10, and 21 days post P. chlamydosporia inoculation. The differential expression of four miRNAs was further validated using qPCR analysis. The P. chlamydosporia endophytic lifestyle in tomato roots included an intricate network of miRNAs and targets. Data provide a first platform of DE tomato miRNAs after P. chlamydosporia colonization. They indicated that several miRNAs are involved in the host response to the fungus, playing important roles for its recognition as a symbiotic microorganism, allowing endophytism by modulating the host defense reaction. Data also indicated that endophytism affects tRNA fragmentation. This is the first study on miRNAs induced by P. chlamydosporia endophytism and related development regulation effects in Solanum lycopersicum.

RocA Is an Accessory Protein to the Virulence-Regulating CovRS Two-Component System in Group A Streptococcus

Regulating gene expression during infection is critical to the ability of pathogens to circumvent the immune response and cause disease. This is true for the group A Streptococcus (GAS), a pathogen that causes both invasive (e.g., necrotizing fasciitis) and noninvasive (e.g., pharyngitis) diseases. The control of virulence (CovRS) two-component system has a major role in regulating GAS virulence factor expression. The regulator of cov (RocA) protein, which is a predicted kinase, functions in an undetermined manner through CovRS to alter gene expression and reduce invasive disease virulence. Here, we show that the ectopic expression of a truncated RocA derivative, harboring the membrane-spanning domains but not the dimerization or HATPase domain, is sufficient to complement a rocA mutant strain. Coupled with a previous bioinformatic study, the data are consistent with RocA being a pseudokinase. RocA reduces the ability of serotype M1 GAS isolates to express capsule and to evade killing in human blood, phenotypes that are not observed for M3 or M18 GAS due to isolates of these serotypes naturally harboring mutant rocA alleles. In addition, we found that varying the RocA concentration attenuates the regulatory activity of Mg²⁺ and the antimicrobial peptide LL-37, which positively and negatively regulate CovS function, respectively. Thus, we propose that RocA is an accessory protein to the CovRS system that influences the ability of GAS to modulate gene expression in response to host factors. A model of how RocA interacts with CovRS, and of the regulatory consequences of such activity, is presented.

Optimized reduced representation bisulfite sequencing reveals tissue-specific mCHH islands in maize

Background

DNA methylation plays important roles in many regulatory processes in plants. It is economically infeasible to profile genome-wide DNA methylation at a single-base resolution in maize, given its genome size of ~2.5 Gb. As an alternative, we adapted region of interest (ROI)-directed reduced representation bisulfite sequencing (RRBS) to survey genome-wide methylation in maize.

Results

We developed a pipeline for selecting restriction enzymes in silico and experimentally showed that, in the maize genome, MseI- and CviQI-digested fragments are precisely enriched in promoters and gene bodies, respectively. We proceeded with comparisons of epigenomes and transcriptomes between shoots and tassels and found that the occurrences of highly methylated, tissue-specific, mCHH islands upstream of transcription start sites (TSSs) were positively correlated with differential gene expression. Furthermore, 5′ regulatory regions between TSS and mCHH islands often contain putative binding sites of known transcription factors (TFs) that regulate the flowering process and the timing of the transition from the vegetative to the reproductive phase. By integrating MNase-seq and siRNA-seq data, we found that regions of mCHH islands accumulate 21nt-siRNAs in a tissue-specific manner, marking the transition to open chromatin, thereby ensuring the accessibility of TFs for tissue-specific gene regulation.

Conclusions

Our ROI-directed RRBS pipeline is eminently applicable to DNA methylation profiling of large genomes. Our results provide novel insights into the tissue-specific epigenomic landscapes in maize, demonstrating that DNA methylation and siRNA and chromatin accessibility constitute a critical, interdependent component that orchestrates the transition from the vegetative to the reproductive phase.

Electronic supplementary material

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

The purplish bifurcate mussel Mytilisepta virgata gene expression atlas reveals a remarkable tissue functional specialization

BACKGROUND

Mytilisepta virgata is a marine mussel commonly found along the coasts of Japan. Although this species has been the subject of occasional studies concerning its ecological role, growth and reproduction, it has been so far almost completely neglected from a genetic and molecular point of view. In the present study we present a high quality de novo assembled transcriptome of the Japanese purplish mussel, which represents the first publicly available collection of expressed sequences for this species.

RESULTS

The assembled transcriptome comprises almost 50,000 contigs, with a N50 statistics of ~1 kilobase and a high estimated completeness based on the rate of BUSCOs identified, standing as one of the most exhaustive sequence resources available for mytiloid bivalves to date. Overall this data, accompanied by gene expression profiles from gills, digestive gland, mantle rim, foot and posterior adductor muscle, presents an accurate snapshot of the great functional specialization of these five tissues in adult mussels.

CONCLUSIONS

We highlight that one of the most striking features of the M. virgata transcriptome is the high abundance and diversification of lectin-like transcripts, which pertain to different gene families and appear to be expressed in particular in the digestive gland and in the gills. Therefore, these two tissues might be selected as preferential targets for the isolation of molecules with interesting carbohydrate-binding properties. In addition, by molecular phylogenomics, we provide solid evidence in support of the classification of M. virgata within the Brachidontinae subfamily. This result is in agreement with the previously proposed hypothesis that the morphological features traditionally used to group Mytilisepta spp. and Septifer spp. within the same clade are inappropriate due to homoplasy.

The NCK and ABI adaptor genes in catfish and their involvement in ESC disease response

Adaptor proteins non-catalytic region of tyrosine kinase (NCK) and Abelson interactor (ABI) are crucial for disease response. NCK1 was identified to be a candidate gene for enteric septicemia of catfish (ESC) disease resistance, and was speculated to play similar roles during ESC and enteropathogenic Escherichia coli (EPEC) pathogenicity. ABI1 was reported as a positional candidate gene for bacterial cold water disease (BCWD) resistance in rainbow trout. In this study, three NCK genes and six ABI genes were identified in the channel catfish (Ictalurus punctatus) genome and blue catfish (I. furcatus) transcriptome, and annotated by domain structures, phylogenetic and syntenic analyses. Their expression patterns were examined in the intestine and liver of catfish after challenge with Edwardsiella ictaluri. In the intestine, NCK1, ABI2a, ABI2b, ABI3a were differentially expressed after E. ictaluri infection. In the liver, NCK2a, NCK2b, ABI1b, ABI2a, ABI2b were significantly upregulated in ESC susceptible fish. In general, the NCK and ABI genes, with exception of ABI3a gene and NCK1 gene, were expressed at higher levels in susceptible fish after infection than in control fish, but were expressed at lower levels in resistant fish than in the control fish. Taken together, these results support the notion that NCK and ABI genes are involved in disease processes facilitating pathogenesis of the E. ictaluri bacteria.

Genome-wide identification and characterization of NB-ARC resistant genes in wheat (Triticum aestivum L.) and their expression during leaf rust infection

NB-ARC domain-containing resistance genes from the wheat genome were identified, characterized and localized on chromosome arms that displayed differential yet positive response during incompatible and compatible leaf rust interactions. Wheat (Triticum aestivum L.) is an important cereal crop; however, its production is affected severely by numerous diseases including rusts. An efficient, cost-effective and ecologically viable approach to control pathogens is through host resistance. In wheat, high numbers of resistance loci are present but only few have been identified and cloned. A comprehensive analysis of the NB-ARC-containing genes in complete wheat genome was accomplished in this study. Complete NB-ARC encoding genes were mined from the Ensembl Plants database to predict 604 NB-ARC containing sequences using the HMM approach. Genome-wide analysis of orthologous clusters in the NB-ARC-containing sequences of wheat and other members of the Poaceae family revealed maximum homology with Oryza sativa indica and Brachypodium distachyon. The identification of overlap between orthologous clusters enabled the elucidation of the function and evolution of resistance proteins. The distributions of the NB-ARC domain-containing sequences were found to be balanced among the three wheat sub-genomes. Wheat chromosome arms 4AL and 7BL had the most NB-ARC domain-containing contigs. The spatio-temporal expression profiling studies exemplified the positive role of these genes in resistant and susceptible wheat plants during incompatible and compatible interaction in response to the leaf rust pathogen Puccinia triticina. Two NB-ARC domain-containing sequences were modelled in silico, cloned and sequenced to analyze their fine structures. The data obtained in this study will augment isolation, characterization and application NB-ARC resistance genes in marker-assisted selection based breeding programs for improving rust resistance in wheat.