An evolutionary view of vanillylamine synthase pAMT, a key enzyme of capsaicinoid biosynthesis pathway in chili pepper

Pungent capsaicinoid is synthesized only in chili pepper (Capsicum spp.). The production of vanillylamine from vanillin is a unique reaction in the capsaicinoid biosynthesis pathway. Although putative aminotransferase (pAMT) has been isolated as the vanillylamine synthase gene, it is unclear how Capsicum acquired pAMT. Here, we present a phylogenetic overview of pAMT and its homologs. The Capsicum genome contained 5 homologs, including pAMT, CaGABA-T1, CaGABA-T3, and two pseudogenes. Phylogenetic analysis indicated that pAMT is a member of the Solanaceae cytoplasmic GABA-Ts. Comparative genome analysis found that multiple copies of GABA-T exist in a specific Solanaceae genomic region, and the cytoplasmic GABA-Ts other than pAMT are located in the region. The cytoplasmic GABA-T was phylogenetically close to pseudo-GABA-T harboring a plastid transit peptide (pseudo-GABA-T3). This suggested that Solanaceae cytoplasmic GABA-Ts occurred via duplication of a chloroplastic GABA-T ancestor and subsequent loss of the plastid transit signal. The cytoplasmic GABA-T may have been translocated from the specific Solanaceae genomic region during Capsicum divergence, resulting in the current pAMT locus. A recombinant protein assay demonstrated that pAMT had higher vanillylamine synthase activity than those of other plant GABA-Ts. pAMT was expressed exclusively in the placental septum of mature green fruit, whereas tomato orthologs SlGABA-T2/4 exhibit a ubiquitous expression pattern in plants. These findings suggested that both the increased catalytic efficiency and transcriptional changes in pAMT may have contributed to establish vanillylamine synthesis in the capsaicinoid biosynthesis pathway. This study provides insights into the establishment of pungency in the evolution of chili peppers.

The phylogeographic journey of a plant species from lowland to highlands during the Pleistocene

Phylogeographic history refers to how species evolve and diversify in response to historical, ecological, and demographic factors. The climate fluctuation during the Pleistocene period marked a crucial time in shaping many species' distribution and genetic structure, particularly those from southern South American grasslands. This work investigated the phylogeographic history of a highland grassland, Petunia altiplana T. Ando & Hashim. (Solanaceae), its diversity, and geographic distribution using a population genomic approach based on RAD-seq data. Our results indicated that, during the Pleistocene, when the grasslands expanded to highlands, the lowland populations of P. altiplana reached the higher open fields, enlarging their geographic distribution. We found that the P. altiplana genetic diversity followed the geographic division into eastern (E) and western (WE) population groups, with a subtle division in the E group regarding the Pelotas River headwater. The results also showed that isolation by distance was the main divergence pattern, with elevation playing a pivotal role in shaping WE and E groups. Our findings indicated that lowland-adapted populations quickly colonized highlands during the late Pleistocene.

Chromosome level genome assembly of endangered medicinal plant Anisodus tanguticus

Anisodus tanguticus is a medicinal herb that belongs to the Anisodus genus of the Solanaceae family. This endangered herb is mainly distributed in Qinghai-Tibet Plateau. In this study, we combined the Illumina short-read, Nanopore long-read and high-throughput chromosome conformation capture (Hi-C) sequencing technologies to de novo assemble the A. tanguticus genome. A high-quality chromosomal-level genome assembly was obtained with a genome size of 1.26 Gb and a contig N50 of 25.07 Mb. Of the draft genome sequences, 97.47% were anchored to 24 pseudochromosomes with a scaffold N50 of 51.28 Mb. In addition, 842.14 Mb of transposable elements occupying 66.70% of the genome assembly were identified and 44,252 protein-coding genes were predicted. The genome assembly of A. tanguticus will provide genetic repertoire to understand the adaptation strategy of Anisodus species in the plateau, which will further promote the conservation of endangered A. tanguticus resources.

A chromosome-level genome assembly of anesthetic drug-producing Anisodus acutangulus provides insights into its evolution and the biosynthesis of tropane alkaloids

Tropane alkaloids (TAs), which are anticholinergic agents, are an essential class of natural compounds, and there is a growing demand for TAs with anesthetic, analgesic, and spasmolytic effects. Anisodus acutangulus (Solanaceae) is a TA-producing plant that was used as an anesthetic in ancient China. In this study, we assembled a high-quality, chromosome-scale genome of A. acutangulus with a contig N50 of 7.4 Mb. A recent whole-genome duplication occurred in A. acutangulus after its divergence from other Solanaceae species, which resulted in the duplication of ADC1 and UGT genes involved in TA biosynthesis. The catalytic activities of H6H enzymes were determined for three Solanaceae plants. On the basis of evolution and co-expressed genes, AaWRKY11 was selected for further analyses, which revealed that its encoded transcription factor promotes TA biosynthesis by activating AaH6H1 expression. These findings provide useful insights into genome evolution related to TA biosynthesis and have potential implications for genetic manipulation of TA-producing plants.

Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. SI is controlled by a single S-locus with multiple haplotypes (S-haplotypes). When the pistil and pollen share the same S-haplotype, the pollen is recognized as self and rejected by the pistil. This review introduces our research on Brassicaceae and Solanaceae SI systems to identify the S-determinants encoded at the S-locus and uncover the mechanisms of self/nonself-discrimination and pollen rejection. The recognition mechanisms of SI systems differ between these families. A self-recognition system is adopted by Brassicaceae, whereas a collaborative nonself-recognition system is used by Solanaceae. Work by our group and subsequent studies indicate that plants have evolved diverse SI systems.

Multiple independent losses of the biosynthetic pathway for two tropane alkaloids in the Solanaceae family

Hyoscyamine and scopolamine (HS), two valuable tropane alkaloids of significant medicinal importance, are found in multiple distantly related lineages within the Solanaceae family. Here we sequence the genomes of three representative species that produce HS from these lineages, and one species that does not produce HS. Our analysis reveals a shared biosynthetic pathway responsible for HS production in the three HS-producing species. We observe a high level of gene collinearity related to HS synthesis across the family in both types of species. By introducing gain-of-function and loss-of-function mutations at key sites, we confirm the reduced/lost or re-activated functions of critical genes involved in HS synthesis in both types of species, respectively. These findings indicate independent and repeated losses of the HS biosynthesis pathway since its origin in the ancestral lineage. Our results hold promise for potential future applications in the artificial engineering of HS biosynthesis in Solanaceae crops.

Abundance of distal repetitive DNA sequences in Capsicum L. (Solanaceae) chromosomes

Chili peppers (Solanaceae family) have great commercial value. They are commercialized in natura and used as spices and for ornamental and medicinal purposes. Although three whole genomes have been published, limited information about satellite DNA sequences, their composition, and genomic distribution has been provided. Here, we exploited the noncoding repetitive fraction, represented by satellite sequences, that tends to accumulate in blocks along chromosomes, especially near the chromosome ends of peppers. Two satellite DNA sequences were identified (CDR-1 and CDR-2), characterized and mapped in silico in three Capsicum genomes (C. annuum, C. chinense, and C. baccatum) using data from the published high-coverage sequencing and repeats finding bioinformatic tools. Localization using FISH in the chromosomes of these species and in two others (C. frutescens and C. chacoense), totaling five species, showed signals adjacent to the rDNA sites. A sequence comparison with existing Solanaceae repeats showed that CDR-1 and CDR-2 have different origins but without homology to rDNA sequences. Satellites occupied subterminal chromosomal regions, sometimes collocated with or adjacent to 35S rDNA sequences. Our results expand knowledge about the diversity of subterminal regions of Capsicum chromosomes, showing different amounts and distributions within and between karyotypes. In addition, these sequences may be useful for future phylogenetic studies.

Genome-Wide Analysis of MBF1 Family Genes in Five Solanaceous Plants and Functional Analysis of SlER24 in Salt Stress

Multiprotein bridging factor 1 (MBF1) is an ancient family of transcription coactivators that play a crucial role in the response of plants to abiotic stress. In this study, we analyzed the genomic data of five Solanaceae plants and identified a total of 21 MBF1 genes. The expansion of MBF1a and MBF1b subfamilies was attributed to whole-genome duplication (WGD), and the expansion of the MBF1c subfamily occurred through transposed duplication (TRD). Collinearity analysis within Solanaceae species revealed collinearity between members of the MBF1a and MBF1b subfamilies, whereas the MBF1c subfamily showed relative independence. The gene expression of SlER24 was induced by sodium chloride (NaCl), polyethylene glycol (PEG), ABA (abscisic acid), and ethrel treatments, with the highest expression observed under NaCl treatment. The overexpression of SlER24 significantly enhanced the salt tolerance of tomato, and the functional deficiency of SlER24 decreased the tolerance of tomato to salt stress. SlER24 enhanced antioxidant enzyme activity to reduce the accumulation of reactive oxygen species (ROS) and alleviated plasma membrane damage under salt stress. SlER24 upregulated the expression levels of salt stress-related genes to enhance salt tolerance in tomato. In conclusion, this study provides basic information for the study of the MBF1 family of Solanaceae under abiotic stress, as well as a reference for the study of other plants.

Anatomical and biophysical basis for graft incompatibility within the Solanaceae

Interspecies grafting is an economically relevant technique that allows beneficial shoot and root combinations from separate species to be combined. One hypothesis for the basis of graft compatibility revolves around taxonomic relatedness. To test how phylogenetic distance affects interspecific graft compatibility within the economically important Solanaceae subfamily, Solanoideae, we characterized the anatomical and biophysical integrity of graft junctions between four species: tomato (Solanum lycopersicum), eggplant (Solanum melongena), pepper (Capsicum annuum), and groundcherry (Physalis pubescens). We analyzed the survival, growth, integrity, and cellular composition of the graft junctions. Utilizing various techniques, we were able to quantitatively assess compatibility among the interspecific grafts. Even though most of our graft combinations could survive, we show that only intrageneric combinations between tomato and eggplant are compatible. Unlike incompatible grafts, the formation of substantial vascular reconnections between tomato and eggplant in the intrageneric heterografts likely contributed to biophysically stable grafts. Furthermore, we identified 10 graft combinations that show delayed incompatibility, providing a useful system to pursue deeper work into graft compatibility. This work provides new evidence that graft compatibility may be limited to intrageneric combinations within the Solanoideae subfamily. Further research amongst additional Solanaceous species can be used to test the extent to which our hypothesis applies to this family.

Genomic and structural basis for evolution of tropane alkaloid biosynthesis

The tropane alkaloids (TAs) cocaine and hyoscyamine have been used medicinally for thousands of years. To understand the evolutionary origins and trajectories of serial biosynthetic enzymes of TAs and especially the characteristic tropane skeletons, we generated the chromosome-level genome assemblies of cocaine-producing Erythroxylum novogranatense (Erythroxylaceae, rosids clade) and hyoscyamine-producing Anisodus acutangulus (Solanaceae, asterids clade). Comparative genomic and phylogenetic analysis suggested that the lack of spermidine synthase/N-methyltransferase (EnSPMT1) in ancestral asterids species contributed to the divergence of polyamine (spermidine or putrescine) methylation in cocaine and hyoscyamine biosynthesis. Molecular docking analysis and key site mutation experiments suggested that ecgonone synthases CYP81AN15 and CYP82M3 adopt different active-site architectures to biosynthesize the same product ecgonone from the same substrate in Erythroxylaceae and Solanaceae. Further synteny analysis showed different evolutionary origins and trajectories of CYP81AN15 and CYP82M3, particularly the emergence of CYP81AN15 through the neofunctionalization of ancient tandem duplication genes. The combination of structural biology and comparative genomic analysis revealed that ecgonone methyltransferase, which is responsible for the biosynthesis of characteristic 2-substituted carboxymethyl group in cocaine, evolved from the tandem copies of salicylic acid methyltransferase by the mutations of critical E216 and S153 residues. Overall, we provided strong evidence for the independent origins of serial TA biosynthetic enzymes on the genomic and structural level, underlying the chemotypic convergence of TAs in phylogenetically distant species.

The SlDOG1 Affect Biosynthesis of Steroidal Glycoalkaloids by Regulating GAME Expression in Tomato

Steroidal alkaloids (SAs) and steroidal glycoalkaloids (SGAs) are common constituents of plant species belonging to the Solanaceae family. However, the molecular mechanism regulating the formation of SAs and SGAs remains unknown. Here, genome-wide association mapping was used to elucidate SA and SGA regulation in tomatoes: a SlGAME5-like glycosyltransferase (Solyc10g085240) and the transcription factor SlDOG1 (Solyc10g085210) were significantly associated with steroidal alkaloid composition. In this study, it was found that rSlGAME5-like can catalyze a variety of substrates for glycosidation and can catalyze SA and flavonol pathways to form O-glucoside and O-galactoside in vitro. The overexpression of SlGAME5-like promoted α-tomatine, hydroxytomatine, and flavonol glycoside accumulation in tomatoes. Furthermore, assessments of natural variation combined with functional analyses identified SlDOG1 as a major determinant of tomato SGA content, which also promoted SA and SGA accumulation via the regulation of GAME gene expression. This study provides new insights into the regulatory mechanisms underlying SGA production in tomatoes.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae)

BACKGROUND AND AIMS

The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.

METHODS

We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.

KEY RESULTS

RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.

CONCLUSIONS

The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

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.

Characterization of the Evolutionary Pressure on Anisodus tanguticus Maxim. with Complete Chloroplast Genome Sequence

Anisodus tanguticus Maxim. (Solanaceae), a traditional endangered Tibetan herb, is endemic to the Qinghai-Tibet Plateau. Here, we report the de novo assembled chloroplast (cp) genome sequences of A. tanguticus (155,765 bp). The cp contains a pair of inverted repeated (IRa and IRb) regions of 25,881 bp that are separated by a large single copy (LSC) region (86,516 bp) and a small single copy SSC (17,487 bp) region. A total of 132 functional genes were annotated in the cp genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Moreover, 199 simple sequence repeats (SSR) and 65 repeat structures were detected. Comparative plastome analyses revealed a conserved gene order and high similarity of protein-coding sequences. The A. tanguticus cp genome exhibits contraction and expansion, which differs from Przewalskia tangutica and other related Solanaceae species. We identified 30 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which may be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Anisodus and family Solanaceae. Analysis of the Ka/Ks ratios of the Hyoscyameae tribe revealed significant positive selection exerted on the cemA, rpoC2, and clpP genes, which suggests that protein metabolism may be an important strategy for A. tanguticus and other species in Hyoscyameae in adapting to the adverse environment on the Qinghai-Tibetan Plateau. Phylogenetic analysis revealed that A. tanguticus clustered closer with Hyoscyamus niger than P. tangutica. Our results provide reliable genetic information for future exploration of the taxonomy and phylogenetic evolution of the Hyoscyameae tribe and related species.

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 evolutionary history of small RNAs in Solanaceae

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

Evolution and development of inflorescences and floral symmetry in Solanaceae

PREMISE

The inflorescences of Solanaceae are unique and complex, which has led to long-standing disputes over floral symmetry mainly due to different interpretations of the cyme-like inflorescence structure. The main disagreements have been over how the phyllomes associated with the flower were arranged relative to the inflorescence axis especially during early flower initiation.

METHODS

Here we investigated the evolution of inflorescences in Solanaceae by analyzing inflorescence structure in the context of phylogeny using ancestral state reconstruction (ASR) to determine the evolutionary transitions between loosely arranged and tightly clustered inflorescences and between monochasial-like and dichasial-like cymes. We also reconstructed two- and three-dimensional models for 12 solanaceous species that represent both inflorescence and phylogenetic diversity in the family.

RESULTS

Our results indicate that the most recent common ancestor of Solanaceae had a loosely arranged and monochasial-like cyme, while tightly clustered inflorescences and dichasial-like cymes were derived. Compared to the known process of scorpioid cyme evolution, Solanaceae achieved their scorpioid cyme-like inflorescences through a previously undescribed way. Along the pedicel, the two flower-preceding prophylls are not in the typical transverse position of dicotyledonous plants; they frequently have axillary buds, and the main inflorescence axis continues in a sympodial fashion. As a result, the plane of symmetry of the flower is 36° from the median, and the inflorescence axis and the two flower-preceding prophylls are symmetrically located along that plane.

CONCLUSIONS

A better understanding of the morphological evolution of solanaceous inflorescence structure helped clarify the floral symmetry of Solanaceae.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

CRISPR gene editing of major domestication traits accelerating breeding for Solanaceae crops improvement

Domestication traits particularly fruit size and plant architecture and flowering are critical in transforming a progenitor's wild stature into a super improved plant. The latest advancements in the CRISPR system, as well as its rapid adoption, are speeding up plant breeding. Solanaceae has a varied range of important crops, with a few model crops, such as tomato and, more recently, groundcherry, serving as a foundation for developing molecular techniques, genome editing tools, and establishing standards for other crops. Domestication traits in agricultural plants are quantified and widely adopted under modern plant breeding to improve small-fruited and bushy crop species like goji berry. The molecular mechanisms of the FW2.2, FW3.2, FW11.3, FAS/CLV3, LC/WUS, SP, SP5G, and CRISPR genome editing technology have been described in detail here. Furthermore, special focus has been placed on CRISPR gene editing achievements for revolutionizing Solanaceae breeding and changing the overall crop landscape. This review seeks to provide a thorough overview of the CRISPR technique's ongoing advancements, particularly in Solanaceae, in terms of domesticated features, future prospects, and regulatory risks. We believe that this vigorous discussion will lead to a broader understanding of CRISPR gene editing as a tool for achieving key breeding goals in other Solanaceae minor crops with significant industrial value.

Using Multiplexed CRISPR/Cas9 for Suppression of Cotton Leaf Curl Virus

In recent decades, Pakistan has suffered a decline in cotton production due to several factors, including insect pests, cotton leaf curl disease (CLCuD), and multiple abiotic stresses. CLCuD is a highly damaging plant disease that seriously limits cotton production in Pakistan. Recently, genome editing through CRISPR/Cas9 has revolutionized plant biology, especially to develop immunity in plants against viral diseases. Here we demonstrate multiplex CRISPR/Cas-mediated genome editing against CLCuD using transient transformation in N. benthamiana plants and cotton seedlings. The genomic sequences of cotton leaf curl viruses (CLCuVs) were obtained from NCBI and the guide RNA (gRNA) were designed to target three regions in the viral genome using CRISPR MultiTargeter. The gRNAs were cloned in pHSE401/pKSE401 containing Cas9 and confirmed through colony PCR, restriction analysis, and sequencing. Confirmed constructs were moved into Agrobacterium and subsequently used for transformation. Agroinfilteration in N. benthamiana revealed delayed symptoms (3–5 days) with improved resistance against CLCuD. In addition, viral titer was also low (20–40%) in infected plants co-infiltrated with Cas9-gRNA, compared to control plants (infected with virus only). Similar results were obtained in cotton seedlings. The results of transient expression in N. benthamiana and cotton seedlings demonstrate the potential of multiplex CRISPR/Cas to develop resistance against CLCuD. Five transgenic plants developed from three experiments showed resistance (60−70%) to CLCuV, out of which two were selected best during evaluation and screening. The technology will help breeding CLCuD-resistant cotton varieties for sustainable cotton production.

Complex evolution of novel red floral color in Petunia

Red flower color has arisen multiple times and is generally associated with hummingbird pollination. The majority of evolutionary transitions to red color proceeded from purple lineages and tend to be genetically simple, almost always involving a few loss-of-function mutations of major phenotypic effect. Here we report on the complex evolution of a novel red floral color in the hummingbird-pollinated Petunia exserta (Solanaceae) from a colorless ancestor. The presence of a red color is remarkable because the genus cannot synthesize red anthocyanins and P. exserta retains a nonfunctional copy of the key MYB transcription factor AN2. We show that moderate upregulation and a shift in tissue specificity of an AN2 paralog, DEEP PURPLE, restores anthocyanin biosynthesis in P. exserta. An essential shift in anthocyanin hydroxylation occurred through rebalancing the expression of three hydroxylating genes. Furthermore, the downregulation of an acyltransferase promotes reddish hues in typically purple pigments by preventing acyl group decoration of anthocyanins. This study presents a rare case of a genetically complex evolutionary transition toward the gain of a novel red color.

Cloning and Expression Analysis of Flavonoid 3', 5'-Hydroxylase Gene from Brunfelsia acuminata

The full-length cDNA sequence of F3′5′H gene from the Brunfelsia acuminata was obtained by RT-PCR and RACE, whose GenBank accession number is JQ678765. The sequence contains a 1521 bp open reading frame, 120 bp 5′UTR and 61 bp 3′UTR, encoding a total of 506 amino acids. The molecular mass of the predicted protein is 56.47 kDa with an estimated pI of 8.78, respectively. Sequence alignment showed that the amino acid sequence of F3′5′H was 91%, 87% and 84% with that of Petunia × hybrida, Nierembergia sp., Solanum tuberosum, respectively. Real-time quantitative PCR analysis showed that the expression of F3′5′H gene was different in petals of different days, which was the highest expression level on day 0 and significantly higher than other days. The results indicated that F3′5′H might play key role in flower color regulation and provide a theoretical reference for blue flower molecular breeding.

Genomic signatures of the evolution of defence against its natural enemies in the poisonous and medicinal plant Datura stramonium (Solanaceae)

Tropane alkaloids and terpenoids are widely used in the medicine and pharmaceutic industry and evolved as chemical defenses against herbivores and pathogens in the annual herb Datura stramonium (Solanaceae). Here, we present the first draft genomes of two plants from contrasting environments of D. stramonium. Using these de novo assemblies, along with other previously published genomes from 11 Solanaceae species, we carried out comparative genomic analyses to provide insights on the genome evolution of D. stramonium within the Solanaceae family, and to elucidate adaptive genomic signatures to biotic and abiotic stresses in this plant. We also studied, in detail, the evolution of four genes of D. stramonium-Putrescine N-methyltransferase, Tropinone reductase I, Tropinone reductase II and Hyoscyamine-6S-dioxygenase-involved in the tropane alkaloid biosynthesis. Our analyses revealed that the genomes of D. stramonium show signatures of expansion, physicochemical divergence and/or positive selection on proteins related to the production of tropane alkaloids, terpenoids, and glycoalkaloids as well as on R defensive genes and other important proteins related with biotic and abiotic pressures such as defense against natural enemies and drought.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1 Mb of the S₂ -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S_6a -haplotype of P. inflata and the S_N -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

Differential Gene Expression with an Emphasis on Floral Organ Size Differences in Natural and Synthetic Polyploids of Nicotiana tabacum (Solanaceae)

Floral organ size, especially the size of the corolla, plays an important role in plant reproduction by facilitating pollination efficiency. Previous studies have outlined a hypothesized organ size pathway. However, the expression and function of many of the genes in the pathway have only been investigated in model diploid species; therefore, it is unknown how these genes interact in polyploid species. Although correlations between ploidy and cell size have been shown in many systems, it is unclear whether there is a difference in cell size between naturally occurring and synthetic polyploids. To address these questions comparing floral organ size and cell size across ploidy, we use natural and synthetic polyploids of Nicotiana tabacum (Solanaceae) as well as their known diploid progenitors. We employ a comparative transcriptomics approach to perform analyses of differential gene expression, focusing on candidate genes that may be involved in floral organ size, both across developmental stages and across accessions. We see differential expression of several known floral organ candidate genes including ARF2, BIG BROTHER, and GASA/GAST1. Results from linear models show that ploidy, cell width, and cell number positively influence corolla tube circumference; however, the effect of cell width varies by ploidy, and diploids have a significantly steeper slope than both natural and synthetic polyploids. These results demonstrate that polyploids have wider cells and that polyploidy significantly increases corolla tube circumference.

Genome-wide identification of small G protein ROPs and their potential roles in Solanaceous family

Small GTPases function as molecular switches to active or inactive signaling cascades via binding or hydrolyzing GTP. A type of plant specific small GTPases, the ROPs are known to be involved in plant growth, development and immunity. We determined whether ROPs are conserved in Solanaceous species and whether they are involved in plant growth, development and resistance against Phytophthora capsisi. In genome-wide screening, a total of 66 ROPs in six Solanaceous species (SolROPs) were identified, including 16 ROPs in Solanum tuberosum L. (potato), 9 in Solanum lycopersicum L. (tomato), 5 in Solanum melongena L. (eggplant), 9 in Capsicum annuum L. (pepper), 13 in Nicotiana benthamiana Domin and 14 in Nicotiana tabacum L. (tobacco). Phylogenetic analysis revealed that 11 AtROPs and 66 SolROPs fall into five distinct clades (I-V) and hence a novel and systematic gene nomenclature was proposed. In addition, a comprehensive expression analysis was performed by making use of an online database. This revealed that ROP genes are differentially expressed during plant growth and development. Moreover, gene expression of SlROP-II.1 in S. lycopersicum could be significantly induced by P. capsici. Subsequently, SlROP-II.1 and its homologues in N. benthamiana and C. annuum (NbROP-II.1 and CaROP-II.1) were selected for functional analysis using virus-induced gene silencing. Infection assays with P. capsici on silenced plants revealed that SlROP-II.1, NbROP-II.1 and CaROP-II.1 play a role in P. capsici resistance, suggesting conserved function of ROP-II clade across different Solanaceous species. In addition, NbROP-II.1 is also involved in regulating plant growth and development. This study signified the diversity of Solanaceous ROPs and their potential roles in plant growth, development and immunity.

Ptr1 evolved convergently with RPS2 and Mr5 to mediate recognition of AvrRpt2 in diverse solanaceous species

The Ptr1 (Pseudomonas tomato race 1) locus in Solanum lycopersicoides confers resistance to strains of Pseudomonas syringae pv. tomato expressing AvrRpt2 and Ralstonia pseudosolanacearum expressing RipBN. Here we describe the identification and phylogenetic analysis of the Ptr1 gene. A single recombinant among 585 F2 plants segregating for the Ptr1 locus was discovered that narrowed the Ptr1 candidates to eight nucleotide-binding leucine-rich repeat protein (NLR)-encoding genes. From analysis of the gene models in the S. lycopersicoides genome sequence and RNA-Seq data, two of the eight genes emerged as the strongest candidates for Ptr1. One of these two candidates was found to encode Ptr1 based on its ability to mediate recognition of AvrRpt2 and RipBN when it was transiently expressed with these effectors in leaves of Nicotiana glutinosa. The ortholog of Ptr1 in tomato and in Solanum pennellii is a pseudogene. However, a functional Ptr1 ortholog exists in Nicotiana benthamiana and potato, and both mediate recognition of AvrRpt2 and RipBN. In apple and Arabidopsis, recognition of AvrRpt2 is mediated by the Mr5 and RPS2 proteins, respectively. Phylogenetic analysis places Ptr1 in a distinct clade compared with Mr5 and RPS2, and it therefore appears to have arisen by convergent evolution for recognition of AvrRpt2.

Evolution of a plant gene cluster in Solanaceae and emergence of metabolic diversity

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Genome-wide identification and characterisation of Aquaporins in Nicotiana tabacum and their relationships with other Solanaceae species

BACKGROUND

Cellular membranes are dynamic structures, continuously adjusting their composition, allowing plants to respond to developmental signals, stresses, and changing environments. To facilitate transmembrane transport of substrates, plant membranes are embedded with both active and passive transporters. Aquaporins (AQPs) constitute a major family of membrane spanning channel proteins that selectively facilitate the passive bidirectional passage of substrates across biological membranes at an astonishing 108 molecules per second. AQPs are the most diversified in the plant kingdom, comprising of five major subfamilies that differ in temporal and spatial gene expression, subcellular protein localisation, substrate specificity, and post-translational regulatory mechanisms; collectively providing a dynamic transportation network spanning the entire plant. Plant AQPs can transport a range of solutes essential for numerous plant processes including, water relations, growth and development, stress responses, root nutrient uptake, and photosynthesis. The ability to manipulate AQPs towards improving plant productivity, is reliant on expanding our insight into the diversity and functional roles of AQPs.

RESULTS

We characterised the AQP family from Nicotiana tabacum (NtAQPs; tobacco), a popular model system capable of scaling from the laboratory to the field. Tobacco is closely related to major economic crops (e.g. tomato, potato, eggplant and peppers) and itself has new commercial applications. Tobacco harbours 76 AQPs making it the second largest characterised AQP family. These fall into five distinct subfamilies, for which we characterised phylogenetic relationships, gene structures, protein sequences, selectivity filter compositions, sub-cellular localisation, and tissue-specific expression. We also identified the AQPs from tobacco's parental genomes (N. sylvestris and N. tomentosiformis), allowing us to characterise the evolutionary history of the NtAQP family. Assigning orthology to tomato and potato AQPs allowed for cross-species comparisons of conservation in protein structures, gene expression, and potential physiological roles.

CONCLUSIONS

This study provides a comprehensive characterisation of the tobacco AQP family, and strengthens the current knowledge of AQP biology. The refined gene/protein models, tissue-specific expression analysis, and cross-species comparisons, provide valuable insight into the evolutionary history and likely physiological roles of NtAQPs and their Solanaceae orthologs. Collectively, these results will support future functional studies and help transfer basic research to applied agriculture.

Tomato yellow leaf curl virus (TYLCV)-resistant tomatoes share molecular mechanisms sustaining resistance with their wild progenitor Solanum habrochaites but not with TYLCV-susceptible tomatoes

The wild tomato species Solanum habrochaites (Sh) has been used as a source for tomato yellow leaf curl virus (TYLCV) resistance in a breeding program to generate a TYLCV-resistant tomato line. Susceptible (S) and resistant (R) lines have been developed through this program. We compared the behavior of R, S and Sh tomato plants upon infection to find out whether the resistant phenotype of R plants originated from Sh. Results showed that mechanisms involving sugar-signaling (i.e., LIN6/HT1), water channels (i.e., TIP1;1), hormone homeostasis (i.e., ABA and SA) and urea accumulation were shared by S. habrochaites and R plants, but not by S. habrochaites and S tomatoes. This finding supports the hypothesis that these mechanisms were introgressed in the R genotype from the wild tomato progenitor during breeding for TYLCV resistance. Hence, identification of genes contributing to resistance to biotic stress from wild tomato species and their introgression into domestic plants ensures tomato supply and food security.

Optimised production of an anti-fungal antibody in Solanaceae hairy roots to develop new formulations against Candida albicans

BACKGROUND

Infections caused by fungi are often refractory to conventional therapies and urgently require the development of novel options, such as immunotherapy. To produce therapeutic antibodies, a plant-based expression platform is an attractive biotechnological strategy compared to mammalian cell cultures. In addition to whole plants, hairy roots (HR) cultures can be used, representing an expression system easy to build up, with indefinite growth while handled under containment conditions.

RESULTS

In this study the production in HR of a recombinant antibody, proved to be a good candidate for human immunotherapy against fungal infections, is reported. Expression and secretion of this antibody, in an engineered single chain (scFvFc) format, by HR from Nicotiana benthamiana and Solanum lycopersicum have been evaluated with the aim of directly using the deriving extract or culture medium against pathogenic fungi. Although both Solanaceae HR showed good expression levels (up to 68 mg/kg), an optimization of rhizosecretion was only obtained for N. benthamiana HR. A preliminary assessment to explain this result highlighted the fact that not only the presence of proteases, but also the chemical characteristics of the growth medium, can influence antibody yield, with implications on recombinant protein production in HR. Finally, the antifungal activity of scFvFc 2G8 antibody produced in N. benthamiana HR was evaluated in Candida albicans growth inhibition assays, evidencing encouraging results.

CONCLUSIONS

Production of this anti-fungal antibody in HR of N. benthamiana and S. lycopersicum elucidated factors affecting pharming in this system and allowed to obtain promising ready-to-use immunotherapeutics against C. albicans.

Ancestral segmental duplication in Solanaceae is responsible for the origin of CRCa-CRCb paralogues in the family

CRABS CLAW (CRC), a member of YABBY transcription factor family, has been previously reported to be principally involved in carpel development across angiosperms, and nectary development in core eudicots. Most of the studies suggest that CRC exists as a single copy gene, except in the Solanaceae where CRC occurs as paralogous pairs-CRCa-CRCb in Solanum lycopersicum, and CRC1-CRC2 in Petunia hybrida. In spite of their crucial role in carpel and nectary development, there is no information about the evolutionary history of the CRC paralogy in Solanaceae and whether the paralogy extends beyond Solanaceae. We analyzed homologues of CRC across angiosperms including genome sequence of fourteen species of Solanaceae available at Sol Genomics Network database, Phytozome and NCBI, to address the questions. Our phylogenetic reconstruction across angiosperms combined with comparative genomic, microsynteny and genome-fractionation analyses across the Solanaceae genomes revealed that (1) the CRCa-CRCb lineage is represented by a single copy in other flowering plants; (2) putative homologues of CRCa and CRCb are present in all the Solanaceae genomes studied; (3) the CRCa-CRCb paralogy in Solanaceae is associated with a large segmental duplication within Solanaceae (perhaps in its common ancestor), and (4) the duplicated segments have undergone different degrees of retention and loss of genes. Also, the CRC gene lineage expanded in Solanaceae following Solanaceae-α hexaploidy event and that two CRC duplicate copies were subsequently retained during the course of evolution. Besides the first detailed description of CRC evolution in Solanaceae, the study identifies potential candidate genes for future functional investigations.

Transcriptome sequencing reveals the effect of biochar improvement on the development of tobacco plants before and after topping

The application of biochar is one of the most useful methods for improving soil quality, which is of the utmost significance for the continuous production of crops. As there are no conclusive studies on the specific effects of biochar application on tobacco quality, this study aimed to improve the yield and quality of tobacco as a model crop for economic and genetic research in southern China, by such application. We used transcriptome sequencing to reveal the effects of applied biochar on tobacco development before and after topping. Our results showed that topping affected carbon and nitrogen metabolism, photosynthesis and secondary metabolism in the tobacco plants, while straw biochar-application to the soil resulted in amino acid and lipid synthesis; additionally, it affected secondary metabolism of the tobacco plants through carbon restoration and hormonal action, before and after topping. In addition to the new insights into the impact of biochar on crops, our findings provide a basis for biochar application measures in tobacco and other crops.

Expression and function of the bHLH genes ALCATRAZ and SPATULA in selected Solanaceae species

The genetic mechanisms underlying fruit development have been identified in Arabidopsis and have been comparatively studied in tomato as a representative of fleshy fruits. However, comparative expression and functional analyses on the bHLH genes downstream the genetic network, ALCATRAZ (ALC) and SPATULA (SPT), which are involved in the formation of the dehiscence zone in Arabidopsis, have not been functionally studied in the Solanaceae. Here, we perform detailed expression and functional studies of ALC/SPT homologs in Nicotiana obtusifolia with capsules, and in Capsicum annuum and Solanum lycopersicum with berries. In Solanaceae, ALC and SPT genes are expressed in leaves, and all floral organs, especially in petal margins, stamens and carpels; however, their expression changes during fruit maturation according to the fruit type. Functional analyses show that downregulation of ALC/SPT genes does not have an effect on gynoecium patterning; however, they have acquired opposite roles in petal expansion and have been co-opted in leaf pigmentation in Solanaceae. In addition, ALC/SPT genes repress lignification in time and space during fruit development in Solanaceae. Altogether, some roles of ALC and SPT genes are different between Brassicaceae and Solanaceae; while the paralogs have undergone some subfunctionalization in the former they are mostly redundant in the latter.

Assessing Anthocyanin Biosynthesis in Solanaceae as a Model Pathway for Secondary Metabolism

Solanaceae have played an important role in elucidating how flower color is specified by the flavonoid biosynthesis pathway (FBP), which produces anthocyanins and other secondary metabolites. With well-established reverse genetics tools and rich genomic resources, Solanaceae provide a robust framework to examine the diversification of this well-studied pathway over short evolutionary timescales and to evaluate the predictability of genetic perturbation on pathway flux. Genomes of eight Solanaceae species, nine related asterids, and four rosids were mined to evaluate variation in copy number of the suite of FBP enzymes involved in anthocyanin biosynthesis. Comparison of annotation sources indicated that the NCBI annotation pipeline generated more and longer FBP annotations on average than genome-specific annotation pipelines. The pattern of diversification of each enzyme among asterids was assessed by phylogenetic analysis, showing that the CHS superfamily encompasses a large paralogous family of ancient and recent duplicates, whereas other FBP enzymes have diversified via recent duplications in particular lineages. Heterologous expression of a pansy F3′5′H gene in tobacco changed flower color from pink to dark purple, demonstrating that anthocyanin production can be predictably modified using reverse genetics. These results suggest that the Solanaceae FBP could be an ideal system to model genotype-to-phenotype interactions for secondary metabolism.

Early consequences of allopolyploidy alter floral evolution in Nicotiana (Solanaceae)

BACKGROUND

Polyploidy has played a major role in angiosperm evolution. Previous studies have examined polyploid phenotypes in comparison to their extant progenitors, but not in context of predicted progenitor phenotypes at allopolyploid origin. In addition, differences in the trends of polyploid versus diploid evolution have not been investigated. We use ancestral character-state reconstructions to estimate progenitor phenotype at allopolyploid origin to determine patterns of polyploid evolution leading to morphology of the extant species. We also compare trends in diploid versus allopolyploid evolution to determine if polyploidy modifies floral evolutionary patterns.

RESULTS

Predicting the ancestral phenotype of a nascent allopolyploid from reconstructions of diploid phenotypes at the time of polyploid formation generates different phenotype predictions than when extant diploid phenotypes are used, the outcome of which can alter conclusions about polyploid evolution; however, most analyses yield the same results. Using ancestral reconstructions of diploid floral phenotypes indicate that young polyploids evolve shorter, wider corolla tubes, but older polyploids and diploids do not show any detectable evolutionary trends. Lability of the traits examined (floral shape, corolla tube length, and corolla tube width) differs across young and older polyploids and diploids. Corolla length is more evolutionarily labile in older polyploids and diploids. Polyploids do not display unique suites of floral characters based on both morphological and color traits, but some suites of characters may be evolving together and seem to have arisen multiple times within Nicotiana, perhaps due to the influence of pollinators.

CONCLUSIONS

Young polyploids display different trends in floral evolution (shorter, wider corolla tubes, which may result in more generalist pollination) than older polyploids and diploids, suggesting that patterns of divergence are impacted by the early consequences of allopolyploidy, perhaps arising from genomic shock and/or subsequent genome stabilization associated with diploidization. Convergent evolution in floral morphology and color in Nicotiana can be consistent with pollinator preferences, suggesting that pollinators may have shaped floral evolution in Nicotiana.

A critical comparison of technologies for a plant genome sequencing project

BACKGROUND

A high-quality genome sequence of any model organism is an essential starting point for genetic and other studies. Older clone-based methods are slow and expensive, whereas faster, cheaper short-read-only assemblies can be incomplete and highly fragmented, which minimizes their usefulness. The last few years have seen the introduction of many new technologies for genome assembly. These new technologies and associated new algorithms are typically benchmarked on microbial genomes or, if they scale appropriately, on larger (e.g., human) genomes. However, plant genomes can be much more repetitive and larger than the human genome, and plant biochemistry often makes obtaining high-quality DNA that is free from contaminants difficult. Reflecting their challenging nature, we observe that plant genome assembly statistics are typically poorer than for vertebrates.

RESULTS

Here, we compare Illumina short read, Pacific Biosciences long read, 10x Genomics linked reads, Dovetail Hi-C, and BioNano Genomics optical maps, singly and combined, in producing high-quality long-range genome assemblies of the potato species Solanum verrucosum. We benchmark the assemblies for completeness and accuracy, as well as DNA compute requirements and sequencing costs.

CONCLUSIONS

The field of genome sequencing and assembly is reaching maturity, and the differences we observe between assemblies are surprisingly small. We expect that our results will be helpful to other genome projects, and that these datasets will be used in benchmarking by assembly algorithm developers.

Comparative analysis of repetitive sequences among species from the potato and the tomato clades

BACKGROUND AND AIMS

The genus Solanum includes important vegetable crops and their wild relatives. Introgression of their useful traits into elite cultivars requires effective recombination between hom(e)ologues, which is partially determined by genome sequence differentiation. In this study we compared the repetitive genome fractions of wild and cultivated species of the potato and tomato clades in a phylogenetic context.

METHODS

Genome skimming followed by a clustering approach was used as implemented in the RepeatExplorer pipeline. Repeat classes were annotated and the sequences of their main domains were compared.

KEY RESULTS

Repeat abundance and genome size were correlated and the larger genomes of species in the tomato clade were found to contain a higher proportion of unclassified elements. Families and lineages of repetitive elements were largely conserved between the clades, but their relative proportions differed. The most abundant repeats were Ty3/Gypsy elements. Striking differences in abundance were found in the highly dynamic Ty3/Gypsy Chromoviruses and Ty1/Copia Tork elements. Within the potato clade, early branching Solanum cardiophyllum showed a divergent repeat profile. There were also contrasts between cultivated and wild potatoes, mostly due to satellite amplification in the cultivated species. Interspersed repeat profiles were very similar among potatoes. The repeat profile of Solanum etuberosum was more similar to that of the potato clade.

CONCLUSIONS

The repeat profiles in Solanum seem to be very similar despite genome differentiation at the level of collinearity. Removal of transposable elements by unequal recombination may have been responsible for structural rearrangements across the tomato clade. Sequence variability in the tomato clade is congruent with clade-specific amplification of repeats after its divergence from S. etuberosum and potatoes. The low differentiation among potato and its wild relatives at the level of interspersed repeats may explain the difficulty in discriminating their genomes by genomic in situ hybridization techniques.

24-nt reproductive phasiRNAs are broadly present in angiosperms

Small RNAs are key regulators in plant growth and development. One subclass, phased siRNAs (phasiRNAs) require a trigger microRNA for their biogenesis. In grasses, two pathways yield abundant phasiRNAs during anther development; miR2275 triggers one class, 24-nt phasiRNAs, coincident with meiosis, while a second class of 21-nt phasiRNAs are present in premeiotic anthers. Here we report that the 24-nt phasiRNA pathway is widely present in flowering plants, indicating that 24-nt reproductive phasiRNAs likely originated with the evolutionary emergence of anthers. Deep comparative genomic analyses demonstrated that this miR2275/24-nt phasiRNA pathway is widely present in eudicots plants, however, it is absent in legumes and in the model plant Arabidopsis, demonstrating a dynamic evolutionary history of this pathway. In Solanaceae species, 24-nt phasiRNAs were observed, but the miR2275 trigger is missing and some loci displaying 12-nt phasing. Both the miR2275-triggered and Solanaceae 24-nt phasiRNAs are enriched in meiotic stages, implicating these phasiRNAs in anther and/or pollen development, a spatiotemporal pattern consistent in all angiosperm lineages that deploy them.

Repeated evolution of a morphological novelty: a phylogenetic analysis of the inflated fruiting calyx in the Physalideae tribe (Solanaceae)

PREMISE OF THE STUDY

The evolution of novel fruit morphologies has been integral to the success of angiosperms. The inflated fruiting calyx, in which the balloon-like calyx swells to completely surround the fruit, has evolved repeatedly across angiosperms and is postulated to aid in protection and dispersal. We investigated the evolution of this trait in the tomatillos and their allies (Physalideae, Solanaceae).

METHODS

The Physalideae phylogeny was estimated using four regions (ITS, LEAFY, trnL-F, waxy) with maximum likelihood (ML) and Bayesian inference. Under the best-fitting ML model of trait evolution, we estimated ancestral states along with the numbers of gains and losses of fruiting calyx accrescence and inflation with Bayesian stochastic mapping. Also, phylogenetic signal in calyx morphology was examined with two metrics (parsimony score and Fritz and Purvis's D).

KEY RESULTS

Based on our well-resolved and densely sampled phylogeny, we infer that calyx evolution has proceeded in a stepwise and directional fashion, from non-accrescent to accrescent to inflated. In total, we inferred 24 gains of accrescence, 24 subsequent transitions to a fully inflated calyx, and only two reversals. Despite this lability, fruiting calyx accrescence and inflation showed strong phylogenetic signal.

CONCLUSIONS

Our phylogeny greatly improves the resolution of Physalideae and highlights the need for taxonomic work. The comparative analyses reveal that the inflated fruiting calyx has evolved many times and that the trajectory toward this phenotype is generally stepwise and irreversible. These results provide a strong foundation for studying the genetic and developmental mechanisms responsible for the repeated origins of this charismatic fruit trait.

Tandem gene duplication and recombination at the AT3 locus in the Solanaceae, a gene essential for capsaicinoid biosynthesis in Capsicum

Capsaicinoids are compounds synthesized exclusively in the genus Capsicum and are responsible for the burning sensation experienced when consuming hot pepper fruits. To date, only one gene, AT3, a member of the BAHD family of acyltransferases, is currently known to have a measurable quantitative effect on capsaicinoid biosynthesis. Multiple AT3 paralogs exist in the Capsicum genome, but their evolutionary relationships have not been characterized well. Recessive alleles at this locus result in absence of capsaicinoids in pepper fruit. To explore the evolution of AT3 in Capsicum and the Solanaceae, we sequenced this gene from diverse Capsicum genotypes and species, along with a number of representative solanaceous taxa. Our results revealed that the coding region of AT3 is highly conserved throughout the family. Further, we uncovered a tandem duplication that predates the diversification of the Solanaceae taxa sampled in this study. This pair of tandem duplications were designated AT3-1 and AT3-2. Sequence alignments showed that the AT3-2 locus, a pseudogene, retains regions of amino acid conservation relative to AT3-1. Gene tree estimation demonstrated that AT3-1 and AT3-2 form well supported, distinct clades. In C. rhomboideum, a non-pungent basal Capsicum species, we describe a recombination event between AT3-1 and AT3-2 that modified the putative active site of AT3-1, also resulting in a frame-shift mutation in the second exon. Our data suggest that duplication of the original AT3 representative, in combination with divergence and pseudogene degeneration, may account for the patterns of sequence divergence and punctuated amino acid conservation observed in this study. Further, an early rearrangement in C. rhomboidium could account for the absence of pungency in this Capsicum species.

Four new species of Capsicum (Solanaceae) from the tropical Andes and an update on the phylogeny of the genus

Four new species of Capsicum (Capsiceae, Solanaceae) from Andean tropical forests in South America are described. Capsicum benoistii Hunz. ex Barboza sp. nov. (incertae sedis) is endemic to a restricted area in south-central Ecuador and is most similar to the more widespread C. geminifolium (Dammer) Hunz. (Colombia, Ecuador, and Peru). Capsicum piuranum Barboza & S. Leiva sp. nov. (Andean clade) is found in northern Peru (Department Piura) and is morphologically most similar to C. caballeroi M. Nee of the Bolivian yungas (Departments Santa Cruz and Cochabamba) but closely related to C. geminifolium and C. lycianthoides Bitter. Capsicum longifolium Barboza & S. Leiva sp. nov. (Andean clade) occurs from northern Peru (Departments Amazonas, Cajamarca, and Piura) to southern Ecuador (Province Zamora-Chinchipe), and is morphologically most similar to C. dimorphum (Miers) Kuntze (Colombia, Ecuador, and Peru). Capsicum neei Barboza & X. Reyes sp. nov. (Bolivian clade) is endemic to southeastern Bolivia (Departments Chuquisaca and Santa Cruz) in the Boliviano-Tucumano Forest, is morphologically most similar to another Bolivian endemic species C. minutiflorum Rusby (Hunz.), and is closely related to C. caballeroi. Complete descriptions, illustrations, distributions and conservation assessments of all new species are given. Chromosome numbers for C. piuranum and C. longifolium are also provided. Three of the new species were included in a new phylogenetic analysis for Capsicum; their positions were strongly resolved within clades previously recognized in the genus.

Evolution by duplication: paleopolyploidy events in plants reconstructed by deciphering the evolutionary history of VOZ transcription factors

BACKGROUND

Facilitated by the rapid progress of sequencing technology, comparative genomic studies in plants have unveiled recurrent whole genome duplication (i.e. polyploidization) events throughout plant evolution. The evolutionary past of plant genes should be analyzed in a background of recurrent polyploidy events in distinctive plant lineages. The Vascular Plant One Zinc-finger (VOZ) gene family encode transcription factors associated with a number of important traits including control of flowering time and photoperiodic pathways, but the evolutionary trajectory of this gene family remains uncharacterized.

RESULTS

In this study, we deciphered the evolutionary history of the VOZ gene family by analyses of 107 VOZ genes in 46 plant genomes using integrated methods: phylogenic reconstruction, Ks-based age estimation and genomic synteny comparisons. By scrutinizing the VOZ gene family phylogeny the core eudicot γ event was well circumscribed, and relics of the precommelinid τ duplication event were detected by incorporating genes from oil palm and banana. The more recent T and ρ polyploidy events, closely coincident with the species diversification in Solanaceae and Poaceae, respectively, were also identified. Other important polyploidy events captured included the "salicoid" event in poplar and willow, the "early legume" and "soybean specific" events in soybean, as well as the recent polyploidy event in Physcomitrella patens. Although a small transcription factor gene family, the evolutionary history of VOZ genes provided an outstanding record of polyploidy events in plants. The evolutionary past of VOZ gene family demonstrated a close correlation with critical plant polyploidy events which generated species diversification and provided answer to Darwin's "abominable mystery".

CONCLUSIONS

We deciphered the evolutionary history of VOZ transcription factor family in plants and ancestral polyploidy events in plants were recapitulated simultaneously. This analysis allowed for the generation of an idealized plant gene tree demonstrating distinctive retention and fractionation patterns following polyploidy events.

Application and development of genome editing technologies to the Solanaceae plants

Genome editing technology using artificial nucleases, including zinc finger nuclease (ZFN), transcription activator-like effector nuclease (TALEN), and clustered regulatory interspaced short palindromic repeats (CRISPR)-Cas9, can mutagenize the target sites of genes of interest. This technology has been successfully applied in several crops, including the Solanaceae plants, such as tomato, potato, tobacco, and petunia. Among the three nucleases, CRISPR-Cas9 is the best for breeding, crop improvement, and the functional analysis of genes of interest, because of its simplicity and high efficiency. Although the technology is useful for reverse genetics, its use in plants is limited due to a lack of regeneration protocols and sequence information. In this review, the present status of genome editing technology in Solanaceae plants is described, and techniques that may improve genome editing technologies are discussed.

The genus Datura L. (Solanaceae) in Mexico and Spain - Ethnobotanical perspective at the interface of medical and illicit uses

ETHNOPHARMACOLOGICAL RELEVANCE

The different species of the genus Datura have been used traditionally by some pre-Columbian civilizations, as well as in medieval rituals linked to magic and witchcraft in both Mexico and Europe. It is also noteworthy the use of different alkaloids obtained from the plants for medicinal purposes in the treatment of various groups of diseases, especially of the respiratory and muscularskeletal systems.

AIM OF THE STUDY

A review of the ethnobotanical uses of the genus Datura in Mexico and Spain has been conducted. We focus on the medicinal and ritualistic uses included in modern ethnobotanical studies, emphasizing the historical knowledge from post-colonial American Codices and medieval European texts. Datura's current social emergency as a drug of recreation and leisure, as well as its link to crimes of sexual abuse is also considered. The work is completed with some notes about the distribution and ecology of the different species and a phytochemical and pharmacological review of Datura alkaloids, necessary to understand their arrival in Europe and the ethnobotanical uses made since then MATERIALS AND METHODS: A literature review and compilation of information on traditional medicinal uses of the genus has been carried out from the main electronic databases. Traditional volumes (codices) have also been consulted in libraries of different institutions. Consultations have been made with the National Toxicological Services of Spain and Mexico for toxicological data.

RESULTS

A total of 118 traditional uses were collected in both territories, 111 medicinal ones to be applied in 76 conditions or symptoms included in 13 pathological groups. Although there are particular medicinal uses in the two countries, we found up to 15 similar uses, of which 80% were previously mentioned in post-Colonial American codices. Applications in the treatment of asthma and rheumatism are also highlighted. Apart from medicinal uses, it is worth noting their cultural and social uses, in the case of Mexico relating to diseases such as being scared, astonishment or falling in love, and in the case of Spain, as a recreational drug and lately, for criminal purposes.

CONCLUSIONS

This review highlights the variety of uses traditionally given to the different species in both territories. The fact that most of the coincident or similar uses in both countries also appear in the classical codices can be found an example of the flow, not only of the plants from America to Europe, but also of their associated information. It is also relevant that particular uses have derived in both countries, reflecting the difference in the cultural factors and traditions linked to rituals and cultural practices. Finally, the significant growth of Datura consumption in recent years as a drug of leisure and recreation, as well as in crimes of sexual submission, should be considered as research of maximum relevance in the field of forensic botany and toxicology.

Evolution of genes associated with gynoecium patterning and fruit development in Solanaceae

BACKGROUND AND AIMS

The genetic basis of fruit development has been extensively studied in Arabidopsis, where major transcription factors controlling valve identity (i.e. FRUITFULL), replum development (i.e. REPLUMLESS) and the differentiation of the dehiscence zones (i.e. SHATTERPROOF, INDEHISCENT and ALCATRAZ) have been identified. This gene regulatory network in other flowering plants is influenced by duplication events during angiosperm diversification. Here we aim to characterize candidate fruit development genes in the Solanaceae and compare them with those of Brassicaceae.

METHODS

ALC/SPT, HEC/IND, RPL and AG/SHP homologues were isolated from publicly available databases and from our own transcriptomes of Brunfelsia australis and Streptosolen jamesonii. Maximum likelihood phylogenetic analyses were performed for each of the gene lineages. Shifts in protein motifs, as well as expression patterns of all identified homologues, are shown in dissected floral organs and fruits in different developmental stages of four Solanaceae species exhibiting different fruit types.

KEY RESULTS

Each gene lineage has undergone different duplication time-points, resulting in very different genetic complements in the Solanaceae when compared with the Brassicaceae. In general, Solanaceae species have more copies of HEC1/2 and RPL than Brassicaceae, have fewer copies of SHP and the same number of copies of AG, ALC and SPT. Solanaceae lack IND orthologues, but have pre-duplication HEC3 homologues. The expression analyses showed opposite expression of SPT and ALC orthologues between dry- and fleshy-fruited species during fruit maturation. Fleshy-fruited species turn off RPL and SPT orthologues during maturation.

CONCLUSIONS

The gynoecium patterning and fruit developmental genetic network in the Brassicaceae cannot be directly extrapolated to the Solanaceae. In Solanaceae ALC, SPT and RPL contribute differently to maturation of dry dehiscent and fleshy fruits, whereas HEC genes are not generally expressed in the gynoecium. RPL genes have broader expression patterns than expected.

Comparative transcriptomics and genomic patterns of discordance in Capsiceae (Solanaceae)

The integration of genomics and phylogenetics allows new insight into the structure of gene tree discordance, the relationships among gene position, gene history, and rate of evolution, as well as the correspondence of gene function, positive selection, and gene ontology enrichment across lineages. We explore these issues using the tribe Capsiceae (Solanaceae), which is comprised of the genera Lycianthes and Capsicum (peppers). In combining the annotated genomes of Capsicum with newly sequenced transcriptomes of four species of Lycianthes and Capsicum, we develop phylogenies for 6747 genes, and construct a backbone species tree using both concordance and explicit phylogenetic network approaches. We quantify phylogenetic discordance among individual gene trees, measure their rates of synonymous and nonsynonymous substitution, and test whether they were positively selected along any branch of the phylogeny. We then map these genes onto the annotated Capsicum genome and test whether rates of evolution, gene history, and gene ontology vary significantly with gene position. We observed substantial discordance among gene trees. A bifurcating species tree placing Capsicum within a paraphyletic Lycianthes was supported over all phylogenetic networks. Rates of synonymous and nonsynonymous substitution varied 41-fold and 130-fold among genes, respectively, and were significantly lower in pericentromeric regions. We found that results of concordance tree analyses vary depending on the subset of genes used, and that genes within the pericentromeric regions only capture a portion of the observed discordance. We identified 787 genes that have been positively selected throughout the diversification history of Capsiceae, and discuss the importance of these genes as targets for investigation of economically important traits in the domesticated peppers.

Plant reference genes for development and stress response studies

Many reference genes are used by different laboratories for gene expression analyses to indicate the relative amount of input RNA/DNA in the experiment. These reference genes are supposed to show least variation among the treatments and with the control sets in a given experiment. However, expression of reference genes varies significantly from one set of experiment to the other. Thus, selection of reference genes depends on the experimental conditions. Sometimes the average expression of two or three reference genes is taken as standard. This review consolidated the details of about 120 genes attempted for normalization during comparative expression analysis in 16 different plants. Plant species included in this review are Arabidopsis thaliana, cotton (Gossypium hirsutum), tobacco (Nicotiana benthamiana and N. tabacum), soybean (Glycine max), rice (Oryza sativa), blueberry (Vaccinium corymbosum), tomato (Solanum lycopersicum), wheat (Triticum aestivum), potato (Solanum tuberosum), sugar cane (Saccharum sp.), carrot (Daucus carota), coffee (Coffea arabica), cucumber (Cucumis sativus), kiwi (Actinidia deliciosa) and grape (Vitis vinifera). The list includes model and cultivated crop plants from both monocot and dicot classes. We have categorized plant-wise the reference genes that have been used for expression analyses in any or all of the four different conditions such as biotic stress, abiotic stress, developmental stages and various organs and tissues, reported till date. This review serves as a guide during the reference gene hunt for gene expression analysis studies.

A novel R3 MYB transcriptional repressor associated with the loss of floral pigmentation in Iochroma

Losses of floral pigmentation represent one of the most common evolutionary transitions in flower color, yet the genetic basis for these changes has been elucidated in only a handful of cases. Here we used crossing studies, bulk-segregant RNA sequencing, phylogenetic analyses and functional tests to identify the gene(s) responsible for the transition to white flowers in Iochroma loxense. Crosses between I. loxense and its blue-flowered sister species, I. cyaneum, suggested that a single locus controls the flower color difference and that the white allele causes a nearly complete loss of pigmentation. Examining sequence variation across phenotypic pools from the crosses, we found that alleles at a novel R3 MYB transcription factor were tightly associated with flower color variation. This gene, which we term MYBL1, falls into a class of MYB transcriptional repressors and, accordingly, higher expression of this gene is associated with downregulation of multiple anthocyanin pigment pathway genes. We confirmed the repressive function of MYBL1 through stable transformation of Nicotiana. The mechanism underlying the evolution of white flowers in I. loxense differs from that uncovered in previous studies, pointing to multiple mechanisms for achieving fixed transitions in flower color intensity.