Identification of reference genes for quantitative RT-PCR analysis of microRNAs and mRNAs in castor bean (Ricinus communis L.) under drought stress

Construct dysregulated miRNA-mRNA interaction networks to conjecture possible pathogenesis for Stomach adenocarcinomas

BACKGROUND

Post-transcriptional regulation of mRNA induced by microRNA is known crucial in tumor occurrence, progression, and metastasis. This study aims at identifying significant miRNA-mRNA axes for stomach adenocarcinomas (STAD).

METHOD

RNA expression profiles were collected from The Cancer Genome Atlas (TCGA) and GEO database for screening differently expressed RNAs and miRNAs (DE-miRNAs/DE-mRNAs). Functional enrichment analysis was conducted with Hiplot and DAVID-mirPath. Connectivity MAP was applied in compounds prediction. MiRNA-mRNA axes were forecasted by TarBase and MiRTarBase. Real-time reverse transcription polymerase chain reaction (RT-qPCR) of stomach specimen verified these miRNA-mRNA pairs. Diagnosis efficacy of miRNA-mRNA interactions was measured by Receiver operation characteristic curve and Decision Curve Analysis. Clinical and survival analysis were also carried out. CIBERSORT and ESTIMATE was employed for immune microenvironment measurement.

RESULT

Totally 228 DE-mRNAs (105 upregulated and 123 downregulated) and 38 DE-miRNAs (22 upregulated and 16 downregulated) were considered significant. TarBase and MiRTarBase identified 18 miRNA-mRNA pairs, 12 of which were verified in RT-qPCR. The network of miR-301a-3p/ELL2 and miR-1-3p/ANXA2 were established and verified in external validation. The model containing all 4 signatures showed better diagnosis ability. Via interacting with M0 macrophage and resting mast cell, these miRNA-mRNA axes may influence tumor microenvironment.

CONCLUSION

This study established a miRNA-mRNA network via bioinformatic analysis and experiment validation for STAD.

Identification of Reference Genes for RT-qPCR Analysis in Gleditsia microphylla under Abiotic Stress and Hormone Treatment

Gleditsia microphylla is an important galactomannan gums source plant with characteristics of drought resistance, barren tolerance, and good adaptability. However, the underlying molecular mechanisms of the biological process are not yet fully understood. Real-time quantitative PCR (RT-qPCR) is an accurate and convenient method to quantify the gene expression level and transcription abundance of suitable reference genes. This study aimed to screen the best internal reference genes in G. microphylla under abiotic stresses, hormone treatments, and different tissues. Based on the transcriptome data, twelve candidate reference genes were selected, and ultimately, nine of them were further evaluated by the geNorm, NormFinder, BestKeeper, and RefFinder algorithms. These results show that TATA-binding protein 1 (TBP1)and Eukaryotic translation initiation factor 4A1 (EIF4A1)were the two most stable reference genes, and glyceraldehyde-3-phosphate dehydrogenase A subunit, chloroplastic (GAPA)and glyceraldehyde-3-phosphate dehydrogenase B subunit, chloroplastic (GAPB)were the two most unstable reference genes across all samples under the given experimental conditions. Meanwhile, the most stable reference genes varied among the different groups and tissues. Therefore, this study suggests that it is better to use a specific reference gene for a particular case rather than using a common reference gene.

Lupeol Accumulation Correlates with Auxin in the Epidermis of Castor

Lupeol, a natural lupane-type pentacyclic triterpene, possesses various pharmacological properties, and its production attracts attention. Significant quantities of lupeol are deposited on the castor aerial organ surface and are easily extractable as a predominant wax constituent. Thus, castor might be considered as a potential bioreactor for the production of lupeol. The lupeol biosynthesis pathway is well known, but how it is regulated remains largely unknown. Among large numbers of castor cultivars, we targeted one accession line (337) with high levels of lupeol on its stem surface and low levels thereof on its hypocotyl surface, implicating that lupeol synthesis is differentially regulated in the two organs. To explore the underlying mechanisms, we did comparative transcriptome analysis of the first internode of 337 stem and the upper hypocotyl. Our results show that large amounts of auxin-related genes are differentially expressed in both parts, implying some possible interactions between auxin and lupeol production. We also found that several auxin-responsive cis-elements are present in promoter regions of HMGR and LUS genes encoding two key enzymes involved in lupeol production. Furthermore, auxin treatments apparently induced the expression levels of RcHMGR and RcLUS. Furthermore, we observed that auxin treatment significantly increased lupeol contents, whereas inhibiting auxin transport led to an opposite phenotype. Our study reveals some relationships between hormone activity and lupeol synthesis and might provide a promising way for improving lupeol yields in castor.

Cassava (Manihot esculenta) defensins: Prospection, structural analysis and tissue-specific expression under biotic/abiotic stresses

Defensins are a prominent family of antimicrobial peptides. They play sophisticated roles in the defense against pathogens in all living organisms, but few works address their expression under different conditions and plant tissues. The present work prospected defensins of Manihot esculenta Crantz, popularly known as cassava. Five defensin candidates (MeDefs) were retrieved from the genome sequences and characterized. Considering chromosome distribution, only MeDef1 and 2 occupy adjacent positions in the same chromosome arm. All 3D structures had antiparallel ß-sheets, an α-helix, and amphipathic residues distributed throughout the peptides with a predominance of cationic surface charge. MeDefs expression was validated by RT-qPCR, including two stress types (biotic: fungus Macrophomina pseudophaseolina, and abiotic: mechanical injury) and a combination of both stresses (fungus+injury) in three different tissues (root, stem, and leaf). For this purpose, ten reference genes (RGs) were tested, and three were chosen to characterize MeDef expression. MeDef3 was up-regulated at roots in all stress situations tested. MeDef1 and MeDef5 were induced in leaves under biotic and abiotic stresses, but not in both stress types simultaneously. Only MeDef2 was down-regulated in the stem tissue also with biotic/abiotic combined stresses. These results indicate that although defensins are known to be responsive to pathogen infection, they may act as preformed defense or, still, have tissue or stress specificities. Aspects of their structure, stability and evolution are also discussed.

Differential responses of thiol metabolism and genes involved in arsenic detoxification in tolerant and sensitive genotypes of bioenergy crop Ricinus communis

Castor, a non-food, dedicated bioenergy crop, has immense potential to be used for phytoremediation/revegetation of heavy metal contaminated sites. In the previous study, we identified arsenate [As(V)]-tolerant (WM) and As(V)-sensitive (GCH 2) genotypes of castor (Ricinus communis L.) with differential accumulation and tolerance of arsenic [As]. The role of thiols in As(V) toxicity and tolerance mechanism in the castor plant is not fully understood. On the one hand, thiol-dependent reduction of As(V) to As(III) by arsenate reductase (AR) makes it capable of reacting with thiol groups of protein leading to disturbed metabolic pathways; on the other hand, reduction of As(V) to arsenite [As(III)] by AR and then complexation of As(III) with phytochelatins (PCs) and compartmentalization of As(III)-PC complex are considered as the major detoxification mechanisms of As(V). In our study, the expression of RcAR increased in leaves and roots of As(V)-tolerant castor genotype WM but decreased in sensitive genotype GCH 2 due to 200 μM As(V) treatment. The activity of glutathione reductase (GR) increased significantly in the tolerant genotype, whereas it remained same in the sensitive genotype. GSH/GSSH ratio declined substantially in the sensitive genotype. The increased expression of phytochelatin synthase 1 isoform 1 (RcPCS1X1) in roots, RcPCS1X2 and metallothionein type 2 (RcMT2) in leaves, and c-type ABC transporter (RcABCC) in roots and leaves of WM was observed, but the expression of these genes declined or remained the same in GCH 2. Overall, our results suggest the essential roles of GR, RcAR, RcPCS1, RcMT2, and RcABCC in the tolerance of WM castor genotype to As(V) toxicity.

Identification of Appropriate Reference Genes for Normalizing miRNA Expression in Citrus Infected by Xanthomonas citri subsp. citri

MicroRNAs (miRNAs) are short noncoding RNA molecules that regulate gene expression at the posttranscriptional level. Reverse transcription-quantitative PCR (RT-qPCR) is one of the most common methods used for quantification of miRNA expression, and the levels of expression are normalized by comparing with reference genes. Thus, the selection of reference genes is critically important for accurate quantification. The present study was intended to identify appropriate miRNA reference genes for normalizing the level of miRNA expression in Citrus sinensis L. Osbeck and Citrus reticulata Blanco infected by Xanthomonas citri subsp. citri, which caused citrus canker disease. Five algorithms (Delta Ct, geNorm, NormFinder, BestKeeper and RefFinder) were used for screening reference genes, and two quantification approaches, poly(A) extension RT-qPCR and stem-loop RT-qPCR, were used to determine the most appropriate method for detecting expression patterns of miRNA. An overall comprehensive ranking output derived from the multi-algorithms showed that poly(A)-tailed miR162-3p/miR472 were the best reference gene combination for miRNA RT-qPCR normalization in citrus canker research. Candidate reference gene expression profiles determined by poly(A) RT-qPCR were more consistent in the two citrus species. To the best of our knowledge, this is the first systematic comparison of two miRNA quantification methods for evaluating reference genes. These results highlight the importance of rigorously assessing candidate reference genes and clarify some contradictory results in miRNA research on citrus.

Selection and validation of castor bean (Ricinus communis) reference genes for quantitative PCR (RT-qPCR) in developing and germinating seeds and expression pattern of four ricin-family genes

This study aims to expand the set of internal control genes used for RT-qPCR experiments with Castor bean (Ricinus communis) seeds by evaluating candidate genes across several seed tissues and developmental stages. Nine reference genes were selected, including actin-11 (ACT11), tubulin alpha-2 (Tα2), elongation factor 1-alpha (EF1-α), protein phosphatase 2A-2 (PP2A2), polyubiquitin-3 (PUB3) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Biological samples consisted of R. communis seeds in 15 stages of maturation and germination. We demonstrate that PP2A2, PUB3 and EF1-α are the most stably expressed genes across the tested conditions and therefore appropriate for RT-qPCR. Subsequently, those reference genes were used for the analysis of the expression of four R. communis ricin-family genes. In developing seeds, the highest ricin expression levels was seen in the nucellus and in the endosperm, whereas in germinating seeds a peak expression occurs 4-6 days after germination. The four tested ricin isoforms exhibited differential expression patterns across tissues and seed developmental stages, which may indicate distinct biological roles for each ricin gene.

Differential responses of growth, photosynthesis, oxidative stress, metals accumulation and NRAMP genes in contrasting Ricinus communis genotypes under arsenic stress

Effect of arsenate [As(V)] on biomass, photosynthetic rate, stomatal conductance, transpiration, oxidative stress, accumulation of As, Fe, Zn, Cu and Mn and expression of NRAMP genes was investigated in As(V) tolerant and sensitive genotypes of bioenergy crop Ricinus communis. As(V) treatments (100 and 200 μM) led to significant reduction in root and leaf biomass, photosynthetic rate, stomatal conductance and transpiration in GCH 2 and GCH 4 genotypes but no significant change or increase was observed in WM and DCH 177 genotypes. No significant difference was observed in hydrogen peroxide content and lipid peroxidation in As(V)-treated tolerant genotypes compared to control, whereas these parameters enhanced significantly in As(V)-treated sensitive genotypes. GCH 2 accumulated around two times As in leaves and showed significant reduction in concentration of Zn and Mn in the leaves and roots due to 200 μM As(V) treatment compared to WM. NRAMP genes are critical for uptake and distribution of essential divalent metal cations, photosynthesis and controlled production of reactive oxygen species in plants. RcNRAMP2, RcNRAMP3 and RcNRAMP5 genes showed differential expression in response to 200 μM As(V) in GCH 2 and WM suggesting that NRAMP genes are associated with differential responses of WM and GCH 2 genotypes to As(V) stress.

Ultrasensitive electroluminescence biosensor for a breast cancer marker microRNA based on target cyclic regeneration and multi-labeled magnetized nanoparticles

An electrochemiluminescent (ECL) biosensor is described for the determination of the breast cancer biomarker microRNA. The method is based on the amplification via target cyclic regeneration through a system of hairpin DNA probes, primers, and Klenow fragment of DNA polymerases combined with CdTe quantum dots (QDs) and gold nanoparticles. The assay is performed by exploiting the luminescence properties of CdTe-QDs and K₂S₂O₈ as a co-reactive agent to increase the ECL signal. It was successfully applied to ECL-based detection of a 20-mer microRNA. The sensor has a linear response in the 0.1 fM to 0.2 pM microRNA concentration range and a detection limit as low as 33 aM. The assay has been applied to the determination of microRNA spiked in serum samples, and recoveries ranged from 94.4 to 100.5%. Graphical abstract A novel electroluminescence biosensor based on the amplification of target cyclic regeneration is described. It is achieved by using a system of hairpin DNA probes, primers, and Klenow fragment of DNA polymerases combined with CdTe QDs and Au NPs, and was successfully applied to microRNA detection.

Tissue-specific transcriptional regulation of seven heavy metal stress-responsive miRNAs and their putative targets in nickel indicator castor bean (R. communis L.) plants

R.communis L. has high capability to accumulate nickel which is a trace nutrient for higher plants and also an environmental contaminant causes toxicity related symptoms at higher concentrations. MicroRNAs (miRNAs) are known to be important modulators of responses against heavy metal stress for detoxification of the metal. In this study, we experimentally measured and validated the transcript levels of the seven heavy metal stress response-related miRNAs and the expression levels of target genes in both leaf and root tissues of R. communis L. subjected to three different concentrations of nickel stress via qRT-PCR quantification. The results demonstrated differential regulations of heavy metal stress-responsive miRNAs and their putative targets in both tissues in same stress treatments. This dynamic regulation suggest that regulatory processes differ between the tissues under nickel stress. Our data suggest that, miR838 was the most responsive to the Ni2+ stress. miR398 target gene Cu-Zn/SOD was found to be up-regulated in both root and leaf tissues. The relations between TCP and expression levels of miR159 and miR319 were also found statistically significant exclusive to leaf tissues. In leaf tissue, changes in miR395 level and its putative target genes, sulphate transporter and sulphate adenyltransferase gene were found in relation whereas, only expression level of sulphate transporter represented a statistically significant relation in root tissue. The sharp decrease in transcript levels of 2r3 myb gene at lower nickel dose suggest to investigate the role of r2r3 myb and the all MYB family members in primary and secondary metabolisms against nickel stress.

Identifying optimal reference genes for the normalization of microRNA expression in cucumber under viral stress

Cucumber green mottle mosaic virus (CGMMV) is an economically important pathogen and causes significant reduction of both yield and quality of cucumber (Cucumis sativus). Currently, there were no satisfied strategies for controlling the disease. A better understanding of microRNA (miRNA) expression related to the regulation of plant-virus interactions and virus resistance would be of great assistance when developing control strategies for CGMMV. However, accurate expression analysis is highly dependent on robust and reliable reference gene used as an internal control for normalization of miRNA expression. Most commonly used reference genes involved in CGMMV-infected cucumber are not universally expressed depending on tissue types and stages of plant development. It is therefore crucial to identify suitable reference genes in investigating the role of miRNA expression. In this study, seven reference genes, including Actin, Tubulin, EF-1α, 18S rRNA, Ubiquitin, GAPDH and Cyclophilin, were evaluated for the most accurate results in analyses using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Gene expression was assayed on cucumber leaves, stems and roots that were collected at different days post inoculation with CGMMV. The expression data were analyzed using algorithms including delta-Ct, geNorm, NormFinder, and BestKeeper as well as the comparative tool RefFinder. The reference genes were subsequently validated using miR159. The results showed that EF-1α and GAPDH were the most reliable reference genes for normalizing miRNA expression in leaf, root and stem samples, while Ubiquitin and EF-1α were the most suitable combination overall.

Transcriptomic analyses reveal complex and interconnected sucrose signaling cascades in developing seeds of castor bean

Seeds are highly specific organs that strongly sink sucrose resources from leaf and stem tissues to trigger seed metabolism and development. In particular, for heterotrophic non-green seeds, the potential molecular mechanism underlying sucrose-driven seed development remains unanswered. Castor bean (Ricinus communis L.), a typical non-green seed, has been considered as a model plant for seed biology study in dicotyledonous plants due to its heterotrophic seeds with persistent endosperms. In the present study, the fast-developing castor bean seeds were treated with exogenous sucrose and mannitol for four hours. The global transcriptomic data were obtained by high-throughput RNA-seq technique, resulting in 468 differentially expressed genes (DGEs). Further analyses revealed that sucrose functioned as both metabolic substrates and signal molecules. Specifically, 73 DGEs involved in carbohydrate and nitrogen metabolism, 42 differentially expressed transcription factors, and 35 DGEs involved in diverse signaling pathways such as auxin, brassinosteroid, ethelyene, cytokinin, gibberellin, and calcium signals, were identified, suggesting that the sucrose signaling pathway might have complex and multi-connected cross-talks with other signals to regulate castor bean seed development. Taken together, this study provides novel data to improve understanding of the potential molecular mechanisms of sucrose in regulating non-green seed development and storage reservoir accumulation during seed development.

Automated identification of reference genes based on RNA-seq data

Background

Gene expression analyses demand appropriate reference genes (RGs) for normalization, in order to obtain reliable assessments. Ideally, RG expression levels should remain constant in all cells, tissues or experimental conditions under study. Housekeeping genes traditionally fulfilled this requirement, but they have been reported to be less invariant than expected; therefore, RGs should be tested and validated for every particular situation. Microarray data have been used to propose new RGs, but only a limited set of model species and conditions are available; on the contrary, RNA-seq experiments are more and more frequent and constitute a new source of candidate RGs.

Results

An automated workflow based on mapped NGS reads has been constructed to obtain highly and invariantly expressed RGs based on a normalized expression in reads per mapped million and the coefficient of variation. This workflow has been tested with Roche/454 reads from reproductive tissues of olive tree (Olea europaea L.), as well as with Illumina paired-end reads from two different accessions of Arabidopsis thaliana and three different human cancers (prostate, small-cell cancer lung and lung adenocarcinoma). Candidate RGs have been proposed for each species and many of them have been previously reported as RGs in literature. Experimental validation of significant RGs in olive tree is provided to support the algorithm.

Conclusion

Regardless sequencing technology, number of replicates, and library sizes, when RNA-seq experiments are designed and performed, the same datasets can be analyzed with our workflow to extract suitable RGs for subsequent PCR validation. Moreover, different subset of experimental conditions can provide different suitable RGs.

Electronic supplementary material

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