Genome-wide analyses reveal the role of noncoding variation in complex traits during rice domestication

Evolutionary comparison of lncRNAs in four cotton species and functional identification of LncR4682-PAS2-KCS19 module in fiber elongation

Long non-coding RNAs (lncRNAs) play an important role in various biological processes in plants. However, there have been few reports on the evolutionary signatures of lncRNAs in closely related cotton species. The lncRNA transcription patterns in two tetraploid cotton species and their putative diploid ancestors were compared in this paper. By performing deep RNA sequencing, we identified 280 429 lncRNAs from 21 tissues in four cotton species. lncRNA transcription evolves more rapidly than mRNAs, and exhibits more severe turnover phenomenon in diploid species compared to that in tetraploid species. Evolutionarily conserved lncRNAs exhibit higher expression levels, and lower tissue specificity compared with species-specific lncRNAs. Remarkably, tissue expression of homologous lncRNAs in Gossypium hirsutum and G. barbadense exhibited similar patterns, suggesting that these lncRNAs may be functionally conserved and selectively maintained during domestication. An orthologous lncRNA, lncR4682, was identified and validated in fibers of G. hirsutum and G. barbadense with the highest conservatism and expression abundance. Through virus-induced gene silencing in upland cotton, we found that lncR4682 and its target genes GHPAS2 and GHKCS19 positively regulated fiber elongation. In summary, the present study provides a systematic analysis of lncRNAs in four closely related cotton species, extending the understanding of transcriptional conservation of lncRNAs across cotton species. In addition, LncR4682-PAS2-KCS19 contributes to cotton fiber elongation by participating in the biosynthesis of very long-chain fatty acids.

LncRNAs: the art of being influential without protein

The plant long noncoding (lnc)RNA field is on the brink of transitioning from large-scale identification of lncRNAs to their functional characterization. Due to the cross-kingdom conservation of interaction types and molecular functions, there is much to be learned from mammalian lncRNA research. Here, we discuss the different molecular processes involving lncRNAs from the regulation of chromatin to splicing. Furthermore, we discuss the lncRNA interactome, which includes proteins, other RNAs, and DNA. We explore and discuss how mammalian lncRNA functionalities could be reflected in similar pathways in plants and hypothesize that several breakthroughs in mammalian research could lead to the discovery of novel plant lncRNA molecular functions. Expanding our knowledge of the biological role of lncRNAs and their multiple applications paves the way for future agricultural applications.

Epigenetic differences between wild and cultivated grapevines highlight the contribution of DNA methylation during crop domestication

The domestication process in grapevines has facilitated the fixation of desired traits. Nowadays, vegetative propagation through cuttings enables easier preservation of these genotypes compared to sexual reproduction. Nonetheless, even with vegetative propagation, various phenotypes are often present within the same vineyard due to the accumulation of somatic mutations. These mutations are not the sole factors influencing phenotype. Alongside somatic variations, epigenetic variation has been proposed as a pivotal player in regulating phenotypic variability acquired during domestication. The emergence of these epialleles might have significantly influenced grapevine domestication over time. This study aims to investigate the impact of domestication on methylation patterns in cultivated grapevines. Reduced-representation bisulfite sequencing was conducted on 18 cultivated and wild accessions. Results revealed that cultivated grapevines exhibited higher methylation levels than their wild counterparts. Differential Methylation Analysis between wild and cultivated grapevines identified a total of 9955 differentially methylated cytosines, of which 78% were hypermethylated in cultivated grapevines. Functional analysis shows that core methylated genes (consistently methylated in both wild and cultivated accessions) are associated with stress response and terpenoid/isoprenoid metabolic processes. Meanwhile, genes with differential methylation are linked to protein targeting to the peroxisome, ethylene regulation, histone modifications, and defense response. Collectively, our results highlight the significant roles that epialleles may have played throughout the domestication history of grapevines.

LncRNA DANA1 promotes drought tolerance and histone deacetylation of drought responsive genes in Arabidopsis

Although many long noncoding RNAs have been discovered in plants, little is known about their biological function and mode of action. Here we show that the drought-induced long intergenic noncoding RNA DANA1 interacts with the L1p/L10e family member protein DANA1-INTERACTING PROTEIN 1 (DIP1) in the cell nucleus of Arabidopsis, and both DANA1 and DIP1 promote plant drought resistance. DANA1 and DIP1 increase histone deacetylase HDA9 binding to the CYP707A1 and CYP707A2 loci. DIP1 further interacts with PWWP3, a member of the PEAT complex that associates with HDA9 and has histone deacetylase activity. Mutation of DANA1 enhances CYP707A1 and CYP707A2 acetylation and expression resulting in impaired drought tolerance, in agreement with dip1 and pwwp3 mutant phenotypes. Our results demonstrate that DANA1 is a positive regulator of drought response and that DANA1 works jointly with the novel chromatin-related factor DIP1 on epigenetic reprogramming of the plant transcriptome during the response to drought.

Selection on the promoter regions plays an important role in complex traits during duck domestication

BACKGROUND

Identifying the key factors that underlie complex traits during domestication is a great challenge for evolutionary and biological studies. In addition to the protein-coding region differences caused by variants, a large number of variants are located in the noncoding regions containing multiple types of regulatory elements. However, the roles of accumulated variants in gene regulatory elements during duck domestication and economic trait improvement are poorly understood.

RESULTS

We constructed a genomics, transcriptomics, and epigenomics map of the duck genome and assessed the evolutionary forces that have been in play across the whole genome during domestication. In total, 304 (42.94%) gene promoters have been specifically selected in Pekin duck among all selected genes. Joint multi-omics analysis reveals that 218 genes (72.01%) with selected promoters are located in open and active chromatin, and 267 genes (87.83%) with selected promoters were highly and differentially expressed in domestic trait-related tissues. One important candidate gene ELOVL3, with a strong signature of differentiation on the core promoter region, is known to regulate fatty acid elongation. Functional experiments showed that the nearly fixed variants in the top selected ELOVL3 promoter in Pekin duck decreased binding ability with HLF and increased gene expression, with the overexpression of ELOVL3 able to increase lipid deposition and unsaturated fatty acid enrichment.

CONCLUSIONS

This study presents genome resequencing, RNA-Seq, Hi-C, and ATAC-Seq data of mallard and Pekin duck, showing that selection of the gene promoter region plays an important role in gene expression and phenotypic changes during domestication and highlights that the variants of the ELOVL3 promoter may have multiple effects on fat and long-chain fatty acid content in ducks.

Identification of novel loci associated with starch content in maize kernels by a genome-wide association study using an enlarged SNP panel

Starch is a major component of cereals, comprising over 70% of dry weight. It serves as a primary carbon source for humans and animals. In addition, starch is an indispensable industrial raw material. While maize (Zea mays) is a key crop and the primary source of starch, the genetic basis for starch content in maize kernels remains poorly understood. In this study, using an enlarged panel, we conducted a genome-wide association study (GWAS) based on best linear unbiased prediction (BLUP) value for starch content of 261 inbred lines across three environments. Compared with previous study, we identified 14 additional significant quantitative trait loci (QTL), encompassed a total of 42 genes, and indicated that increased marker density contributes to improved statistical power. By integrating gene expression profiling, Gene Ontology (GO) enrichment and haplotype analysis, several potential target genes that may play a role in regulating starch content in maize kernels have been identified. Notably, we found that ZmAPC4, associated with the significant SNP chr4.S_175584318, which encodes a WD40 repeat-like superfamily protein and is highly expressed in maize endosperm, might be a crucial regulator of maize kernel starch synthesis. Out of the 261 inbred lines analyzed, they were categorized into four haplotypes. Remarkably, it was observed that the inbred lines harboring hap4 demonstrated the highest starch content compared to the other haplotypes. Additionally, as a significant achievement, we have developed molecular markers that effectively differentiate maize inbred lines based on their starch content. Overall, our study provides valuable insights into the genetic basis of starch content and the molecular markers can be useful in breeding programs aimed at developing maize varieties with high starch content, thereby improving breeding efficiency.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01437-6.

Emerging roles of long non-coding RNAs in regulating agriculturally important seed traits

Seeds have enormous economic importance as source of calories, nutrition, edible oil, and biofuels. Therefore, seed traits like seed size and shape, weight, micronutrient content, oil content, quality, post-harvest sprouting, etc., are some of the main targets in crop improvement. Designing the strategies for their improvement benefits heavily from understanding the regulatory aspects of seed development. Recent studies indicate that long non-coding RNAs (lncRNAs) are one of the important regulators of seed development. They played a significant role in crop domestication by influencing seed traits. LncRNAs are conventionally defined as non-coding RNAs greater than 200 bp in length but lacking protein coding potential. Here we highlight the emerging pieces of evidence of lncRNA-mediated regulation of seed development through diverse mechanisms, for instance, by acting as target mimics or precursors of regulatory small RNAs or through chromatin remodeling and post-transcriptional repression. We also enumerate the insights from high-throughput transcriptomic studies from developing seeds of cereal, oilseed, biofuel, and pulse crops. We highlight the lncRNA candidates and lncRNA-mediated regulatory networks regulating seed development and related agronomic traits. Further, we discuss the potential of lncRNAs for improvement of agriculturally important seed traits through marker-assisted breeding and/or transgenic approaches.

qCTB7 positively regulates cold tolerance at booting stage in rice

KEY MESSAGE

LOC_Os07g07690 on qCTB7 is associated with cold tolerance at the booting stage in rice, and analysis of transgenic plants demonstrated that qCTB7 influenced cold tolerance by altering the morphology and cytoarchitecture of anthers and pollen. Cold tolerance at the booting stage (CTB) in rice can significantly affect yield in high-latitude regions. Although several CTB genes have been isolated, their ability to induce cold tolerance is insufficient to ensure adequate rice yields in cold regions at high latitudes. Here, we identified the PHD-finger domain-containing protein gene qCTB7 using QTL-seq and linkage analysis through systematic measurement of CTB differences and the spike fertility of the Longjing31 and Longdao3 cultivars, resulting in the derivation of 1570 F2 progeny under cold stress. We then characterized the function of qCTB7 in rice. It was found that overexpression of qCTB7 promoted CTB and the same yield as Longdao3 under normal growing conditions while the phenotype of qctb7 knockout showed anther and pollen failure under cold stress. When subjected to cold stress, the germination of qctb7 pollen on the stigma was reduced, resulting in lower spike fertility. These findings indicate that qCTB7 regulates the appearance, morphology, and cytoarchitecture of the anthers and pollen. Three SNPs in the promoter region and coding region of qCTB7 were identified as recognition signals for CTB in rice and could assist breeding efforts to improve cold tolerance for rice production in high latitudes.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

Identification of LncRNAs and Functional Analysis of ceRNA Related to Fatty Acid Synthesis during Flax Seed Development

Flax is a flowering plant cultivated for its oil and contains various unsaturated fatty acids. Linseed oil is known as the "deep-sea fish oil" of plants, and is beneficial to brain and blood lipids, among other positive effects. Long non-coding RNAs (lncRNAs) play an important role in plant growth and development. There are not many studies assessing how lncRNAs are related to the fatty acid synthesis of flax. The relative oil contents of the seeds of the variety Heiya NO.14 (for fiber) and the variety Macbeth (for oil) were determined at 5 day, 10 day, 20 day, and 30 day after flowering. We found that 10-20 day is an important period for ALA accumulation in the Macbeth variety. The strand-specific transcriptome data were analyzed at these four time points, and a series of lncRNAs related to flax seed development were screened. A competing endogenous RNA (ceRNA) network was constructed and the accuracy of the network was verified using qRT-PCR. MSTRG.20631.1 could act with miR156 on the same target, squamosa promoter-binding-like protein (SPL), to influence fatty acid biosynthesis through a gluconeogenesis-related pathway during flax seed development. This study provides a theoretical basis for future studies assessing the potential functions of lncRNAs during seed development.

Functional variation in the non-coding genome: molecular implications for food security

The growing world population, in combination with the anticipated effects of climate change, is pressuring food security. Plants display an impressive arsenal of cellular mechanisms conferring resilience to adverse environmental conditions, and humans rely on these mechanisms for stable food production. The elucidation of the molecular basis of the mechanisms used by plants to achieve resilience promises knowledge-based approaches to enhance food security. DNA sequence polymorphisms can reveal genomic regions that are linked to beneficial traits of plants. However, our ability to interpret how a given DNA sequence polymorphism confers a fitness advantage at the molecular level often remains poor. A key factor is that these polymorphisms largely localize to the enigmatic non-coding genome. Here, we review the functional impact of sequence variations in the non-coding genome on plant biology in the context of crop breeding and agricultural traits. We focus on examples of non-coding with particularly convincing functional support. Our survey combines findings that are consistent with the view that the non-coding genome contributes to cellular mechanisms assisting many plant traits. Understanding how DNA sequence polymorphisms in the non-coding genome shape plant traits at the molecular level offers a largely unexplored reservoir of solutions to address future challenges in plant growth and resilience.

LncRNA expression analysis by comparative transcriptomics among closely related poplars and their regulatory roles in response to salt stress

Long noncoding RNAs (lncRNAs) play crucial roles in regulating key biological processes; however, our knowledge of lncRNAs' roles in plant adaptive evolution is still limited. Here, we determined the divergence of conserved lncRNAs in closely related poplar species that were either tolerant or sensitive to salt stress by comparative transcriptome analysis. Among the 34,363 identified lncRNAs, approximately 3% were shared among poplar species with conserved sequences but diversified in their function, copy number, originating genomic region and expression patterns. Further cluster analysis revealed that the conserved lncRNAs showed more similar expression patterns within salt-tolerant poplars (P. euphratica and P. pruinosa) than between salt-tolerant and salt-sensitive poplars. Among these lncRNAs, the antisense lncRNA lncERF024 was induced by salt and differentiated expression between salt-sensitive and salt-tolerant poplars. Overexpression of lncERF024 in P. alba var. pyramidalis enhanced poplar tolerance to salt stress. Furthermore, RNA pull-down and RNA-seq analysis showed that numerous candidate genes or proteins associated with stress response and photosynthesis might be involved in salt resistance in PeulncERF024-OE poplars. Altogether, our study provided novel insight into how the diversification of lncRNA expression contributes to plant adaptation traits and showed that lncERF024 may be involved in the regulation both of gene expression and protein function conferring salt tolerance in Populus.

Combined population transcriptomic and genomic analysis reveals cis-regulatory differentiation of non-coding RNAs in maize

Long intergenic non-coding RNA (lincRNA), cis-acting expression quantitative trait locus (cis-eQTL), maize, regulatory evolution. The law of genetic variation during domestication explains the evolutionary mechanism and provides a theoretical basis for improving existing varieties of maize. Previous studies focused on exploiting regulatory variations controlling the expression of protein-coding genes rather than of non-protein-coding genes. Here, we examined the genetic and evolutionary features of long non-coding RNAs from intergenic regions (long intergenic non-coding RNAs, lincRNAs) using population-scale transcriptome data and identified 1168 lincRNAs with cis-acting expression quantitative trait loci (cis-eQTLs). We found that lincRNAs are more likely to be regulated by cis-eQTLs, which exert stronger effects than the protein-coding genes. During maize domestication and improvement, upregulated alleles of lincRNAs, which originated from both standing variation and new mutation, accumulate more frequently and show larger effect sizes than the coding genes. A stronger signature of genetic differentiation was observed in their regulatory regions compared to those of randomly sampled lincRNAs. In addition, we found that cis-regulatory differentiation of lincRNAs is related to the sequence conservation of lincRNA transcripts. Non-conserved lincRNAs more tend to gain upregulated alleles and show a stronger relationship with selected traits than conserved lincRNAs between maize and its wild relatives. Our findings in maize improve the understanding of cis-regulatory variation in lincRNA genes during domestication and improvement and provide an effective approach for prioritizing candidates for further investigation.

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.

CAN OF SPINACH, a novel long non-coding RNA, affects iron deficiency responses in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are RNA molecules with functions independent of any protein-coding potential. A whole transcriptome (RNA-seq) study of Arabidopsis shoots under iron sufficient and deficient conditions was carried out to determine the genes that are iron-regulated in the shoots. We identified two previously unannotated transcripts on chromosome 1 that are significantly iron-regulated. We have called this iron-regulated lncRNA, CAN OF SPINACH (COS). cos mutants have altered iron levels in leaves and seeds. Despite the low iron levels in the leaves, cos mutants have higher chlorophyll levels than WT plants. Moreover, cos mutants have abnormal development during iron deficiency. Roots of cos mutants are longer than those of WT plants, when grown on iron deficient medium. In addition, cos mutant plants accumulate singlet oxygen during iron deficiency. The mechanism through which COS affects iron deficiency responses is unclear, but small regions of sequence similarity to several genes involved in iron deficiency responses occur in COS, and small RNAs from these regions have been detected. We hypothesize that COS is required for normal adaptation to iron deficiency conditions.

OsWRKY115 on qCT7 links to cold tolerance in rice

qCT7, a novel QTL for increasing seedling cold tolerance in rice, was fine-mapped to a 70.9-kb region on chromosome 7, and key OsWRKY115 was identified in transgenic plants. Cold stress caused by underground cold-water irrigation seriously limits rice productivity. We systemically measured the cold-responsive traits of 2,570 F2 individuals derived from two widely cultivated rice cultivars, Kong-Yu-131 and Dong-Nong-422, to identify the major genomic regions associated with cold tolerance. A novel major QTL, qCT7, was mapped on chromosome 7 associated with the cold tolerance and survival, using whole-genome re-sequencing with bulked segregant analysis. Local QTL linkage analysis with F₂ and fine mapping with recombinant plant revealed a 70.9-kb core region on qCT7 encoding 13 protein-coding genes. Only the LOC_Os07g27670 expression level encoding the OsWRKY115 transcription factor on the locus was specifically induced by cold stress in the cold-tolerant cultivar. Moreover, haplotype analysis and the KASP8 marker indicated that OsWRKY115 was significantly associated with cold tolerance. Overexpression and knockout of OsWRKY115 significantly affected cold tolerance in seedlings. Our experiments identified OsWRKY115 as a novel regulatory gene associated with cold response in rice, and the Kong-Yu-131 allele with specific cold-induced expression may be an important molecular variant.

Cereal Endosperms: Development and Storage Product Accumulation

The persistent triploid endosperms of cereal crops are the most important source of human food and animal feed. The development of cereal endosperms progresses through coenocytic nuclear division, cellularization, aleurone and starchy endosperm differentiation, and storage product accumulation. In the past few decades, the cell biological processes involved in endosperm formation in most cereals have been described. Molecular genetic studies performed in recent years led to the identification of the genes underlying endosperm differentiation, regulatory network governing storage product accumulation, and epigenetic mechanism underlying imprinted gene expression. In this article, we outline recent progress in this area and propose hypothetical models to illustrate machineries that control aleurone and starchy endosperm differentiation, sugar loading, and storage product accumulations. A future challenge in this area is to decipher the molecular mechanisms underlying coenocytic nuclear division, endosperm cellularization, and programmed cell death.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences.

Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Compensation of Wild Plants Weakens the Effects of Crop-Wild Gene Flow on Wild Rice Populations

Crop-wild gene flow may alter the fitness of the recipient i.e., crop-wild hybrids, then potentially impact wild populations, especially for the gene flow carrying selective advantageous crop alleles, such as transgenes conferring insect resistance. Given the continuous crop-wild gene flow since crop domestication and the occasionally stressful environments, the extant wild populations of most crops are still "wild." One interpretation for this phenomenon is that wild populations have the mechanism buffered for the effects of crop alleles. However, solid evidence for this has been scarce. We used wild rice (Oryza rufipogon) and transgenic (Bt/CpTI) rice (O. sativa) as a crop-wild gene flow model and established cultivated, wild, and F7 hybrid rice populations under four levels of insect (Chilo suppressalis) pressure. Then, we measured the trait performance of the plants and estimated fitness to test the compensatory response of relatively high fitness compared to the level of insect damage. The performance of all plants varied with the insect pressure level; wild plants had higher insect-tolerance that was expressed as over- or equal-compensatory responses to insect damage, whereas crop and hybrids exhibited under-compensatory responses. The higher compensation resulted in a better performance of wild rice under insect pressure where transgenes conferring insect resistance had a somewhat beneficial effect. Remarkable hybrid vigour and the benefit effect of transgenes increased the fitness of hybrids together, but this joint effect was weakened by the compensation of wild plants. These results suggest that compensation to environmental stress may reduce the potential impacts of crop alleles on wild plants, thereby it is a mechanism maintaining the "wild" characteristics of wild populations under the scenario of continuous crop-wild gene flow.

Functional examination of lncRNAs in allotetraploid Gossypium hirsutum

BACKGROUND

An evolutionary model using diploid and allotetraploid cotton species identified 80 % of non-coding transcripts in allotetraploid cotton as being uniquely activated in comparison with its diploid ancestors. The function of the lncRNAs activated in allotetraploid cotton remain largely unknown.

RESULTS

We employed transcriptome analysis to examine the relationship between the lncRNAs and mRNAs of protein coding genes (PCGs) in cotton leaf tissue under abiotic stresses. LncRNA expression was preferentially associated with that of the flanking PCGs. Selected highly-expressed lncRNA candidates (n = 111) were subjected to a functional screening pilot test in which virus-induced gene silencing was integrated with abiotic stress treatment. From this low-throughput screen, we obtained candidate lncRNAs relating to plant height and tolerance to drought and other abiotic stresses.

CONCLUSIONS

Low-throughput screen is an effective method to find functional lncRNA for further study. LncRNAs were more active in abiotic stresses than PCG expression, especially temperature stress. LncRNA XLOC107738 may take a cis-regulatory role in response to environmental stimuli. The degree to which lncRNAs are constitutively expressed may impact expression patterns and functions on the individual gene level rather than in genome-wide aggregate.

Identification of long noncoding natural antisense transcripts (lncNATs) correlated with drought stress response in wild rice (Oryza nivara)

Background

Wild rice, including Oryza nivara and Oryza rufipogon, which are considered as the ancestors of Asian cultivated rice (Oryza sativa), possess high genetic diversity and serve as a crucial resource for breeding novel cultivars of cultivated rice. Although rice domestication related traits, such as seed shattering and plant architecture, have been intensively studied at the phenotypic and genomic levels, further investigation is needed to understand the molecular basis of phenotypic differences between cultivated and wild rice. Drought stress is one of the most severe abiotic stresses affecting rice growth and production. Adaptation to drought stress involves a cascade of genes and regulatory factors that form complex networks. O. nivara inhabits swampy areas with a seasonally dry climate, which is an ideal material to discover drought tolerance alleles. Long noncoding natural antisense transcripts (lncNATs), a class of long noncoding RNAs (lncRNAs), regulate the corresponding sense transcripts and play an important role in plant growth and development. However, the contribution of lncNATs to drought stress response in wild rice remains largely unknown.

Results

Here, we conducted strand-specific RNA sequencing (ssRNA-seq) analysis of Nipponbare (O. sativa) and two O. nivara accessions (BJ89 and BJ278) to determine the role of lncNATs in drought stress response in wild rice. A total of 1246 lncRNAs were identified, including 1091 coding–noncoding NAT pairs, of which 50 were expressed only in Nipponbare, and 77 were expressed only in BJ89 and/or BJ278. Of the 1091 coding–noncoding NAT pairs, 240 were differentially expressed between control and drought stress conditions. Among these 240 NAT pairs, 12 were detected only in Nipponbare, and 187 were detected uniquely in O. nivara. Furthermore, 10 of the 240 coding–noncoding NAT pairs were correlated with genes enriched in stress responsive GO terms; among these, nine pairs were uniquely found in O. nivara, and one pair was shared between O. nivara and Nipponbare.

Conclusion

We identified lncNATs associated with drought stress response in cultivated rice and O. nivara. These results will improve our understanding of the function of lncNATs in drought tolerance and accelerate rice breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07754-4.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

Genome-Wide Analysis Identified a Set of Conserved lncRNAs Associated with Domestication-Related Traits in Rice

Crop domestication, which gives rise to a number of desirable agronomic traits, represents a typical model system of plant evolution. Numerous genomic evidence has proven that noncoding RNAs such as microRNAs and phasiRNAs, as well as protein-coding genes, are selected during crop domestication. However, limited data shows plant long noncoding RNAs (lncRNAs) are also involved in this biological process. In this study, we performed strand-specific RNA sequencing of cultivated rice Oryza sativa ssp. japonica and O. sativa ssp. indica, and their wild progenitor O. rufipogon. We identified a total of 8528 lncRNAs, including 4072 lncRNAs in O. rufipogon, 2091 lncRNAs in japonica rice, and 2365 lncRNAs in indica rice. The lncRNAs expressed in wild rice were revealed to be shorter in length and had fewer exon numbers when compared with lncRNAs from cultivated rice. We also identified a number of conserved lncRNAs in the wild and cultivated rice. The functional study demonstrated that several of these conserved lncRNAs are associated with domestication-related traits in rice. Our findings revealed the feature and conservation of lncRNAs during rice domestication and will further promote functional studies of lncRNAs in rice.

Mutation of a major CG methylase alters genome-wide lncRNA expression in rice

Plant long non-coding RNAs (lncRNAs) function in diverse biological processes, and lncRNA expression is under epigenetic regulation, including by cytosine DNA methylation. However, it remains unclear whether 5-methylcytosine (5mC) plays a similar role in different sequence contexts (CG, CHG, and CHH). In this study, we characterized and compared the profiles of genome-wide lncRNA profiles (including long intergenic non-coding RNAs [lincRNAs] and long noncoding natural antisense transcripts [lncNATs]) of a null mutant of the rice DNA methyltransferase 1, OsMET1-2 (designated OsMET1-2-/-) and its isogenic wild type (OsMET1-2+/+). The En/Spm transposable element (TE) family, which was heavily methylated in OsMET1-2+/+, was transcriptionally de-repressed in OsMET1-2-/- due to genome-wide erasure of CG methylation, and this led to abundant production of specific lncRNAs. In addition, RdDM-mediated CHH hypermethylation was increased in the 5'-upstream genomic regions of lncRNAs in OsMET1-2-/-. The positive correlation between the expression of lincRNAs and that of their proximal protein-coding genes was also analyzed. Our study shows that CG methylation negatively regulates the TE-related expression of lncRNA and demonstrates that CHH methylation is also involved in the regulation of lncRNA expression.

Full-length annotation with multistrategy RNA-seq uncovers transcriptional regulation of lncRNAs in cotton

Long noncoding RNAs (lncRNAs) are crucial factors during plant development and environmental responses. To build an accurate atlas of lncRNAs in the diploid cotton Gossypium arboreum, we combined Isoform-sequencing, strand-specific RNA-seq (ssRNA-seq), and cap analysis gene expression (CAGE-seq) with PolyA-seq and compiled a pipeline named plant full-length lncRNA to integrate multi-strategy RNA-seq data. In total, 9,240 lncRNAs from 21 tissue samples were identified. 4,405 and 4,805 lncRNA transcripts were supported by CAGE-seq and PolyA-seq, respectively, among which 6.7% and 7.2% had multiple transcription start sites (TSSs) and transcription termination sites (TTSs). We revealed that alternative usage of TSS and TTS of lncRNAs occurs pervasively during plant growth. Besides, we uncovered that many lncRNAs act in cis to regulate adjacent protein-coding genes (PCGs). It was especially interesting to observe 64 cases wherein the lncRNAs were involved in the TSS alternative usage of PCGs. We identified lncRNAs that are coexpressed with ovule- and fiber development-associated PCGs, or linked to GWAS single-nucleotide polymorphisms. We mapped the genome-wide binding sites of two lncRNAs with chromatin isolation by RNA purification sequencing. We also validated the transcriptional regulatory role of lnc-Ga13g0352 via virus-induced gene suppression assay, indicating that this lncRNA might act as a dual-functional regulator that either activates or inhibits the transcription of target genes.

Comparative analysis of the transcriptomes of two rice subspecies during domestication

Two subspecies of rice, Oryza sativa ssp. indica and O. sativa ssp. japonica, with reproductive isolation and differences in morphology and phenotypic differences, were established during the process of rice domestication. To understand how domestication has changed the transcriptomes of the two rice subspecies and given rise to the phenotypic differences, we obtained approximately 700 Gb RNA-Seq data from 26 indica and 25 japonica accessions, and identified 97,005 transcribed fragments and 4579 novel transcriptionally active regions. The two rice subspecies had significantly different gene expression profiles, we identified 1,357 (3.3% in all genes) differentially expressed genes (DEGs) between indica and japonica rice. Combining existing gene function studies, it is found that some of these differential genes are related to the differentiation of the two subspecies, such as grain shape and cold tolerance, etc. Functional annotation of these DEGs indicates that they are involved in cell wall biosynthesis and reproductive processes. Furthermore, compared with the non-DEGs, the DEGs from both subspecies had more 5'flanking regions with low polymorphism to divergence ratios, indicating a stronger positive selection pressure on the regulation of the DEGs. This study improves our understanding of the rice genome by comparatively analyzing the transcriptomes of indica and japonica rice and identifies DEGs those may be responsible for the reproductive isolation and phenotypic differences between the two rice subspecies.

Interpopulation differences of retroduplication variations (RDVs) in rice retrogenes and their phenotypic correlations

Retroduplication variation (RDV), a type of retrocopy polymorphism, is considered to have essential biological significance, but its effect on gene function and species phenotype is still poorly understood. To this end, we analyzed the retrocopies and RDVs in 3,010 rice genomes. We calculated the RDV frequencies in the genome of each rice population; detected the mutated, ancestral and expressed retrogenes in rice genomes; and analyzed their RDV influence on rice phenotypic traits. Collectively, 73 RDVs were identified, and 14 RDVs in ancestral retrogenes can significantly affect rice phenotypes. Our research reveals that RDV plays an important role in rice migration, domestication and evolution. We think that RDV is a good molecular breeding marker candidate. To our knowledge, this is the first study on the relationship between retrogene function, expression, RDV and species phenotype.

Long Non-Coding RNAs, the Dark Matter: An Emerging Regulatory Component in Plants

Long non-coding RNAs (lncRNAs) are pervasive transcripts of longer than 200 nucleotides and indiscernible coding potential. lncRNAs are implicated as key regulatory molecules in various fundamental biological processes at transcriptional, post-transcriptional, and epigenetic levels. Advances in computational and experimental approaches have identified numerous lncRNAs in plants. lncRNAs have been found to act as prime mediators in plant growth, development, and tolerance to stresses. This review summarizes the current research status of lncRNAs in planta, their classification based on genomic context, their mechanism of action, and specific bioinformatics tools and resources for their identification and characterization. Our overarching goal is to summarize recent progress on understanding the regulatory role of lncRNAs in plant developmental processes such as flowering time, reproductive growth, and abiotic stresses. We also review the role of lncRNA in nutrient stress and the ability to improve biotic stress tolerance in plants. Given the pivotal role of lncRNAs in various biological processes, their functional characterization in agriculturally essential crop plants is crucial for bridging the gap between phenotype and genotype.

Genetic variations associated with long noncoding RNAs

Genetic variations, including single nucleotide polymorphisms (SNPs) and structural variations, are widely distributed in the genome, including the long noncoding RNA (lncRNA) regions. The changes at locus might produce numerous effects in a variety of aspects. Multiple bioinformatics resources and tools were also developed for systematically dealing with genetic variations associated with lncRNAs. Moreover, correlation of the genetic variations in lncRNAs with immune disease, cancers, and other disease as well as development process were all included for discussion. In this essay, we summarized how and in what aspects these changes would affect lncRNA functions.

Genome-wide identification and characterization of long non-coding RNAs related to grain yield in foxtail millet [Setaria italica (L.) P. Beauv.]

Background

Long noncoding RNAs (lncRNAs) have been reported to play critical roles in diverse growth and development processes in plants. However, the systematic identification and characterization of lncRNAs in foxtail millet is nearly blank.

Results

In this study, we performed high-throughput sequencing of young spikelets from four foxtail millet varieties in different yield levels at booting stage. As a result, a total of 12,378 novel lncRNAs were identified, and 70 were commonly significantly differentially expressed in comparisons between high-yield varieties and conventional varieties, suggesting that they involved in yield formation and regulation in foxtail millet. Functional analysis revealed that among the 70 significantly differentially expressed lncRNAs, 67 could transcriptionally modulate target genes in cis and in trans. Moreover, 18 lncRNAs related to grain yield in foxtail millet were predicted to function as miRNA target mimics and regulate gene expression by competing for the interaction between miRNAs and their target mRNAs.

Conclusion

Our results will provide materials for elucidation of the molecular mechanisms of lncRNAs participate in yield regulation, and will contribute to high yield foxtail millet breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07272-9.

Molecular Control and Application of Male Fertility for Two-Line Hybrid Rice Breeding

The significance of the climate change may involve enhancement of plant growth as well as utilization of the environmental alterations in male fertility (MF) regulation via male sterility (MS) systems. We described that MS systems provide a fundamental platform for improvement in agriculture production and have been explicated for creating bulk germplasm of the two-line hybrids (EGMS) in rice as compared to the three-line, to gain production sustainability and exploit its immense potential. Environmental alterations such as photoperiod and/or temperature and humidity regulate MS in EGMS lines via genetic and epigenetic changes, regulation of the noncoding RNAs, and RNA-metabolism including the transcriptional factors (TFs) implication. Herein, this article enlightens a deep understanding of the molecular control of MF in EGMS lines and exploring the regulatory driving forces that function efficiently during plant adaption under a changing environment. We highlighted a possible solution in obtaining more stable hybrids through apomixis (single-line system) for seed production.

Identifying Risk Genes and Interpreting Pathogenesis for Parkinson's Disease by a Multiomics Analysis

Genome-wide association studies (GWAS) have identified tens of genetic variants associated with Parkinson’s disease (PD). Nevertheless, the genes or DNA elements that affect traits through these genetic variations are usually undiscovered. This study was the first to combine meta-analysis GWAS data and expression data to identify PD risk genes. Four known genes, CRHR1, KANSL1, NSF and LRRC37A, and two new risk genes, STX4 and BST1, were identified. Among them, CRHR1 is a known drug target, indicating that hydrocortisone may become a potential drug for the treatment of PD. Furthermore, the potential pathogenesis of CRHR1 and LRRC37A was explored by applying DNA methylation (DNAm) data, indicating a pathogenesis whereby the effect of a genetic variant on PD is mediated by genetic regulation of transcription through DNAm. Overall, this research identified the risk genes and pathogenesis that affect PD through genetic variants, which has significance for the diagnosis and treatment of PD.

Conservation recommendations for Oryza rufipogon Griff. in China based on genetic diversity analysis

Over the past 30 years, human disturbance and habitat fragmentation have severely endangered the survival of common wild rice (Oryza rufipogon Griff.) in China. A better understanding of the genetic structure of O. rufipogon populations will therefore be useful for the development of conservation strategies. We examined the diversity and genetic structure of natural O. rufipogon populations at the national, provincial, and local levels using simple sequence repeat (SSR) markers. Twenty representative populations from sites across China showed high levels of genetic variability, and approximately 44% of the total genetic variation was among populations. At the local level, we studied fourteen populations in Guangxi Province and four populations in Jiangxi Province. Populations from similar ecosystems showed less genetic differentiation, and local environmental conditions rather than geographic distance appeared to have influenced gene flow during population genetic evolution. We identified a triangular area, including northern Hainan, southern Guangdong, and southwestern Guangxi, as the genetic diversity center of O. rufipogon in China, and we proposed that this area should be given priority during the development of ex situ and in situ conservation strategies. Populations from less common ecosystem types should also be given priority for in situ conservation.