Genomes of the Banyan Tree and Pollinator Wasp Provide Insights into Fig-Wasp Coevolution.

Banyan trees are distinguished by their extraordinary aerial roots. The Ficus genus includes species that have evolved a species-specific mutualism system with wasp pollinators. We sequenced genomes of the Chinese banyan tree, F. microcarpa, and a species lacking aerial roots, F. hispida, and one wasp genome coevolving with F. microcarpa, Eupristina verticillata. Comparative analysis of the two Ficus genomes revealed dynamic karyotype variation associated with adaptive evolution. Copy number expansion of auxin-related genes from duplications and elevated auxin production are associated with aerial root development in F. microcarpa. A male-specific AGAMOUS paralog, FhAG2, was identified as a candidate gene for sex determination in F. hispida. Population genomic analyses of Ficus species revealed genomic signatures of morphological and physiological coadaptation with their pollinators involving terpenoid- and benzenoid-derived compounds. These three genomes offer insights into and genomic resources for investigating the geneses of aerial roots, monoecy and dioecy, and codiversification in a symbiotic system.

Transcriptome dynamics provides insights into divergences of the photosynthesis pathway between Saccharum officinarum and Saccharum spontaneum.

Saccharum spontaneum and Saccharum officinarum contributed to the genetic background of modern sugarcane cultivars. Saccharum spontaneum has shown a higher net photosynthetic rate and lower soluble sugar than S. officinarum. Here, we analyzed 198 RNA-sequencing samples to investigate the molecular mechanisms for the divergences of photosynthesis and sugar accumulation between the two Saccharum species. We constructed gene co-expression networks based on differentially expressed genes (DEGs) both for leaf developmental gradients and diurnal rhythm. Our results suggested that the divergence of sugar accumulation may be attributed to the enrichment of major carbohydrate metabolism and the oxidative pentose phosphate pathway. Compared with S. officinarum, S. spontaneum DEGs showed a high enrichment of photosynthesis and contained more complex regulation of photosynthesis-related genes. Noticeably, S. spontaneum lacked gene interactions with sulfur assimilation stimulated by photorespiration. In S. spontaneum, core genes related to clock and photorespiration displayed a sensitive regulation by the diurnal rhythm and phase-shift. Small subunit of Rubisco (RBCS) displayed higher expression in the source tissues of S. spontaneum. Additionally, it was more sensitive under a diurnal rhythm, and had more complex gene networks than that in S. officinarum. This indicates that the differential regulation of RBCS Rubisco contributed to photosynthesis capacity divergence in both Saccharum species.

Haplotype-resolved genome of Mimosa bimucronata revealed insights into leaf movement and nitrogen fixation.

BACKGROUND: Mimosa bimucronata originates from tropical America and exhibits distinctive leaf movement characterized by a relative slow speed. Additionally, this species possesses the ability to fix nitrogen. Despite these intriguing traits, comprehensive studies have been hindered by the lack of genomic resources for M. bimucronata. RESULTS: To unravel the intricacies of leaf movement and nitrogen fixation, we successfully assembled a high-quality, haplotype-resolved, reference genome at the chromosome level, spanning 648 Mb and anchored in 13 pseudochromosomes. A total of 32,146 protein-coding genes were annotated. In particular, haplotype A was annotated with 31,035 protein-coding genes, and haplotype B with 31,440 protein-coding genes. Structural variations (SVs) and allele specific expression (ASE) analyses uncovered the potential role of structural variants in leaf movement and nitrogen fixation in M. bimucronata. Two whole-genome duplication (WGD) events were detected, that occurred ~ 2.9 and ~ 73.5 million years ago. Transcriptome and co-expression network analyses revealed the involvement of aquaporins (AQPs) and Ca2+-related ion channel genes in leaf movement. Moreover, we also identified nodulation-related genes and analyzed the structure and evolution of the key gene NIN in the process of symbiotic nitrogen fixation (SNF). CONCLUSION: The detailed comparative genomic and transcriptomic analyses provided insights into the mechanisms governing leaf movement and nitrogen fixation in M. bimucronata. This research yielded genomic resources and provided an important reference for functional genomic studies of M. bimucronata and other legume species.

Transcriptome and small RNA analysis unveils novel insights into the C4 gene regulation in sugarcane.

MAIN CONCLUSION: This study reveals miRNA indirect regulation of C4 genes in sugarcane through transcription factors, highlighting potential key regulators like SsHAM3a. C4 photosynthesis is crucial for the high productivity and biomass of sugarcane, however, the miRNA regulation of C4 genes in sugarcane remains elusive. We have identified 384 miRNAs along the leaf gradients, including 293 known miRNAs and 91 novel miRNAs. Among these, 86 unique miRNAs exhibited differential expression patterns, and we identified 3511 potential expressed targets of these differentially expressed miRNAs (DEmiRNAs). Analyses using Pearson correlation coefficient (PCC) and Gene Ontology (GO) enrichment revealed that targets of miRNAs with positive correlations are integral to chlorophyll-related photosynthetic processes. In contrast, negatively correlated pairs are primarily associated with metabolic functions. It is worth noting that no C4 genes were predicted as targets of DEmiRNAs. Our application of weighted gene co-expression network analysis (WGCNA) led to a gene regulatory network (GRN) suggesting miRNAs might indirectly regulate C4 genes via transcription factors (TFs). The GRAS TF SsHAM3a emerged as a potential regulator of C4 genes, targeted by miR171y and miR171am, and exhibiting a negative correlation with miRNA expression along the leaf gradient. This study sheds light on the complex involvement of miRNAs in regulating C4 genes, offering a foundation for future research into enhancing sugarcane's photosynthetic efficiency.

Upland rice genomic signatures of adaptation to drought resistance and navigation to molecular design breeding.

Upland rice is a distinctive drought-aerobic ecotype of cultivated rice highly resistant to drought stress. However, the genetic and genomic basis for the drought-aerobic adaptation of upland rice remains largely unclear due to the lack of genomic resources. In this study, we identified 25 typical upland rice accessions and assembled a high-quality genome of one of the typical upland rice varieties, IRAT109, comprising 384 Mb with a contig N50 of 19.6 Mb. Phylogenetic analysis revealed upland and lowland rice have distinct ecotype differentiation within the japonica subgroup. Comparative genomic analyses revealed that adaptive differentiation of lowland and upland rice is likely attributable to the natural variation of many genes in promoter regions, formation of specific genes in upland rice, and expansion of gene families. We revealed differentiated gene expression patterns in the leaves and roots of the two ecotypes and found that lignin synthesis mediated by the phenylpropane pathway plays an important role in the adaptive differentiation of upland and lowland rice. We identified 28 selective sweeps that occurred during domestication and validated that the qRT9 gene in selective regions can positively regulate drought resistance in rice. Eighty key genes closely associated with drought resistance were appraised for their appreciable potential in drought resistance breeding. Our study enhances the understanding of the adaptation of upland rice and provides a genome navigation map of drought resistance breeding, which will facilitate the breeding of drought-resistant rice and the "blue revolution" in agriculture.

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.

Gene Regulatory Network Controlling Flower Development in Spinach (Spinacia oleracea L.).

Spinach (Spinacia oleracea L.) is a dioecious, diploid, wind-pollinated crop cultivated worldwide. Sex determination plays an important role in spinach breeding. Hence, this study aimed to understand the differences in sexual differentiation and floral organ development of dioecious flowers, as well as the differences in the regulatory mechanisms of floral organ development of dioecious and monoecious flowers. We compared transcriptional-level differences between different genders and identified differentially expressed genes (DEGs) related to spinach floral development, as well as sex-biased genes to investigate the flower development mechanisms in spinach. In this study, 9189 DEGs were identified among the different genders. DEG analysis showed the participation of four main transcription factor families, MIKC_MADS, MYB, NAC, and bHLH, in spinach flower development. In our key findings, abscisic acid (ABA) and gibberellic acid (GA) signal transduction pathways play major roles in male flower development, while auxin regulates both male and female flower development. By constructing a gene regulatory network (GRN) for floral organ development, core transcription factors (TFs) controlling organ initiation and growth were discovered. This analysis of the development of female, male, and monoecious flowers in spinach provides new insights into the molecular mechanisms of floral organ development and sexual differentiation in dioecious and monoecious plants in spinach.

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

Nelumbo lutea (American lotus), which differs from Nelumbo nucifera (Asian lotus) morphologically, is one of the two remaining species in the basal eudicot family Nelumbonaceae. Here, we assembled the 843-Mb genome of American lotus into eight pseudochromosomes containing 31 382 protein-coding genes. Comparative analyses revealed conserved synteny without large chromosomal rearrangements between the genomes of American and Asian lotus and identified 29 533 structural variants (SVs). Carotenoid and anthocyanin pigments determine the yellow and red petal colors of American and Asian lotus, respectively. The structural genes encoding enzymes of the carotenoid and anthocyanin biosynthesis pathways were conserved between two species but differed in expression. We detected SVs caused by repetitive sequence expansion or contraction among the anthocyanin biosynthesis regulatory MYB genes. Further transient overexpression of candidate NnMYB5 induced anthocyanin accumulation in lotus petals. Alternative oxidase (AOX), uncoupling proteins (UCPs), and sugar metabolism and transportation contributed to carpel thermogenesis. Carpels produce heat with sugars transported from leaves as the main substrates, because there was weak tonoplast sugar transporter (TST) activity, and with SWEETs were highly expressed during thermogenesis. Cell proliferation-related activities were particularly enhanced in the warmer carpels compared with stamens during the cold night before blooming, which suggested that thermogenesis plays an important role in flower protogyny. Population genomic analyses revealed deep divergence between American and Asian lotus, and independent domestication affecting seed, rhizome, and flower traits. Our findings provide a high-quality reference genome of American lotus for exploring the genetic divergence and variation between two species and revealed possible genomic bases for petal color, carpel thermogenesis and domestication in lotus.

Responses of roots and rhizosphere of female papaya to the exogenous application of GA3.

Exogenous GAs have an indeterminate effect on root development. Our current study used female papaya to reveal how the roots and rhizosphere respond to the exogenous application of GA3 by investigating the transcriptome profile in roots, metabolic profile and microbial community in both roots and rhizosphere of GA3-treated and control female papaya. The results demonstrated that exogenous GA3 treatment enhanced female papaya lateral root development, which gave plants physical advantages of water and nutrient uptake. In addition, it was likely that GA3 spraying in papaya shoot apices increased the level of auxin, which was transported to roots by CpPIN1, where auxin upregulated CpLBD16 and repressed CpBP to promote the lateral root initiation and development. In papaya roots, corresponding transporters (CpTMT3, CpNRT1:2, CpPHT1;4, CpINT2, CpCOPT2, CpABCB11, CpNIP4;1) were upregulated and excretion transporters were downregulated such as CpNAXT1 for water and nutrients uptake with exogenous GA3 application. Moreover, in GA3-treated papaya roots, CpALS3 and CpMYB62 were downregulated, indicating a stronger abiotic resistance to aluminum toxic and phosphate starvation. On the other hand, BRs and JAs, which involve in defense responses, were enriched in the roots and rhizosphere of GA3-treated papayas. The upregulation of the two hormones might result in the reduction of pathogens in roots and rhizosphere such as Colletotrichum and Verticillium. GA3-treated female papaya increased the abundance of beneficial bacteria species including Mycobacterium, Mitsuaria, and Actinophytocola, but decreased that of the genera Candidatus and Bryobacter for that it required less nitrate. Overall, the roots and rhizosphere of female papaya positively respond to exogenous application of GA3 to promote development and stress tolerance. Treatment of female papaya with GA3 might result in the promotion of lateral root formation and development by upregulating CpLBD16 and downregulating CpBP. GA3-treated papaya roots exhibited feedback control of brassinolide and jasmonate signaling in root development and defense. These findings revealed complex response to a growth hormone treatment in papaya roots and rhizosphere and will lead to investigations on the impact of other plant hormones on belowground development in papaya.

Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection.

BACKGROUND: Diseases are the major factor affecting the quality and yield of sugarcane during its growth and development. However, our knowledge about the factors regulating disease responses remain limited. The present study focuses on identifying genes regulating transcriptional mechanisms responsible for resistance to leaf scald caused by Xanthomonas albilineans in S. spontaneum and S. officinarum. RESULTS: After inoculation of the two sugarcane varieties SES208 (S. spontaneum) and LA Purple (S. officinarum) with Xanthomonas albilineans, SES208 exhibited significantly greater resistance to leaf scald caused by X. albilineans than did LA Purple. Using transcriptome analysis, we identified a total of 4323 and 1755 differentially expressed genes (DEGs) in inoculated samples of SES208 and LA Purple, respectively. Significantly, 262 DEGs were specifically identified in SES208 that were enriched for KEGG pathway terms such as plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction. Furthermore, we built a transcriptional regulatory co-expression network that specifically identified 16 and 25 hub genes in SES208 that were enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction. All of these essential genes might be significantly involved in resistance-regulating responses in SES208 after X. albilineans inoculation. In addition, we found allele-specific expression in SES208 that was associated with the resistance phenotype of SES208 when infected by X. albilineans. After infection with X. albilineans, a great number of DEGs associated with the KEGG pathways 'phenylpropanoid biosynthesis' and 'flavonoid biosynthesis' exhibited significant expression changes in SES208 compared to LA Purple that might contribute to superior leaf scald resistance in SES208. CONCLUSIONS: We provided the first systematical transcriptome map that the higher resistance of SES208 is associated with and elicited by the rapid activation of multiple clusters of defense response genes after infection by X. albilineans and not merely due to changes in the expression of genes generically associated with stress resistance. These results will serve as the foundation for further understanding of the molecular mechanisms of resistance against X. albilineans in S. spontaneum.

Reward history modulates visual attention in an avian model.

Attention can be biased towards previously reward-associated stimuli even when they are task-irrelevant and physically non-salient, although studies of reward-modulated attention have been largely limited to primate (including human and nonhuman) models. Birds have been shown to have the capacity to discriminate reward and spatial cues in a manner similar to primates, but whether reward history involuntarily affects their attention in the same way remains unclear. We adapted a spatial cueing paradigm with differential rewards to investigate how reward modulates the allocation of attention in peafowl (Pavo cristatus). The birds were required to locate and peck a target on a computer screen that was preceded by a high-value or low-value color cue that was uninformative with respect to the location of the upcoming target. All birds exhibited a validity effect (performance enhanced on valid compared to invalid cue), and an interaction effect between value and validity was evident at the group level, being particularly pronounced in the birds with the greatest amount of reward training. The time course of reward learning was conspicuously incremental, phenomenologically slower compared to primates. Our findings suggest a similar influence of reward history on attention across phylogeny despite a significant difference in neuroanatomy.

The evolutionary origin and domestication history of goldfish (Carassius auratus).

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.

Genome-Wide Identification and Expression Pattern Profiling of the Aquaporin Gene Family in Papaya (Carica papaya L.).

Aquaporins (AQPs) are mainly responsible for the transportation of water and other small molecules such as CO2 and H2O2, and they perform diverse functions in plant growth, in development, and under stress conditions. They are also active participants in cell signal transduction in plants. However, little is known about AQP diversity, biological functions, and protein characteristics in papaya. To better understand the structure and function of CpAQPs in papaya, a total of 29 CpAQPs were identified and classified into five subfamilies. Analysis of gene structure and conserved motifs revealed that CpAQPs exhibited a degree of conservation, with some differentiation among subfamilies. The predicted interaction network showed that the PIP subfamily had the strongest protein interactions within the subfamily, while the SIP subfamily showed extensive interaction with members of the PIP, TIP, NIP, and XIP subfamilies. Furthermore, the analysis of CpAQPs' promoters revealed a large number of cis-elements participating in light, hormone, and stress responses. CpAQPs exhibited different expression patterns in various tissues and under different stress conditions. Collectively, these results provided a foundation for further functional investigations of CpAQPs in ripening, as well as leaf, flower, fruit, and seed development. They also shed light on the potential roles of CpAQP genes in response to environmental factors, offering valuable insights into their biological functions in papaya.

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum.

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.

Evolutionary expansion and functional divergence of sugar transporters in Saccharum (S. spontaneum and S. officinarum).

The sugar transporter (ST) family is considered to be the most important gene family for sugar accumulation, but limited information about the ST family in the important sugar-yielding crop Saccharum is available due to its complex genetic background. Here, 105 ST genes were identified and clustered into eight subfamilies in Saccharum spontaneum. Comparative genomics revealed that tandem duplication events contributed to ST gene expansions of two subfamilies, PLT and STP, in S. spontaneum, indicating an early evolutionary step towards high sugar content in Saccharum. The analyses of expression patterns were based on four large datasets with a total of 226 RNA sequencing samples from S. spontaneum and Saccharum officinarum. The results clearly demonstrated 50 ST genes had different spatiotemporal expression patterns in leaf tissues, 10 STs were specifically expressed in the stem, and 10 STs responded to the diurnal rhythm. Heterologous expression experiments in the defective yeast strain EBY.VW4000 indicated STP13, pGlcT2, VGT3, and TMT4 are the STs with most affinity for glucose/fructose and SUT1_T1 has the highest affinity to sucrose. Furthermore, metabolomics analysis suggested STP7 is a sugar starvation-induced gene and STP13 has a function in retrieving sugar in senescent tissues. PLT11, PLT11_T1, TMT3, and TMT4 contributed to breaking the limitations of the storage sink. SUT1, SUT1_T1, PLT11, TMT4, pGlcT2, and VGT3 responded for different functions in these two Saccharum species. This study demonstrated the evolutionary expansion and functional divergence of the ST gene family and will enable the further investigation of the molecular mechanism of sugar metabolism in Saccharum.

The genome of Shanputao (Vitis amurensis) provides a new insight into cold tolerance of grapevine.

Vitis amurensis (Shanputao) is the most cold tolerant Vitis species and so is of great interest to grape breeders and producers in areas with low winter temperatures. Here, we report its high-quality, chromosome-level genome assembly based on a combination of sequence data from Illumina and PacBio platforms, BioNano optical mapping and high-throughput chromosome conformation Capture (Hi-C) mapping. The 604.56-Mb genome contains 32 885 protein-coding genes. Shanputao was found to share a common ancestor with PN40024 (V. vinifera) approximately 2.17-2.91 million years ago, and gene expansion observed in Shanputao might contribute to the enhancement of cold tolerance. Transcriptome analysis revealed 17 genes involved in cold signal transduction, suggesting that there was a different response mechanism to chilling temperature and freezing conditions. Furthermore, a genome-wide association study uncovered a phosphoglycerate kinase gene that may contribute to the freezing resistance of buds in the winter. The Shanputao genome sequence not only represents a valuable resource for grape breeders, but also is important for clarifying the molecular mechanisms involved in cold tolerance.

Rambutan genome revealed gene networks for spine formation and aril development.

Rambutan is a popular tropical fruit known for its exotic appearance, has long flexible spines on shells, extraordinary aril growth, desirable nutrition, and a favorable taste. The genome of an elite rambutan cultivar Baoyan 7 was assembled into 328 Mb in 16 pseudo-chromosomes. Comparative genomics analysis between rambutan and lychee revealed that rambutan chromosomes 8 and 12 are collinear with lychee chromosome 1, which resulted in a chromosome fission event in rambutan (n = 16) or a fusion event in lychee (n = 15) after their divergence from a common ancestor 15.7 million years ago. Root development genes played a crucial role in spine development, such as endoplasmic reticulum pathway genes, jasmonic acid response genes, vascular bundle development genes, and K+ transport genes. Aril development was regulated by D-class genes (STK and SHP1), plant hormone and phenylpropanoid biosynthesis genes, and sugar metabolism genes. The lower rate of male sterility of hermaphroditic flowers appears to be regulated by MYB24. Population genomic analyses revealed genes in selective sweeps during domestication that are related to fruit morphology and environment stress response. These findings enhance our understanding of spine and aril development and provide genomic resources for rambutan improvement.

Gene regulation network analyses of pistil development in papaya.

BACKGROUND: The pistil is an essential part of flowers that functions in the differentiation of the sexes and reproduction in plants. The stigma on the pistil can accept pollen to allow fertilization and seed development. Papaya (Carica papaya L.) is a dioecious plant, where female flowers exhibit normal pistil, while the male flowers exhibit aborted pistil at a late stage of pistil development. RESULTS: The developmental stages of papaya pistil were analyzed after first dividing it into slices representing the primordium stage 1 (S1), the pre-meiotic stages S2, post-meiotic stage S3, and the mitotic stage S4. The SS scoring algorithm analysis of genes preferentially expressed at different stages revealed differentially expressed genes between male and female flowers. A transcription factor regulatory network for each stage based on the genes that are differentially expressed between male and female flowers was constructed. Some transcription factors related to pistil development were revealed based on the analysis of regulatory networks such as CpAGL11, CpHEC2, and CpSUPL. Based on the specific expression of genes, constructed a gene regulatory subnetwork with CpAGL11-CpSUPL-CpHEC2 functioning as the core. Analysis of the functionally enriched terms in this network reveals several differentially expressed genes related to auxin/ brassinosteroid signal transduction in the plant hormone signal transduction pathway. At the same time, significant differences in the expression of auxin and brassinosteroid synthesis-related genes between male and female flowers at different developmental stages were detected. CONCLUSIONS: The pistil abortion of papaya might be caused by the lack of expression or decreased expression of some transcription factors and hormone-related genes, affecting hormone signal transduction or hormone biosynthesis. Analysis of aborted and normally developing pistil in papaya provided new insights into the molecular mechanism of pistil development and sex differentiation in dioecious papaya.

Transcriptome analyses shed light on floral organ morphogenesis and bract color formation in Bougainvillea.

BACKGROUND: Bougainvillea is a popular ornamental plant with brilliant color and long flowering periods. It is widely distributed in the tropics and subtropics. The primary ornamental part of the plant is its colorful and unusual bracts, rich in the stable pigment betalain. The developmental mechanism of the bracts is not clear, and the pathway of betalain biosynthesis is well characterized in Bougainvillea. RESULTS: At the whole-genome level, we found 23,469 protein-coding genes by assembling the RNA-Seq and Iso-Seq data of floral and leaf tissues. Genome evolution analysis revealed that Bougainvillea is related to spinach; the two diverged approximately 52.7 million years ago (MYA). Transcriptome analysis of floral organs revealed that flower development of Bougainvillea was regulated by the ABCE flower development genes; A-class, B-class, and E-class genes exhibited high expression levels in bracts. Eight key genes of the betalain biosynthetic pathway were identified by homologous alignment, all of which were upregulated concurrently with bract development and betalain accumulation during the bract initiation stage of development. We found 47 genes specifically expressed in stamens, including seven highly expressed genes belonging to the pentose and glucuronate interconversion pathways. BgSEP2b, BgSWEET11, and BgRD22 are hub genes and interacted with many transcription factors and genes in the carpel co-expression network. CONCLUSIONS: We assembled protein-coding genes of Bougainvilea, identified the floral development genes, and constructed the gene co-expression network of petal, stamens, and carpel. Our results provide fundamental information about the mechanism of flower development and pigment accumulation in Bougainvillea, and will facilitate breeding of cultivars with high ornamental value.

The genomic architecture of the sex-determining region and sex-related metabolic variation in Ginkgobiloba.

Sex differences and evolutionary differences are critical biological issues. Ginkgo is an ancient lineage of dioecious gymnosperms with special value for studying the mechanism of sex determination in plants. However, the major genetic basic underlying sex chromosomes remains to be uncovered. In this study, we identify the sex-determining region of Ginkgo and locate it to the area from megabases 48 to 75 on chromosome 2. We find that the male sex-determining region of Ginkgo contains more than 200 genes, including four MADS-box genes, demonstrating that the Ginkgo sex determination system is of the XY type. We also find that genetic sex differences result in specialized flavonoid metabolism and regulation in each sex. These findings establish a foundation for revealing the molecular mechanism of sexual dimorphism and promoting the development of the Ginkgo industry.