Comparison of ddRADseq and EUChip60K SNP genotyping systems for population genetics and genomic selection in Eucalyptus dunnii (Maiden).

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Co-Expression Networks in Sunflower: Harnessing the Power of Multi-Study Transcriptomic Public Data to Identify and Categorize Candidate Genes for Fungal Resistance.

Fungal plant diseases are a major threat to food security worldwide. Current efforts to identify and list loci involved in different biological processes are more complicated than originally thought, even when complete genome assemblies are available. Despite numerous experimental and computational efforts to characterize gene functions in plants, about ~40% of protein-coding genes in the model plant Arabidopsis thaliana L. are still not categorized in the Gene Ontology (GO) Biological Process (BP) annotation. In non-model organisms, such as sunflower (Helianthus annuus L.), the number of BP term annotations is far fewer, ~22%. In the current study, we performed gene co-expression network analysis using eight terabytes of public transcriptome datasets and expression-based functional prediction to categorize and identify loci involved in the response to fungal pathogens. We were able to construct a reference gene network of healthy green tissue (GreenGCN) and a gene network of healthy and stressed root tissues (RootGCN). Both networks achieved robust, high-quality scores on the metrics of guilt-by-association and selective constraints versus gene connectivity. We were able to identify eight modules enriched in defense functions, of which two out of the three modules in the RootGCN were also conserved in the GreenGCN, suggesting similar defense-related expression patterns. We identified 16 WRKY genes involved in defense related functions and 65 previously uncharacterized loci now linked to defense response. In addition, we identified and classified 122 loci previously identified within QTLs or near candidate loci reported in GWAS studies of disease resistance in sunflower linked to defense response. All in all, we have implemented a valuable strategy to better describe genes within specific biological processes.

Double Digest Restriction-Site Associated DNA Sequencing (ddRADseq) Technology.

Double digest restriction-site associated DNA sequencing (ddRADseq) technology combines genome reduced representation by digestion with two restriction enzymes and next generation sequencing (NGS) to obtain thousands of markers (SNP, SSR, and InDels) and genotype tens to hundreds of samples simultaneously. In this chapter, we describe a 96-plex derived ddRADseq protocol that can be set up to obtain different depth of coverage per locus and can be exploited to model and non-model plant species.

ddRADseq-mediated detection of genetic variants in sugarcane.

KEY MESSAGE: The article presents an optimization of the key parameters for the identification of SNPs in sugarcane using a GBS protocol based on two Illumina NextSeq and NovaSeq platforms. Sugarcane (Saccharum sp.), a world-wide known feedstock for sugar production, bioethanol, and energy, has an extremely complex genome, being highly polyploid and aneuploid. A double-digestion restriction site-associated DNA sequencing protocol (ddRADseq) was tested in four commercial sugarcane hybrids and one high-fibre biotype for the detection of single nucleotide polymorphisms (SNPs). In this work we tested two Illumina sequencing platforms, read size (70 vs. 150 bp), different sequencing coverage per individual (medium and high coverage), and single-reads versus paired-end reads. We also explored different variant calling strategies (with and without reference genome) and filtering schemes [combining two minor allele frequencies (MAFs) with three depth of coverage thresholds]. For the discovery of a large number of novel SNPs in sugarcane, we recommend longer size and paired-end reads, medium sequencing coverage per individual and Illumina platform NovaSeq6000 for a cost-effective approach, and filter parameters of lower MAF and higher depth coverages thresholds. Although the de novo analysis retrieved more SNPs, the reference-based method allows downstream characterization of variants. For the two best performing matrices, the number of SNPs per chromosome correlated positively with chromosome length, demonstrating the presence of variants throughout the genome. Multivariate comparisons, with both matrices, showed closer relationships among commercial hybrids than with the high-fibre biotype. Functional analysis of the SNPs demonstrated that more than half of them landed within regulatory regions, whereas the other half affected coding, intergenic and intronic regions. Allelic distances values were lower than 0.07 when analysing two replicated genotypes, confirming the protocol robustness.

Genome-Wide Association Studies in Sunflower: Towards Sclerotinia sclerotiorum and Diaporthe/Phomopsis Resistance Breeding.

Diseases caused by necrotrophic fungi, such as the cosmopolitan Sclerotinia sclerotiorum and the Diaporthe/Phomopsis complex, are among the most destructive diseases of sunflower worldwide. The lack of complete resistance combined with the inefficiency of chemical control makes assisted breeding the best strategy for disease control. In this work, we present an integrated genome-wide association (GWA) study investigating the response of a diverse panel of sunflower inbred lines to both pathogens. Phenotypic data for Sclerotinia head rot (SHR) consisted of five disease descriptors (disease incidence, DI; disease severity, DS; area under the disease progress curve for DI, AUDPCI, and DS, AUDPCS; and incubation period, IP). Two disease descriptors (DI and DS) were evaluated for two manifestations of Diaporthe/Phomopsis: Phomopsis stem canker (PSC) and Phomopsis head rot (PHR). In addition, a principal component (PC) analysis was used to derive transformed phenotypes as inputs to a univariate GWA (PC-GWA). Genotypic data comprised a panel of 4269 single nucleotide polymorphisms (SNP), generated via genotyping-by-sequencing. The GWA analysis revealed 24 unique marker-trait associations for SHR, 19 unique marker-trait associations for Diaporthe/Phomopsis diseases, and 7 markers associated with PC1 and PC2. No common markers were found for the response to the two pathogens. Nevertheless, epistatic interactions were identified between markers significantly associated with the response to S. sclerotiorum and Diaporthe/Phomopsis. This suggests that, while the main determinants of resistance may differ for the two pathogens, there could be an underlying common genetic basis. The exploration of regions physically close to the associated markers yielded 364 genes, of which 19 were predicted as putative disease resistance genes. This work presents the first simultaneous evaluation of two manifestations of Diaporthe/Phomopsis in sunflower, and undertakes a comprehensive GWA study by integrating PSC, PHR, and SHR data. The multiple regions identified, and their exploration to identify candidate genes, contribute not only to the understanding of the genetic basis of resistance, but also to the development of tools for assisted breeding.

Morphological and genetic diversity of maize landraces along an altitudinal gradient in the Southern Andes.

Maize (Zea mays ssp. mays) is a major cereal crop worldwide and is traditionally or commercially cultivated almost all over the Americas. The North-Western Argentina (NWA) region constitutes one of the main diversity hotspots of the Southern Andes, with contrasting landscapes and a large number of landraces. Despite the extensive collections performed by the "Banco Activo de Germoplasma INTA Pergamino, Argentina" (BAP), most of them have not been characterized yet. Here we report the morphological and molecular evaluation of 30 accessions collected from NWA, along an altitudinal gradient between 1120 and 2950 meters above sea level (masl). Assessment of morphological variation in a common garden allowed the discrimination of two groups, which differed mainly in endosperm type and overall plant size. Although the groups retrieved by the molecular analyses were not consistent with morphological clusters, they showed a clear pattern of altitudinal structuring. Affinities among accessions were not in accordance with racial assignments. Overall, our results revealed that there are two maize gene pools co-existing in NWA, probably resulting from various waves of maize introduction in pre-Columbian times as well as from the adoption of modern varieties by local farmers. In conclusion, the NWA maize landraces preserved at the BAP possess high morphological and molecular variability. Our results highlight their potential as a source of diversity for increasing the genetic basis of breeding programs and provide useful information to guide future sampling and conservation efforts.

On-field phenotypic evaluation of sunflower populations for broad-spectrum resistance to Verticillium leaf mottle and wilt.

Sunflower Verticillium Wilt and Leaf Mottle (SVW), caused by Verticillium dahliae (Kleb.; Vd), is a soil-borne disease affecting sunflower worldwide. A single dominant locus, known as V1, was formerly effective in controlling North-American Vd races, whereas races from Argentina, Europe and an emerging race from USA overcome its resistance. This emphasizes the need for identifying broad-spectrum genetic resistance (BSR) sources. Here we characterize two sunflower mapping populations (MPs) for SVW resistance: a biparental MP and the association MP from the National Institute of Agricultural Technology (INTA), under field growing conditions. Nine field-trials (FTs) were conducted in highly infested fields in the most SVW-affected region of Argentina. Several disease descriptors (DDs), including incidence and severity, were scored across four phenological stages. Generalized linear models were fitted according to the nature of each variable, adjusting mean phenotypes for inbred lines across and within FTs. Comparison of these responses allowed the identification of novel BSR sources. Furthermore, we present the first report of SVW resistance heritability, with estimates ranging from 35 to 45% for DDs related to disease incidence and severity, respectively. This study constitutes the largest SVW resistance characterization reported to date in sunflower, identifying valuable genetic resources for BSR-breeding to cope with a pathogen of increasing importance worldwide.

Plastome genomics in South American maize landraces: chloroplast lineages parallel the geographical structuring of nuclear gene pools.

BACKGROUND AND AIMS: The number of plastome sequences has increased exponentially during the last decade. However, there is still little knowledge of the levels and distribution of intraspecific variation. The aims of this study were to estimate plastome diversity within Zea mays and analyse the distribution of haplotypes in connection with the landrace groups previously delimited for South American maize based on nuclear markers. METHODS: We obtained the complete plastomes of 30 South American maize landraces and three teosintes by means of next-generation sequencing (NGS) and used them in combination with data from public repositories. After quality filtering, the curated data were employed to search for single-nucleotide polymorphisms, indels and chloroplast simple sequence repeats. Exact permutational contingency tests were performed to assess associations between plastome and nuclear variation. Network and Bayesian phylogenetic analyses were used to infer evolutionary relationships among haplotypes. KEY RESULTS: Our analyses identified a total of 124 polymorphic plastome loci, with the intergenic regions psbE-rps18, petN-rpoB, trnL_UAG-ndhF and rpoC2-atpI exhibiting the highest marker densities. Although restricted in number, these markers allowed the discrimination of 27 haplotypes in a total of 51 Zea mays individuals. Andean and lowland South American landraces differed significantly in haplotype distribution. However, overall differentiation patterns were not informative with respect to subspecies diversification, as evidenced by the scattered distribution of maize and teosinte plastomes in both the network and Bayesian phylogenetic reconstructions. CONCLUSIONS: Knowledge of intraspecific plastome variation provides the framework for a more comprehensive understanding of evolutionary processes at low taxonomic levels and may become increasingly important for future plant barcoding efforts. Whole-plastome sequencing provided useful variability to contribute to maize phylogeographic studies. The structuring of haplotype diversity in the maize landraces examined here clearly reflects the distinction between the Andean and South American lowland gene pools previously inferred based on nuclear markers.

Exploring sunflower responses to Sclerotinia head rot at early stages of infection using RNA-seq analysis.

Sclerotinia head rot (SHR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating sunflower crop diseases. Despite its worldwide occurrence, the genetic determinants of plant resistance are still largely unknown. Here, we investigated the Sclerotinia-sunflower pathosystem by analysing temporal changes in gene expression in one susceptible and two tolerant inbred lines (IL) inoculated with the pathogen under field conditions. Differential expression analysis showed little overlapping among ILs, suggesting genotype-specific control of cell defense responses possibly related to differences in disease resistance strategies. Functional enrichment assessments yielded a similar pattern. However, all three ILs altered the expression of genes involved in the cellular redox state and cell wall remodeling, in agreement with current knowledge about the initiation of plant immune responses. Remarkably, the over-representation of long non-coding RNAs (lncRNA) was another common feature among ILs. Our findings highlight the diversity of transcriptional responses to SHR within sunflower breeding lines and provide evidence of lncRNAs playing a significant role at early stages of defense.

Genetic Diversity, Population Structure and Linkage Disequilibrium Assessment among International Sunflower Breeding Collections.

Sunflower germplasm collections are valuable resources for broadening the genetic base of commercial hybrids and ameliorate the risk of climate events. Nowadays, the most studied worldwide sunflower pre-breeding collections belong to INTA (Argentina), INRA (France), and USDA-UBC (United States of America-Canada). In this work, we assess the amount and distribution of genetic diversity (GD) available within and between these collections to estimate the distribution pattern of global diversity. A mixed genotyping strategy was implemented, by combining proprietary genotyping-by-sequencing data with public whole-genome-sequencing data, to generate an integrative 11,834-common single nucleotide polymorphism matrix including the three breeding collections. In general, the GD estimates obtained were moderate. An analysis of molecular variance provided evidence of population structure between breeding collections. However, the optimal number of subpopulations, studied via discriminant analysis of principal components (K = 12), the bayesian STRUCTURE algorithm (K = 6) and distance-based methods (K = 9) remains unclear, since no single unifying characteristic is apparent for any of the inferred groups. Different overall patterns of linkage disequilibrium (LD) were observed across chromosomes, with Chr10, Chr17, Chr5, and Chr2 showing the highest LD. This work represents the largest and most comprehensive inter-breeding collection analysis of genomic diversity for cultivated sunflower conducted to date.

Exploring gene networks in two sunflower lines with contrasting leaf senescence phenotype using a system biology approach.

BACKGROUND: Leaf senescence is a complex process, controlled by multiple genetic and environmental variables. In sunflower, leaf senescence is triggered abruptly following anthesis thereby limiting the capacity of plants to keep their green leaf area during grain filling, which subsequently has a strong impact on crop yield. Recently, we performed a selection of contrasting sunflower inbred lines for the progress of leaf senescence through a physiological, cytological and molecular approach. Here we present a large scale transcriptomic analysis using RNA-seq and its integration with metabolic profiles for two contrasting sunflower inbred lines, R453 and B481-6 (early and delayed senescence respectively), with the aim of identifying metabolic pathways associated to leaf senescence. RESULTS: Gene expression profiles revealed a higher number of differentially expressed genes, as well as, higher expression levels in R453, providing evidence for early activation of the senescence program in this line. Metabolic pathways associated with sugars and nutrient recycling were differentially regulated between the lines. Additionally, we identified transcription factors acting as hubs in the co-expression networks; some previously reported as senescence-associated genes in model species but many are novel candidate genes. CONCLUSIONS: Understanding the onset and the progress of the senescence process in crops and the identification of these new candidate genes will likely prove highly useful for different management strategies to mitigate the impact of senescence on crop yield. Functional characterization of candidate genes will help to develop molecular tools for biotechnological applications in breeding crop yield.

Transcript and metabolic adjustments triggered by drought in Ilex paraguariensis leaves.

MAIN CONCLUSION: Abscisic acid is involved in the drought response of Ilex paraguariensis. Acclimation includes root growth stimulation, stomatal closure, osmotic adjustment, photoprotection, and regulation of nonstructural carbohydrates and amino acid metabolisms. Ilex paraguariensis (yerba mate) is cultivated in the subtropical region of South America, where the occurrence of drought episodes limit yield. To explore the mechanisms that allow I. paraguariensis to overcome dehydration, we investigated (1) how gene expression varied between water-stressed and non-stressed plants and (2) in what way the modulation of gene expression was linked to physiological status and metabolite composition. A total of 4920 differentially expressed transcripts were obtained through RNA-Seq after water deprivation. Drought induced the expression of several transcripts involved in the ABA-signalling pathway. Stomatal closure and leaf osmotic adjustments were promoted to minimize water loss, and these responses were accompanied by a high transcriptional remodeling of stress perception, signalling and transcriptional regulation, the photoprotective and antioxidant systems, and other stress-responsive genes. Simultaneously, significant changes in metabolite contents were detected. Glutamine, phenylalanine, isomaltose, fucose, and malate levels were shown to be positively correlated with dehydration. Principal component analysis showed differences in the metabolic profiles of control and stressed leaves. These results provide a comprehensive overview of how I. paraguariensis responds to dehydration at transcriptional and metabolomic levels and provide further characterization of the molecular mechanisms associated with drought response in perennial subtropical species.

Snakin-1 affects reactive oxygen species and ascorbic acid levels and hormone balance in potato.

Snakin-1 is a member of the Solanum tuberosum Snakin/GASA family. We previously demonstrated that Snakin-1 is involved in plant defense to pathogens as well as in plant growth and development, but its mechanism of action has not been completely elucidated yet. Here, we showed that leaves of Snakin-1 silenced potato transgenic plants exhibited increased levels of reactive oxygen species and significantly reduced content of ascorbic acid. Furthermore, Snakin-1 silencing enhanced salicylic acid content in accordance with an increased expression of SA-inducible PRs genes. Interestingly, gibberellic acid levels were also enhanced and transcriptome analysis revealed that a large number of genes related to sterol biosynthesis were downregulated in these silenced lines. Moreover, we demonstrated that Snakin-1 directly interacts with StDIM/DWF1, an enzyme involved in plant sterols biosynthesis. Additionally, the analysis of the expression pattern of PStSN1::GUS in potato showed that Snakin-1 is present mainly in young tissues associated with active growth and cell division zones. Our comprehensive analysis of Snakin-1 silenced lines demonstrated for the first time in potato that Snakin-1 plays a role in redox balance and participates in a complex crosstalk among different hormones.

De novo transcriptome sequencing and SSR markers development for Cedrela balansae C.DC., a native tree species of northwest Argentina.

The endangered Cedrela balansae C.DC. (Meliaceae) is a high-value timber species with great potential for forest plantations that inhabits the tropical forests in Northwestern Argentina.Research on this species is scarce because of the limited genetic and genomic information available. Here, we explored the transcriptome of C. balansae using 454 GS FLX Titanium next-generation sequencing (NGS) technology. Following de novo assembling, we identified 27,111 non-redundant unigenes longer than 200 bp, and considered these transcripts for further downstream analysis. The functional annotation was performed searching the 27,111 unigenes against the NR-Protein and the Interproscan databases. This analysis revealed 26,977 genes with homology in at least one of the Database analyzed. Furthermore, 7,774 unigenes in 142 different active biological pathways in C. balansae were identified with the KEGG database. Moreover, after in silico analyses, we detected 2,663 simple sequence repeats (SSRs) markers. A subset of 70 SSRs related to important "stress tolerance" traits based on functional annotation evidence, were selected for wet PCR-validation in C. balansae and other Cedrela species inhabiting in northwest and northeast of Argentina (C. fissilis, C. saltensis and C. angustifolia). Successful transferability was between 77% and 93% and thanks to this study, 32 polymorphic functional SSRs for all analyzed Cedrela species are now available. The gene catalog and molecular markers obtained here represent a starting point for further research, which will assist genetic breeding programs in the Cedrela genus and will contribute to identifying key populations for its preservation.

The first complete genomic structure of Butyrivibrio fibrisolvens and its chromid.

Butyrivibrio fibrisolvens forms part of the gastrointestinal microbiome of ruminants and other mammals, including humans. Indeed, it is one of the most common bacteria found in the rumen and plays an important role in ruminal fermentation of polysaccharides, yet, to date, there is no closed reference genome published for this species in any ruminant animal. We successfully assembled the nearly complete genome sequence of B. fibrisolvens strain INBov1 isolated from cow rumen using Illumina paired-end reads, 454 Roche single-end and mate pair sequencing technology. Additionally, we constructed an optical restriction map of this strain to aid in scaffold ordering and positioning, and completed the first genomic structure of this species. Moreover, we identified and assembled the first chromid of this species (pINBov266). The INBov1 genome encodes a large set of genes involved in the cellulolytic process but lacks key genes. This seems to indicate that B. fibrisolvens plays an important role in ruminal cellulolytic processes, but does not have autonomous cellulolytic capacity. When searching for genes involved in the biohydrogenation of unsaturated fatty acids, no linoleate isomerase gene was found in this strain. INBov1 does encode oleate hydratase genes known to participate in the hydrogenation of oleic acids. Furthermore, INBov1 contains an enolase gene, which has been recently determined to participate in the synthesis of conjugated linoleic acids. This work confirms the presence of a novel chromid in B. fibrisolvens and provides a new potential reference genome sequence for this species, providing new insight into its role in biohydrogenation and carbohydrate degradation.

Main and epistatic QTL analyses for Sclerotinia Head Rot resistance in sunflower.

Sclerotinia Head Rot (SHR), a disease caused by Sclerotinia sclerotiorum, is one of the most limiting factors in sunflower production. In this study, we identified genomic loci associated with resistance to SHR to support the development of assisted breeding strategies. We genotyped 114 Recombinant Inbred Lines (RILs) along with their parental lines (PAC2 -partially resistant-and RHA266 -susceptible-) by using a 384 single nucleotide polymorphism (SNP) Illumina Oligo Pool Assay to saturate a sunflower genetic map. Subsequently, we tested these lines for SHR resistance using assisted inoculations with S. sclerotiorum ascospores. We also conducted a randomized complete-block assays with three replicates to visually score disease incidence (DI), disease severity (DS), disease intensity (DInt) and incubation period (IP) through four field trials (2010-2014). We finally assessed main effect quantitative trait loci (M-QTLs) and epistatic QTLs (E-QTLs) by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively. As a result of this study, the improved map incorporates 61 new SNPs over candidate genes. We detected a broad range of narrow sense heritability (h2) values (1.86-59.9%) as well as 36 M-QTLs and 13 E-QTLs along 14 linkage groups (LGs). On LG1, LG10, and LG15, we repeatedly detected QTLs across field trials; which emphasizes their putative effectiveness against SHR. In all selected variables, most of the identified QTLs showed high determination coefficients, associated with moderate to high heritability values. Using markers shared with previous Sclerotinia resistance studies, we compared the QTL locations in LG1, LG2, LG8, LG10, LG11, LG15 and LG16. This study constitutes the largest report of QTLs for SHR resistance in sunflower. Further studies focusing on the regions in LG1, LG10, and LG15 harboring the detected QTLs are necessary to identify causal alleles and contribute to unraveling the complex genetic basis governing the resistance.

Integration of transcriptomic and metabolic data reveals hub transcription factors involved in drought stress response in sunflower (Helianthus annuus L.).

KEY MESSAGE: By integration of transcriptional and metabolic profiles we identified pathways and hubs transcription factors regulated during drought conditions in sunflower, useful for applications in molecular and/or biotechnological breeding. Drought is one of the most important environmental stresses that effects crop productivity in many agricultural regions. Sunflower is tolerant to drought conditions but the mechanisms involved in this tolerance remain unclear at the molecular level. The aim of this study was to characterize and integrate transcriptional and metabolic pathways related to drought stress in sunflower plants, by using a system biology approach. Our results showed a delay in plant senescence with an increase in the expression level of photosynthesis related genes as well as higher levels of sugars, osmoprotectant amino acids and ionic nutrients under drought conditions. In addition, we identified transcription factors that were upregulated during drought conditions and that may act as hubs in the transcriptional network. Many of these transcription factors belong to families implicated in the drought response in model species. The integration of transcriptomic and metabolomic data in this study, together with physiological measurements, has improved our understanding of the biological responses during droughts and contributes to elucidate the molecular mechanisms involved under this environmental condition. These findings will provide useful biotechnological tools to improve stress tolerance while maintaining crop yield under restricted water availability.

ATGC transcriptomics: a web-based application to integrate, explore and analyze de novo transcriptomic data.

BACKGROUND: In the last years, applications based on massively parallelized RNA sequencing (RNA-seq) have become valuable approaches for studying non-model species, e.g., without a fully sequenced genome. RNA-seq is a useful tool for detecting novel transcripts and genetic variations and for evaluating differential gene expression by digital measurements. The large and complex datasets resulting from functional genomic experiments represent a challenge in data processing, management, and analysis. This problem is especially significant for small research groups working with non-model species. RESULTS: We developed a web-based application, called ATGC transcriptomics, with a flexible and adaptable interface that allows users to work with new generation sequencing (NGS) transcriptomic analysis results using an ontology-driven database. This new application simplifies data exploration, visualization, and integration for a better comprehension of the results. CONCLUSIONS: ATGC transcriptomics provides access to non-expert computer users and small research groups to a scalable storage option and simple data integration, including database administration and management. The software is freely available under the terms of GNU public license at http://atgcinta.sourceforge.net .

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

BACKGROUND: In 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. RESULTS: Weighted 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. CONCLUSIONS: The 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.

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.