MicroRNA and mRNA expression profiling analysis revealed the regulation of plant height in Gossypium hirsutum

Comparative Microscopic, Transcriptome and IAA Content Analyses Reveal the Stem Growth Variations in Two Cultivars Ilex verticillata

Ilex verticillata is not only an excellent ornamental tree species for courtyards, but it is also a popular bonsai tree. 'Oosterwijk' and 'Red sprite' are two varieties of Ilex verticillata. The former has a long stem with few branches, while the latter has a short stem. In order to explain the stem growth differences between the two cultivars 'Oosterwijk' and 'Red sprite', determination of the microstructure, transcriptome sequence and IAA content was carried out. The results showed that the xylem thickness, vessel area and vessel number of 'Oosterwijk' were larger than in 'Red sprite'. In addition, our analysis revealed that the differentially expressed genes which were enriched in phenylpropanoid biosynthesis; phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis in the black and tan modules of the two varieties. We found that AST, HCT and bHLH 94 may be key genes in the formation of shoot difference. Moreover, we found that the IAA content and auxin-related DEGs GH3.6, GH3, ATRP5, IAA27, SAUR36-like, GH3.6-like and AIP 10A5-like may play important roles in the formation of shoot differences. In summary, these results indicated that stem growth variations of 'Oosterwijk' and 'Red sprite' were associated with DEGs related to phenylpropanoid biosynthesis, phenylalanine metabolism and phenylalanine, tyrosine and tryptophan biosynthesis, as well as auxin content and DEGs related to the auxin signaling pathway.

Regulatory roles of miRNAs associated with the aging pathway in tree vegetative phase changes

The transition from the vegetative juvenile phase to the adult phase is a crucial event in the life cycle of flowering plants, with flowering being the most important milestone. While the regulatory pathways of flowering have been well established in model plants such as Arabidopsis and a few crops, the flowering regulation pathways in perennial forest trees remain poorly understood. This paper summarizes the regulation of flowering time by miR156 and miR172, which are the main members of the aging pathway, and also presents new information on the role of miR159 and miR169. These two microRNAs interact with miR156 and miR172 to jointly regulate flowering time in forest trees. Overall, this review sheds light on the complex regulatory mechanisms underlying flowering time in forest trees and provides insights into potential targets for manipulating the flowering time of these economically and ecologically important species.

Integration of miRNA and mRNA analysis reveals the role of ribosome in to anti-artificial aging in sweetcorn

Sweetcorn is a kind of maize with high sugar content and has poor seed aging tolerance, which seriously limits its production. However, few studies have explored the artificial aging (AA) tolerance by miRNA-mRNA integration analysis in sweetcorn. Here, we characterized the physiological, biochemical and transcriptomic changes of two contrasting lines K62 and K107 treated with AA during time series. Both the germination indexes and antioxidant enzymes showed significant difference between two lines. The MDA content of AA-tolerant genotype K62 was significantly lower than that of K107 on the fourth and sixth day. Subsequently, 157 differentially expressed miRNAs (DEMIs) and 8878 differentially expressed mRNAs (DEMs) were identified by RNA-seq analysis under aging stress. The "ribosome" and "peroxisome" pathways were enriched to respond to aging stress, genes for both large units and small ribosomal subunits were significantly upregulated expressed and higher translation efficiency might exist in K62. Thirteen pairs of miRNA-target genes were obtained, and 8 miRNA-mRNA pairs might involve in ribosome protein and translation process. Our results elucidate the mechanism of sweetcorn response to AA at miRNA-mRNA level, and provide a new insight into sweetcorn AA response to stress.

Mapping of dwarfing QTL of Ari1327, a semi-dwarf mutant of upland cotton

BACKGROUND

Upland Cotton (Gossypium hirsutum L.) has few cotton varieties suitable for mechanical harvesting. The plant height of the cultivar is one of the key features that need to modify. Hence, this study was planned to locate the QTL for plant height in a 60Co γ treated upland cotton semi-dwarf mutant Ari1327.

RESULTS

Interestingly, bulk segregant analysis (BSA) and genotyping by sequencing (GBS) methods exhibited that candidate QTL was co-located in the region of 5.80-9.66 Mb at D01 chromosome in two F2 populations. Using three InDel markers to genotype a population of 1241 individuals confirmed that the offspring's phenotype is consistent with the genotype. Comparative analysis of RNA-seq between the mutant and wild variety exhibited that Gh_D01G0592 was identified as the source of dwarfness from 200 genes. In addition, it was also revealed that the appropriate use of partial separation markers in QTL mapping can escalate linkage information.

CONCLUSIONS

Overwhelmingly, the results will provide the basis to reveal the function of candidate genes and the utilization of excellent dwarf genetic resources in the future.

Profiling of MicroRNAs Involved in Mepiquat Chloride-Mediated Inhibition of Internode Elongation in Cotton (Gossypium hirsutum L.) Seedlings

Mepiquat chloride (MC) is the most important plant growth retardant that is widely used in cotton (Gossypium hirsutum L.) production to suppress excessive vegetative growth and improve plant architecture. MicroRNAs (miRNAs) are important gene expression regulators that control plant growth and development. However, miRNA-mediated post-transcriptional regulation in MC-induced growth inhibition remains unclear. In this study, the dynamic expression profiles of miRNAs responsive to MC in cotton internodes were investigated. A total of 508 known miRNAs belonging to 197 families and five novel miRNAs were identified. Among them, 104 miRNAs were differentially expressed at 48, 72, or 96 h post MC treatment compared with the control (0 h); majority of them were highly conserved miRNAs. The number of differentially expressed miRNAs increased with time after treatment. The expression of 14 known miRNAs was continuously suppressed, whereas 12 known miRNAs and one novel miRNA were continuously induced by MC. The expression patterns of the nine differentially expressed miRNAs were verified using qRT-PCR. The targets of the known and novel miRNAs were predicted. Four conserved and six novel targets were validated using the RLM-5' RACE assay. This study revealed that miRNAs play crucial regulatory roles in the MC-induced inhibition of internode elongation. It can improve our understanding of post-transcriptional gene regulation in MC-mediated growth inhibition and could potentially facilitate the breeding of dwarf cotton.

Identification of miRNA-eQTLs in maize mature leaf by GWAS

Background

MiRNAs play essential roles in plant development and response to biotic and abiotic stresses through interaction with their target genes. The expression level of miRNAs shows great variations among different plant accessions, developmental stages, and tissues. Little is known about the content within the plant genome contributing to the variations in plants. This study aims to identify miRNA expression-related quantitative trait loci (miR-QTLs) in the maize genome.

Results

The miRNA expression level from next generation sequencing (NGS) small RNA libraries derived from mature leaf samples of the maize panel (200 maize lines) was estimated as phenotypes, and maize Hapmap v3.2.1 was chosen as the genotype for the genome-wide association study (GWAS). A total of four significant miR-eQTLs were identified contributing to miR156k-5p, miR159a-3p, miR390a-5p and miR396e-5p, and all of them are trans-eQTLs. In addition, a strong positive coexpression of miRNA was found among five miRNA families. Investigation of the effects of these miRNAs on the expression levels and target genes provided evidence that miRNAs control the expression of their targets by suppression and enhancement.

Conclusions

These identified significant miR-eQTLs contribute to the diversity of miRNA expression in the maize penal at the developmental stages of mature leaves in maize, and the positive and negative regulation between miRNA and its target genes has also been uncovered.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Identification of Grafting-Responsive MicroRNAs Associated with Growth Regulation in Pecan [Carya illinoinensis (Wangenh.) K. Koch]

Pecan [Carya illinoinensis (Wangenh.) K. Koch] is an economically important nut tree and grafting is often used for clonal propagation of cultivars. However, there is a lack of research on the effects of rootstocks on scions, which are meaningful targets for directed breeding of pecan grafts. MicroRNAs (miRNAs) play an important role in many biological processes, but the mechanism underlying the involvement of miRNAs in grafting-conferred physiological changes is unclear. To identify the grafting-responsive miRNAs that may be involved in the regulation of growth in grafted pecan, six small RNA libraries were constructed from the phloem of two groups of grafts with significantly different growth performance on short and tall rootstocks. A total of 441 conserved miRNAs belonging to 42 miRNA families and 603 novel miRNAs were identified. Among the identified miRNAs, 24 (seven conserved and 17 novel) were significantly differentially expressed by the different grafts, implying that they might be responsive to grafting and potentially involved in the regulation of graft growth. Ninety-five target genes were predicted for the differentially expressed miRNAs; gene annotation was available for 33 of these. Analysis of their targets suggested that the miRNAs may regulate auxin transport, cell activity, and inorganic phosphate (Pi) acquisition, and thereby, mediate pecan graft growth. Use of the recently-published pecan genome enabled identification of a substantial population of miRNAs, which are now available for further research. We also identified the grafting-responsive miRNAs and their potential roles in pecan graft growth, providing a basis for research on long-distance regulation in grafted pecan.

Transcriptome sequencing reveals hotspot mutation regions and dwarfing mechanisms in wheat mutants induced by γ-ray irradiation and EMS

Induced mutation is an important approach for creating novel plant germplasms. The introduction of dwarf or semi-dwarf genes into wheat has led to great advancements in yield improvement. In this study, four elite dwarf wheat mutants, named dm1-dm4, induced from γ-ray irradiation or ethyl methanesulfonate (EMS) mutagenesis, were used to identify transcriptome variations and dwarfing mechanisms. The results showed that the hotspot regions of mutations distributed on the chromosomes were consistent among the four mutant lines and these regions were mainly located around the 50, 360 and 400 Mb positions of chromosome 1A and the distal regions of chromosomes 2A and 2BL. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that 'protein processing in endoplasmic reticulum' was the most common significantly enriched pathway based on the differentially expressed genes (DEGs) between wildtype (WT) and the mutants. Notably, 18 out of 20 genes involved in this process encode heat shock proteins (HSPs). The results implied that HSPs might participate in wheat dwarfism response and function in the dwarfism process through protein folding and/or degradation. Moreover, seven genes in dm4 involved in modulating auxin levels were down-regulated and dm4 was more sensitive to auxin treatment compared with WT, indicating the important roles of auxin in regulation of dwarf phenotype in dm4. This study not only identified transcriptome sequence variation induced by physical and chemical mutagenesis but also revealed potential dwarfing mechanisms in the wheat mutant lines.

Long non-coding RNAs and their potential functions in Ligon-lintless-1 mutant cotton during fiber development

BACKGROUND

Long non-coding RNAs (LncRNAs) are part of genes, which are not translated into proteins and play a vital role in plant growth and development. Nevertheless, the presence of LncRNAs and how they functions in Ligon-lintless-1 mutant during the early cessation of cotton fiber development are still not well understood. In order to investigate the function of LncRNAs in cotton fiber development, it is necessary and important to identify LncRNAs and their potential roles in fiber cell development.

RESULTS

In this work, we identified 18,333 LncRNAs, with the proportion of long intergenic noncoding RNAs (LincRNAs) (91.5%) and anti-sense LncRNAs (8.5%), all transcribed from Ligon-lintless-1 (Li1) and wild-type (WT). Expression differences were detected between Ligon-lintless-1 and wild-type at 0 and 8 DPA (day post anthesis). Pathway analysis and Gene Ontology based on differentially expressed LncRNAs on target genes, indicated fatty acid biosynthesis and fatty acid elongation being integral to lack of fiber in mutant cotton. The result of RNA-seq and RT-qPCR clearly singles out two potential LncRNAs, LNC_001237 and LNC_017085, to be highly down-regulated in the mutant cotton. The two LncRNAs were found to be destabilized or repressed by ghr-miR2950. Both RNA-seq analysis and RT-qPCR results in Ligon-lintless-1 mutant and wild-type may provide strong evidence of LNC_001237, LNC_017085 and ghr-miR2950 being integral molecular elements participating in various pathways of cotton fiber development.

CONCLUSION

The results of this study provide fundamental evidence for the better understanding of LncRNAs regulatory role in the molecular pathways governing cotton fiber development. Further research on designing and transforming LncRNAs will help not only in the understanding of their functions but will also in the improvement of fiber quality.

Joint RNA-Seq and miRNA Profiling Analyses to Reveal Molecular Mechanisms in Regulating Thickness of Pod Canopy in Brassica napus

Oilseed rape (Brassica napus) is the second largest oilseed crop worldwide. As an architecture component of B. napus, thickness of pod canopy (TPC) plays an important role in yield formation, especially under high-density cultivation conditions. However, the mechanisms underlying the regulation of TPC remain unclear. RNA and microRNA (miRNA) profiling of two groups of B. napus lines with significantly different TPC at the bolting with a tiny bud stage revealed differential expressions of numerous genes involved in nitrogen-related pathways. Expression of several nitrogen-related response genes, including ASP5, ASP2, ASN3, ATCYSC1, PAL2, APT2, CRTISO, and COX15, was dramatically changed in the thick TPC lines compared to those in the thin TPC lines. Differentially expressed miRNAs also included many involved in nitrogen-related pathways. Expression of most target genes was negatively associated with corresponding miRNAs, such as miR159, miR6029, and miR827. In addition, 12 (including miR319, miR845, and miR158) differentially expressed miRNAs between two plant tissues sampled (stem apex and flower bud) were identified, implying that they might have roles in determining overall plant architecture. These results suggest that nitrogen signaling may play a pivotal role in regulating TPC in B. napus.

Identification and profiling of upland cotton microRNAs at fiber initiation stage under exogenous IAA application

Background

Cotton is the most essential textile crop worldwide, and phytohormones are critical for cotton fiber development. One example is the role of auxin in fiber initiation, but we know little molecular basis. MicroRNAs (miRNAs) have a significant function in cotton development; nevertheless their role in fiber initiation remains unclear. Here, exogenous IAA was applied to cotton plant before anthesis. Utilizing small RNA sequencing, the mechanism underlying miRNA-mediated regulation of fiber initiation under exogenous IAA treatment was investigated.

Results

With exogenous IAA application, the endogenous IAA and GA contents of IAA treated (IT) ovules were higher than control (CK) ovules at the fiber initiation stage, while endogenous ABA content was lower in IT than CK. Using scanning electron microscopy, we found the fiber number and size were significantly promoted in IT at 0 DPA. Fiber quality analysis showed that fiber length, uniformity, strength, elongation, and micronaire of IT were higher than CK, though not statistically significant, while lint percent was significantly higher in IT. We generated six small RNA libraries using − 3, 0, and 3 DPA ovules of IT and CK, and identified 58 known miRNAs and 83 novel miRNAs together with the target genes. The differential expressed miRNAs number between IT and CK at − 3, 0, 3 DPA was 34, 16 and 24, respectively. Gene ontology and KEGG pathway enrichment analyses for the target genes of the miRNAs expressed in a differential manner showed that they were significantly enriched in 30 terms and 8 pathways. QRT-PCR for those identified miRNAs and the target genes related to phytohormones and fiber development was performed, and results suggested a potential role of these miRNAs in fiber initiation.

Conclusions

The exogenous IAA application affected the relative phytohormone contents in ovule and promoted fiber initiation in cotton. Identification and profiling of miRNAs and their targets at the fiber initiation stage provided insights for miRNAs’ regulation function of fiber initiation. These findings not only shed light on the regulatory network of fiber growth but also offer clues for cotton fiber amelioration strategies in cotton.

Electronic supplementary material

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

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

Global transcriptome and weighted gene co-expression network analyses reveal hybrid-specific modules and candidate genes related to plant height development in maize

Weighted gene co-expression network analysis was explored to find key hub genes involved in plant height regulation. Plant height, an important trait for maize breeding because of its close relatedness to lodging resistance and yield, has been reported to be determined by multiple qualitative and quantitative genes. However, few genes related to plant height have been characterized in maize. Herein, three different maize hybrids, with extremely distinct plant height, which were further classified into low (L), middle (M) and high (H) group, were selected for RNA sequencing at three key developmental stages, namely, jointing stage (I), big flare period (II) and tasseling stage (III). Intriguingly, transcriptome profiles for hybrids ranging from low to high group exhibited significantly similarity in both jointing stage and big flare period. However, remarkably larger differentially expressed genes could be detected between hybrid from low to either middle or high group in tasseling stage. These results were repeatedly observed in both phenotyping and gene ontology enrichment analysis, indicating that transition from big flare period to tasseling stage plays a critical role in determination of plant height. Furthermore, weighted gene co-expression network analysis was explored to find key hub genes involved in plant height regulation. Hundreds of candidate genes, encoding various transcription factors, and regulators involved in internode cell regulation and cell wall synthesis were identified in our network. More importantly, great majority of candidates were correlated to either metabolism or signaling pathway of several plant phytohormones. Particularly, numerous functionally characterized genes in gibberellic acid as well as brassinosteroids signaling transduction pathways were also discovered, suggesting their critical roles in plant height regulation. The present study could provide a modestly comprehensive insight into networks for regulation of plant height in maize.

Genome-wide analysis of poplar SAUR gene family and expression profiles under cold, polyethylene glycol and indole-3-acetic acid treatments

Small auxin-up RNA (SAUR) proteins play an important role in the regulation of plant growth and development. Here, we identified 105 SAUR genes and comprehensively analyzed them in Populus trichocarpa. Based on the phylogenetic relationships, the PtSAURs were classified into ten subfamilies. Of the 105 PtSAURs, 100 were randomly distributed along the nineteen chromosomes, while the remaining genes were located along unassigned scafoolds. These genes mainly evolved through segmental duplications. In total, 94 PtSAURs contained no introns, and each group had a similar conserved motif structure. A promoter analysis revealed various cis-elements related to growth, development and stress responses, and a synteny analysis established orthologous relationships among SAURs in Arabidopsis, rice, grape and poplar. The qRT-PCR and tissue expression analyses indicated that PtSAURs show different expression levels in various tissues in response to different treatments. PtSAUR53 was located on the nuclear and plasma membrane by conducting subcellular localization analysis. This study provides a comprehensive overview of poplar SAUR proteins and a foundation for further investigations for functional analysis of SAURs in poplar growth and development. At the same time, it will be valuable to further study the poplar SAUR genes to reveal their biological effects.

A genome-wide analysis of the small auxin-up RNA (SAUR) gene family in cotton

Background

Small auxin-up RNA (SAUR) gene family is the largest family of early auxin response genes in higher plants, which have been implicated in the regulation of multiple biological processes. However, no comprehensive analysis of SAUR genes has been reported in cotton (Gossypium spp.).

Results

In the study, we identified 145, 97, 214, and 176 SAUR homologous genes in the sequenced genomes of G. raimondii, G. arboreum, G. hirsutum, and G. barbadense, respectively. A phylogenetic analysis revealed that the SAUR genes can be classified into 10 groups. A further analysis of chromosomal locations and gene duplications showed that tandem duplication and segmental duplication events contributed to the expansion of the SAUR gene family in cotton. An exon-intron organization and motif analysis revealed the conservation of SAUR-specific domains, and the auxin responsive elements existed in most of the upstream sequences. The expression levels of 16 GhSAUR genes in response to an exogenous application of IAA were determined by a quantitative RT-PCR analysis. The genome-wide RNA-seq data and qRT-PCR analysis of selected SAUR genes in developing fibers revealed their differential expressions. The physical mapping showed that 20 SAUR genes were co-localized with fiber length quantitative trait locus (QTL) hotspots. Single nucleotide polymorphisms (SNPs) were detected for 12 of these 20 genes between G. hirsutum and G. barbadense, but no SNPs were identified between two backcross inbred lines with differing fiber lengths derived from a cross between the two cultivated tetraploids.

Conclusions

This study provides an important piece of genomic information for the SAUR genes in cotton and lays a solid foundation for elucidating the functions of SAUR genes in auxin signaling pathways to regulate cotton growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4224-2) contains supplementary material, which is available to authorized users.

Integrated analysis of mRNA and miRNA expression profiling in rice backcrossed progenies (BC2F12) with different plant height

Inter-specific hybridization and backcrossing commonly occur in plants. The use of progeny generated from inter-specific hybridization and backcrossing has been developed as a novel model system to explore gene expression divergence. The present study investigated the analysis of gene expression and miRNA regulation in backcrossed introgression lines constructed from cultivated and wild rice. High-throughput sequencing was used to compare gene and miRNA expression profiles in three progeny lines (L1710, L1817 and L1730), with different plant heights resulting from the backcrossing of introgression lines (BC₂F₁₂) and their parents (O. sativa and O. longistaminata). A total of 25,387 to 26,139 mRNAs and 379 to 419 miRNAs were obtained in these rice lines. More differentially expressed genes and miRNAs were detected in progeny/O. longistaminata comparison groups than in progeny/O. sativa comparison groups. Approximately 80% of the genes and miRNAs showed expression level dominance to O. sativa, indicating that three progeny lines were closer to the recurrent parent, which might be influenced by their parental genome dosage. Approximately 16% to 64% of the differentially expressed miRNAs possessing coherent target genes were predicted, and many of these miRNAs regulated multiple target genes. Most genes were up-regulated in progeny lines compared with their parents, but down-regulated in the higher plant height line in the comparison groups among the three progeny lines. Moreover, certain genes related to cell walls and plant hormones might play crucial roles in the plant height variations of the three progeny lines. Taken together, these results provided valuable information on the molecular mechanisms of hybrid backcrossing and plant height variations based on the gene and miRNA expression levels in the three progeny lines.

Suppressive effect of microRNA319 expression on rice plant height

KEY MESSAGE: miR319 was identified as a dwarf-inducing gene from Shiokari and its dwarf near isogenic line, and its transgenic rice showed a reduced plant height. This finding reveals the potential application of miR319 in future molecular breeding. It is well known that microRNAs (miRNAs) play important roles in plant physiology, especially in development and stress responses. However, little is known about the role of miRNAs in plant height. In this study, the rice cultivar Shiokari and its dwarf near isogenic line Shiokari-d6 were analysed to identify and characterize plant height-associated miRNAs. This anatomic and morphological investigation revealed that the major cause of the shorter height of Shiokari-d6 is the significantly dis-elongated internodes, particularly the second internode and those underneath it. The results of miRNA microarray profiling and real-time RT-PCR indicated that miR319 is expressed at a significantly higher level in Shiokari-d6 than in Shiokari. Transgenic rice overexpressing miR319 in Oryza sativa L. cv. Tainung 67 generated through Agrobacterium-mediated transformation had a stable dwarf phenotype regardless of whether the plants were from the T1 or T2 generation. We also found that the internodes of miR319-overexpressing rice are shortened, particularly the third internode and those underneath it. Furthermore, we identified three putative miR319 target genes that were previously uncharacterized with expression levels that were negatively correlated with the expression of miR319. In conclusion, miR319 is the first miRNA proposed to be involved in plant height regulation, and its function may influence the elongation of internodes, which leads to decreased plant height.

Genome-wide identification and characterization of microRNAs differentially expressed in fibers in a cotton phytochrome A1 RNAi line

Cotton fiber is an important commodity throughout the world. Fiber property determines fiber quality and commercial values. Previous studies showed that silencing phytochrome A1 gene (PHYA1) by RNA interference in Upland cotton (Gossypium hirsutum L. cv. Coker 312) had generated PHYA1 RNAi lines with simultaneous improvements in fiber quality (longer, stronger and finer fiber) and other key agronomic traits. Characterization of the altered molecular processes in these RNAi genotypes and its wild-type controls is a great interest to better understand the PHYA1 RNAi phenotypes. In this study, a total of 77 conserved miRNAs belonging to 61 families were examined in a PHYA1 RNAi line and its parental Coker 312 genotype by using multiplex sequencing. Of these miRNAs, seven (miR7503, miR7514, miR399c, miR399d, miR160, miR169b, and miR2950) were found to be differentially expressed in PHYA1 RNAi cotton. The target genes of these differentially expressed miRNAs were involved in the metabolism and signaling pathways of phytohormones, which included Gibberellin, Auxin and Abscisic Acid. The expression of several MYB transcription factors was also affected by miRNAs in RNAi cotton. In addition, 35 novel miRNAs (novel miR1-novel miR35) were identified in fibers for the first time in this study. Target genes of vast majority of these novel miRNAs were also predicted. Of these, nine novel miRNAs (novel-miR1, 2, 16, 19, 26, 27, 28, 31 and 32) were targeted to cytochrome P450-like TATA box binding protein (TBP). The qRT-PCR confirmed expression levels of several differentially regulated miRNAs. Expression patterns of four miRNAs-targets pairs were also examined via RNA deep sequencing. Together, the results imply that the regulation of miRNA expression might confer to the phenotype of the PHYA1 RNAi line(s) with improved fiber quality.