Adaptive evolution and functional innovation of Populus-specific recently evolved microRNAs

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

Profiling of Known and Novel microRNAs in an Oleaginous Crop Native to the Amazon Basin, Sacha Inchi (Plukenetia volubilis), Through smallRNA-Seq

BACKGROUND

MicroRNAs (miRNAs) play crucial roles in regulating tissue-specific gene expression and plant development. This study explores the identification and functional characterization of miRNAs in Plukenetia volubilis (sacha inchi), an economically and nutritionally significant crop native to the Amazon basin, across three organs: root, stem, and leaf.

METHODS

Small RNA libraries were sequenced on the Illumina Novaseq 6000 platform, yielding high-quality reads that facilitated the discovery of known and novel miRNAs using miRDeep-P.

RESULTS

A total of 277 miRNAs were identified, comprising 71 conserved and 206 novel miRNAs, across root, stem, and leaf tissues. In addition, differential expression analysis using DESeq2 identified distinct miRNAs exhibiting tissue-specific regulation. Notably, novel miRNAs like novel_1, novel_88, and novel_189 showed significant roles in processes such as auxin signaling, lignin biosynthesis, and stress response. Functional enrichment analysis of miRNA target genes revealed pathways related to hormonal regulation, structural reinforcement, and environmental adaptation, highlighting tissue-specific functions. The Principal Component Analysis and PERMANOVA confirmed clear segregation of miRNA expression profiles among tissues, underlining organ-specific regulation. Differential expression patterns emphasized unique regulatory roles in each organ: roots prioritized stress response and nutrient uptake, leaves focused on photosynthesis and UV protection, and stems contributed to structural integrity and nutrient transport, suggesting evolutionary adaptations in P. volubilis.

CONCLUSIONS

This study identified novel miRNA-mediated networks that regulate developmental and adaptive processes in P. volubilis, underscoring its molecular adaptations for resilience and productivity. By characterizing both conserved and novel miRNAs, the findings lay a foundation for genetic improvement and molecular breeding strategies aimed at enhancing agronomic traits, stress tolerance, and the production of bioactive compounds.

The miR6445-NAC029 module regulates drought tolerance by regulating the expression of glutathione S-transferase U23 and reactive oxygen species scavenging in Populus

MicroRNAs are essential in plant development and stress resistance, but their specific roles in drought stress require further investigation. Here, we have uncovered that a Populus-specific microRNAs (miRNA), miR6445, targeting NAC (NAM, ATAF, and CUC) family genes, is involved in regulating drought tolerance of poplar. The expression level of miR6445 was significantly upregulated under drought stress; concomitantly, seven targeted NAC genes showed significant downregulation. Silencing the expression of miR6445 by short tandem target mimic technology significantly decreased the drought tolerance in poplar. Furthermore, 5' RACE experiments confirmed that miR6445 directly targeted NAC029. The overexpression lines of PtrNAC029 (OE-NAC029) showed increased sensitivity to drought compared with knockout lines (Crispr-NAC029), consistent with the drought-sensitive phenotype observed in miR6445-silenced strains. PtrNAC029 was further verified to directly bind to the promoters of glutathione S-transferase U23 (GSTU23) and inhibit its expression. Both Crispr-NAC029 and PtrGSTU23 overexpressing plants showed higher levels of PtrGSTU23 transcript and GST activity while accumulating less reactive oxygen species (ROS). Moreover, poplars overexpressing GSTU23 demonstrated enhanced drought tolerance. Taken together, our research reveals the crucial role of the miR6445-NAC029-GSTU23 module in enhancing poplar drought tolerance by regulating ROS homeostasis. This finding provides new molecular targets for improving the drought resistance of trees.

Biogenesis, Mode of Action and the Interactions of Plant Non-Coding RNAs

The central dogma of genetics, which outlines the flow of genetic information from DNA to RNA to protein, has long been the guiding principle in molecular biology. In fact, more than three-quarters of the RNAs produced by transcription of the plant genome are not translated into proteins, and these RNAs directly serve as non-coding RNAs in the regulation of plant life activities at the molecular level. The breakthroughs in high-throughput transcriptome sequencing technology and the establishment and improvement of non-coding RNA experiments have now led to the discovery and confirmation of the biogenesis, mechanisms, and synergistic effects of non-coding RNAs. These non-coding RNAs are now predicted to play important roles in the regulation of gene expression and responses to stress and evolution. In this review, we focus on the synthesis, and mechanisms of non-coding RNAs, and we discuss their impact on gene regulation in plants.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

Coordinating Diverse Functions of miRNA and lncRNA in Fleshy Fruit

Non-coding RNAs play vital roles in the diverse biological processes of plants, and they are becoming key topics in horticulture research. In particular, miRNAs and long non-coding RNAs (lncRNAs) are receiving increased attention in fruit crops. Recent studies in horticulture research provide both genetic and molecular evidence that miRNAs and lncRNAs regulate biological function and stress responses during fruit development. Here, we summarize multiple regulatory modules of miRNAs and lncRNAs and their biological roles in fruit sets and stress responses, which would guide the development of molecular breeding techniques on horticultural crops.

Ascorbate-Glutathione Cycle Genes Families in Euphorbiaceae: Characterization and Evolutionary Analysis

Ascorbate peroxidase (APX), Monodehydroascorbate Reductase (MDAR), Dehydroascorbate Reductase (DHAR) and Glutathione Reductase (GR) enzymes participate in the ascorbate-glutathione cycle, which exerts a central role in the antioxidant metabolism in plants. Despite the importance of this antioxidant system in different signal transduction networks related to development and response to environmental stresses, the pathway has not yet been comprehensively characterized in many crop plants. Among different eudicotyledons, the Euphorbiaceae family is particularly diverse with some species highly tolerant to drought. Here the APX, MDAR, DHAR, and GR genes in Ricinus communis, Jatropha curcas, Manihot esculenta, and Hevea brasiliensis were identified and characterized. The comprehensive phylogenetic and genomic analyses allowed the classification of the genes into different classes, equivalent to cytosolic, peroxisomal, chloroplastic, and mitochondrial enzymes, and revealed the duplication events that contribute to the expansion of these families within plant genomes. Due to the high drought stress tolerance of Ricinus communis, the expression patterns of ascorbate-glutathione cycle genes in response to drought were also analyzed in leaves and roots, indicating a differential expression during the stress. Altogether, these data contributed to the characterization of the expression pattern and evolutionary analysis of these genes, filling the gap in the proposed functions of core components of the antioxidant mechanism during stress response in an economically relevant group of plants.

Identification and Comparative Analysis of Conserved and Species-Specific microRNAs in Four Populus Sections

The conservation and diversity of microRNA (miRNA) families provide insights into the evolution of miRNA genes. However, there are few studies to explore the miRNA genes at the genus level in plants. Here, we identified 1194 miRNA loci in four Populus species P. deltoides, P. euphratica, P. tremula, and P. trichocarpa from Aigeiros, Turanga, Populus, and Tacamahaca sections, respectively, by combining de novo and homolog-based approaches. Our results indicated that a similar number of miRNA loci exist in each species (296–301 miRNA loci). Among the identified 143 miRNA families, 68 families are shared by the studied four species, and 31 families are species-specific, which might be related to local adaptation. Additionally, multiple miRNA-related single nucleotide polymorphisms (SNPs) were found, indicating that polymorphisms in pre-miRNA hairpins were likely to affect miRNA biogenesis. This study expanded the breadth and depth of miRNA annotations and provided valuable resources for further exploring the diversity and function of poplar miRNAs.

Association Study and Mendelian Randomization Analysis Reveal Effects of the Genetic Interaction Between PtoMIR403b and PtoGT31B-1 on Wood Formation in Populus tomentosa

MicroRNAs (miRNAs), important posttranscriptional regulators of gene expression, play a crucial role in plant growth and development. A single miRNA can regulate numerous target genes, making the determination of its function and interaction with targets challenging. We identified PtomiR403b target to PtoGT31B-1, which encodes a galactosyltransferase responsible for the biosynthesis of cell wall polysaccharides. We performed an association study and epistasis and Mendelian randomization (MR) analyses to explore how the genetic interaction between PtoMIR403b and its target PtoGT31B-1 underlies wood formation. Single nucleotide polymorphism (SNP)-based association studies identified 25 significant associations (P < 0.01, Q < 0.05), and PtoMIR403b and PtoGT31B-1 were associated with five traits, suggesting a role for PtomiR403b and PtoGT31B-1 in wood formation. Epistasis analysis identified 93 significant pairwise epistatic associations with 10 wood formation traits, and 37.89% of the SNP-SNP pairs indicated interactions between PtoMIR403b and PtoGT31B-1. We performed an MR analysis to demonstrate the causality of the relationships between SNPs in PtoMIR403b and wood property traits and that PtoMIR403b modulates wood formation by regulating expression of PtoGT31B-1. Therefore, our findings will facilitate dissection of the functions and interactions with miRNA-targets.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Synonymous mutation in Growth Regulating Factor 15 of miR396a target sites enhances photosynthetic efficiency and heat tolerance in poplar

Heat stress damages plant tissues and induces multiple adaptive responses. Complex and spatiotemporally specific interactions among transcription factors (TFs), microRNAs (miRNAs), and their targets play crucial roles in regulating stress responses. To explore these interactions and to identify regulatory networks in perennial woody plants subjected to heat stress, we integrated time-course RNA-seq, small RNA-seq, degradome sequencing, weighted gene correlation network analysis, and multi-gene association approaches in poplar. Results from Populus trichocarpa enabled us to construct a three-layer, highly interwoven regulatory network involving 15 TFs, 45 miRNAs, and 77 photosynthetic genes. Candidate gene association studies in a population of P. tomentosa identified 114 significant associations and 696 epistatic SNP–SNP pairs that were linked to 29 photosynthetic and growth traits (P&lt;0.0001, q&lt;0.05). We also identified miR396a and its target, Growth-Regulating Factor 15 (GRF15) as an important regulatory module in the heat-stress response. Transgenic plants of hybrid poplar (P. alba × P. glandulosa) overexpressing a GRF15 mRNA lacking the miR396a target sites exhibited enhanced heat tolerance and photosynthetic efficiency compared to wild-type plants. Together, our observations demonstrate that GRF15 plays a crucial role in responding to heat stress, and they highlight the power of this new, multifaceted approach for identifying regulatory nodes in plants.

Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene

Background

Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays.

Results

A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3’UTR at a unique binding site.

Conclusions

Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07536-y.

PhasiRNAs in Plants: Their Biogenesis, Genic Sources, and Roles in Stress Responses, Development, and Reproduction

Phased secondary small interfering RNAs (phasiRNAs) constitute a major category of small RNAs in plants, but most of their functions are still poorly defined. Some phasiRNAs, known as trans-acting siRNAs, are known to target complementary mRNAs for degradation and to function in development. However, the targets or biological roles of other phasiRNAs remain speculative. New insights into phasiRNA biogenesis, their conservation, and their variation across the flowering plants continue to emerge due to the increased availability of plant genomic sequences, deeper and more sophisticated sequencing approaches, and improvements in computational biology and biochemical/molecular/genetic analyses. In this review, we survey recent progress in phasiRNA biology, with a particular focus on two classes associated with male reproduction: 21-nucleotide (accumulate early in anther ontogeny) and 24-nucloetide (produced in somatic cells during meiosis) phasiRNAs. We describe phasiRNA biogenesis, function, and evolution and define the unanswered questions that represent topics for future research.

Evolution and Functional Divergence of the Fructokinase Gene Family in Populus

New kinase has emerged throughout evolution, but how new kinase evolve while maintaining their functions and acquiring new functions remains unclear. Fructokinase (FRK), the gateway kinase to fructose metabolism, plays essential roles in plant development, and stress tolerance. Here, we explored the evolution of FRK gene family in 20 plant species (from green algae to angiosperms) and their functional roles in Populus. We identified 125 putative FRK genes in the 20 plant species with an average of 6 members per species. Phylogenetic analysis separated these 125 genes into 8 clades including 3 conserved clades and 5 specific clades, the 5 of which only exist in green algae or angiosperms. Evolutionary analysis revealed that FRK genes in ancient land plants have the largest number of functional domains with the longest amino acid sequences, and the length of FRK genes became shorter during the transition to vascular plants. This was accompanied by loss, acquisition, and diversification of functional domains. In Populus, segmental duplication appears to be the main mechanism for the expansion of FRK genes. Specially, most FRK genes duplicated in salicoids are regulated by Populus-specific microRNAs. Furthermore, compared with common FRKs, Populus-specific FRKs have showed higher expression specificity and are associated with fewer growth and wood property traits, which suggests that these FRKs may have undergone functional divergence. Our study explores the specific roles of FRKs in the Populus genome and provides new insights for functional investigation of this gene family.

Genetic interactions among Pto-miR319 family members and their targets influence growth and wood properties in Populus tomentosa

MicroRNAs (miRNAs) play crucial roles in all aspects of plant growth and development, but the genetic interactions of miRNAs and their target genes in woody plants are largely unknown. Here, we integrated association genetics and expression profiling to decipher the allelic variations and interactions of the Pto-MIR319 family of miRNAs and 12 putative Pto-miR319 target genes related to wood formation in 435 unrelated individuals of Populus tomentosa Carrière (Chinese white poplar). Expression pattern analysis showed that among all pairings between expressions of pre-miRNA of Pto-MIR319 members and targets, 70.0% showed negative correlation of expression levels (r = - 0.944 to 0.674, P < 0.01) in eight tissues and organs of poplar, suggesting that Pto-miR319 may participate in the regulatory network of wood formation. Single SNP-based association studies identified 137 significant associations (P < 0.01, Q < 0.1), representing 126 unique SNPs from Pto-MIR319 members and their targets, with 10 tree growth traits, revealing that these genetic factors have common roles related to wood formation. Epistasis analysis uncovered 105 significant SNP-SNP associations (P < 0.01) influencing the 10 traits, demonstrating the close genetic interactions between Pto-MIR319 family members and the 12 Pto-miR319 target genes. Notably, one common SNP, in the precursor region of Pto-MIR319e, affected the stability of Pto-MIR319e's secondary structure by altering the stem-loop structure and minimum free energy, contributing to variations in the expression of Pto-MIR319e and Pto-miR319e target genes. This study enriches the understanding of the functions of miR319 family miRNAs in poplar and exemplifies a feasible approach to exploring the genetic effects underlying miRNA-mRNA interactions related to complex traits in trees.

Functional Similarity and Difference among Bra-MIR319 Family in Plant Development

miR319 was the first plant miRNA discovered via forward genetic mutation screening. In this study, we found that miR319 family members had similar sequences but different expression patterns in Brassica campestris and Arabidopsis thaliana. RT-PCR analysis revealed that Bra-MIR319a and Bra-MIR319c had similar expression patterns and were widely expressed in plant development, whereas Bra-MIR319b could only be detected in stems. The overexpression of each Bra-MIR319 family member in Arabidopsis could inhibit cell division and function in leaf and petal morphogenesis. Bra-miR319a formed a new regulatory relationship after whole genome triplication, and Bra-MIR319a overexpressing in Arabidopsis led to the degradation of pollen content and affected the formation of intine, thereby causing pollen abortion. Our results suggest that Bra-MIR319 family members have functional similarity and difference in plant development.

Populus trichocarpa clade A PP2C protein phosphatases: their stress-induced expression patterns, interactions in core abscisic acid signaling, and potential for regulation of growth and development

Overexpression of the poplar PP2C protein phosphatase gene PtrHAB2 resulted in increased tree height and altered leaf morphology and phyllotaxy, implicating PP2C phosphatases as growth regulators functioning under favorable conditions. We identified and studied Populus trichocarpa genes, PtrHAB1 through PtrHAB15, belonging to the clade A PP2C family of protein phosphatases known to regulate abscisic acid (ABA) signaling. PtrHAB1 through PtrHAB3 and PtrHAB12 through PtrHAB15 were the most highly expressed genes under non-stress conditions. The poplar PP2C genes were differentially regulated by drought treatments. Expression of PtrHAB1 through PtrHAB3 was unchanged or downregulated in response to drought, while all other PtrHAB genes were weakly to strongly upregulated in response to drought stress treatments. Yeast two-hybrid assays involving seven ABA receptor proteins (PtrRCAR) against 12 PtrHAB proteins detected 51 interactions involving eight PP2Cs and all PtrRCAR proteins with 22 interactions requiring the addition of ABA. PtrHAB2, PtrHAB12, PtrHAB13 and PtrHAB14 also interacted with the sucrose non-fermenting related kinase 2 proteins PtrSnRK2.10 and PtrSnRK2.11, supporting conservation of a SnRK2 signaling cascade regulated by PP2C in poplar. Additionally, PtrHAB2, PtrHAB12, PtrHAB13 and PtrHAB14 interacted with the mitogen-activated protein kinase protein PtrMPK7. Due to its interactions with PtrSnRK2 and PtrMPK7 proteins, and its reduced expression during drought stress, PtrHAB2 was overexpressed in poplar to test its potential as a growth regulator under non-stress conditions. 35S::PtrHAB2 transgenics exhibited increased growth rate for a majority of transgenic events and alterations in leaf phyllotaxy and morphology. These results indicate that PP2Cs have additional roles which extend beyond canonical ABA signaling, possibly coordinating plant growth and development in response to environmental conditions.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

Conserved noncoding sequences conserve biological networks and influence genome evolution

Comparative genomics approaches have identified numerous conserved cis-regulatory sequences near genes in plant genomes. Despite the identification of these conserved noncoding sequences (CNSs), our knowledge of their functional importance and selection remains limited. Here, we used a combination of DNA methylome analysis, microarray expression analyses, and functional annotation to study these sequences in the model tree Populus trichocarpa. Methylation in CG contexts and non-CG contexts was lower in CNSs, particularly CNSs in the 5'-upstream regions of genes, compared with other sites in the genome. We observed that CNSs are enriched in genes with transcription and binding functions, and this also associated with syntenic genes and those from whole-genome duplications, suggesting that cis-regulatory sequences play a key role in genome evolution. We detected a significant positive correlation between CNS number and protein interactions, suggesting that CNSs may have roles in the evolution and maintenance of biological networks. The divergence of CNSs indicates that duplication-degeneration-complementation drives the subfunctionalization of a proportion of duplicated genes from whole-genome duplication. Furthermore, population genomics confirmed that most CNSs are under strong purifying selection and only a small subset of CNSs shows evidence of adaptive evolution. These findings provide a foundation for future studies exploring these key genomic features in the maintenance of biological networks, local adaptation, and transcription.

Association Genetics in Populus Reveal the Allelic Interactions of Pto-MIR167a and Its Targets in Wood Formation

MicroRNAs (miRNAs) play crucial regulatory roles in plant growth and development by interacting with RNA molecules, including messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs); however, the genetic networks of miRNAs and their targets influencing the phenotypes of perennial trees remain to be investigated. Here, we integrated expression profiling and association analysis of underlying physiology and expression traits to dissect the allelic variations and genetic interactions of Pto-MIR167a and its targets, sponge lncRNA ARFRL, and Pto-ARF8, in 435 unrelated individuals of Populus tomentosa. Tissue-specific expression analysis in eight tissues, including stem, leaf, root, and shoot apex, revealed negative correlations between Pto-MIR167a and lncRNA ARFRL and Pto-ARF8 (r = -0.60 and -0.61, respectively, P < 0.01), and a positive correlation between sponge lncRNA ARFRL and Pto-ARF8 (r = 0.90, P < 0.01), indicating their potential regulatory roles in tree growth and wood formation. Single nucleotide polymorphism (SNP)-based association studies detected 53 significant associations (P < 0.01, Q < 0.1) representing 41 unique SNPs from the three genes and six traits, suggesting their potential roles in wood formation. Epistasis uncovered 88 pairwise interactions for 10 traits, which provided substantial evidence for genetic interactions among Pto-MIR167a, lncRNA ARFRL, and Pto-ARF8. Using gene expression-based association mapping, we also examined SNPs within the three genes that influence phenotypes by regulating the expression of Pto-ARF8. Interestingly, SNPs in the precursor region of Pto-MIR167a altered its secondary structure stability and transcription, thereby affecting the expression of its targets. In summary, we elucidated the genetic interactions between Pto-MIR167a and its targets, sponge lncRNA ARFRL, and Pto-ARF8, in tree growth and wood formation, and provide a feasible method for further investigation of multi-factor genetic networks influencing phenotypic variation in the population genetics of trees.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.