Selection and validation of suitable reference genes for miRNA expression normalization by quantitative RT-PCR in citrus somatic embryogenic and adult tissues

The miR159a-DUO1 module regulates pollen development by modulating auxin biosynthesis and starch metabolism in citrus

Achieving seedlessness in citrus varieties is one of the important objectives of citrus breeding. Male sterility associated with abnormal pollen development is an important factor in seedlessness. However, our understanding of the regulatory mechanism underlying the seedlessness phenotype in citrus is still limited. Here, we determined that the miR159a-DUO1 module played an important role in regulating pollen development in citrus, which further indirectly modulated seed development and fruit size. Both the overexpression of csi-miR159a and the knocking out of DUO1 in Hong Kong kumquat (Fortunella hindsii) resulted in small and seedless fruit phenotypes. Moreover, pollen was severely aborted in both transgenic lines, with arrested pollen mitotic I and abnormal pollen starch metabolism. Through additional cross-pollination experiments, DUO1 was proven to be the key target gene for miR159a to regulate male sterility in citrus. Based on DNA affinity purification sequencing (DAP-seq), RNA-seq, and verified interaction assays, YUC2/YUC6, SS4 and STP8 were identified as downstream target genes of DUO1, those were all positively regulated by DUO1. In transgenic F. hindsii lines, the miR159a-DUO1 module down-regulated the expression of YUC2/YUC6, which decreased indoleacetic acid (IAA) levels and modulated auxin signaling to repress pollen mitotic I. The miR159a-DUO1 module reduced the expression of the starch synthesis gene SS4 and sugar transport gene STP8 to disrupt starch metabolism in pollen. Overall, this work reveals a new mechanism by which the miR159a-DUO1 module regulates pollen development and elucidates the molecular regulatory network underlying male sterility in citrus.

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L

BACKGROUND

Theobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao.

RESULTS

We describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3-5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues.

CONCLUSION

The Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

Screening and validation of optimal miRNA reference genes in different developing stages and tissues of Lilium henryi Baker

Dynamic miRNA detection using the qRT-PCR technique requires appropriate reference genes to ensure data reliability. Previous studies have screened internal reference genes in plants during embryonic development and various stress treatment, involving relatively few tissues and organs. There is no relevant miRNA study in Lilium henryi Baker and limited research on the optimal miRNA reference genes in lilies, such as 5S, 18S, U6 and Actin. Twelve genes were selected as candidate reference genes whose expression stability was analyzed in petals at different developmental stages and other tissues using various algorithms, such as geNorm, NormFinder, BestKeeper, and Delta CT. The results revealed that the optimal combination of reference genes for Lilium henryi Baker petals at different developmental stages was osa-miR166m and osa-miR166a-3p, while that for different tissues of Lilium henryi Baker was osa-miR166g-3p and osa-miR166a-3p.Four important genes related to growth and development regulation, namely, osa-miR156a, osa-miR395b, osa-miR396a-3p, and osa-miR396a-5p, were selected for validation. The findings of the present study could contribute to future investigations onmiRNA expression and the related functions in Lilium henryi Baker while providing important references for the normalization of the miRNA expression in other varieties of lily.

Identifying the Soybean microRNAs Related to Phytophthora sojae Based on RNA Sequencing and Bioinformatics Analysis

Phytophthora root rot in soybeans is caused by a pathogen called Phytophthora sojae (P. sojae), which results in a significant decrease in soybean production within affected regions. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that play a key post-transcriptional regulatory role in eukaryotes. In this paper, the miRNAs that respond to P. sojae were analyzed from the gene level to complement the study of molecular resistance mechanisms in soybean. The study utilized high-throughput sequencing of soybean data to predict miRNAs that respond to P. sojae, analyze their specific functions, and verify regulatory relationships using qRT-PCR. The results showed that the miRNAs in soybean respond to P. sojae infection. MiRNAs can be transcribed independently, suggesting the presence of transcription factor binding sites in the promoter regions. Additionally, we performed an evolutionary analysis on conserved miRNAs that respond to P. sojae. Finally, we investigated the regulatory relationships among miRNAs, genes, and transcription factors, and identified five regulatory patterns. These findings lay the groundwork for future studies on the evolution of miRNAs responsive to P. sojae.

Differential Expression of miRNAs Involved in Response to Candidatus Liberibacter asiaticus Infection in Mexican Lime at Early and Late Stages of Huanglongbing Disease

Huanglongbing (HLB) is one of the most destructive diseases threatening citriculture worldwide. This disease has been associated with α-proteobacteria species, namely Candidatus Liberibacter. Due to the unculturable nature of the causal agent, it has been difficult to mitigate the disease, and nowadays a cure is not available. MicroRNAs (miRNAs) are key regulators of gene expression, playing an essential role in abiotic and biotic stress in plants including antibacterial responses. However, knowledge derived from non-model systems including Candidatus Liberibacter asiaticus (CLas)-citrus pathosystem remains largely unknown. In this study, small RNA profiles from Mexican lime (Citrus aurantifolia) plants infected with CLas at asymptomatic and symptomatic stages were generated by sRNA-Seq, and miRNAs were obtained with ShortStack software. A total of 46 miRNAs, including 29 known miRNAs and 17 novel miRNAs, were identified in Mexican lime. Among them, six miRNAs were deregulated in the asymptomatic stage, highlighting the up regulation of two new miRNAs. Meanwhile, eight miRNAs were differentially expressed in the symptomatic stage of the disease. The target genes of miRNAs were related to protein modification, transcription factors, and enzyme-coding genes. Our results provide new insights into miRNA-mediated regulation in C. aurantifolia in response to CLas infection. This information will be useful to understand molecular mechanisms behind the defense and pathogenesis of HLB.

miR156 regulates somatic embryogenesis by modulating starch accumulation in citrus

Somatic embryogenesis (SE) is a major regeneration approach for in vitro cultured tissues of plants, including citrus. However, SE capability is difficult to maintain, and recalcitrance to SE has become a major obstacle to plant biotechnology. We previously reported that miR156-SPL modules regulate SE in citrus callus. However, the downstream regulatory pathway of the miR156-SPL module in SE remains unclear. In this study, we found that transcription factors CsAGL15 and CsFUS3 bind to the CsMIR156A promoter and activate its expression. Suppression of csi-miR156a function leads to up-regulation of four target genes, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (CsSPL) genes, and reduction of SE efficiency. In the short tandem target mimic (STTM)-miR156a overexpression callus (MIM156), the number of amyloplasts and starch content were significantly reduced, and genes involved in starch synthesis and transport were down-regulated. csi-miR172d was down-regulated, whereas the target genes, CsTOE1.1 and CsTOE1.2, which inhibit the expression of starch biosynthesis genes, were up-regulated. In our working model, CsAGL15 and CsFUS3 activate csi-miR156a, which represses CsSPLs and further regulates csi-miR172d and CsTOEs, thus altering starch accumulation in callus cells and regulating SE in citrus. This study elucidates the pathway of miR156-SPLs and miR172-TOEs-mediated regulation of SE, and provides new insights into enhancing SE capability in citrus.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

miR-434 and miR-242 have a potential role in heat stress response in rainbow trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are being extensively studied as they function as key metabolic regulators which play a role in the heat stress response. However, the role of miRNAs in heat stress remains uncertain and many new miRNAs have not yet been discovered. In a previous study, we performed high-throughput sequencing of differentially expressed miRNAs identified on exposing rainbow trout (Oncorhynchus mykiss) to heat stress (18 vs. 24°C), which led to the identification of two novel miRNAs, temporarily named novel miR-434 and -242. The differential expression level of these miRNAs was extremely significant (P < 0.01); we analysed target gene mRNA transcripts by bioinformatics software (miRanda). We found novel miR-434 and -242 were predicted to regulate the transcripts of heat shock 70-kDa protein 4-like (HSPA4L) and calreticulin (CRT), respectively, by bioinformatics software. Here our core objective was to validate if HSPA4L and CRT are indeed the target genes of novel miR-434 and -242, respectively, and for this purpose we used the dual-luciferase reporter assay system. Target gene sequences were synthesized and cloned into a dual-luciferase vector. To better understand the function of the target genes, we combined the previous sequencing results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. We found that novel miR-434 regulated HSPA4L expression by binding to a putative binding site in the 3'-UTR of HSPA4L, and luciferase activity inhibition was observed. In contrast, novel miR-242 was not involved in regulating CRT expression. To conclude, we believe our results should serve as a foundation for future studies aiming to comprehensively understand the mechanisms used by rainbow trout to cope with heat stress.

Selection and Optimization of Reference Genes for MicroRNA Expression Normalization by qRT-PCR in Chinese Cedar (Cryptomeria fortunei) under Multiple Stresses

MicroRNA (miRNA) expression analysis is very important for investigating its functions. To date, no research on reference genes (RGs) for miRNAs in gymnosperms, including Cryptomeria fortunei, has been reported. Here, ten miRNAs (i.e., pab-miR159a, cln-miR162, cas-miR166d, pab-miR395b, ppt-miR894, cln-miR6725, novel1, novel6, novel14 and novel16) and three common RGs (U6, 5S and 18S) were selected as candidate RGs. qRT-PCR was used to analyse their expressions in C. fortunei under various experimental conditions, including multiple stresses (cold, heat, drought, salt, abscisic acid and gibberellin) and in various tissues (roots, stems, tender needles, cones and seeds). Four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) were employed to assess the stability of candidate RG expression; the geometric mean and RefFinder program were used to comprehensively evaluate RG stability. According to the results, novel16, cln-miR6725, novel1 and U6 were the most stable RGs for studying C. fortunei miRNA expression. In addition, the expression of three target miRNAs (aly-miR164c-5p, aly-miR168a-5p and smo-miR396) was examined to verify that the selected RGs are suitable for miRNA expression normalisation. This study may aid further investigations of miRNA expression/function in the response of C. fortunei to abiotic stress and provides an important basis for the standardisation of miRNA expression in other gymnosperm species.

MicroRNAs Involved in Regulatory Cytoplasmic Male Sterility by Analysis RNA-seq and Small RNA-seq in Soybean

Cytoplasmic male sterility (CMS) is an important plant characteristic for exploiting heterosis to enhance crop traits during breeding. However, the CMS regulatory network remains unclear in plants, even though researchers have attempted to isolate genes associated with CMS. In this study, we performed high-throughput sequencing and degradome analyses to identify microRNAs (miRNAs) and their targets in a soybean CMS line (JLCMS9A) and its maintainer line (JLCMS9B). Additionally, the differentially expressed genes during reproductive development were identified using RNA-seq data. A total of 280 miRNAs matched soybean miRNA sequences in miRBase, including mature miRNAs and pre-miRNAs. Of the 280 miRNAs, 30, 23, and 21 belonged to the miR166, miR156, and miR171 families, respectively. Moreover, 410 novel low-abundant miRNAs were identified in the JLCMS9A and JLCMS9B flower buds. Furthermore, 303 and 462 target genes unique to JLCMS9A and JLCMS9B, respectively, as well as 782 common targets were predicted based on the degradome analysis. Target genes differentially expressed between the CMS line and the maintainer line were revealed by an RNA-seq analysis. Moreover, all target genes were annotated with diverse functions related to biological processes, cellular components, and molecular functions, including transcriptional regulation, the nucleus, meristem maintenance, meristem initiation, cell differentiation, auxin-activated signaling, plant ovule development, and anther development. Finally, a network was built based on the interactions. Analyses of the miRNA, degradome, and transcriptome datasets generated in this study provided a comprehensive overview of the reproductive development of a CMS soybean line. The data presented herein represent useful information for soybean hybrid breeding. Furthermore, the study results indicate that miRNAs might contribute to the soybean CMS regulatory network by modulating the expression of CMS-related genes. These findings lay the foundation for future studies on the molecular mechanisms underlying soybean CMS.

Identification and validation of miRNA reference genes in poplar under pathogen stress

Quantitative real time polymerase chain reaction (qRT-PCR) is a common method to analyze gene expression. Due to differences in RNA quantity, quality, and reverse transcription efficiency between qRT-PCR samples, reference genes are used as internal standards to normalize gene expression. However, few universal genes, especially miRNAs, have been identified as reference so far. Therefore, it is essential to identify reference genes that can be used across various experimental conditions, stress treatments, or tissues. In this study, 14 microRNAs (miRNAs) and 5.8S rRNA were assessed for expression stability in poplar trees infected with canker pathogen. Using geNorm, NormFinder and Bestkeeper reference gene analysis programs, we found that miR156g and miR156a exhibited stable expression throughout the infection process. miR156g, miR156a and 5.8S rRNA were then tested as internal standards to measure the expression of miR1447 and miR171c, and the results were compared to small RNA sequencing (RNA-seq) data. We found that when miR156a and 5.8S rRNA were used as the reference gene, the expression of miR1447 and miR171c were consistent with the small RNA-seq expression profiles. Therefore, miR156a was the most stable miRNAs examined in this study, and could be used as a reference gene in poplar under canker pathogen stress, which should enable comprehensive comparisons of miRNAs expression and avoid the bias caused by different length between detected miRNAs and traditional reference genes. The present study has expanded the miRNA reference genes available for gene expression studies in trees under biotic stress.

Next-Generation Sequencing Identification and Characterization of MicroRNAs in Dwarfed Citrus Trees Infected With Citrus Dwarfing Viroid in High-Density Plantings

Citrus dwarfing viroid (CDVd) induces stunting on sweet orange trees [Citrus sinensis (L.) Osbeck], propagated on trifoliate orange rootstock [Citrus trifoliata (L.), syn. Poncirus trifoliata (L.) Raf.]. MicroRNAs (miRNAs) are a class of non-coding small RNAs (sRNAs) that play important roles in the regulation of tree gene expression. To identify miRNAs in dwarfed citrus trees, grown in high-density plantings, and their response to CDVd infection, sRNA next-generation sequencing was performed on CDVd-infected and non-infected controls. A total of 1,290 and 628 miRNAs were identified in stem and root tissues, respectively, and among those, 60 were conserved in each of these two tissue types. Three conserved miRNAs (csi-miR479, csi-miR171b, and csi-miR156) were significantly downregulated (adjusted p-value < 0.05) in the stems of CDVd-infected trees compared to the non-infected controls. The three stem downregulated miRNAs are known to be involved in various physiological and developmental processes some of which may be related to the characteristic dwarfed phenotype displayed by CDVd-infected C. sinensis on C. trifoliata rootstock field trees. Only one miRNA (csi-miR535) was significantly downregulated in CDVd-infected roots and it was predicted to target genes controlling a wide range of cellular functions. Reverse transcription quantitative polymerase chain reaction analysis performed on selected miRNA targets validated the negative correlation between the expression levels of these targets and their corresponding miRNAs in CDVd-infected trees. Our results indicate that CDVd-responsive plant miRNAs play a role in regulating important citrus growth and developmental processes that may participate in the cellular changes leading to the observed citrus dwarf phenotype.

Selection of reference genes for gene expression analysis in Liriodendron hybrids' somatic embryogenesis and germinative tissues

The differential expression of genes is crucial for plant somatic embryogenesis (SE), and the accurate quantification of gene expression levels relies on choosing appropriate reference genes. To select the most suitable reference genes for SE studies, 10 commonly used reference genes were examined in synchronized somatic embryogenic and subsequent germinative cultures of Liriodendron hybrids by using quantitative real-time reverse transcription PCR. Four popular normalization algorithms: geNorm, NormFinder, Bestkeeper and Delta-Ct were used to select and validate the suitable reference genes. The results showed that elongation factor 1-gamma, histone H1 linker protein, glyceraldehyde-3-phosphate dehydrogenase and α-tubulin were suitable for SE tissues, while elongation factor 1-gamma and actin were best for the germinative organ tissues. Our work will benefit future studies of gene expression and functional analyses of SE in Liriodendron hybrids. It is also serves as a guide of reference gene selection in early embryonic gene expression analyses for other woody plant species.

Identification and function prediction of iron-deficiency-responsive microRNAs in citrus leaves

Iron is a critical micronutrient for growth and development of plants and its deficiency limiting the crop productivity. MicroRNAs (miRNAs) play vital roles in adaptation of plants to various nutrient deficiencies. However, the role of miRNAs and their target genes related to Fe-deficiency is limited. In this study, we identified Fe-deficiency-responsive miRNAs from citrus. In Fe-deficiency conditions, about 50 and 31 miRNAs were up-regulated and down-regulated, respectively. The differently expressed miRNAs might play critical roles in contributing the Fe-deficiency tolerance in citrus plants. The miRNAs-mediated Fe-deficiency tolerance in citrus plants might related to the enhanced stress tolerance by decreased expression of miR172; regulation of S homeostasis by decreased expression of miR395; inhibition of plant growth by increased expression of miR319 and miR477; regulation of Cu homeostasis as well as activation of Cu/Zn superoxide dismutase activity due to decreased expression of miR398 and miR408 and regulation of lignin accumulation by decreased expression of miR397 and miR408. The identified miRNAs in present study laid a foundation to understand the Fe-deficiency adaptive mechanisms in citrus plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02669-z.

Optimization of reference genes for qRT-PCR analysis of microRNA expression under abiotic stress conditions in sweetpotato

Sweetpotato (Ipomoea batatas. L) is an important food crop, harvested for its nutrient-rich tuberous roots. Drought and salt stresses are two major factors limiting the sweetpotato production. Since microRNAs (miRNAs) are well known to play crucial roles in regulation of plant stress responses, quantitative profiling of miRNA expression under stress conditions will facilitate identification and genetic manipulation of novel miRNAs to improve stress tolerance. Real-time quantitative reverse transcription PCR (qRT-PCR) is a commonly used tool for this purpose, but not without challenges. Although stem-loop and poly(A)-tail modified qRT-PCR methods were developed for characterizing miRNA expression, accurate profiling of miRNAs is still difficult in many plant species because of a lack of reliable reference genes for normalizing miRNA transcripts. To identify reference genes that are suitable for normalizing miRNA expression in sweetpotato, the expression stability of eight candidate miRNAs and two commonly used reference genes were tested in 96 samples involving four tissues and two cultivars under drought and salt stress treatments. Data analysis using the geNorm, NormFinder and Bestkeeper algorithms demonstrated that miRn60, miR482, and their combination were reliable references. We further validated the reference genes by expression analysis of the well-characterized miR319 and miR156 that regulate drought and salt stress responses, respectively. The reference genes identified in this study will facilitate future miRNA analysis under abiotic stress conditions in sweetpotato.

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.

Evaluation of reference genes for quantitative real-time PCR analysis of messenger RNAs and microRNAs in rainbow trout Oncorhynchus mykiss under heat stress

We assessed the expression stability of several messenger (m)RNAs and micro (mi)RNAs from liver and head kidney of rainbow trout Oncorhynchus mykiss using high-throughput RNA sequencing (RNA-seq) and miRNA-seq data. Additionally, four commonly used reference genes and one small non-coding RNA (u6) were also selected to identify ideal reference mRNAs and miRNAs for quantitative real-time (qrt)-PCR analysis of heat stress responses. GeNorm, NormFinder, BestKeeper and comparative ΔC_t were employed for analysis of qrt-PCR data to systematically assess the expression stability of candidate mRNAs and miRNAs and stability was ranked using geometric means. β-actin and ef1-α were the most stably expressed reference mRNAs in liver and head kidney, respectively and ssa-mir-26a-5p and ssa-mir-462b-5p were the most stably expressed miRNAs in these tissues. This is the first identification of appropriate reference mRNAs and miRNAs for qrt-PCR analysis of O. mykiss under heat stress.

Selection of suitable candidate genes for miRNA expression normalization in Yellow River Carp (Cyprinus carpio. var)

Yellow River carp is widely cultivated in the world due to its economic value in aquaculture, and the faster growth of females compared to males. It is believed that microRNAs (miRNA) are involved in gonadal differentiation and development. qPCR is the most preferred method for miRNA functional analysis. Reliable reference genes for normalization in qRT-PCR are the key to ensuring the accuracy of this method. The aim of present research was to evaluate as well as identify the efficacy of reference genes for miRNA expression using qRT-PCR in Yellow River carp. Nine ncRNAs (miR-101, miR-23a, let7a, miR-26a, miR-146a, miR-451, U6, 5S, and 18S) were chosen and tested in four sample sets: (1) different tissues in adult carp, (2) different tissues in juvenile carp, (3) different early developmental stages of carp, and (4) different developmental stages of carp gonads. The stability and suitability values were calculated using NormFinder, geNorm, and BestKeeper software. The results showed that 5S was a suitable reference gene in different tissues of adult and juvenile carp. The genes 5S, 18S, and U6 were the most stable reference genes in the early developmental stages of carp. Let-7a and miR-23a were considered as the suitable reference genes in the development of gonads. All these reference genes were subsequently validated using miR-430. The results showed that genes 5S and 18S were the most suitable reference genes to normalize miRNA expression under normal growth conditions in early different developmental stages. The genes Let-7a, and miR-23a were the most suitable in different developmental stages. The present study is the first comprehensive study of the stability of miRNA reference genes in Yellow River carp, providing valuable as well as basic data for investigating more accurate miRNA expression during gonadal differentiation and development of carp.

Genome-wide comprehensive analysis of transcriptomes and small RNAs offers insights into the molecular mechanism of alkaline stress tolerance in a citrus rootstock

Alkaline stress has serious-negative effects on citrus production. Ziyang xiangcheng (Citrus junos Sieb. ex Tanaka) (Cj) is a rootstock that is tolerant to alkaline stress and iron deficiency. Trifoliate orange (Poncirus trifoliata (L.) Raf.) (Pt), the most widely used rootstock in China, is sensitive to alkaline stress. To investigate the molecular mechanism underlying the tolerance of Cj to alkaline stress, next-generation sequencing was employed to profile the root transcriptomes and small RNAs of Cj and Pt seedlings that were cultured in nutrient solutions along a three pH gradient. This two-level regulation data set provides a system-level view of molecular events with a precise resolution. The data suggest that the auxin pathway may play a central role in the inhibitory effect of alkaline stress on root growth and that the regulation of auxin homeostasis under alkaline stress is important for the adaptation of citrus to alkaline stress. Moreover, the jasmonate (JA) pathway exhibits the opposite response to alkaline stress in Cj and Pt and may contribute to the differences in the alkaline stress tolerance and iron acquisition between Cj and Pt. The dataset provides a wealth of genomic resources and new clues to further study the mechanisms underlying alkaline stress resistance in Cj.

Identification of appropriate reference genes for RT-qPCR analysis in Juglans regia L

Reverse transcription quantitative real-time PCR (RT-qPCR) is a popular adopted technique to detect gene expression, and the selection of appropriate reference genes is crucial for data normalization. In the present study, seven candidate reference genes were screened to evaluate their expression stability in various flower buds, leaf buds, tissues and cultivars of the English walnut (Juglans regia L.) based on four algorithms (geNorm, Normfinder, Bestkeeper and RefFinder). The results demonstrated that TUA, EF1 and TUB were appropriate reference genes for flower buds at different stages of female flower buds differentiation; TUB and 18S rRNA were best for leaf buds at different stages of female flower buds differentiation; TUB and TUA were suitable for different cultivars; and ACT2, 18S rRNA and GAPDH were useful for different tissues. Moreover, the expression of ACT was not stable among different flower buds, leaf buds and cultivars. The stability of reference genes were confirmed through the analysis of the expression of SPL18 gene. These results will contribute to a reliable normalization of gene expression in J. regia.

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

Regulation mechanism of microRNA in plant response to abiotic stress and breeding

microRNAs (miRNAs) in plants are a class of small RNAs consisting of approximately 21-24 nucleotides. The mature miRNA binds to the target mRNA through the formation of a miRNA-induced silencing complex (MIRISC), and cleaves or inhibits translation, thereby achieving negative regulation of the target gene. Based on miRNA plays an important role in regulating plant gene expression, studies on the prediction, identification, function and evolution of plant miRNAs have been carried out. In addition, many researches prove that miRNAs are also involved in many kinds of abiotic and biotic stress, under abiotic stress, plants can express some miRNA, and act on stress-related target genes, which can make plants adapt to stress in physiological response. In this review, the synthetic pathway and mechanism of plant miRNA are briefly described, and we discuss the biological functions and regulatory mechanisms of miRNAs responding to abiotic stresses including low temperature, salt, drought stress and breeding to lay the foundation for further exploring the mechanism of action of miRNAs in stress resistance of plant. And analyze its utilization prospects in plant stress resistance research.

Selection and validation of suitable reference genes for qRT-PCR analysis in pear leaf tissues under distinct training systems

Training systems generally alter tree architecture, which modulates light microclimate within the canopy, for the purpose of improving photosynthetic efficiency and fruit quality. Gene expression quantification is one of the most important methods for exploring the molecular mechanisms underlying the influence of training systems on pear photosynthesis, and suitable reference genes for gene expression normalization are a prerequisite for this method. In this study, the expression stability of nine common and four novel candidate genes were evaluated in 14 different pear leaf samples in two training systems, including those at four developmental stages (training_period) and from different parts of the trees (training_space), using two distinct algorithms, geNorm and NormFinder. Our results revealed that SKD1 (Suppressor of K⁺Transport Growth Defect1)/ YLS8 (Yellow Leaf Specific 8) and ARM (Armadillo) were the most stable single reference genes for the ‘training_period’ and ‘training_space’ subsets, respectively, although these single genes were not as stable as the optimal pairs of reference genes, SKD1+YLS8 and ARM+YLS8, respectively. Furthermore, the expression levels of the PpsAPX (Ascorbate peroxidase) gene showed that the arbitrary use of reference genes without previous testing could lead to misinterpretation of data. This work constitutes the first systematic analysis regarding the selection of superior reference genes in training system studies, facilitating the elucidation of gene function in pear and providing valuable information for similar studies in other higher plants.

Selection and Validation of Reference Genes for mRNA Expression by Quantitative Real-Time PCR Analysis in Neolamarckia cadamba

Neolamarckia cadamba is an economically-important fast-growing tree species in South China and Southeast Asia. As a prerequisite first step for future gene expression studies, we have identified and characterized a series of stable reference genes that can be used as controls for quantitative real time PCR (qRT-PCR) expression analysis in this study. The expression stability of 15 candidate reference genes in various tissues and mature leaves under different conditions was evaluated using four different algorithms, i.e., geNorm, NormFinder, BestKeeper and RefFinder. Our results showed that SAMDC was the most stable of the selected reference genes across the set of all samples, mature leaves at different photosynthetic cycles and under drought stress, whereas RPL10A had the most stable expression in various tissues. PGK and RPS25 were considered the most suitable reference for mature leaves at different developmental stages and under cold treatment, respectively. Additionally, the gene expression profiles of sucrose transporter 4 (NcSUT4), and 9-cis-epoxycarotenoid dioxygenase 3 (NcNCED3) were used to confirm the validity of candidate reference genes. Collectively, our study is the first report to validate the optimal reference genes for normalization under various conditions in N. cadamba and will benefit the future discovery of gene function in this species.

miR156-SPL modules regulate induction of somatic embryogenesis in citrus callus

miR156 is a highly conserved plant miRNA and has been extensively studied because of its versatile roles in plant development. Here, we report a novel role of miR156 in regulating somatic embryogenesis (SE) in citrus, one of the most widely cultivated fruit crops in the world. SE is an important means of in vitro regeneration, but over the course of long-term sub-culturing there is always a decline in the SE potential of the preserved citrus embryogenic callus, and this represents a key obstacle for citrus biotechnology. In this study, the SE competence of citrus callus of wild kumquat (Fortunella hindsii) was significantly enhanced by either overexpression of csi-miR156a or by individual knock-down of the two target genes, CsSPL3 and CsSPL14, indicating that the effect of miR156-SPL modules was established during the initial phases of SE induction. Biological processes that might promote SE in response to miR156 overexpression were explored using RNA-seq, and mainly included hormone signaling pathways, stress responses, DNA methylation, and the cell cycle. CsAKIN10 was identified as interacting protein of CsSPL14. Our results provide insights into the regulatory pathway through which miR156-SPL modules enhance the SE potential of citrus callus, and provide a theoretical basis for improvement of plant SE competence.

Identification and Expression of miRNAs Related to Female Flower Induction in Walnut (Juglans regia L.)

Flower induction is an essential stage in walnut (Juglans regia L.) trees, directly affecting yield, yield stability, fruit quality and commodity value. The objective of this study was to identify miRNAs related to female flower induction via high-throughput sequencing and bioinformatics analysis. A total of 123 miRNAs were identified including 51 known miRNAs and 72 novel miRNAs. Differential expression was observed in 19 of the known miRNAs and 34 of the novel miRNAs. Twelve miRNAs were confirmed by RT-qPCR. A total of 1339 target genes were predicted for the differentially expressed miRNAs. The functions of 616 of those target genes had been previously annotated. The target genes of the differentially expressed miRNAs included: (i) floral homeotic protein APETALA 2 (AP2) and ethylene-responsive transcription factor RAP2-7 which were targeted by jre-miRn69; (ii) squamosa promoter-binding protein 1 (SPB1) and various SPLs (squamosa promoter-binding-like protein) which were targeted by jre-miR157a-5p; (iii) various hormone response factors which were targeted by jre-miR160a-5p (ARF18) and jre-miR167a-5p (ARF8) and (iv) transcription factor SCL6 which was targeted by jre-miR171b-3p, jre-miRn46 and jre-miRn49. The KEGG pathway analysis of the target genes indicated that the differentially expressed miRNAs were mainly enriched to ubiquitin mediated proteolysis, RNA degradation and various carbohydrate metabolism pathways. Many miRNAs were detected in J. regia during female flower induction. Some miRNAs (jre-miR157a-5p, jre-miR160a-5p, jre-miR167a-5p, miR171b-3p jre-miRn69 and jre-miRn49) were involved in female flower induction. The results of this experiment will contribute valuable information for further research about the function of miRNAs in flower induction of J. regia and other fruit trees.

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.

Selection of reference genes for miRNA qRT-PCR under abiotic stress in grapevine

Grapevine is among the fruit crops with high economic value, and because of the economic losses caused by abiotic stresses, the stress resistance of Vitis vinifera has become an increasingly important research area. Among the mechanisms responding to environmental stresses, the role of miRNA has received much attention recently. qRT-PCR is a powerful method for miRNA quantitation, but the accuracy of the method strongly depends on the appropriate reference genes. To determine the most suitable reference genes for grapevine miRNA qRT-PCR, 15 genes were chosen as candidate reference genes. After eliminating 6 candidate reference genes with unsatisfactory amplification efficiency, the expression stability of the remaining candidate reference genes under salinity, cold and drought was analysed using four algorithms, geNorm, NormFinder, deltaCt and Bestkeeper. The results indicated that U6 snRNA was the most suitable reference gene under salinity and cold stresses; whereas miR168 was the best for drought stress. The best reference gene sets for salinity, cold and drought stresses were miR160e + miR164a, miR160e + miR168 and ACT + UBQ + GAPDH, respectively. The selected reference genes or gene sets were verified using miR319 or miR408 as the target gene.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Identification of novel miRNAs and miRNA expression profiling in embryogenic tissues of Picea balfouriana treated by 6-benzylaminopurine

Here, we compared miRNA expression profiles in embryonic cell cultures of the conifer Picea balfouriana following application of the synthetic cytokinin 6-benzylaminopurine (6-BAP). We used next-generation sequencing to analyze three libraries of small RNAs from the treated embryogenic cell cultures and generated 24,000,000 raw reads from each of the libraries. Over 70 differentially regulated micro RNA (miRNA) families (≥2 fold change in expression) were identified between pairs of treatments. A quantitative analysis showed that miR3633 and miR1026 were upregulated in tissues with the highest embryogenic ability. These two miRNAs were predicted to target genes encoding receptor-like protein kinase and GAMYB transcription factors, respectively. In one library, miR1160, miR5638, miR1315, and miR5225 were downregulated. These four miRNAs were predicted to target genes encoding APETALA2, calmodulin-binding protein, and calcium-dependent protein kinase transcription factors. The expression patterns of the miRNAs and their targets were negatively correlated. Approximately 181 potentially novel P. balfouriana miRNAs were predicted from the three libraries, and seven were validated during the quantitative analysis. This study is the first report of differential miRNA regulation in tissues treated with 6-BAP during somatic embryogenesis. The differentially expressed miRNAs will be of value for investigating the mechanisms of embryogenic processes that are responsive to 6-BAP in P. balfouriana.

microRNAs as reference genes for quantitative PCR in cotton

Cotton (Gossypium hirsutum) is the most important non-food plant in the world. Studies concerning the fiber quality and plant fitness of cotton at molecular level depend on high sensitive and reproducible gene-expression assays. However, only a few reports have described genes suitable for normalizing gene expression data. In this study, we report for the first time that microRNAs (miRNAs) are reliable reference genes (RGs) for cotton gene expression normalization in quantitative real-time reverse transcription (RT)-PCR. The stability of cotton miRNAs was assayed in root, stem, leaf and flower samples from three different cultivars [FiberMax (FM966), Delta Opal (DO) and Cedro] and under conditions of biotic stress caused by infection with Cotton leafroll dwarf virus (CLRDV). The stability of mRNAs already described as reference genes in cotton was also assessed. The geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to select the best reference genes. In 8 of the 12 sets tested, miRNAs (miR172, 168 and 390) were found to be the best RGs. To validate the best selected RGs, miR159, miR164, miR2118, miR2910, miR3476, GhDCL2 and GhDCL4 expression levels were evaluated under biotic stress conditions, and miR164 and a putative myo-inositol oxigenase gene (GhMIOX) were assessed in leaves and flowers. The RGs selected in this work proved to be excellent reference genes in the two cases studied. Our results support the use of miRNAs as reference genes for miRNA and protein-coding genes.

Integrated mRNA and microRNA analysis identifies genes and small miRNA molecules associated with transcriptional and post-transcriptional-level responses to both drought stress and re-watering treatment in tobacco

BACKGROUND

Drought stress is one of the most severe problem limited agricultural productivity worldwide. It has been reported that plants response to drought-stress by sophisticated mechanisms at both transcriptional and post-transcriptional levels. However, the precise molecular mechanisms governing the responses of tobacco leaves to drought stress and water status are not well understood. To identify genes and miRNAs involved in drought-stress responses in tobacco, we performed both mRNA and small RNA sequencing on tobacco leaf samples from the following three treatments: untreated-control (CL), drought stress (DL), and re-watering (WL).

RESULTS

In total, we identified 798 differentially expressed genes (DEGs) between the DL and CL (DL vs. CL) treatments and identified 571 DEGs between the WL and DL (WL vs. DL) treatments. Further analysis revealed 443 overlapping DEGs between the DL vs. CL and WL vs. DL comparisons, and, strikingly, all of these genes exhibited opposing expression trends between these two comparisons, strongly suggesting that these overlapping DEGs are somehow involved in the responses of tobacco leaves to drought stress. Functional annotation analysis showed significant up-regulation of genes annotated to be involved in responses to stimulus and stress, (e.g., late embryogenesis abundant proteins and heat-shock proteins) antioxidant defense (e.g., peroxidases and glutathione S-transferases), down regulation of genes related to the cell cycle pathway, and photosynthesis processes. We also found 69 and 56 transcription factors (TFs) among the DEGs in, respectively, the DL vs. CL and the WL vs. DL comparisons. In addition, small RNA sequencing revealed 63 known microRNAs (miRNA) from 32 families and 368 novel miRNA candidates in tobacco. We also found that five known miRNA families (miR398, miR390, miR162, miR166, and miR168) showed differential regulation under drought conditions. Analysis to identify negative correlations between the differentially expressed miRNAs (DEMs) and DEGs revealed 92 mRNA-miRNA interactions between CL and DL plants, and 32 mRNA-miRNA interactions between DL and WL plants.

CONCLUSIONS

This study provides a global view of the transcriptional and the post-transcriptional responses of tobacco under drought stress and re-watering conditions. Our results establish an empirical foundation that should prove valuable for further investigations into the molecular mechanisms through which tobacco, and plants more generally, respond to drought stress at multiple molecular genetic levels.

Identification of Appropriate Reference Genes for Normalization of miRNA Expression in Grafted Watermelon Plants under Different Nutrient Stresses

Watermelon (Citrullus lanatus) is a globally important crop belonging to the family Cucurbitaceae. The grafting technique is commonly used to improve its tolerance to stress, as well as to enhance its nutrient uptake and utilization. It is believed that miRNA is most likely involved in its nutrient-starvation response as a graft-transportable signal. The quantitative real-time reverse transcriptase polymerase chain reaction is the preferred method for miRNA functional analysis, in which reliable reference genes for normalization are crucial to ensure the accuracy. The purpose of this study was to select appropriate reference genes in scion (watermelon) and rootstocks (squash and bottle gourd) of grafted watermelon plants under normal growth conditions and nutrient stresses (nitrogen and phosphorus starvation). Under nutrient starvation, geNorm identified miR167c and miR167f as two most stable genes in both watermelon leaves and squash roots. miR166b was recommended by both geNorm and NormFinder as the best reference in bottle gourd roots under nutrient limitation. Expression of a new Cucurbitaceae miRNA, miR85, was used to validate the reliability of candidate reference genes under nutrient starvation. Moreover, by comparing several target genes expression in qRT-PCR analysis with those in RNA-seq data, miR166b and miR167c were proved to be the most suitable reference genes to normalize miRNA expression under normal growth condition in scion and rootstock tissues, respectively. This study represents the first comprehensive survey of the stability of miRNA reference genes in Cucurbitaceae and provides valuable information for investigating more accurate miRNA expression involving grafted watermelon plants.

Genome-scale mRNA and small RNA transcriptomic insights into initiation of citrus apomixis

Nucellar embryony (NE) is an adventitious form of apomixis common in citrus, wherein asexual embryos initiate directly from nucellar cells surrounding the embryo sac. NE enables the fixation of desirable agronomic traits and the production of clonal offspring of virus-free rootstock, but impedes progress in hybrid breeding. In spite of the great importance of NE in citrus breeding and commercial production, little is understood about the underlying molecular mechanisms. In this study, the stages of nucellar embryo initiation (NEI) were determined for two polyembryonic citrus cultivars via histological observation. To explore the genes and regulatory pathways involved in NEI, we performed mRNA-seq and sRNA-seq analyses of ovules immediately prior to and at stages during NEI in the two pairs of cultivars. A total of 305 differentially expressed genes (DEGs) were identified between the poly- and monoembryonic ovules. Gene ontology (GO) analysis revealed that several processes are significantly enriched based on DEGs. In particular, response to stress, and especially response to oxidative stress, was over-represented in polyembryonic ovules. Nearly 150 miRNAs, comprising ~90 conserved and ~60 novel miRNAs, were identified in the ovules of either cultivar pair. Only two differentially expressed miRNAs (DEMs) were identified, of which the novel miRN23-5p was repressed whereas the targets accumulated in the polyembryonic ovules. This integrated study on the transcriptional and post-transcriptional regulatory profiles between poly- and monoembryonic citrus ovules provides new insights into the mechanism of NE, which should contribute to revealing the regulatory mechanisms of plant apomixis.

Selection of suitable reference genes for quantitative RT-PCR normalization in the halophyte Halostachys caspica under salt and drought stress

The plants are always subjected to various environmental stress, because of plant sessile growth. qRT-PCR is a sensitive and reliable technology, and the normalization of target gene expression with suitable reference genes is very important for obtaining accurate data. Halostachys caspica is an extremely salt-tolerant halophyte belonging to Chenopodiaceae and a good candidate to explore the stress-physiological and molecular mechanism. To get truly the expression profiles of coding genes and miRNAs in H. caspica in response to salt and drought stress using qRT-PCR, suitable reference genes need to be confirmed. In this study, 10 candidate genes including ACT, UBC10, UBC13, TUB2, TUB3, EF1α, 5S rRNA, tRNA, U6 and miR1436 from H. caspica are chosen, and among them, the former nine are commonly used as internal control genes, and miR1436 with high sequence copies is no significant difference expression in high salinity-treated and untreated small RNA libraries of this species. The three softwares are used to analyze expression stability. The results showed that EF1α and TUB3 were the most stable under salt and drought stress, respectively, and UBC10 was the most constant aross all the samples with the both stressed combination. This work will benefit deep studies on abiotic tolerance in H. caspica.

High-throughput sequencing and degradome analysis reveal altered expression of miRNAs and their targets in a male-sterile cybrid pummelo (Citrus grandis)

Background

G1 + HBP is a male sterile cybrid line with nuclear genome from Hirado Buntan pummelo (C. grandis Osbeck) (HBP) and mitochondrial genome from “Guoqing No.1” (G1, Satsuma mandarin), which provides a good opportunity to study male sterility and nuclear-cytoplasmic cross talk in citrus. High-throughput sRNA and degradome sequencing were applied to identify miRNAs and their targets in G1 + HBP and its fertile type HBP during reproductive development.

Results

A total of 184 known miRNAs, 22 novel miRNAs and 86 target genes were identified. Some of the targets are transcription factors involved in floral development, such as auxin response factors (ARFs), SQUAMOSA promoter binding protein box (SBP-box), MYB, basic region-leucine zipper (bZIP), APETALA2 (AP2) and transport inhibitor response 1 (TIR1). Eight target genes were confirmed to be sliced by corresponding miRNAs using 5’ RACE technology. Based on the sequencing abundance, 42 differentially expressed miRNAs between sterile line G1 + HBP and fertile line HBP were identified. Differential expression of miRNAs and their target genes between two lines was validated by quantitative RT-PCR, and reciprocal expression patterns between some miRNAs and their targets were demonstrated. The regulatory mechanism of miR167a was investigated by yeast one-hybrid and dual-luciferase assays that one dehydrate responsive element binding (DREB) transcription factor binds to miR167a promoter and transcriptionally repress miR167 expression.

Conclusion

Our study reveals the altered expression of miRNAs and their target genes in a male sterile line of pummelo and highlights that miRNA regulatory network may be involved in floral bud development and cytoplasmic male sterility in citrus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2882-0) contains supplementary material, which is available to authorized users.

Validation of Suitable Reference Genes for Assessing Gene Expression of MicroRNAs in Lonicera japonica

MicroRNAs (miRNAs), which play crucial regulatory roles in plant secondary metabolism and responses to the environment, could be developed as promising biomarkers for different varieties and production areas of herbal medicines. However, limited information is available for miRNAs from Lonicera japonica, which is widely used in East Asian countries owing to various pharmaceutically active secondary metabolites. Selection of suitable reference genes for quantification of target miRNA expression through quantitative real-time (qRT)-PCR is important for elucidating the molecular mechanisms of secondary metabolic regulation in different tissues and varieties of L. japonica. For precise normalization of gene expression data in L. japonica, 16 candidate miRNAs were examined in three tissues, as well as 21 cultivated varieties collected from 16 production areas, using GeNorm, NormFinder, and RefFinder algorithms. Our results revealed combination of u534122 and u3868172 as the best reference genes across all samples. Their specificity was confirmed by detecting the cycling threshold (C t) value ranges in different varieties of L. japonica collected from diverse production areas, suggesting the use of these two reference miRNAs is sufficient for accurate transcript normalization with different tissues, varieties, and production areas. To our knowledge, this is the first report on validation of reference miRNAs in honeysuckle (Lonicera spp.). Restuls from this study can further facilitate discovery of functional regulatory miRNAs in different varieties of L. japonica.

Selection and Verification of Candidate Reference Genes for Mature MicroRNA Expression by Quantitative RT-PCR in the Tea Plant (Camellia sinensis)

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is a rapid and sensitive method for analyzing microRNA (miRNA) expression. However, accurate qRT-PCR results depend on the selection of reliable reference genes as internal positive controls. To date, few studies have identified reliable reference genes for differential expression analysis of miRNAs among tissues, and among experimental conditions in plants. In this study, three miRNAs and four non-coding small RNAs (ncRNA) were selected as reference candidates, and the stability of their expression was evaluated among different tissues and under different experimental conditions in the tea plant (Camellia sinensis) using the geNorm and NormFinder programs. It was shown that miR159a was the best single reference gene in the bud to the fifth leaf, 5S rRNA was the most suitable gene in different organs, miR6149 was the most stable gene when the leaves were attacked by Ectropis oblique and U4, miR5368n and miR159a were the best genes when the leaves were treated by methyl jasmonate (MeJA), salicylic acid (SA) and abscisic acid (ABA), respectively. Our results provide suitable reference genes for future investigations on miRNA functions in tea plants.

Identification of miRNAs and their targets by high-throughput sequencing and degradome analysis in cytoplasmic male-sterile line NJCMS1A and its maintainer NJCMS1B of soybean

BACKGROUND

Cytoplasmic male sterility (CMS) provides crucial breeding materials that facilitate hybrid seed production in various crops, and thus plays an important role in the study of hybrid vigor (heterosis), in plants. However, the CMS regulatory network in soybean remains unclear. MicroRNAs (miRNAs) play crucial roles in flower and pollen development by targeting genes that regulate their expression in plants. To identify the miRNAs and their targets that exist in the soybean CMS line NJCMS1A and its maintainer NJCMS1B, high-throughput sequencing and degradome analysis were conducted in this study.

RESULTS

Two small RNA libraries were constructed from the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. A total of 105 new miRNAs present on the other arm of known pre-miRNAs, 23 new miRNA members, 158 novel miRNAs and 160 high-confidence soybean miRNAs were identified using high-throughput sequencing. Among the identified miRNAs, 101 differentially expressed miRNAs with greater than two-fold changes between NJCMS1A and NJCMS1B were discovered. The different expression levels of selected miRNAs were confirmed by stem-loop quantitative real-time PCR. A degradome analysis showed that 856 targets were predicted to be targeted by 296 miRNAs, including a squamosa promoter-binding protein-like transcription factor family protein, a pentatricopeptide repeat-containing protein, and an auxin response factor, which were previously shown to be involved in floral organ or anther development in plants. Additionally, some targets, including a MADS-box transcription factor, NADP-dependent isocitrate dehydrogenase and NADH-ubiquinone oxidoreductase 24 kDa subunit, were identified, and they may have some relationship with the programmed cell death, reactive oxygen species accumulation and energy deficiencies, which might lead to soybean male sterility.

CONCLUSIONS

The present study is the first to use deep sequencing technology to identify miRNAs and their targets in the flower buds of the soybean CMS line NJCMS1A and its maintainer NJCMS1B. The results revealed that the miRNAs might participate in flower and pollen development, which could facilitate our understanding of the molecular mechanisms behind CMS in soybean.

Selection of Reference Genes for Normalization of MicroRNA Expression by RT-qPCR in Sugarcane Buds under Cold Stress

Sugarcane, accounting for 80% of world's sugar, originates in the tropics but is cultivated mainly in the subtropics. Therefore, chilling injury frequently occurs and results in serious losses. Recent studies in various plant species have established microRNAs as key elements in the post-transcriptional regulation of response to biotic and abiotic stresses including cold stress. Though, its accuracy is largely influenced by the use of reference gene for normalization, quantitative PCR is undoubtedly a popular method used for identification of microRNAs. For identifying the most suitable reference genes for normalizing miRNAs expression in sugarcane under cold stress, 13 candidates among 17 were investigated using four algorithms: geNorm, NormFinder, deltaCt, and Bestkeeper, and four candidates were excluded because of unsatisfactory efficiency and specificity. Verification was carried out using cold-related genes miR319 and miR393 in cold-tolerant and sensitive cultivars. The results suggested that miR171/18S rRNA and miR171/miR5059 were the best reference gene sets for normalization for miRNA RT-qPCR, followed by the single miR171 and 18S rRNA. These results can aid research on miRNA responses during sugarcane stress, and the development of sugarcane tolerant to cold stress. This study is the first report concerning the reference gene selection of miRNA RT-qPCR in sugarcane.

Comparative Analysis of miRNAs and Their Target Transcripts between a Spontaneous Late-Ripening Sweet Orange Mutant and Its Wild-Type Using Small RNA and Degradome Sequencing

Fruit ripening in citrus is not well-understood at the molecular level. Knowledge of the regulatory mechanism of citrus fruit ripening at the post-transcriptional level in particular is lacking. Here, we comparatively analyzed the miRNAs and their target genes in a spontaneous late-ripening mutant, "Fengwan" sweet orange (MT) (Citrus sinensis L. Osbeck), and its wild-type counterpart ("Fengjie 72-1," WT). Using high-throughput sequencing of small RNAs and RNA degradome tags, we identified 107 known and 21 novel miRNAs, as well as 225 target genes. A total of 24 miRNAs (16 known miRNAs and 8 novel miRNAs) were shown to be differentially expressed between MT and WT. The expression pattern of several key miRNAs and their target genes during citrus fruit development and ripening stages was examined. Csi-miR156k, csi-miR159, and csi-miR166d suppressed specific transcription factors (GAMYBs, SPLs, and ATHBs) that are supposed to be important regulators involved in citrus fruit development and ripening. In the present study, miRNA-mediated silencing of target genes was found under complicated and sensitive regulation in citrus fruit. The identification of miRNAs and their target genes provide new clues for future investigation of mechanisms that regulate citrus fruit ripening.

Identification of miRNAs involved in fruit ripening in Cavendish bananas by deep sequencing

BACKGROUND

MicroRNAs (miRNAs) are a family of non-coding small RNAs that play an important regulatory role in various biological processes. Previous studies have reported that miRNAs are closely related to the ripening process in model plants. However, the miRNAs that are closely involved in the banana fruit ripening process remain unknown.

METHODS

Here, we investigated the miRNA populations from banana fruits in response to ethylene or 1-MCP treatment using a deep sequencing approach and bioinformatics analysis combined with quantitative RT-PCR validation.

RESULTS

A total of 125 known miRNAs and 26 novel miRNAs were identified from three libraries. MiRNA profiling of bananas in response to ethylene treatment compared with 1-MCP treatment showed differential expression of 82 miRNAs. Furthermore, the differentially expressed miRNAs were predicted to target a total of 815 target genes. Interestingly, some targets were annotated as transcription factors and other functional proteins closely involved in the development and the ripening process in other plant species. Analysis by qRT-PCR validated the contrasting expression patterns between several miRNAs and their target genes.

CONCLUSIONS

The miRNAome of the banana fruit in response to ethylene or 1-MCP treatment were identified by high-throughput sequencing. A total of 82 differentially expressed miRNAs were found to be closely associated with the ripening process. The miRNA target genes encode transcription factors and other functional proteins, including SPL, APETALA2, EIN3, E3 ubiquitin ligase, β-galactosidase, and β-glucosidase. These findings provide valuable information for further functional research of the miRNAs involved in banana fruit ripening.

Comprehensive selection of reference genes for quantitative gene expression analysis during seed development in Brassica napus

MicroRNAs have higher expression stability than protein-coding genes in B. napus seeds and are therefore good reference genes for miRNA and mRNA RT-qPCR analysis. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) has become the "gold standard" to gain insight into function of genes. However, the accuracy of the technique depends on appropriate reference genes for quantification analysis in different experimental conditions. Accumulation of microRNAs (miRNAs) has also been studied by RT-qPCR, but there are no reference genes currently validated for normalization of Brassica napus miRNA expression data. In this study, we selected 43 B. napus miRNAs and 18 previously validated mRNA reference genes. The expression stability of the candidate reference genes was evaluated in different tissue samples (stages of seed development, flowers, and leaves) using geNorm, NormFinder, and RefFinder analysis. The best-ranked reference genes for expression studies during seed development (miR167-1_2, miR11-1, miR159-1 and miR168-1) were used to asses the expression of miR03-1. Since candidate miRNAs showed higher expression stability than protein-coding genes in most of the tested conditions, the expression profile of DGAT1 gene was compared when normalized by the four most stable miRNAs reference genes and by the four most stable mRNA reference genes. The expected expression pattern of DGAT1 during seed development was achieved with the use of miRNA as reference genes. In conclusion, the most stable miRNA reference genes can be employed in the normalization of RT-qPCR quantification of miRNAs and protein-coding genes. This work is the first to perform a comprehensive survey of the stability of miRNA reference genes in B. napus and provides guidelines to obtain more accurate RT-qPCR results in B. napus seeds studies.

Genomewide analysis of small RNAs in nonembryogenic and embryogenic tissues of citrus: microRNA- and siRNA-mediated transcript cleavage involved in somatic embryogenesis

Somatic embryogenesis (SE) is a process of somatic cells becoming dedifferentiated and generating embryos. SE has been widely used in biotechnology as a powerful way of regeneration and a model system for studying plant embryogenesis, but the controlling mechanisms of SE are far from clear. Here, we show the genomewide profiles of miRNAs/siRNAs and their target genes in nonembryogenic and embryogenic tissues of 'Valencia' sweet orange. By high-throughput sequencing (HTS) of small RNAs and RNA degradome tags, we identified 50 known and 45 novel miRNAs, 130 miniature inverted-repeat transposable elements (MITEs) derived, 94 other and 235 phased small interfering RNAs (siRNAs), as well as 203 target genes. The majority of the abundantly expressed miRNAs/siRNAs exhibit lower expression levels in embryogenic callus (EC) or during SE process than in nonembryogenic callus (NEC), which is supposed to derepress the target genes that are involved in development and stress response, thus to activate the biological processes required for cell differentiation. However, the conserved csi-miR156a/b, miR164b and 171c directed suppression of specific transcription factors (TFs) are supposed to inactivate the postembryonic growth thus to maintain normal SE. In this study, miRNA- and siRNA-mediated silencing of target genes was found under sophisticated regulation in citrus SE system; the enhancement effect of specific conserved miRNAs on SE was discussed, providing new clues for future investigation of mechanisms that control SE.

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

An efficient method for miRNA detection and localization in crop plants

microRNAs are a class of non-coding small RNAs (sRNAs) that are important regulators of gene expression at the post-transcriptional level by mRNA cleavage or translation inhibition. Another class of sRNAs are siRNAs, which also regulate gene expression but by causing DNA methylation. This epigenetic regulatory role has been observed for some miRNAs as well. The use of sRNAs allows the development of biotechnological applications in plants. To develop these types of applications, and to better understand the natural roles they play, it is important to be able to detect and to localize these sRNAs at the plant tissue level. Sometimes, in crop plants this can be challenging. Therefore, we developed a tissue printing hybridization protocol for easy and efficient detection of sRNAs and demonstrate this by the analysis of the spatio-temporal expression patterns of the miRNAs miR159 and miR164 in fruits of various crop plants. Moreover, we show the possibility to also detect the expression of miRNAs in fruit juice using a dot blot hybridization approach.

Identification of miRNAs involved in pear fruit development and quality

Background

MicroRNAs (miRNAs) are a class of small, endogenous RNAs that take part in regulating genes through mediating gene expressions at the post-transcriptional level in plants. Previous studies have reported miRNA identification in various plants ranging from model plants to perennial fruit trees. However, the role of miRNAs in pear (Pyrus bretschneideri) fruit development is not clear. Here, we investigated the miRNA profiles of pear fruits from different time stages during development with Illumina HiSeq 2000 platform and bioinformatics analysis. Quantitative real-time PCR was used to validate the expression levels of miRNAs.

Results

Both conserved and species-specific miRNAs in pear have been identified in this study. Total reads, ranging from 19,030,925 to 25,576,773, were obtained from six small RNA libraries constructed for different stages of fruit development after flowering. Comparative profiling showed that an average of 90 miRNAs was expressed with significant differences between various developmental stages. KEGG pathway analysis on 2,216 target genes of 188 known miRNAs and 1,127 target genes of 184 novel miRNAs showed that miRNAs are widely involved in the regulation of fruit development. Among these, a total of eleven miRNAs putatively participate in the pathway of lignin biosynthesis, nine miRNAs were identified to take part in sugar and acid metabolism, and MiR160 was identified to regulate auxin response factor.

Conclusion

Comparative analysis of miRNAomes during pear fruit development is presented, and miRNAs were proved to be widely involved in the regulation of fruit development and formation of fruit quality, for example through lignin synthesis, sugar and acid metabolism, and hormone signaling. Combined with computational analysis and experimental confirmation, the research contributes valuable information for further functional research of microRNA in fruit development for pear and other species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-953) contains supplementary material, which is available to authorized users.

Genome-wide comparison of microRNAs and their targeted transcripts among leaf, flower and fruit of sweet orange

BACKGROUND

In plants, microRNAs (miRNAs) regulate gene expression mainly at the post-transcriptional level. Previous studies have demonstrated that miRNA-mediated gene silencing pathways play vital roles in plant development. Here, we used a high-throughput sequencing approach to characterize the miRNAs and their targeted transcripts in the leaf, flower and fruit of sweet orange.

RESULTS

A total of 183 known miRNAs and 38 novel miRNAs were identified. An in-house script was used to identify all potential secondary siRNAs derived from miRNA-targeted transcripts using sRNA and degradome sequencing data. Genome mapping revealed that these miRNAs were evenly distributed across the genome with several small clusters, and 69 pre-miRNAs were co-localized with simple sequence repeats (SSRs). Noticeably, the loop size of pre-miR396c was influenced by the repeat number of CUU unit. The expression pattern of miRNAs among different tissues and developmental stages were further investigated by both qRT-PCR and RNA gel blotting. Interestingly, Csi-miR164 was highly expressed in fruit ripening stage, and was validated to target a NAC transcription factor. This study depicts a global picture of miRNAs and their target genes in the genome of sweet orange, and focused on the comparison among leaf, flower and fruit tissues.

CONCLUSIONS

This study provides a global view of miRNAs and their target genes in different tissue of sweet orange, and focused on the identification of miRNA involved in the regulation of fruit ripening. The results of this study lay a foundation for unraveling key regulators of orange fruit development and ripening on post-transcriptional level.