The Sol Genomics Network (solgenomics.net): growing tomatoes using Perl

Genetic and functional traits limit the success of colonisation by arbuscular mycorrhizal fungi in a tomato wild relative

To understand whether domestication had an impact on susceptibility and responsiveness to arbuscular mycorrhizal fungi (AMF) in tomato (Solanum lycopersicum), we investigated two tomato cultivars ("M82" and "Moneymaker") and a panel of wild relatives including S. neorickii, S. habrochaites and S. pennellii encompassing the whole Lycopersicon clade. Most genotypes revealed good AM colonisation levels when inoculated with the AMF Funneliformis mosseae. By contrast, both S. pennellii accessions analysed showed a very low colonisation, but with normal arbuscule morphology, and a negative response in terms of root and shoot biomass. This behaviour was independent of fungal identity and environmental conditions. Genomic and transcriptomic analyses revealed in S. pennellii the lack of genes identified within QTLs for AM colonisation, a limited transcriptional reprogramming upon mycorrhization and a differential regulation of strigolactones and AM-related genes compared to tomato. Donor plants experiments indicated that the AMF could represent a cost for S. pennellii: F. mosseae could extensively colonise the root only when it was part of a mycorrhizal network, but a higher mycorrhization led to a higher inhibition of plant growth. These results suggest that genetics and functional traits of S. pennellii are responsible for the limited extent of AMF colonisation.

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

Pathogen elicitor peptide (pep), systemin, and their receptors in tomato: sequence analysis sheds light on standing disagreements about biotic stress signaling components

Peps are endogenous damage-associated polypeptides that evoke defense responses in plants. Like other damage-associated molecular patterns, Pep signals are transduced by receptors. PEPRs are the receptors that transduce Pep danger signals. This paper identifies new putative Peps in the Solanaceae (including Solanum spp., Nicotiana spp., and Petunia spp.) and Coffea and explores their properties. Using these newly identified Peps we derive sequence logos that present a refinement of the current understanding of the importance of specific residues in the Pep signaling molecules in Solanaceae, including several arginines, prolines that restrict peptide's conformations, and C-terminal asparagine. We examine the degree of disorder in Pep, which is likely important to the mechanism of Pep perception. This work also calls into question some of the evolutionary relationships between Peps in Solanaceae and specific Arabidopsis Peps published in previous literature, culminating in a conclusion that SlPep should not be named SlPep6 due to the lack of conservation of protein sequences in AtPROPEP6 and SlPROPEP, and that SlPep probably does not have two receptors in tomato, based on phylogenetic analysis. Our analyses advance understanding of the Pep signaling system in Solanaceae.

A novel tomato interspecific (Solanum lycopersicum var. cerasiforme and Solanum pimpinellifolium) MAGIC population facilitates trait association and candidate gene discovery in untapped exotic germplasm

We developed a novel eight-way tomato multiparental advanced generation intercross (MAGIC) population to improve the accessibility of tomato relatives genetic resources to geneticists and breeders. The interspecific tomato MAGIC population (ToMAGIC) was obtained by intercrossing four accessions each of Solanum lycopersicum var. cerasiforme and Solanum pimpinellifolium, which are the weedy relative and the ancestor of cultivated tomato, respectively. The eight exotic ToMAGIC founders were selected based on a representation of the genetic diversity and geographical distribution of the two taxa. The resulting MAGIC population comprises 354 lines, which were genotyped using a new 12k tomato single primer enrichment technology panel and yielded 6488 high-quality single-nucleotide polymorphism (SNPs). The genotyping data revealed a high degree of homozygosity, an absence of genetic structure, and a balanced representation of the founder genomes. To evaluate the potential of the ToMAGIC population, a proof of concept was conducted by phenotyping it for fruit size, plant pigmentation, leaf morphology, and earliness. Genome-wide association studies identified strong associations for the studied traits, pinpointing both previously identified and novel candidate genes near or within the linkage disequilibrium blocks. Domesticated alleles for fruit size were recessive and were found, at low frequencies, in wild/ancestral populations. Our findings demonstrate that the newly developed ToMAGIC population is a valuable resource for genetic research in tomato, offering significant potential for identifying new genes that govern key traits in tomato. ToMAGIC lines displaying a pyramiding of traits of interest could have direct applicability for integration into breeding pipelines providing untapped variation for tomato breeding.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.)

BACKGROUND

N6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy.

RESULTS

In this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance.

CONCLUSIONS

This study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.

Multi-year field trials provide a massive repository of trait data on a highly diverse population of tomato and uncover novel determinants of tomato productivity

Tomato (Solanum lycopersicum) is a prominent fruit with rich genetic resources for crop improvement. By using a phenotype-guided screen of over 7900 tomato accessions from around the world, we identified new associations for complex traits such as fruit weight and total soluble solids (Brix). Here, we present the phenotypic data from several years of trials. To illustrate the power of this dataset we use two case studies. First, evaluation of color revealed allelic variation in phytoene synthase 1 that resulted in differently colored or even bicolored fruit. Secondly, in view of the negative relationship between fruit weight and Brix, we pre-selected a subset of the collection that includes high and low Brix values in each category of fruit size. Genome-wide association analysis allowed us to detect novel loci associated with total soluble solid content and fruit weight. In addition, we developed eight F2 biparental intraspecific populations. Furthermore, by taking a phenotype-guided approach we were able to isolate individuals with high Brix values that were not compromised in terms of yield. In addition, the demonstration of novel results despite the high number of previous genome-wide association studies of these traits in tomato suggests that adoption of a phenotype-guided pre-selection of germplasm may represent a useful strategy for finding target genes for breeding.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

Allele mining through TILLING and EcoTILLING approaches in vegetable crops

MAIN CONCLUSION

The present review illustrates a comprehensive overview of the allele mining for genetic improvement in vegetable crops, and allele exploration methods and their utilization in various applications related to pre-breeding of economically important traits in vegetable crops. Vegetable crops have numerous wild descendants, ancestors and terrestrial races that could be exploited to develop high-yielding and climate-resilient varieties resistant/tolerant to biotic and abiotic stresses. To further boost the genetic potential of economic traits, the available genomic tools must be targeted and re-opened for exploitation of novel alleles from genetic stocks by the discovery of beneficial alleles from wild relatives and their introgression to cultivated types. This capability would be useful for giving plant breeders direct access to critical alleles that confer higher production, improve bioactive compounds, increase water and nutrient productivity as well as biotic and abiotic stress resilience. Allele mining is a new sophisticated technique for dissecting naturally occurring allelic variants in candidate genes that influence important traits which could be used for genetic improvement of vegetable crops. Target-induced local lesions in genomes (TILLINGs) is a sensitive mutation detection avenue in functional genomics, particularly wherein genome sequence information is limited or not available. Population exposure to chemical mutagens and the absence of selectivity lead to TILLING and EcoTILLING. EcoTILLING may lead to natural induction of SNPs and InDels. It is anticipated that as TILLING is used for vegetable crops improvement in the near future, indirect benefits will become apparent. Therefore, in this review we have highlighted the up-to-date information on allele mining for genetic enhancement in vegetable crops and methods of allele exploration and their use in pre-breeding for improvement of economic traits.

Solanum tuberosum Microtuber Development under Darkness Unveiled through RNAseq Transcriptomic Analysis

Potato microtuber (MT) development through in vitro techniques are ideal propagules for producing high quality potato plants. MT formation is influenced by several factors, i.e., photoperiod, sucrose, hormones, and osmotic stress. We have previously developed a protocol of MT induction in medium with sucrose (8% w/v), gelrite (6g/L), and 2iP as cytokinin under darkness. To understand the molecular mechanisms involved, we performed a transcriptome-wide analysis. Here we show that 1715 up- and 1624 down-regulated genes were involved in this biological process. Through the protein-protein interaction (PPI) network analyses performed in the STRING database (v11.5), we found 299 genes tightly associated in 14 clusters. Two major clusters of up-regulated proteins fundamental for life growth and development were found: 29 ribosomal proteins (RPs) interacting with 6 PEBP family members and 117 cell cycle (CC) proteins. The PPI network of up-regulated transcription factors (TFs) revealed that at least six TFs-MYB43, TSF, bZIP27, bZIP43, HAT4 and WOX9-may be involved during MTs development. The PPI network of down-regulated genes revealed a cluster of 83 proteins involved in light and photosynthesis, 110 in response to hormone, 74 in hormone mediate signaling pathway and 22 related to aging.

Breedbase: a digital ecosystem for modern plant breeding

Modern breeding methods integrate next-generation sequencing and phenomics to identify plants with the best characteristics and greatest genetic merit for use as parents in subsequent breeding cycles to ultimately create improved cultivars able to sustain high adoption rates by farmers. This data-driven approach hinges on strong foundations in data management, quality control, and analytics. Of crucial importance is a central database able to (1) track breeding materials, (2) store experimental evaluations, (3) record phenotypic measurements using consistent ontologies, (4) store genotypic information, and (5) implement algorithms for analysis, prediction, and selection decisions. Because of the complexity of the breeding process, breeding databases also tend to be complex, difficult, and expensive to implement and maintain. Here, we present a breeding database system, Breedbase (https://breedbase.org/, last accessed 4/18/2022). Originally initiated as Cassavabase (https://cassavabase.org/, last accessed 4/18/2022) with the NextGen Cassava project (https://www.nextgencassava.org/, last accessed 4/18/2022), and later developed into a crop-agnostic system, it is presently used by dozens of different crops and projects. The system is web based and is available as open source software. It is available on GitHub (https://github.com/solgenomics/, last accessed 4/18/2022) and packaged in a Docker image for deployment (https://hub.docker.com/u/breedbase, last accessed 4/18/2022). The Breedbase system enables breeding programs to better manage and leverage their data for decision making within a fully integrated digital ecosystem.

Comprehensive Analysis of N6-Methyladenosine Regulatory Genes from Citrus grandis and Expression Profilings in the Fruits of “Huajuhong” (C. grandis “Tomentosa”) during Various Development Stages

Citrus grandis “Tomentosa” (“Huajuhong”) is a famous traditional Chinese medicine. The aim of the present study is to provide a comprehensive characterization of the m6A regulatory genes from C. grandis, and examine their expression patterns in fruits of C. grandis “Tomentosa” during various developmental stages. A total of 26 N6-methyladenosine (m6A) regulatory proteins were identified from the genome of C. grandis, which were distributed across nine chromosomes in C. grandis. Phylogenetic relationships revealed that all m6A regulatory genes were divided into groups of m6A writers, erasers, and readers. The m6A writer groups included CgMTA, CgMTB, and CgMTC three MTs (methyltransferases), one CgVIR (virilizer), one CgHAKAI (E3 ubiquitin ligase HAKAI), and one CgFIP37 (FKBP interacting protein 37). Moreover, 10 CgALKBH (α-ketoglutarate-dependent dioxygenase homolog) members (numbered from CgALKBH1 to CgALKBH10) and 10 CgECT (C-terminal region) members (numbered from CgECT1 to CgECT10) in C. grandis were identified as m6A erasers and readers, respectively. The domain structures and motif architectures among the groups of m6A writers, erasers, and readers were diverse. Cis-acting elements in the promoters of the 26 m6A regulatory genes predicted that the abscisic acid-responsive (ABA) element (ABRE) was present on the promoters of 19 genes. In addition, the expression profiles of all m6A regulatory genes were examined in the fruits of two varieties of C. grandis “Tomentosa” during different growth stages to give basic hints for further investigation of the function of the N6-methyladenosine regulatory genes in C. grandis “Tomentosa”.

Comparative transcriptomic analysis of apple and peach fruits: insights into fruit type specification

Fruits represent key evolutionary innovations in angiosperms and exhibit diverse types adapted for seed dissemination. However, the mechanisms that underlie fruit type diversity are not understood. The Rosaceae family comprises many different fruit types, including 'pome' and 'drupe' fruits, and hence is an excellent family for investigating the genetic basis of fruit type specification. Using comparative transcriptomics, we investigated the molecular events that correlate with pome (apple) and drupe (peach) fleshy fruit development, focusing on the earliest stages of fruit initiation. We identified PI and TM6, MADS box genes whose expression negatively correlates with fruit flesh-forming tissues irrespective of fruit type. In addition, the MADS box gene FBP9 is expressed in fruit-forming tissues in both species, and was lost multiple times in the genomes of dry-fruit-forming eudicots including Arabidopsis. Network analysis reveals co-expression between FBP9 and photosynthesis genes in both apple and peach, suggesting that FBP9 and photosynthesis may both promote fleshy fruit development. The large transcriptomic datasets at the earliest stages of pome and drupe fruit development provide rich resources for comparative studies, and the work provides important insights into fruit-type specification.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium species.

DeepLRR: An Online Webserver for Leucine-Rich-Repeat Containing Protein Characterization Based on Deep Learning

Members of the leucine-rich repeat (LRR) superfamily play critical roles in multiple biological processes. As the LRR unit sequence is highly variable, accurately predicting the number and location of LRR units in proteins is a highly challenging task in the field of bioinformatics. Existing methods still need to be improved, especially when it comes to similarity-based methods. We introduce our DeepLRR method based on a convolutional neural network (CNN) model and LRR features to predict the number and location of LRR units in proteins. We compared DeepLRR with six existing methods using a dataset containing 572 LRR proteins and it outperformed all of them when it comes to overall F1 score. In addition, DeepLRR has integrated identifying plant disease-resistance proteins (NLR, LRR-RLK, LRR-RLP) and non-canonical domains. With DeepLRR, 223, 191 and 183 LRR-RLK genes in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa ssp. Japonica) and tomato (Solanum lycopersicum) genomes were re-annotated, respectively. Chromosome mapping and gene cluster analysis revealed that 24.2% (54/223), 29.8% (57/191) and 16.9% (31/183) of LRR-RLK genes formed gene cluster structures in Arabidopsis, rice and tomato, respectively. Finally, we explored the evolutionary relationship and domain composition of LRR-RLK genes in each plant and distributions of known receptor and co-receptor pairs. This provides a new perspective for the identification of potential receptors and co-receptors.

Gene Expression Analysis of Microtubers of Potato Solanum tuberosum L. Induced in Cytokinin Containing Medium and Osmotic Stress

Potato microtuber productions through in vitro techniques are ideal propagules for producing high quality seed potatoes. Microtuber development is influenced by several factors, i.e., high content sucrose and cytokinins are among them. To understand a molecular mechanism of microtuberization using osmotic stress and cytokinin signaling will help us to elucidate this process. We demonstrate in this work a rapid and efficient protocol for microtuber development and gene expression analysis. Medium with high content of sucrose and gelrite supplemented with 2iP as cytokinin under darkness condition produced the higher quantity and quality of microtubers. Gene expression analysis of genes involved in the two-component signaling system (StHK1), cytokinin signaling, (StHK3, StHP4, StRR1) homeodomains (WUSCHEL, POTH1, BEL5), auxin signaling, ARF5, carbon metabolism (TPI, TIM), protein synthesis, NAC5 and a morphogenetic regulator of tuberization (POTH15) was performed by qPCR real time. Differential gene expression was observed during microtuber development. Gene regulation of two component and cytokinin signaling is taking place during this developmental process, yielding more microtubers. Further analysis of each component is required to elucidate it.

Genome-Wide Investigation of N6-Methyladenosine Regulatory Genes and Their Roles in Tea (Camellia sinensis) Leaves During Withering Process

N6-methyladenosine (m6A), one of the internal modifications of RNA molecules, can directly influence RNA abundance and function without altering the nucleotide sequence, and plays a pivotal role in response to diverse environmental stresses. The precise m6A regulatory mechanism comprises three types of components, namely, m6A writers, erasers, and readers. To date, the research focusing on m6A regulatory genes in plant kingdom is still in its infancy. Here, a total of 34 m6A regulatory genes were identified from the chromosome-scale genome of tea plants. The expansion of m6A regulatory genes was driven mainly by whole-genome duplication (WGD) and segmental duplication, and the duplicated gene pairs evolved through purifying selection. Gene structure analysis revealed that the sequence variation contributed to the functional diversification of m6A regulatory genes. Expression pattern analysis showed that most m6A regulatory genes were differentially expressed under environmental stresses and tea-withering stage. These observations indicated that m6A regulatory genes play essential roles in response to environmental stresses and tea-withering stage. We also found that RNA methylation and DNA methylation formed a negative feedback by interacting with each other's methylation regulatory genes. This study provided a foundation for understanding the m6A-mediated regulatory mechanism in tea plants under environmental stresses and tea-withering stage.

High density genotype storage for plant breeding in the Chado schema of Breedbase

Modern breeding programs routinely use genome-wide information for selecting individuals to advance. The large volumes of genotypic information required present a challenge for data storage and query efficiency. Major use cases require genotyping data to be linked with trait phenotyping data. In contrast to phenotyping data that are often stored in relational database schemas, next-generation genotyping data are traditionally stored in non-relational storage systems due to their extremely large scope. This study presents a novel data model implemented in Breedbase (https://breedbase.org/) for uniting relational phenotyping data and non-relational genotyping data within the open-source PostgreSQL database engine. Breedbase is an open-source, web-database designed to manage all of a breeder's informatics needs: management of field experiments, phenotypic and genotypic data collection and storage, and statistical analyses. The genotyping data is stored in a PostgreSQL data-type known as binary JavaScript Object Notation (JSONb), where the JSON structures closely follow the Variant Call Format (VCF) data model. The Breedbase genotyping data model can handle different ploidy levels, structural variants, and any genotype encoded in VCF. JSONb is both compressed and indexed, resulting in a space and time efficient system. Furthermore, file caching maximizes data retrieval performance. Integration of all breeding data within the Chado database schema retains referential integrity that may be lost when genotyping and phenotyping data are stored in separate systems. Benchmarking demonstrates that the system is fast enough for computation of a genomic relationship matrix (GRM) and genome wide association study (GWAS) for datasets involving 1,325 diploid Zea mays, 314 triploid Musa acuminata, and 924 diploid Manihot esculenta samples genotyped with 955,690, 142,119, and 287,952 genotype-by-sequencing (GBS) markers, respectively.

EXPath 2.0: An Updated Database for Integrating High-Throughput Gene Expression Data with Biological Pathways

Co-expressed genes tend to have regulatory relationships and participate in similar biological processes. Construction of gene correlation networks from microarray or RNA-seq expression data has been widely applied to study transcriptional regulatory mechanisms and metabolic pathways under specific conditions. Furthermore, since transcription factors (TFs) are critical regulators of gene expression, it is worth investigating TFs on the promoters of co-expressed genes. Although co-expressed genes and their related metabolic pathways can be easily identified from previous resources, such as EXPath and EXPath Tool, this information is not simultaneously available to identify their regulatory TFs. EXPath 2.0 is an updated database for the investigation of regulatory mechanisms in various plant metabolic pathways with 1,881 microarray and 978 RNA-seq samples. There are six significant improvements in EXPath 2.0: (i) the number of species has been extended from three to six to include Arabidopsis, rice, maize, Medicago, soybean and tomato; (ii) gene expression at various developmental stages have been added; (iii) construction of correlation networks according to a group of genes is available; (iv) hierarchical figures of the enriched Gene Ontology (GO) terms are accessible; (v) promoter analysis of genes in a metabolic pathway or correlation network is provided; and (vi) user's gene expression data can be uploaded and analyzed. Thus, EXPath 2.0 is an updated platform for investigating gene expression profiles and metabolic pathways under specific conditions. It facilitates users to access the regulatory mechanisms of plant biological processes. The new version is available at http://EXPath.itps.ncku.edu.tw.

The res (restored cell structure by salinity) tomato mutant reveals the role of the DEAD-box RNA helicase SlDEAD39 in plant development and salt response

Increasing evidences highlight the importance of DEAD-box RNA helicases in plant development and stress responses. In a previous study, we characterized the tomato res mutant (restored cell structure by salinity), showing chlorosis and development alterations that reverted under salt-stress conditions. Map-based cloning demonstrates that RES gene encodes SlDEAD39, a chloroplast-targeted DEAD-box RNA helicase. Constitutive expression of SlDEAD39 complements the res mutation, while the silencing lines had a similar phenotype than res mutant, which is also reverted under salinity. Functional analysis of res mutant proved SlDEAD39 is involved in the in vivo processing of the chloroplast, 23S rRNA, at the hidden break-B site, a feature also supported by in vitro binding experiments of the protein. In addition, our results show that other genes coding for chloroplast-targeted DEAD-box proteins are induced by salt-stress, which might explain the rescue of the res mutant phenotype. Interestingly, salinity restored the phenotype of res adult plants by increasing their sugar content and fruit yield. Together, these results propose an unprecedented role of a DEAD-box RNA helicase in regulating plant development and stress response through the proper ribosome and chloroplast functioning, which, in turn, represents a potential target to improve salt tolerance in tomato crops.

Proteomic analysis and interactions network in leaves of mycorrhizal and nonmycorrhizal sorghum plants under water deficit

For understanding the water deficit stress mechanism in sorghum, we conducted a physiological and proteomic analysis in the leaves of Sorghum bicolor L. Moench (a drought tolerant crop model) of non-colonized and colonized plants with a consortium of arbuscular mycorrhizal fungi. Physiological results indicate that mycorrhizal fungi association enhances growth and photosynthesis in plants, under normal and water deficit conditions. 2D-electrophoresis profiles revealed 51 differentially accumulated proteins in response to water deficit, of which HPLC/MS successfully identified 49. Bioinformatics analysis of protein–protein interactions revealed the participation of different metabolic pathways in nonmycorrhizal compared to mycorrhizal sorghum plants under water deficit. In noninoculated plants, the altered proteins are related to protein synthesis and folding (50S ribosomal protein L1, 30S ribosomal protein S10, Nascent polypeptide-associated complex subunit alpha), coupled with multiple signal transduction pathways, guanine nucleotide-binding beta subunit (Rack1) and peptidyl-prolyl-cis-trans isomerase (ROC4). In contrast, in mycorrhizal plants, proteins related to energy metabolism (ATP synthase-24kDa, ATP synthase β), carbon metabolism (malate dehydrogenase, triosephosphate isomerase, sucrose-phosphatase), oxidative phosphorylation (mitochondrial-processing peptidase) and sulfur metabolism (thiosulfate/3-mercaptopyruvate sulfurtransferase) were found. Our results provide a set of proteins of different metabolic pathways involved in water deficit produced by sorghum plants alone or associated with a consortium of arbuscular mycorrhizal fungi isolated from the tropical rain forest Los Tuxtlas Veracruz, México.

Molecular breeding approaches for production of disease-resilient commercially important tobacco

Tobacco is one of the most widely cultivated nonfood cash crops, a source of income, model organism for plant molecular research, a natural pesticide and of pharmaceutical importance. First domesticated in South Americas, the modern-day tobacco (Nicotiana tabacum) is now cultivated in more than 125 countries to generate revenues worth billions of dollars each year. However, the production of this crop is highly threatened by the global presence of devastating infectious agents, which cause huge fiscal loss. These threats have been battled through breeding for acquiring disease resilience in tobacco plants, first, via conventional and now with the use of modern molecular breeding approaches. For efficacy and precision, the characterization of the genetic components underlying disease resistance is the key tool in tobacco for resistance breeding programs. The past few decades have witnessed significant progress in resilience breeding through advanced molecular techniques. The current review discusses history of tobacco breeding since its time of origin till date, highlighting the most widely used techniques and recent advances in molecular research and strategies for resistance breeding. In addition, we narrate the budding possibilities for the future. This review will provide a comprehensive and valuable information for the tobacco growers and researchers to deal with the destructive infectious diseases.

De novo whole-genome assembly in Chrysanthemum seticuspe, a model species of Chrysanthemums, and its application to genetic and gene discovery analysis

Cultivated chrysanthemum (Chrysanthemum morifolium Ramat.) is one of the most economically important ornamental crops grown worldwide. It has a complex hexaploid genome (2n = 6x = 54) and large genome size. The diploid Chrysanthemum seticuspe is often used as a model of cultivated chrysanthemum, since the two species are closely related. To expand our knowledge of the cultivated chrysanthemum, we here performed de novo whole-genome assembly in C. seticuspe using the Illumina sequencing platform. XMRS10, a C. seticuspe accession developed by five generations of self-crossing from a self-compatible strain, AEV2, was used for genome sequencing. The 2.72 Gb of assembled sequences (CSE_r1.0), consisting of 354,212 scaffolds, covered 89.0% of the 3.06 Gb C. seticuspe genome estimated by k-mer analysis. The N50 length of scaffolds was 44,741 bp. For protein-encoding genes, 71,057 annotated genes were deduced (CSE_r1.1_cds). Next, based on the assembled genome sequences, we performed linkage map construction, gene discovery and comparative analyses for C. seticuspe and cultivated chrysanthemum. The generated C. seticuspe linkage map revealed skewed regions in segregation on the AEV2 genome. In gene discovery analysis, candidate flowering-related genes were newly found in CSE_r1.1_cds. Moreover, single nucleotide polymorphism identification and annotation on the C. × morifolium genome showed that the C. seticuspe genome was applicable to genetic analysis in cultivated chrysanthemums. The genome sequences assembled herein are expected to contribute to future chrysanthemum studies. In addition, our approach demonstrated the usefulness of short-read genome assembly and the importance of choosing an appropriate next genome sequencing technology based on the purpose of the post-genome analysis.

Agrobacterium tumefaciens-Mediated Transformation of Tomato

Tomato is both an important food crop and serves as a model plant species that is used for various research investigations including understanding gene function. Transformation is commonly utilized to facilitate these investigations in combination with all the extensive genetic and genomic resources available for tomato. The transformation protocol routinely used in our laboratory has been applied to many different tomato genotypes and relies on Agrobacterium tumefaciens infection of young cotyledon sections. We have used vector systems for overexpression, RNA interference for gene silencing, and CRISPR/Cas9 for genome editing. Vectors used to design gene constructs contained selectable marker genes that conferred resistance to kanamycin, hygromycin, and the herbicide component, bialaphos. The protocol we follow for Agrobacterium-mediated transformation of both cultivated and wild species of tomato is detailed in this chapter.

Regulation of Root Angle and Gravitropism

Regulation of plant root angle is critical for obtaining nutrients and water and is an important trait for plant breeding. A plant’s final, long-term root angle is the net result of a complex series of decisions made by a root tip in response to changes in nutrient availability, impediments, the gravity vector and other stimuli. When a root tip is displaced from the gravity vector, the short-term process of gravitropism results in rapid reorientation of the root toward the vertical. Here, we explore both short- and long-term regulation of root growth angle, using natural variation in tomato to identify shared and separate genetic features of the two responses. Mapping of expression quantitative trait loci mapping and leveraging natural variation between and within species including Arabidopsis suggest a role for PURPLE ACID PHOSPHATASE 27 and CELL DIVISION CYCLE 73 in determining root angle.

TM8 represses developmental timing in Nicotiana benthamiana and has functionally diversified in angiosperms

BACKGROUND

MADS-box genes are key regulators of plant reproductive development and members of most lineages of this gene family have been extensively studied. However, the function and diversification of the ancient TM8 lineage remains elusive to date. The available data suggest a possible function in flower development in tomato and fast evolution through numerous gene loss events in flowering plants.

RESULTS

We show the broad conservation of TM8 within angiosperms and find that in contrast to other MADS-box gene lineages, no gene duplicates have been retained after major whole genome duplication events. Through knock-down of NbTM8 by virus induced gene silencing in Nicotiana benthamiana, we show that NbTM8 represses miR172 together with another MADS-box gene, SHORT VEGETATIVE PHASE (NbSVP). In the closely related species Petunia hybrida, PhTM8 is not expressed under the conditions we investigated and consistent with this, a knock-out mutant did not show a phenotype. Finally, we generated transgenic tomato plants in which TM8 was silenced or ectopically expressed, but these plants did not display a clear phenotype. Therefore, no clear function could be confirmed for Solanum lycopersium.

CONCLUSIONS

While the presence of TM8 is generally conserved, it remains difficult to propose a general function in angiosperms. Based on all the available data to date, supplemented with our own results, TM8 function seems to have diversified quickly throughout angiosperms and acts as repressor of miR172 in Nicotiana benthamiana, together with NbSVP.

The coupling of transcriptome and proteome adaptation during development and heat stress response of tomato pollen

BACKGROUND

Pollen development is central for plant reproduction and is assisted by changes of the transcriptome and proteome. At the same time, pollen development and viability is largely sensitive to stress, particularly to elevated temperatures. The transcriptomic and proteomic changes during pollen development and of different stages in response to elevated temperature was targeted to define the underlying molecular principles.

RESULTS

The analysis of the transcriptome and proteome of Solanum lycopersicum pollen at tetrad, post-meiotic and mature stage before and after heat stress yielded a decline of the transcriptome but an increase of the proteome size throughout pollen development. Comparison of the transcriptome and proteome led to the discovery of two modes defined as direct and delayed translation. Here, genes of distinct functional processes are under the control of direct and delayed translation. The response of pollen to elevated temperature occurs rather at proteome, but not as drastic at the transcriptome level. Heat shock proteins, proteasome subunits, ribosomal proteins and eukaryotic initiation factors are most affected. On the example of heat shock proteins we demonstrate a decoupling of transcript and protein levels as well as a distinct regulation between the developmental stages.

CONCLUSIONS

The transcriptome and proteome of developing pollen undergo drastic changes in composition and quantity. Changes at the proteome level are a result of two modes assigned as direct and delayed translation. The response of pollen to elevated temperature is mainly regulated at the proteome level, whereby proteins related to synthesis and degradation of proteins are most responsive and might play a central role in the heat stress response of pollen.

Alternative splicing in tomato pollen in response to heat stress

Alternative splicing (AS) is a key control mechanism influencing signal response cascades in different developmental stages and under stress conditions. In this study, we examined heat stress (HS)-induced AS in the heat sensitive pollen tissue of two tomato cultivars. To obtain the entire spectrum of HS-related AS, samples taken directly after HS and after recovery were combined and analysed by RNA-seq. For nearly 9,200 genes per cultivar, we observed at least one AS event under HS. In comparison to control, for one cultivar we observed 76% more genes with intron retention (IR) or exon skipping (ES) under HS. Furthermore, 2,343 genes had at least one transcript with IR or ES accumulated under HS in both cultivars. These genes are involved in biological processes like protein folding, gene expression and heat response. Transcriptome assembly of these genes revealed that most of the alternative spliced transcripts possess truncated coding sequences resulting in partial or total loss of functional domains. Moreover, 141 HS specific and 22 HS repressed transcripts were identified. Further on, we propose AS as layer of stress response regulating constitutively expressed genes under HS by isoform abundance.

Evolution of Eukaryal and Archaeal Pseudouridine Synthase Pus10

In archaea, pseudouridine (Ψ) synthase Pus10 modifies uridine (U) to Ψ at positions 54 and 55 of tRNA. In contrast, Pus10 is not found in bacteria, where modifications at those two positions are carried out by TrmA (U54 to m⁵U54) and TruB (U55 to Ψ55). Many eukaryotes have an apparent redundancy; their genomes contain orthologs of archaeal Pus10 and bacterial TrmA and TruB. Although eukaryal Pus10 genes share a conserved catalytic domain with archaeal Pus10 genes, their biological roles are not clear for the two reasons. First, experimental evidence suggests that human Pus10 participates in apoptosis induced by the tumor necrosis factor-related apoptosis-inducing ligand. Whether the function of human Pus10 is in place or in addition to of Ψ synthesis in tRNA is unknown. Second, Pus10 is found in earlier evolutionary branches of fungi (such as chytrid Batrachochytrium) but is absent in all dikaryon fungi surveyed (Ascomycetes and Basidiomycetes). We did a comprehensive analysis of sequenced genomes and found that orthologs of Pus10, TrmA, and TruB were present in all the animals, plants, and protozoa surveyed. This indicates that the common eukaryotic ancestor possesses all the three genes. Next, we examined 116 archaeal and eukaryotic Pus10 protein sequences to find that Pus10 existed as a single copy gene in all the surveyed genomes despite ancestral whole genome duplications had occurred. This indicates a possible deleterious gene dosage effect. Our results suggest that functional redundancy result in gene loss or neofunctionalization in different evolutionary lineages.

SolCyc: a database hub at the Sol Genomics Network (SGN) for the manual curation of metabolic networks in Solanum and Nicotiana specific databases

Database URL

https://solgenomics.net/tools/solcyc/.

NbEXPA1, an α-expansin, is plasmodesmata-specific and a novel host factor for potyviral infection

Plasmodesmata (PD), unique to the plant kingdom, are structurally complex microchannels that cross the cell wall to establish symplastic communication between neighbouring cells. Viral intercellular movement occurs through PD. To better understand the involvement of PD in viral infection, we conducted a quantitative proteomic study on the PD-enriched fraction from Nicotiana benthamiana leaves in response to infection by Turnip mosaic virus (TuMV). We report the identification of a total of 1070 PD protein candidates, of which 100 (≥2-fold increase) and 48 (≥2-fold reduction) are significantly differentially accumulated in the PD-enriched fraction, when compared with protein levels in the corresponding healthy control. Among the differentially accumulated PD protein candidates, we show that an α-expansin designated NbEXPA1, a cell wall loosening protein, is PD-specific. TuMV infection downregulates NbEXPA1 mRNA expression and protein accumulation. We further demonstrate that NbEXPA1 is recruited to the viral replication complex via the interaction with NIb, the only RNA-dependent RNA polymerase of TuMV. Silencing of NbEXPA1 inhibits plant growth and TuMV infection, whereas overexpression of NbEXPA1 promotes viral replication and intercellular movement. These data suggest that NbEXPA1 is a host factor for potyviral infection. This study not only generates a PD-proteome dataset that is useful in future studies to expound PD biology and PD-mediated virus-host interactions but also characterizes NbEXPA1 as the first PD-specific cell wall loosening protein and its essential role in potyviral infection.

Gene editing in tomatoes

Tomato is an effective model plant species because it possesses the qualities necessary for genetic and functional studies, but is also a food crop making what is learned more translatable for crop improvement when compared with other non-food crop models. The availability of genome sequences for many genotypes and amenability to transformation methodologies (Agrobacterium-mediated, direct DNA uptake via protoplasts, biolistics) make tomato the perfect platform to study the application of gene-editing technologies. This review includes information related to tomato transformation methodology, one of the necessary requirements for gene editing, along with the status of site-directed mutagenesis by TALENs (transcription activator-like effector nucleases) and CRISPR/Cas (clustered regularly interspaced short palindromic repeats and CRISPR-associated Proteins). In addition to the reports on proof-of-concept experiments to demonstrate the feasibility of gene editing in tomato, there are many reports that show the power of these technologies for modification of traits, such as fruit characteristics (ripening, size, and parthenocarpy), pathogen susceptibility, architecture (plant and inflorescence), and metabolic engineering. Also highlighted in this review are reports on the application of a recent CRISPR technology called base editing that allows the modification of one base pair in a gene sequence and a strategy that takes advantage of a geminivirus replicon for delivery of DNA repair template.

Transcriptome profiling of tobacco (Nicotiana tabacum) pollen and pollen tubes

BACKGROUND

Pollen tube growth is essential for plant reproduction and represents a widely employed model to investigate polarized cell expansion, a process important for plant morphogenesis and development. Cellular and regulatory mechanisms underlying pollen tube elongation are under intense investigation, which stands to greatly benefit from a comprehensive understanding of global gene expression profiles in pollen and pollen tubes. Here, RNA sequencing technology was applied to de novo assemble a Nicotiana tabacum male gametophytic transcriptome and to compare transcriptome profiles at two different stages of gametophyte development: mature pollen grains (MPG) and pollen tubes grown for six hours in vitro (PT6).

RESULTS

De novo assembly of data obtained by 454 sequencing of a normalized cDNA library representing tobacco pollen and pollen tube mRNA (pooled mRNA isolated from mature pollen grains [MPG] and from pollen tubes grown in vitro for 3 [PT3] or 6 [PT6] hours) resulted in the identification of 78,364 unigenes. Among these unigenes, which mapped to 24,933 entries in the Sol Genomics Network (SGN) N. tabacum unigene database, 24,672 were predicted to represent full length cDNAs. In addition, quantitative analyses of data obtained by Illumina sequencing of two separate non-normalized MPG and PT6 cDNA libraries showed that 8979 unigenes were differentially expressed (differentially expressed unigenes: DEGs) between these two developmental stages at a FDR q-value of <0.0001. Interestingly, whereas most of these DEGs were downregulated in PT6, the minor fraction of DEGs upregulated in PT6 was enriched for GO (gene ontology) functions in pollen tube growth or fertilization.

CONCLUSIONS

A major output of our study is the development of two different high-quality databases representing the tobacco male gametophytic transcriptome and containing encompassing information about global changes in gene expression after pollen germination. Quantitative analyses of these databases 1) indicated that roughly 30% of all tobacco genes are expressed in the male gametophyte, and 2) support previous observations suggesting a global reduction of transcription after pollen germination. Interestingly, a small number of genes, many of which predicted to function in pollen tube growth or fertilization, were found to be upregulated in elongating pollen tubes despite globally reduced transcription.

DES-TOMATO: A Knowledge Exploration System Focused On Tomato Species

Tomato is the most economically important horticultural crop used as a model to study plant biology and particularly fruit development. Knowledge obtained from tomato research initiated improvements in tomato and, being transferrable to other such economically important crops, has led to a surge of tomato-related research and published literature. We developed DES-TOMATO knowledgebase (KB) for exploration of information related to tomato. Information exploration is enabled through terms from 26 dictionaries and combination of these terms. To illustrate the utility of DES-TOMATO, we provide several examples how one can efficiently use this KB to retrieve known or potentially novel information. DES-TOMATO is free for academic and nonprofit users and can be accessed at http://cbrc.kaust.edu.sa/des_tomato/, using any of the mainstream web browsers, including Firefox, Safari and Chrome.

In silico analysis reveals widespread presence of three gene families, MAPK, MAPKK and MAPKKK, of the MAPK cascade from crop plants of Solanaceae in comparison to the distantly-related syntenic species from Rubiaceae, coffee

Mitogen-activated protein kinases (MAPKs) are an important family of genes which play roles in vital plant processes, and they also help in coping against various kinds of environmental stresses including abiotic as well as biotic factors. The advancement of genomics calls for the annotation, identification, and detailed processing of the essential gene families in plants in order to provide insights into the importance of their central roles as well as for providing the basis for making their growth vigorous even under stressed conditions and, ultimately, to benefit from them by foreseeing the potential threats to their growth. In the current study, MAPK, MAPKK, and MAPKKK families of the MAPK cascade were identified and reported from five different agriculturally and economically important crop species of the Solanaceae and Rubiaceae families based on conserved signature motifs aligned throughout the members of the families under this gene superfamily. Genes reported from the species after strict filtering were: 89, tomato; 108, potato; 63, eggplant; 79, pepper; 64, coffee. These MAPKs were found to be randomly distributed throughout the genome on the chromosomes of the respective species. Various characteristics of the identified genes were studied including gene structure, gene and coding sequence length, protein length, isoelectric point, molecular weight, and subcellular localization. Moreover, maximum likelihood test of phylogeny was conducted on the retrieved sequences for the three MAPK cascade families to determine their homologous relationships which were also analyzed quantitatively by heat plots.

Genome-Wide Prediction of Metabolic Enzymes, Pathways, and Gene Clusters in Plants

Plant metabolism underpins many traits of ecological and agronomic importance. Plants produce numerous compounds to cope with their environments but the biosynthetic pathways for most of these compounds have not yet been elucidated. To engineer and improve metabolic traits, we need comprehensive and accurate knowledge of the organization and regulation of plant metabolism at the genome scale. Here, we present a computational pipeline to identify metabolic enzymes, pathways, and gene clusters from a sequenced genome. Using this pipeline, we generated metabolic pathway databases for 22 species and identified metabolic gene clusters from 18 species. This unified resource can be used to conduct a wide array of comparative studies of plant metabolism. Using the resource, we discovered a widespread occurrence of metabolic gene clusters in plants: 11,969 clusters from 18 species. The prevalence of metabolic gene clusters offers an intriguing possibility of an untapped source for uncovering new metabolite biosynthesis pathways. For example, more than 1,700 clusters contain enzymes that could generate a specialized metabolite scaffold (signature enzymes) and enzymes that modify the scaffold (tailoring enzymes). In four species with sufficient gene expression data, we identified 43 highly coexpressed clusters that contain signature and tailoring enzymes, of which eight were characterized previously to be functional pathways. Finally, we identified patterns of genome organization that implicate local gene duplication and, to a lesser extent, single gene transposition as having played roles in the evolution of plant metabolic gene clusters.

Label-free deep shotgun proteomics reveals protein dynamics during tomato fruit tissues development

Current innovations in mass-spectrometry-based technologies allow deep coverage of protein expression. Despite its immense value and in contrast to transcriptomics, only a handful of studies in crop plants engaged with global proteome assays. Here, we present large-scale shotgun proteomics profiling of tomato fruit across two key tissues and five developmental stages. A total of 7738 individual protein groups were identified and reliably measured at least in one of the analyzed tissues or stages. The depth of our assay enabled identification of 61 differentially expressed transcription factors, including renowned ripening-related regulators and elements of ethylene signaling. Significantly, we measured proteins involved in 83% of all predicted enzymatic reactions in the tomato metabolic network. Hence, proteins representing almost the complete set of reactions in major metabolic pathways were identified, including the cytosolic and plastidic isoprenoid and the phenylpropanoid pathways. Furthermore, the data allowed us to discern between protein isoforms according to expression patterns, which is most significant in light of the weak transcript-protein expression correspondence. Finally, visualization of changes in protein abundance associated with a particular process provided us with a unique view of skin and flesh tissues in developing fruit. This study adds a new dimension to the existing genomic, transcriptomic and metabolomic resources. It is therefore likely to promote translational and post-translational research in tomato and additional species, which is presently focused on transcription.

Oil palm drought inducible DREB1 induced expression of DRE/CRT- and non-DRE/CRT-containing genes in lowland transgenic tomato under cold and PEG treatments

Dehydration-responsive element binding (DREB) transcription factor plays an important role in controlling the expression of abiotic stress responsive genes. An intronless oil palm EgDREB1 was isolated and confirmed to be a nuclear localized protein. Electrophoretic mobility shift and yeast one-hybrid assays validated its ability to interact with DRE/CRT motif. Its close evolutionary relation to the dicot NtDREB2 suggests a universal regulatory role. In order to determine its involvement in abiotic stress response, functional characterization was performed in oil palm seedlings subjected to different levels of drought severity and in EgDREB1 transgenic tomato seedlings treated by abiotic stresses. Its expression in roots and leaves was compared with several antioxidant genes using quantitative real-time PCR. Early accumulation of EgDREB1 in oil palm roots under mild drought suggests possible involvement in the initiation of signaling communication from root to shoot. Ectopic expression of EgDREB1 in T1 transgenic tomato seedlings enhanced expression of DRE/CRT and non-DRE/CRT containing genes, including tomato peroxidase (LePOD), ascorbate peroxidase (LeAPX), catalase (LeCAT), superoxide dismutase (LeSOD), glutathione reductase (LeGR), glutathione peroxidase (LeGP), heat shock protein 70 (LeHSP70), late embryogenesis abundant (LeLEA), metallothionine type 2 (LeMET2), delta 1-pyrroline-5- carboxylate synthetase (LePCS), ABA-aldehyde oxidase (LeAAO) and 9-cis- Epoxycarotenoid dioxygenase (LeECD) under PEG treatment and cold stress (4 °C). Altogether, these findings suggest that EgDREB1 is a functional regulator in enhancing tolerance to drought and cold stress.

Plant Genome DataBase Japan (PGDBj)

A portal website that integrates a variety of information related to genomes of model and crop plants from databases (DBs) and the literature was generated. This website, named the Plant Genome DataBase Japan (PGDBj, http://pgdbj. jp/en/ ), is comprised of three component DBs and a cross-search engine which provides a seamless search over their contents. One of the three component DBs is the Ortholog DB, which provides gene cluster information based on the amino acid sequence similarity. Over 1,000,000 amino acid sequences of 40 Viridiplantae species were collected from the public DNA DBs, and plant genome DBs such as TAIR and RAP-DB were subjected to reciprocal BLAST searches for clustering. Another component DB is the Plant Resource DB for genomic- and bio-resources. This DB also integrates the SABRE DB, which provides cDNA and genome sequence resources maintained in the RIKEN BioResource Center and National BioResource Projects Japan. The third component DB of PGDBj is the DNA Marker DB, which manually or automatically collects curated information on DNA markers, quantitative trait loci (QTL), and related genetic linkage maps, from the literature and external DBs. By combining these component DBs and a cross-search engine, PGDBj serves as a useful platform to study genetic systems for both fundamental and applied researches for a wide range of plant species.

Agroinfiltration contributes to VP1 recombinant protein degradation

There is a growing interest in applying tobacco agroinfiltration for recombinant protein production in a plant based system. However, in such a system, the action of proteases might compromise recombinant protein production. Protease sensitivity of model recombinant foot-and-mouth disease (FMD) virus P1-polyprotein (P1) and VP1 (viral capsid protein 1) as well as E. coli glutathione reductase (GOR) were investigated. Recombinant VP1 was more severely degraded when treated with the serine protease trypsin than when treated with the cysteine protease papain. Cathepsin L- and B-like as well as legumain proteolytic activities were elevated in agroinfiltrated tobacco tissues and recombinant VP1 was degraded when incubated with such a protease-containing tobacco extract. In silico analysis revealed potential protease cleavage sites within the P1, VP1 and GOR sequences. The interaction modeling of the single VP1 protein with the proteases papain and trypsin showed greater proximity to proteolytic active sites compared to modeling with the entire P1-polyprotein fusion complex. Several plant transcripts with differential expression were detected 24 hr post-agroinfiltration when the RNA-seq technology was applied to identify changed protease transcripts using the recently available tobacco draft genome. Three candidate genes were identified coding for proteases which included the Responsive-to-Desiccation-21 (RD21) gene and genes for coding vacuolar processing enzymes 1a (NbVPE1a) and 1b (NbVPE1b). The data demonstrates that the tested recombinant proteins are sensitive to protease action and agroinfiltration induces the expression of potential proteases that can compromise recombinant protein production.

Indel Group in Genomes (IGG) Molecular Genetic Markers

Genetic markers are essential when developing or working with genetically variable populations. Indel Group in Genomes (IGG) markers are primer pairs that amplify single-locus sequences that differ in size for two or more alleles. They are attractive for their ease of use for rapid genotyping and their codominant nature. Here, we describe a heuristic algorithm that uses a k-mer-based approach to search two or more genome sequences to locate polymorphic regions suitable for designing candidate IGG marker primers. As input to the IGG pipeline software, the user provides genome sequences and the desired amplicon sizes and size differences. Primer sequences flanking polymorphic insertions/deletions are produced as output. IGG marker files for three sets of genomes, Solanum lycopersicum/Solanum pennellii, Arabidopsis (Arabidopsis thaliana) Columbia-0/Landsberg erecta-0 accessions, and S. lycopersicum/S. pennellii/Solanum tuberosum (three-way polymorphic) are included.

Intron gain by tandem genomic duplication: a novel case in a potato gene encoding RNA-dependent RNA polymerase

The origin and subsequent accumulation of spliceosomal introns are prominent events in the evolution of eukaryotic gene structure. However, the mechanisms underlying intron gain remain unclear because there are few proven cases of recently gained introns. In an RNA-dependent RNA polymerase (RdRp) gene, we found that a tandem duplication occurred after the divergence of potato and its wild relatives among other Solanum plants. The duplicated sequence crosses the intron-exon boundary of the first intron and the second exon. A new intron was detected at this duplicated region, and it includes a small previously exonic segment of the upstream copy of the duplicated sequence and the intronic segment of the downstream copy of the duplicated sequence. The donor site of this new intron was directly obtained from the small previously exonic segment. Most of the splicing signals were inherited directly from the parental intron/exon structure, including a putative branch site, the polypyrimidine tract, the 3' splicing site, two putative exonic splicing enhancers, and the GC contents differed between the intron and exon. In the widely cited model of intron gain by tandem genomic duplication, the duplication of an AGGT-containing exonic segment provides the GT and AG splicing sites for the new intron. Our results illustrate that the tandem duplication model of intron gain should be diverse in terms of obtaining the proper splicing signals.

Untargeted Metabolomics of Tomato Plants after Root-Knot Nematode Infestation

After 2 months from the infestation of tomato plants with the root-knot nematode (RKN) Meloidogyne incognita, we performed a gas chromatography-mass spectrometry untargeted fingerprint analysis for the identification of characteristic metabolites and biomarkers. Principal component analysis, and orthogonal projections to latent structures discriminant analysis suggested dramatic local changes of the plant metabolome. In the case of tomato leaves, β-alanine, phenylalanine, and melibiose were induced in response to RKN stimuli, while ribose, glycerol, myristic acid, and palmitic acid were reduced. For tomato stems, upregulated metabolites were ribose, sucrose, fructose, and glucose, while fumaric acid and glycine were downregulated. The variation in molecular strategies to the infestation of RKNs may play an important role in how Solanum lycopersicum and other plants adapt to nematode parasitic stress.

Gr and hp-1 tomato mutants unveil unprecedented interactions between arbuscular mycorrhizal symbiosis and fruit ripening

Systemic responses to an arbuscular mycorrhizal fungus reveal opposite phenological patterns in two tomato ripening mutants depending whether ethylene or light reception is involved. The availability of tomato ripening mutants has revealed many aspects of the genetics behind fleshy fruit ripening, plant hormones and light signal reception. Since previous analyses revealed that arbuscular mycorrhizal symbiosis influences tomato berry ripening, we wanted to test the hypothesis that an interplay might occur between root symbiosis and fruit ripening. With this aim, we screened seven tomato mutants affected in the ripening process for their responsiveness to the arbuscular mycorrhizal fungus Funneliformis mosseae. Following their phenological responses we selected two mutants for a deeper analysis: Green ripe (Gr), deficient in fruit ethylene perception and high-pigment-1 (hp-1), displaying enhanced light signal perception throughout the plant. We investigated the putative interactions between ripening processes, mycorrhizal establishment and systemic effects using biochemical and gene expression tools. Our experiments showed that both mutants, notwithstanding a normal mycorrhizal phenotype at root level, exhibit altered arbuscule functionality. Furthermore, in contrast to wild type, mycorrhization did not lead to a higher phosphate concentration in berries of both mutants. These results suggest that the mutations considered interfere with arbuscular mycorrhiza inducing systemic changes in plant phenology and fruits metabolism. We hypothesize a cross talk mechanism between AM and ripening processes that involves genes related to ethylene and light signaling.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

Dynamic Metabolic Regulation by a Chromosome Segment from a Wild Relative During Fruit Development in a Tomato Introgression Line, IL8-3

We performed comparative metabolome and transcriptome analyses throughout fruit development using the tomato cultivar M82 and its near-isogenic line IL8-3, with interesting and useful traits such as a high content of soluble solids. Marked differences between M82 and IL8-3 were found not only in ripe fruits but also at 20 days after flowering (DAF) in the hierarchical clustering analysis of the metabolome, whereas patterns were similar between the two genotypes at 10 and 30 DAF. Our metabolome analysis conclusively showed that 20 DAF is an important stage of fruit metabolism and that the Solanum pennellii introgressed region in IL8-3 plays a key role in metabolic changes at this stage. Carbohydrate and amino acid metabolism were found to be promoted in IL8-3 at 20 DAF and the ripening stage, respectively, whereas transcriptome analysis showed no marked differences between the two genotypes, indicating that dynamic metabolic regulation at 20 DAF and the ripening stage was controlled by relatively few genes. The transcript levels of the cell wall invertase (LIN6) and sucrose synthase (TOMSSF) genes in starch and sucrose metabolic pathway and that of the glutamate synthase (SlGOGAT) gene in the amino acid metabolic pathway in IL8-3 fruit were higher than those in M82, and SlGOGAT expression was enhanced under high-sugar conditions. The results suggest that the promotion of carbohydrate metabolism by LIN6 and TOMSSF in IL8-3 fruit at 20 DAF affects SlGOGAT expression and amino acid accumulation via higher sugar concentration at the late stage of fruit development.

A DUF-246 family glycosyltransferase-like gene affects male fertility and the biosynthesis of pectic arabinogalactans

BACKGROUND

Pectins are a group of structurally complex plant cell wall polysaccharides whose biosynthesis and function remain poorly understood. The pectic polysaccharide rhamnogalacturonan-I (RG-I) has two types of arabinogalactan side chains, type-I and type-II arabinogalactans. To date few enzymes involved in the biosynthesis of pectin have been described. Here we report the identification of a highly conserved putative glycosyltransferase encoding gene, Pectic ArabinoGalactan synthesis-Related (PAGR), affecting the biosynthesis of RG-I arabinogalactans and critical for pollen tube growth.

RESULTS

T-DNA insertions in PAGR were identified in Arabidopsis thaliana and were found to segregate at a 1:1 ratio of heterozygotes to wild type. We were unable to isolate homozygous pagr mutants as pagr mutant alleles were not transmitted via pollen. In vitro pollen germination assays revealed reduced rates of pollen tube formation in pollen from pagr heterozygotes. To characterize a loss-of-function phenotype for PAGR, the Nicotiana benthamiana orthologs, NbPAGR-A and B, were transiently silenced using Virus Induced Gene Silencing. NbPAGR-silenced plants exhibited reduced internode and petiole expansion. Cell wall materials from NbPAGR-silenced plants had reduced galactose content compared to the control. Immunological and linkage analyses support that RG-I has reduced type-I arabinogalactan content and reduced branching of the RG-I backbone in NbPAGR-silenced plants. Arabidopsis lines overexpressing PAGR exhibit pleiotropic developmental phenotypes and the loss of apical dominance as well as an increase in RG-I type-II arabinogalactan content.

CONCLUSIONS

Together, results support a function for PAGR in the biosynthesis of RG-I arabinogalactans and illustrate the essential roles of these polysaccharides in vegetative and reproductive plant growth.

Arbuscular mycorrhiza Symbiosis Induces a Major Transcriptional Reprogramming of the Potato SWEET Sugar Transporter Family

Biotrophic microbes feeding on plants must obtain carbon from their hosts without killing the cells. The symbiotic Arbuscular mycorrhizal (AM) fungi colonizing plant roots do so by inducing major transcriptional changes in the host that ultimately also reprogram the whole carbon partitioning of the plant. AM fungi obtain carbohydrates from the root cortex apoplast, in particular from the periarbuscular space that surrounds arbuscules. However, the mechanisms by which cortical cells export sugars into the apoplast for fungal nutrition are unknown. Recently a novel type of sugar transporter, the SWEET, able to perform not only uptake but also efflux from cells was identified. Plant SWEETs have been shown to be involved in the feeding of pathogenic microbes and are, therefore, good candidates to play a similar role in symbiotic associations. Here we have carried out the first phylogenetic and expression analyses of the potato SWEET family and investigated its role during mycorrhiza symbiosis. The potato genome contains 35 SWEETs that cluster into the same four clades defined in Arabidopsis. Colonization of potato roots by the AM fungus Rhizophagus irregularis imposes major transcriptional rewiring of the SWEET family involving, only in roots, changes in 22 of the 35 members. None of the SWEETs showed mycorrhiza-exclusive induction and most of the 12 induced genes belong to the putative hexose transporters of clade I and II, while only two are putative sucrose transporters from clade III. In contrast, most of the repressed transcripts (10) corresponded to clade III SWEETs. Promoter-reporter assays for three of the induced genes, each from one cluster, showed re-localization of expression to arbuscule-containing cells, supporting a role for SWEETs in the supply of sugars at biotrophic interfaces. The complex transcriptional regulation of SWEETs in roots in response to AM fungal colonization supports a model in which symplastic sucrose in cortical cells could be cleaved in the cytoplasm by sucrose synthases or cytoplasmic invertases and effluxed as glucose, but also directly exported as sucrose and then converted into glucose and fructose by cell wall-bound invertases. Precise biochemical, physiological and molecular analyses are now required to profile the role of each potato SWEET in the arbuscular mycorrhizal symbiosis.