Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

Transposable elements (TEs) significantly drive dynamic changes that characterize genome evolution. However, understanding the variability associated with TE insertions among different cultivars remains challenging. The pomegranate (Punica granatum L.) has yet to be extensively studied regarding the roles of TEs in the diversification of cultivars. Herein, we explored the genome distribution of TEs and its potential functional implications among four pomegranate cultivars, ‘Bhagwa’, ‘Dabenzi’, ‘Taishanhong’ and ‘Tunisia’, whose genome sequences are available. A total of 8404 full-length TEs were isolated. The content of TEs varied among the cultivars, ranging from 41.67% of ‘Taishanhong’ to 52.45% of ‘Bhagwa’. In all cultivars, the Gypsy superfamily of retrotransposons accounted for a larger genome proportion than the Copia superfamily. Seventy-three full-length TEs were found at the same genomic loci in all four cultivars. By contrast, 947, 297, 311, and 874 TEs were found exclusively in ‘Bhagwa’, ‘Dabenzi’, ‘Taishanhong’, and ‘Tunisia’ cultivars, respectively. Phylogenetic clustering based on the presence of TE insertions in specific loci reflected the geographic origins of the cultivars. The insertion time profiles of LTR-REs were studied in the four cultivars. Shared elements across the four cultivars exhibited, on average, a more ancient insertion date than those exclusive to three, two, or one cultivars. The majority of TEs were located within 1000 bp from the nearest gene. This localization was observed for 57% of DNA TEs and 55% of long-terminal repeat retrotransposons (LTR-RE). More than 10% of TEs resulted inserted within genes. Concerning DNA TEs, 3.91% of insertions occurred in introns, while 2.42% occurred in exons. As to LTR-REs, 4% of insertions occurred in exons and 1.98% in introns. Functional analysis of the genes lying close to TEs was performed to infer if differences in TE insertion can affect the fruit quality. Two TE insertions were found close to two genes encoding 4-coumarate--CoA ligase, an enzyme involved in the phenylpropanoid pathway. Moreover, a TIR/Mariner element was found within the exon of a gene encoding anthocyanidin reductase in the ‘Tunisia’ genotype, crucial in the biosynthesis of flavan-3-ols and proanthocyanidins, strictly correlated with the nutraceutical properties of pomegranate. Although functional and metabolomic studies are essential to elucidate the consequences of TE insertions, these results contribute to advancing our comprehension of the role of TEs in pomegranate genomics, providing insights for crop breeding.

Ascorbic acid (AsA) is a key antioxidant and nutrient in plants, regulating reactive oxygen species (ROS) levels and maintaining cellular redox homeostasis. The AsA recycling pathway sustains AsA pools by restoring its oxidized forms, ensuring intracellular balance. Among the enzymes involved, monodehydroascorbate reductase (MDHAR) is important for the regeneration of AsA from monodehydroascorbate. In this study, we analyzed the four MDHAR paralogs in Lactuca sativa using CRISPR/Cas9 to determine whether disruption of individual MDHAR genes could alter AsA levels in lettuce leaves. Unexpectedly, none of the knockouts caused long-term changes in leaf AsA content. Transcriptomic analyses at 14 and 28 days showed minimal effects on AsA recycling or biosynthesis genes, except MDHAR genes. However, several other genes indirectly implicated in AsA regulation displayed differential expression in all mutants compared to the wild type, suggesting the presence of a complex regulatory network. In particular, genes encoding transcription factors (TFs), such as mTERF15, COL9, UPBEAT1, NAC28, and NAC42, were differentially regulated in all MDHAR mutants compared to the wild type at 28 days. These findings indicate that, although AsA content remains unchanged, MDHAR single knockouts alter expression of other genes through which the plants may indirectly compensate to maintain redox homeostasis.

The rapid development of high-throughput sequencing technologies and the integration of omics sciences have markedly advanced plant genomics over the past decade