Alternative splicing triggered by the insertion of a CACTA transposon attenuates LsGLK and leads to the development of pale-green leaves in lettuce

SWATH-MS based proteomics reveals the role of photosynthesis related proteins and secondary metabolic pathways in the colored leaves of sweet olive (Osmanthus fragrans)

Colored leaves, a notable horticultural trait, have high research and ornamental value. The evergreen sweet olive (Osmanthus fragrans), one of the top ten traditional flowers in China, has been cultivated for more than two thousand years. However, in recent years, an increasing number of O. fragrans cultivars with colored leaves have been cultivated for their ornamental value. To study the molecular mechanism underlying the observed changes in leaf color, we selected O. fragrans 'Yinbi Shuanghui' (Y), which has yellow-white leaves, and O. fragrans 'Sijigui' (S), which has green leaves, as materials. Pigment content measurement showed that the chlorophyll, carotenoid and anthocyanin contents in Y were lower than in S. According to the SWATH-MS sequencing results, a total of 3,959 proteins were quantitatively identified, 1,300 of which were differentially expressed proteins (DEPs), including 782 up-regulated and 518 down-regulated proteins in Y compared to S. Functional enrichment analysis of DEPs revealed that down-regulated expression of photosynthesis related proteins may lead to the inhibition of chlorophyll synthesis in Y, this may be the main cause of leaf color change. Moreover, a protein interaction prediction model also showed that proteins such as PetC, PsbO, PsbP, and PsbQ were key proteins in the interaction network, and the up-regulated proteins participating in the anthocyanin and carotenoid pathways may be related to the formation of yellow-white leaves. Taken together, our findings represent the first SWATH-MS-based proteomic report on colored leaf O. fragrans and reveal that chlorophyll synthesis and secondary metabolism pathways contribute to the changes in leaf color.

LsRTDv1, a reference transcript dataset for accurate transcript-specific expression analysis in lettuce

Accurate quantification of gene and transcript-specific expression, with the underlying knowledge of precise transcript isoforms, is crucial to understanding many biological processes. Analysis of RNA sequencing data has benefited from the development of alignment-free algorithms which enhance the precision and speed of expression analysis. However, such algorithms require a reference transcriptome. Here we generate a reference transcript dataset (LsRTDv1) for lettuce (cv. Saladin), combining long- and short-read sequencing with publicly available transcriptome annotations, and filtering to keep only transcripts with high-confidence splice junctions and transcriptional start and end sites. LsRTDv1 identifies novel genes (mostly long non-coding RNAs) and increases the number of transcript isoforms per gene in the lettuce genome from 1.4 to 2.7. We show that LsRTDv1 significantly increases the mapping rate of RNA-seq data from a lettuce time-series experiment (mock- and Botrytis cinerea-inoculated) and enables detection of genes that are differentially alternatively spliced in response to infection as well as transcript-specific expression changes. LsRTDv1 is a valuable resource for investigation of transcriptional and alternative splicing regulation in lettuce.

The pleiotropic functions of GOLDEN2-LIKE transcription factors in plants

The regulation of gene expression is crucial for biological plant growth and development, with transcription factors (TFs) serving as key switches in this regulatory mechanism. GOLDEN2-LIKE (GLK) TFs are a class of functionally partially redundant nuclear TFs belonging to the GARP superfamily of MYB TFs that play a key role in regulating genes related to photosynthesis and chloroplast biogenesis. Here, we summarized the current knowledge of the pleiotropic roles of GLKs in plants. In addition to their primary functions of controlling chloroplast biogenesis and function maintenance, GLKs have been proven to regulate the photomorphogenesis of seedlings, metabolite synthesis, flowering time, leaf senescence, and response to biotic and abiotic stress, ultimately contributing to crop yield. This review will provide a comprehensive understanding of the biological functions of GLKs and serve as a reference for future theoretical and applied studies of GLKs.

PE (Prickly Eggplant) encoding a cytokinin-activating enzyme responsible for the formation of prickles in eggplant

Eggplant is one of the most important vegetables worldwide, with some varieties displaying prickles. These prickles, present on the leaves, stems, and fruit calyxes, posing challenges during cultivation, harvesting, and transportation, making them an undesirable agronomic trait. However, the genetic mechanisms underlying prickle morphogenesis in eggplant remain poorly understood, impeding genetic improvements. In this study, genetic analyses revealed that prickle morphogenesis is governed by a single dominant nuclear gene, termed PE (Prickly Eggplant). Subsequent bulk segregant RNA-sequencing (BSR-seq) and linkage analysis preliminarily mapped PE to chromosome 6. This locus was then fine mapped to a 9233 bp interval in a segregating population of 1109 plants, harboring only one candidate gene, SmLOG1, which encodes a LONELY GUY (LOG)-family cytokinin biosynthetic enzyme. Expression analyses via transcriptome and qRT-PCR demonstrate that SmLOG1 is predominantly expressed in immature prickles. CRISPR-Cas9 knockout experiments targeting SmLOG1 in prickly parental line 'PI 381159' abolished prickles across all tissues, confirming its critical role in prickle morphogenesis. Sequence analysis of SmLOG1 pinpointed variations solely within the non-coding region. We developed a cleaved amplified polymorphic sequences (CAPS) marker from a distinct SNP located at -735-bp within the SmLOG1 promoter, finding significant association with prickle variation in 190 eggplant germplasms. These findings enhance our understanding of the molecular mechanisms governing prickle development in eggplant and facilitate the use of marker-assisted selection (MAS) for breeding prickleless cultivars.

Deciphering transcriptomic signatures explaining the phenotypic plasticity of nonheading lettuce genotypes under artificial light conditions

Elucidating the mechanisms and pathways involved in genotype-environment (G×E) interactions and phenotypic plasticity is critical for improving plant growth. Controlled environment agricultural systems allow growers to modulate the environment for particular genotypes. In this study, we evaluated the effects of interactions among 14 genotypes and four artificial light environments on leaf lettuce phenotypes and dissected the underlying molecular mechanism via transcriptome-based modeling. Variations in morphological traits and phytochemical concentrations in response to artificial light treatments revealed significant G×E interactions. The appropriate genotype and artificial light combinations for maximizing phenotypic expression were determined on the basis of a joint regression analysis and the additive main effect and multiplicative interaction model for these G×E interactions. Transcriptome-based regression modeling explained approximately 50%-90% of the G×E variations. Further analyzes indicated Red Lettuce Leaves 4 (RLL4) regulates UV-B and blue light signaling through the effects of the HY5-MBW pathway on flavonoid biosynthesis and contributes to natural variations in the light-responsive plasticity of lettuce traits. Our study represents an important step toward elucidating the phenotypic variations due to G×E interactions in nonheading lettuce under artificial light conditions.

An insertion of transposon in DcNAP inverted its function in the ethylene pathway to delay petal senescence in carnation (Dianthus caryophyllus L.)

Petal senescence is the final stage of flower development. Transcriptional regulation plays key roles in this process. However, whether and how post-transcriptional regulation involved is still largely unknown. Here, we identified an ethylene-induced NAC family transcription factor DcNAP in carnation (Dianthus caryophyllus L.). One allele, DcNAP-dTdic1, has an insertion of a dTdic1 transposon in its second exon. The dTdic1 transposon disrupts the structure of DcNAP and causes alternative splicing, which transcribes multiple domain-deleted variants (DcNAP2 and others). Conversely, the wild type allele DcNAP transcribes DcNAP1 encoding an intact NAC domain. Silencing DcNAP1 delays and overexpressing DcNAP1 accelerates petal senescence in carnation, while silencing and overexpressing DcNAP2 have the opposite effects, respectively. Further, DcNAP2 could interact with DcNAP1 and interfere the binding and activation activity of DcNAP1 to the promoters of its downstream target ethylene biosynthesis genes DcACS1 and DcACO1. Lastly, ethylene signalling core transcriptional factor DcEIL3-1 can activate the expression of DcNAP1 and DcNAP2 in the same way by binding their promoters. In summary, we discovered a novel mechanism by which DcNAP regulates carnation petal senescence at the post-transcriptional level. It may also provide a useful strategy to manipulate the NAC domains of NAC transcription factors for crop genetic improvement.

Post-transcriptional and translational regulation of plant gene expression by transposons

Transposons are mobile DNA sequences that can move within the genome and integrate in new genomic locations. They are widespread in eukaryotes and prokaryotes and can influence gene expression when landing within or nearby a gene. Although transposon-induced regulation of gene expression at the transcriptional level has been extensively studied, there has been less focus on regulation at the post-transcriptional and translational levels. Recent studies in maize (Zea mays) and other plant species suggest that transposon insertions can affect RNA processing, RNA stability, protein translation and protein stability. We will describe the diverse mechanisms by which transposons can influence gene expression at the post-transcriptional and translational levels, and discuss the interactions between these mechanisms.

Advances in cis-element- and natural variation-mediated transcriptional regulation and applications in gene editing of major crops

Transcriptional regulation is crucial to control of gene expression. Both spatio-temporal expression patterns and expression levels of genes are determined by the interaction between cis-acting elements and trans-acting factors. Numerous studies have focused on the trans-acting factors that mediate transcriptional regulatory networks. However, cis-acting elements, such as enhancers, silencers, transposons, and natural variations in the genome, are also vital for gene expression regulation and could be utilized by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9)-mediated gene editing to improve crop quality and yield. In this review, we discuss current understanding of cis-element-mediated transcriptional regulation in major crops, including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays), as well as the latest advancements in gene editing techniques and their applications in crops to highlight prospective strategies for crop breeding.

RsGLK2.1-RsNF-YA9a module positively regulates the chlorophyll biosynthesis by activating RsHEMA2 in green taproot of radish

Radish (Raphanus sativus L.) is an economically important and widely cultivated root vegetable crop. The coloration of the green skin and green flesh is an important trait influencing the nutrition and flavor quality in fruit radish. GOLDEN2-LIKEs (GLKs) play critically important roles in plastid development and chlorophyll biosynthesis in plants. However, the molecular mechanism underlying chlorophyll biosynthesis still remain elusive in green fruit radish taproot. Herein, the RsGLK2.1 gene exhibited higher expression level in taproot with a green skin (GS) and green flesh (GF) than that in taproot of the white or red radish genotypes. RsGLK2.1 is a nuclear transcription factor that has intrinsic transcriptional activation activity. Overexpression of RsGLK2.1 increased the total chlorophyll content of 20.68%-45.84% in radish leaves. Knockout of the RsGLK2.1 gene via CRISPR/Cas9 technology resulted in a significant decrease in the chlorophyll content. Overexpression of the RsGLK2.1 gene could restore the phenotype of the glk1glk2 mutant Arabidopsis. RsGLK2.1 was participated in regulating the chlorophyll biosynthesis by directly binding to the promoter of RsHEMA2 and activating its transcription. The interaction of RsNF-YA9a with RsGLK2.1 increased the transcriptional activity of the downstream gene RsHEMA2 under the light condition rather than the dark condition, indicating that both of them regulate the chlorophyll biosynthesis in a light-dependent manner of radish. Overall, these results provided insights into the molecular framework of the RsGLK2.1-RsNF-YA9a module, and could facilitate dissecting the regulatory mechanism underlying chlorophyll biosynthesis in green taproot of radish, and genetic improvement of quality traits in fruit radish breeding programs.

The genetic architecture of prolificacy in maize revealed by association mapping and bulk segregant analysis

KEY MESSAGE

Here, we revealed maize prolificacy highly correlated with domestication and identified a causal gene ZmEN1 located in one novel QTL qGEN261 that regulating maize prolificacy by using multiple-mapping methods. The development of maize prolificacy (EN) is crucial for enhancing yield and breeding specialty varieties. To achieve this goal, we employed a genome-wide association study (GWAS) to analyze the genetic architecture of EN in maize. Using 492 inbred lines with a wide range of EN variability, our results demonstrated significant differences in genetic, environmental, and interaction effects. The broad-sense heritability (H2) of EN was 0.60. Through GWAS, we identified 527 significant single nucleotide polymorphisms (SNPs), involved 290 quantitative trait loci (QTL) and 806 genes. Of these SNPs, 18 and 509 were classified as major effect loci and minor loci, respectively. In addition, we performed a bulk segregant analysis (BSA) in an F2 population constructed by a few-ears line Zheng58 and a multi-ears line 647. Our BSA results identified one significant QTL, qBEN1. Importantly, combining the GWAS and BSA, four co-located QTL, involving six genes, were identified. Three of them were expressed in vegetative meristem, shoot tip, internode and tip of ear primordium, with ZmEN1, encodes an unknown auxin-like protein, having the highest expression level in these tissues. It suggested that ZmEN1 plays a crucial role in promoting axillary bud and tillering to encourage the formation of prolificacy. Haplotype analysis of ZmEN1 revealed significant differences between different haplotypes, with inbred lines carrying hap6 having more EN. Overall, this is the first report about using GWAS and BSA to dissect the genetic architecture of EN in maize, which can be valuable for breeding specialty maize varieties and improving maize yield.

Genetic determinants of lettuce resistance to drop caused by Sclerotinia minor identified through genome-wide association mapping frequently co-locate with loci regulating anthocyanin content

KEY MESSAGE

GWAS identified 19 QTLs for resistance to Sclerotinia minor, 11 of them co-locating with red leaf color. Lower disease incidence was observed in red and dark red accessions. Lettuce (Lactuca sativa L.), one of the most economically important vegetables grown primarily in moderate climates around the world, is susceptible to many diseases including lettuce drop caused by the soilborne fungus Sclerotinia minor. Complete resistance to S. minor has not been identified in cultivated lettuce or its wild relatives. We conducted five experiments over 4 years with the diversity panel of almost 500 lettuce accessions to evaluate their response to the pathogen in an artificially infested field. The lowest disease incidence (DI) was observed in cultivars Eruption, Infantry, and Annapolis (median DI of 12.1-17.5%), while the highest DI was recorded for cultivars Reine des Glaces, Wayahead, and line FL. 43007 (median DI of 81.0-95.2%). Overall, significantly lower DI was observed in red and dark red accessions compared to those with a lower anthocyanin content. Genome-wide association mapping identified 19 QTLs for resistance to S. minor, 21 for the presence of red leaf color or its variations caused by the anthocyanin content, and one for the green color intensity. Eleven of the QTLs for disease resistance were located within 10 Mb of the loci associated with red color or anthocyanin content identified in this diversity panel. The frequent, non-random co-location of QTLs, together with the lower DI observed in red and dark red accessions suggests that lettuce interaction with S. minor may be partly influenced by anthocyanins. We have identified RLL2 and ANS, the genes of the anthocyanin biosynthesis pathway that co-locate with resistance QTLs, as candidates for functional studies to ascertain the involvement of anthocyanins in lettuce resistance against S. minor. Resistance QTLs closely linked with QTLs for anthocyanin content could be used to develop lettuce with a relatively high partial resistance and red color, while those not associated with anthocyanins could be used to develop partially resistant cultivars of green color.

LettuceGDB: The community database for lettuce genetics and omics

As a globally popular leafy vegetable and a representative plant of the Asteraceae family, lettuce has great economic and academic significance. In the last decade, high-throughput sequencing, phenotyping, and other multi-omics data in lettuce have accumulated on a large scale, thus increasing the demand for an integrative lettuce database. Here, we report the establishment of a comprehensive lettuce database, LettuceGDB (https://www.lettucegdb.com/). As an omics data hub, the current LettuceGDB includes two reference genomes with detailed annotations; re-sequencing data from over 1000 lettuce varieties; a collection of more than 1300 worldwide germplasms and millions of accompanying phenotypic records obtained with manual and cutting-edge phenomics technologies; re-analyses of 256 RNA sequencing datasets; a complete miRNAome; extensive metabolite information for representative varieties and wild relatives; epigenetic data on the genome-wide chromatin accessibility landscape; and various lettuce research papers published in the last decade. Five hierarchically accessible functions (Genome, Genotype, Germplasm, Phenotype, and O-Omics) have been developed with a user-friendly interface to enable convenient data access. Eight built-in tools (Assembly Converter, Search Gene, BLAST, JBrowse, Primer Design, Gene Annotation, Tissue Expression, Literature, and Data) are available for data downloading and browsing, functional gene exploration, and experimental practice. A community forum is also available for information sharing, and a summary of current research progress on different aspects of lettuce is included. We believe that LettuceGDB can be a comprehensive functional database amenable to data mining and database-driven exploration, useful for both scientific research and lettuce breeding.

LsNRL4 enhances photosynthesis and decreases leaf angles in lettuce

Lettuce (Lactuca sativa) is one of the most important vegetables worldwide and an ideal plant for producing protein drugs. Both well-functioning chloroplasts that perform robust photosynthesis and small leaf angles that enable dense planting are essential for high yields. In this study, we used an F2 population derived from a cross between a lettuce cultivar with pale-green leaves and large leaf angles to a cultivar with dark-green leaves and small leaf angles to clone LsNRL4, which encodes an NPH3/RPT2-Like (NRL) protein. Unlike other NRL proteins in lettuce, the LsNRL4 lacks the BTB domain. Knockout mutants engineered using CRISPR/Cas9 and transgenic lines overexpressing LsNRL4 verified that LsNRL4 contributes to chloroplast development, photosynthesis and leaf angle. The LsNRL4 gene was not present in the parent with pale-green leaves and enlarged leaf angles. Loss of LsNRL4 results in the enlargement of chloroplasts, decreases in the amount of cellular space allocated to chloroplasts and defects in secondary cell wall biosynthesis in lamina joints. Overexpressing LsNRL4 significantly improved photosynthesis and decreased leaf angles. Indeed, the plant architecture of the overexpressing lines is ideal for dense planting. In summary, we identified a novel NRL gene that enhances photosynthesis and influences plant architecture. Our study provides new approaches for the breeding of lettuce that can be grown in dense planting in the open field or in modern plant factories. LsNRL4 homologues may also be used in other crops to increase photosynthesis and improve plant architecture.

Combination analysis of single-molecule long-read and Illumina sequencing provides insights into the anthocyanin accumulation mechanism in an ornamental grass, Pennisetum setaceum cv. Rubrum

Combination analysis of single-molecule long-read and Illumina sequencing provide full-length transcriptome information and shed new light on the anthocyanin accumulation mechanism of Pennisetum setaceum cv. 'Rubrum'. Pennisetum setaceum cv. 'Rubrum' is an ornamental grass with purple leaves widely used in landscaping. However, the current next-generation sequencing (NGS) transcriptome information of this species is not satisfactory due to the difficulties in obtaining full-length transcripts. Furthermore, the molecular mechanisms of anthocyanin accumulation in P. setaceum have not been thoroughly studied. In this study, we used PacBio full-length transcriptome sequencing (SMRT) combined with NGS technology to build and improve the transcriptomic datasets and reveal the molecular mechanism of anthocyanin accumulation in P. setaceum cv. 'Rubrum'. Therefore, 280,413 full-length non-chimeric reads sequences were obtained using the SMRT technology. We obtained 97,450 high-quality non-redundant transcripts and identified 5352 alternative splicing events. In addition, 93,066 open reading frames (ORFs), including 57,457 full ORFs and 2910 long non-coding RNA (lncRNAs) were screened out. Furthermore, 10,795 differentially expressed genes were identified using NGS. We also explored key genes, synthesis pathways, and detected lncRNA involved in anthocyanin accumulation, providing new insights into anthocyanin accumulation in P. setaceum cv. 'Rubrum'. To our best knowledge, we provided the full-length transcriptome information of P. setaceum cv. 'Rubrum' for the first time. The results of this study will provide baseline information for gene function studies and pave the way for future P. setaceum cv. 'Rubrum' breeding projects.