Ruth A. Heinz

Transcriptomic identification of candidate genes involved in sunflower responses to chilling and salt stresses based on cDNA microarray analysis

BackgroundConsidering that sunflower production is expanding to arid regions, tolerance to abiotic stresses as drought, low temperatures and salinity arises as one of the main constrains nowadays. Differential organ-specific sunflower ESTs (expressed sequence tags) were previously generated by a subtractive hybridization method that included a considerable number of putative abiotic stress associated sequences. The objective of this work is to analyze concerted gene expression profiles of organ-specific ESTs by fluorescence microarray assay, in response to high sodium chloride concentration and chilling treatments with the aim to identify and follow up candidate genes for early responses to abiotic stress in sunflower.ResultsAbiotic-related expressed genes were the target of this characterization through a gene expression analysis using an organ-specific cDNA fluorescence microarray approach in response to high salinity and low temperatures. The experiment included three independent replicates from leaf samples. We analyzed 317 unigenes previously isolated from differential organ-specific cDNA libraries from leaf, stem and flower at R1 and R4 developmental stage. A statistical analysis based on mean comparison by ANOVA and ordination by Principal Component Analysis allowed the detection of 80 candidate genes for either salinity and/or chilling stresses. Out of them, 50 genes were up or down regulated under both stresses, supporting common regulatory mechanisms and general responses to chilling and salinity. Interestingly 15 and 12 sequences were up regulated or down regulated specifically in one stress but not in the other, respectively. These genes are potentially involved in different regulatory mechanisms including transcription/translation/protein degradation/protein folding/ROS production or ROS-scavenging. Differential gene expression patterns were confirmed by qRT-PCR for 12.5% of the microarray candidate sequences.ConclusionEighty genes isolated from organ-specific cDNA libraries were identified as candidate genes for sunflower early response to low temperatures and salinity. Microarray profiling of chilling and NaCl-treated sunflower leaves revealed dynamic changes in transcript abundance, including transcription factors, defense/stress related proteins, and effectors of homeostasis, all of which highlight the complexity of both stress responses. This study not only allowed the identification of common transcriptional changes to both stress conditions but also lead to the detection of stress-specific genes not previously reported in sunflower. This is the first organ-specific cDNA fluorescence microarray study addressing a simultaneous evaluation of concerted transcriptional changes in response to chilling and salinity stress in cultivated sunflower.

Comparison of predictive methods and biological validation for qPCR reference genes in sunflower leaf senescence transcript analysis

The selection and validation of reference genes constitute a key point for gene expression analysis based on qPCR, requiring efficient normalization approaches. In this work, the expression profiles of eight genes were evaluated to identify novel reference genes for transcriptional studies associated to the senescence process in sunflower. Three alternative strategies were applied for the evaluation of gene expression stability in leaves of different ages and exposed to different treatments affecting the senescence process: algorithms implemented in geNorm, BestKeeper software, and the fitting of a statistical linear mixed model (LMModel). The results show that geNorm suggested the use of all combined genes, although identifying α-TUB1 as the most stable expressing gene. BestKeeper revealed α-TUB and β-TUB as stable genes, scoring β-TUB as the most stable one. The statistical LMModel identified α-TUB, actin, PEP, and EF-1α as stable genes in this order. The model-based approximation allows not only the estimation of systematic changes in gene expression, but also the identification of sources of random variation through the estimation of variance components, considering the experimental design applied. Validation of α-TUB and EF-1α as reference genes for expression studies of three sunflower senescence associated genes showed that the first one was more stable for the assayed conditions. We conclude that, when biological replicates are available, LMModel allows a more reliable selection under the assayed conditions. This study represents the first analysis of identification and validation of genuine reference genes for use as internal control in qPCR expression studies in sunflower, experimentally validated throughout six different controlled leaf senescence conditions.

Identification of Single Nucleotide Polymorphisms and analysis of Linkage Disequilibrium in sunflower elite inbred lines using the candidate gene approach

BackgroundAssociation analysis is a powerful tool to identify gene loci that may contribute to phenotypic variation. This includes the estimation of nucleotide diversity, the assessment of linkage disequilibrium structure (LD) and the evaluation of selection processes. Trait mapping by allele association requires a high-density map, which could be obtained by the addition of Single Nucleotide Polymorphisms (SNPs) and short insertion and/or deletions (indels) to SSR and AFLP genetic maps. Nucleotide diversity analysis of randomly selected candidate regions is a promising approach for the success of association analysis and fine mapping in the sunflower genome. Moreover, knowledge of the distance over which LD persists, in agronomically meaningful sunflower accessions, is important to establish the density of markers and the experimental design for association analysis.ResultsA set of 28 candidate genes related to biotic and abiotic stresses were studied in 19 sunflower inbred lines. A total of 14,348 bp of sequence alignment was analyzed per individual. In average, 1 SNP was found per 69 nucleotides and 38 indels were identified in the complete data set. The mean nucleotide polymorphism was moderate (θ = 0.0056), as expected for inbred materials. The number of haplotypes per region ranged from 1 to 9 (mean = 3.54 ± 1.88). Model-based population structure analysis allowed detection of admixed individuals within the set of accessions examined. Two putative gene pools were identified (G1 and G2), with a large proportion of the inbred lines being assigned to one of them (G1). Consistent with the absence of population sub-structuring, LD for G1 decayed more rapidly (r2 = 0.48 at 643 bp; trend line, pooled data) than the LD trend line for the entire set of 19 individuals (r2 = 0.64 for the same distance).ConclusionKnowledge about the patterns of diversity and the genetic relationships between breeding materials could be an invaluable aid in crop improvement strategies. The relatively high frequency of SNPs within the elite inbred lines studied here, along with the predicted extent of LD over distances of 100 kbp (r2~0.1) suggest that high resolution association mapping in sunflower could be achieved with marker densities lower than those usually reported in the literature.

Differential representation of sunflower ESTs in enriched organ-specific cDNA libraries in a small scale sequencing project

BackgroundSubtractive hybridization methods are valuable tools for identifying differentially regulated genes in a given tissue avoiding redundant sequencing of clones representing the same expressed genes, maximizing detection of low abundant transcripts and thus, affecting the efficiency and cost effectiveness of small scale cDNA sequencing projects aimed to the specific identification of useful genes for breeding purposes. The objective of this work is to evaluate alternative strategies to high-throughput sequencing projects for the identification of novel genes differentially expressed in sunflower as a source of organ-specific genetic markers that can be functionally associated to important traits.ResultsDifferential organ-specific ESTs were generated from leaf, stem, root and flower bud at two developmental stages (R1 and R4). The use of different sources of RNA as tester and driver cDNA for the construction of differential libraries was evaluated as a tool for detection of rare or low abundant transcripts. Organ-specificity ranged from 75 to 100% of non-redundant sequences in the different cDNA libraries. Sequence redundancy varied according to the target and driver cDNA used in each case. The R4 flower cDNA library was the less redundant library with 62% of unique sequences. Out of a total of 919 sequences that were edited and annotated, 318 were non-redundant sequences. Comparison against sequences in public databases showed that 60% of non-redundant sequences showed significant similarity to known sequences. The number of predicted novel genes varied among the different cDNA libraries, ranging from 56% in the R4 flower to 16 % in the R1 flower bud library. Comparison with sunflower ESTs on public databases showed that 197 of non-redundant sequences (60%) did not exhibit significant similarity to previously reported sunflower ESTs. This approach helped to successfully isolate a significant number of new reported sequences putatively related to responses to important agronomic traits and key regulatory and physiological genes.ConclusionsThe application of suppressed subtracted hybridization technology not only enabled the cost effective isolation of differentially expressed sequences but it also allowed the identification of novel sequences in sunflower from a relative small number of analyzed sequences when compared to major sequencing projects.

Association mapping in sunflower for sclerotinia head rot resistance

BackgroundSclerotinia Head Rot (SHR) is one of the most damaging diseases of sunflower in Europe, Argentina, and USA, causing average yield reductions of 10 to 20 %, but leading to total production loss under favorable environmental conditions for the pathogen. Association Mapping (AM) is a promising choice for Quantitative Trait Locus (QTL) mapping, as it detects relationships between phenotypic variation and gene polymorphisms in existing germplasm without development of mapping populations. This article reports the identification of QTL for resistance to SHR based on candidate gene AM.ResultsA collection of 94 sunflower inbred lines were tested for SHR under field conditions using assisted inoculation with the fungal pathogen Sclerotinia sclerotiorum. Given that no biological mechanisms or biochemical pathways have been clearly identified for SHR, 43 candidate genes were selected based on previous transcript profiling studies in sunflower and Brassica napus infected with S. sclerotiorum. Associations among SHR incidence and haplotype polymorphisms in 16 candidate genes were tested using Mixed Linear Models (MLM) that account for population structure and kinship relationships. This approach allowed detection of a significant association between the candidate gene HaRIC_B and SHR incidence (P < 0.01), accounting for a SHR incidence reduction of about 20 %.ConclusionsThese results suggest that AM will be useful in dissecting other complex traits in sunflower, thus providing a valuable tool to assist in crop breeding.

Development, characterization and experimental validation of a cultivated sunflower (Helianthus annuus L.) gene expression oligonucleotide microarray.

Oligonucleotide-based microarrays with accurate gene coverage represent a key strategy for transcriptional studies in orphan species such as sunflower, H. annuus L., which lacks full genome sequences. The goal of this study was the development and functional annotation of a comprehensive sunflower unigene collection and the design and validation of a custom sunflower oligonucleotide-based microarray. A large scale EST (>130,000 ESTs) curation, assembly and sequence annotation was performed using Blast2GO (www.blast2go.de). The EST assembly comprises 41,013 putative transcripts (12,924 contigs and 28,089 singletons). The resulting Sunflower Unigen Resource (SUR version 1.0) was used to design an oligonucleotide-based Agilent microarray for cultivated sunflower. This microarray includes a total of 42,326 features: 1,417 Agilent controls, 74 control probes for sunflower replicated 10 times (740 controls) and 40,169 different non-control probes. Microarray performance was validated using a model experiment examining the induction of senescence by water deficit. Pre-processing and differential expression analysis of Agilent microarrays was performed using the Bioconductor limma package. The analyses based on p-values calculated by eBayes (p<0.01) allowed the detection of 558 differentially expressed genes between water stress and control conditions; from these, ten genes were further validated by qPCR. Over-represented ontologies were identified using FatiScan in the Babelomics suite. This work generated a curated and trustable sunflower unigene collection, and a custom, validated sunflower oligonucleotide-based microarray using Agilent technology. Both the curated unigene collection and the validated oligonucleotide microarray provide key resources for sunflower genome analysis, transcriptional studies, and molecular breeding for crop improvement.

Integrating transcriptomic and metabolomic analysis to understand natural leaf senescence in sunflower.

Leaf senescence is a complex process, which has dramatic consequences on crop yield. In sunflower, gap between potential and actual yields reveals the economic impact of senescence. Indeed, sunflower plants are incapable of maintaining their green leaf area over sustained periods. This study characterizes the leaf senescence process in sunflower through a systems biology approach integrating transcriptomic and metabolomic analyses: plants being grown under both glasshouse and field conditions. Our results revealed a correspondence between profile changes detected at the molecular, biochemical and physiological level throughout the progression of leaf senescence measured at different plant developmental stages. Early metabolic changes were detected prior to anthesis and before the onset of the first senescence symptoms, with more pronounced changes observed when physiological and molecular variables were assessed under field conditions. During leaf development, photosynthetic activity and cell growth processes decreased, whereas sucrose, fatty acid, nucleotide and amino acid metabolisms increased. Pathways related to nutrient recycling processes were also up-regulated. Members of the NAC, AP2-EREBP, HB, bZIP and MYB transcription factor families showed high expression levels, and their expression level was highly correlated, suggesting their involvement in sunflower senescence. The results of this study thus contribute to the elucidation of the molecular mechanisms involved in the onset and progression of leaf senescence in sunflower leaves as well as to the identification of candidate genes involved in this process.

Single nucleotide polymorphism genotyping by heteroduplex analysis in sunflower (Helianthus annuus L.)

Single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) are increasingly used for cultivar identification, construction of genetic maps, genetic diversity assessment, association mapping and marker-assisted breeding. Although there are several highly sensitive methods for the detection of polymorphisms, most of them are often beyond the budget of medium-throughput academic laboratories or seed companies. Heteroduplex analysis by enzymatic cleavage (CEL1CH) or denaturing high-performance liquid chromatography (dHPLC) has been successfully used to examine genetic variation in several plant and animal species. In this work, we assess and compare the performance of both methods in sunflower by genotyping SNPs from a set of 24 selected polymorphic candidate genes. The CEL1CH method allowed us to accurately detect allele differences in 10 out of 24 regions using an in-house prepared CEL1 enzyme (celery single strand endonuclease 1, Apium graveolens L.). Similarly, a total of 11 regions were successfully optimized for dHPLC analysis. As a scaling-up approach, both strategies were tested to genotype either 42 SNPs/indels in 22 sunflower accessions from the local germplasm bank or 33 SNPs/indels in 90 recombinant inbred lines (RILs) for genetic mapping purposes. Summarizing, a total of 601 genotypes were efficiently analyzed either with CEL1CH (110) or dHPCL (491). In conclusion, CEL1CH and dHPLC proved to be robust, complementary methods, allowing medium-scale laboratories to scale up the number of both SNPs and individuals to be included in genetic studies and targeted germplasm diversity characterization (EcoTILLING).

Cloning and functional characterization of a fructan 1-exohydrolase (1-FEH) in the cold tolerant Patagonian species Bromus pictus.

Fructans are fructose polymers synthesized in a wide range of species such as bacteria, fungi and plants. Fructans are synthesized by fructosyltransferases (FTs) and depolymerized by fructan exohydrolases (FEHs). Bromus pictus is a graminean decaploid species from the Patagonian region of Argentina, which accumulates large amounts of fructans even at temperate temperatures. The first gene isolated from B. pictus fructan metabolism was a putative sucrose:fructan 6-fructosyltransferase (6-SFT). Here, a complete cDNA of the first fructan exohydrolase (FEH) from B. pictus (Bp1-FEHa) was isolated using RT-PCR strategies. The Bp1-FEHa encoding gene is present as a single copy in B. pictus genome. Functional characterization in Pichia pastoris confirmed Bp1-FEHa is a fructan exohydrolase with predominant activity towards β-(2-1) linkages. Its expression was analyzed in different leaf sections, showing the highest expression levels in the second section of the sheath and the tip of the blade. Bp1-FEHa expression was studied along with FEH and FT activities and fructan accumulation profile in response to chilling conditions during a 7-day time course experiment. Bp1-FEHa expression and FEH activity followed a similar pattern in response to low temperatures, especially in basal sections of the sheaths. In these sections the FEH and FT activities were particularly high and they were significantly correlated to fructan accumulation profile, along with cold treatment.

Network and biosignature analysis for the integration of transcriptomic and metabolomic data to characterize leaf senescence process in sunflower

BackgroundIn recent years, high throughput technologies have led to an increase of datasets from omics disciplines allowing the understanding of the complex regulatory networks associated with biological processes. Leaf senescence is a complex mechanism controlled by multiple genetic and environmental variables, which has a strong impact on crop yield. Transcription factors (TFs) are key proteins in the regulation of gene expression, regulating different signaling pathways; their function is crucial for triggering and/or regulating different aspects of the leaf senescence process. The study of TF interactions and their integration with metabolic profiles under different developmental conditions, especially for a non-model organism such as sunflower, will open new insights into the details of gene regulation of leaf senescence.ResultsWeighted Gene Correlation Network Analysis (WGCNA) and BioSignature Discoverer (BioSD, Gnosis Data Analysis, Heraklion, Greece) were used to integrate transcriptomic and metabolomic data. WGCNA allowed the detection of 10 metabolites and 13 TFs whereas BioSD allowed the detection of 1 metabolite and 6 TFs as potential biomarkers. The comparative analysis demonstrated that three transcription factors were detected through both methodologies, highlighting them as potentially robust biomarkers associated with leaf senescence in sunflower.ConclusionsThe complementary use of network and BioSignature Discoverer analysis of transcriptomic and metabolomic data provided a useful tool for identifying candidate genes and metabolites which may have a role during the triggering and development of the leaf senescence process. The WGCNA tool allowed us to design and test a hypothetical network in order to infer relationships across selected transcription factor and metabolite candidate biomarkers involved in leaf senescence, whereas BioSignature Discoverer selected transcripts and metabolites which discriminate between different ages of sunflower plants. The methodology presented here would help to elucidate and predict novel networks and potential biomarkers of leaf senescence in sunflower.