Comparative analyses of American and Asian lotus genomes reveal insights into petal color, carpel thermogenesis and domestication

Genome-wide development of simple sequence repeat (SSR) markers at 2-Mb intervals in lotus (Nelumbo Adans.)

BACKGROUND

Despite the rapid advancement of high-throughput sequencing, simple sequence repeats (SSRs) remain indispensable molecular markers for various applied and research tasks owing to their cost-effectiveness and ease of use. However, existing SSR markers cannot meet the growing demand for research on lotus (Nelumbo Adans.) given their scarcity and weak connections to the lotus genome.

METHODS

Using whole-genome resequencing, active SSR loci were identified throughout the genomes of eight typical Asian lotus. After that, high polymorphism SSR molecular markers were mined from each 2n + 0.5 Mb site on each chromosome (e.g., Chr.1-2.5, 4.5, 6.5 Mb) through four steps: online primer design, primer pair evaluation, agarose gel electrophoresis testing using six Asian lotus, one American lotus, and two their hybrids, and DNA sequence alignment. Finally, the polymerase chain reaction (PCR) efficiency of several SSR markers was validated in 20 Asian temperate lotus, eight Asian tropical lotus, and one American lotus.

RESULTS

A total of 463 SSR markers were developed based on each 2n + 0.5 Mb site of the eight lotus chromosomes (totaling 821.29 Mb). These markers were evenly distributed throughout the lotus genome at a density of 1 SSR per 1.76 Mb. The chromosomal locations of the SSR markers were determined precisely, and the specificity of the primer pairs for each site was verified by sequencing the PCR products. We further provided a set of genome-wide SSR loci, covering 129 per Mb, identified from eight representative Asian lotus, allowing other researchers to independently discover specific SSR markers for particular experiments.

CONCLUSION

These SSR markers, which have a density of 1 SSR marker per 1.76 Mb in this study, will act as a bridge connecting lotus phenotypes with the genome. This work reveals a novel and convenient strategy for developing highly polymorphic SSR markers at any location throughout the lotus genome, and it sheds light on the development of SSR molecular markers in other plant species.

A transposon DNA/hAT-Ac insertion promotes the formation of yellow tepals in lotus (Nelumbo)

Yellow tepal is a unique trait of the American lotus (Nelumbo lutea), and all yellow lotus cultivars in the market possess genetic material from the American lotus. However, the formation of yellow tepals in lotus and the genetic mechanism of their formation remain unclear. In this study, we identified a transposon DNA/hAT-Ac, located within the promoter region of an R2R3-MYB transcription factor, MYB12, by comparing the insertion patterns of transposons in the genomes of American and Asian lotus (Nelumbo nucifera). The transposon was found exclusively in yellow lotus cultivars and not in red or white lotus. The insertion of DNA/hAT-Ac facilitated the specific expression of MYB12 in the yellow lotus tepals. Transient expression in lotus tepals, dual-luciferase, and yeast one-hybrid assays demonstrated that MYB12 promotes the accumulation of carotenoids and flavonols by activating the expression of genes involved in carotenoid and flavonols biosynthesis, and it directly binds to the promoters of PSY and FLS. Our results indicated that the transposon DNA/hAT-Ac-mediated specific expression of MYB12 is crucial for the formation of yellow tepals in lotus, and the findings provide a theoretical basis for the breeding of yellow lotus cultivars.

NnSBE1 encodes a starch branching enzyme involved in starch biosynthesis in lotus seeds

Lotus seed starch holds vast potential for utilization across various industries, with its content and structure directly influencing the commercial value of lotus seeds. However, there has been limited information available on the molecular mechanisms underlying lotus seed starch biosynthesis. In this study, three starch branching enzyme homologs were identified in the lotus genome, designated as NnSBE1 to NnSBE3, which possess conserved CBM_48 and α_Aamy domains. Among them, NnSBE1 exhibited predominant expression, with abundant transcript levels observed in lotus seeds and flower-related organs. Expression of NnSBE1 remained consistently up-regulated in lotus cotyledons from 6 to 21 days after pollination. Additionally, a C2H2-type finger protein encoding gene, NnLOL1, co-expressed with NnSBE1 in lotus cotyledons. As a seed-predominantly expressed transcription factor, NnLOL1 was confirmed to activate NnSBE1 expression. Transient overexpression of NnSBE1 in lotus cotyledons resulted in a significant increase in both amylopectin and total starch content compared to the control. Furthermore, multiple variation sites within the NnSBE1 gene gave rise to diverse haplotypes between seed-lotus and other lotus varieties. These findings contribute to our understanding of the regulation mechanisms involved in lotus seed starch biosynthesis, offering valuable theoretical insights for the genetic improvement of lotus seed starch by molecular breeding strategies.

Integration of Mitoflash and Time-Series Transcriptomics Facilitates Energy Dynamics Tracking and Substrate Supply Analysis of Floral Thermogenesis in Lotus

The high biosynthetic and energetic demands of floral thermogenesis render thermogenic plants the ideal systems to characterize energy metabolism in plants, but real-time tracking of energy metabolism in plant cells remains challenging. In this study, a new method was developed for tracking the mitochondrial energy metabolism at the single mitochondria level by real-time imaging of mitochondrial superoxide production (i.e., mitoflash). Using this method, we observed the increased mitoflash frequencies in the receptacles of Nelumbo nucifera Gaertn. at the thermogenic stages. This increase, combined with the higher expression of antioxidant response-related genes identified through time-series transcriptomics at the same stages, shows us a new regulatory mechanism for plant redox balance. Furthermore, we found that the upregulation of respiratory metabolism-related genes during the thermogenic stages not only correlates with changes in mitoflash frequency but also underscores the critical roles of these pathways in ensuring adequate substrate supply for thermogenesis. Metabolite analysis revealed that sugars are likely one of the substrates for thermogenesis and may be transported over long distances by sugar transporters. Taken together, our findings demonstrate that mitoflash is a reliable tool for tracking energy metabolism in thermogenic plants and contributes to our understanding of the regulatory mechanisms underlying floral thermogenesis.

Gene expression and metabolite levels converge in the thermogenic spadix of skunk cabbage

The inflorescence (spadix) of skunk cabbage (Symplocarpus renifolius) is strongly thermogenic and can regulate its temperature at around 23 °C even when the ambient temperature drops below freezing. To elucidate the mechanisms underlying developmentally controlled thermogenesis and thermoregulation in skunk cabbage, we conducted a comprehensive transcriptome and metabolome analysis across 3 developmental stages of spadix development. Our RNA-seq analysis revealed distinct groups of expressed genes, with selenium-binding protein 1/methanethiol oxidase (SBP1/MTO) exhibiting the highest levels in thermogenic florets. Notably, the expression of alternative oxidase (AOX) was consistently high from the prethermogenic stage through the thermogenic stage in the florets. Metabolome analysis showed that alterations in nucleotide levels correspond with the developmentally controlled and tissue-specific thermogenesis of skunk cabbage, evident by a substantial increase in AMP levels in thermogenic florets. Our study also reveals that hydrogen sulfide, a product of SBP1/MTO, inhibits cytochrome c oxidase (COX)-mediated mitochondrial respiration, while AOX-mediated respiration remains relatively unaffected. Specifically, at lower temperatures, the inhibitory effect of hydrogen sulfide on COX-mediated respiration increases, promoting a shift toward the dominance of AOX-mediated respiration. Finally, despite the differential regulation of genes and metabolites throughout spadix development, we observed a convergence of gene expression and metabolite accumulation patterns during thermogenesis. This synchrony may play a key role in developmentally regulated thermogenesis. Moreover, such convergence during the thermogenic stage in the spadix may provide a solid molecular basis for thermoregulation in skunk cabbage.

Evaluation and validation of suitable reference genes for quantitative real-time PCR analysis in lotus (Nelumbo nucifera Gaertn.)

The qRT-PCR technique has been regarded as an important tool for assessing gene expression diversity. Selection of appropriate reference genes is essential for validating deviation and obtaining reliable and accurate results. Lotus (Nelumbo nucifera Gaertn) is a common aquatic plant with important aesthetic, commercial, and cultural values. Twelve candidate genes, which are typically used as reference genes for qRT-PCR in other plants, were selected for this study. These candidate reference genes were cloned with, specific primers designed based on published sequences. In particular, the expression level of each gene was examined in different tissues and growth stages of Lotus. Notably, the expression stability of these candidate genes was assessed using the software programs geNorm and NormFinder. As a result, the most efficient reference genes for rootstock expansion were TBP and UBQ. In addition, TBP and EF-1α were the most efficient reference genes in various floral tissues, while ACT and GAPDH were the most stable genes at all developmental stages of the seed. CYP and GAPDH were the best reference genes at different stages of leaf development, but TUA was the least stable. Meanwhile, the gene expression profile of NnEXPA was analyzed to confirm the validity of the findings. It was concluded that, TBP and GAPDH were identified as the best reference genes. The results of this study may help researchers to select appropriate reference genes and thus obtain credible results for further quantitative RT-qPCR gene expression analyses in Lotus.

NnSUS1 encodes a sucrose synthase involved in sugar accumulation in lotus seed cotyledons

Fresh lotus seeds are gaining favor with consumers for their crunchy texture and natural sweetness. However, the intricacies of sugar accumulation in lotus seeds remain elusive, which greatly hinders the quality improvement of fresh lotus seeds. This study endeavors to elucidate this mechanism by identifying and characterizing the sucrose synthase (SUS) gene family in lotus. Comprising five distinct members, namely NnSUS1 to NnSUS5, each gene within this family features a C-terminal glycosyl transferase1 (GT1) domain. Among them, NnSUS1 is the predominately expressed gene, showing high transcript abundance in the floral organs and cotyledons. NnSUS1 was continuously up-regulated from 6 to 18 days after pollination (DAP) in lotus cotyledons. Furthermore, NnSUS1 demonstrates co-expression relationships with numerous genes involved in starch and sucrose metabolism. To investigate the function of NnSUS1, a transient overexpression system was established in lotus cotyledons, which confirmed the gene's contribution to sugar accumulation. Specifically, transient overexpression of NnSUS1 in seed cotyledons leads to a significant increase in the levels of total soluble sugar, including sucrose and fructose. These findings provide valuable theoretical insights for improving sugar content in lotus seeds through molecular breeding methods.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

FLOWERING LOCUS T genes control floral induction in lotus

The transition to flowering is a vital process in the lotus life cycle that significantly impacts its ornamental value and seed production. However, the molecular basis of floral transition in lotus remains largely unknown. Here, eight homologous FLOWERING LOCUS T (FT) genes were initially characterized in lotus, which were designated as NnFT1-NnFT8. All of these genes were found to possess the conserved PEBP domain and exhibited high transcript levels in both lotus leaves and floral organs. The proNnFT:β-glucuronidase (GUS) assay exhibited GUS staining in the vascular tissues of leaves. Furthermore, subcellular localization revealed that NnFT proteins were present in various cellular organelles, including the nucleus, cytoplasm, and endoplasmic reticulum. Overexpression of two NnFT homologs, NnFT2 and NnFT3, rescued the late flowering phenotype in the Arabidopsis ft-10 mutant, indicating the stimulative roles of NnFTs in floral induction. Moreover, NnFTs demonstrated interactions with a bZIP transcription factor, FLOWERING LOCUS D (NnFD), both in vitro and in vivo. These findings will not only deepen our understanding of the regulatory mechanism underlying lotus floral transition, but also provide valuable genetic resources for creating new lotus varieties with extended blooming periods using molecular strategies in the future.

Exploring domestication pattern in lotus: insights from dispensable genome assembly

Lotus (Nelumbo nucifera Gaertn.), an important aquatic plant in horticulture and ecosystems, has been cultivated for more than 7000 years and domesticated into three different subgroups: flower lotus, rhizome lotus, and seed lotus. To explore the domesticated regions of each subgroup, re-sequencing data of 371 lotus accessions collected from the public database were aligned to the genome of 'China-Antique (CA)'. Unmapped reads were used to build the dispensable genome of each subgroup using a metagenome-like assembly strategy. More than 27 Mb of the dispensable genome in these three subgroups and the wild group was assembled, of which 11,761 genes were annotated. Some of the contigs in the dispensable genome were similar to the genomic segments of other lotus accessions other than 'CA'. The annotated genes in each subgroup played essential roles in specific developmental processes. Dissection of selective signals in three cultivated subgroups also demonstrated that subgroup-specific metabolic pathways, such as the brassinosteroids metabolism enrichment in FL, associated with these selected genes in each subgroup and the contigs in dispensable genome nearly located in the domesticated regions of each subgroup, respectively. Our data presented a valuable resource for facilitating lotus genomic studies, complemented the helpful information to the reference genome, and shed light on the selective signals of domesticated subgroups.

Transcription factor NnMYB5 controls petal color by regulating GLUTATHIONE S-TRANSFERASE2 in Nelumbo nucifera

Lotus (Nelumbo spp.) is an important aquatic ornamental genus in the family Nelumbonaceae comprising only 2 species: Nelumbo lutea with yellow flowers and Nelumbo nucifera with red or white flowers. The petal color variations between these 2 species have previously been associated with the potential activities of FLAVONOL SYNTHASE (FLS) and MYB5. However, the underlying genetic mechanisms of flower color divergence within the N. nucifera species remain unclear. Here, quantitative trait locus mapping led to the identification of MYB5, a candidate gene controlling petal color in N. nucifera. Genotyping of 213 natural lotus accessions revealed an 80 kb presence/absence variant (PAV) of the NnMYB5 gene that is associated with petal color variation. Transcriptome analysis, dual-luciferase, and yeast 1-hybrid assays showed that NnMYB5 could directly activate the anthocyanin transporter gene GLUTATHIONE S-TRANSFERASE2 (NnGST2). Heterologous expression of NnGST2 in Arabidopsis (Arabidopsis thaliana) and its overexpression in lotus petals induced anthocyanin accumulation. Deletion of the 80 kb PAV within NnMYB5 inactivated NnGST2 expression and blocked anthocyanin accumulation in white N. nucifera petals. In contrast, the anthocyanin deficiency of N. lutea occurred due to pseudogenized NlMYB5 alleles. Our results establish a regulatory link between NnMYB5 and NnGST2 in petal anthocyanin accumulation and demonstrate the independent mechanisms controlling flower coloration in Nelumbo.

Analysis of bZIP gene family in lotus (Nelumbo) and functional study of NnbZIP36 in regulating anthocyanin synthesis

BACKGROUND

The basic leucine zipper (bZIP) family is a predominant group of transcription factors in plants, involved in regulating plant growth, development, and response to stressors. Additionally, the bZIP gene family has a key role in anthocyanin production. Despite the significant role of bZIP genes in plants, their potential contribution in lotus remains understudied.

RESULTS

A total of 124 bZIP genes (59 NnbZIPs and 65 NlbZIPs) were identified from genomes of two lotus species. These genes were classified into 13 groups according to the grouping principle of the Arabidopsis bZIP gene family. Analysis of promoter cis-acting elements indicated that most bZIP gene family members in lotus are associated with response to abiotic stresses. The promoters of some bZIP genes contain MYB binding sites that regulate anthocyanin synthesis. We examined the anthocyanin content of the petals from three different colored lotus, combined with transcriptome data analysis and qRT-PCR results, showing that the expression trends of NnbZIP36 and the homologous gene NlbZIP38 were significantly correlated with the anthocyanin content in lotus petals. Furthermore, we found that overexpression of NnbZIP36 in Arabidopsis promoted anthocyanin accumulation by upregulating the expression of genes (4CL, CHI, CHS, F3H, F3'H, DFR, ANS and UF3GT) related to anthocyanin synthesis.

CONCLUSIONS

Our study enhances the understanding of the bZIP gene family in lotus and provides evidence for the role of NnbZIP36 in regulating anthocyanin synthesis. This study also sets the stage for future investigations into the mechanism by which the bZIP gene family regulates anthocyanin biosynthesis in lotus.

Metabolite profiling and screening of callus browning-related genes in lotus (Nelumbo nucifera)

Callus browning is a major drawback to lotus callus proliferation and regeneration. However, the underlying mechanism of its formation remains largely unknown. Herein, we aimed to explore the metabolic and molecular basis of lotus callus browning by combining histological staining, high-throughput metabolomics, and transcriptomic assays for lotus callus at three browning stages. Histological stained brown callus cross sections displayed severe cell death symptoms, accompanied by an obvious accumulation of polyphenols and lignified materials. Widely targeted metabolomics revealed extensively decreased accumulation of most detected flavonoids and benzylisoquinoline alkaloids (BIAs), as well as a few phenolic acids, amino acids and their derivatives in callus with browning symptoms. Conversely, the contents of most detected tannins were significantly increased. Subsequent comparative transcriptomics identified a set of differentially expressed genes (DEGs) associated with the biosynthesis and regulation of flavonoids and BIAs in lotus. Notably, callus browning was coupled with significantly up-regulated expression of two polyphenol oxidase (PPO) and 17 peroxidase (POD) encoding genes, while the expression of ethylene associated genes remained at marginal levels. These results suggest that lotus callus browning is primarily controlled at the level of metabolism, wherein the oxidation of flavonoids and BIAs is crucially decisive.

Molecular characterization of SrSTP14, a sugar transporter from thermogenic skunk cabbage, and its possible role in developing pollen

In floral thermogenesis, sugars play an important role not only as energy providers but also as growth and development facilitators. Yet, the mechanisms underlying sugar translocation and transport in thermogenic plants remain to be studied. Asian skunk cabbage (Symplocarpus renifolius) is a species that can produce durable and intense heat in its reproductive organ, the spadix. Significant morphological and developmental changes in the stamen are well-characterized in this plant. In this study, we focused on the sugar transporters (STPs), SrSTP1 and SrSTP14, whose genes were identified by RNA-seq as the upregulated STPs during thermogenesis. Real-time PCR confirmed that mRNA expression of both STP genes was increased from the pre-thermogenic to the thermogenic stage in the spadix, where it is predominantly expressed in the stamen. SrSTP1 and SrSTP14 complemented the growth defects of a hexose transporter-deficient yeast strain, EBY4000, on media containing 0.02, 0.2, and 2% (w/v) glucose and galactose. Using a recently developed transient expression system in skunk cabbage leaf protoplasts, we revealed that SrSTP1 and SrSTP14-GFP fusion proteins were mainly localized to the plasma membrane. To dig further into the functional analysis of SrSTPs, tissue-specific localization of SrSTPs was investigated by in situ hybridization. Using probes for SrSTP14, mRNA expression was observed in the microspores within the developing anther at the thermogenic female stage. These results indicate that SrSTP1 and SrSTP14 transport hexoses (e.g., glucose and galactose) at the plasma membrane and suggest that SrSTP14 may play a role in pollen development through the uptake of hexoses into pollen precursor cells.

Comprehensive genome-wide identification and functional characterization of MAPK cascade gene families in Nelumbo

Mitogen-activated protein kinase (MAPK) cascade signaling pathway plays pivotal roles in various plant biological processes. However, systematic study of MAPK cascade gene families is yet to be conducted in lotus. Herein, 198 putative MAPK genes, including 152 MAP3Ks, 15 MKKs, and 31 MPKs genes were identified in Nelumbo. Segmental duplication was identified as the predominant factor driving MAPK cascade gene family expansion in lotus. MAPK cascade genes in N. nucifera and N. lutea shared high degree of sequence homologies, with 84, 9, and 19 homologous MAP3K, MKK, and MPK gene pairs being detected between the two species, respectively, with most genes predominantly undergoing purifying selection. Gene expression profiling indicated that NnMAPK cascade genes were extensively involved in plant development and submergence stress response. Co-expression analysis revealed potential interaction between transcription factors (TFs) and NnMAPK cascade genes in various biological processes. NnMKK showed predicted interactions with multiple NnMAP3K or NnMPK proteins, which suggested that functional diversity of MAPK cascade genes could be as a result of their complex protein interaction mechanisms. This first systematic analysis of MAPK cascade families in lotus provides deeper insights into their evolutionary dynamics and functional properties, which potentially could be crucial for lotus genetic improvement.

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Genome-wide study and functional characterization elucidates the potential association of late embryogenesis abundant (LEA) genes with lotus seed development

Late embryogenesis abundant (LEA) proteins are extremely hydrophilic proteins imperatively associated with plant growth and development, as well as cell protection from abiotic stress. However, the genome-wide characterization of LEA gene family remains limited, especially in aquatic species such as lotus (Nelumbo spp.). Here, 57 putative LEA genes, including 28 NnLEAs and 29 NlLEAs were identified in the N.nucifera and N.lutea genomes, respectively. A total of 27 homologous LEA gene pairs were identified, indicating high degree of sequence homologies between the two Nelumbo species. Secondary structure prediction indicated high prevalence of alpha (α) helix structure among LEA proteins in the LEA_1, LEA_4, and SMP groups. Screening of putative promoter cis-elements revealed that NnLEA genes were involved in diverse biological processes. Most NnLEA genes were predominantly expressed in the late cotyledons and plumules development stages, suggesting their potential vital roles in lotus seed maturation. In addition, genes co-expressed with NnLEAs were involved in ABA signaling, seed maturation, and development processes. Overall, this study provides new insights for the in-depth understanding of the functions of NnLEA proteins in lotus seed development, and could act as a useful reference for the molecular breeding of seeds with prolonged lifespan.

Metabolic profile and transcriptome reveal the mystery of petal blotch formation in rose

BACKGROUND

Petal blotch is a unique ornamental trait in angiosperm families, and blotch in rose petal is rare and has great esthetic value. However, the cause of the formation of petal blotch in rose is still unclear. The influence of key enzyme genes and regulatory genes in the pigment synthesis pathways needs to be explored and clarified.

RESULTS

In this study, the rose cultivar 'Sunset Babylon Eyes' with rose-red to dark red blotch at the base of petal was selected as the experimental material. The HPLC-DAD and UPLC-TQ-MS analyses indicated that only cyanidin 3,5-O-diglucoside (Cy3G5G) contributed to the blotch pigmentation of 'Sunset Babylon Eyes', and the amounts of Cy3G5G varied at different developmental stages. Only flavonols but no flavone were found in blotch and non-blotch parts. As a consequence, kaempferol and its derivatives as well as quercetin and its derivatives may act as background colors during flower developmental stages. Despite of the differences in composition, the total content of carotenoids in blotch and non-blotch parts were similar, and carotenoids may just make the petals show a brighter color. Transcriptomic data, quantitative real-time PCR and promoter sequence analyses indicated that RC7G0058400 (F3'H), RC6G0470600 (DFR) and RC7G0212200 (ANS) may be the key enzyme genes for the early formation and color deepening of blotch at later stages. As for two transcription factor, RC7G0019000 (MYB) and RC1G0363600 (WRKY) may bind to the promoters of critical enzyme genes, or RC1G0363600 (WRKY) may bind to the promoter of RC7G0019000 (MYB) to activate the anthocyanin accumulation in blotch parts of 'Sunset Babylon Eyes'.

CONCLUSIONS

Our findings provide a theoretical basis for the understanding of the chemical and molecular mechanism for the formation of petal blotch in rose.

Integrated metabolite profiling and transcriptome analysis reveal candidate genes involved in the formation of yellow Nelumbo nucifera

As a worldwide major ornamental flower and a edible plant, lotus (Nelumbo nucifera) is also used as medicine and tea beverage. Here, transcriptome and metabolites of yellow (MLQS) and white (YGB) lotus cultivars during five key flower coloration stages were profiled. 2014 differentially expressed genes were detected with 11 carotenoids in lotus were identified for the first time. Then, regulatory networks between and within functional modules was reconstructed, and the correlation between module-metabolites and gene-metabolites was conducted within 3 core modules. 18 candidate genes related to the formation of yellow flower were screened out and a gene regulatory model for the flower color difference between MLQS and YGB were speculated as follows: The substrate competition between F3'H and F3'5'H and substrate specificity of FLS, together with differential expression of CCD4a and CCD4b were contribute to the differences in flavonoids and carotenoids. Besides, UGT73C2, UGT91C1-2 and SGTase, and regulation of UGTs by transcription factors PLATZ, MADS, NAC031, and MYB308 may also play a role in the upstream regulation. The following verification results indicated that functional differences existed in the coding sequences of NnCCD4b and promoters of NnCCD4a of MLQS and YGB. In all, this study preliminarily reveals the mechanism of yellow flower coloration in lotus and provides new ideas for the study of complex ornamental characters of other plants.

Genome-Wide Association Study for Screening and Identifying Potential Shin Color Loci in Ducks

Shin color diversity is a widespread phenomenon in birds. In this study, ducks were assessed to identify candidate genes for yellow, black, and spotted tibiae. For this purpose, we performed whole-genome resequencing of an F2 population consisting of 275 ducks crossed between Runzhou crested-white ducks and Cherry Valley ducks. We obtained 12.6 Mb of single nucleotide polymorphism (SNP) data, and the three shin colors were subsequently genotyped. Genome-wide association studies (GWASs) were performed to identify candidate and potential SNPs for the three shin colors. According to the results, 2947 and 3451 significant SNPs were associated with black and yellow shins, respectively, and six potential SNPs were associated with spotted shins. Based on the SNP annotations, the MITF, EDNRB2, POU family members, and the SLC superfamily were the candidate genes regulating pigmentation. In addition, the isoforms of EDNRB2, TYR, TYRP1, and MITF-M were significantly different between the black and yellow tibiae. MITF and EDNRB2 may have synergistic roles in the regulation of melanin synthesis, and their mutations may lead to phenotypic differences in the melanin deposition between individuals. This study provides new insights into the genetic factors that may influence tibia color diversity in birds.

Genome-Wide Analysis Identifies Candidate Genes Encoding Beak Color of Duck

Beak color diversity is a broadly occurring phenomenon in birds. Here, we used ducks to identify candidate genes for yellow, black, and spotted beaks. For this, an F2 population consisting of 275 ducks was genotyped using whole genome resequencing containing 12.6 M single-nucleotide polymorphisms (SNPs) and three beak colors. Genome-wide association studies (GWAS) was used to identify the candidate and potential SNPs for three beak colors in ducks (yellow, spotted, and black). The results showed that 2753 significant SNPs were associated with black beaks, 7462 with yellow, and 17 potential SNPs with spotted beaks. Based on SNP annotation, MITF, EDNRB2, members of the POU family, and the SLC superfamily were the candidate genes regulating pigmentation. Meanwhile, isoforms MITF-M and EDNRB2 were significantly different between black and yellow beaks. MITF and EDNRB2 likely play a synergistic role in the regulation of melanin synthesis, and their mutations contribute to phenotypic differences in beak melanin deposition among individuals. This study provides new insights into genetic factors that may influence the diversity of beak color.

Genome-Wide Analysis Identifies Candidate Genes Encoding Feather Color in Ducks

Comparative population genomics and genome-wide association studies (GWAS) offer opportunities to discover human-driven detectable signatures within the genome. From the point of view of evolutionary biology, the identification of genes associated with the domestication of traits is of interest for the elucidation of the selection of these traits. To this end, an F2 population of ducks, consisting of 275 ducks, was genotyped using a whole genome re-sequence containing 12.6 Mb single nucleotide polymorphisms (SNPs) and four plumage colors. GWAS was used to identify the candidate and potential SNPs of four plumage colors in ducks (white, spot, grey, and black plumage). In addition, FST and genetic diversity (π ratio) were used to screen signals of the selective sweep, which relate to the four plumage colors. Major genomic regions associated with white, spotted, and black feathers overlapped with their candidate selection regions, whereas no such overlap was observed with grey plumage. In addition, MITF and EDNRB2 are functional candidate genes that contribute to white and black plumage due to their indirect involvement in the melanogenesis pathway. This study provides new insights into the genetic factors that may influence the diversity of plumage color.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Comprehensive Analysis for GRF Transcription Factors in Sacred Lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) is an aquatic perennial plant with essential food, ornamental, and pharmacological value. Growth-regulating factor (GRF) is a transcription factor (TF) family that plays an important role in regulating the growth and development of plants. In this study, a comprehensive analysis of the GRF family in N. nucifera was performed, and its role in N. nucifera development was studied. A total of eight GRF genes were identified in the N. nucifera genome. Phylogenetic analysis divided the 38 GRF genes into six clades, while the NuGRFs only contained five clades. The analyses of gene structures, motifs, and cis-acting regulatory elements of the GRF gene family were performed. In addition, the chromosome location and collinearity were analyzed. The expression pattern based on transcriptomic data and real-time reverse transcription-quantitative PCR (qRT-PCR) revealed that the GRF genes were expressed in multiple organs and were abundant in actively growing tissues, and the expression levels decreased as the age of N. nucifera increased. Then, 3D structures of the NuGRF proteins were predicted by homology modeling. Finally, the subcellular localization of GRF1 was ascertained in the tobacco leaf through a vector. Therefore, this study provides a comprehensive overview of the GRF TF family in N. nucifera.