Genetical genomics of quality related traits in potato tubers using proteomics.

BACKGROUND: Recent advances in ~omics technologies such as transcriptomics, metabolomics and proteomics along with genotypic profiling have permitted the genetic dissection of complex traits such as quality traits in non-model species. To get more insight into the genetic factors underlying variation in quality traits related to carbohydrate and starch metabolism and cold sweetening, we determined the protein content and composition in potato tubers using 2D-gel electrophoresis in a diploid potato mapping population. Upon analyzing we made sure that the proteins from the patatin family were excluded to ensure a better representation of the other proteins. RESULTS: We subsequently performed pQTL analyses for all other proteins with a sufficient representation in the population and established a relationship between proteins and 26 potato tuber quality traits (e.g. flesh colour, enzymatic discoloration) by co-localization on the genetic map and a direct correlation study of protein abundances and phenotypic traits. Over 1643 unique protein spots were detected in total over the two harvests. We were able to map pQTLs for over 300 different protein spots some of which co-localized with traits such as starch content and cold sweetening. pQTLs were observed on every chromosome although not evenly distributed over the chromosomes. The largest number of pQTLs was found for chromosome 8 and the lowest for chromosome number 10. For some 20 protein spots multiple QTLs were observed. CONCLUSIONS: From this analysis, hotspot areas for protein QTLs were identified on chromosomes three, five, eight and nine. The hotspot on chromosome 3 coincided with a QTL previously identified for total protein content and had more than 23 pQTLs in the region from 70 to 80 cM. Some of the co-localizing protein spots associated with some of the most interesting tuber quality traits were identified, albeit far less than we had anticipated at the onset of the experiments.

Fine mapping of Restorer-of-fertility in pepper (Capsicum annuum L.) identified a candidate gene encoding a pentatricopeptide repeat (PPR)-containing protein.

KEY MESSAGE: Using fine mapping techniques, the genomic region co-segregating with Restorer - of - fertility ( Rf ) in pepper was delimited to a region of 821 kb in length. A PPR gene in this region, CaPPR6 , was identified as a strong candidate for Rf based on expression pattern and characteristics of encoding sequence. Cytoplasmic-genic male sterility (CGMS) has been used for the efficient production of hybrid seeds in peppers (Capsicum annuum L.). Although the mitochondrial candidate genes that might be responsible for cytoplasmic male sterility (CMS) have been identified, the nuclear Restorer-of-fertility (Rf) gene has not been isolated. To identify the genomic region co-segregating with Rf in pepper, we performed fine mapping using an Rf-segregating population consisting of 1068 F2 individuals, based on BSA-AFLP and a comparative mapping approach. Through six cycles of chromosome walking, the co-segregating region harboring the Rf locus was delimited to be within 821 kb of sequence. Prediction of expressed genes in this region based on transcription analysis revealed four candidate genes. Among these, CaPPR6 encodes a pentatricopeptide repeat (PPR) protein with PPR motifs that are repeated 14 times. Characterization of the CaPPR6 protein sequence, based on alignment with other homologs, showed that CaPPR6 is a typical Rf-like (RFL) gene reported to have undergone diversifying selection during evolution. A marker developed from a sequence near CaPPR6 showed a higher prediction rate of the Rf phenotype than those of previously developed markers when applied to a panel of breeding lines of diverse origin. These results suggest that CaPPR6 is a strong candidate for the Rf gene in pepper.

Integration of multi-omics data for prediction of phenotypic traits using random forest.

BACKGROUND: In order to find genetic and metabolic pathways related to phenotypic traits of interest, we analyzed gene expression data, metabolite data obtained with GC-MS and LC-MS, proteomics data and a selected set of tuber quality phenotypic data from a diploid segregating mapping population of potato. In this study we present an approach to integrate these ~ omics data sets for the purpose of predicting phenotypic traits. This gives us networks of relatively small sets of interrelated ~ omics variables that can predict, with higher accuracy, a quality trait of interest. RESULTS: We used Random Forest regression for integrating multiple ~ omics data for prediction of four quality traits of potato: tuber flesh colour, DSC onset, tuber shape and enzymatic discoloration. For tuber flesh colour beta-carotene hydroxylase and zeaxanthin epoxidase were ranked first and forty-fourth respectively both of which have previously been associated with flesh colour in potato tubers. Combining all the significant genes, LC-peaks, GC-peaks and proteins, the variation explained was 75 %, only slightly more than what gene expression or LC-MS data explain by themselves which indicates that there are correlations among the variables across data sets. For tuber shape regressed on the gene expression, LC-MS, GC-MS and proteomics data sets separately, only gene expression data was found to explain significant variation. For DSC onset, we found 12 significant gene expression, 5 metabolite levels (GC) and 2 proteins that are associated with the trait. Using those 19 significant variables, the variation explained was 45 %. Expression QTL (eQTL) analyses showed many associations with genomic regions in chromosome 2 with also the highest explained variation compared to other chromosomes. Transcriptomics and metabolomics analysis on enzymatic discoloration after 5 min resulted in 420 significant genes and 8 significant LC metabolites, among which two were putatively identified as caffeoylquinic acid methyl ester and tyrosine. CONCLUSIONS: In this study, we made a strategy for selecting and integrating multiple ~ omics data using random forest method and selected representative individual peaks for networks based on eQTL, mQTL or pQTL information. Network analysis was done to interpret how a particular trait is associated with gene expression, metabolite and protein data.

Down regulation of StGA3ox genes in potato results in altered GA content and affect plant and tuber growth characteristics.

GA biosynthesis and catabolism has been shown to play an important role in regulating tuberization in potato. Active GAs are inactivated in the stolon tips shortly after induction to tuberization. Overexpression of a GA inactivation gene results in an earlier tuberization phenotype, while reducing expression of the same gene results in delayed tuberization. In addition, overexpression of genes involved in GA biosynthesis results in delayed tuberization, while decreased expression of those genes results in earlied tuberization. The final step in GA biosynthesis is catalysed by StGA3ox1 and StGA3ox2 activity, that convert inactive forms of GA into active GA1 and GA4. In this study we cloned StGA3ox2 gene in an RNAi construct and used this construct to transform potato plants. The StGA3ox2 silenced plants were smaller and had shorter internodes. In addition, we assayed the concentrations of various GAs in the transgenic plants and showed an altered GA content. No difference was observed on the time point of tuber initiation. However, the transgenic clones had increased number of tubers with the same yield, resulting in smaller average tuber weight. In addition, we cloned the promoter of StGA3ox2 to direct expression of the GUS reporter gene to visualize the sites of GA biosynthesis in the potato plant. Finally, we discuss how changes of several GA levels can have an impact on shoot, stolon and tuber development, as well as the possible mechanisms that mediate feed-forward and feed-back regulation loops in the GA biosynthetic pathway in potato.

Naturally occurring allele diversity allows potato cultivation in northern latitudes.

Potato (Solanum tuberosum L.) originates from the Andes and evolved short-day-dependent tuber formation as a vegetative propagation strategy. Here we describe the identification of a central regulator underlying a major-effect quantitative trait locus for plant maturity and initiation of tuber development. We show that this gene belongs to the family of DOF (DNA-binding with one finger) transcription factors and regulates tuberization and plant life cycle length, by acting as a mediator between the circadian clock and the StSP6A mobile tuberization signal. We also show that natural allelic variants evade post-translational light regulation, allowing cultivation outside the geographical centre of origin of potato. Potato is a member of the Solanaceae family and is one of the world's most important food crops. This annual plant originates from the Andean regions of South America. Potato develops tubers from underground stems called stolons. Its equatorial origin makes potato essentially short-day dependent for tuberization and potato will not make tubers in the long-day conditions of spring and summer in the northern latitudes. When introduced in temperate zones, wild material will form tubers in the course of the autumnal shortening of day-length. Thus, one of the first selected traits in potato leading to a European potato type is likely to have been long-day acclimation for tuberization. Potato breeders can exploit the naturally occurring variation in tuberization onset and life cycle length, allowing varietal breeding for different latitudes, harvest times and markets.

The PIN family of proteins in potato and their putative role in tuberization.

The PIN family of trans-membrane proteins mediates auxin efflux throughout the plant and during various phases of plant development. In Arabidopsis thaliana, the PIN family comprised of 8 members, divided into "short" and "long" PINs according to the length of the hydrophilic domain of the protein. Based on sequence homology using the recently published potato genome sequence (Solanum tuberosum group Phureja) we identified ten annotated potato StPIN genes. Mining the publicly available gene expression data, we constructed a catalog tissue specificity of StPIN gene expression, focusing on the process of tuberization. A total of four StPIN genes exhibited increased expression 4 days after tuber induction, prior to the onset of stolon swelling. For two PIN genes, StPIN4 and StPIN2, promoter sequences were cloned and fused to the GUS reporter protein to study tissue specificity in more detail. StPIN4 promoter driven GUS staining was detected in the flower stigma, in the flower style, below the ovary and petals, in the root tips, in the vascular tissue of the stolons and in the tuber parenchyma cells. StPIN2 promoter driven GUS staining was detected in flower buds, in the vascular tissue of the swelling stolons and in the storage parenchyma of the growing tubers. Based on our results, we postulate a role for the StPINs in redistributing auxin in the swelling stolon during early events in tuber development.

The effects of auxin and strigolactones on tuber initiation and stolon architecture in potato.

Various transcriptional networks and plant hormones have been implicated in controlling different aspects of potato tuber formation. Due to its broad impact on many plant developmental processes, a role for auxin in tuber initiation has been suggested but never fully resolved. Here, auxin concentrations were measured throughout the plant prior to and during the process of tuber formation. Auxin levels increase dramatically in the stolon prior to tuberization and remain relatively high during subsequent tuber growth, suggesting a promoting role for auxin in tuber formation. Furthermore, in vitro tuberization experiments showed higher levels of tuber formation from axillary buds of explants where the auxin source (stolon tip) had been removed. This phenotype could be rescued by application of auxin on the ablated stolon tips. In addition, a synthetic strigolactone analogue applied on the basal part of the stolon resulted in fewer tubers. The experiments indicate that a system for the production and directional transport of auxin exists in stolons and acts synergistically with strigolactones to control the outgrowth of the axillary stolon buds, similar to the control of above-ground shoot branching.

Organ specificity and transcriptional control of metabolic routes revealed by expression QTL profiling of source--sink tissues in a segregating potato population.

BACKGROUND: With the completion of genome sequences belonging to some of the major crop plants, new challenges arise to utilize this data for crop improvement and increased food security. The field of genetical genomics has the potential to identify genes displaying heritable differential expression associated to important phenotypic traits. Here we describe the identification of expression QTLs (eQTLs) in two different potato tissues of a segregating potato population and query the potato genome sequence to differentiate between cis- and trans-acting eQTLs in relation to gene subfunctionalization. RESULTS: Leaf and tuber samples were analysed and screened for the presence of conserved and tissue dependent eQTLs. Expression QTLs present in both tissues are predominantly cis-acting whilst for tissue specific QTLs, the percentage of trans-acting QTLs increases. Tissue dependent eQTLs were assigned to functional classes and visualized in metabolic pathways. We identified a potential regulatory network on chromosome 10 involving genes crucial for maintaining circadian rhythms and controlling clock output genes. In addition, we show that the type of genetic material screened and sampling strategy applied, can have a high impact on the output of genetical genomics studies. CONCLUSIONS: Identification of tissue dependent regulatory networks based on mapped differential expression not only gives us insight in tissue dependent gene subfunctionalization but brings new insights into key biological processes and delivers targets for future haplotyping and genetic marker development.

Data integration and network reconstruction with ~omics data using Random Forest regression in potato.

In the post-genomic era, high-throughput technologies have led to data collection in fields like transcriptomics, metabolomics and proteomics and, as a result, large amounts of data have become available. However, the integration of these ~omics data sets in relation to phenotypic traits is still problematic in order to advance crop breeding. We have obtained population-wide gene expression and metabolite (LC-MS) data from tubers of a diploid potato population and present a novel approach to study the various ~omics datasets to allow the construction of networks integrating gene expression, metabolites and phenotypic traits. We used Random Forest regression to select subsets of the metabolites and transcripts which show association with potato tuber flesh color and enzymatic discoloration. Network reconstruction has led to the integration of known and uncharacterized metabolites with genes associated with the carotenoid biosynthesis pathway. We show that this approach enables the construction of meaningful networks with regard to known and unknown components and metabolite pathways.

Identification of alleles of carotenoid pathway genes important for zeaxanthin accumulation in potato tubers.

We have investigated the genetics and molecular biology of orange flesh colour in potato (Solanum tuberosum L.). To this end the natural diversity in three genes of the carotenoid pathway was assessed by SNP analyses. Association analysis was performed between SNP haplotypes and flesh colour phenotypes in diploid and tetraploid potato genotypes. We observed that among eleven beta-carotene hydroxylase 2 (Chy2) alleles only one dominant allele has a major effect, changing white into yellow flesh colour. In contrast, none of the lycopene epsilon cyclase (Lcye) alleles seemed to have a large effect on flesh colour. Analysis of zeaxanthin epoxidase (Zep) alleles showed that all (diploid) genotypes with orange tuber flesh were homozygous for one specific Zep allele. This Zep allele showed a reduced level of expression. The complete genomic sequence of the recessive Zep allele, including the promoter, was determined, and compared with the sequence of other Zep alleles. The most striking difference was the presence of a non-LTR retrotransposon sequence in intron 1 of the recessive Zep allele, which was absent in all other Zep alleles investigated. We hypothesise that the presence of this large sequence in intron 1 caused the lower expression level, resulting in reduced Zep activity and accumulation of zeaxanthin. Only genotypes combining presence of the dominant Chy2 allele with homozygosity for the recessive Zep allele produced orange-fleshed tubers that accumulated large amounts of zeaxanthin.

From QTL to candidate gene: genetical genomics of simple and complex traits in potato using a pooling strategy.

BACKGROUND: Utilization of the natural genetic variation in traditional breeding programs remains a major challenge in crop plants. The identification of candidate genes underlying, or associated with, phenotypic trait QTLs is desired for effective marker assisted breeding. With the advent of high throughput -omics technologies, screening of entire populations for association of gene expression with targeted traits is becoming feasible but remains costly. Here we present the identification of novel candidate genes for different potato tuber quality traits by employing a pooling approach reducing the number of hybridizations needed. Extreme genotypes for a quantitative trait are collected and the RNA from contrasting bulks is then profiled with the aim of finding differentially expressed genes. RESULTS: We have successfully implemented the pooling strategy for potato quality traits and identified candidate genes associated with potato tuber flesh color and tuber cooking type. Elevated expression level of a dominant allele of the beta-carotene hydroxylase (bch) gene was associated with yellow flesh color through mapping of the gene under a major QTL for flesh color on chromosome 3. For a second trait, a candidate gene with homology to a tyrosine-lysine rich protein (TLRP) was identified based on allele specificity of the probe on the microarray. TLRP was mapped on chromosome 9 in close proximity to a QTL for potato cooking type strengthening its significance as a candidate gene. Furthermore, we have performed a profiling experiment targeting a polygenic trait, by pooling individual genotypes based both on phenotypic and marker data, allowing the identification of candidate genes associated with the two different linkage groups. CONCLUSIONS: A pooling approach for RNA-profiling with the aim of identifying novel candidate genes associated with tuber quality traits was successfully implemented. The identified candidate genes for tuber flesh color (bch) and cooking type (tlrp) can provide useful markers for breeding schemes in the future. Strengths and limitations of the approach are discussed.

Genes driving potato tuber initiation and growth: identification based on transcriptional changes using the POCI array.

The increasing amount of available expressed gene sequence data makes whole-transcriptome analysis of certain crop species possible. Potato currently has the second largest number of publicly available expressed sequence tag (EST) sequences among the Solanaceae. Most of these ESTs, plus other proprietary sequences, were combined and used to generate a unigene assembly. The set of 246,182 sequences produced 46,345 unigenes, which were used to design a 44K 60-mer oligo array (Potato Oligo Chip Initiative: POCI). In this study, we attempt to identify genes controlling and driving the process of tuber initiation and growth by implementing large-scale transcriptional changes using the newly developed POCI array. Major gene expression profiles could be identified exhibiting differential expression at key developmental stages. These profiles were associated with functional roles in cell division and growth. A subset of genes involved in the regulation of the cell cycle, based on their Gene Ontology classification, exhibit a clear transient upregulation at tuber onset indicating increased cell division during these stages. The POCI array allows the study of potato gene expression on a much broader level than previously possible and will greatly enhance analysis of transcriptional control mechanisms in a wide range of potato research areas. POCI sequence and annotation data are publicly available through the POCI database ( http://pgrc.ipk-gatersleben.de/poci ).

StGA2ox1 is induced prior to stolon swelling and controls GA levels during potato tuber development.

The formation and growth of a potato (Solanum tuberosum) tuber is a complex process regulated by different environmental signals and plant hormones. In particular, the action of gibberellins (GAs) has been implicated in different aspects of potato tuber formation. Here we report on the isolation and functional analysis of a potato GA 2-oxidase gene (StGA2ox1) and its role in tuber formation. StGA2ox1 is upregulated during the early stages of potato tuber development prior to visible swelling and is predominantly expressed in the subapical region of the stolon and growing tuber. 35S-over-expression transformants exhibit a dwarf phenotype, reduced stolon growth and earlier in vitro tuberization. Transgenic plants with reduced expression levels of StGA2ox1 showed normal plant growth, an altered stolon swelling phenotype and delayed in vitro tuberization. Tubers of the StGA2ox1 suppression clones contain increased levels of GA20, indicating altered GA metabolism. We propose a role for StGA2ox1 in early tuber initiation by modifying GA levels in the subapical stolon region at the onset of tuberization, thereby facilitating normal tuber development and growth.

Unravelling enzymatic discoloration in potato through a combined approach of candidate genes, QTL, and expression analysis.

Enzymatic discoloration (ED) of potato tubers was investigated in an attempt to unravel the underlying genetic factors. Both enzyme and substrate concentration have been reported to influence the degree of discoloration and as such this trait can be regarded as polygenic. The diploid mapping population C x E, consisting of 249 individuals, was assayed for the degree of ED and levels of chlorogenic acid and tyrosine. Using this data, Quantitative Trait Locus (QTL) analysis was performed. Three QTLs for ED have been found on parental chromosomes C3, C8, E1, and E8. For chlorogenic acid a QTL has been identified on C2 and for tyrosine levels, a QTL has been detected on C8. None of the QTLs overlap, indicating the absence of genetic correlations between these components underlying ED, in contrast to earlier reports in literature. An obvious candidate gene for the QTL for ED on Chromosome 8 is polyphenol oxidase (PPO), which was previously mapped on chromosome 8. With gene-specific primers for PPO gene POT32 a CAPS marker was developed. Three different alleles (POT32-1, -2, and -3) could be discriminated. The segregating POT32 alleles were used to map the POT32 CAPS marker and QTL analysis was redone, showing that POT32 coincides with the QTL peak. A clear correlation between allele combinations and degree of discoloration was observed. In addition, analysis of POT32 gene expression in a subset of genotypes indicated a correlation between the level of gene expression and allele composition. On average, genotypes having two copies of allele 1 had both the highest degree of discoloration as well as the highest level of POT32 gene expression.

Tuber on a chip: differential gene expression during potato tuber development.

Potato tuber development has proven to be a valuable model system for studying underground sink organ formation. Research on this topic has led to the identification of many genes involved in this complex process and has aided in the unravelling of the mechanisms underlying starch synthesis. However, less attention has been paid to the biochemical pathways of other important metabolites or to the changing metabolic fluxes occurring during potato tuber development. In this paper, we describe the construction of a potato complementary DNA (cDNA) microarray specifically designed for genes involved in processes related to tuber development and tuber quality traits. We present expression profiles of 1315 cDNAs during tuber development where the predominant profiles were strong up- and down-regulation. Gene expression profiles showing transient increases or decreases were less abundantly represented and followed more moderate changes, mainly during tuber initiation. In addition to the confirmation of gene expression patterns during tuber development, many novel differentially expressed genes were identified and are considered as candidate genes for direct involvement in potato tuber development. A detailed analysis of starch metabolism genes provided a unique overview of expression changes during tuber development. Characteristic expression profiles were often clearly different between gene family members. A link between differential gene expression during tuber development and potato tissue specificity is described. This dataset provides a firm basis for the identification of key regulatory genes in a number of metabolic pathways that may provide researchers with new tools to achieve breeding goals for use in industrial applications.