Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

Rate of double reduction and genetic variability in yield, quality, and senescence related traits in tetraploid potato (Solanum tuberosum L.)

The amount of genetic variability is the foundation for genetic change in any plant breeding program, and the amount of double reduction can influence genetic gain and the amount of future genetic diversity in polyploid species. Our study investigates these factors using variance components analysis on a dataset comprising 13,131 potato breeding lines and phenotypic data from Scandinavian environments spanning 17 years (2003 to 2021). Pedigree information was used in quantitative genetic models to estimate additive genetic variance and the relative importance of additive and non-additive genetic variance. We used two models, a baseline model (M1) without effects due to specific combining ability (SCA) and M2 (including SCA due to interaction between parental genomes). Two cross-validation (CV) schemes [5-Fold and leave-one-breeding-cycle-out (LBCO)] were used to evaluate the prediction ability (PA) of each model. We estimated the rate of double reduction phenomenon (DRP) by determining the rate best fitting the data using a marginal likelihood approach. Our findings showed a wide range of variation in different traits, with very large proportion of additive genetic variance in dry matter content (DMC), but intermediate additive genetic variance for relative yield (RY), germination (GR), and withering (WNG). All traits showed modest non-additive genetic variance. Furthermore, genotype x environment interaction played a significant role in trait variability but is still much smaller than the additive genetic variance. After using different DRP rates, we found that a model with a 0.05 DRP rate provided the best fit to the data. Heritability estimates indicated a strong genetic basis for DMC, while other traits showed more moderate heritability, which shows contributions from both additive and interaction factors. Model comparison by 5-Fold CV and LBCO and the log likelihood ratio test (LRT) highlighted the importance of considering SCA when capturing trait variability. In 5-Fold CV, PA ranged from 0.296 to 0.812 in M1 and 0.300 to 0.813 in M2. Under LBCO CV, PA ranged from 0.180 to 0.726 in M1 and 0.180 to 0.728 in M2. However, an increase in PA in Model 2, which incorporates SCA, compared to Model 1, can be attributed to the inclusion of SCA effects. Furthermore, the LRT results indicated a highly significant difference between the models. CV and LRT suggest the need for genetic models that account for both additive and SCA effects. Our analysis also showed that genotype x environment interactions should be accounted for in order to maximize the accuracy of predicted breeding values of tetraploid potato clones. The rate of double reductions was small and insignificant.

Global identification and regulatory network analysis reveal the significant roles of lncRNAs during anther and pollen development in Arabidopsis

KEY MESSAGE

A high-throughput sequencing identified 1283 lncRNAs in anthers at different stages in Arabidopsis and their relationship with protein-coding genes and miRNAs during anther and pollen development were analyzed. Long non-coding RNAs (lncRNAs) are important regulatory molecules involved in various biological processes. However, their roles in male reproductive development and interactions with miRNAs remained elusive. In this study, a high-throughput sequencing of anthers at different developmental stages in Arabidopsis identified 1283 lncRNAs including 524 differentially expressed lncRNAs (DELs). Most of these DELs exhibited positive correlations with the expression patterns of adjacent protein-coding genes. Weighted gene co-expression network analysis (WGCNA) revealed that protein-coding genes targeted by DELs in four modules related to the tetrad stage were associated with functions such as pollen wall formation, pollen germination, or pollen tube growth, respectively. Furthermore, five, 10, and 11 lncRNAs were predicted as miRNAs' endogenous target mimics (eTMs), precursors, and natural antisense transcripts of pri-miRNA, respectively. Remarkably, the lncRNA, host gene of ath-miR167a (ath-miR167aHG), predicted as the precursor of miR167a, was selected for function validation. Its overexpression resulted in the up-regulation of miR167a and the subsequent down-regulation of miR167a's target genes ARF6 and ARF8, demonstrating a functional interaction between ath-miR167aHG and miR167a. The transgenic plants showed delayed flowering, shorter filaments, abnormal anther dehiscence, and undeveloped siliques ultimately, suggesting a role of ath-miR167aHG in male reproductive development. Collectively, our research shed new light on the functions of lncRNAs in male reproductive development and uncovered the unique interactions between lncRNAs and miRNAs.

RNA sequencing and functional analysis uncover key long non-coding RNAs involved in regulating pollen fertility during the process of gametocidal action in wheat

Gametocidal (Gc) chromosomes have been widely utilized in genetic breeding due to their ability to induce chromosomal breakage and eliminate gametes that lack them. Long noncoding RNAs (lncRNAs) have various functional mechanisms in regulating pollen and anther development; however, their regulatory contributions to Gc action are still unknown. Here, we identified 2824 differentially expressed lncRNAs (DE-lncRNAs) from the anther tissues of Triticum aestivum cv. Chinese Spring (CS) and Chinese Spring-Gc 3C chromosome monosomic addition line (CS-3C) through sequencing. In this study, we predicted 161 target mRNAs for 145 DE-lncRNAs, including 104 cis-regulatory, 60 trans-regulatory, and three both cis-regulatory and trans-regulatory manner. Combined with our previous miRNA sequencing data, 241 DE-lncRNAs functioned as potential endogenous target mimics (eTMs) for 84 differentially expressed microRNAs (DE-miRNAs, including 12 novel miRNAs). The results of transient transformation in tobacco leaves indicated that L006278 could bind to MTCONS_00006277, which encoded a calcineurin CBL-interacting protein kinase 19-like, and suppress its expression. Furthermore, L117735 could function as an eTM for tae-miR9657b-3p, and L056972 could function as an eTM for gc-m2240-5p. To explore the function of lncRNAs in the process of Gc action, we transformed L006278, an up-regulated lncRNA in CS-3C, into rice to analyze its effect on pollen fertility. Overexpression of L006278 led to a reduction in rice pollen fertility. Overall, our findings indicate that lncRNAs can contribute to the regulation of pollen fertility during the process of Gc action by regulating the expression levels of target mRNAs and acting as eTMs for certain key miRNAs.

Construction of a Dataset for All Expressed Transcripts for Alzheimer's Disease Research

Accurate identification and functional annotation of splicing isoforms and non-coding RNAs (lncRNAs), alongside full-length protein-encoding transcripts, are critical for understanding gene (mis)regulation and metabolic reprogramming in Alzheimer's disease (AD). This study aims to provide a comprehensive and accurate transcriptome resource to improve existing AD transcript databases. Background/Objectives: Gene mis-regulation and metabolic reprogramming play a key role in AD, yet existing transcript databases lack accurate and comprehensive identification of splicing isoforms and lncRNAs. This study aims to generate a refined transcriptome dataset, expanding the understanding of AD onset and progression. Methods: Publicly available RNA-seq data from pre-AD and AD tissues were utilized. Advanced bioinformatics tools were applied to assemble and annotate full-length transcripts, including splicing isoforms and lncRNAs, with an emphasis on correcting errors and enhancing annotation accuracy. Results: A significantly improved transcriptome dataset was generated, which includes detailed annotations of splicing isoforms and lncRNAs. This dataset expands the scope of existing AD transcript databases and provides new insights into the molecular mechanisms underlying AD. The findings demonstrate that the refined dataset captures more relevant details about AD progression compared to publicly available data. Conclusions: The newly developed transcriptome resource and the associated analysis tools offer a valuable contribution to AD research, providing deeper insights into the disease's molecular mechanisms. This work supports future research into gene regulation and metabolic reprogramming in AD and serves as a foundation for exploring novel therapeutic targets.

Elucidating the Functional Roles of Long Non-Coding RNAs in Alzheimer's Disease

Alzheimer's disease (AD) is a multifaceted neurodegenerative disorder characterized by cognitive decline and neuronal loss, representing a most challenging health issue. We present a computational analysis of transcriptomic data of AD tissues vs. healthy controls, focused on the elucidation of functional roles played by long non-coding RNAs (lncRNAs) throughout the AD progression. We first assembled our own lncRNA transcripts from the raw RNA-Seq data generated from 527 samples of the dorsolateral prefrontal cortex, resulting in the identification of 31,574 novel lncRNA genes. Based on co-expression analyses between mRNAs and lncRNAs, a co-expression network was constructed. Maximal subnetworks with dense connections were identified as functional clusters. Pathway enrichment analyses were conducted over mRNAs and lncRNAs in each cluster, which served as the basis for the inference of functional roles played by lncRNAs involved in each of the key steps in an AD development model that we have previously built based on transcriptomic data of protein-encoding genes. Detailed information is presented about the functional roles of lncRNAs in activities related to stress response, reprogrammed metabolism, cell polarity, and development. Our analyses also revealed that lncRNAs have the discerning power to distinguish between AD samples of each stage and healthy controls. This study represents the first of its kind.

Transcriptional Modulation during Photomorphogenesis in Rice Seedlings

Light is one of the most important factors regulating plant gene expression patterns, metabolism, physiology, growth, and development. To explore how light may induce or alter transcript splicing, we conducted RNA-Seq-based transcriptome analyses by comparing the samples harvested as etiolated seedlings grown under continuous dark conditions vs. the light-treated green seedlings. The study aims to reveal differentially regulated protein-coding genes and novel long noncoding RNAs (lncRNAs), their light-induced alternative splicing, and their association with biological pathways. We identified 14,766 differentially expressed genes, of which 4369 genes showed alternative splicing. We observed that genes mapped to the plastid-localized methyl-erythritol-phosphate (MEP) pathway were light-upregulated compared to the cytosolic mevalonate (MVA) pathway genes. Many of these genes also undergo splicing. These pathways provide crucial metabolite precursors for the biosynthesis of secondary metabolic compounds needed for chloroplast biogenesis, the establishment of a successful photosynthetic apparatus, and photomorphogenesis. In the chromosome-wide survey of the light-induced transcriptome, we observed intron retention as the most predominant splicing event. In addition, we identified 1709 novel lncRNA transcripts in our transcriptome data. This study provides insights on light-regulated gene expression and alternative splicing in rice.

Lectin Receptor-Like Protein Kinase OsNRFG6 is Required for Embryo Sac Development and Fertilization in Neo-Tetraploid Rice

Great yield-enhancing prospects of autotetraploid rice was restricted by various polyploidy-induced reproductive dysfunction. To surmount these challenges, our group has generated a series of valuable fertile tetraploid lines (denoted as neo-tetraploid rice) through 20-year efforts. With this context, a G-type lectin receptor-like kinase, OsNRFG6, was identified as a pivotal factor associated with reproductive regulation in neo-tetraploid rice. Nevertheless, it is still elusive about a comprehensive understanding of its precise functional roles and underlying molecular mechanisms during reproduction of neo-tetraploid rice. Here, we demonstrated that OsNRFG6 executed a constitutive expression pattern and encoded proteins localizing in perinucleus and endoplasmic reticulum. Subsequently, four independent mutant lines of OsNRFG6 within neo-tetraploid rice background were further identified, all displaying low seed-setting rate due to abortive embryo sacs and defective double fertilization. RNA-seq and RT-qPCR revealed a significant down-regulation of OsNRFG6 and female reproductive genes such as OsMEL1 and LOG in ovaries prior to and post-fertilization, attributing this effect to OsNRFG6 mutation. Furthermore, through yeast-two hybrids, bimolecular fluorescence complementation assays, and luciferase complementation imaging assays, it was determined that OsNRFG6 could interact with itself and two female reproductive proteins (LOG and OsDES1) to form protein complexes. These results elucidate the reproductive functions and molecular pathway governed by OsNRFG6 in regulating fertility of neo-tetraploid rice, offering insights into molecular understanding of fertility improvement in polyploid rice.

Seed development-related genes contribute to high yield heterosis in integrated utilization of elite autotetraploid and neo-tetraploid rice

Introduction

Autotetraploid rice holds high resistance to abiotic stress and substantial promise for yield increase, but it could not be commercially used because of low fertility. Thus, our team developed neo-tetraploid rice with high fertility and hybrid vigor when crossed with indica autotetraploid rice. Despite these advances, the molecular mechanisms underlying this heterosis remain poorly understood.

Methods

An elite indica autotetraploid rice line (HD11) was used to cross with neo-tetraploid rice, and 34 hybrids were obtained to evaluate agronomic traits related to yield. WE-CLSM, RNA-seq, and CRISPR/Cas9 were employed to observe endosperm structure and identify candidate genes from two represent hybrids.

Results and discussion

These hybrids showed high seed setting and an approximately 55% increase in 1000-grain weight, some of which achieved grain yields comparable to those of the diploid rice variety. The endosperm observations indicated that the starch grains in the hybrids were more compact than those in paternal lines. A total of 119 seed heterosis related genes (SHRGs) with different expressions were identified, which might contribute to high 1000-grain weight heterosis in neo-tetraploid hybrids. Among them, 12 genes had been found to regulate grain weight formation, including OsFl3, ONAC023, OsNAC024, ONAC025, ONAC026, RAG2, FLO4, FLO11, OsISA1, OsNF-YB1, NF-YC12, and OsYUC9. Haplotype analyses of these 12 genes revealed the various effects on grain weight among different haplotypes. The hybrids could polymerize more dominant haplotypes of above grain weight regulators than any homozygous cultivar. Moreover, two SHRGs (OsFl3 and SHRG2) mutants displayed a significant reduction in 1000-grain weight and an increase in grain chalkiness, indicating that OsFl3 and SHRG2 positively regulate grain weight. Our research has identified a valuable indica autotetraploid germplasm for generating strong yield heterosis in combination with neo-tetraploid lines and gaining molecular insights into the regulatory processes of heterosis in tetraploid rice.

Update on functional analysis of long non-coding RNAs in common crops

With the rapid advances in next-generation sequencing technology, numerous non-protein-coding transcripts have been identified, including long noncoding RNAs (lncRNAs), which are functional RNAs comprising more than 200 nucleotides. Although lncRNA-mediated regulatory processes have been extensively investigated in animals, there has been considerably less research on plant lncRNAs. Nevertheless, multiple studies on major crops showed lncRNAs are involved in crucial processes, including growth and development, reproduction, and stress responses. This review summarizes the progress in the research on lncRNA roles in several major crops, presents key strategies for exploring lncRNAs in crops, and discusses current challenges and future prospects. The insights provided in this review will enhance our comprehension of lncRNA functions in crops, with potential implications for improving crop genetics and breeding.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Alleviated lead toxicity in rice plant by co-augmented action of genome doubling and TiO2 nanoparticles on gene expression, cytological and physiological changes

Lead is a very toxic and futile heavy metal for rice plants because of its injurious effects on plant growth and metabolic processes. Polyploidy or whole genome doubling increases the ability of plants to withstand biotic and abiotic stress. Considering the beneficial effects of nanoparticles and tetraploid rice, this research was conducted to examine the effectiveness of tetraploid and titanium dioxide nanoparticles (TiO2 NPs) in mitigating the toxic effects of lead. A diploid (E22-2x) and it's tetraploid (T-42) rice line were treated with Pb (200 μM) and TiO2 NPs (15 mg L-1). Lead toxicity dramatically reduced shoot length (16 % and 4 %) and root length (17 % and 9 %), biological yield (55 % and 36 %), and photosynthetic activity, as evidenced by lower levels of chlorophyll a and b (30 % and 9 %) in E-22 and T-42 rice cultivars compared to the control rice plants, respectively. Furthermore, lead toxicity amplified the levels of reactive oxygen species (ROS), such as malondialdehyde and H2O2, while decreasing activities of all antioxidant enzymes, such as superoxidase, peroxidase, and glutathione predominately in the diploid cultivar. Transmission electron microscopy and semi-thin section observations revealed that Pb-treated cells in E22-2x had more cell abnormalities than T-42, such as irregularly shaped mitochondria, cell wall, and reduced root cell size. Polyploidy and TiO2 reduced Pb uptake in rice cultivars and expression levels of metal transporter genes such as OsHMA9 and OsNRAMP5. According to the findings, genome doubling alleviates Pb toxicity by reducing Pb accumulation, ROS, and cell damage. Tetraploid rice can withstand the toxic effect of Pb better than diploid rice, and TiO2 NPs can alleviate the toxic impact of Pb. Our study findings act as a roadmap for future research endeavours, directing the focus toward risk management and assessing long-term impacts to balance environmental sustainability and agricultural growth.

Transcription Factors and Their Regulatory Roles in the Male Gametophyte Development of Flowering Plants

Male gametophyte development in plants relies on the functions of numerous genes, whose expression is regulated by transcription factors (TFs), non-coding RNAs, hormones, and diverse environmental stresses. Several excellent reviews are available that address the genes and enzymes associated with male gametophyte development, especially pollen wall formation. Growing evidence from genetic studies, transcriptome analysis, and gene-by-gene studies suggests that TFs coordinate with epigenetic machinery to regulate the expression of these genes and enzymes for the sequential male gametophyte development. However, very little summarization has been performed to comprehensively review their intricate regulatory roles and discuss their downstream targets and upstream regulators in this unique process. In the present review, we highlight the research progress on the regulatory roles of TF families in the male gametophyte development of flowering plants. The transcriptional regulation, epigenetic control, and other regulators of TFs involved in male gametophyte development are also addressed.

The applications of CRISPR/Cas-mediated microRNA and lncRNA editing in plant biology: shaping the future of plant non-coding RNA research

MAIN CONCLUSION

CRISPR/Cas technology has greatly facilitated plant non-coding RNA (ncRNA) biology research, establishing itself as a promising tool for ncRNA functional characterization and ncRNA-mediated plant improvement. Throughout the last decade, the promising genome editing tool clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated proteins (Cas; CRISPR/Cas) has allowed unprecedented advances in the field of plant functional genomics and crop improvement. Even though CRISPR/Cas-mediated genome editing system has been widely used to elucidate the biological significance of a number of plant protein-coding genes, this technology has been barely applied in the functional analysis of those non-coding RNAs (ncRNAs) that modulate gene expression, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). Nevertheless, compelling findings indicate that CRISPR/Cas-based ncRNA editing has remarkable potential for deciphering the biological roles of ncRNAs in plants, as well as for plant breeding. For instance, it has been demonstrated that CRISPR/Cas tool could overcome the challenges associated with other approaches employed in functional genomic studies (e.g., incomplete knockdown and off-target activity). Thus, in this review article, we discuss the current status and progress of CRISPR/Cas-mediated ncRNA editing in plant science in order to provide novel prospects for further assessment and validation of the biological activities of plant ncRNAs and to enhance the development of ncRNA-centered protocols for crop improvement.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice

BACKGROUND

Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice.

RESULTS

H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21.

CONCLUSION

Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

CircPCBL: Identification of Plant CircRNAs with a CNN-BiGRU-GLT Model

Circular RNAs (circRNAs), which are produced post-splicing of pre-mRNAs, are strongly linked to the emergence of several tumor types. The initial stage in conducting follow-up studies involves identifying circRNAs. Currently, animals are the primary target of most established circRNA recognition technologies. However, the sequence features of plant circRNAs differ from those of animal circRNAs, making it impossible to detect plant circRNAs. For example, there are non-GT/AG splicing signals at circRNA junction sites and few reverse complementary sequences and repetitive elements in the flanking intron sequences of plant circRNAs. In addition, there have been few studies on circRNAs in plants, and thus it is urgent to create a plant-specific method for identifying circRNAs. In this study, we propose CircPCBL, a deep-learning approach that only uses raw sequences to distinguish between circRNAs found in plants and other lncRNAs. CircPCBL comprises two separate detectors: a CNN-BiGRU detector and a GLT detector. The CNN-BiGRU detector takes in the one-hot encoding of the RNA sequence as the input, while the GLT detector uses k-mer (k = 1 - 4) features. The output matrices of the two submodels are then concatenated and ultimately pass through a fully connected layer to produce the final output. To verify the generalization performance of the model, we evaluated CircPCBL using several datasets, and the results revealed that it had an F1 of 85.40% on the validation dataset composed of six different plants species and 85.88%, 75.87%, and 86.83% on the three cross-species independent test sets composed of Cucumis sativus, Populus trichocarpa, and Gossypium raimondii, respectively. With an accuracy of 90.9% and 90%, respectively, CircPCBL successfully predicted ten of the eleven circRNAs of experimentally reported Poncirus trifoliata and nine of the ten lncRNAs of rice on the real set. CircPCBL could potentially contribute to the identification of circRNAs in plants. In addition, it is remarkable that CircPCBL also achieved an average accuracy of 94.08% on the human datasets, which is also an excellent result, implying its potential application in animal datasets. Ultimately, CircPCBL is available as a web server, from which the data and source code can also be downloaded free of charge.

Research progress on the roles of lncRNAs in plant development and stress responses

Long non-coding RNAs (lncRNAs) are RNAs of more than 200 nucleotides in length that are not (or very rarely) translated into proteins. In eukaryotes, lncRNAs regulate gene expression at the transcriptional, post-transcriptional, and epigenetic levels. lncRNAs are categorized according to their genomic position and molecular mechanism. This review summarized the characteristics and mechanisms of plant lncRNAs involved in vegetative growth, reproduction, and stress responses. Our discussion and model provide a theoretical basis for further studies of lncRNAs in plant breeding.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

Identification and Functional Analysis of Long Non-Coding RNA (lncRNA) in Response to Seed Aging in Rice

Many lncRNAs have been shown to play a vital role in aging processes. However, how lncRNAs regulate seed aging remains unknown. In this study, we performed whole transcriptome strand-specific RNA sequencing of samples from rice embryos, analyzed the differences in expression of rice seed lncRNAs before and after artificial aging treatment (AAT), and systematically screened 6002 rice lncRNAs. During the AAT period, the expression levels of most lncRNAs (454) were downregulated and only four were upregulated among the 458 differentially expressed lncRNAs (DELs). Cis- or trans-regulated target genes of the four upregulated lncRNAs were mainly related to base repair, while 454 downregulated lncRNAs were related to plant-pathogen interaction, plant hormones, energy metabolism, and secondary metabolism. The pathways of DEL target genes were similar with those of differentially expressed mRNAs (DEGs). A competing endogenous RNA (ceRNA) network composed of 34 lncRNAs, 24 microRNAs (miRNA), and 161 mRNAs was obtained. The cDNA sequence of lncRNA LNC_037529 was obtained by rapid amplification of cDNA ends (RACE) cloning with a total length of 1325 bp, a conserved 5' end, and a non-conserved 3' end. Together, our findings indicate that genome-wide selection for lncRNA downregulation was an important mechanism for rice seed aging. LncRNAs can be used as markers of seed aging in rice. These findings provide a future path to decipher the underlying mechanism associated with lncRNAs in seed aging.

Role of long non-coding RNAs in rice reproductive development

Rice is a staple crop, feeding over half of the global population. The future demand of population growth and climate change requires substantial rice improvement. Recent advances in rice genomics have highlighted the vital role of the non-coding part of the genome. The protein-coding regions account for only a tiny portion of the eukaryotic genome, and most of the genomic regions transcribe copious amounts of non-coding RNAs. Of these, the long non-coding RNAs, including linear non-coding RNAs (lncRNAs) and circular non-coding RNAs (circRNAs), have been shown to play critical roles in various developmental processes by regulating the expression of genes and functions of proteins at transcriptional, post-transcriptional and post-translational levels. With the advances in next-generation sequencing technologies, a substantial number of long non-coding RNAs have been found to be expressed in plant reproductive organs in a cell- and tissue-specific manner suggesting their reproductive development-related functions. Accumulating evidence points towards the critical role of these non-coding RNAs in flowering, anther, and pollen development, ovule and seed development and photoperiod and temperature regulation of male fertility. In this mini review, we provide a brief overview of the role of the linear and circular long non-coding RNAs in rice reproductive development and control of fertility and crop yield.

Regulatory dynamics of gene expression in the developing male gametophyte of Arabidopsis

Sexual reproduction in angiosperms requires the production and delivery of two male gametes by a three-celled haploid male gametophyte. This demands synchronized gene expression in a short developmental window to ensure double fertilization and seed set. While transcriptomic changes in developing pollen are known for Arabidopsis, no studies have integrated RNA and proteomic data in this model. Further, the role of alternative splicing has not been fully addressed, yet post-transcriptional and post-translational regulation may have a key role in gene expression dynamics during microgametogenesis. We have refined and substantially updated global transcriptomic and proteomic changes in developing pollen for two Arabidopsis accessions. Despite the superiority of RNA-seq over microarray-based platforms, we demonstrate high reproducibility and comparability. We identify thousands of long non-coding RNAs as potential regulators of pollen development, hundreds of changes in alternative splicing and provide insight into mRNA translation rate and storage in developing pollen. Our analysis delivers an integrated perspective of gene expression dynamics in developing Arabidopsis pollen and a foundation for studying the role of alternative splicing in this model.

Population transcriptomic analysis identifies the comprehensive lncRNAs landscape of spike in wheat (Triticum aestivum L.)

BACKGROUND

Long noncoding RNAs (lncRNAs) are emerging as the important regulators involving in growth and development as well as stress response in plants. However, current lncRNA studies were mainly performed at the individual level and the significance of it is not well understood in wheat.

RESULTS

In this study, the lncRNA landscape of wheat spike was characterized through analysing a total of 186 spike RNA-seq datasets from 93 wheat genotypes. A total of 35,913 lncRNAs as well as 1,619 lncRNA-mRNA pairs comprised of 443 lncRNAs and 464 mRNAs were obtained. Compared to coding genes, these lncRNAs displayed rather low conservation among wheat and other gramineous species. Based on re-sequencing data, the genetic variations of these lncRNA were investigated and obvious genetic bottleneck were found on them during wheat domestication process. Furthermore, 122 lncRNAs were found to act as ceRNA to regulate endogenous competition. Finally, association and co-localization analysis of the candidate lncRNA-mRNA pairs identified 170 lncRNAs and 167 target mRNAs significantly associated with spike-related traits, including lncRNA.127690.1/TraesCS2A02G518500.1 (PMEI) and lncRNA.104854.1/TraesCS6A02G050300.1 (ATG5) associated with heading date and spike length, respectively.

CONCLUSIONS

This study reported the lncRNA landscape of wheat spike through the population transcriptome analysis, which not only contribute to better understand the wheat evolution from the perspective of lncRNA, but also lay the foundation for revealing roles of lncRNA playing in spike development.

A Long Noncoding RNA Derived from lncRNA-mRNA Networks Modulates Seed Vigor

The discovery of long noncoding RNAs (lncRNAs) has filled a great gap in our understanding of posttranscriptional gene regulation in a variety of biological processes related to plant stress responses. However, systematic analyses of the lncRNAs expressed in rice seeds that germinate under cold stress have been elusive. In this study, we performed strand-specific whole transcriptome sequencing in germinated rice seeds under cold stress and normal temperature. A total of 6258 putative lncRNAs were identified and expressed in a stage-specific manner compared to mRNA. By investigating the targets of differentially expressed (DE) lncRNAs of LT-I (phase I of low temperature)/NT-I (phase I of normal temperature), it was shown that the auxin-activated signaling pathway was significantly enriched, and twenty-three protein-coding genes with most of the members of the SAUR family located in chromosome 9 were identified as the candidate target genes that may interact with five lncRNAs. A seed vigor-related lncRNA, SVR, which interplays with the members of the SAUR gene family in cis was eventually identified. The CRISPR/Cas 9 engineered mutations in SVR cause delay of germination. The findings provided new insights into the connection between lncRNAs and the auxin-activated signaling pathway in the regulation of rice seed vigor.

Thermo-Sensitive Genic Male Sterile Lines of Neo-Tetraploid Rice Developed through Gene Editing Technology Revealed High Levels of Hybrid Vigor

Neo-tetraploid rice, which developed from the progenies of autotetraploid hybrid by our research group, is a useful germplasm with high fertility and strong heterosis when they crossed with other autotetraploid rice lines. The CRISPR/Cas9-mediated TMS5 gene editing system has been widely used in diploid rice, but there are few reports in tetraploid rice. Here, we used CRISPR/Cas9 technology to edit the TMS5 gene, which is a temperature sensitive gene controlling the fertility in diploid rice, in neo-tetraploid rice to develop male sterile lines. Two mutant lines, H2s and H3s, were developed from the gene editing and displayed characteristics of thermo-sensitive genic male sterility. The daily mean temperatures of 23 °C to 26 °C were found to be critical for sterility (restrictive temperature) in H2s and H3s under both controlled (growth chambers) and natural growing conditions (field). Cytological observation showed the anther dysplasia appeared later in H2s and H3s than that of the TMS5 mutant of diploid rice (E285s) under the same conditions. Then these mutant lines, H2s and H3s, were crossed with tetraploid rice to generate F₁ hybrids, which exhibited obvious advantages for effective number of panicles, total grains and seed setting. The high levels of hybrids heterosis were maintained for several generations that can save seed cost. Our research provides an effective way of developing thermo-sensitive genic male sterility (TGMS) lines of tetraploid rice using gene editing, which will accelerate the utilization of polyploid heterosis.

RNA-Seq with a novel glabrous-ZM24fl reveals some key lncRNAs and the associated targets in fiber initiation of cotton

BACKGROUND

Cotton fiber is an important natural resource for textile industry and an excellent model for cell biology study. Application of glabrous mutant cotton and high-throughput sequencing facilitates the identification of key genes and pathways for fiber development and cell differentiation and elongation. LncRNA is a type of ncRNA with more than 200 nt in length and functions in the ways of chromatin modification, transcriptional and post-transcriptional modification, and so on. However, the detailed lncRNA and associated mechanisms for fiber initiation are still unclear in cotton.

RESULTS

In this study, we used a novel glabrous mutant ZM24fl, which is endowed with higher somatic embryogenesis, and functions as an ideal receptor for cotton genetic transformation. Combined with the high-throughput sequencing, fatty acid pathway and some transcription factors such as MYB, ERF and bHLH families were identified the important roles in fiber initiation; furthermore, 3,288 lncRNAs were identified, and some differentially expressed lncRNAs were also analyzed. From the comparisons of ZM24_0 DPA vs ZM24_-2 DPA and fl_0 DPA vs ZM24_0 DPA, one common lncRNA MSTRG 2723.1 was found that function upstream of fatty acid metabolism, MBY25-mediating pathway, and pectin metabolism to regulate fiber initiation. In addition, other lncRNAs MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1 were also showed potential important roles in fiber development; and the co-expression analysis between lncRNAs and targets showed the distinct models of different lncRNAs and complicated interaction between lncRNAs in fiber development of cotton.

CONCLUSIONS

From the above results, a key lncRNA MSTRG 2723.1 was identified that might mediate some key genes transcription of fatty acid metabolism, MYB25-mediating pathway, and pectin metabolism to regulate fiber initiation of ZM24 cultivar. Co-expression analysis implied that some other important lncRNAs (e.g., MSTRG 3390.1, MSTRG 48719.1, and MSTRG 31176.1) were also showed the different regulatory model and interaction between them, which proposes some valuable clues for the lncRNAs associated mechanisms in fiber development.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers

Nitrate (NO₃ ^-) and ammonium (NH₄ ⁺) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO₃ ^- and NH₄NO₃ treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH₄ ⁺ and NH₄NO₃ treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

The emerging role of small RNAs in ovule development, a kind of magic

In plants, small RNAs have been recognized as key genetic and epigenetic regulators of development. Small RNAs are usually 20 to 30 nucleotides in length and they control, in a sequence specific manner, the transcriptional or post-transcriptional expression of genes. In this review, we present a comprehensive overview of the most recent findings about the function of small RNAs in ovule development, including megasporogenesis and megagametogenesis, both in sexual and apomictic plants. We discuss recent studies on the role of miRNAs, siRNAs and trans-acting RNAs (ta-siRNAs) in early female germline differentiation. The mechanistic complexity and unique regulatory features are reviewed, and possible directions for future research are provided.

Integrated Analysis of lncRNA and mRNA Reveals Novel Insights into Wool Bending in Zhongwei Goat

Simple Summary

The high-quality lambskin of the Chinese Zhongwei goat has a high economic value. The quality of lamb skin is mainly affected by the curvature of the wool, which is regulated by the growth and development of hair follicles. In this study, the expression profiles of long non-coding RNAs (lncRNAs) of 45-day-old and 108-day-old Zhongwei goats were constructed by the Ribo Zero RNA sequencing. A total of 60 differential lncRNAs and 352 differential mRNAs were identified. The functional annotation of differential lncRNAs target genes showed that they were mainly enriched in PI3K-Akt and Arachidonic acid metabolic signaling pathways. In combination with qRT-PCR and WGCNA results, we speculate that LOC102172600 and LOC102191729 might affect hair follicle development and wool curvature by regulating the target genes. These results provide new insights into the potential role of lncRNA in regulating wool bending.

Abstract

Chinese Zhongwei goat is a rare and precious fur breed as its lamb fur is a well-known fur product. Wool bending of lamb fur of the Zhongwei goat is its most striking feature. However, the curvature of the wool decreases gradually with growth, which significantly affects its quality and economic value. The mechanism regulating the phenotypic changes of hair bending is still unclear. In the present study, the skin tissues of Zhongwei goats at 45 days (curving wool) and 108 days (slight-curving wool) after birth were taken as the research objects, and the expression profiling of long non-coding RNAs (lncRNAs) and mRNAs were analyzed based on the Ribo Zero RNA sequencing (RNA-seq) method. In total, 46,013 mRNAs and 13,549 lncRNAs were identified, of which 352 were differentially expressed mRNAs and 60 were. lncRNAs. Functional enrichment analysis of the target genes of lncRNAs were mainly enriched in PI3K-Akt, Arachidonic acid metabolic, cAMP, Wnt, and other signaling pathways. The qRT-PCR results of eight selected lncRNAs and target genes were consistent with the sequencing result, which indicated our data were reliable. Through the analysis of the weighted gene co-expression network, 13 co-expression modules were identified. The turquoise module contained a large number of differential expressed lncRNAs, which were mainly enriched in the PI3K-Akt signaling pathway and cAMP signaling pathway. The predicted LOC102172600 and LOC102191729 might affect the development of hair follicles and the curvature of wool by regulating the target genes. Our study provides novel insights into the potential roles of lncRNAs in the regulation of wool bending. In addition, the study offers a theoretical basis for further study of goat wool growth, so as to be a guidance and reference for breeding and improvement in the future.

License to Regulate: Noncoding RNA Special Agents in Plant Meiosis and Reproduction

The complexity of the subcellular processes that take place during meiosis requires a significant remodeling of cellular metabolism and dynamic changes in the organization of chromosomes and the cytoskeleton. Recently, investigations of meiotic transcriptomes have revealed additional noncoding RNA factors (ncRNAs) that directly or indirectly influence the course of meiosis. Plant meiosis is the point at which almost all known noncoding RNA-dependent regulatory pathways meet to influence diverse processes related to cell functioning and division. ncRNAs have been shown to prevent transposon reactivation, create germline-specific DNA methylation patterns, and affect the expression of meiosis-specific genes. They can also influence chromosome-level processes, including the stimulation of chromosome condensation, the definition of centromeric chromatin, and perhaps even the regulation of meiotic recombination. In many cases, our understanding of the mechanisms underlying these processes remains limited. In this review, we will examine how the different functions of each type of ncRNA have been adopted in plants, devoting attention to both well-studied examples and other possible functions about which we can only speculate for now. We will also briefly discuss the most important challenges in the investigation of ncRNAs in plant meiosis.

A Genome Doubling Event Reshapes Rice Morphology and Products by Modulating Chromatin Signatures and Gene Expression Profiling

Evolutionarily, polyploidy represents a smart method for adjusting agronomically important in crops through impacts on genomic abundance and chromatin condensation. Autopolyploids have a relatively concise genetic background with great diversity and provide an ideal system to understand genetic and epigenetic mechanisms attributed to the genome-dosage effect. However, whether and how genome duplication events during autopolyploidization impact chromatin signatures are less understood in crops. To address it, we generated an autotetraploid rice line from a diploid progenitor, Oryza sativa ssp. indica 93-11. Using transposase-accessible chromatin sequencing, we found that autopolyploids lead to a higher number of accessible chromatin regions (ACRs) in euchromatin, most of which encode protein-coding genes. As expected, the profiling of ACR densities supported that the effect of ACRs on transcriptional gene activities relies on their positions in the rice genome, regardless of genome doubling. However, we noticed that genome duplication favors genic ACRs as the main drivers of transcriptional changes. In addition, we probed intricate crosstalk among various kinds of epigenetic marks and expression patterns of ACR-associated gene expression in both diploid and autotetraploid rice plants by integrating multiple-omics analyses, including chromatin immunoprecipitation sequencing and RNA-seq. Our data suggested that the combination of H3K36me2 and H3K36me3 may be associated with dynamic perturbation of ACRs introduced by autopolyploidization. As a consequence, we found that numerous metabolites were stimulated by genome doubling. Collectively, our findings suggest that autotetraploids reshape rice morphology and products by modulating chromatin signatures and transcriptional profiling, resulting in a pragmatic means of crop genetic improvement.

Identification of stable pollen development related reference genes for accurate qRT-PCR analysis and morphological variations in autotetraploid and diploid rice

Autotetraploid rice exhibited hybrid vigor and greater genetic variation compared to diploid rice, but low pollen fertility is a major hindrance for its utilization. Our previous analysis revealed that large number of pollen fertility genes were exhibited down-regulation in autotetraploid rice. Hence, it is of utmost importance to reveal the expression patterns of pollen fertility genes with high accuracy. To find stable reference genes for autotetraploid rice, we compared the pollen development stages between diploid and autotetraploid rice, and 14 candidate genes were selected based on transcriptome analysis to evaluate their expression levels. Autotetraploid rice (i.e. Taichung65-4x) displayed lower seed set (40.40%) and higher percentage of abnormalities during the pollen development process than its diploid counterpart. To detect the candidate reference genes for pollen development of autotetraploid and diploid rice, we used five different algorithms, including NormFinder, BestKeeper, ΔCt method, geNorm and Re-Finder to evaluate their expression patterns stability. Consequently, we identified two genes, Cytochrome b5 and CPI, as the best candidate reference genes for qRT-PCR normalization in autotetraploid and diploid rice during pre-meiosis, meiosis, single microspore and bicellular pollen development stages. However, Cytochrome b5 was found to be the most stably expressed gene during different pollen development stages in autotetraploid rice. The results of our study provide a platform for subsequent gene expression analyses in autotetraploid rice, which could also be used in other polyploid plants.

A New Way of Rice Breeding: Polyploid Rice Breeding

Polyploid rice, first discovered by Japanese scientist Eiiti Nakamori in 1933, has a history of nearly 90 years. In the following years, polyploid rice studies have mainly focused on innovations in breeding theory, induction technology and the creation of new germplasm, the analysis of agronomic traits and nutritional components, the study of gametophyte development and reproduction characteristics, DNA methylation modification and gene expression regulation, distant hybridization and utilization among subspecies, species and genomes. In recent years, PMeS lines and neo-tetraploid rice lines with stable high seed setting rate characteristics have been successively selected, breaking through the bottleneck of low seed setting rate of polyploid rice. Following, a series of theoretical and applied studies on high seed setting rate tetraploid rice were carried out. This has pushed research on polyploid rice to a new stage, opening new prospects for polyploid rice breeding.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered “junk”, DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these “junk” RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2n = 6x = 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of “junk” DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10–30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.

Fertile Tetraploids: New Resources for Future Rice Breeding?

Ploidy manipulation is an efficient technique for the development of novel phenotypes in plant breeding. However, in rice (Oryza sativa L.), severe seed sterility has been considered a barrier preventing cultivation of autotetraploids since the 1930s. Recently, a series of studies identified two fertile autotetraploids, identified herein as the PMeS (Polyploid Meiosis Stability) and Neo-Tetraploid lines. Here, we summarize their characteristics, focusing on the recovery of seed fertility, and discuss potential future directions of study in this area, providing a comprehensive understanding of current progress in the study of fertile tetraploid rice, a classical, but promising, concept for rice breeding.