Summary: Naturally coloured cotton (NCC) fibres need little or no dyeing process in textile industry to low‐carbon emission and are environment‐friendly. Proanthocyanidins (PAs) and their derivatives were considered as the main components causing fibre coloration and made NCCs very popular and healthy, but the monotonous fibre colours greatly limit the wide application of NCCs. Here a G. hirsutum empurpled mutant (HS2) caused by T‐DNA insertion is found to enhance the anthocyanidins biosynthesis and accumulate anthocyanidins in the whole plant. HPLC and LC/MS‐ESI analysis confirmed the anthocyanidins methylation and peonidin, petunidin and malvidin formation are blocked. The deficiency of GhOMT1 in HS2 was associated with the activation of the anthocyanidin biosynthesis and the altered components of anthocyanidins. The transcripts of key genes in anthocyanidin biosynthesis pathway are significantly up‐regulated in HS2, while transcripts of the genes for transport and decoration were at similar levels as in WT. To investigate the potential mechanism of GhOMT1 deficiency in cotton fibre coloration, HS2 mutant was crossed with NCCs. Surprisingly, offsprings of HS2 and NCCs enhanced PAs biosynthesis and increased PAs levels in their fibres from the accumulated anthocyanidins through up‐regulated GhANR and GhLAR. As expected, multiple novel lines with improved fibre colours including orange red and navy blue were produced in their generations. Based on this work, a new strategy for breeding diversified NCCs was brought out by promoting PA biosynthesis. This work will help shed light on mechanisms of PA biosynthesis and bring out potential molecular breeding strategy to increase PA levels in NCCs.

Cotton Verticillium wilt, caused by the notorious fungal phytopathogen Verticillium dahliae (V. dahliae), is a destructive soil-borne vascular disease and severely decreases cotton yield and quality worldwide. Transcriptional and post-transcriptional regulation of genes responsive to V. dahliae are crucial for V. dahliae tolerance in plants. However, the specific microRNAs (miRNAs) and the miRNA/target gene crosstalk involved in cotton resistance to Verticillium wilt remain largely limited. To investigate the roles of regulatory RNAs under V. dahliae induction in upland cotton, mRNA and small RNA libraries were constructed from mocked and infected roots of two upland cotton cultivars with the V. dahliae-sensitive cultivar Jimian 11 (J11) and the V. dahliae-tolerant cultivar Zhongzhimian 2 (Z2). A comparative transcriptome analysis revealed 8330 transcripts were differentially expressed under V. dahliae stress and associated with several specific biological processes. Moreover, small RNA sequencing identified a total of 383 miRNAs, including 330 unique conserved miRNAs and 53 novel miRNAs. Analysis of the regulatory network involved in the response to V. dahliae stress revealed 31 differentially expressed miRNA–mRNA pairs, and the up-regulation of GhmiR395 and down-regulation of GhmiR165 were possibly involved in the response to V. dahliae by regulating sulfur assimilation through the GhmiR395-APS1/3 module and the establishment of the vascular pattern and secondary cell wall formation through GhmiR165-REV module, respectively. The integrative analysis of mRNA and miRNA expression profiles from upland cotton lays the foundation for further investigation of regulatory mechanisms of resistance to Verticillium wilt in cotton and other crops.

Background: Cold stress is one of the primary environmental factors that affect plant growth and productivity, especially for crops like Brassica napus that live through cold seasons. Till recently, although a number of genes and pathways involved in B. napus cold response have been revealed by independent studies, a genome-wide identification of the key regulators and the regulatory networks is still lack. In this study, we investigated the transcriptomes of cold stressed semi-winter and winter type rapeseeds in short day condition, mainly with the purpose to systematically identify the functional conserved transcription factors (TFs) in cold response of B. napus. Results: Global modulation of gene expression was observed in both the semi-winter type line (158A) and the winter type line (SGDH284) rapeseeds, in response to a seven-day chilling stress in short-day condition. Function analysis of differentially expressed genes (DEGs) revealed enhanced stresses response mechanisms and inhibited photosynthesis in both lines, as well as a more extensive inhibition of some primary biological processes in the semi-winter type line. Over 400 TFs were differentially expressed in response to cold stress, including 56 of them showed high similarity to the known cold response TFs and were consistently regulated in 158A and SGDH284, as well as 25 TFs which targets were over-represented in the total DEGs. A further investigation based on their interactions indicated the critical roles of several TFs in cold response of B. napus. Conclusion: In summary, our results revealed the alteration of gene expression in cold stressed semi-winter and winter ecotype B. napus lines and provided a valuable collection of candidate key regulators involved in B. napus response to cold stress, which could expand our understanding of plant stress response and benefit the future improvement of the breed of rapeseeds.

Arsenic (As) contamination in subsoil and groundwater is a big problem, especially in many South-East Asian countries. As a staple crop growing under flooded condition in these areas, rice (Oryza sativa L.) becomes a big threat to human health through the food chain since As is highly accumulated in grains. Thus, reducing As accumulation in rice through molecular breeding and identification of rice varieties with low As content are the pressing issues. However, the current understanding on the molecular mechanism of As stress response is still limited for rice. In this study, we performed a comprehensive search for the As-responsive small RNAs (sRNAs) of rice. Briefly, 4,762 and 18,152 sRNAs were identified to be highly activated under As stress in roots and shoots respectively, while 14,603 and 8,308 sRNAs were intensively repressed by As treatment in roots and shoots, respectively. A number of the As-responsive sRNAs found their loci on tRNAs, rRNAs or long non-coding RNAs (lncRNAs). Interestingly, these loci preferentially distributed on the 5ʹ halves of the tRNA, rRNA or lncRNA precursors. Among the above-identified As-responsive sRNAs, 252 Argonaute 1 (AGO1)-enriched sRNAs were extracted for target identification, resulting in 200 pairs of sRNA–protein-coding target interactions. Many targets are functionally involved in the development, stress response, reproduction, or lipid metabolism. Additionally, 56 lncRNAs were discovered to be targeted by nine AGO1-enriched sRNAs, indicating the potential involvement of these lncRNAs in As signaling. Taken together, our results could expand the understanding on the non-coding RNA-mediated As stress response in rice.

A great mass of long noncoding RNAs (lncRNAs) have been identified in mouse genome and increasing evidences in the last decades have revealed their crucial roles in diverse biological processes. Nevertheless, the biological roles of lncRNAs in the mouse retina remains largely unknown due to the lack of a comprehensive annotation of lncRNAs expressed in the retina. In this study, we applied the long-reads sequencing strategy to unravel the transcriptomes of developing mouse retinas and identified a total of 940 intergenic lncRNAs (lincRNAs) in embryonic and neonatal retinas, including about 13% of them were transcribed from unannotated gene loci. Subsequent analysis revealed that function of lincRNAs expressed in mouse retinas were closely related to the physiological roles of this tissue, including 90 lincRNAs that were differentially expressed after the functional loss of key regulators of retinal ganglion cell (RGC) differentiation. In situ hybridization results demonstrated the enrichment of three class IV POU-homeobox genes adjacent lincRNAs (linc-3a, linc-3b and linc-3c) in ganglion cell layer and indicated they were potentially RGC-specific. In summary, this study systematically annotated the lincRNAs expressed in embryonic and neonatal mouse retinas and implied their crucial regulatory roles in retinal development such as RGC differentiation.

The wild cotton species Gossypium stocksii produces a brown fiber that provides a valuable resource for the color improvement of naturally colored cotton (NCC) fiber. However, the biochemical basis and molecular mechanism of its fiber pigmentation remain unclear. Herein, we analyzed the dynamics of proanthocyanidins (PAs) accumulation in developing the fiber of G. stocksii, which suggested a similar role of PAs and/or their derivatives in the fiber coloration of G. stocksii. In addition, comparative transcriptomics analyses revealed that the PA biosynthetic genes were expressed at higher levels and for a longer period in developing fibers of G. stocksii than G. arboreum (white fiber), and the transcription factors, such as TT8, possibly played crucial regulatory roles in regulating the PA branch genes. Moreover, we found that the anthocyanidin reductase (ANR) was expressed at a higher level than the leucoanthocyanidin reductases (LARs) and significantly upregulated during fiber elongation, suggesting a major role of ANR in PA synthesis in G. stocksii fiber. In summary, this work revealed the accumulation of PAs and the expression enhancement of PA biosynthetic genes in developing fibers of G. stocksii. We believe this work will help our understanding of the molecular mechanisms of cotton fiber coloration and further promote the future breeding of novel NCCs.