Differential gene expression among three sex types reveals a MALE STERILITY 1 (CpMS1) for sex differentiation in papaya

Comparative proteomic analysis of papaya bud flowers reveals metabolic signatures and pathways driving hermaphrodite development

Papaya (Carica papaya) is a trioecious species with female, male, and hermaphrodite plants. Given the sex segregation, selecting hermaphroditic plants is vital for orchard establishment due to their greater commercial value. However, selecting hermaphrodite plants through sexing is laborious and costly. Moreover, environmental stressors can exacerbate the issue by potentially inducing abnormal flower development, thus affecting fruit quality. Despite these challenges, the molecular mechanisms governing sex development in papaya remain poorly understood. Thus, this study aimed to identify proteins associated with sex development in female and hermaphrodite flowers of papaya through comparative proteomic analysis. Proteins from flower buds at the early and late developmental stages of three papaya genotypes (UENF-CALIMAN 01, JS12, and Sunrise Solo 72/12) were studied via proteomic analysis via the combination of the shotgun method and nanoESI-HDMSE technology. In buds at an early stage of development, 496 (35.9%) proteins exhibited significantly different abundances between sexes for the SS72/12 genotype, 139 (10%) for the JS12 genotype, and 165 (11.9%) for the UC-01 genotype. At the final stage of development, there were 181 (13.5%) for SS72/12, 113 (8.4%) for JS12, and 125 (9.1%) for UC-01. The large group of differentially accumulated proteins (DAPs) between the sexes was related to metabolism, as shown by the observation of only the proteins that exhibited the same pattern of accumulation in the three genotypes. Specifically, carbohydrate metabolism proteins were up-regulated in hermaphrodite flower buds early in development, while those linked to monosaccharide and amino acid metabolism increased during late development. Enrichment of sporopollenin and phenylpropanoid biosynthesis pathways characterizes hermaphrodite samples across developmental stages, with predicted protein interactions highlighting the crucial role of phenylpropanoids in sporopollenin biosynthesis for pollen wall formation. Most of the DAPs played key roles in pectin, cellulose, and lignin synthesis and were essential for cell wall formation and male flower structure development, notably in the pollen coat. These findings suggest that hermaphrodite flowers require more energy for development, likely due to complex pollen wall formation. Overall, these insights illuminate the molecular mechanisms of papaya floral development, revealing complex regulatory networks and energetic demands in the formation of male reproductive structures.

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.

Characterization and analysis of the promoter region of monodehydroascorbate reductase 4 (CpMDAR4) in papaya

Differential spatial and temporal expression patterns due to regulatory cis-elements and two different isoforms are detected among CpMDAR4 alleles in papaya. The aim of this research was to study the effects of cis-element differences between the X, Y and Y^h alleles on the expression of CpMDAR4, a potential candidate gene for sex differentiation in papaya, using a transcriptional reporter system in a model species Arabidopsis thaliana. Possible effects of a retrotransposon insertion in the Y and Y^h alleles on the transcription and expression of CpMDAR4 alleles in papaya flowers were also examined. When comparing promoters and cis-regulatory elements among genes in the non-recombining region of the sex chromosomes, paired genes exhibited differences. Our results showed that differences in the promoter sequences of the CpMDAR4 alleles drove the expression of a reporter gene to different flower tissues in Arabidopsis. β-glucuronidase staining analysis of T₂ and T₃ lines for constructs containing 5' deletions of native Y and Y^h allele promoters showed the loss of specific expression of the reporter gene in the anthers, confirming the existence and location of cis-regulatory element POLLEN1LELAT52. The expression analysis of CpMDAR4 alleles in papaya flowers also showed that all alleles are actively expressed in different flower tissues, with the existence of a shorter truncated isoform, with unknown function, for the Y and Y^h alleles due to an LTR-RT insertion in the Y and Y^h chromosomes. The observed expression patterns in Arabidopsis thaliana flowers and the expression patterns of CpMDAR4 alleles in papaya flowers suggest that MDAR4 might have a role on development of reproductive organs in papaya, and that it constitutes an important candidate for sex differentiation.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

Integrative Analyses of Transcriptomics and Metabolomics in Sex Differentiation of Mulberry Flowers

Sex determination and sex differentiation of plants are important physiological processes of plant development. Mulberry (Morus indica L.) is an important economic tree being cultivated in sericulture countries, and mulberry leaf is commonly used for sericulture. The transcriptomic and metabolomic differences between the staminate flowers (SFs) and pistillate flowers (PFs) of mulberry were investigated by RNA sequencing and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). Overall, we uncovered 4,230 genes and 209 metabolites are significantly differentially expressed between the SFs and PFs of mulberry. The combined transcriptomic and metabolomic analysis revealed these differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) are involved in flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism, and these detected DEGs and DEMs may be associated with sex differentiation of mulberry through the regulation of the enrichment pathways, such as the MAPK pathway, flavonoid biosynthesis, galactose metabolism, plant–pathogen interaction, and starch and sucrose metabolism. This study will provide a rich source for the analysis of the molecular mechanism of mulberry sex differentiation processes.

Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower

Sex determination and differentiation is an important biological process for unisexual flower development. Spinach is a model plant to study the mechanism of sex determination and differentiation of dioecious plant. Till now, little is known about spinach sex determination and differentiation mechanism. MicroRNAs are key factors in flower development. Herein, small RNA sequencing was performed to explore the roles of microRNAs in spinach sex determination and differentiation. As a result, 92 known and 3402 novel microRNAs were identified in 18 spinach female and male flower samples. 74 differentially expressed microRNAs were identified between female and male flowers, including 20 female-biased and 48 male-biased expression microRNAs. Target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Among the differentially expressed microRNAs between FNS and M03, 55 microRNAs were found to reside in sex chromosome; one of them, sol-miR2550n, was functionally studied via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n was a male-promoting factor and validating the reliability of our small RNA sequencing data. Conclusively, this work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

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.

The transcriptional landscape of the giant freshwater prawn: Embryonic development and early sexual differentiation mechanisms

The giant freshwater prawn pjMacrobrachium rosenbergii is one of the best studied species in aquaculture. However, the transcriptional changes associated with embryonic development and the sexual differentiation mechanism of M. rosenbergii remain to be elucidated. To characterize the embryonic development of this prawn and to determine whether differential expression and differential splicing play roles in the early sexual differentiation of M. rosenbergii, we profiled five developmental days of male and female embryos by RNA sequencing. We identified modules of co-expressed genes representing waves of transcription that correspond to physiological processes in early embryonic development (such as the maternal-to-zygotic transition) up to preparation for life outside the egg (development of muscles, cuticle etc.). Additionally, we found that hundreds of genes are differentially expressed between sexes, most of them uncharacterized, suggesting that the sex differentiation mechanism of M. rosenbergii might contain clade-specific elements. The resulting first-of-a-kind transcriptional map of embryonic development of male and female M. rosenbergii will guide future studies to reveal the roles of specific genes and splicing isoforms in the embryonic development and sexual differentiation process of M. rosenbergii.

Sex-linked gene expression and the emergence of hermaphrodites in Carica papaya

PREMISE

One evolutionary path from hermaphroditism to dioecy is via a gynodioecious intermediate. The evolution of dioecy may also coincide with the formation of sex chromosomes that possess sex-determining loci that are physically linked in a region of suppressed recombination. Dioecious papaya (Carica papaya) has an XY chromosome system, where the presence of a Y chromosome determines maleness. However, in cultivation, papaya is gynodioecious, due to the conversion of the male Y chromosome to a hermaphroditic Y^h chromosome during its domestication.

METHODS

We investigated gene expression linked to the X, Y, and Y^h chromosomes at different floral developmental stages to identify differentially expressed genes that may be involved in the sexual transition of males to hermaphrodites.

RESULTS

We identified 309 sex-biased genes found on the sex chromosomes, most of which are found in the pseudoautosomal regions. Female (XX) expression in the sex-determining region was almost double that of X-linked expression in males (XY) and hermaphrodites (XY^h ), which rules out dosage compensation for most sex-linked genes; although, an analysis of hemizygous X-linked loci found evidence of partial dosage compensation. Furthermore, we identified a candidate gene associated with sex determination and the transition to hermaphroditism, a homolog of the MADS-box protein SHORT VEGETATIVE PHASE.

CONCLUSIONS

We identified a pattern of partial dosage compensation for hemizygous genes located in the papaya sex-determining region. Furthermore, we propose that loss-of-expression of the Y-linked SHORT VEGETATIVE PHASE homolog facilitated the transition from males to hermaphrodites in papaya.