Deciphering the evolutionary development of the "Chinese lantern" within Solanaceae.

MAIN CONCLUSION: The key genetic variation underlying the evo-devo of ICS in Solanaceae may be further pinpointed using an integrated strategy of forward and reverse genetics studies under the framework of phylogeny. The calyx of Physalis remains persistent throughout fruit development. Post-flowering, the fruiting calyx is inflated rapidly to encapsulate the berry, giving rise to a "Chinese lantern" structure called inflated calyx syndrome (ICS). It is unclear how this novelty arises. Over the past 2 decades, the role of MADS-box genes in the evolutionary development (evo-devo) of ICS has mainly been investigated within Solanaceae. In this review, we analyze the main achievements, challenges, and new progress. ICS acts as a source for fruit development, provides a microenvironment to protect fruit development, and assists in long-distance fruit dispersal. ICS is a typical post-floral trait, and the onset of its development is triggered by specific developmental signals that coincide with fertilization. These signals can be replaced by exogenous gibberellin and cytokinin application. MPF2-like heterotopic expression and MBP21-like loss have been proposed to be two essential evolutionary events for ICS origin, and manipulating the related MADS-box genes has been shown to affect the ICS size, sepal organ identity, and/or male fertility, but not completely disrupt ICS. Therefore, the core genes or key links in the ICS biosynthesis pathways may have undergone secondary mutations during evolution, or they have not yet been pinpointed. Recently, we have made some encouraging progress in acquiring lantern mutants in Physalis floridana. In addition to technological innovation, we propose an integrated strategy to further analyze the evo-devo mechanisms of ICS in Solanaceae using forward and reverse genetics studies under the framework of phylogeny.

Unraveling the Genetic Control of Pigment Accumulation in Physalis Fruits.

Physalis pubescens and Physalis alkekengi, members of the Physalis genus, are valued for their delicious and medicinal fruits as well as their different ripened fruit colors-golden for P. pubescens and scarlet for P. alkekengi. This study aimed to elucidate the pigment composition and genetic mechanisms during fruit maturation in these species. Fruit samples were collected at four development stages, analyzed using spectrophotometry and high-performance liquid chromatography (HPLC), and complemented with transcriptome sequencing to assess gene expression related to pigment biosynthesis. ß-carotene was identified as the dominant pigment in P. pubescens, contrasting with P. alkekengi, which contained both lycopene and ß-carotene. The carotenoid biosynthesis pathway was central to fruit pigmentation in both species. Key genes pf02G043370 and pf06G178980 in P. pubescens, and TRINITY_DN20150_c1_g3, TRINITY_DN10183_c0_g1, and TRINITY_DN23805_c0_g3 in P. alkekengi were associated with carotenoid production. Notably, the MYB-related and bHLH transcription factors (TFs) regulated zeta-carotene isomerase and ß-hydroxylase activities in P. pubescens with the MYB-related TF showing dual regulatory roles. In P. alkekengi, six TF families-bHLH, HSF, WRKY, M-type MADS, AP2, and NAC-were implicated in controlling carotenoid synthesis enzymes. Our findings highlight the intricate regulatory network governing pigmentation and provide insights into Physalis germplasm's genetic improvement and conservation.

First draft reference genome and annotation of the alternative oil species Physaria fendleri.

In the wake of increasing demand for renewable energy sources, plant-based sources including alternative oilseeds have come to the forefront of interest. Hydroxy fatty acids (HFAs), produced in a few oilseed species, are important chemical feed stocks for industrial applications. An integrated approach was taken to assemble the first draft genome of the alternative HFA producer Physaria fendleri (n = 6), an outcrossing species with high heterozygosity. Both de novo transcriptome assemblies and genome assemblies were produced with public and generated sequencing reads. Resulting intermediate assemblies were then scaffolded and patched with multiple data sources, followed by super-scaffolding onto a masked genome of Camelina laxa (n = 6). Despite a current lack of available resources for the physical mapping of genomic scaffolds of P. fendleri, topography of the genome with respect to repeat and gene content was preserved at the scaffold level and not significantly lost via super-scaffolding. Read representation, gene and genome completion statistics, and annotation results illustrated the creation of a functional draft genome and a tool for future research on alternative oil species.

Characterization and expression analysis of a thaumatin-like protein PpTLP1 from ground cherry Physalis pubescens.

Ground cherry, Physalis pubescens, is mainly cultivated as a fruit worldwide and popularly used as a food supplement and traditional Chinese medicine. Plants are challenged by external environmental stress and can initiate resistance to the stress through the regulation of pathogenesis-related (PR) proteins. Among PR proteins, PR-5, a thaumatin-like protein (TLP), was identified in many plants and found to be able to enhance stress resistance. However, PR-5 in ground cherry is not characterized and its expression is yet to be understood. In this study, a PR-5 protein PpTLP1 in P. pubescens was firstly identified. Analysis of the amino acid sequences revealed that PpTLP1 was highly similar to PR-NP24 identified in tomato with a difference in only one amino acid. Expression analysis indicated that the PpTLP1 gene was highly expressed in leaf while the PpTLP1 protein was tissue-specifically accumulated in cherry exocarp. Furthermore, the down-regulation of PpTLP1 in ground cherry was induced by NaCl treatment while the up-regulation was promoted by the infection of Sclerotinia sclerotiorum and Botrytis cinerea. This study will provide a new plant resource containing a TLP in Physalis genus and a novel insight for the improvement of postharvest management of ground cherry and other Solanaceae plants.

Identification of P450 Candidates Associated with the Biosynthesis of Physalin-Class Compounds in Physalis angulata.

The Physalis genus has long been used as traditional medicine in the treatment of various diseases. Physalins, the characteristic class of compounds in this genus, are major bioactive constituents. To date, the biogenesis of physalins remains largely unknown, except for the recently established knowledge that 24-methyldesmosterol is a precursor of physalin. To identify the genes encoding P450s that are putatively involved in converting 24-methyldesmosterol to physalins, a total of 306 P450-encoding unigenes were retrieved from our recently constructed P. angulata transcriptome. Extensive phylogenetic analysis proposed 21 P450s that might participate in physalin biosynthesis. To validate the candidates, we developed a virus-induced gene silencing (VIGS) system for P. angulata, and four P450 candidates were selected for the VIGS experiments. The reduction in the transcripts of the four P450 candidates by VIGS all led to decreased levels of physalin-class compounds in the P. angulata leaves. Thus, this study provides a number of P450 candidates that are likely associated with the biosynthesis of physalin-class compounds, forming a strong basis to reveal the unknown physalin biosynthetic pathway in the future.

Establishing Physalis as a Solanaceae model system enables genetic reevaluation of the inflated calyx syndrome.

The highly diverse Solanaceae family contains several widely studied models and crop species. Fully exploring, appreciating, and exploiting this diversity requires additional model systems. Particularly promising are orphan fruit crops in the genus Physalis, which occupy a key evolutionary position in the Solanaceae and capture understudied variation in traits such as inflorescence complexity, fruit ripening and metabolites, disease and insect resistance, self-compatibility, and most notable, the striking inflated calyx syndrome (ICS), an evolutionary novelty found across angiosperms where sepals grow exceptionally large to encapsulate fruits in a protective husk. We recently developed transformation and genome editing in Physalis grisea (groundcherry). However, to systematically explore and unlock the potential of this and related Physalis as genetic systems, high-quality genome assemblies are needed. Here, we present chromosome-scale references for P. grisea and its close relative Physalis pruinosa and use these resources to study natural and engineered variations in floral traits. We first rapidly identified a natural structural variant in a bHLH gene that causes petal color variation. Further, and against expectations, we found that CRISPR-Cas9-targeted mutagenesis of 11 MADS-box genes, including purported essential regulators of ICS, had no effect on inflation. In a forward genetics screen, we identified huskless, which lacks ICS due to mutation of an AP2-like gene that causes sepals and petals to merge into a single whorl of mixed identity. These resources and findings elevate Physalis to a new Solanaceae model system and establish a paradigm in the search for factors driving ICS.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae.

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Natural infection of cape gooseberry (Physalis peruviana) by the potyvirus Tamarillo leaf malformation virus (TaLMV).

Tamarillo leaf malformation virus (TaLMV) is a potyvirus first discovered in cape gooseberry fields in Eastern, and South-western Antioquia. This virus is responsible for a very damaging disease that has resulted in significant reductions in yields and cultivated area for this crop in Colombia. Tamarillo is frequently co-cultivated with other solanaceous plants but no evidence for cross-pathogenicity of TaLMV has been found until now. In this work, we report a natural infection of cape gooseberry (Physalis peruviana L.) by TaLMV. Infection by TaLMV was detected by RNAseq screening of cape gooseberry fields and confirmed by RT-qPCR and Sanger sequencing. The sequenced genome is 99.3% identical to previously sequenced TaLMV isolates, and evidence suggests that it can accumulate at high loads in this new reported host. RT-qPCR analysis indicates that TaLMV is already widely distributed, can naturally infect other solanaceous hosts and may become an emerging threat to the cape gooseberry agroindustry, the second most important exotic fruit export in Colombia. Keywords: high-throughput sequencing; plant virology; Potyviridae; RT-qPCR; Solanaceae.

A lineage-specific arginine in POS1 is required for fruit size control in Physaleae (Solanaceae) via gene co-option.

Solanaceae have important economic value mainly due to their edible fruits. Physalis organ size 1/cytokinin response factor 3 (POS1/CRF3), a unique gene in Solanaceae, is involved in fruit size variation in Physalis but not in Solanum. However, the underlying mechanisms remain elusive. Here, we found that POS1/CRF3 was likely created via the fusion of CRF7 and CRF8 duplicates. Multiple genetic manipulations revealed that only POS1 and Capsicum POS1 (CaPOS1) functioned in fruit size control via the positive regulation of cell expansion. Comparative studies in a phylogenetic framework showed the directional enhancement of POS1-like expression in the flowers and fruits of Physaleae and the specific gain of certain interacting proteins associated with cell expansion by POS1 and CaPOS1. A lineage-specific single nucleotide polymorphism (SNP) caused the 68th amino acid histidine in the POS1 orthologs of non-Physaleae (Nicotiana and Solanum) to change to arginine in Physaleae (Physalis and Capsicum). Substituting the arginine in Physaleae POS1-like by histidine completely abolished their function in the fruits and the protein-protein interaction (PPI) with calreticulin-3. Transcriptomic comparison revealed the potential downstream pathways of POS1, including the brassinosteroid biosynthesis pathway. However, POS1-like may have functioned ancestrally in abiotic stress within Solanaceae. Our work demonstrated that heterometric expression and a SNP caused a single amino acid change to establish new PPIs, which contributed to the co-option of POS1 in multiple regulatory pathways to regulate cell expansion and thus fruit size in Physaleae. These results provide new insights into fruit morphological evolution and fruit yield control.

The draft chromosome-level genome assembly of tetraploid ground cherry (Prunus fruticosa Pall.) from long reads.

Cherries are stone fruits and belong to the economically important plant family of Rosaceae with worldwide cultivation of different species. The ground cherry, Prunus fruticosa Pall., is an ancestor of cultivated sour cherry, an important tetraploid cherry species. Here, we present a long read chromosome-level draft genome assembly and related plastid sequences using the Oxford Nanopore Technology PromethION platform and R10.3 pore type. We generated a final consensus genome sequence of 366 Mb comprising eight chromosomes. The N50 scaffold was ~44 Mb with the longest chromosome being 66.5 Mb. The chloroplast and mitochondrial genomes were 158,217 bp and 383,281 bp long, which is in accordance with previously published plastid sequences. This is the first report of the genome of ground cherry (P. fruticosa) sequenced by long read technology only. The datasets obtained from this study provide a foundation for future breeding, molecular and evolutionary analysis in Prunus studies.

Multiple and integrated functions of floral C-class MADS-box genes in flower and fruit development of Physalis floridana.

KEY MESSAGE: This work reveals potentially multiple and integrated roles in flower and fruit development of floral C-class MADS-box genes in Physalis. The Physalis fruit features a morphological novelty, the Chinese lantern. Floral C-class MADS-domain AGAMOUS-like (AG-like) proteins can interact with the identified regulators of this novel structure. However, the developmental role of the floral C-class genes is unknown in Physalis. Here, we characterized two AG-like genes from Physalis floridana, designated PFAG1 and PFAG2. The two paralogous genes shared around 61.0% of sequence identity and had similar expression domains, with different expression levels in the floral and berry development. However, the genes had distinct expression patterns in leaf and calyx development. Protein-protein interaction analyses revealed that PFAG1 and PFAG2 could commonly or specifically dimerize with certain floral MADS-domain proteins as well as non-MADS-domain proteins involved in various floral developmental processes. Gene downregulation analyses demonstrated that PFAG1 may repress PFAG2, but PFAG2 did not affect PFAG1. Downregulating PFAG1 led to incomplete floral homeotic variation in the stamens and carpels, and alteration of petal coloration pattern, while downregulating PFAG2 did not result in any floral homeotic variation. PFAG1 affected pollen maturation, while PFAG2 affected female fertility. However, simultaneously downregulating PFAG1 and PFAG2 caused loss of the complete C-function, indicating that the two PFAG genes interact to determine the identity and functionality of androecia and gynoecia organs. Their potential roles in regulating fruit size and the Chinese lantern are also discussed. Our results reveal functional divergence of floral C-class MADS-box genes in Physalis, demonstrating that they may play multiple and integrated roles in flower and fruit development.

Anther Culture of the Gametophytic Self-Incompatible Species Physalis ixocarpa Brot.

Here we present an optimized protocol for in vitro embryo formation and plant regeneration through anther culture of the Mexican husk tomato (Physalis ixocarpa Brot.). This protocol relies on the application of an anther thermal shock at a specific developmental stage prior to the in vitro culture, ensures embryo formation from anthers without callus formation, and allows spending less time to regenerate doubled haploid complete plants. This protocol has been used for different cultivars of Physalis ixocarpa (Chapingo, Rendidora, Puebla, Arandaz, Manzano, Tamazula, Salamanca, and Milpero), and also for two wild-type accessions, all of them cultivated in Mexico. Chapingo cultivar responded with the highest percentage of androgenesis on the embryo induction medium (EIM).

Physalis floridana CRABS CLAW mediates neofunctionalization of GLOBOSA genes in carpel development.

Floral B-function MADS-box genes, such as GLOBOSA (GLO), function in corolla and stamen organ identity specification. The functions of these genes outside these floral whorls are rarely reported. DOLL1 is a GLO gene controlling corolla and androecium organ identity. In this study we found that, in Physalis floridana double-layered-lantern 1 (doll1) mutant pollinated with wild-type pollen, fruit set was extremely low, indicating that doll1 females are dysfunctional. Stigma and style structure, stigma receptivity, pollen tube guidance, and embryo sac development were also impaired in doll1. P. floridana CRABS CLAW (PFCRC), predominantly expressed in carpels, was repressed in doll1 native carpels. Loss-of-function of PFCRC altered carpel meristem determinacy, carpel closure, and ovule number, and the resultant 'pistil' consisted of multiple spirally-arranged dorsiventral carpels occasionally with 1-2 naked ovules on the margin and trichomes at each mutated carpel tip, implying an alteration of carpel organ identity. Regulatory and genetic interactions between B-class MADS-box genes and PFCRC were revealed in a context-dependent manner in floral development. Our work reveals a new role for the B-function genes in carpel and ovule development via regulating PFCRC, providing a new understanding of genetic regulatory networks between MADS-domain and CRC transcription factors in mediating carpel organ specification, functionality, and origin.

Optimization of the genotyping-by-sequencing SNP calling for diversity analysis in cape gooseberry (Physalis peruviana L.) and related taxa.

A robust Genotyping-By-Sequencing (GBS) pipeline platform was examined to provide accurate discovery of Single Nucleotide Polymorphisms (SNPs) in a cape gooseberry (Physalis peruviana L.) and related taxa germplasm collection. A total of 176 accessions representing, wild, weedy, and commercial cultivars as well as related taxa from the Colombian germplasm bank and other world repositories were screened using GBS. The pipeline parameters mnLCov of 0.5 and a mnScov of 0.7, tomato and potato genomes, and cape gooseberry transcriptome for read alignments, were selected to better assess diversity and population structure in cape gooseberry and related taxa. A total of 7,425 SNPs, derived from P. peruviana common tags (unique 64 bp sequences shared between selected species), were used. Within P. peruviana, five subpopulations with a high genetic diversity and allele fixation (HE: 0.35 to 0.36 and FIS: -0.11 to -0.01, respectively) were detected. Conversely, low genetic differentiation (FST: 0.01 to 0.05) was also observed, indicating a high gene flow among subpopulations. These results contribute to the establishment of adequate conservation and breeding strategies for Cape gooseberry and closely related Physalis species.

Transcriptomic variation of the flower-fruit transition in Physalis and Solanum.

MAIN CONCLUSION: Gene expression variations in response to fertilization between Physalis and Solanum might play essential roles in species divergence and fruit evolution. Fertilization triggers variation in fruit development and morphology. The Chinese lantern, a morphological novelty derived from the calyx, is formed upon fertilization in Physalis but is not observed in Solanum. The underlying genetic variations are largely unknown. Here, we documented the developmental and morphological differences in the flower and fruit between Physalis floridana and Solanum pimpinellifolium and then evaluated both the transcript sequence variation and gene expression at the transcriptomic level at fertilization between the two species. In Physalis transcriptomic analysis, 468 unigenes were identified as differentially expressed genes (DEGs) that were strongly regulated by fertilization across 3 years. In comparison with tomato, 14,536 strict single-copy orthologous gene pairs were identified between P. floridana and S. pimpinellifolium in the flower-fruit transcriptome. Nine types of gene variations with specific GO-enriched patterns were identified, covering 58.82% orthologous gene pairs that were DEGs in either trend or dosage at the flower-fruit transition between the two species, which could adequately distinguish Solanum and Physalis, implying that differential gene expression at fertilization might play essential roles during the divergence and fruit evolution of Solanum-Physalis. Virus-induced gene silencing analyses revealed the developmental roles of some transcription factor genes in fertility, Chinese lantern development, and fruit weight control in Physalis. This study presents the first floral transcriptomic resource of Physalis, and reveals some candidate genetic variations accounting for the early fruit developmental evolution in Physalis in comparison to Solanum.

Complete chloroplast genomes of four Physalis species (Solanaceae): lights into genome structure, comparative analysis, and phylogenetic relationships.

BACKGROUND: Physalis L. is a genus of herbaceous plants of the family Solanaceae, which has important medicinal, edible, and ornamental values. The morphological characteristics of Physalis species are similar, and it is difficult to rapidly and accurately distinguish them based only on morphological characteristics. At present, the species classification and phylogeny of Physalis are still controversial. In this study, the complete chloroplast (cp) genomes of four Physalis species (Physalis angulata, P. alkekengi var. franchetii, P. minima and P. pubescens) were sequenced, and the first comprehensive cp genome analysis of Physalis was performed, which included the previously published cp genome sequence of Physalis peruviana. RESULTS: The Physalis cp genomes exhibited typical quadripartite and circular structures, and were relatively conserved in their structure and gene synteny. However, the Physalis cp genomes showed obvious variations at four regional boundaries, especially those of the inverted repeat and the large single-copy regions. The cp genomes' lengths ranged from 156,578 bp to 157,007 bp. A total of 114 different genes, 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes, were observed in four new sequenced Physalis cp genomes. Differences in repeat sequences and simple sequence repeats were detected among the Physalis cp genomes. Phylogenetic relationships among 36 species of 11 genera of Solanaceae based on their cp genomes placed Physalis in the middle and upper part of the phylogenetic tree, with a monophyletic evolution having a 100% bootstrap value. CONCLUSION: Our results enrich the data on the cp genomes of the genus Physalis. The availability of these cp genomes will provide abundant information for further species identification, increase the taxonomic and phylogenetic resolution of Physalis, and assist in the investigation and utilization of Physalis plants.

Transcriptome-wide identification of microRNAs and functional insights inferred from microRNA-target pairs in Physalis angulata L.

Physalis angulata L., a member of the family Solanaceae, is widely used as the folk medicine in various countries. Continuous research efforts are devoted to the discovery of the effective medicinal ingredients from Physalis angulata. However, due to the limited resources of genome and transcriptome sequencing data, only a few studies have been performed at the gene regulatory level. In this study, the transcriptomes of five organs (roots, stems, leaves, flowers and fruits) of Physalis angulata were reported. Based on the transcriptome assembly containing 196,117 unique transcripts, a total of 17,556 SSRs (simple sequence repeats) were identified, which could be useful RNA-based barcoding for discrimination of the plants closely relative to Physalis angulata. Additionally, 24 transcripts were discovered to be the potential microRNA (miRNA) precursors which encode a total of 31 distinct mature miRNAs. Some of these precursors showed organ-specific expression patterns. Target prediction revealed 116 miRNA-target pairs, involving 31 miRNAs and 83 target transcripts in Physalis angulata. Taken together, our results could serve as the data resource for in-depth studies on the molecular regulatory mechanisms related to the production of medicinal ingredients in Physalis angulata.

An Extension of the Kimura Two-Parameter Model to the Natural Evolutionary Process.

Accurate estimates of genetic difference are required for research in evolutionary biology. Here we extend the Kimura two-parameter (K2P) model by considering gaps (insertions and/or deletions) and introduce a new measure for estimating genetic difference between two nucleotide sequences in terms of nucleotide changes that have occurred during the evolutionary process. Using the nuclear ribosomal DNA internal transcribed spacer 2 region from the genus Physalis, we demonstrate that species identification and phylogenetic studies strongly depend on evolutionary models. It is especially noteworthy that the use of different models affects the degree of overlap between intraspecific and interspecific genetic differences. We observe that the percentage of interspecific sequence pairs with values less than the maximum intraspecific genetic difference is 43.2% for the K2P model which is calculated by removing gap sites across all sequences, 22.7% for the K2P model which is calculated by removing gap sites for sequence pairs, and 16.9% for our model which is calculated without removing gap sites. Additionally, the numbers of sequence pairs with interspecific genetic differences of zero are 50 for the K2P model and 29 for our model. The genetic difference measure based on the K2P model, compared to our model, overestimates 21 sequence pairs that are not originally identical. These results indicate the importance of estimating genetic differences under the model of sequence evolution that includes insertions and deletions in addition to substitutions.

Quantification of withaferin-A and withanolide-A in diploid (n = 12) and tetraploid cytotypes (n = 24) of "Rassbhary", Physalis angulata L.

During the present study an analytical method based on reverse-phase high performance liquid chromatography photodiode array detection method was developed for simultaneous determination of withaferin-A and withanolide-A in plant parts of two cytotypes (diploid n = 12 & tetraploid n = 24) of Physalis angulata. All the plant parts were extracted in different solvent solutions i.e., acidic [HCl] methanol (i.e., methanol containing 0.3% of HCl), methanol, n-hexane, chloroform. Both the compounds were comparatively analysed. The results revealed that tetraploid cytotype (n = 24) showed the higher composition of both the reference compounds. The method is simple, rapid and provides better resolution can be easily applied to the quantitative analyses of withanolides in plant matrices.

Exon junction complex (EJC) core genes play multiple developmental roles in Physalis floridana.

KEY MESSAGE: Molecular and functional characterization of four gene families of the Physalis exon junction complex (EJC) core improved our understanding of the evolution and function of EJC core genes in plants. The exon junction complex (EJC) plays significant roles in posttranscriptional regulation of genes in eukaryotes. However, its developmental roles in plants are poorly known. We characterized four EJC core genes from Physalis floridana that were named PFMAGO, PFY14, PFeIF4AIII and PFBTZ. They shared a similar phylogenetic topology and were expressed in all examined organs. PFMAGO, PFY14 and PFeIF4AIII were localized in both the nucleus and cytoplasm while PFBTZ was mainly localized in the cytoplasm. No protein homodimerization was observed, but they could form heterodimers excluding the PFY14-PFBTZ heterodimerization. Virus-induced gene silencing (VIGS) of PFMAGO or PFY14 aborted pollen development and resulted in low plant survival due to a leaf-blight-like phenotype in the shoot apex. Carpel functionality was also impaired in the PFY14 knockdowns, whereas pollen maturation was uniquely affected in PFBTZ-VIGS plants. Once PFeIF4AIII was strongly downregulated, plant survival was reduced via a decomposing root collar after flowering and Chinese lantern morphology was distorted. The expression of Physalis orthologous genes in the DYT1-TDF1-AMS-bHLH91 regulatory cascade that is associated with pollen maturation was significantly downregulated in PFMAGO-, PFY14- and PFBTZ-VIGS flowers. Intron-retention in the transcripts of P. floridana dysfunctional tapetum1 (PFDYT1) occurred in these mutated flowers. Additionally, the expression level of WRKY genes in defense-related pathways in the shoot apex of PFMAGO- or PFY14-VIGS plants and in the root collar of PFeIF4AIII-VIGS plants was significantly downregulated. Taken together, the Physalis EJC core genes play multiple roles including a conserved role in male fertility and newly discovered roles in Chinese lantern development, carpel functionality and defense-related processes. These data increase our understanding of the evolution and functions of EJC core genes in plants.