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.