SciELO - Scientific Electronic Library Online

 
vol.53 issue2Tropicalization of canola (Brassica napus L.): commercial hybrids show potential for cultivation in the Brazilian CerradoGermination and development of M1 seedlings of two Selliera radicans Cav. accessions subjected to gamma radiation author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

  • Have no cited articlesCited by SciELO

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista de la Facultad de Ciencias Agrarias. Universidad Nacional de Cuyo

On-line version ISSN 1853-8665

Rev. Fac. Cienc. Agrar., Univ. Nac. Cuyo vol.53 no.2 Mendoza Dec. 2021

 

Original article

Prioritization of vigor QTL-associated genes for future genome-directed Vitis breeding

Priorización de genes asociados a QTLs de vigor para futuros planes de mejoramiento dirigido en Vitis

Inés Hugalde1  2  * 

Marcos Paolinelli1  3  * 

Cecilia B. Agüero2 

Summaira Riaz2 

Sebastián Gómez Talquenca1 

M. Andrew Walker2 

Hernán Vila1 

1 Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agrope cuaria Mendoza. Luján de Cuyo (5507). Mendoza. Argentina. hugalde.ines@inta.gob.ar

2 University of California Davis. Department of Viticulture and Enology. Davis (CA 95616). California. USA.

3 Centro Científico Tecnológico-Consejo Nacional de Investigaciones Científicas y Técnicas de Argentina (CCT-CONICET). Mendoza (5500). Mendoza. Argentina.

Abstract

Vigor control in grapevine may become especially important under climate change. A better understanding of gene-phenotype relationships is required in order to exploit plant genomics for breeding purposes. This research aims to use quantitative trait loci (QTLs) for vigor identified in the progeny from a cross of Ramsey (Vitis champinii) × Riparia Gloire (V. riparia). Genes located 700 kb up and downstream from each QTL position were interro gated for functional enrichment through ShinyGO online tool, based on the gene ontology annotation of Vitis vinifera PN40024. Key biological processes like phloem and xylem development, cell cycle, response to hormones, amino acid transport, tissue development, sugar metabolism, nitrogen transport, and stress/immune responses, showed functional enrichment. Integral response to light and auxin might be required for fine molecular tuning of vegetative growth in Vitis. Fifty out of 1318 candidate genes were prioritized, reducing their amount to a manageable number of candidates genes for further directed breeding strategies.

Keywords: Biological processes; Gene ontology; QTL map; Vigor; Gene prioritization

Resumen

El control del vigor en vid es un factor de gran importancia en el contexto de cambio climático actual. Es necesario desarrollar una mejor comprensión de la relación gen-fe notipo para la asistencia del mejoramiento vegetal mediante genómica. En este trabajo se utilizó la información obtenida de un mapeo de QTLs para vigor realizado con la progenie de Ramsey (Vitis champinii) × Riparia Gloire (V. riparia). De acuerdo a la anotación onto lógica de Vitis vinifera PN40024, los genes ubicados 700 kb por encima y por debajo de cada marcador fueron interrogados para determinar enriquecimiento funcional a través de la herramienta online ShinyGO. Distintos procesos clave en la definición de vigor mostraron enriquecimientos altamente significativos, tales como desarrollo tisular (floema y xilema), ciclo celular, transporte de aminoácidos (nitrógeno) metabolismo de azúcares y respuesta a hormonas, a estrés e inmune. La respuesta integral a la luz y la producción de auxinas parece modular molecularmente el crecimiento vegetativo en Vitis.

Palabras clave: Proceso biológico; Ontología génica; Mapeo de QTL; Vigor; Priorización de genes

“We respectfully remember Dr. Hernán Vila, who lost his battle against COVID-19 in October 2020 and dedicate this article to his memory”

Introduction

The understanding of genotype-phenotype associations constitute a major challenge for plant scientists. The causal association between variation of a certain trait and genotypic differences is the foundation for developing targeted strategies to be used in molecular breeding 1.

Vigor is considered as the propensity to assimilate, store, and/or use non-structural carbohydrates to produce large canopies. It is associated with high metabolism and rapid shoot growth 17. The main processes that account for vigor are: 1- carbon assimilation or photosynthesis, 2- cell growth, by division and expansion, and 3- leaf area development, for the reception of CO2 and light. All these complex processes contribute to what is known as a highly quantitative trait.

There is increasing interest in understanding the genetic basis of grapevine vigor and biomass production, given its substantial impact on yield, water management, plant health, and fruit quality. Advances in understanding the genetic bases of complex traits through genetic mapping and quantitative trait locus (QTL) analysis have linked certain complex phenotypes to specific regions of chromosomes, and helped to identify the action, number, and precise location of these regions. The Ramsey (Vitis champinii) × Riparia Gloire (V. riparia) linkage map was developed by Lowe and Walker (2006) with the goal of mapping traits like biotic resistance, drought tolerance and vigor. We recently reported QTLs for vigor in this progeny that are associated with leaf area, specific leaf area, field canopy biomass and several partitioning indices 12. However, the genes behind these quantitative traits have not been identified. Here, we analyze the function of genes surrounding the QTLs based on gene location and the Gene Ontology (GO) annotation of the V. vinifera genome. GO is a computational representation of the functions of protein and non-coding RNA mole cules produced by genes from many different organisms. It can be used for interpreting large-scale molecular biology experiments to gain insight into the structure, function, and dynamics of an organism. This approach was proposed by Correa et al. (2014) when iden tifying QTLs and candidate genes associated with cluster architecture in grapevines, and later by Ritcher et al. (2019), when approaching cluster architecture and revealing genes enriched by physical projection of the QTLs onto the PN40024 reference genome sequence.

In this paper, we communicate the results of a GO functional enrichment to prioritize candidate genes responsible for vigor that interrogated 1318 genes obtained from the QTL mapping of the Ramsey x Riparia Gloire progeny. We also attempt to elucidate the putative role of genes in terms of physiological and biochemical processes and identify priority genes for further targeted breeding approaches.

Materials and methods

Six highly significant QTLs for vigor components from the Ramsey x Riparia Gloire map were selected based on their LOD values above the genome wide significance level 12. The BLAST server against the V. vinifera PN40024 genome (https://urgi.versailles.inra.fr/blast/) was then used to locate the hybridization regions of PCR primers targeting the simple sequence repeat (SSR) markers linked to the selected QTLs (Table 1, page 30-31).

Table 1: Prioritized genes for vigor based on Gene Ontology enrichment analysis and hypothetical function of predicted proteins. SSR markers associated with selected QTLs are shown with chromosome number. Tabla 1: Genes priorizados para vigor basados en enriquecimiento por ontología génica y sus respectivas funciones proteicas hipotéticas. Los marcadores SSR asociados con los QTLs elegidos se identifican por número de cromosoma. 

The coordinates 700 kb upstream and downstream of primer hybridization on the reference sequence were used to identify the genes within the region flanked by those limits. The 700 kb limit was established considering that candidate genes should be no farther than 3-4 cM from the marker, since the V. vinifera genome has been estimated at about 475 Mb and 1 cM equals 200-300 kb 3. For that purpose, a script was developed to search for annotated genes flanking the QTL regions in the annotation file of the V. vinifera PN40024 genome.

The annotated genes found with the script were then used to interrogate for functional enrichment through ShinyGO online tool version 0.61 8. For the set of genes contained in the QTL regions associated with vigor, the occurrence of associated gene ontology terms (GO terms) was statistically evaluated for overrepresentation through an hypergeometric test retrieving those GOterms with statistical significance (False Discovery Rate < 0.05) when comparing their percentage of occurrence with the percentage of each GOterm in the whole annotated genome. Integral plots that associate gene locations in chromosomes with enriched GO categories were performed with shinyCircos 22.

Results

In this work, GO functional enrichment allowed the identification of genes involved in key processes related to vigor and prioritization reduced the number of candidate genes from 1318 to 50.

Our results show that associations between overrepresented GO terms and vigor helped to rank candidate genes, based on their putative function. Phloem and xylem development, cell cycle, response to hormones, tissue development, aminoacid and nitrogen transport, sugar metabolism and immune responses, all showed functional enrichment (Figure 1 page 32; Table 1, page 30-31).

Lines connect genes to enriched GO terms (ranging from FDR 0.034 in black, to FDR 1x10-9 in gold).

Las líneas conetan los genes con los términos de enriquecimiento desde FDR 0,034 en negro hasta FDR 1x10-9 en dorado).

Figure 1: Figura 1: Chromosome location of vigor-associated genes identified through GO functional enrichment. Genes and SSR markers showing GO categories enrichments (FDR<0.05) are indicated on chromosomes of V vinifera. Ubicación cromosómica de genes asociados a vigor identificados por enriquecimiento funcional GO (FDR<0,05) y marcadores SSR en los cromosomas de V vinifera. 

On chromosome 1, most identified genes encode for aminoacid transmembrane transporters (Figure 1 page 32; Table 1, page 30-31). In addition, the single gene found on chromosome 3 influences the enrichment for functions related to transport of organic acids and nitrogen compounds.

Two important transcription factors (TFs) related to photomorphogenesis (TF 104879018) and auxins (TF 104879021) were identified on chromosome 4 (Figure 1 page 32; Table 1, page 30-31). TFs are especially interesting because they control the expression of several genes. Genes in chromosome 10 are involved in biotic and abiotic stress responses (Figure 1 page 32; Table 1, page 30-31). Another transcription factor related to this processes, identified in Table 1 (page 30-31) as 100243518, appears significant in this chro mosome. On chromosomes 14 and 19, we found the most varied group of genes in terms of function (Figure 1 page 32; Table 1, page 30-31), including, once again, genes encoding for nitrogen transport, that are key for growth, cell cycle regulation and development 13,19. Many genes encode for proteins involved in phloem and xylem development. This was the function with major enrichment (Figure 1 page 32; Table 1, page 30-31). It represents a key aspect for growth, as transport of water, nutrients, proteins and sugars is vital for the plant to develop.

Discussion

Although the analysis used the V. vinifera PN40024 genome, due to its superior charac terization and annotation, sequence alignment with V. riparia Gloire 9 of chromosomes with significant QTLs averaged 96,77 % (data not shown). Also, Liang et al. (2019) deter mined a 97.34-97.65% identity through whole genome comparison of V. vinifera PN40024 with two V. riparia accessions.

As mentioned, on chromosome 1, most identified genes encode for aminoacid transmem brane transporters (Figure 1; Table 1, page 30-31). Transmembrane aminoacid transport is associated with enhanced growth and high rates of protein synthesis 20. Nitrogen can be taken up from soil in various forms, being one of them amino acids, and is of considerable importance in vigor control, yield and berry quality 7. The fact that these genes are closely located, suggests some kind of structural regulation. It has been observed that root elon gation and enlargement in the rootstock 110R, partly depend on transcriptomic regulations of sugar and protein transporter genes SWEET and NRT1/PTR in roots. This was found to facilitate carbohydrate and nitrogen accumulation, providing essential energy to 110R roots under drought 21.

In relation to the two important TFs related to photomorphogenesis and auxins identified on chromosome 4, they are particularly intriguing given that they control the expression of several genes. The particularity of having both TFs in the same chromosomal region increases the probability of a synergistic response through coregulation. Indeed, the coor dination between TFs in response to light and auxins was well established in Arabidopsis thaliana by Halliday, et al. (2009), who showed that light regulates growth of distant tissues from the site of light exposure through auxin production. Something similar was observed in A. thaliana by Hornitschek, et al. (2012), where phytochrome interacting factors 4 and 5 (PIF4 and PIF5) regulated elongation growth by controlling the expression of genes that encode for auxin biosynthesis and signaling. Interestingly, TF 104879018 in our study is homologous to HY5 of A. thaliana, which is found downstream in the signaling cascade of PIF1/PIF310. HY5 promotes growth, especially through photosynthesis induction and higher nutrient uptake by roots. The other prioritized TF 104879021, is a homolog to the auxin-responsive protein IAA28 of A. thaliana, which plays a role in regulation of lateral root growth. In grapevines, kaolin, a mineral that reflects radiation from the leaf surface, produced an increase in IAA content. This treatment also caused higher values of stomatal conductance, net CO2 assimilation rate, intrinsic water use efficiency, and a slight decrease in ABA 5. These results might be supported by the same mechanism connecting growth hormones and light interception. Therefore, vigor in grapevines may partly depend on promotion of photosynthesis, lateral root development and nitrate uptake, and these processes may be associated through the expression of genes downstream-regulated by TF 104879018 and TF 104879021.

Regarding the genes and TF found in chromosome 10, involved in biotic and abiotic stress responses, some authors have shown significant correlation between vigor and tolerance/ susceptibility to diseases that could induce different defense responses in the host plant 2. Vigorous plants may have developed stronger immune responses to defend themselves.

Considering our findings in chromosomes 14 and 19, many genes encode for proteins involved in phloem and xylem development. This function showed the major enrichment (Figure 1 page 32; Table 1, page 30-31), representing a key aspect for growth, as transport of water, nutrients, proteins and sugars is vital for the plant to develop. These functions are tightly correlated to auxins and soluble carbohydrates seasonal dynamics, since cambium activity and xylem/phloem development respond to this signaling in woody species. It has been observed that IAA and soluble carbohydrate dynamics directly affect xylem and phloem formation in trees 6. In addition, gibberellins increase carbon allocation to different organs by inducing accumulation of non-structural carbohydrates in leaves, enhancement of phloem area and expression of sugar transporters 16.

Our results may lead to deeper gene selection strategies, aiming at choosing a smaller number of genes. Candidates with the smallest enrichment FDR values that are associated with QTLs that explain the highest percentages of variation, constitute interesting targets. For quantitative characters, where positive feedbacks cause large effects, significant explanatory effects from 10% to 20%, may result in impressive effects. This strategy considers that both approaches, QTL mapping and GO enrichment, work at different levels. While QTLs indentify regions on chromosomes containing genes encoding for a certain trait, the enrichment process takes into account all the genes involved in a pathway related to that trait. In processes where numerous genes are involved, a meaningful change should include many of them. Conse quently, transcription factors could produce phenotype differences even at higher FDR values. In our work, this last effort selected 16 genes and 4 TRFs as plausible candidates for further breeding studies (Table 1, page 30-31). Further functional genomic studies should weigh the importance of these selected genes on the final phenotype.

Conclusions

In conclusion, this analysis allowed the detection of plausible candidate genes encoding for the components of key processes governing vegetative growth in Vitis. The analysis allowed the reduction of candidate gene number based on marker proximity and functional enrichment, clearly demonstrating a suitable shortcut for target-directed genome-guided breeding strategies. This approach is particularly useful when a map is not densely saturated.

Two TFs, which potentially enhance growth by relating light response to hormone acti vation, and then to photosynthesis and morphogenesis, are strong candidates for targeted breeding. The nitrogen transport encoding genes would allow this light/hormone promoted growth by facilitating amino acid/protein synthesis and transport. Phloem and xylem related genes would also be part of this process by enabling water and nutrient transport. All these functions need to be tightly correlated, since auxins and soluble carbohydrates seasonal dynamics play key roles in tissue growth, cambium activity and xylem/phloem development. Gene characterization in individuals of the Ramsey × Riparia Gloire progeny will be the topic for future research.

References

1. Bargsten, J. W.; Nap, J. P.; Sanchez-Perez, G. F.; Aalt, D. J. van Dijk. 2014. Prioritization of candidate genes in QTL regions based on associations between traits and biological processes. BMC Plant Biol 14: 330 https://doi.org/10.1186/s12870-014-0330-3 [ Links ]

2. Calonnec A.; Burie, J. B.; Langlais, M.; Guyader, S.; Saint-Jean, S.; Sache, I.; Tivoli, B. 2013. Impacts of plant growth and architecture on pathogen processes and their consequences for epidemic behaviour. Eur J Plant Pathol 135: 479-497. doi: 10.1007/s10658-012-0111-5 [ Links ]

3. Cipriani, G.; Di Gaspero, G.; Canaguier, A.; Jusseaume, J.; Tassin, J.; Lemainque, A.; Thareau, V.; Adam-Blondon, A. F.; Testolin, R. 2011. Molecular linkage maps: strategies, resources and achievements. In: Adam-Blondon, AF; Martinez Zapater, JM; Kole, C; (eds). Genetics, genomics and breeding of grapes. Science Publishers and CRC Press. 111-136. doi: 10.1201/b10948 [ Links ]

4. Correa, J.; Mamani, M.; Muñoz-Espinoza, C.; Laborie, D.; Muñoz, C.; Pinto, M.; Hinrichsen, P. 2014. Heritability and identification of QTLs and underlying candidate genes associated with the architecture of the grapevine cluster (Vitis vinifera L.). Theoretical and applied genetics. 127(5): 1143-1162. doi: https://doi.org/10.1007/s00122-014-2286-y [ Links ]

5. Dinis, L. T.; Bernardo, S.; Luzio, A.; Pinto, G.; Meijón, M.; Pintó-Marijuan, M.; Cotado, A.; Correia, C.; Moutinho Pereira, J. 2018. Kaolin modulates ABA and IAA dynamics and physiology of grapevine under Mediterranean summer stress. J Plant Physiol. 220: 181-192. doi: 10.1016/j.jplph.2017.11.007 [ Links ]

6. Fajstavr, M.; Paschová, Z.; Giagli, K.; Vavrčík, H.; Gryc, V.; Urban, J. 2018. Auxin (IAA) and soluble carbohydrate seasonal dynamics monitored during xylogenesis and phloemogenesis in Scots pine. Forest. 11: 553-562. doi: 10.3832/ifor2734-011 [ Links ]

7. Gatti, M.; Squeri, C.; Garavani, A.; Frioni, T.; Dosso, P.; Diti, I.; Poni, S. 2019. Effects of variable rate nitrogen application on cv Barbera performance: yield and grape composition. Am J Enol Vitic. 70:188-200. doi: 10.5344/ajev.2019.18072 [ Links ]

8. Ge, S. X.; Jung, D.; Yao, R. 2020. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics. 36: 2628-2629. doi: 10.1093/bioinformatics/btz931 [ Links ]

9. Girollet, N.; Rubio, B.; López-Roques, C.; Valière, S.; Ollat, N.; Bert, P. F. 2019. De novo phased assembly of the Vitis riparia grape genome. Sci Data. 6: 127. doi: 10.1038/s41597-019-0133-3 [ Links ]

10. Halliday, K. J.; Martínez-García, J. F.; Josse, E. M. 2009. Integration of light and auxin signaling. Cold Spring Harb Perspect Biol. 1:a001586. doi: 10.1101/cshperspect.a001586 [ Links ]

11. Hornitschek, P.; Kohnen, M. V.; Lorrain, S.; Rougemont, J.; Ljung, K.; López-Vidriero, I.; Franco-Zorrilla, J. M.; Solano, R.; Trevisan, M.; Pradervand, S.; Xenarios, I.; Fankhauser, C. 2012. Phytochrome interacting factors 4 and 5 control seedling growth in changing light conditions by directly controlling auxin signaling. Plant J. 71: 699-711. doi: 10.1111/j.1365-313X.2012.05033.x [ Links ]

12. Hugalde, I. P.; Riaz, S.; Agüero, C. B.; Vila, H.; Gomez Talquenca, S.; Walker, M. A. 2019. Studying growth and vigor as quantitative traits in grapevine populations. In: Maia, RT; de Araújo Campos, M; (eds). Integrated View of Population Genetics. IntechOPen. doi: 10.5772/intechopen.82537 [ Links ]

13. Lang, C. P.; Merkt, N.; Zörb, C. 2018. Different nitrogen (N) forms affect responses to N form and N supply of rootstocks and grafted grapevines. Plant Sci. 277: 311-321. doi: 10.1016/j. plantsci.2018.10.004 [ Links ]

14. Liang, Z.; Duan, S.; Sheng, J.; Zhu, S.; Ni, X.; Shao, J.; Liu, C.; Nick, P.; Du, F.; Fan, P.; Mao, R.; Zhu, Y.; Deng, W.; Yang, M.; Huang, H.; Liu, Y.; Ding, Y.; Liu, X.; Jiang, J.; Zhu, Y.; Li, S.; He, X.; Chen, W.; Dong, Y. 2019. Wholegenome resequencing of 472 Vitis accessions for grapevine diversity and demographic history analyses. Nat Commun. 10: 1190. doi: 10.1038/s41467-019-09135-8 [ Links ]

15. Lowe, K. M.; Walker, M. A. 2006. Genetic linkage map of the interspecific grape rootstock cross Ramsey (Vitis champinii) × Riparia Gloire (Vitis riparia). Theor Appl Genet. 112: 1582- 1592. doi: 10.1007/s00122-006-0264-8 [ Links ]

16. Murcia, G.; Pontin, M.; Reinoso, H.; Baraldi, R.; Bertazza, G.; Gómez-Talquenca, S.; Bottini, R.; Piccoli, P. N. 2016. ABA and GA3 increase carbon allocation in different organs of grapevine plants by inducing accumulation of non-structural carbohydrates in leaves, enhancement of phloem area and expression of sugar transporters. Physiol Plant. 156: 323-337. doi: 10.1111/ppl.12390 [ Links ]

17. Rebolledo, M. C.; Dingkuhn, M.; Courtois, B.; Gibon, Y.; Clément-Vidal, A.; Cruz, D. F.; Duitama, J.; Lorieux, M.; Luquet, D. 2015. Phenotypic and genetic dissection of component traits for early vigour in rice using plant growth modelling, sugar content analyses and association mapping. J Exp Bot. 66: 5555-5566. doi: 10.1093/jxb/erv258 [ Links ]

18. Richter, R.; Gabriel, D.; Rist, F.; Töpfer, R.; Zyprian, E. 2019. Identification of co-located QTLs and genomic regions affecting grapevine cluster architecture. Theor Appl Genet. 132: 1159-1177. https://doi. org/10.1007/s00122-018-3269-1 [ Links ]

19. Tegeder, M.; Masclaux-Daubresse, C. 2018. Source and sink mechanisms of nitrogen transport and use. New Phytol. 217: 35-53. doi: 10.1111/nph.14876 [ Links ]

20. Wipf, D.; Ludewig, U.; Tegeder, M.; Rentsch, D.; Koch, W.; Frommer, W. B. 2002. Conservation of amino acid transporters in fungi, plants and animals. Trends Biochem Sci. 27: 139-147. doi: 10.1016/ S0968-0004(01)02054-0 [ Links ]

21. Yıldırım, K.; Yağcı, A.; Sucu, S.; Tunç, S. 2018. Responses of grapevine rootstocks to drought through altered root system architecture and root transcriptomic regulations. Plant Physiol Bioch. 127: 256-268. doi: 10.1016/j.plaphy.2018.03.034 [ Links ]

22. Yu, Y.; Ouyang, Y.; Yao, W. 2018. ShinyCircos: an R/Shiny application for interactive creation of circos plot. Bioinformatics. 34: 1229-1231. doi: 10.1093/bioinformatics/btx763. [ Links ]

Received: May 11, 2021; Accepted: October 12, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License