Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Image-Based Quantitative Analysis of Epidermal Morphology in Wild Potato Leaves

The epidermal leaf patterns of plants exhibit remarkable diversity in cell shapes, sizes, and arrangements, driven by environmental interactions that lead to significant adaptive changes even among closely related species. The Solanaceae family, known for its high diversity of adaptive epidermal structures, has traditionally been studied using qualitative phenotypic descriptions. To advance this, we developed a workflow combining multi-scale computer vision, image processing, and data analysis to extract digital descriptors for leaf epidermal cell morphology. Applied to nine wild potato species, this workflow quantified key morphological parameters, identifying descriptors for trichomes, stomata, and pavement cells, and revealing interdependencies among these traits. Principal component analysis (PCA) highlighted two main axes, accounting for 45% and 21% of variance, corresponding to features such as guard cell shape, trichome length, stomatal density, and trichome density. These axes aligned well with the historical and geographical origins of the species, separating southern from Central American species, and forming distinct clusters for monophyletic groups. This workflow thus establishes a quantitative foundation for investigating leaf epidermal cell morphology within phylogenetic and geographic contexts.

Metabolic and transcriptomic analyses of nectaries reveal differences in the mechanism of nectar production between monocots (Ananas comosus) and dicots (Nicotiana tabacum)

BACKGROUND

Nectar is offered by numerous flowering plants to attract pollinators. To date, the production and secretion of nectar have been analyzed mainly in eudicots, particularly rosids such as Arabidopsis. However, due to the enormous diversity of flowering plants, further research on other plant species, especially monocots, is needed. Ananas comosus (monocot) is an economically important species that is ideal for such analyses because it produces easily accessible nectar in sufficient quantities. In addition, the analyses were also carried out with Nicotiana tabacum (dicot, asterids) for comparison.

RESULTS

We performed transcriptome sequencing (RNA-Seq) analyses of the nectaries of Ananas comosus and Nicotiana tabacum, to test whether the mechanisms described for nectar production and secretion in Arabidopsis are also present in these plant species. The focus of these analyses is on carbohydrate metabolism and transport (e.g., sucrose-phosphate synthases, invertases, sucrose synthases, SWEETs and further sugar transporters). In addition, the metabolites were analyzed in the nectar, nectaries and leaves of both plant species to address the question of whether concentration gradients for different metabolites exist between the nectaries and nectar The nectar of N. tabacum contains large amounts of glucose, fructose and sucrose, and the sucrose concentration in the nectar appears to be similar to the sucrose concentration in the nectaries. Nectar production and secretion in this species closely resemble corresponding processes in some other dicots, including sucrose synthesis in nectaries and sucrose secretion by SWEET9. The nectar of A. comosus also contains large amounts of glucose, fructose and sucrose and in this species the sucrose concentration in the nectar appears to be higher than the sucrose concentration in the nectaries. Furthermore, orthologs of SWEET9 generally appear to be absent in A. comosus and other monocots. Therefore, sucrose export by SWEETs from nectaries into nectar can be excluded; rather, other mechanisms, such as active sugar export or exocytosis, are more likely.

CONCLUSION

The mechanisms of nectar production and secretion in N. tabacum appear to be largely similar to those in other dicots, whereas in the monocotyledonous species A. comosus, different synthesis and transport processes are involved.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

Genetic analysis of heat tolerance in hot pepper: insights from comprehensive phenotyping and QTL mapping

High temperatures present a formidable challenge to the cultivation of hot pepper, profoundly impacting not only vegetative growth but also leading to flower and fruit abscission, thereby causing a significant reduction in yield. To unravel the intricate genetic mechanisms governing heat tolerance in hot pepper, an F2 population was developed through the crossing of two distinct genotypes exhibiting contrasting heat tolerance characteristics: DLS-161-1 (heat tolerant) and DChBL-240 (heat susceptible). The F2 population, along with the parental lines, was subjected to comprehensive phenotyping encompassing diverse morphological, physiological, and biochemical heat-related traits under high temperature conditions (with maximum temperature ranging from 31 to 46.5°C and minimum temperature from 15.4 to 30.5°C). Leveraging the Illumina Nova Seq-6000 platform, Double digest restriction-site associated DNA sequencing (ddRAD-seq) was employed to generate 67.215 Gb data, with subsequent alignment of 218.93 million processed reads against the reference genome of Capsicum annuum. Subsequent variant calling and ordering resulted in 5806 polymorphic SNP markers grouped into 12 LGs. Further QTL analysis identified 64 QTLs with LOD values ranging from 2.517 to 11.170 and explained phenotypic variance ranging from 4.05 to 19.39%. Among them, 21 QTLs, explaining more than 10% phenotypic variance, were identified as major QTLs controlling 9 morphological, 3 physiological, and 2 biochemical traits. Interestingly, several QTLs governing distinct parameters were found to be colocalized, suggesting either a profound correlation between the QTLs regulating these traits or their significant genomic proximity. In addition to the QTLs, we also identified 368380 SSR loci within the identified QTL regions, dinucleotides being the most abundant type (211,381). These findings provide valuable insights into the genetics of heat tolerance in hot peppers. The identified QTLs and SSR markers offer opportunities to develop heat-tolerant varieties, ensuring better crop performance under high-temperature conditions.

The LEA gene family in tomato and its wild relatives: genome-wide identification, structural characterization, expression profiling, and role of SlLEA6 in drought stress

BACKGROUND

Late embryogenesis abundant (LEA) proteins are widely distributed in higher plants and play crucial roles in regulating plant growth and development processes and resisting abiotic stress. Cultivated tomato (Solanum lycopersicum) is an important vegetable crop worldwide; however, its growth, development, yield, and quality are currently severely constrained by abiotic stressors. In contrast, wild tomato species are more tolerant to abiotic stress and can grow normally in extreme environments. The main objective of this study was to identify, characterize, and perform gene expression analysis of LEA protein families from cultivated and wild tomato species to mine candidate genes and determine their potential role in abiotic stress tolerance in tomatoes.

RESULTS

Total 60, 69, 65, and 60 LEA genes were identified in S. lycopersicum, Solanum pimpinellifolium, Solanum pennellii, and Solanum lycopersicoides, respectively. Characterization results showed that these genes could be divided into eight clusters, with the LEA_2 cluster having the most members. Most LEA genes had few introns and were non-randomly distributed on chromosomes; the promoter regions contained numerous cis-acting regulatory elements related to abiotic stress tolerance and phytohormone responses. Evolutionary analysis showed that LEA genes were highly conserved and that the segmental duplication event played an important role in evolution of the LEA gene family. Transcription and expression pattern analyses revealed different regulatory patterns of LEA genes between cultivated and wild tomato species under normal conditions. Certain S. lycopersicum LEA (SlLEA) genes showed similar expression patterns and played specific roles under different abiotic stress and phytohormone treatments. Gene ontology and protein interaction analyses showed that most LEA genes acted in response to abiotic stimuli and water deficit. Five SlLEA proteins were found to interact with 11 S. lycopersicum WRKY proteins involved in development or resistance to stress. Virus-induced gene silencing of SlLEA6 affected the antioxidant and reactive oxygen species defense systems, increased the degree of cellular damage, and reduced drought resistance in S. lycopersicum.

CONCLUSION

These findings provide comprehensive information on LEA proteins in cultivated and wild tomato species and their possible functions under different abiotic and phytohormone stresses. The study systematically broadens our current understanding of LEA proteins and candidate genes and provides a theoretical basis for future functional studies aimed at improving stress resistance in tomato.

Keeping time in the dark: Potato diel and circadian rhythmic gene expression reveals tissue-specific circadian clocks

The circadian clock is an internal molecular oscillator and coordinates numerous physiological processes through regulation of molecular pathways. Tissue-specific clocks connected by mobile signals have previously been found to run at different speeds in Arabidopsis thaliana tissues. However, tissue variation in circadian clocks in crop species is unknown. In this study, leaf and tuber global gene expression in cultivated potato under cycling and constant environmental conditions was profiled. In addition, we used a circadian-regulated luciferase reporter construct to study tuber gene expression rhythms. Diel and circadian expression patterns were present among 17.9% and 5.6% of the expressed genes in the tuber. Over 500 genes displayed differential tissue specific diel phases. Intriguingly, few core circadian clock genes had circadian expression patterns, while all such genes were circadian rhythmic in cultivated tomato leaves. Furthermore, robust diel and circadian transcriptional rhythms were observed among detached tubers. Our results suggest alternative regulatory mechanisms and/or clock composition is present in potato, as well as the presence of tissue-specific independent circadian clocks. We have provided the first evidence of a functional circadian clock in below-ground storage organs, holding important implications for other storage root and tuberous crops.

Genome evolution and diversity of wild and cultivated potatoes

Potato (Solanum tuberosum L.) is the world's most important non-cereal food crop, and the vast majority of commercially grown cultivars are highly heterozygous tetraploids. Advances in diploid hybrid breeding based on true seeds have the potential to revolutionize future potato breeding and production^1-4. So far, relatively few studies have examined the genome evolution and diversity of wild and cultivated landrace potatoes, which limits the application of their diversity in potato breeding. Here we assemble 44 high-quality diploid potato genomes from 24 wild and 20 cultivated accessions that are representative of Solanum section Petota, the tuber-bearing clade, as well as 2 genomes from the neighbouring section, Etuberosum. Extensive discordance of phylogenomic relationships suggests the complexity of potato evolution. We find that the potato genome substantially expanded its repertoire of disease-resistance genes when compared with closely related seed-propagated solanaceous crops, indicative of the effect of tuber-based propagation strategies on the evolution of the potato genome. We discover a transcription factor that determines tuber identity and interacts with the mobile tuberization inductive signal SP6A. We also identify 561,433 high-confidence structural variants and construct a map of large inversions, which provides insights for improving inbred lines and precluding potential linkage drag, as exemplified by a 5.8-Mb inversion that is associated with carotenoid content in tubers. This study will accelerate hybrid potato breeding and enrich our understanding of the evolution and biology of potato as a global staple food crop.

Allelic variation for alpha-Glucan Water Dikinase is associated with starch phosphate content in tetraploid potato

Association analysis resulted in the identification of specific StGWD alleles causing either an increase or decrease in starch phosphate content which was verified in diploid and tetraploid potato mapping populations. Potatoes are grown for various purposes like French fries, table potatoes, crisps and for their starch. One of the most important aspects of potato starch is that it contains a high amount of phosphate ester groups which are considered to be important for providing improved functionalization after derivatization processes. Little is known about the variation in phosphate content as such in different potato varieties and thus we studied the genetic diversity for this trait. From other studies it was clear that the phosphate content is controlled by a quantitative trait locus (QTL) underlying the candidate gene α-Glucan Water Dikinase (StGWD) on chromosome 5. We performed direct amplicon sequencing of this gene by Sanger sequencing. Sequences of two StGWD amplicons from a global collection of 398 commercial cultivars and progenitor lines were used to identify 16 different haplotypes. By assigning tag SNPs to these haplotypes, each of the four alleles present in a cultivar could be deduced and linked to a phosphate content. A high value for intra-individual heterozygosity was observed (Ho = 0.765). The average number of different haplotypes per individual (Ai) was 3.1. Pedigree analysis confirmed that the haplotypes are identical-by-descent (IBD) and offered insight in the breeding history of elite potato germplasm. Haplotypes originating from introgression of wild potato accessions carrying resistance genes could be traced. Furthermore, association analysis resulted in the identification of specific StGWD alleles causing either an increase or decrease in starch phosphate content varying from 12 nmol PO₄/mg starch to 38 nmol PO₄/mg starch. These allele effects were verified in diploid and tetraploid mapping populations and offer possibilities to breed and select for this trait.

EMBRYO SAC DEVELOPMENT 1 affects seed setting rate in rice by controlling embryo sac development

Seed setting rate is one of the critical factors that determine rice yield. Grain formation is a complex biological process, whose molecular mechanism is yet to be improved. Here we investigated the function of an OVATE family protein, Embryo Sac Development 1 (ESD1), in the regulation of seed setting rate in rice (Oryza sativa) by examining its loss-of-function mutants generated via clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated9 (Cas9) technology. ESD1 was predominantly expressed at Stage 6 of panicle development, especially in the ovules. esd1 mutants displayed reduced seed setting rates with normal stamen development and pollen tube growth but abnormal pistil group. Investigation of embryo sacs revealed that during the mitosis of functional megaspores, some egg cells degraded during differentiation in esd1 mutants, thereby hindering subsequent fertilization process and reducing seed setting rate. In addition, the transcriptional level of O. sativa anaphase-promoting complex 6, a reported embryo sac developing gene, was significantly reduced in esd1 mutants. These results support that ESD1 is an important modulator of ESD and seed setting rate in rice. Together, this finding demonstrates that ESD1 positively regulates the seed setting rate by controlling ESD in rice and has implications for the improvement of rice yield.

A high-quality genome assembly of the eggplant provides insights into the molecular basis of disease resistance and chlorogenic acid synthesis

The eggplant (Solanum melongena L.) is one of the most important Solanaceae crops, ranking third for total production and economic value in its genus. Herein, we report a high-quality, chromosome-scale eggplant reference genome sequence of 1155.8 Mb, with an N50 of 93.9 Mb, which was assembled by combining PacBio long reads and Hi-C sequencing data. Repetitive sequences occupied 70.1% of the assembly length, and 35,018 high-confidence protein-coding genes were annotated based on multiple sources. Comparative analysis revealed 646 species-specific families and 364 positive selection genes, conferring distinguishing traits on the eggplant. We performed genome-wide comparative identification of disease resistance genes and discovered an expanded gene family of bacterial spot resistance in eggplant and pepper, but not in tomato and potato. The genes involved in chlorogenic acid synthesis were comprehensively characterized. Highly similar chromosomal distribution patterns of polyphenol oxidase genes were observed in the eggplant, tomato, and potato genomes. The eggplant reference genome sequence will not only facilitate evolutionary studies of the Solanaceae but also facilitate their breeding and improvement.

Genome-Wide Identification and Expression Analysis of the Dof Gene Family in Medicago sativa L. Under Various Abiotic Stresses

The Dof transcription factor is a plant-specific transcriptional regulator that plays important roles in plant development and acts as a mediator in plant external stress responses. However, Dofs have previously been identified in several plants but not in alfalfa (Medicago sativa L.), one of the most widely cultivated forage legumes. In the present study, a total of 40 MsDof genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, and expression patterns under abscisic acid (ABA), cold, heat, drought and salt stresses of these Dof genes were comprehensively analyzed. The Dof genes family in alfalfa could be classified into eight classes. Gene ontology (GO) and tissue-specific analysis indicated that most MsDof genes may be involved in biological functions during plant growth. Moreover, the expression profiles and quantitative real-time PCR analysis indicated that eight candidate abiotic tolerance genes were induced in response to four abiotic stresses. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole genome level, which would provide new information on the Dof family in alfalfa.

Evolution of a guarded decoy protease and its receptor in solanaceous plants

Rcr3 is a secreted protease of tomato that is targeted by fungal effector Avr2, a secreted protease inhibitor of the fungal pathogen Cladosporium fulvum. The Avr2-Rcr3 complex is recognized by receptor-like protein Cf-2, triggering hypersensitive cell death (HR) and disease resistance. Avr2 also targets Rcr3 paralog Pip1, which is not required for Avr2 recognition but contributes to basal resistance. Thus, Rcr3 acts as a guarded decoy in this interaction, trapping the fungus into a recognition event. Here we show that Rcr3 evolved > 50 million years ago (Mya), whereas Cf-2 evolved <6Mya by co-opting the pre-existing Rcr3 in the Solanum genus. Ancient Rcr3 homologs present in tomato, potato, eggplants, pepper, petunia and tobacco can be inhibited by Avr2 with the exception of tobacco Rcr3. Four variant residues in Rcr3 promote Avr2 inhibition, but the Rcr3 that co-evolved with Cf-2 lacks three of these residues, indicating that the Rcr3 co-receptor is suboptimal for Avr2 binding. Pepper Rcr3 triggers HR with Cf-2 and Avr2 when engineered for enhanced inhibition by Avr2. Nicotiana benthamiana (Nb) is a natural null mutant carrying Rcr3 and Pip1 alleles with deleterious frame-shift mutations. Resurrected NbRcr3 and NbPip1 alleles were active proteases and further NbRcr3 engineering facilitated Avr2 inhibition, uncoupled from HR signalling. The evolution of a receptor co-opting a conserved pathogen target contrasts with other indirect pathogen recognition mechanisms.

A universal chromosome identification system for maize and wild Zea species

Maize was one of the first eukaryotic species in which individual chromosomes can be identified cytologically, which made maize one of the oldest models for genetics and cytogenetics research. Nevertheless, consistent identification of all 10 chromosomes from different maize lines as well as from wild Zea species remains a challenge. We developed a new technique for maize chromosome identification based on fluorescence in situ hybridization (FISH). We developed two oligonucleotide-based probes that hybridize to 24 chromosomal regions. Individual maize chromosomes show distinct FISH signal patterns, which allow universal identification of all chromosomes from different Zea species. We developed karyotypes from three Zea mays subspecies and two additional wild Zea species based on individually identified chromosomes. A paracentric inversion was discovered on the long arm of chromosome 4 in Z. nicaraguensis and Z. luxurians based on modifications of the FISH signal patterns. Chromosomes from these two species also showed distinct distribution patterns of terminal knobs compared with other Zea species. These results support that Z. nicaraguensis and Z. luxurians are closely related species.

Evolutionary Rate Heterogeneity and Functional Divergence of Orthologous Genes in Pyrus

Negatively selected genes (NSGs) and positively selected genes (PSGs) are the two types of most nuclear protein-coding genes in organisms. However, the evolutionary rates and characteristics of different types of genes have been rarely understood. In the present study, we investigate the rates of synonymous substitution (Ks) and the rates of non-synonymous substitution (Ka) by comparing the orthologous genes of two sequenced Pyrus species, Pyrus bretschneideri and Pyrus communis. Subsequently, we compared the evolutionary rates, gene structures, and expression profiles during different fruit development between PSGs and NSGs. Compared with the NSGs, the PSGs have fewer exons, shorter gene length, lower synonymous substitution rates and have higher evolutionary rates. Remarkably, gene expression patterns between two Pyrus species fruit indicated functional divergence for most of the orthologous genes derived from a common ancestor, and subfunctionalization for some of them. Overall, the present study shows that PSGs differs from NSGs not only under environmental selective pressure (Ka/Ks), but also in their structural, functional, and evolutionary properties. Additionally, our resulting data provides important insights for the evolution and highlights the diversification of orthologous genes in two Pyrus species.

Environmental and epigenetic regulation of Rider retrotransposons in tomato

Transposable elements in crop plants are the powerful drivers of phenotypic variation that has been selected during domestication and breeding programs. In tomato, transpositions of the LTR (long terminal repeat) retrotransposon family Rider have contributed to various phenotypes of agronomical interest, such as fruit shape and colour. However, the mechanisms regulating Rider activity are largely unknown. We have developed a bioinformatics pipeline for the functional annotation of retrotransposons containing LTRs and defined all full-length Rider elements in the tomato genome. Subsequently, we showed that accumulation of Rider transcripts and transposition intermediates in the form of extrachromosomal DNA is triggered by drought stress and relies on abscisic acid signalling. We provide evidence that residual activity of Rider is controlled by epigenetic mechanisms involving siRNAs and the RNA-dependent DNA methylation pathway. Finally, we demonstrate the broad distribution of Rider-like elements in other plant species, including crops. Our work identifies Rider as an environment-responsive element and a potential source of genetic and epigenetic variation in plants.

Transposons are major constituents of plant genomes and represent a powerful source of internal genetic and epigenetic variation. For example, domestication of maize has been facilitated by a dramatic change in plant architecture, the consequence of a transposition event. Insertion of transposons near genes often confers quantitative phenotypic variation linked to changes in transcriptional patterns, as documented for blood oranges and grapes. In tomato, the most widely grown fruit crop and model for fleshy fruit biology, occurrences of several beneficial traits related to fruit shape and plant architecture are due to the activity of the transposon family Rider. While Rider represents a unique endogenous source of genetic and epigenetic variation, mechanisms regulating Rider activity remain unexplored. By achieving experimentally-controlled activation of the Rider family, we shed light on the regulation of these transposons by drought stress, signalling by phytohormones, as well as epigenetic pathways. Furthermore, we reveal the presence of Rider-like elements in other economically important crops such as rapeseed, beetroot and quinoa. This suggests that drought-inducible Rider activation could be further harnessed to generate genetic and epigenetic variation for crop breeding, and highlights the potential of transposon-directed mutagenesis for crop improvement.

Nucleotide sequence polymorphism in the RFL-PPR genes of potato

Cytoplasmic male sterility (CMS) is widely used for hybrid seed production in cultivated Solanaceae species. However, there is very limited information about CMS-Rf genetic systems in potato (Solanum tuberosum). Studying the CMS-Rf systems in potato is both of theoretical and practical significance due to the emergence of a new revolutionary strategy of reinventing potato as adiploid inbred line-based crop to develop F₁ hybrid seed potato breeding (Lindhout et al. 2011; Jansky et al. 2016). To search for potato Rf gene candidates, the comparative genetic approach was applied. Based on similarity to petunia Rf-PPR592 gene, 38 fragments were identified in five loci of the whole-genome nucleotide sequence of the accession DM 1-3 516 R44 S. tuberosum Phureja group (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The putative encoded mitochondrial proteins have 589-597 amino acid residues, similarto RF-PPR proteins of petunia and chili pepper and contain 14 or 15 PPR motifs. Primers have been developed flanking the most variable 782-865 bp regions of the selected loci, and polymorphism of the cloned fragments has been investigated in a subset of nine potato genotypes. The amplified fragments included seven or eight PPR motifs and lacked introns. The SNP frequencies ranged from 7.0 to 19.8% depending on the locus, while the ratio of nonsynonymous to synonymous substitutions varied between 0.9 and 2.1. Positions 1, 3 and 6 were the most variable in the studied PPR motifs. Our results demonstrated that the analysed sequences belong to the RFL-PPR gene subfamily and may be considered as Rf gene candidates in potato.

Comparative Transcriptome Analysis in Eggplant Reveals Selection Trends during Eggplant Domestication

Eggplant (Solanum melongena L.) is an economically and nutritionally important fruit crop of the Solanaceae family, which was domesticated in India and southern China. However, the genome regions subjected to selective sweeps in eggplant remain unknown. In the present study, we performed comparative transcriptome analysis of cultivated and wild eggplant species with emphasis on the selection pattern during domestication. In total, 44,073 (S. sisymbriifolium) to 58,677 (S. melongena cultivar S58) unigenes were generated for the six eggplant accessions with total lengths of 36.6-46 Mb. The orthologous genes were assessed using the ratio of nonsynonymous (K a) to synonymous (K s) nucleotide substitutions to characterize selective patterns during eggplant domestication. We identified 19 genes under positive selection across the phylogeny that were classified into four groups. The gene (OG12205) under positive selection was possibly associated with fruit-related traits in eggplant, which may have resulted from human manipulation. Eight positive selected genes were potentially involved in stress tolerance or disease resistance, suggesting that environmental changes and biotic stresses were important selective pressures in eggplant domestication. Taken together, our results shed light on the effects of artificial and natural selection on the transcriptomes of eggplant and its wild relatives. Identification of the selected genes will facilitate the understanding of genetic architecture of domesticated-related traits and provide resources for resistant breeding in eggplant.

Association Study Reveals Novel Genes Related to Yield and Quality of Fruit in Cape Gooseberry (Physalis peruviana L.)

Association mapping has been proposed as an efficient approach to assist plant breeding programs to investigate the genetic basis of agronomic traits. In this study, we evaluated 18 traits related to yield, (FWP, NF, FWI, and FWII), fruit size-shape (FP, FA, MW, WMH, MH, HMW, DI, FSI, FSII, OVO, OBO), and fruit quality (FIR, CF, and SST), in a diverse collection of 100 accessions of Physalis peruviana including wild, landraces, and anther culture derived lines. We identified seven accessions with suitable traits: fruit weight per plant (FWP) > 7,000 g/plant and cracked fruits (CF) < 4%, to be used as parents in cape gooseberry breeding program. In addition, the accessions were also characterized using Genotyping By Sequencing (GBS). We discovered 27,982 and 36,142 informative SNP markers based on the alignment against the two cape gooseberry references transcriptomes. Besides, 30,344 SNPs were identified based on alignment to the tomato reference genome. Genetic structure analysis showed that the population could be divided into two or three sub-groups, corresponding to landraces-anther culture and wild accessions for K = 2 and wild, landraces, and anther culture plants for K = 3. Association analysis was carried out using a Mixed Linear Model (MLM) and 34 SNP markers were significantly associated. These results reveal the basis of the genetic control of important agronomic traits and may facilitate marker-based breeding in P. peruviana.

Comparative Oligo-FISH Mapping: An Efficient and Powerful Methodology To Reveal Karyotypic and Chromosomal Evolution

Developing the karyotype of a eukaryotic species relies on identification of individual chromosomes, which has been a major challenge for most nonmodel plant and animal species. We developed a novel chromosome identification system by selecting and labeling oligonucleotides (oligos) located in specific regions on every chromosome. We selected a set of 54,672 oligos (45 nt) based on single copy DNA sequences in the potato genome. These oligos generated 26 distinct FISH signals that can be used as a "bar code" or "banding pattern" to uniquely label each of the 12 chromosomes from both diploid and polyploid (4× and 6×) potato species. Remarkably, the same bar code can be used to identify the 12 homeologous chromosomes among distantly related Solanum species, including tomato and eggplant. Accurate karyotypes based on individually identified chromosomes were established in six Solanum species that have diverged for >15 MY. These six species have maintained a similar karyotype; however, modifications to the FISH signal bar code led to the discovery of two reciprocal chromosomal translocations in Solanum etuberosum and S. caripense We also validated these translocations by oligo-based chromosome painting. We demonstrate that the oligo-based FISH techniques are powerful new tools for chromosome identification and karyotyping research, especially for nonmodel plant species.

Parallel evolution of the POQR prolyl oligo peptidase gene conferring plant quantitative disease resistance

Plant pathogens with a broad host range are able to infect plant lineages that diverged over 100 million years ago. They exert similar and recurring constraints on the evolution of unrelated plant populations. Plants generally respond with quantitative disease resistance (QDR), a form of immunity relying on complex genetic determinants. In most cases, the molecular determinants of QDR and how they evolve is unknown. Here we identify in Arabidopsis thaliana a gene mediating QDR against Sclerotinia sclerotiorum, agent of the white mold disease, and provide evidence of its convergent evolution in multiple plant species. Using genome wide association mapping in A. thaliana, we associated the gene encoding the POQR prolyl-oligopeptidase with QDR against S. sclerotiorum. Loss of this gene compromised QDR against S. sclerotiorum but not against a bacterial pathogen. Natural diversity analysis associated POQR sequence with QDR. Remarkably, the same amino acid changes occurred after independent duplications of POQR in ancestors of multiple plant species, including A. thaliana and tomato. Genome-scale expression analyses revealed that parallel divergence in gene expression upon S. sclerotiorum infection is a frequent pattern in genes, such as POQR, that duplicated both in A. thaliana and tomato. Our study identifies a previously uncharacterized gene mediating QDR against S. sclerotiorum. It shows that some QDR determinants are conserved in distantly related plants and have emerged through the repeated use of similar genetic polymorphisms at different evolutionary time scales.

Conservation and Dispersion of Genes Conferring Resistance to Tomato Begomoviruses between Tomato and Pepper Genomes

In the present climate change scenario, controlling plant disease through exploitation of host plant resistance could contribute toward the sustainable crop production and global food security. In this respect, the identification of new sources of resistance and utilization of genetic diversity within the species may help in the generation of cultivars with improved disease resistance. Begomoviruses namely, Tomato yellow leaf curl virus (TYLCV) and Chilli leaf curl virus (ChLCV) are known to cause major yield losses in several economically important crop plants of the family Solanaceae. Though co-occurrence, association and synergistic interactions among these viruses in the host plants is reported, whether orthologous genetic loci in related host plants could be responsible for conferring resistance to these viruses has not been investigated yet. Several loci including Ty1, Ty2, Ty3, Ty4, and ty5 have been reported to confer resistance to leaf curl viruses in tomato. Here, we examined the pepper orthologous markers, corresponding to these QTL regions, for polymorphism between ChLCV susceptible and resistant genotypes of pepper. Further, to examine if the polymorphic markers are segregating with the disease resistance, Bulk Segregant Analysis (BSA) was performed on F₂ population derived from crosses between resistant and susceptible lines. However, none of the markers showed polymorphism in BSA suggesting that the tested markers are not linked to genes/QTLs responsible for conferring resistance to ChLCV in the selected genotypes. In silico analysis was performed to study the synteny and collinearity of genes located within these QTL regions in tomato and pepper genomes, which revealed that more than 60% genes located in Ty2 and Ty4, 13.71% genes in Ty1, 23.07% in Ty3, and 44.77% genes located within ty5 QTL region in tomato are conserved in pepper genome. However, despite such a high conservation in gene content, the linkage relationship in these regions seems to be greatly affected by gross rearrangements in both the species.

Microbial taxonomy in the era of OMICS: application of DNA sequences, computational tools and techniques

The current prokaryotic taxonomy classifies phenotypically and genotypically diverse microorganisms using a polyphasic approach. With advances in the next-generation sequencing technologies and computational tools for analysis of genomes, the traditional polyphasic method is complemented with genomic data to delineate and classify bacterial genera and species as an alternative to cumbersome and error-prone laboratory tests. This review discusses the applications of sequence-based tools and techniques for bacterial classification and provides a scheme for more robust and reproducible bacterial classification based on genomic data. The present review highlights promising tools and techniques such as ortho-Average Nucleotide Identity, Genome to Genome Distance Calculator and Multi Locus Sequence Analysis, which can be validly employed for characterizing novel microorganisms and assessing phylogenetic relationships. In addition, the review discusses the possibility of employing metagenomic data to assess the phylogenetic associations of uncultured microorganisms. Through this article, we present a review of genomic approaches that can be included in the scheme of taxonomy of bacteria and archaea based on computational and in silico advances to boost the credibility of taxonomic classification in this genomic era.

Dissecting Tissue-Specific Transcriptomic Responses from Leaf and Roots under Salt Stress in Petunia hybrida Mitchell

One of the primary objectives of plant biotechnology is to increase resistance to abiotic stresses, such as salinity. Salinity is a major abiotic stress and increasing crop resistant to salt continues to the present day as a major challenge. Salt stress disturbs cellular environment leading to protein misfolding, affecting normal plant growth and causing agricultural losses worldwide. The advent of state-of-the-art technologies such as high throughput mRNA sequencing (RNA-seq) has revolutionized whole-transcriptome analysis by allowing, with high precision, to measure changes in gene expression. In this work, we used tissue-specific RNA-seq to gain insight into the Petunia hybrida transcriptional responses under NaCl stress using a controlled hydroponic system. Roots and leaves samples were taken from a continuum of 48 h of acute 150 mM NaCl. This analysis revealed a set of tissue and time point specific differentially expressed genes, such as genes related to transport, signal transduction, ion homeostasis as well as novel and undescribed genes, such as Peaxi162Scf00003g04130 and Peaxi162Scf00589g00323 expressed only in roots under salt stress. In this work, we identified early and late expressed genes in response to salt stress while providing a core of differentially express genes across all time points and tissues, including the trehalose-6-phosphate synthase 1 (TPS1), a glycosyltransferase reported in salt tolerance in other species. To test the function of the novel petunia TPS1 allele, we cloned and showed that TPS1 is a functional plant gene capable of complementing the trehalose biosynthesis pathway in a yeast tps1 mutant. The list of candidate genes to enhance salt tolerance provided in this work constitutes a major effort to better understand the detrimental effects of salinity in petunia with direct implications for other economically important Solanaceous species.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

Distinct Patterns of Gene Gain and Loss: Diverse Evolutionary Modes of NBS-Encoding Genes in Three Solanaceae Crop Species

Plant resistance conferred by nucleotide binding site (NBS)-encoding resistance genes plays a key role in the defense against various pathogens throughout the entire plant life cycle. However, comparative analyses for the systematic evaluation and determination of the evolutionary modes of NBS-encoding genes among Solanaceae species are rare. In this study, 447, 255, and 306 NBS-encoding genes were identified from the genomes of potato, tomato, and pepper, respectively. These genes usually clustered as tandem arrays on chromosomes; few existed as singletons. Phylogenetic analysis indicated that three subclasses [TNLs (TIR-NBS-LRR), CNLs (CC-NBS-LRR), and RNLs (RPW8-NBS-LRR)] each formed a monophyletic clade and were distinguished by unique exon/intron structures and amino acid motif sequences. By comparing phylogenetic and systematic relationships, we inferred that the NBS-encoding genes in the present genomes of potato, tomato, and pepper were derived from 150 CNL, 22 TNL, and 4 RNL ancestral genes, and underwent independent gene loss and duplication events after speciation. The NBS-encoding genes therefore exhibit diverse and dynamic evolutionary patterns in the three Solanaceae species, giving rise to the discrepant gene numbers observed today. Potato shows a “consistent expansion” pattern, tomato exhibits a pattern of “first expansion and then contraction,” and pepper presents a “shrinking” pattern. The earlier expansion of CNLs in the common ancestor led to the dominance of this subclass in gene numbers. However, RNLs remained at low copy numbers due to their specific functions. Along the evolutionary process of NBS-encoding genes in Solanaceae, species-specific tandem duplications contributed the most to gene expansions.

Fast, scalable prediction of deleterious noncoding variants from functional and population genomic data

Many genetic variants that influence phenotypes of interest are located outside of protein-coding genes, yet existing methods for identifying such variants have poor predictive power. Here we introduce a new computational method, called LINSIGHT, that substantially improves the prediction of noncoding nucleotide sites at which mutations are likely to have deleterious fitness consequences, and which, therefore, are likely to be phenotypically important. LINSIGHT combines a generalized linear model for functional genomic data with a probabilistic model of molecular evolution. The method is fast and highly scalable, enabling it to exploit the 'big data' available in modern genomics. We show that LINSIGHT outperforms the best available methods in identifying human noncoding variants associated with inherited diseases. In addition, we apply LINSIGHT to an atlas of human enhancers and show that the fitness consequences at enhancers depend on cell type, tissue specificity, and constraints at associated promoters.

Phytochrome Interacting Factors (PIFs) in Solanum lycopersicum: Diversity, Evolutionary History and Expression Profiling during Different Developmental Processes

Although the importance of light for tomato plant yield and edible fruit quality is well known, the PHYTOCHROME INTERACTING FACTORS (PIFs), main components of phytochrome-mediated light signal transduction, have been studied almost exclusively in Arabidopsis thaliana. Here, the diversity, evolution and expression profile of PIF gene subfamily in Solanum lycopersicum was characterized. Eight tomato PIF loci were identified, named SlPIF1a, SlPIF1b, SlPIF3, SlPIF4, SlPIF7a, SlPIF7b, SlPIF8a and SlPIF8b. The duplication of SlPIF1, SlPIF7 and SlPIF8 genes were dated and temporally coincided with the whole-genome triplication event that preceded tomato and potato divergence. Different patterns of mRNA accumulation in response to light treatments were observed during seedling deetiolation, dark-induced senescence, diel cycle and fruit ripening. SlPIF4 showed similar expression profile as that reported for A. thaliana homologs, indicating an evolutionary conserved function of PIF4 clade. A comprehensive analysis of the evolutionary and transcriptional data allowed proposing that duplicated SlPIFs have undergone sub- and neofunctionalization at mRNA level, pinpointing the importance of transcriptional regulation for the maintenance of duplicated genes. Altogether, the results indicate that genome polyploidization and functional divergence have played a major role in diversification of the Solanum PIF gene subfamily.

Genome-wide analysis of Dof transcription factors reveals functional characteristics during development and response to biotic stresses in pepper

The DNA-binding with one zinc finger proteins (Dofs) are a plant-specific family of transcription factors. The Dofs are involved in a variety of biological processes such as phytohormone production, seed development, and environmental adaptation. Dofs have been previously identified in several plants, but not in pepper. We identified 33 putative Dof genes in pepper (CaDofs). To gain an overview of the CaDofs, we analyzed phylogenetic relationships, protein motifs, and evolutionary history. We divided the 33 CaDofs, containing 25 motifs, into four major groups distributed on eight chromosomes. We discovered an expansion of the CaDofs dated to a recent duplication event. Segmental duplication that occurred before the speciation of the Solanaceae lineages was predominant among the CaDofs. The global gene-expression profiling of the CaDofs by RNA-seq analysis showed distinct temporal and pathogen-specific variation during development and response to biotic stresses (two TMV strains, PepMoV, and Phytophthora capsici), suggesting functional diversity among the CaDofs. These results will provide the useful clues into the responses of Dofs in biotic stresses and promote a better understanding of their multiple function in pepper and other species.

The first de novo transcriptome of pepino (Solanum muricatum): assembly, comprehensive analysis and comparison with the closely related species S. caripense, potato and tomato

BACKGROUND

Solanum sect. Basarthrum is phylogenetically very close to potatoes (Solanum sect. Petota) and tomatoes (Solanum sect. Lycopersicon), two groups with great economic importance, and for which Solanum sect. Basarthrum represents a tertiary gene pool for breeding. This section includes the important regional cultigen, the pepino (Solanum muricatum), and several wild species. Among the wild species, S. caripense is prominent due to its major involvement in the origin of pepino and its wide geographical distribution. Despite the value of the pepino as an emerging crop, and the potential for gene transfer from both the pepino and S. caripense to potatoes and tomatoes, there has been virtually no genomic study of these species.

RESULTS

Using Illumina HiSeq 2000, RNA-Seq was performed with a pool of three tissues (young leaf, flowers in pre-anthesis and mature fruits) from S. muricatum and S. caripense, generating almost 111,000,000 reads among the two species. A high quality de novo transcriptome was assembled from S. muricatum clean reads resulting in 75,832 unigenes with an average length of 704 bp. These unigenes were functionally annotated based on similarity of public databases. We used Blast2GO, to conduct an exhaustive study of the gene ontology, including GO terms, EC numbers and KEGG pathways. Pepino unigenes were compared to both potato and tomato genomes in order to determine their estimated relative position, and to infer gene prediction models. Candidate genes related to traits of interest in other Solanaceae were evaluated by presence or absence and compared with S. caripense transcripts. In addition, by studying five genes, the phylogeny of pepino and five other members of the family, Solanaceae, were studied. The comparison of S. caripense reads against S. muricatum assembled transcripts resulted in thousands of intra- and interspecific nucleotide-level variants. In addition, more than 1000 SSRs were identified in the pepino transcriptome.

CONCLUSIONS

This study represents the first genomic resource for the pepino. We suggest that the data will be useful not only for improvement of the pepino, but also for potato and tomato breeding and gene transfer. The high quality of the transcriptome presented here also facilitates comparative studies in the genus Solanum. The accurate transcript annotation will enable us to figure out the gene function of particular traits of interest. The high number of markers (SSR and nucleotide-level variants) obtained will be useful for breeding programs, as well as studies of synteny, diversity evolution, and phylogeny.

Identification of key genes involved in root development of tomato using expressed sequence tag analysis

Root system of plants are actually fascinating structures, not only critical for plant development, but also important for storage and conduction. Due to its agronomic importance, identification of genes involved in root development has been a subject of intense study. Tomato is the one of the most consumed vegetables in the world. Tomato has been used as model system for dicot plants because of its small genome, well-established transformation techniques and well-constructed physical map. The present study is targeted to identify of root specific genes expressed temporally and also gene(s) involved in lateral root and profuse root development. A total of 890 ESTs were identified from five EST libraries constructed using SSH approach which included temporal gene regulation (early and late) and genes involved in morphogenetic traits (lateral and profuse rooting). One hundred sixty-one unique ESTs identified from various libraries were categorized based on their putative functions and deposited in NCBI-dbEST database. In addition, 36 ESTs were selected for validation of their expression by RT-PCR. The present findings will help in shedding light to the unexplored developmental process of root growth in tomato and plant in general.

Identification and phylogenetic analysis of late embryogenesis abundant proteins family in tomato (Solanum lycopersicum)

This study provided a comparative genomic analysis of the LEA gene family, and these may provide valuable information for their functional investigations in the future. Late embryogenesis abundant (LEA) proteins are a group of proteins that accumulate in response to cellular dehydration in many organisms. Here, we identified 27 LEA genes in tomato. A strong correlation between phylogeny, gene structure, and motif composition was found. The predicted SlLEA genes were non-randomly distributed within their chromosomes, and segmental and tandem duplications were probably important for their expansion. Many cis-elements potentially mediating transcription in response to abiotic stress were also found in the 1,000 bp upstream sequence of the promoter region. An additional intragenic recombination played an important role in the evolution of SlLEA genes. Selection analysis also identified some significant site-specific constraints that acted on the evolution of most LEA paralogs. Expression analysis using both microarray data and quantitative real-time PCR indicated that SlLEA genes were widely expressed in various tissues, and that a few members responded to some abiotic stresses. Our study provides useful information on the LEA genes in tomato and will facilitate their further characterization to better understand their functions.

Physalis floridana Cell Number Regulator1 encodes a cell membrane-anchored modulator of cell cycle and negatively controls fruit size

Physalis species show a significant variation in berry size; however, the underlying molecular basis is unknown. In this work, we showed that cell division difference in the ovaries might contribute to the ultimate berry size variation within Physalis species, and that mRNA abundance of Physalis floridana Cell Number Regulator1 (PfCNR1), the putative orthologue of the tomato fruit weight 2.2 (FW2.2), was negatively correlated with cell division in the ovaries. Moreover, heterochronic expression variation of the PfCNR1 genes in the ovaries concomitantly correlated with berry weight variation within Physalis species. In transgenic Physalis, multiple organ sizes could be negatively controlled by altering PfCNR1 levels, and cell division instead of cell expansion was primarily affected. PfCNR1 was shown to be anchored in the plasma membrane and to interact with PfAG2 (an AGAMOUS-like protein determining ovary identity). The expression of PfCYCD2;1, a putative orthologue of the mitosis-specific gene CyclinD2;1 in the cell cycle was negatively correlated with the PfCNR1 mRNA levels. PfAG2 was found to selectively bind to the CArG-box in the PfCYCD2;1 promoter and to repress PfCYCD2;1 expression, thus suggesting a PfAG2-mediated pathway for PfCNR1 to regulate cell division. The interaction of PfCNR1 with PfAG2 enhanced the repression of PfCYCD2;1 expression. The nuclear import of PfAG2 was essential in the proposed pathway. Our data provide new insights into the developmental pathways of a cell membrane-anchored protein that modulates cell division and governs organ size determination. This study also sheds light on the link between organ identity and organ growth in plants.

Development of fruit color in Solanaceae: a story of two biosynthetic pathways

This review highlights the major differences between the regulation of two important pathways namely anthocyanin and carotenoid pathways, responsible for fruit color generation in Solanaceae mediated by transcription factors (TFs). The anthocyanin pathway is regulated by a common set of TFs (MYB, MYC and WD40) belonging to specific families of DNA-binding proteins. Their regulation is aimed at controlling the type and amount of pigments produced and the physiological conditions (like pH) at which they are finally stored. In the carotenoid pathway, the color diversity depends on the quantity of pigment produced and the point where the pathway is arrested. TFs in the latter case are accordingly found to influence the sequestration and degradation of these pigments, which determines their final concentration in the tissue. TFs (phytochrome interacting factors, MADS-BOX, HB-ZIP and B-ZIP) also regulate important rate-determining steps, which decide the direction in which the pathway proceeds and the point at which it is terminated. In the absence of a clear pattern of TF-mediated regulation, it is suggested that the carotenoid pathway is more significantly influenced by other regulatory methods which need to be explored. It is expected that common factors affecting these pathways are the ones acting much before the initiation of the biosynthesis of respective pigments.

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

BACKGROUND AND AIMS

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

METHODS

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

KEY RESULTS

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

CONCLUSIONS

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.

Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis

BACKGROUND

Eggplant (Solanum melongena L.) and turkey berry (S. torvum Sw.), a wild ally of eggplant with promising multi-disease resistance traits, are of great economic, medicinal and genetic importance, but genomic resources for these species are lacking. In the present study, we sequenced the transcriptomes of eggplant and turkey berry to accelerate research on these two non-model species.

RESULTS

We built comprehensive, high-quality de novo transcriptome assemblies of the two Leptostemonum clade Solanum species from short-read RNA-Sequencing data. We obtained 34,174 unigenes for eggplant and 38,185 unigenes for turkey berry. Functional annotations based on sequence similarity to known plant datasets revealed a distribution of functional categories for both species very similar to that of tomato. Comparison of eggplant, turkey berry and another 11 plant proteomes resulted in 276 high-confidence single-copy orthologous groups, reasonable phylogenetic tree inferences and reliable divergence time estimations. From these data, it appears that eggplant and its wild Leptostemonum clade relative turkey berry split from each other in the late Miocene, ~6.66 million years ago, and that Leptostemonum split from the Potatoe clade in the middle Miocene, ~15.75 million years ago. Furthermore, 621 and 815 plant resistance genes were identified in eggplant and turkey berry respectively, indicating the variation of disease resistance genes between them.

CONCLUSIONS

This study provides a comprehensive transcriptome resource for two Leptostemonum clade Solanum species and insight into their evolutionary history and biological characteristics. These resources establish a foundation for further investigations of eggplant biology and for agricultural improvement of this important vegetable. More generally, we show that RNA-Seq is a fast, reliable and cost-effective method for assessing genome evolution in non-model species.

A comprehensive phylogenetic and structural analysis of the carcinoembryonic antigen (CEA) gene family

The carcinoembryonic antigen (CEA) gene family belongs to the immunoglobulin (Ig) superfamily and codes for a vast number of glycoproteins that differ greatly both in amino acid composition and function. The CEA family is divided into two groups, the carcinoembryonic antigen-related cell adhesion molecules (CEACAMs) and the pregnancy-specific glycoproteins. The CEA family members are implicated in pleiotropic (patho)physiological functions including cell-cell adhesion, pregnancy, immunity, neovascularization, regulation of insulin homeostasis, and carcinogenesis. In general, the CEA-encoded proteins are composed of an extracellular region with Ig variable and constant-like domains and a cytoplasmic region containing signaling motifs. Of particular interest, the well-studied human and mouse CEA genes are arranged in clusters in a single chromosome. Taking into account this characteristic, we made an effort to reconstruct the evolutionary history of the CEA gene family. Toward this end, the publicly available genomes were searched extensively for CEA homologs. The domain organization of the retrieved protein sequences was analyzed, and, subsequently, comprehensive phylogenetic analyses of the entire length CEA homologous proteins were performed. A series of evolutionarily conserved amino acid residues, functionally important, were identified. The relative positioning of these residues on the modeled tertiary structure of novel CEA protein domains revealed that they are, also, spatially conserved. Furthermore, the chromosomal arrangement of CEA genes was examined, and it was found that the CEA genes are preserved in terms of position, transcriptional orientation, and number in all species under investigation.

A phylogenetic framework for evolutionary study of the nightshades (Solanaceae): a dated 1000-tip tree

BACKGROUND

The Solanaceae is a plant family of great economic importance. Despite a wealth of phylogenetic work on individual clades and a deep knowledge of particular cultivated species such as tomato and potato, a robust evolutionary framework with a dated molecular phylogeny for the family is still lacking. Here we investigate molecular divergence times for Solanaceae using a densely-sampled species-level phylogeny. We also review the fossil record of the family to derive robust calibration points, and estimate a chronogram using an uncorrelated relaxed molecular clock.

RESULTS

Our densely-sampled phylogeny shows strong support for all previously identified clades of Solanaceae and strongly supported relationships between the major clades, particularly within Solanum. The Tomato clade is shown to be sister to section Petota, and the Regmandra clade is the first branching member of the Potato clade. The minimum age estimates for major splits within the family provided here correspond well with results from previous studies, indicating splits between tomato and potato around 8 Million years ago (Ma) with a 95% highest posterior density (HPD) 7-10 Ma, Solanum and Capsicum c. 19 Ma (95% HPD 17-21), and Solanum and Nicotiana c. 24 Ma (95% HPD 23-26).

CONCLUSIONS

Our large time-calibrated phylogeny provides a significant step towards completing a fully sampled species-level phylogeny for Solanaceae, and provides age estimates for the whole family. The chronogram now includes 40% of known species and all but two monotypic genera, and is one of the best sampled angiosperm family phylogenies both in terms of taxon sampling and resolution published thus far. The increased resolution in the chronogram combined with the large increase in species sampling will provide much needed data for the examination of many biological questions using Solanaceae as a model system.

Identification of immunity related genes to study the Physalis peruviana--Fusarium oxysporum pathosystem

The Cape gooseberry (Physalisperuviana L) is an Andean exotic fruit with high nutritional value and appealing medicinal properties. However, its cultivation faces important phytosanitary problems mainly due to pathogens like Fusarium oxysporum, Cercosporaphysalidis and Alternaria spp. Here we used the Cape gooseberry foliar transcriptome to search for proteins that encode conserved domains related to plant immunity including: NBS (Nucleotide Binding Site), CC (Coiled-Coil), TIR (Toll/Interleukin-1 Receptor). We identified 74 immunity related gene candidates in P. peruviana which have the typical resistance gene (R-gene) architecture, 17 Receptor like kinase (RLKs) candidates related to PAMP-Triggered Immunity (PTI), eight (TIR-NBS-LRR, or TNL) and nine (CC–NBS-LRR, or CNL) candidates related to Effector-Triggered Immunity (ETI) genes among others. These candidate genes were categorized by molecular function (98%), biological process (85%) and cellular component (79%) using gene ontology. Some of the most interesting predicted roles were those associated with binding and transferase activity. We designed 94 primers pairs from the 74 immunity-related genes (IRGs) to amplify the corresponding genomic regions on six genotypes that included resistant and susceptible materials. From these, we selected 17 single band amplicons and sequenced them in 14 F. oxysporum resistant and susceptible genotypes. Sequence polymorphisms were analyzed through preliminary candidate gene association, which allowed the detection of one SNP at the PpIRG-63 marker revealing a nonsynonymous mutation in the predicted LRR domain suggesting functional roles for resistance.

Linking the potato genome to the conserved ortholog set (COS) markers

BACKGROUND

Conserved ortholog set (COS) markers are an important functional genomics resource that has greatly improved orthology detection in Asterid species. A comprehensive list of these markers is available at Sol Genomics Network (http://solgenomics.net/) and many of these have been placed on the genetic maps of a number of solanaceous species.

RESULTS

We amplified over 300 COS markers from eight potato accessions involving two diploid landraces of Solanum tuberosum Andigenum group (formerly classified as S. goniocalyx, S. phureja), and a dihaploid clone derived from a modern tetraploid cultivar of S. tuberosum and the wild species S. berthaultii, S. chomatophilum, and S. paucissectum. By BLASTn (Basic Local Alignment Search Tool of the NCBI, National Center for Biotechnology Information) algorithm we mapped the DNA sequences of these markers into the potato genome sequence. Additionally, we mapped a subset of these markers genetically in potato and present a comparison between the physical and genetic locations of these markers in potato and in comparison with the genetic location in tomato. We found that most of the COS markers are single-copy in the reference genome of potato and that the genetic location in tomato and physical location in potato sequence are mostly in agreement. However, we did find some COS markers that are present in multiple copies and those that map in unexpected locations. Sequence comparisons between species show that some of these markers may be paralogs.

CONCLUSIONS

The sequence-based physical map becomes helpful in identification of markers for traits of interest thereby reducing the number of markers to be tested for applications like marker assisted selection, diversity, and phylogenetic studies.

Genome Microscale Heterogeneity among Wild Potatoes Revealed by Diversity Arrays Technology Marker Sequences

Tuber-bearing potato species possess several genes that can be exploited to improve the genetic background of the cultivated potato Solanum tuberosum. Among them, S. bulbocastanum and S. commersonii are well known for their strong resistance to environmental stresses. However, scant information is available for these species in terms of genome organization, gene function, and regulatory networks. Consequently, genomic tools to assist breeding are meager, and efficient exploitation of these species has been limited so far. In this paper, we employed the reference genome sequences from cultivated potato and tomato and a collection of sequences of 1,423 potato Diversity Arrays Technology (DArT) markers that show polymorphic representation across the genomes of S. bulbocastanum and/or S. commersonii genotypes. Our results highlighted microscale genome sequence heterogeneity that may play a significant role in functional and structural divergence between related species. Our analytical approach provides knowledge of genome structural and sequence variability that could not be detected by transcriptome and proteome approaches.

Conservation and functional element discovery in 20 angiosperm plant genomes

Here, we describe the construction of a phylogenetically deep, whole-genome alignment of 20 flowering plants, along with an analysis of plant genome conservation. Each included angiosperm genome was aligned to a reference genome, Arabidopsis thaliana, using the LASTZ/MULTIZ paradigm and tools from the University of California-Santa Cruz Genome Browser source code. In addition to the multiple alignment, we created a local genome browser displaying multiple tracks of newly generated genome annotation, as well as annotation sourced from published data of other research groups. An investigation into A. thaliana gene features present in the aligned A. lyrata genome revealed better conservation of start codons, stop codons, and splice sites within our alignments (51% of features from A. thaliana conserved without interruption in A. lyrata) when compared with previous publicly available plant pairwise alignments (34% of features conserved). The detailed view of conservation across angiosperms revealed not only high coding-sequence conservation but also a large set of previously uncharacterized intergenic conservation. From this, we annotated the collection of conserved features, revealing dozens of putative noncoding RNAs, including some with recorded small RNA expression. Comparing conservation between kingdoms revealed a faster decay of vertebrate genome features when compared with angiosperm genomes. Finally, conserved sequences were searched for folding RNA features, including but not limited to noncoding RNA (ncRNA) genes. Among these, we highlight a double hairpin in the 5'-untranslated region (5'-UTR) of the PRIN2 gene and a putative ncRNA with homology targeting the LAF3 protein.

Tight genetic linkage of prezygotic barrier loci creates a multifunctional speciation island in Petunia

Most flowering plants depend on animal vectors for pollination and seed dispersal. Differential pollinator preferences lead to premating isolation and thus reduced gene flow between interbreeding plant populations. Sets of floral traits, adapted to attract specific pollinator guilds, are called pollination syndromes. Shifts in pollination syndromes have occurred surprisingly frequently, considering that they must involve coordinated changes in multiple genes affecting multiple floral traits. Although the identification of individual genes specifying single pollination syndrome traits is in progress in many species, little is known about the genetic architecture of coadapted pollination syndrome traits and how they are embedded within the genome. Here we describe the tight genetic linkage of loci specifying five major pollination syndrome traits in the genus Petunia: visible color, UV absorption, floral scent production, pistil length, and stamen length. Comparison with other Solanaceae indicates that, in P. exserta and P. axillaris, loci specifying these floral traits have specifically become clustered into a multifunctional "speciation island". Such an arrangement promotes linkage disequilibrium and avoids the dissolution of pollination syndromes by recombination. We suggest that tight genetic linkage provides a mechanism for rapid switches between distinct pollination syndromes in response to changes in pollinator availabilities.

Wild relatives of the eggplant (Solanum melongena L.: Solanaceae): new understanding of species names in a complex group

BACKGROUND

The common or brinjal eggplant (Solanum melongena L.) belongs to the Leptostemonum Clade (the "spiny" solanums) of the species-rich genus Solanum (Solanaceae). Unlike most of the genus, the eggplant and its relatives are from the Old World; most eggplant wild relatives are from Africa. An informal system for naming eggplant wild relatives largely based on crossing and other biosystematics data has been in use for approximately a decade. This system recognises several forms of two broadly conceived species, S. incanum L. and S. melongena. Recent morphological and molecular work has shown that species-level differences exist between these entities, and a new species-level nomenclature has been identified as necessary for plant breeders and for the maintenance of accurately named germplasm.

METHODOLOGY/PRINCIPAL FINDINGS

We examined herbarium specimens from throughout the wild species ranges as part of a larger revision of the spiny solanums of Africa. Based on these morphological and molecular studies, we delimited species in the group to which the common eggplant belongs and constructed identification keys for the group. We also examined the monophyly of the group considered as the eggplant relatives by previous authors.

CONCLUSIONS/SIGNIFICANCE

WE RECOGNISE TEN SPECIES IN THIS GROUP: S. aureitomentosum Bitter, S. campylacanthum A.Rich., S. cerasiferum Dunal, S. incanum L., S. insanum L., S. lichtensteinii Willd., S. linnaeanum Hepper & P.-M.L.Jaeger, S. melongena L., S. rigidum Lam. and S. umtuma Voronts. & S.Knapp. We review the history of naming and provide keys and character lists for all species. Ploidy level differences have not been investigated in the eggplant wild relatives; we identify this as a priority for improvement of crop wild relative use in breeding. The application of species-level names to these entities will help focus new collecting efforts for brinjal eggplant improvement and help facilitate information exchange.

Overview of tomato (Solanum lycopersicum) candidate pathogen recognition genes reveals important Solanum R locus dynamics

To investigate the genome-wide spatial arrangement of R loci, a complete catalogue of tomato (Solanum lycopersicum) and potato (Solanum tuberosum) nucleotide-binding site (NBS) NBS, receptor-like protein (RLP) and receptor-like kinase (RLK) gene repertories was generated. Candidate pathogen recognition genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution. NBS genes frequently occur in clusters of related gene copies that also include RLP or RLK genes. This scenario is compatible with the existence of selective pressures optimizing coordinated transcription. A number of duplication events associated with lineage-specific evolution were discovered. These findings suggest that different evolutionary mechanisms shaped pathogen recognition gene cluster architecture to expand and to modulate the defence repertoire. Analysis of pathogen recognition gene clusters associated with documented resistance function allowed the identification of adaptive divergence events and the reconstruction of the evolution history of these loci. Differences in candidate pathogen recognition gene number and organization were found between tomato and potato. Most candidate pathogen recognition gene orthologues were distributed at less than perfectly matching positions, suggesting an ongoing lineage-specific rearrangement. Indeed, a local expansion of Toll/Interleukin-1 receptor (TIR)-NBS-leucine-rich repeat (LRR) (TNL) genes in the potato genome was evident. Taken together, these findings have implications for improved understanding of the mechanisms of molecular adaptive selection at Solanum R loci.

Differential control of ethylene responses by GREEN-RIPE and GREEN-RIPE LIKE1 provides evidence for distinct ethylene signaling modules in tomato

The factors that mediate specific responses to the plant hormone ethylene are not fully defined. In particular, it is not known how signaling at the receptor complex can control distinct subsets of ethylene responses. Mutations at the Green-ripe (Gr) and reversion to ethylene sensitivity1 (rte1) loci, which encode homologous proteins of unknown function, influence ethylene responses in tomato (Solanum lycopersicum) and Arabidopsis (Arabidopsis thaliana), respectively. In Arabidopsis, AtRTE1 is required for function of the ETR1 ethylene receptor and acts predominantly through this receptor via direct protein-protein interaction. While most eudicot families including the Brassicaceae possess a single gene that is closely related to AtRTE1, we report that members of the Solanaceae family contain two phylogenetically distinct genes defined by GR and GREEN-RIPE LIKE1 (GRL1), creating the possibility of subfunctionalization. We also show that SlGR and SlGRL1 are differentially expressed in tomato tissues and encode proteins predominantly localized to the Golgi. A combination of overexpression in tomato and complementation of the rte1-3 mutant allele indicates that SlGR and SlGRL1 influence distinct but overlapping ethylene responses. Overexpression of SlGRL1 in the Gr mutant background provides evidence for the existence of different ethylene signaling modules in tomato that are influenced by GR, GRL1, or both. In addition, overexpression of AtRTE1 in tomato leads to reduced ethylene responsiveness in a subset of tissues but does not mimic the Gr mutant phenotype. Together, these data reveal species-specific heterogeneity in the control of ethylene responses mediated by members of the GR/RTE1 family.

Role of an esterase in flavor volatile variation within the tomato clade

Tomato flavor is dependent upon a complex mixture of volatiles including multiple acetate esters. Red-fruited species of the tomato clade accumulate a relatively low content of acetate esters in comparison with the green-fruited species. We show that the difference in volatile ester content between the red- and green-fruited species is associated with insertion of a retrotransposon adjacent to the most enzymatically active member of a family of esterases. This insertion causes higher expression of the esterase, resulting in the reduced levels of multiple esters that are negatively correlated with human preferences for tomato. The insertion was evolutionarily fixed in the red-fruited species, suggesting that high expression of the esterase and consequent low ester content may provide an adaptive advantage in the ancestor of the red-fruited species. These results illustrate at a molecular level how closely related species exhibit major differences in volatile production by altering a volatile-associated catabolic activity.

A genome-wide analysis of the RNA helicase gene family in Solanum lycopersicum

Helicases belong to a class of molecular motor proteins that are found in yeast, animals, and plants. The helicase family is divided into three subfamilies, including the DEAD-box, DEAH-box and DExD/H-box helicases, which are classified based on variations within a common motif, known as motif II. The RNA helicases are involved in every step of RNA metabolism, including nuclear transcription, pre-mRNA splicing, ribosome biogenesis, nucleocytoplasmic transport, translation, RNA decay, and organellar gene expression. The RNA helicase protein family plays a crucial role in plant growth and development as well as in response to biotic and abiotic stresses. However, unlike Arabidopsis, no detailed information regarding the RNA helicase family is currently available for tomato (Solanum lycopersicum) due to a limited number of whole-genome sequences. In this study, we identified a total of 157 RNA helicase genes in the tomato genome. According to the structural features of the motif II region, we classified the tomato RNA helicase genes into DEAD-box, DEAH-box and DExD/H-box helicase genes. But there are 27 RNA helicases not belonging to this three subfamilies, we called that "other helicase". We mapped the 157 RNA helicase genes onto the tomato chromosomes, which range from chr01 to chr12. Microarray and expressed sequence tag data showed that many of these RNA helicase proteins may be involved in diverse biological processes and responses to various stresses. To our knowledge, this is the first report of a genome-wide analysis of the tomato RNA helicase gene family. This study provides valuable information for understanding the classification and putative functions of the RNA helicase gene family in Solanaceae.

Fractionation, stability, and isolate-specificity of QTL for resistance to Phytophthora infestans in cultivated tomato (Solanum lycopersicum)

Cultivated tomato (Solanum lycopersicum) is susceptible to late blight, a major disease caused by Phytophthora infestans, but quantitative resistance exists in the wild tomato species S. habrochaites. Previously, we mapped several quantitative trait loci (QTL) from S. habrochaites and then introgressed each individually into S. lycopersicum. Near-isogenic lines (NILs) were developed, each containing a single introgressed QTL on chromosome 5 or 11. NILs were used to create two recombinant sub-NIL populations, one for each target chromosome region, for higher-resolution mapping. The sub-NIL populations were evaluated for foliar and stem resistance to P. infestans in replicated field experiments over two years, and in replicated growth chamber experiments for resistance to three California isolates. Each of the original single QTL on chromosomes 5 and 11 fractionated into between two and six QTL for both foliar and stem resistance, indicating a complex genetic architecture. The majority of QTL from the field experiments were detected in multiple locations or years, and two of the seven QTL detected in growth chambers were co-located with QTL detected in field experiments, indicating stability of some QTL across environments. QTL that confer foliar and stem resistance frequently co-localized, suggesting that pleiotropy and/or tightly linked genes control the trait phenotypes. Other QTL exhibited isolate-specificity and QTL × environment interactions. Map-based comparisons between QTL mapped in this study and Solanaceae resistance genes/QTL detected in other published studies revealed multiple cases of co-location, suggesting conservation of gene function.