Identifying conserved and novel microRNAs in developing seeds of Brassica napus using deep sequencing

Identification of conserved miRNAs and their targets in Jatropha curcas: an in silico approach

BACKGROUND

MicroRNAs (miRNAs) are small endogenous RNAs with an approximate length of 18-22 nucleotides and involved in the regulation of gene expression in transcriptional or post-transcriptional levels. They were found to be associated with leaf morphogenesis, flowering time, vegetative phase change, and response to environmental cues in plants, where they act as a critical regulatory factor. The nature of high conservancy of plant miRNAs within the plant species made it possible to detect the conserved miRNAs by computational approaches. Expressed Sequence Tags (EST) based comparative genomic approaches provide advantages over wet lab approaches as it is convenient, easy to carry out and less time consuming. EST-based in silico approach can unravel new conserved miRNAs in plants, even when the complete genome sequence is not available.

RESULTS

To identify the novel miRNAs, a total of 46,865 ESTs from Jatropha curcas were searched for homology to all available 6746 mature miRNAs of plant eudicotyledons. Finally, we ended up with 12 novel miRNAs in Jatropha that range from 18 to 19 nucleotides where their respective precursor miRNAs had 54.11-71.76% (A + U) content. The putative miRNAs belong to 12 individual miRNA family and most of them have higher (A + U) content ranging from 47.36 to 77.77% than their respective miRNA homologs. Many of the target genes by the newly identified miRNAs were associated with plant growth and development, stress response, defense and hormone signaling, and oil synthesis pathways.

CONCLUSION

These findings have the potential to speed up miRNA identification and expand our understanding of miRNA functions in J. curcas.

Multiple Functions of MiRNAs in Brassica napus L

The worldwide climate changes every year due to global warming, waterlogging, drought, salinity, pests, and pathogens, impeding crop productivity. Brassica napus is one of the most important oil crops in the world, and rapeseed oil is considered one of the most health-beneficial edible vegetable oils. Recently, miRNAs have been found and confirmed to control the expression of targets under disruptive environmental conditions. The mechanism is through the formation of the silencing complex that mediates post-transcriptional gene silencing, which pairs the target mRNA and target cleavage and/or translation inhibition. However, the functional role of miRNAs and targets in B. napus is still not clarified. This review focuses on the current knowledge of miRNAs concerning development regulation and biotic and abiotic stress responses in B. napus. Moreover, more strategies for miRNA manipulation in plants are discussed, along with future perspectives, and the enormous amount of transcriptome data available provides cues for miRNA functions in B. napus. Finally, the construction of the miRNA regulatory network can lead to the significant development of climate change-tolerant B. napus through miRNA manipulation.

Clone and Function Verification of the OPR gene in Brassica napus Related to Linoleic Acid Synthesis

BACKGROUND

Fatty acid composition and content affect rapeseed oil quality. Fatty acid synthesis-related genes in rapeseed have been studied globally by researchers. Nevertheless, rapeseed oil is mainly composed of seven different fatty acids (FA), and each fatty acid was regulated by different genes. Furthermore, different FA affect each other, which needs continuous and in-depth research to obtain more clear results in Brassica napus.

RESULTS

In this paper, broad-scale miRNA expression profiles were constructed and 21 differentially expressed miRNAs were detected. GO enrichment analysis showed that most up-regulated proteins were involved in transcription factor activity and catalytic activity. KEGG pathway enrichment analysis indicated that 20 pathways involving 36 target genes were enriched, of which the bna00592 pathway may be involved in fatty acid metabolism. The results were verified using a quantitative real-time PCR (RT-qPCR) analysis, we found that the target gene of bna-miR156b > c > g was the OPR (12-oxo-phytodienoic acid reductase). Four copies of OPR gene were found, and the over-expression vectors (pCAMBIA1300-35 s-OPR and pCAMBIA1300-RNAi-OPR) were constructed to verify their functions. In T1 and T2 generation, the content of linoleic acid (LA) increased significantly in OE but deceased in OPRi.

CONCLUSIONS

This is the first study to provide four copies of the OPR gene that regulates LA metabolism, can be used for the molecular mechanism of LA and optimizing fatty acid profiles in oilseed for breeding programs.

Identification of miRNA-mRNA Regulatory Modules Involved in Lipid Metabolism and Seed Development in a Woody Oil Tree (Camellia oleifera)

Tea oil camellia (Camellia oleifera), an important woody oil tree, is a source of seed oil of high nutritional and medicinal value that is widely planted in southern China. However, there is no report on the identification of the miRNAs involved in lipid metabolism and seed development in the high- and low-oil cultivars of tea oil camellia. Thus, we explored the roles of miRNAs in the key periods of oil formation and accumulation in the seeds of tea oil camellia and identified miRNA–mRNA regulatory modules involved in lipid metabolism and seed development. Sixteen small RNA libraries for four development stages of seed oil biosynthesis in high- and low-oil cultivars were constructed. A total of 196 miRNAs, including 156 known miRNAs from 35 families, and 40 novel miRNAs were identified, and 55 significantly differentially expressed miRNAs were found, which included 34 upregulated miRNAs, and 21 downregulated miRNAs. An integrated analysis of the miRNA and mRNA transcriptome sequence data revealed that 10 miRNA–mRNA regulatory modules were related to lipid metabolism; for example, the regulatory modules of ath-miR858b–MYB82/MYB3/MYB44 repressed seed oil biosynthesis, and a regulation module of csi-miR166e-5p–S-ACP-DES6 was involved in the formation and accumulation of oleic acid. A total of 23 miRNA–mRNA regulatory modules were involved in the regulation of the seed size, such as the regulatory module of hpe-miR162a_L-2–ARF19, involved in early seed development. A total of 12 miRNA–mRNA regulatory modules regulating growth and development were identified, such as the regulatory modules of han-miR156a_L+1–SPL4/SBP2, promoting early seed development. The expression changes of six miRNAs and their target genes were validated using quantitative real-time PCR, and the targeting relationship of the cpa-miR393_R-1–AFB2 regulatory module was verified by luciferase assays. These data provide important theoretical values and a scientific basis for the genetic improvement of new cultivars of tea oil camellia in the future.

Comprehensive In Silico Characterization and Expression Profiling of TCP Gene Family in Rapeseed

TCP proteins are plant-specific transcription factors that have multipurpose roles in plant developmental procedures and stress responses. Therefore, a genome-wide analysis was performed to categorize the TCP genes in the rapeseed genome. In this study, a total of 80 BnTCP genes were identified in the rapeseed genome and grouped into two main classes (PCF and CYC/TB1) according to phylogenetic analysis. The universal evolutionary analysis uncovered that BnTCP genes had experienced segmental duplications and positive selection pressure. Gene structure and conserved motif examination presented that Class I and Class II have diverse intron-exon patterns and motifs numbers. Overall, nine conserved motifs were identified and varied from 2 to 7 in all TCP genes; and some of them were gene-specific. Mainly, Class II (PCF and CYC/TB1) possessed diverse structures compared to Class I. We identified four hormone- and four stress-related responsive cis-elements in the promoter regions. Moreover, 32 bna-miRNAs from 14 families were found to be targeting 21 BnTCPs genes. Gene ontology enrichment analysis presented that the BnTCP genes were primarily related to RNA/DNA binding, metabolic processes, transcriptional regulatory activities, etc. Transcriptome-based tissue-specific expression analysis showed that only a few genes (mainly BnTCP9, BnTCP22, BnTCP25, BnTCP48, BnTCP52, BnTCP60, BnTCP66, and BnTCP74) presented higher expression in root, stem, leaf, flower, seeds, and silique among all tested tissues. Likewise, qRT-PCR-based expression analysis exhibited that BnTCP36, BnTCP39, BnTCP53, BnTCP59, and BnTCP60 showed higher expression at certain time points under various hormones and abiotic stress conditions but not by drought and MeJA. Our results opened the new groundwork for future understanding of the intricate mechanisms of BnTCP in various developmental processes and abiotic stress signaling pathways in rapeseed.

Identification of miRNAs and their target genes in genic male sterility lines in Brassica napus by small RNA sequencing

BACKGROUND

Brassica napus is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been used successfully for rapeseed hybrid production in China. MicroRNAs (miRNAs) play crucial regulatory roles in various plant growth, development, and stress response processes. However, reports on miRNAs that regulate the pollen development of GMS lines in B. napus are few.

RESULTS

In this study, 12 small RNA and transcriptome libraries were constructed and sequenced for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines, namely, "6251AB" and "6284AB". At the same time, 12 small RNA and transcriptome libraries were also constructed and sequenced for the flower buds from the fertile and sterile lines of two dominant GMS (DGMS) lines, namely, "4001AB" and "4006AB". Based on the results, 46 known miRNAs, 27 novel miRNAs on the other arm of known pre-miRNAs, and 44 new conserved miRNAs were identified. Thirty-five pairs of novel miRNA-3p/miRNA-5p were found. Among all the identified miRNAs, fifteen differentially expressed miRNAs with over 1.5-fold change between flower buds of sterile and fertile lines were identified, including six differentially expressed miRNAs between "4001A" and "4001B", two differentially expressed miRNAs between "4006A" and "4006B", four differentially expressed miRNAs between "6251A" and "6251B", and ten differentially expressed miRNAs between "6284A" and "6284B". The correlation analysis of small RNA and transcriptome sequencing was conducted. And 257 candidate target genes were predicted for the 15 differentially expressed miRNAs. The results of 5' modified RACE indicated that BnaA09g48720D, BnaA09g11120D, and BnaCnng51960D were cleaved by bna-miR398a-3p, bna-miR158-3p and bna-miR159a, respectively. Among the differentially expressed miRNAs, miR159 was chosen to analyze its function. Overexpression of bna-miR159 in Arabidopsis resulted in decreased seed setting rate, and shortened siliques, illustrating that miR159 may regulate the fertility and silique development in rapeseed.

CONCLUSIONS

Our findings provide an overview of miRNAs that are potentially involved in GMS and pollen development. New information on miRNAs and their related target genes are provided to exploit the GMS mechanism and reveal the miRNA networks in B. napus.

Genome-Wide Characterization of Glutathione Peroxidase (GPX) Gene Family in Rapeseed (Brassica napus L.) Revealed Their Role in Multiple Abiotic Stress Response and Hormone Signaling

Plant glutathione peroxidases (GPXs) are the main enzymes in the antioxidant defense system that sustain H₂O₂ homeostasis and normalize plant reaction to abiotic stress conditions. To understand the major roles of the GPX gene family in rapeseed (Brassica napus L.), for the first time, a genome-wide study identified 25 BnGPX genes in the rapeseed genome. The phylogenetic analysis discovered that GPX genes were grouped into four major groups (Group I-Group IV) from rapeseed and three closely interrelated plant species. The universal investigation uncovered that the BnGPXs gene experienced segmental duplications and positive selection pressure. Gene structure and motifs examination recommended that most of the BnGPX genes demonstrated a comparatively well-maintained exon-intron and motifs arrangement within the identical group. Likewise, we recognized five hormones-, four stress-, and numerous light-reactive cis-elements in the promoters of BnGPXs. Five putative bna-miRNAs from two families were also prophesied, targeting six BnGPXs genes. Gene ontology annotation results proved the main role of BnGPXs in antioxidant defense systems, ROS, and response to stress stimulus. Several BnGPXs genes revealed boosted expression profiles in many developmental tissues/organs, i.e., root, seed, leaf, stem, flower, and silique. The qRT-PCR based expression profiling exhibited that two genes (BnGPX21 and BnGPX23) were suggestively up-regulated against different hormones (ABA, IAA, and MeJA) and abiotic stress (salinity, cold, waterlogging, and drought) treatments. In short, our discoveries provide a basis for additional functional studies on the BnGPX genes in future rapeseed breeding programs.

Genome-Wide Analysis and Expression Profile of Superoxide Dismutase (SOD) Gene Family in Rapeseed (Brassica napus L.) under Different Hormones and Abiotic Stress Conditions

Superoxide dismutase (SOD) is an important enzyme that acts as the first line of protection in the plant antioxidant defense system, involved in eliminating reactive oxygen species (ROS) under harsh environmental conditions. Nevertheless, the SOD gene family was yet to be reported in rapeseed (Brassica napus L.). Thus, a genome-wide investigation was carried out to identify the rapeseed SOD genes. The present study recognized 31 BnSOD genes in the rapeseed genome, including 14 BnCSDs, 11 BnFSDs, and six BnMSDs. Phylogenetic analysis revealed that SOD genes from rapeseed and other closely related plant species were clustered into three groups based on the binding domain with high bootstrap values. The systemic analysis exposed that BnSODs experienced segmental duplications. Gene structure and motif analysis specified that most of the BnSOD genes displayed a relatively well-maintained exon-intron and motif configuration within the same group. Moreover, we identified five hormones and four stress- and several light-responsive cis-elements in the promoters of BnSODs. Thirty putative bna-miRNAs from seven families were also predicted, targeting 13 BnSODs. Gene ontology annotation outcomes confirm the BnSODs role under different stress stimuli, cellular oxidant detoxification processes, metal ion binding activities, SOD activity, and different cellular components. Twelve BnSOD genes exhibited higher expression profiles in numerous developmental tissues, i.e., root, leaf, stem, and silique. The qRT-PCR based expression profiling showed that eight genes (BnCSD1, BnCSD3, BnCSD14, BnFSD4, BnFSD5, BnFSD6, BnMSD2, and BnMSD10) were significantly up-regulated under different hormones (ABA, GA, IAA, and KT) and abiotic stress (salinity, cold, waterlogging, and drought) treatments. The predicted 3D structures discovered comparable conserved BnSOD protein structures. In short, our findings deliver a foundation for additional functional investigations on the BnSOD genes in rapeseed breeding programs.

Genome-wide analysis and expression patterns of lipid phospholipid phospholipase gene family in Brassica napus L

BACKGROUND

Lipid phosphate phosphatases (LPP) are critical for regulating the production and degradation of phosphatidic acid (PA), an essential signaling molecule under stress conditions. Thus far, the LPP family genes have not been reported in rapeseed (Brassica napus L.).

RESULTS

In this study, a genome-wide analysis was carried out to identify LPP family genes in rapeseed that respond to different stress conditions. Eleven BnLPPs genes were identified in the rapeseed genome. Based on phylogenetic and synteny analysis, BnLPPs were classified into four groups (Group I-Group IV). Gene structure and conserved motif analysis showed that similar intron/exon and motifs patterns occur in the same group. By evaluating cis-elements in the promoters, we recognized six hormone- and seven stress-responsive elements. Further, six putative miRNAs were identified targeting three BnLPP genes. Gene ontology analysis disclosed that BnLPP genes were closely associated with phosphatase/hydrolase activity, membrane parts, phosphorus metabolic process, and dephosphorylation. The qRT-PCR based expression profiles of BnLPP genes varied in different tissues/organs. Likewise, several gene expression were significantly up-regulated under NaCl, PEG, cold, ABA, GA, IAA, and KT treatments.

CONCLUSIONS

This is the first report to describe the comprehensive genome-wide analysis of the rapeseed LPP gene family. We identified different phytohormones and abiotic stress-associated genes that could help in enlightening the plant tolerance against phytohormones and abiotic stresses. The findings unlocked new gaps for the functional verification of the BnLPP gene family during stresses, leading to rapeseed improvement.

Catalase (CAT) Gene Family in Rapeseed (Brassica napus L.): Genome-Wide Analysis, Identification, and Expression Pattern in Response to Multiple Hormones and Abiotic Stress Conditions

Catalase (CAT) is an antioxidant enzyme expressed by the CAT gene family and exists in almost all aerobic organisms. Environmental stresses induce the generation of reactive oxygen species (ROS) that eventually hinder plant growth and development. The CAT enzyme translates the hydrogen peroxide (H₂O₂) to water (H₂O) and reduce the ROS levels to shelter the cells’ death. So far, the CAT gene family has not been reported in rapeseed (Brassica napus L.). Therefore, a genome-wide comprehensive analysis was conducted to classify the CAT genes in the rapeseed genome. The current study identified 14 BnCAT genes in the rapeseed genome. Based on phylogenetic and synteny analysis, the BnCATs belong to four groups (Groups I–IV). A gene structure and conserved motif analysis showed that Group I, Group II, and Group IV possess almost the same intron/exon pattern, and an equal number of motifs, while Group III contains diverse structures and contain 15 motifs. By analyzing the cis-elements in the promoters, we identified five hormone-correlated responsive elements and four stress-related responsive elements. Further, six putative bna-miRNAs were also identified, targeting three genes (BnCAT4, BnCAT6, and BnCAT8). Gene ontology (GO) enrichment analysis showed that the BnCAT genes were largely related to cellular organelles, ROS response, stimulus response, stress response, and antioxidant enzymes. Almost 10 BnCAT genes showed higher expression levels in different tissues, i.e., root, leaf, stem, and silique. The expression analysis showed that BnCAT1–BnCAT3 and BnCAT11–BnCAT13 were significantly upregulated by cold, salinity, abscisic acid (ABA), and gibberellic acid (GA) treatment, but not by drought and methyl jasmonate (MeJA). Notably, most of the genes were upregulated by waterlogging stress, except BnCAT6, BnCAT9, and BnCAT10. Our results opened new windows for future investigations and provided insights into the CAT family genes in rapeseed.

Comprehensive In Silico Characterization and Expression Profiling of TCP Gene Family in Rapeseed

TCP proteins are plant-specific transcription factors that have multipurpose roles in plant developmental procedures and stress responses. Therefore, a genome-wide analysis was performed to categorize the TCP genes in the rapeseed genome. In this study, a total of 80 BnTCP genes were identified in the rapeseed genome and grouped into two main classes (PCF and CYC/TB1) according to phylogenetic analysis. The universal evolutionary analysis uncovered that BnTCP genes had experienced segmental duplications and positive selection pressure. Gene structure and conserved motif examination presented that Class I and Class II have diverse intron-exon patterns and motifs numbers. Overall, nine conserved motifs were identified and varied from 2 to 7 in all TCP genes; and some of them were gene-specific. Mainly, Class II (PCF and CYC/TB1) possessed diverse structures compared to Class I. We identified four hormone- and four stress-related responsive cis-elements in the promoter regions. Moreover, 32 bna-miRNAs from 14 families were found to be targeting 21 BnTCPs genes. Gene ontology enrichment analysis presented that the BnTCP genes were primarily related to RNA/DNA binding, metabolic processes, transcriptional regulatory activities, etc. Transcriptome-based tissue-specific expression analysis showed that only a few genes (mainly BnTCP9, BnTCP22, BnTCP25, BnTCP48, BnTCP52, BnTCP60, BnTCP66, and BnTCP74) presented higher expression in root, stem, leaf, flower, seeds, and silique among all tested tissues. Likewise, qRT-PCR-based expression analysis exhibited that BnTCP36, BnTCP39, BnTCP53, BnTCP59, and BnTCP60 showed higher expression at certain time points under various hormones and abiotic stress conditions but not by drought and MeJA. Our results opened the new groundwork for future understanding of the intricate mechanisms of BnTCP in various developmental processes and abiotic stress signaling pathways in rapeseed.

Identification and Characterization of microRNAs in the Developing Seed of Linseed Flax (Linum usitatissimum L.)

Seed development plays an important role during the life cycle of plants. Linseed flax is an oil crop and the seed is a key organ for fatty acids synthesis and storage. So it is important to understand the molecular mechanism of fatty acid biosynthesis during seed development. In this study, four small RNA libraries from early seeds at 5, 10, 20 and 30 days after flowering (DAF) were constructed and used for high-throughput sequencing to identify microRNAs (miRNAs). A total of 235 miRNAs including 114 known conserved miRNAs and 121 novel miRNAs were identified. The expression patterns of these miRNAs in the four libraries were investigated by bioinformatics and quantitative real-time polymerase chain reaction (qPCR) analysis. It was found that several miRNAs, including Lus-miRNA156a was significantly correlated with seed development process. In order to confirm the actual biological function of Lus-miRNA156a, over-expression vector was constructed and transformed to Arabidopsis. The phenotypes of homozygous transgenic lines showed decreasing of oil content and most of the fatty acid content in seeds as well as late flowering time. The results provided a clue that miRNA156a participating the fatty acid biosynthesis pathway and the detailed molecular mechanism of how it regulates the pathway needs to be further investigated.

Comparative Analysis of miRNA Expression Profiles Between Heat-Tolerant and Heat-Sensitive Genotypes of Flowering Chinese Cabbage Under Heat Stress Using High-Throughput Sequencing

Heat stress disturbs cellular homeostasis, thus usually impairs yield of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). MicroRNAs (miRNAs) play a significant role in plant responses to different stresses by modulating gene expression at the post-transcriptional level. However, the roles that miRNAs and their target genes may play in heat tolerance of flowering Chinese cabbage remain poorly characterized. The current study sequenced six small RNA libraries generated from leaf tissues of flowering Chinese cabbage collected at 0, 6, and 12 h after 38 °C heat treatment, and identified 49 putative novel miRNAs and 43 known miRNAs that differentially expressed between heat-tolerant and heat-sensitive flowering Chinese cabbage. Among them, 14 novel and nine known miRNAs differentially expressed only in the heat-tolerant genotype under heat-stress, therefore, their target genes including disease resistance protein TAO1-like, RPS6, reticuline oxidase-like protein, etc. might play important roles in enhancing heat-tolerance. Gene Ontology (GO) analysis revealed that targets of these differentially expressed miRNAs may play key roles in responses to temperature stimulus, cell part, cellular process, cell, membrane, biological regulation, binding, and catalytic activities. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified their important functions in signal transduction, environmental adaptation, global and overview maps, as well as in stress adaptation and in MAPK signaling pathways such as cell death. These findings provide insight into the functions of the miRNAs in heat stress tolerance of flowering Chinese cabbage.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

Identifying High Confidence microRNAs in the Developing Seeds of Jatropha curcas

MicroRNAs (miRNAs) are endogenously short noncoding regulatory RNAs implicated in plant development and physiology. Nine small RNA (sRNA) libraries from three typical seed developmental stages (young, intermediate, and mature) were generated by deep sequencing to identify the miRNAs of J. curcas, a potential oilseed crop for the production of renewable oil. Strict criteria were adopted to identify 93 high confidence miRNAs including 48 conserved miRNAs and 45 novel miRNAs. Target genes of these miRNAs were involved in a broad range of physiological functions, including gene expression regulation, primary & secondary metabolism, growth & development, signal transduction, and stress response. About one third (29 out of 93) miRNAs showed significant changes in expression levels during the seed developmental process, indicating that the miRNAs might regulate its targets by their changes of transcription levels in seed development. However, most miRNAs were found differentially expressed in the late stage of seed development, suggesting that miRNAs play more important roles in the stage when seed accumulating organic matters and suffering dehydration stress. This study presents the first large scale identification of high confidence miRNAs in the developing seeds of J. curcas.

Modulation of miRNA expression in natural populations of A. thaliana along a wide altitudinal gradient of Indian Himalayas

Plant populations growing along an altitudinal gradient are exposed to different environmental conditions. They are excellent resources to study regulatory mechanisms adopted by plants to respond to different environmental stresses. Regulation by miRNA is one of such strategies. Here, we report how different miRNAs are preferentially expressed in the three natural populations of A. thaliana originating from a wide altitudinal range. The expression level of miRNAs was mostly governed by temperature and radiation. Majority of the identified miRNAs expressed commonly in the three populations. However, 30 miRNAs expressed significantly at different level between the low and the high altitude populations. Most of these miRNAs regulate the genes associated with different developmental processes, abiotic stresses including UV, cold, secondary metabolites, etc. Further, the expression of miR397 and miR858 involved in lignin biosynthesis and regulation of secondary metabolites respectively, may be regulated by light intensity. A few miRNAs expressed at increasing level with the increase in the altitude of the site indicating environment driven tight regulation of these miRNAs. Further, several novel miRNAs and isomiR diversity specific to the Himalayas are reported which might have an adaptive advantage. To the best of our knowledge, this is the first report on miRNA expression from natural plant populations.

Physiological studies and genome-wide microRNA profiling of cold-stressed Brassica napus

Temperature extremes, including cold, adversely impact plant growth and development. Plant responses to cold stress (CS) are regulated at both transcriptional and post-transcriptional levels. MicroRNAs (miRNAs), small non-coding RNAs, are known to be involved in post-transcriptional regulation of various developmental processes and metal stress in Brassica napus L. (canola), however, their role in response to CS is largely unknown. In this study, changes in various physiological parameters and endogenous abundance of miRNAs were characterized in spring canola seedlings (DH12075) exposed to 4 °C for 0-48 h. Cold stress induced electrolyte leakage, increased the levels of malondialdheyde and antioxidant enzymes and reduced photosynthetic efficiency. Using small RNA sequencing, 70 known and 126 novel miRNAs were identified in CS leaf tissues and among these, 25 known and 104 novel miRNAs were differentially expressed. Quantitative real-time (qRT) PCR analysis of eight selected miRNAs confirmed their CS responsiveness. Furthermore, the expression of six out of eight miRNAs exhibited an opposite trend in a winter variety of canola, 'Mendel', when compared to 'DH12075'. This first study on the B. napus miRNAome provides a framework for further functional analysis of these miRNAs and their targets in response to CS which may contribute towards the future development of cold resilient crops.

Integrated mRNA, sRNA, and degradome sequencing reveal oilseed rape complex responses to Sclerotinia sclerotiorum (Lib.) infection

Sclerotinia stem rot (SSR), caused by the fungal pathogen Sclerotinia sclerotiorum, is a devastating disease resulting in yield losses and decreases in seed quality in oilseed rape (Brassica napus) worldwide. However, the molecular mechanisms underlying the response of oilseed rape to S. sclerotiorum infection at the transcriptional and post-transcriptional levels are poorly understood. Here, we used an integrated omics approach (transcriptome, sRNAome, and degradome sequencing) on the Illumina platform to compare the RNA expression and post-transcriptional profiles of oilseed rape plants inoculated or not with S. sclerotiorum. In total, 7,065 differentially expressed genes (DEGs) compared with the mock-inoculated control at 48 hours post inoculation were identified. These DEGs were associated with protein kinases, signal transduction, transcription factors, hormones, pathogenesis-related proteins, secondary metabolism, and transport. In the sRNA-Seq analysis, 77 known and 176 novel miRNAs were identified; however, only 10 known and 41 novel miRNAs were differentially expressed between the samples inoculated or not with S. sclerotiorum. Degradome sequencing predicted 80 cleavage sites with 64 miRNAs. Integrated mRNA, sRNA and degradome sequencing analysis reveal oilseed rape complex responses to S. sclerotiorum infection. This study provides a global view of miRNA and mRNA expression profiles in oilseed rape following S. sclerotiorum infection.

Integrated microRNA and mRNA expression profiling reveals a complex network regulating pomegranate (Punica granatum L.) seed hardness

The breeding of new soft-seeded pomegranate cultivars provides new products for the market and increases farmers' incomes, yet the genetic architecture mediating seed hardness is largely unknown. Here, the seed hardness and hundred-seed weights of 26 cultivars were determined in 2 successive years. We conducted miRNA and mRNA sequencing to analyse the seeds of two varieties of Punica granatum: soft-seeded Tunisia and hard-seeded Sanbai, at 60 and 120 d after flowering. Seed hardness was strongly positively correlated with hundred-seed weight. We detected 25 and 12 differentially expressed miRNA-mRNA pairs with negative regulatory relationships between the two genotypes at 60 and 120 d after flowering, respectively. These miRNA-mRNA pairs mainly regulated seed hardness by altering cell wall structure. Transcription factors including NAC1, WRKY and MYC, which are involved in seed hardness, were targeted by differentially expressed mdm-miR164e and mdm-miR172b. Thus, seed hardness is the result of a complex biological process regulated by a miRNA-mRNA network in pomegranate. These results will help us understand the complexity of seed hardness and help to elucidate the miRNA-mediated molecular mechanisms that contribute to seed hardness in pomegranate.

Identification of microRNAs involved in lipid biosynthesis and seed size in developing sea buckthorn seeds using high-throughput sequencing

Sea buckthorn is a plant of medicinal and nutritional importance owing in part to the high levels of essential fatty acids, linoleic (up to 42%) and α-linolenic (up to 39%) acids in the seed oil. Sea buckthorn can produce seeds either via the sexual pathway or by apomixis. The seed development and maturation programs are critically dependent on miRNAs. To understand miRNA-mediated regulation of sea buckthorn seed development, eight small RNA libraries were constructed for deep sequencing from developing seeds of a low oil content line ‘SJ1’ and a high oil content line ‘XE3’. High-throughput sequencing identified 137 known miRNA from 27 families and 264 novel miRNAs. The potential targets of the identified miRNAs were predicted based on sequence homology. Nineteen (four known and 15 novel) and 22 (six known and 16 novel) miRNAs were found to be involved in lipid biosynthesis and seed size, respectively. An integrated analysis of mRNA and miRNA transcriptome and qRT-PCR identified some key miRNAs and their targets (miR164d-ARF2, miR168b-Δ9D, novelmiRNA-108-ACC, novelmiRNA-23-GPD1, novelmiRNA-58-DGAT1, and novelmiRNA-191-DGAT2) potentially involved in seed size and lipid biosynthesis of sea buckthorn seed. These results indicate the potential importance of miRNAs in regulating lipid biosynthesis and seed size in sea buckthorn.

Integration of omics approaches to understand oil/protein content during seed development in oilseed crops

Oilseed crops, especially soybean (Glycine max) and canola/rapeseed (Brassica napus), produce seeds that are rich in both proteins and oils and that are major sources of energy and nutrition worldwide. Most of the nutritional content in the seed is accumulated in the embryo during the seed filling stages of seed development. Understanding the metabolic pathways that are active during seed filling and how they are regulated are essential prerequisites to crop improvement. In this review, we summarize various omics studies of soybean and canola/rapeseed during seed filling, with emphasis on oil and protein traits, to gain a systems-level understanding of seed development. Currently, most (80-85%) of the soybean and rapeseed reference genomes have been sequenced (950 and 850 megabases, respectively). Parallel to these efforts, extensive omics datasets from different seed filling stages have become available. Transcriptome and proteome studies have detected preponderance of starch metabolism and glycolysis enzymes to be the possible cause of higher oil in B. napus compared to other crops. Small RNAome studies performed during the seed filling stages have revealed miRNA-mediated regulation of transcription factors, with the suggestion that this interaction could be responsible for transitioning the seeds from embryogenesis to maturation. In addition, progress made in dissecting the regulation of de novo fatty acid synthesis and protein storage pathways is described. Advances in high-throughput omics and comprehensive tissue-specific analyses make this an exciting time to attempt knowledge-driven investigation of complex regulatory pathways.

Identification and profiling of conserved and novel microRNAs involved in oil and oleic acid production during embryogenesis in Carya cathayensis Sarg

MicroRNAs (miRNAs) are important regulators of plant development and fruit formation. Mature embryos of hickory (Carya cathayensis Sarg.) nuts contain more than 70% oil (comprising 90% unsaturated fatty acids), along with a substantial amount of oleic acid. To understand the roles of miRNAs involved in oil and oleic acid production during hickory embryogenesis, three small RNA libraries from different stages of embryogenesis were constructed. Deep sequencing of these three libraries identified 95 conserved miRNAs with 19 miRNA*s, 7 novel miRNAs (as well as their corresponding miRNA*s), and 26 potentially novel miRNAs. The analysis identified 15 miRNAs involved in oil and oleic acid production that are differentially expressed during embryogenesis in hickory. Among them, nine miRNA sequences, including eight conserved and one novel, were confirmed by qRT-PCR. In addition, 145 target genes of the novel miRNAs were predicted using a bioinformatic approach. Our results provide a framework for better understanding the roles of miRNAs during embryogenesis in hickory.

Small RNA changes in synthetic Brassica napus

Small RNAs and microRNAs were found to vary extensively in synthetic Brassica napus and subsequent generations, accompanied by the activation of transposable elements in response to hybridization and polyploidization. Resynthesizing B. napus by hybridization and chromosome doubling provides an approach to create novel polyploids and increases the usable genetic variability in oilseed rape. Although many studies have shown that small RNAs (sRNAs) act as important factor during hybridization and polyploidization in plants, much less is known on how sRNAs change in synthetic B. napus, particularly in subsequent generations after formation. We performed high-throughput sequencing of sRNAs in S1-S4 generations of synthetic B. napus and in the homozygous B. oleracea and B. rapa parent lines. We found that the number of small RNAs (sRNAs) and microRNAs (miRNAs) doubled in synthetic B. napus relative to the parents. The proportions of common sRNAs detected varied from the S1 to S4 generations, suggesting sRNAs are unstable in synthetic B. napus. The majority of miRNAs (67.2 %) were non-additively expressed in the synthesized Brassica allotetraploid, and 33.3 % of miRNAs were novel in the resynthesized B. napus. The percentage of miRNAs derived from transposable elements (TEs) also increased, indicating transposon activation and increased transposon-associated miRNA production in response to hybridization and polyploidization. The number of target genes for each miRNA in the synthesized Brassica allotetraploid was doubled relative to the parents, enhancing the complexity of gene expression regulation. The potential roles of miRNAs and their targets are discussed. Our data demonstrate generational changes in sRNAs and miRNAs in synthesized B. napus.

Tight regulation of the interaction between Brassica napus and Sclerotinia sclerotiorum at the microRNA level

MicroRNAs (miRNAs) are multifunctional non-coding short nucleotide molecules. Nevertheless, the role of miRNAs in the interactions between plants and necrotrophic pathogens is largely unknown. Here, we report the identification of the miRNA repertoire of the economically important oil crop oilseed rape (Brassica napus) and those involved in interacting with its most devastating necrotrophic pathogen Sclerotinia sclerotiorum. We identified 280 B. napus miRNA candidates, including 53 novel candidates and 227 canonical members or variants of known miRNA families, by high-throughput deep sequencing of small RNAs from both normal and S. sclerotiorum-inoculated leaves. Target genes of 15 novel candidates and 222 known miRNAs were further identified by sequencing of degradomes from the two types of samples. MiRNA microarray analysis revealed that 68 miRNAs were differentially expressed between S. sclerotiorum-inoculated and uninoculated leaves. A set of these miRNAs target genes involved in plant defense to S. sclerotiorum and/or other pathogens such as nucleotide binding site-leucine-rich repeat (NBS-LRR) R genes and nitric oxygen and reactive oxygen species related genes. Additionally, three miRNAs target AGO1 and AGO2, key components of post-transcriptional gene silencing (PTGS). Expression of several viral PTGS suppressors reduced resistance to S. sclerotiorum. Arabidopsis mutants of AGO1 and AGO2 exhibited reduced resistance while transgenic lines over-expressing AGO1 displayed increased resistance to S. sclerotiorum in an AGO1 expression level-dependent manner. Moreover, transient over-expression of miRNAs targeting AGO1 and AGO2 decreased resistance to S. sclerotiorum in oilseed rape. Our results demonstrate that the interactions between B. napus and S. sclerotiorum are tightly regulated at miRNA level and probably involve PTGS.

Molecular and systems approaches towards drought-tolerant canola crops

1169 I. 1170 II. 1170 III. 1172 IV. 1176 V. 1181 VI. 1182 1183 References 1183 SUMMARY: Modern agriculture is facing multiple challenges including the necessity for a substantial increase in production to meet the needs of a burgeoning human population. Water shortage is a deleterious consequence of both population growth and climate change and is one of the most severe factors limiting global crop productivity. Brassica species, particularly canola varieties, are cultivated worldwide for edible oil, animal feed, and biodiesel, and suffer dramatic yield loss upon drought stress. The recent release of the Brassica napus genome supplies essential genetic information to facilitate identification of drought-related genes and provides new information for agricultural improvement in this species. Here we summarize current knowledge regarding drought responses of canola, including physiological and -omics effects of drought. We further discuss knowledge gained through translational biology based on discoveries in the closely related reference species Arabidopsis thaliana and through genetic strategies such as genome-wide association studies and analysis of natural variation. Knowledge of drought tolerance/resistance responses in canola together with research outcomes arising from new technologies and methodologies will inform novel strategies for improvement of drought tolerance and yield in this and other important crop species.

Dysregulation of miRNA isoform level at 5' end in Alzheimer's disease

Alzheimer's disease (AD) is the most common form of dementia, whose mechanism is still not yet fully understood. A miRNA-based signature method, commonly according to the changes of expression levels, is widely used for AD analysis in previous studies. Recently, miRNA isoforms called as isomiR variants, which is considered to play important biological roles, have been demonstrated as the applications of high throughput sequencing platforms. Here, we presented an entropy-based model to detect the miRNA isoform level at the 5' end, and found many miRNAs with significant changes of isoform levels between the early stage and the late stage of AD by the application of this model to the public data. The statistical significance of the overlap between isoform-level changed miRNAs and AD related miRNAs extracted from HMDD2 supports that these miRNA isoforms are not degradation products. Based on the most common isomiR seed analysis of isoform-level changed AD related miRNAs, the predicted targets are also found to be enriched for genes involved in transcriptional regulation and the nervous system. After comparing with the expression level based method, we detected that changes of 5' isoform levels are more stable than those of expression levels for AD related miRNA detecting.

Identification and Expression Analyses of miRNAs from Two Contrasting Flower Color Cultivars of Canna by Deep Sequencing

miRNAs are endogenous small RNA (sRNA) that play critical roles in plant development processes. Canna is an ornamental plant belonging to family Cannaceae. Here, we report for the first time the identification and differential expression of miRNAs in two contrasting flower color cultivars of Canna, Tropical sunrise and Red president. A total of 313 known miRNAs belonging to 78 miRNA families were identified from both the cultivars. Thirty one miRNAs (17 miRNA families) were specific to Tropical sunrise and 43 miRNAs (10 miRNA families) were specific to Red president. Thirty two and 18 putative new miRNAs were identified from Tropical sunrise and Red president, respectively. One hundred and nine miRNAs were differentially expressed in the two cultivars targeting 1343 genes. Among these, 16 miRNAs families targeting60 genes were involved in flower development related traits and five miRNA families targeting five genes were involved in phenyl propanoid and pigment metabolic processes. We further validated the expression analysis of a few miRNA and their target genes by qRT-PCR. Transcription factors were the major miRNA targets identified. Target validation of a few randomly selected miRNAs by RLM-RACE was performed but was successful with only miR162. These findings will help in understanding flower development processes, particularly the color development in Canna.

Identification of miRNAs Affecting the Establishment of Brassica Alboglabra Seedling

MicroRNAs (miRNAs) are important for plant development including seed formation, dormancy, and germination, as well as seedling establishment. The Brassica vegetable seedling establishment stage influences the development of high quality seedlings, but also affects the nutrient content of sprouts. Chinese kale (Brassica alboglabra) seedlings at different growth stages were used to construct two small-RNA (sRNA) libraries. We comprehensively analyzed the miRNAs in 2- and 9-day-old seedlings. An average of 11,722,490 clean reads were generated after removing low-quality reads and adapter contaminants. The results revealed that 37.65 and 26.69% of the sRNAs in 2- and 9-day-old seedlings, respectively, were 24 nt long. In total, 254 known mature miRNA sequences from 228 miRNA families and 343 novel miRNAs were identified. Of these miRNAs, 224 were differentially expressed between the two analyzed libraries. The most abundant miRNAs identified by sequence homology were miR156, miR167, and miR157, each with more than 100,000 sequenced reads. Compared with the expression levels in 2-day-old seedlings, MiR8154 and miR390 were the most up- and down-regulated miRNAs respectively in 9-day-old seedlings. Gene ontology enrichment analysis of the differentially expressed-miRNA target genes affecting biological processes revealed that most genes were in the "regulation of transcription" category. Additionally, the expression patterns of some miRNAs and target genes were validated by quantitative real-time polymerase chain reaction. We determined that development-associated miRNAs (e.g., bal-miR156/157/159/166/167/172/396), were highly-expressed during seedling-establishment stage, as were stress-related (bal-miR408) and metabolism-related (bal-miR826) miRNAs. Combined with the low level of targets SPL9 and AP2, it was concluded that miR156-SPL9 and miR172-AP modules play key roles during the B. alboglabra seedling establishment stage.

Identification of Rapeseed MicroRNAs Involved in Early Stage Seed Germination under Salt and Drought Stresses

Drought and salinity are severe and wide-ranging abiotic stresses that substantially affect crop germination, development and productivity, and seed germination is the first critical step in plant growth and development. To comprehensively investigate small-RNA targets and improve our understanding of miRNA-mediated post-transcriptional regulation networks during Brassica napus seed imbibition under drought and salt stresses, we constructed three small-RNA libraries from B. napus variety ZS11 embryos exposed to salt (200 mM NaCl, denoted "S"), drought (200 g L(-1) PEG-6000, denoted "D"), and distilled water (denoted "CK") during imbibition and sequenced them using an Illumina Genome Analyzer. A total of 11,528,557, 12,080,081, and 12,315,608 raw reads were obtained from the CK, D, and S libraries, respectively. Further analysis identified 85 known miRNAs belonging to 31 miRNA families and 882 novel miRNAs among the three libraries. Comparison of the D and CK libraries revealed significant down-regulation of six miRNA families, miR156, miR169, miR860, miR399, miR171, and miR395, whereas only miR172 was significantly up-regulated. In contrast, comparison of the S library with the CK library showed significant down-regulation of only two miRNA families: miRNA393 and miRNA399. Putative targets for 336, 376, and 340 novel miRNAs were successfully predicted in the CK, D, and S libraries, respectively, and 271 miRNA families and 20 target gene families [including disease resistance protein (DIRP), drought-responsive family protein (DRRP), early responsive to dehydration stress protein (ERD), stress-responsive alpha-beta barrel domain protein (SRAP), and salt tolerance homolog2 (STH2)] were confirmed as being core miRNAs and genes involved in the seed imbibition response to salt and drought stresses. The sequencing results were partially validated by quantitative RT-PCR for both conserved and novel miRNAs as well as the predicted target genes. Our data suggest that diverse and complex miRNAs are involved in seed imbibition, indicating that miRNAs are involved in plant hormone regulation, and may play important roles during seed germination under salt- or drought-stress conditions.

Identification of microRNAs Actively Involved in Fatty Acid Biosynthesis in Developing Brassica napus Seeds Using High-Throughput Sequencing

Seed development has a critical role during the spermatophyte life cycle. In Brassica napus, a major oil crop, fatty acids are synthesized and stored in specific tissues during embryogenesis, and understanding the molecular mechanism underlying fatty acid biosynthesis during seed development is an important research goal. In this study, we constructed three small RNA libraries from early seeds at 14, 21, and 28 days after flowering (DAF) and used high-throughput sequencing to examine microRNA (miRNA) expression. A total of 85 known miRNAs from 30 families and 1160 novel miRNAs were identified, of which 24, including 5 known and 19 novel miRNAs, were found to be involved in fatty acid biosynthesis.bna-miR156b, bna-miR156c, bna-miR156g, novel_mir_1706, novel_mir_1407, novel_mir_173, and novel_mir_104 were significantly down-regulated at 21 DAF and 28 DAF, whereas bna-miR159, novel_mir_1081, novel_mir_19 and novel_mir_555 were significantly up-regulated. In addition, we found that some miRNAs regulate functional genes that are directly involved in fatty acid biosynthesis and that other miRNAs regulate the process of fatty acid biosynthesis by acting on a large number of transcription factors. The miRNAs and their corresponding predicted targets were partially validated by quantitative RT-PCR. Our data suggest that diverse and complex miRNAs are involved in the seed development process and that miRNAs play important roles in fatty acid biosynthesis during seed development.

Identification, evolution, and expression partitioning of miRNAs in allopolyploid Brassica napus

The recently published genome of Brassica napus offers for the first time the opportunity to gain insights into the genomic organization and the evolution of miRNAs in oilseed rape. In this study, 12 small RNA libraries from two B. napus cultivars (Tapidor and Ningyou7) and their four double-haploid lines were sequenced, employing the newly sequenced B. napus genome, together with genomes of its progenitors Brassica rapa and Brassica oleracea. A total of 645 miRNAs including 280 conserved and 365 novel miRNAs were identified. Comparative analysis revealed a high level of genomic conservation of MIRNAs (75.9%) between the subgenomes of B. napus and its two progenitors' genomes, and MIRNA lost/gain events (133) occurred in B. napus after its speciation. Furthermore, significant partitioning of miRNA expressions between the two subgenomes in B. napus was detected. The data of degradome sequencing, miRNA-mediated cleavage, and expression analyses support specific interactions between miRNAs and their targets in the modulation of diverse physiological processes in roots and leaves, as well as in biosynthesis of, for example, glucosinolates and lipids in oilseed rape. These data provide a first genome-wide view on the origin, evolution, and genomic organization of B. napus MIRNAs.

ASR5 is involved in the regulation of miRNA expression in rice

The work describes an ASR knockdown transcriptomic analysis by deep sequencing of rice root seedlings and the transactivation of ASR cis-acting elements in the upstream region of a MIR gene. MicroRNAs are key regulators of gene expression that guide post-transcriptional control of plant development and responses to environmental stresses. ASR (ABA, Stress and Ripening) proteins are plant-specific transcription factors with key roles in different biological processes. In rice, ASR proteins have been suggested to participate in the regulation of stress response genes. This work describes the transcriptomic analysis by deep sequencing two libraries, comparing miRNA abundance from the roots of transgenic ASR5 knockdown rice seedlings with that of the roots of wild-type non-transformed rice seedlings. Members of 59 miRNA families were detected, and 276 mature miRNAs were identified. Our analysis detected 112 miRNAs that were differentially expressed between the two libraries. A predicted inverse correlation between miR167abc and its target gene (LOC_Os07g29820) was confirmed using RT-qPCR. Protoplast transactivation assays showed that ASR5 is able to recognize binding sites upstream of the MIR167a gene and drive its expression in vivo. Together, our data establish a comparative study of miRNAome profiles and is the first study to suggest the involvement of ASR proteins in miRNA gene regulation.

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.

Altered Fruit and Seed Development of Transgenic Rapeseed (Brassica napus) Over-Expressing MicroRNA394

Fruit and seed development in plants is a complex biological process mainly involved in input and biosynthesis of many storage compounds such as proteins and oils. Although the basic biochemical pathways for production of the storage metabolites in plants are well characterized, their regulatory mechanisms are not fully understood. In this study, we functionally identified rapeseed (Brassica napus) miR394 with its target gene Brassica napus LEAF CURLING RESPONSIVENESS (BnLCR) to dissect a role of miR394 during the fruit and seed development. Transgenic rapeseed plants over-expressing miR394 under the control of the cauliflower mosaic virus 35S promoter were generated. miR394 over-expression plants exhibited a delayed flowering time and enlarged size of plants, leaf blade, pods and seed body, but developed seeds with higher contents of protein and glucosinolates (GLS) and lower levels of oil accumulation as compared to wild-type. Over-expression of miR394 altered the fatty acid (FA) composition by increasing several FA species such as C16:0 and C18:0 and unsaturated species of C20:1 and C22:1 but lowering C18:3. This change was accompanied by induction of genes coding for transcription factors of FA synthesis including LEAFY COTYLEDON1 (BnLEC1), BnLEC2, and FUSCA3 (FUS3). Because the phytohormone auxin plays a crucial role in fruit development and seed patterning, the DR5-GUS reporter was used for monitoring the auxin response in Arabidopsis siliques and demonstrated that the DR5 gene was strongly expressed. These results suggest that BnmiR394 is involved in rapeseed fruit and seed development.

Identification of microRNAs and transcript targets in Camelina sativa by deep sequencing and computational methods

Camelina sativa is an annual oilseed crop that is under intensive development for renewable resources of biofuels and industrial oils. MicroRNAs, or miRNAs, are endogenously encoded small RNAs that play key roles in diverse plant biological processes. Here, we conducted deep sequencing on small RNA libraries prepared from camelina leaves, flower buds and two stages of developing seeds corresponding to initial and peak storage products accumulation. Computational analyses identified 207 known miRNAs belonging to 63 families, as well as 5 novel miRNAs. These miRNAs, especially members of the miRNA families, varied greatly in different tissues and developmental stages. The predicted miRNA target genes are involved in a broad range of physiological functions including lipid metabolism. This report is the first step toward elucidating roles of miRNAs in C. sativa and will provide additional tools to improve this oilseed crop for biofuels and biomaterials.

Entropy-based model for miRNA isoform analysis

MiRNAs have been widely studied due to their important post-transcriptional regulatory roles in gene expression. Many reports have demonstrated the evidence of miRNA isoform products (isomiRs) in high-throughput small RNA sequencing data. However, the biological function involved in these molecules is still not well investigated. Here, we developed a Shannon entropy-based model to estimate isomiR expression profiles of high-throughput small RNA sequencing data extracted from miRBase webserver. By using the Kolmogorov-Smirnov statistical test (KS test), we demonstrated that the 5p and 3p miRNAs present more variants than the single arm miRNAs. We also found that the isomiR variant, except the 3’ isomiR variant, is strongly correlated with Minimum Free Energy (MFE) of pre-miRNA, suggesting the intrinsic feature of pre-miRNA should be one of the important factors for the miRNA regulation. The functional enrichment analysis showed that the miRNAs with high variation, particularly the 5’ end variation, are enriched in a set of critical functions, supporting these molecules should not be randomly produced. Our results provide a probabilistic framework for miRNA isoforms analysis, and give functional insights into pre-miRNA processing.

Transcriptome analysis of canola (Brassica napus) under salt stress at the germination stage

Canola (Brassica napus) is one of the most important oil crops in the world. However, its yield has been constrained by salt stress. In this study, transcriptome profiles were explored using Digital Gene Expression (DGE) at 0, 3, 12 and 24 hours after H2O (control) and NaCl treatments on B. napus roots at the germination stage. Comparisons of gene-expression between the control and the treatment were conducted after tag-mapping to the sequenced Brassica rapa genome. The differentially expressed genes during the time course of salt stress were focused on, and 163 genes were identified to be differentially expressed at all the time points. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that some of the genes were involved in proline metabolism, inositol metabolism, carbohydrate metabolic processes and oxidation-reduction processes and may play vital roles in the salt-stress response at the germination stage. Thus, this study provides new candidate salt stress responding genes, which may function in novel putative nodes in the molecular pathways of salt stress resistance.

Identification of miRNAs and their targets in transgenic Brassica napus and its acceptor (Westar) by high-throughput sequencing and degradome analysis

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that play many roles in plant growth, development, and the stress response.

A genome-wide identification of the miRNAome in response to salinity stress in date palm (Phoenix dactylifera L.)

Although date palm is relatively salt-tolerant, little is known about the underlying molecular mechanisms that contribute to its salt tolerance. Only recently, investigators have uncovered microRNA-mediated post-transcriptional gene regulation, which is critical for typical plant development and adaptation to stress conditions such as salinity. To identify conserved and novel miRNAs in date palm and to characterize miRNAs that could play a role in salt tolerance, we have generated sRNA libraries from the leaves and roots of NaCl-treated and untreated seedlings of date palm. Deep sequencing of these four sRNA libraries yielded approximately 251 million reads. The bioinformatics analysis has identified 153 homologs of conserved miRNAs, 89 miRNA variants, and 180 putative novel miRNAs in date palm. Expression profiles under salinity revealed differential regulation of some miRNAs in date palm. In leaves, 54 of the identified miRNAs were significantly affected and the majority (70%) of them were upregulated, whereas in roots, 25 of the identified miRNAs were significantly affected and 76% of them were upregulated by the salinity stress. The salt-responsiveness of some of these miRNAs was further validated using semi-quantitative PCR (qPCR). Some of the predicted targets for the identified miRNA include genes with known functions in plant salt tolerance, such as potassium channel AKT2-like proteins, vacuolar protein sorting-associated protein, calcium-dependent and mitogen-activated proteins. As one of the first cultivated trees in the world that can tolerate a wide range of abiotic stresses, date palm contains a large population of conserved and non-conserved miRNAs that function at the post-transcriptional level. This study provided insights into miRNA-mediated gene expression that are important for adaptation to salinity in date palms.

The elucidation of stress memory inheritance in Brassica rapa plants

Plants are able to maintain the memory of stress exposure throughout their ontogenesis and faithfully propagate it into the next generation. Recent evidence argues for the epigenetic nature of this phenomenon. Small RNAs (smRNAs) are one of the vital epigenetic factors because they can both affect gene expression at the place of their generation and maintain non-cell-autonomous gene regulation. Here, we have made an attempt to decipher the contribution of smRNAs to the heat-shock-induced transgenerational inheritance in Brassica rapa plants using sequencing technology. To do this, we have generated comprehensive profiles of a transcriptome and a small RNAome (smRNAome) from somatic and reproductive tissues of stressed plants and their untreated progeny. We have demonstrated that the highest tissue-specific alterations in the transcriptome and smRNAome profile are detected in tissues that were not directly exposed to stress, namely, in the endosperm and pollen. Importantly, we have revealed that the progeny of stressed plants exhibit the highest fluctuations at the smRNAome level but not at the transcriptome level. Additionally, we have uncovered the existence of heat-inducible and transgenerationally transmitted tRNA-derived small RNA fragments in plants. Finally, we suggest that miR168 and braAGO1 are involved in the stress-induced transgenerational inheritance in plants.

Identification and characterization of microRNAs in oilseed rape (Brassica napus) responsive to infection with the pathogenic fungus Verticillium longisporum using Brassica AA (Brassica rapa) and CC (Brassica oleracea) as reference genomes

Verticillium longisporum, a soil-borne pathogenic fungus, causes vascular disease in oilseed rape (Brassica napus). We proposed that plant microRNAs (miRNAs) are involved in the plant-V. longisporum interaction. To identify oilseed rape miRNAs, we deep-sequenced two small RNA libraries made from V. longisporum infected/noninfected roots and employed Brassica rapa and Brassica oleracea genomes as references for miRNA prediction and characterization. We identified 893 B. napus miRNAs representing 360 conserved and 533 novel miRNAs, and mapped 429 and 464 miRNAs to the AA and CC genomes, respectively. Microsynteny analysis with the conserved miRNAs and their flanking protein coding sequences revealed 137 AA-CC genome syntenic miRNA pairs and 61 AA and 42 CC genome-unique miRNAs. Sixty-two miRNAs were responsive to the V. longisporum infection. We present data for specific interactions and simultaneously reciprocal changes in the expression levels of the miRNAs and their targets in the infected roots. We demonstrate that miRNAs are involved in the plant-fungus interaction and that miRNA168-Argonaute 1 (AGO1) expression modulation might act as a key regulatory module in a compatible plant-V. longisporum interaction. Our results suggest that V. longisporum may have evolved a virulence mechanism by interference with plant miRNAs to reprogram plant gene expression and achieve infection.

Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding

The transcript pool of a plant part, under any given condition, is a collection of mRNAs that will pave the way for a biochemical reaction of the plant to stimuli. Over the past decades, transcriptome study has advanced from Northern blotting to RNA sequencing (RNA-seq), through other techniques, of which real-time quantitative polymerase chain reaction (PCR) and microarray are the most significant ones. The questions being addressed by such studies have also matured from a solitary process to expression atlas and marker-assisted genetic enhancement. Not only genes and their networks involved in various developmental processes of plant parts have been elucidated, but also stress tolerant genes have been highlighted. The transcriptome of a plant with altered expression of a target gene has given information about the downstream genes. Marker information has been used for breeding improved varieties. Fortunately, the data generated by transcriptome analysis has been made freely available for ample utilization and comparison. The review discusses this wide variety of transcriptome data being generated in plants, which includes developmental stages, abiotic and biotic stress, effect of altered gene expression, as well as comparative transcriptomics, with a special emphasis on microarray and RNA-seq. Such data can be used to determine the regulatory gene networks, which can subsequently be utilized for generating improved plant varieties.

A genome-wide perspective of miRNAome in response to high temperature, salinity and drought stresses in Brassica juncea (Czern) L

Micro RNAs (miRNAs) are involved in diverse biological processes including adaptive response towards abiotic stresses. To unravel small RNAs and more specifically miRNAs that can potentially regulate determinants of abiotic stress tolerance, next generation sequencing of B. juncea seedlings subjected to high temperature, high salt and drought conditions was carried out. With the help of UEA sRNA workbench software package, 51 conserved miRNAs belonging to 30 miRNA families were identified. As there was limited genomic information available for B. juncea, we generated and assembled its genome sequence at a low coverage. Using the generated sequence and other publically available Brassica genomic/transcriptomic resources as mapping reference, 126 novel (not reported in any plant species) were discovered for the first time in B. juncea. Further analysis also revealed existence of 32 and 37 star sequences for conserved and novel miRNAs, respectively. The expression of selected conserved and novel miRNAs under conditions of different abiotic stresses was revalidated through universal TaqMan based real time PCR. Putative targets of identified conserved and novel miRNAs were predicted in B. rapa to gain insights into functional roles manifested by B. juncea miRNAs. Furthermore, SPL2-like, ARF17-like and a NAC domain containing protein were experimentally validated as targets of miR156, miR160 and miR164 respectively. Investigation of gene ontologies linked with targets of known and novel miRNAs forecasted their involvement in various biological functions.

Identifying microRNAs and transcript targets in Jatropha seeds

MicroRNAs, or miRNAs, are endogenously encoded small RNAs that play a key role in diverse plant biological processes. Jatropha curcas L. has received significant attention as a potential oilseed crop for the production of renewable oil. Here, a sRNA library of mature seeds and three mRNA libraries from three different seed development stages were generated by deep sequencing to identify and characterize the miRNAs and pre-miRNAs of J. curcas. Computational analysis was used for the identification of 180 conserved miRNAs and 41 precursors (pre-miRNAs) as well as 16 novel pre-miRNAs. The predicted miRNA target genes are involved in a broad range of physiological functions, including cellular structure, nuclear function, translation, transport, hormone synthesis, defense, and lipid metabolism. Some pre-miRNA and miRNA targets vary in abundance between the three stages of seed development. A search for sequences that produce siRNA was performed, and the results indicated that J. curcas siRNAs play a role in nuclear functions, transport, catalytic processes and disease resistance. This study presents the first large scale identification of J. curcas miRNAs and their targets in mature seeds based on deep sequencing, and it contributes to a functional understanding of these miRNAs.

Identifying microRNAs and transcript targets in Jatropha seeds

MicroRNAs, or miRNAs, are endogenously encoded small RNAs that play a key role in diverse plant biological processes. Jatropha curcas L. has received significant attention as a potential oilseed crop for the production of renewable oil. Here, a sRNA library of mature seeds and three mRNA libraries from three different seed development stages were generated by deep sequencing to identify and characterize the miRNAs and pre-miRNAs of J. curcas. Computational analysis was used for the identification of 180 conserved miRNAs and 41 precursors (pre-miRNAs) as well as 16 novel pre-miRNAs. The predicted miRNA target genes are involved in a broad range of physiological functions, including cellular structure, nuclear function, translation, transport, hormone synthesis, defense, and lipid metabolism. Some pre-miRNA and miRNA targets vary in abundance between the three stages of seed development. A search for sequences that produce siRNA was performed, and the results indicated that J. curcas siRNAs play a role in nuclear functions, transport, catalytic processes and disease resistance. This study presents the first large scale identification of J. curcas miRNAs and their targets in mature seeds based on deep sequencing, and it contributes to a functional understanding of these miRNAs.

Applications and challenges of next-generation sequencing in Brassica species

Next-generation sequencing (NGS) produces numerous (often millions) short DNA sequence reads, typically varying between 25 and 400 bp in length, at a relatively low cost and in a short time. This revolutionary technology is being increasingly applied in whole-genome, transcriptome, epigenome and small RNA sequencing, molecular marker and gene discovery, comparative and evolutionary genomics, and association studies. The Brassica genus comprises some of the most agro-economically important crops, providing abundant vegetables, condiments, fodder, oil and medicinal products. Many Brassica species have undergone the process of polyploidization, which makes their genomes exceptionally complex and can create difficulties in genomics research. NGS injects new vigor into Brassica research, yet also faces specific challenges in the analysis of complex crop genomes and traits. In this article, we review the advantages and limitations of different NGS technologies and their applications and challenges, using Brassica as an advanced model system for agronomically important, polyploid crops. Specifically, we focus on the use of NGS for genome resequencing, transcriptome sequencing, development of single-nucleotide polymorphism markers, and identification of novel microRNAs and their targets. We present trends and advances in NGS technology in relation to Brassica crop improvement, with wide application for sophisticated genomics research into agronomically important polyploid crops.

High-throughput sequencing identification of novel and conserved miRNAs in the Brassica oleracea leaves

Background

Plant microRNAs are short (~21 nt) non-coding molecules that regulate gene expression by targeting the mRNA cleavage or protein translation inhibition. In this manner, they play many important roles in the cells of living organisms. One of the plant species in which the entire set of miRNAs has not been yet completely identified is Brassica oleracea var. capitata (cabbage). For this reason and for the economic and nutritional importance of this food crop, high-throughput small RNAs sequencing has been performed to discover the novel and conserved miRNAs in mature cabbage leaves.

Results

In this study, raw reads generated from three small RNA libraries were bioinformatically processed and further analyzed to select sequences homologous to known B. oleracea and other plant miRNAs. As a result of this analysis, 261 conserved miRNAs (belonging to 62 families) have been discovered. MIR169, MIR167 and MIR166 were the largest miRNA families, while the highest abundance molecules were miR167, miR166, miR168c and miR157a. Among the generated sequencing reads, miRNAs* were also found, such as the miR162c*, miR160a* and miR157a*. The unannotated tags were used in the prediction and evaluation of novel miRNAs, which resulted in the 26 potential miRNAs proposal. The expressions of 13 selected miRNAs were analyzed by northern blot hybridization. The target prediction and annotation for identified miRNAs were performed, according to which discovered molecules may target mRNAs encoding several potential proteins – e.g., transcription factors, polypeptides that regulate hormone stimuli and abiotic stress response, and molecules participating in transport and cell communication. Additionally, KEGG maps analysis suggested that the miRNAs in cabbage are involved in important processing pathways, including glycolysis, glycerolipid metabolism, flavonoid biosynthesis and oxidative phosphorylation.

Conclusions

Conclusively, for the first time, the large set of miRNAs was identified in mature cabbage leaves. Potential targets designation for these miRNAs may suggest their essential role in many plants primary biological processes. Presented study not only supplements the knowledge about B. oleracea miRNAs, but additionally it may be used in other research concerning the improvement of the cabbage cultivation.

Identification of potential miRNAs and their targets in Vriesea carinata (Poales, Bromeliaceae)

The miRNAs play important roles in regulation of gene expression at the post-transcriptional level. A small RNA and RNA-seq of libraries were constructed to identify miRNAs in Vriesea carinata, a native bromeliad species from Brazilian Atlantic Rainforest. Illumina technology was used to perform high throughput sequencing and data was analyzed using bioinformatics tools. We obtained 2,191,509 mature miRNAs sequences representing 54 conserved families in plant species. Further analysis allowed the prediction of secondary structures for 19 conserved and 16 novel miRNAs. Potential targets were predicted from pre-miRNAs by sequence homology and validated using RTqPCR approach. This study provides the first identification of miRNAs and their potential targets of a bromeliad species.