SunUp and Sunset genomes revealed impact of particle bombardment mediated transformation and domestication history in papaya

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion

BACKGROUND

Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions.

RESULTS

Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP".

CONCLUSION

The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

Sexy ways: approaches to studying plant sex chromosomes

Sex chromosomes have evolved in many plant species with separate sexes. Current plant research is shifting from examining the structure of sex chromosomes to exploring their functional aspects. New studies are progressively unveiling the specific genetic and epigenetic mechanisms responsible for shaping distinct sexes in plants. While the fundamental methods of molecular biology and genomics are generally employed for the analysis of sex chromosomes, it is often necessary to modify classical procedures not only to simplify and expedite analyses but sometimes to make them possible at all. In this review, we demonstrate how, at the level of structural and functional genetics, cytogenetics, and bioinformatics, it is essential to adapt established procedures for sex chromosome analysis.

Molecular mechanism of vivipary as revealed by the genomes of viviparous mangroves and non-viviparous relatives

Vivipary is a prominent feature of mangroves, allowing seeds to complete germination while attached to the mother plant, and equips propagules to endure and flourish in challenging coastal intertidal wetlands. However, vivipary-associated genetic mechanisms remain largely elusive. Genomes of two viviparous mangrove species and a non-viviparous inland relative were sequenced and assembled at the chromosome level. Comparative genomic analyses between viviparous and non-viviparous genomes revealed that DELAY OF GERMINATION 1 (DOG1) family genes (DFGs), the proteins from which are crucial for seed dormancy, germination, and reserve accumulation, are either lost or dysfunctional in the entire lineage of true viviparous mangroves but are present and functional in their inland, non-viviparous relatives. Transcriptome dynamics at key stages of vivipary further highlighted the roles of phytohormonal homeostasis, proteins stored in mature seeds, and proanthocyanidins in vivipary under conditions lacking DFGs. Population genomic analyses elucidate dynamics of syntenic regions surrounding the missing DFGs. Our findings demonstrated the genetic foundation of constitutive vivipary in Rhizophoraceae mangroves.

Multi-omic applications for understanding and enhancing tropical fruit flavour

Consumer trends towards nutrient-rich foods are contributing to global increasing demand for tropical fruit. However, commercial cultivars in the breeding pipeline that are tailored to meet market demand are at risk of possessing reduced fruit flavour qualities. This stems from recurrent prioritised selection for superior agronomic traits and not fruit flavour, which may in turn reduce consumer satisfaction. There is realisation that fruit quality traits, inclusive of flavour, must be equally selected for; but currently, there are limited tools and resources available to select for fruit flavour traits, particularly in tropical fruit species. Although sugars, acids, and volatile organic compounds are known to define fruit flavour, the specific combinations of these, that result in defined consumer preferences, remain unknown for many tropical fruit species. To define and include fruit flavour preferences in selective breeding, it is vital to determine the metabolites that underpin them. Then, objective quantitative analysis may be implemented instead of solely relying on human sensory panels. This may lead to the development of selective genetic markers through integrated omics approaches that target biosynthetic pathways of flavour active compounds. In this review, we explore progress in the development of tools to be able to strategically define and select for consumer-preferred flavour profiles in the breeding of new cultivars of tropical fruit species.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

Carica papaya L. sex chromosome review and physical mapping of the serk 2, svp-like and mdar 4 sequences

Physical mapping evidences the chromosome organization and structure. Despite the data about plant cytogenomics, physical mapping has been conducted from single-copy and/or low-copy genes for few species. Carica papaya cytogenomics has been accomplished from BAC-FISH and repeatome sequences. We aimed to map the serk 2, svp-like and mdar 4 sequences in C. papaya. The sequences were amplified and the amplicons sequenced, showing similarity in relation to serk 2, svp-like and mdar 4 genes. Carica papaya diploidy was confirmed and the mitotic chromosomes characterized. The chromosome 1 exhibited the secondary constriction pericentromeric to the centromere of the long arm. So, we concluded that it is the sex chromosomes. serk 2 was mapped in the long arm interstitial portion of the sex chromosomes, and the interphase nuclei showed two fluorescence signals. Considering these results and the sequencing data from the C. papaya sex chromosomes, svp-like and mdar 4 genes were mapped in the interstitial region of the sex chromosome long arm. Both sequences showed only one fluorescence signal in the interphase nuclei. The procedure adopted here can be reproduced for other single-copy and/or low-copy genes, allowing the construction of cytogenetic maps. In addition, we revisited the cytogenomics data about C. papaya sex chromosomes, presenting a revised point of view about the structure and evolution to these chromosomes.

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

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.

Identification and Validation of Key Genes Related to Preferred Flavour Profiles in Australian Commercial Papaya (Carica papaya L.)

Commercial papaya varieties grown in Australia vary greatly in taste and aroma. Previous profiling has identified undesirable 'off tastes' in existing varieties, discouraging a portion of the population from consuming papayas. Our focus on enhancing preferred flavours led to an exploration of the genetic mechanisms and biosynthesis pathways that underlie these desired taste profiles. To identify genes associated with consumer-preferred flavours, we conducted whole RNA sequencing and de novo genome assembly on papaya varieties RB1 (known for its sweet flavour and floral aroma) and 1B (less favoured due to its bitter taste and musty aroma) at both ripe and unripe stages. In total, 180,368 transcripts were generated, and 118 transcripts related to flavours were differentially expressed between the two varieties at the ripe stage. Five genes (cpBGH3B, cpPFP, cpSUS, cpGES and cpLIS) were validated through qPCR and significantly differentially expressed. These genes are suggested to play key roles in sucrose metabolism and aromatic compound production pathways, holding promise for future selective breeding strategies. Further exploration will involve assessing their potential across broader germplasm and various growth environments.

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.

Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata)

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

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.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

Microparticle-mediated CRISPR DNA delivery for genome editing in poplar

The use of CRISPR/Cas9 is currently the method of choice for precise genome engineering in plants, including in the biomass crop poplar. The most commonly used method for delivering CRISPR/Cas9 and its components in poplar is via Agrobacterium-mediated transformation, that besides the desired gene-editing event also results in stable T-DNA integration. Here we explore the delivery of the gene-editing reagents via DNA-coated microparticle bombardment into the model tree Populus tremula x P. alba to evaluate its potential for developing transgene-free, gene-edited trees, as well as its potential for integrating donor DNA at specific target sites. Using an optimized transformation method, which favors the regeneration of plants that transiently express the genes on the delivered donor DNA, we regenerated gene-edited plants that are free of the Cas9 and the antibiotic resistance-encoding transgenes. In addition, we report the frequent integration of donor DNA fragments at the Cas9-induced double-strand break, opening opportunities toward targeted gene insertions.

A perspective from the EU: unintended genetic changes in plants caused by NGT-their relevance for a comprehensive molecular characterisation and risk assessment

Several regions in the world are currently holding discussions in regard to the regulation of new genomic techniques (NGTs) and their application in agriculture. The European Commission, for instance, is proposing the introduction of specific regulation for NGT plants. Various questions need to be answered including e.g., the extent to which NGT-induced intended and unintended genetic modifications must be subjected to a mandatory risk assessment as part of an approval procedure. This review mostly focuses on findings in regard to unintended genetic changes that can be caused by the application of NGTs. More specifically, the review deals with the application of the nuclease CRISPR/Cas, which is currently the most important tool for developing NGT plants, and its potential to introduce double strand breaks (DSBs) at a targeted DNA sequence. For this purpose, we identified the differences in comparison to non-targeted mutagenesis methods used in conventional breeding. The review concludes that unintended genetic changes caused by NGT processes are relevant to risk assessment. Due to the technical characteristics of NGTs, the sites of the unintended changes, their genomic context and their frequency (in regard to specific sites) mean that the resulting gene combinations (intended or unintended) may be unlikely to occur with conventional methods. This, in turn, implies that the biological effects (phenotypes) can also be different and may cause risks to health and the environment. Therefore, we conclude that the assessment of intended as well as unintended genetic changes should be part of a mandatory comprehensive molecular characterisation and risk assessment of NGT plants that are meant for environmental releases or for market authorisation.

Sex determination in papaya: Current status and perspectives

Papaya (Carica papaya L.) is an economically significant plant that produces fruit consumed worldwide due to its organoleptic characteristics. Since their commercial production, papaya fruits have faced several problems, such as pests, which have been partly resolved using transgenic varieties. Nevertheless, a principal challenge in this cultivation is the plant's sex determination. The sex issue in papaya is complex because papaya flowers can bear three sex forms: male, female, and hermaphrodite, which affects their fruit production, shape, and yield. Fruits from hermaphrodite plants are preferred more by consumers than female ones, and male plants rarely produce fruits without commercial value. Chromosomes are responsible for sex determination in papaya, denoted as XY for male, XX for female, and XYh for hermaphrodite. However, genes related to sex have been reported but are not conclusive. Factors such as the environment, hormones, and genetic and epigenetic background can also affect sex expression. Therefore, in this review, we will discuss recent research on the sex of papaya, from reported genes to date, their biology, and sexing approaches using molecular markers and their advantages.

Genomic and Transcriptomic Insights into the Evolution and Divergence of MIKC-Type MADS-Box Genes in Carica papaya

MIKC-type MADS-box genes, also known as type II genes, play a crucial role in regulating the formation of floral organs and reproductive development in plants. However, the genome-wide identification and characterization of type II genes as well as a transcriptomic survey of their potential roles in Carica papaya remain unresolved. Here, we identified and characterized 24 type II genes in the C. papaya genome, and investigated their evolutional scenario and potential roles with a widespread expression profile. The type II genes were divided into thirteen subclades, and gene loss events likely occurred in papaya, as evidenced by the contracted member size of most subclades. Gene duplication mainly contributed to MIKC-type gene formation in papaya, and the duplicated gene pairs displayed prevalent expression divergence, implying the evolutionary significance of gene duplication in shaping the diversity of type II genes in papaya. A large-scale transcriptome analysis of 152 samples indicated that different subclasses of these genes showed distinct expression patterns in various tissues, biotic stress response, and abiotic stress response, reflecting their divergent functions. The hub-network of male and female flowers and qRT-PCR suggested that TT16-3 and AGL8 participated in male flower development and seed germination. Overall, this study provides valuable insights into the evolution and functions of MIKC-type genes in C. papaya.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Genome-Wide Analysis of WRKY and NAC Transcription Factors in Carica papaya L. and Their Possible Role in the Loss of Drought Tolerance by Recent Cultivars through the Domestication of Their Wild Ancestors

A genome-wide analysis for two families of key transcription factors (TF; WRKY and NAC) involved in drought response revealed 46 WRKY and 66 NAC members of the Carica papaya genome. A phylogenetic analysis grouped the CpWRKY proteins into three groups (I, II a, b, c, d, e and III), while the CpNAC proteins were clustered into 15 groups. The conserved domains, chromosomal localization and promoter cis-acting elements were also analyzed. In addition, from a previous transcriptome study of two contrasting genotypes in response to 14 days of water deficit stress (WDS), we found that 29 of the 46 CpWRKYs genes and 25 of the 66 CpNACs genes were differentially expressed in response to the WDS. In the present paper, the native wild genotype (WG) (collected in its center of origin) consistently showed a higher expression (transcripts per million; TPM and fold change; FC) than the commercial genotype (CG) in almost all the members of the CpWRKY and CpNAC gene families. To corroborate this, we selected CpWRKY50 and CpNAC83.1 for further evaluation by RT-qPCR. Consistently, the WG showed higher relative expression levels (REL) after 14 days of WDS than the CG, in both the leaves and roots. The results suggest that the CpWRKY and CpNAC TF families are important for drought tolerance in this species. The results may also suggest that, during the domestication process, the ability of the native (wild) C. papaya genotypes to respond to drought (including the overexpression of the CpWRKY and CpNAC genes) was somehow reduced in the current commercial genotypes.

Strategies for delivery of CRISPR/Cas-mediated genome editing to obtain edited plants directly without transgene integration

Increased understanding of plant genetics and the development of powerful and easier-to-use gene editing tools over the past century have revolutionized humankind's ability to deliver precise genotypes in crops. Plant transformation techniques are well developed for making transgenic varieties in certain crops and model organisms, yet reagent delivery and plant regeneration remain key bottlenecks to applying the technology of gene editing to most crops. Typical plant transformation protocols to produce transgenic, genetically modified (GM) varieties rely on transgenes, chemical selection, and tissue culture. Typical protocols to make gene edited (GE) varieties also use transgenes, even though these may be undesirable in the final crop product. In some crops, the transgenes are routinely segregated away during meiosis by performing crosses, and thus only a minor concern. In other crops, particularly those propagated vegetatively, complex hybrids, or crops with long generation times, such crosses are impractical or impossible. This review highlights diverse strategies to deliver CRISPR/Cas gene editing reagents to regenerable plant cells and to recover edited plants without unwanted integration of transgenes. Some examples include delivering DNA-free gene editing reagents such as ribonucleoproteins or mRNA, relying on reagent expression from non-integrated DNA, using novel delivery mechanisms such as viruses or nanoparticles, using unconventional selection methods to avoid integration of transgenes, and/or avoiding tissue culture altogether. These methods are advancing rapidly and already enabling crop scientists to make use of the precision of CRISPR gene editing tools.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Engineered biocontainable RNA virus vectors for non-transgenic genome editing across crop species and genotypes

CRISPR/Cas genome-editing tools provide unprecedented opportunities for basic plant biology research and crop breeding. However, the lack of robust delivery methods has limited the widespread adoption of these revolutionary technologies in plant science. Here, we report an efficient, non-transgenic CRISPR/Cas delivery platform based on the engineered tomato spotted wilt virus (TSWV), an RNA virus with a host range of over 1000 plant species. We eliminated viral elements essential for insect transmission to liberate genome space for accommodating large genetic cargoes without sacrificing the ability to infect plant hosts. The resulting non-insect-transmissible viral vectors enabled effective and stable in planta delivery of Cas12a and Cas9 nucleases as well as adenine and cytosine base editors. In systemically infected plant tissues, the deconstructed TSWV-derived vectors induced efficient somatic gene mutations and base conversions in multiple crop species with little genotype dependency. Plants with heritable, bi-allelic mutations could be readily regenerated by culturing the virus-infected tissues in vitro without antibiotic selection. Moreover, we showed that antiviral treatment with ribavirin during tissue culture cleared the viral vectors in 100% of regenerated plants and further augmented the recovery of heritable mutations. Because many plants are recalcitrant to stable transformation, the viral delivery system developed in this work provides a promising tool to overcome gene delivery bottlenecks for genome editing in various crop species and elite varieties.

Comparative transcriptomics and genomic analyses reveal differential gene expression related to Colletotrichum brevisporum resistance in papaya (Carica papaya L.)

Colletotrichum brevisporum is an important causal pathogen of anthracnose that seriously affects the fruit quality and yield of papaya (Carica papaya L.). Although many genes and biological processes involved in anthracnose resistance have been reported in other species, the molecular mechanisms involved in the response or resistance to anthracnose in post-harvest papaya fruits remain unclear. In this study, we compared transcriptome changes in the post-harvest fruits of the anthracnose-susceptible papaya cultivar Y61 and the anthracnose-resistant cultivar G20 following C. brevisporum inoculation. More differentially expressed genes (DEGs) and differentially expressed long non-coding RNAs (DElnRNAs) were identified in G20 than in Y61, especially at 24 h post-inoculation (hpi), suggesting a prompt activation of defense responses in G20 in the first 24 h after C. brevisporum inoculation. These DEGs were mainly enriched in plant-pathogen interaction, phenylpropanoid biosynthesis/metabolism, and peroxisome and flavonoid biosynthesis pathways in both cultivars. However, in the first 24 hpi, the number of DEGs related to anthracnose resistance was greater in G20 than in Y61, and changes in their expression levels were faster in G20 than in Y61. We also identified a candidate anthracnose-resistant gene cluster, which consisted of 12 genes, 11 in G20 and Y61, in response to C. brevisporum inoculation. Moreover, 529 resistance gene analogs were identified in papaya genome, most of which responded to C. brevisporum inoculation and were genetically different between papaya cultivars and wild-type populations. The total expression dose of the resistance gene analogs may help papaya resist C. brevisporum infection. This study revealed the mechanisms underlying different anthracnose resistance between the anthracnose-resistant and anthracnose-susceptible cultivars based on gene expression, and identified some potential anthracnose resistance-related candidate genes/major regulatory factors. Our findings provided potential targets for developing novel genetic strategies to overcome anthracnose in papaya.

Genomic Analysis of LEA Genes in Carica papaya and Insight into Lineage-Specific Family Evolution in Brassicales

Late embryogenesis abundant (LEA) proteins comprise a diverse superfamily involved in plant development and stress responses. This study presents a first genome-wide analysis of LEA genes in papaya (Carica papaya L., Caricaceae), an economically important tree fruit crop widely cultivated in the tropics and subtropics. A total of 28 members were identified from the papaya genome, which belong to eight families with defined Pfam domains, i.e., LEA_1 (3), LEA_2 (4), LEA_3 (5), LEA_4 (5), LEA_5 (2), LEA_6 (2), DHN (4), and SMP (3). The family numbers are comparable to those present in Ricinus communis (Euphorbiaceae, 28) and Moringa oleifera (Moringaceae, 29), but relatively less than that found in Moringa oleifera (Cleomaceae, 39) and Arabidopsis thaliana (Brassicaceae, 51), implying lineage-specific evolution in Brassicales. Indeed, best-reciprocal-hit-based sequence comparison and synteny analysis revealed the presence of 29 orthogroups, and significant gene expansion in Tarenaya and Arabidopsis was mainly contributed by whole-genome duplications that occurred sometime after their split with the papaya. Though a role of transposed duplication was also observed, tandem duplication was shown to be a key contributor in gene expansion of most species examined. Further comparative analyses of exon-intron structures and protein motifs supported fast evolution of this special superfamily, especially in Arabidopsis. Transcriptional profiling revealed diverse expression patterns of CpLEA genes over various tissues and different stages of developmental fruit. Moreover, the transcript level of most genes appeared to be significantly regulated by drought, cold, and salt stresses, corresponding to the presence of cis-acting elements associated with stress response in their promoter regions. These findings not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further functional analysis of LEA genes in papaya.