- Research article
- Open Access
Association mapping in sunflower (Helianthus annuus L.) reveals independent control of apical vs. basal branching
BMC Plant Biologyvolume 15, Article number: 84 (2015)
Shoot branching is an important determinant of plant architecture and influences various aspects of growth and development. Selection on branching has also played an important role in the domestication of crop plants, including sunflower (Helianthus annuus L.). Here, we describe an investigation of the genetic basis of variation in branching in sunflower via association mapping in a diverse collection of cultivated sunflower lines.
Detailed phenotypic analyses revealed extensive variation in the extent and type of branching within the focal population. After correcting for population structure and kinship, association analyses were performed using a genome-wide collection of SNPs to identify genomic regions that influence a variety of branching-related traits. This work resulted in the identification of multiple previously unidentified genomic regions that contribute to variation in branching. Genomic regions that were associated with apical and mid-apical branching were generally distinct from those associated with basal and mid-basal branching. Homologs of known branching genes from other study systems (i.e., Arabidopsis, rice, pea, and petunia) were also identified from the draft assembly of the sunflower genome and their map positions were compared to those of associations identified herein. Numerous candidate branching genes were found to map in close proximity to significant branching associations.
In sunflower, variation in branching is genetically complex and overall branching patterns (i.e., apical vs. basal) were found to be influenced by distinct genomic regions. Moreover, numerous candidate branching genes mapped in close proximity to significant branching associations. Although the sunflower genome exhibits localized islands of elevated linkage disequilibrium (LD), these non-random associations are known to decay rapidly elsewhere. The subset of candidate genes that co-localized with significant associations in regions of low LD represents the most promising target for future functional analyses.
Shoot branching is a major determinant of plant architecture and plays an important role in the adaptation of plants to their environment. Variation in branching helps plants compete with their neighbors and also offers protection against herbivory [1-4]. Shoot branching can also affect developmental phenotypes such as flowering time and reproductive success . Moreover, this trait is an important component of the so-called “domestication syndrome” , with many crops exhibiting reduced branching (i.e., increased apical dominance) relative to their wild progenitors.
In cultivated sunflower (Helianthus annuus L.), selection during domestication resulted in the production of an apically dominant, unbranched growth form that differs markedly from its highly branched wild progenitor (common sunflower; also H. annuus) [7-10]. During the transition to hybrid breeding in the mid-20th century, however, recessive branching was reintroduced to the sunflower gene pool to produce male-fertile restorer (R) lines that can be crossed with unbranched, cytoplasmic male-sterile (i.e., female; A) lines to produce unbranched, fully fertile hybrids. The apical branching of the R-lines is desirable because it provides fertile pollen for a longer time period, resulting in a longer window for pollination and hybrid production . The modern-day cultivated sunflower gene pool thus exhibits substantial variation in plant architecture, making it an ideal system to study the genetics of branching.
In general terms, branching is initiated from axillary meristems in leaf axils on the primary shoot. These meristems give rise to axillary buds which remain dormant or grow out into a branch that can be influenced by environmental conditions or developmental signals such as hormones . Three phytohormones (auxin, cytokinin [CK], and strigolactone [SL]) and genes associated with their homeostasis and signaling are thought to be largely responsible for the regulation of branching [12-15]. Bud outgrowth is inhibited by basipetal transport of auxin produced at the shoot apical meristem. SL is a carotenoid-derived phytohormone that also inhibits bud outgrowth. It is produced in the roots and transported acropetally in the stem [12,16]. In contrast to auxin and SL, CKs are locally synthesized in the bud and promote the outgrowth of the axillary bud. Ultimately, cross talk by these phytohormone related pathways regulates branching . Additionally, genes related to gibberellic acid (GA) and polyamine metabolism, and genes encoding transcription factors, at least one MAP kinase, and cytochrome P450 all play important roles in axillary bud initiation and branch growth [15-18].
In crosses between cultivated and wild sunflower, branching has been found to be a genetically complex trait influenced by numerous small effect loci distributed throughout the genome [19,20]. Classical genetic analyses in cultivated sunflower have, however, revealed the existence of loci with major effects on both apical and basal branching [11,21-23]. More recently, quantitative trait locus (QTL) mapping has been used to localize the recessive apical branching of restorer lines to a region (known as the B locus) on the upper half of linkage group (LG) 10 [24,25]. The unbranched phenotype characteristic of female lines and hybrids is thought to be controlled by the dominant B allele, while the branched R-lines are homozygous for the recessive b allele. While traditional QTL analyses have proven to be useful in identifying genomic regions that influence plant architecture in sunflower, this general approach suffers some limitations. Most notably, the use of biparental populations only enables the analysis of two alleles per gene, and also limits the number of recombination events, thereby providing relatively limited genetic resolution .
Association mapping (also known as LD mapping) has emerged as an alternative to QTL mapping for investigating the genetic basis of quantitative traits . Because it involves the analysis of a diverse collection of more or less unrelated individuals, association mapping allows for the simultaneous evaluation of the effects of multiple haplotypes across diverse genetic backgrounds. Moreover, because association populations typically capture numerous generations of historical recombination, this approach provides much higher resolution than is possible with a family-based mapping population. Herein, we report the results of a detailed analysis of variation in branching in an association mapping population that captures nearly 90% of the allelic diversity present within the cultivated sunflower gene pool [28,29]. We evaluated this population for various branching-related traits at three different locations and tested for genetic associations across the genome using genotypic data derived from a high-density SNP array [30,31]. We also identified candidate genes involved in hormonal or transcriptional regulation that mapped in close proximity to significant associations.
Development of the association mapping population
The development and initial characterization of the sunflower association mapping population utilized herein are described by Mandel et al. [28,29]. Briefly, a diverse collection of cultivated sunflower lines was obtained from the USDA North Central Regional Plant Introduction Station (NCRPIS; Ames, IA, USA) and from the French National Institute for Agricultural Research (INRA; France). These lines were genotyped with simple-sequence repeat (SSR) markers distributed across all 17 LGs and the resulting data were used to identify hierarchical subsets of lines that captured maximum diversity . The present study was based on the same subset of 288 lines employed by Mandel et al. , which differed slightly from the core 288 described by Mandel et al.  due to limited seed availability for some lines. This population, which is available for distribution from the Germplasm Resources Information Network (GRIN) of the National Plant Germplasm System, is known as UGA-SAM1. Of the full set of 288 lines in this population, only 271 were included in our final analyses due to germination difficulties and plant loss during the growing season. These lines capture nearly 90% of the allelic diversity and include lines that are oil and confectionery types from the two major heterotic groups in cultivated sunflower as well as select open-pollinated varieties (OPVs) and land races. Many of the lines were advanced via single-seed descent for one or two cycles to reduce residual heterozygosity prior to the start of this experiment.
Field design and phenotypic analysis of branching traits
All 288 lines were planted in two replicates using an alpha lattice at three different locations during the spring of 2010. The three locations were: the UGA Plant Sciences Farm, (Oconee County, GA, USA), the NCRPIS at Iowa State University (Ames, IA, USA), and the University of British Columbia’s Botany Gardens (Vancouver, BC, Canada) . Individuals of each line (3–4 individuals per replicate per location) were evaluated at the R9 reproductive stage, which represents physiological maturity .
At maturity, each plant was divided into four quarters (1st, 2nd, 3rd, and 4th, from top to bottom) and the number of primary branches was counted in each quarter. Primary branches were recorded as present if they were longer than 2 cm and had developed a terminal inflorescence. The numbers of primary branches in each quarter were then used to estimate the extent and type of branching on each plant. Apical branching was estimated as the number of primary branches in the 1st quarter, mid-apical branching was estimated as total number of branches in the 1st and 2nd quarters, and so forth (Figure 1). If a particular quarter did not branch, the leaf axils were examined to determine if axillary bud initiation had occurred. If more than 50% of the nodes in a particular quarter displayed initiation, it was recorded as branch initiation. In addition, the presence or absence of secondary branches was recorded.
For each line, branch numbers and lengths were averaged within replicates prior to analysis. For secondary branching, scores of 0 for absence or 1 for presence were assigned and these values were subsequently averaged, as above. Using PROC GLM in SAS (ver. 9.3; SAS Institute, Cary NC), we found a significant genotype x environment (G x E) interaction (P < 0.001). Therefore, all subsequent analyses were performed separately by location, which also made possible the identification of loci that appear to be susceptible to environmental effects (i.e., those that were significant in at least one, but not all locations). In each location, genotypes were treated as fixed effects and block and replicates as random effects. For association mapping of the various branching traits, we used least-squares means (LS means) since block and rep effects were found to be statistically significant (P < 0.05).
Correlations between various branching traits and among locations were determined using Spearman’s rank correlation coefficients (ρ) in JMP (ver. 9; SAS Institute), and corrected for multiple tests using the sequential Bonferroni correction . Principal component analysis (PCA) was performed to visualize the various branching types (e.g., apical vs. basal branching) within the SAM association mapping population using the pca function implemented in the FactoMineR package ver. 1.16  available in the R statistical computing language (ver. 3.1.0) .
As described by Mandel et al. , total DNA was extracted from pooled leaf tissue from four individuals of each line using a CTAB extraction protocol . These samples were then genotyped using an Infinium SNP array (Illumina, San Diego, CA) at the Emory University Biomarker Service Center. This array was designed to target polymorphic SNPs from across the sunflower genome. Details about the development of this array have been previously provided by Bachlava et al. . Genome studio ver. 2011.1 (Illumina) was used to make SNP calls and map positions were assigned based on the consensus genetic map of Bowers et al. . Of a total of 9,480 SNPs on the array, 5,788 genetically mapped polymorphic SNPs could be reliably scored as apparently single copy loci in our population . Further, only SNPs with a minor allele frequency of ≥10% (5,359) was used for association mapping analyses . The previously identified B locus was not directly used as a marker in this study since its presence or absence was not known for all the lines used in this study. However, markers spanning this region (as determined based on marker position on the sunflower consensus map; Bowers et al. ) were included in our analyses.
Association analyses were performed using TASSEL ver. 3.0 . Because such analyses are prone to false positives due to unrecognized kinship and/or population structure , three different mixed models were employed. The first model corrected for kinship (K, estimated using SPAGeDi)  only. The second and third models corrected for both kinship and population structure. In the Q + K model, population structure (Q) was estimated using STRUCTURE ver. 2.2 . In the P + K model, population structure (P) was estimated via a principal coordinate analysis (PCoA) using GenAlEx ver. 6.41 . Details of the underlying SPAGeDi, STRUCTURE, and PCoA analyses can be found in Mandel et al. .
Following the association analyses, Q-Q (quantile-quantile) plots were constructed for each of the three models and compared to the results of a naïve model to select the most appropriate model for analysis. This was performed for apical (1st quarter), mid (2nd and 3rd quarter) and basal (4th quarter) branching. For association mapping of apical branching, the total number of branches in the first quarter from lines that displayed apical, mid-apical, and whole plant branching were included in the analysis. Similarly for all other branching types, the number of branches from all lines that displayed branching in the respective quarter were included in our analyses. Because non-independence of linked markers can result in highly conservative significance thresholds , we set a threshold -log(P) value of 3.6 (alpha=0.05, P=0.00025, log 1/P=3.60) to identify significant associations using a multiple testing correction method that accounts for correlation among markers while also controlling the type I error rate .
Identification and mapping of candidate branching genes
In order to identify putative sunflower orthologs of branching genes identified in other plant species (i.e., Arabidopsis, rice, pea, and petunia), we searched the literature for genes involved in axillary meristem initiation and outgrowth. These genes included transcription factors such as REVOLUTA (REV), LATERAL SUPPRESSOR (LAS) and REGULATORS OF AXILLARY MERSITEMS (RAX 1, 2, and 3) from Arabidopsis that have been shown to play important roles during initiation of axillary meristem and bud formation [44-46]. In addition, genes associated with homeostasis and signaling of phytohormones and growth regulators such as auxin, CK, SL, GA, and polyamines were included [12-15,17,18]. Many other genes involved in branch outgrowth that encode transcription factors, cytochrome P450, MAP KINASE KINASE 7 (MAPKK7), arabinogalactan proteins, and other DNA binding proteins were also included. For genes that were a part of a multi-gene family, only the genes that have been implicated in branching have been included. This resulted in the identification of 48 candidate genes (Additional file 1).
Once identified, the sequence of each candidate gene was searched against v0.1 of the draft sunflower genome assembly (http://www.sunflowergenome.org/) using tblastx with an E-value threshold of 1e-6. Since sunflower has undergone at least three whole genome duplication events  and is also likely to have experienced numerous segmental duplications, up to eight sunflower blast hits were considered for each of the candidate genes to allow for multiple possible homologs. Map positions of as many of these genes as possible were then determined by comparison to genetically mapped contigs from the whole genome assembly using the whole genome shotgun (WGS) sequence-based sunflower genetic map . Genetic positions from the WGS map were translated into map positions on the Bowers et al.  consensus genetic map using common markers. This allowed us to determine if any of the candidate genes mapped in close proximity to SNPs associated with the branching traits. We used a window size of 2.5 cM to determine co-localization since the SNPs in question had been previously ordered using multiple mapping populations differing from the population used to order the WGS map. As a result, an uncertainty of several cM in the exact genetic position of markers remained after map integration. The sunflower contigs containing the candidate genes (based on the blast results) and their position on LGs are listed in Additional file 2. To determine if candidate genes and significant associations occurred in regions of high or low LD, r 2 values were computed using the diversity panel between all SNP pairs within 2.5 cMs of each other. An overall r 2 value was then computed for each SNP by averaging the individual pairwise r 2 values.
Association mapping of branching patterns
Phenotypic analyses revealed extensive variation in the type and extent of branching within the association mapping population (Table 1, Figure 2, Additional file 3). Depending on the location, 89–102 lines exhibited whole-plant branching, while 70–110 lines exhibited no branching. The numbers of lines that exhibited other types of branching are listed in Table 1. There was high overlap of 77 lines exhibiting whole plant branching at all three locations (Table 2). However, all other branching patterns exhibited extensive variation across locations also displayed as a heat map (Figure 2). There was a greater overlap of lines exhibiting basal branching in GA and BC compared to IA, where many did not branch. Collectively these data illustrate that the environment can exert a strong influence on branching patterns. Interestingly, unbranched lines did not display any axillary bud initiation.
As expected, the correlation analyses revealed positive correlations amongst branching types and across locations (Table 3). However, basal branching tended to exhibit lower correlations overall, both within and across locations. This suggests that basal branching may be under different genetic control, and might also be more subject to environmental influences than apical and mid branching. In the PCA (Figure 3), the first dimension accounted for 77-83% (depending on location) of the underlying variation, and largely reflected variation in the extent of branching (i.e., the number of branches produced). The second dimension captured 12-16% of the phenotypic variation and primarily reflected differences in apical vs. basal branching.
The P + K model, which corrected for both kinship and population structure (using PCoA results) appeared to be the most conservative model across traits and showed the lowest tendency toward false positives (Additional file 4). This model was thus selected for all subsequent association analyses. Our analyses detected significant associations on multiple LGs (Table 4). In total, SNPs on 14 of the 17 sunflower LGs groups (all but LGs 7, 11, and 14) showed significant associations with various branching types. The Manhattan plots (Figure 4) illustrate our results for the various branching types at all three locations. Overall, the largest number of significant SNPs could be found on LG 10 in a region that was mostly associated with branching in the apical and mid regions of the plant. The broad peak of associations that is visible on the upper portion of LG 10 (referred to herein as 10a; Table 4; Figure 4) corresponds to the so-called B locus. In general, we were able to identify SNPs associated with basal and mid-basal branching which were distinct from apical, mid-apical, and mid branching (Table 4; Figure 4). For example, certain LGs (such as LG 4 in GA and BC and LG 8 in GA and IA) were associated with basal and mid-basal branching whereas SNPs near the bottom of LG 6 (in GA, IA and BC) and 13 (IA and BC) were associated with branching in the apical and mid regions. Only one secondary branching association was observed on LG 10 (in BC). The results for whole plant branching are not presented here, as that information was analyzed by Mandel et al. .
Identification of candidate genes associated with various branching phenotypes
A total of 48 candidate branching genes were identified from other species (Additional file 1), and homologs to 39 of these genes could be identified in sunflower. This attrition is likely due to the use of an incomplete, draft assembly of the genome. We allowed up to eight of the best blast hits for each gene which resulted in 278 homologs. Of these, 153 (corresponding to 38 genes) could be placed on the sunflower genetic map (Additional files 2 and 5). Among these 153 blast hits, 92 unique contigs were identified. In many cases, multiple genes (belonging to gene families) were found within the same sunflower contig.
In comparing the map positions of these genes to the locations of SNPs exhibiting significant associations, we identified 13 potential candidate branching genes (from the larger set of 92 unique genes) that mapped in close proximity to various branching traits (Figure 5). Genes of interest involved in hormone related pathways included an auxin biosynthetic gene (YUCCA; YUC), an auxin induced gene (AUXIN/INDOLE-3-ACETIC ACID; Aux/IAA), and genes involved in CK biosynthesis (ISOPENTYLTRANSFERASE; IPT) and degradation (CK OXIDASE/DEHYDROGENASE; CKX), GA catabolism (GIBBERELLIN 2-OXIDASES; GA2ox), and SL biosynthesis (DWARF27; D27). Several transcription factors also co-localized with significant branching associations, including CUP-SHAPED COTYLEDON (CUC), LATERAL SHOOT INDUCING FACTOR (LIF), BRANCHED (BRC), as well as a histone methyltransferase, SET DOMAIN GROUP8 (SDG8). Of these, eight candidate genes were in regions of elevated LD (Figures 5 and 6; Additional file 6), including all genes on LG 10, one on the lower half of LG 13, and one on LG 16. The remaining genes on LGs 4, 5, 6, and the upper half of 13 were present in regions of lower LD (i.e., r 2 < 0.20).
The identification of loci influencing specific branching patterns has the potential to facilitate the manipulation of plant architecture, which can influence yield and seed/fruit quality in crops [15,49-52]. Here, we have identified numerous distinct loci that influence apical vs. basal branching in sunflower. It is important to note, however, that the occurrence of variation in branching patterns within lines across the three locations suggests that environmental variation also plays an important role in determining sunflower branching architecture. This conclusion is consistent with the results of previous studies on the effect of external factors on branch formation. For example, increased planting density can result in a suppression of branching [53,54], and photoperiod is also known to affect branching [55-58]. Other environmental factors influencing branching include light levels and quality [59-61], plant nutrition status, and availability of nutrients [62,63]. It seems likely that some combination of these factors influenced our results across locations.
The formation of a branch involves two developmental processes: the initiation of an axillary bud and its elongation into a branch [5,64]. As noted above, genes such as REV, LAS, and RAX in Arabidopsis and their orthologs in rice and tomato play a role during axillary meristem initiation and bud formation [44-46,65-67]. Because the non-branching lines in our study also did not exhibit axillary bud initiation, we were unable to separate initiation from outgrowth and were thus unable to identify loci contributing specifically to one process or the other. The association mapping population did, however, exhibit substantial diversity in branching patterns which facilitated the identification of multiple branching-related loci that had not been previously identified in sunflower. These included associations on LGs 1, 2b, 5, 12, 14, 15, and 16 not previously identified via traditional QTL mapping [19,20,68] or via association mapping .
Importantly, our results also revealed that apical and basal branching are under largely independent genetic control (Figure 4; Table 4). Previous studies have identified a small number of loci that influence either apical or basal branching [22,23] and suggested that these traits are controlled by different loci, but our results indicate that the effects of branching-related loci are largely compartmentalized in sunflower. In fact, with the exception of loci on LGs 6a, 10a, and 10b, all three of which are associated with whole plant branching, the remaining 20 significant associations had primarily apical or basal effects. As mentioned above, the broad peak of associations on LG 10a (Figure 4) corresponds to the so-called B locus . This locus is responsible for the branching that was reintroduced into the sunflower gene pool to extend the flowering time in R lines, and is present as a large haplotypic block within many of these lines , resulting in a sizable island of elevated LD (Figure 6). Similar instances of elevated LD, mostly associated with past episodes of selection during the evolution of cultivated sunflower, are also visible elsewhere in the genome . While multiple novel branching loci were identified in this study, it is possible that segregation of the B locus masked the effect of other loci contributing to the variation observed within this population. The analysis of a population that varies for branching while being fixed for the recessive, unbranched b/b genotype at this locus has the potential to shed light on this possibility.
We further identified candidate branching genes that are co-localized with significant branching associations. Similar approaches using association mapping have been used to identify positional candidates for important agronomic traits in crops such as rice (e.g., [69,70]) and maize (e.g., ), to identify candidate genes associated with flowering time, pathogen resistance, and tolerance to salinity in Arabidopsis (e.g., [72,73]). Of course, the resolution afforded by association mapping varies across the genome due to localized variation in the extent of LD. We thus focused on the identification of candidate genes that co-localized with branching associations in regions of low LD (i.e., r 2 < 0.20). These genes included homologs of AtYUC that were closely linked to mid-basal and basal branching on LG 4 and mid-apical and mid-basal branching on LG 6; a homolog of CKX, which co-localized with apical and mid-branching on LG 13; and homologs of GA2ox5 and GA2ox6 in that same region of LG 13. In addition to these hormone-related genes, a homolog of SDG8 that was associated with mid-basal branching co-localized on LG 5. SDG8 is thought to epigenetically regulate other branching genes, and loss of this gene in Arabidopsis exhibit increased shoot branching . The above genes are thus excellent candidates for further functional characterization. In maize, the primary determinant of branched vs. unbranched is teosinte branched1 (tb1), which is a transcription factor within the CYCLOIDEA (CYC)/TB1 subfamily of the TEOSINTE-BRANCHED1/CYCLOIDEA/PCF (TCP) transcription factor family [75-77]. While homologs of this gene are known to influence branching in other species such as rice  and Arabidopsis , we found no evidence for a similar role of CYC-like genes in sunflower (see also Chapman et al. ).
In this study, we identified numerous loci with significant effects on branching in sunflower, many of which had variable effects across environments. This includes multiple genomic regions that had not previously been implicated in branching. Interestingly, the majority of these loci primarily affected either apical or basal branching, as opposed to influencing branching at the whole-plant level. We also identified a collection of branching-related candidate genes that co-localized with significant association in regions of low LD, providing us with a pool of promising candidates for future functional validation.
- A lines:
Cytoplasmic male-sterile female lines
- Aux/IAA :
- BRC :
- CKX :
- CUC :
- CYC :
- D27 :
- GA2ox :
Germplasm Resources Information Network
French National Institute for Agricultural Research
- IPT :
- LAS :
- LIF :
LATERAL SHOOT INDUCING FACTOR
- LS Means:
- MAPKK7 :
MAP KINASE KINASE 7
North Central Regional Plant Introduction Station
Population structure estimated using principal coordinate analysis
Principal component analysis
Principal coordinate analysis
Population structure estimated using STRUCTURE
- Q-Q plots:
Quantitative trait locus
- R lines:
Male-fertile restorer lines
- RAX :
REGULATORS OF AXILLARY MERSITEMS
- REV :
- SDG8 :
SET DOMAIN GROUP8
- tb1 :
teosinte branched 1
Whole genome shotgun
Irwin DL, Aarssen LW. Effects of nutrient level on cost and benefit of apical dominance in Epilobium ciliatum. Am Midl Nat. 1996;136:14–28.
Jaremo J, Nilsson P, Tuomi J. Plant compensatory growth: herbivory or competition? Oikos. 1996;77:238–47.
Fay PA, Throop HL. Branching responses in Silphium Integrifolium (Asteraceae) following mechanical or gall damage to apical meristems and neighbor removal. Am J Bot. 2005;92:954–9.
Evers JB, Van der Krol AR, Vos J, Struik PC. Understanding shoot branching by modelling form and function. Trends Plant Sci. 2011;16:464–7.
Ward SP, Leyser O. Shoot branching. Curr Opin Plant Biol. 2004;7:73–8.
Hammer K. Das Domestikationssyndrom. Kulturpflanze. 1984;32:11–34.
Heiser CB. The sunflower among the North American Indians. Proc Am Philos Soc. 1951;95:432–48.
Heiser CB, Smith DM, Clevenger SB, Martin WC. The North American sunflowers (Helianthus). Mem Torrey Bot Club. 1969;22:1–218.
Rogers CE, Thompson TE, Seiler GJ. Sunflower species of the United States. In National Sunflower Association, Bismarck, North Dakota; 1982:75.
Rieseberg LH, Seiler GJ. Molecular evidence and the origin and development of the domesticated sunflower (Helianthus annuus, Asteraceae). Econ Bot. 1990;44:79–91.
Fick GN, Miller JF. Sunflower breeding. In: AA S, editor. Sunflower technology and production. Madison: American Society of Agronomy; 1997. p. 395–440.
Ongaro V, Leyser O. Hormonal control of shoot branching. J Exp Bot. 2008;59:67–74.
Ferguson BJ, Beveridge CA. Roles for auxin, cytokinin, and strigolactone in regulating shoot branching. Plant Physiol. 2009;149:1929–44.
Shimizu-Sato S, Tanaka M, Mori H. Auxin-cytokinin interactions in the control of shoot branching. Plant Mol Biol. 2009;69:429–35.
Yaish MWF, Guevara DR, El-Kereamy A, Rothstein SJ. Axillary shoot branching in plants. In: Pua EC, Davey MR, editors. Plant developmental biology–biotechnological perspectives:. Volume 1. Berlin, Heidelberg: Springer Berlin Heidelberg; 2010. p. 37–52.
McSteen P. Hormonal regulation of branching in grasses. Plant Physiol. 2009;149:46–55.
Lo SF, Yang SY, Chen KT, Hsing YI, Zeevaart JA, Chen LJ, et al. A novel class of gibberellin 2-oxidases control semidwarfism, tillering, and root development in rice. Plant Cell. 2008;20:2603–18.
Ge C, Cui X, Wang Y, Hu Y, Fu Z, Zhang D, et al. BUD2, encoding an S-adenosylmethionine decarboxylase, is required for Arabidopsis growth and development. Cell Res. 2006;16:446–56.
Burke JM, Tang S, Knapp SJ, Rieseberg LH. Genetic analysis of sunflower domestication. Genetics. 2002;161:1257–67.
Wills DM, Burke JM. Quantitative trait locus analysis of the early domestication of sunflower. Genetics. 2007;176:2589–99.
Putt ED. Recessive branching in sunflowers. Crop Sci. 1964;4:444–5.
Gentzbittel L, Mestries E, Mouzeyar S, Mazeyrat F, Badaoui S, Vear F, et al. A composite map of expressed sequences and phenotypic traits of the sunflower (Helianthus annuus L.) genome. Theor Appl Genet. 1999;99:218–34.
Hockett EA, Knowles PF. Inheritance of branching in sunflowers, Helianthus annuus L. Crop Sci. 1970;10:432–6.
Tang S, Leon A, Bridges WC, Knapp SJ. Quantitative trait loci for genetically correlated seed traits are tightly linked to branching and pericarp pigment loci in sunflower. Crop Sci. 2006;46:721–34.
Tang S, Yu JK, Slabaugh B, Shintani K, Knapp J. Simple sequence repeat map of the sunflower genome. Theor Appl Genet. 2002;105:1124–36.
Zhu C, Gore M, Buckler ES, Yu J. Status and prospects of association mapping in plants. Plant Genome. 2008;1:5–20.
Yu J, Buckler ES. Genetic association mapping and genome organization of maize. Curr Opin Biotechnol. 2006;17:155–60.
Mandel JR, Dechaine JM, Marek LF, Burke JM. Genetic diversity and population structure in cultivated sunflower and a comparison to its wild progenitor, Helianthus annuus L. Theor Appl Genet. 2011;123:693–704.
Mandel JR, Nambeesan S, Bowers JE, Marek LF, Ebert D, Rieseberg LH, et al. Association mapping and the genomic consequences of selection in sunflower. PLoS Genet. 2013;9:e1003378. 10.1371/journal.pgen.1003378.
Bachlava E, Taylor CA, Tang S, Bowers JE, Mandel JR, Burke JM, et al. SNP discovery and development of a high-density genotyping array for sunflower. PLoS One. 2012;7:e29814. 10.1371/journal.pone.0029814.
Bowers JE, Bachlava E, Brunick RL, Rieseberg LH, Knapp SJ, Burke JM. Development of a 10,000 locus genetic map of the sunflower genome based on multiple crosses. Genes Genomes Genet. 2012;2:721–9.
Schneiter AA, Miller JF. Description of sunflower growth stages. Crop Sci. 1981;21:901–3.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
Le S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008;25:1–18.
R Core Team: R: A language and environment for statistical computing. 2014.
Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.
Pritchard JK, Stephens M, Rosenberg N, Donnelly P. Association mapping in structured populations. Am J Hum Genet. 2000;67:170–81.
Hardy OJ, Vekemans X. SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002;2:618–20.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
Peakall R, Smouse PE. Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.
De Silva H, Ball R. Linkage disequilibrium mapping concepts. In: Oraguzie NC, Rikkerink EHA, Gardiner SE SH, editors. Association mapping in plants. New York: Springer; 2007. p. 103–32.
Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol. 2008;32:361–9.
Greb T, Clarenz O, Schafer E, Muller D, Herrero R, Schmitz G, et al. Molecular analysis of the LATERAL SUPPRESSOR gene in Arabidopsis reveals a conserved control mechanism for axillary meristem formation. Genes Dev. 2003;17:1175–87.
Otsuga D, DeGuzman B, Prigge MJ, Drews GN, Clark SE. REVOLUTA regulates meristem initiation at lateral positions. Plant J. 2001;25:223–36.
Muller D, Schmitz G, Theres K. Blind homologous R2R3 Myb genes control the pattern of lateral meristem initiation in Arabidopsis. Plant Cell. 2006;18:586–97.
Barker MS, Kane NC, Matvienko M, Kozik A, Michelmore RW, Knapp SJ, et al. Multiple paleopolyploidizations during the evolution of the Compositae reveal parallel patterns of duplicate gene retention after millions of years. Mol Biol Evol. 2008;25:2445–55.
Renaut S, Grassa CJ, Yeaman S, Moyers BT, Lai Z, Kane NC, et al. Genomic islands of divergence are not affected by geography of speciation in sunflowers. Nat Commun. 2013;4:1827.
Ledger SE, Janssen BJ, Karunairetnam S, Wang T, Snowden KC. Modified CAROTENOID CLEAVAGE DIOXYGENASE8 expression correlates with altered branching in kiwifruit (Actinidia chinensis). New Phytol. 2010;188:803–13.
Cieslak M, Seleznyova AN, Hanan J. A functional-structural kiwifruit vine model integrating architecture, carbon dynamics and effects of the environment. Ann Bot. 2011;107:747–64.
Lauri PA, Costes E, Regnard JL, Brun L, Simon S, Monney P, et al. Does knowledge on fruit tree architecture and its implications for orchard management improve horticultural sustainability? An overview of recent advances in the apple. Acta Hortic. 2009;817:243–50.
Upadyayula N, Da Silva HS, Bohn MO, Rocheford TR. Genetic and QTL analysis of maize tassel and ear inflorescence architecture. Theor Appl Genet. 2006;112:592–606.
Casal JJ, Sanchez RA, Deregibus VA. The effect of plant density on tillering: The involvement of R/FR ratio and the proportion of radiation intercepted per plant. Environ Exp Bot. 1986;26:365–71.
Aguilar-Martínez JA, Poza-Carrión C, Cubas P. Arabidopsis BRANCHED1 acts as an integrator of branching signals within axillary buds. Plant Cell. 2007;19:458–72.
Arumingtyas EL, Floyd R, Gregory M, Murfet I. Branching in Pisum: inheritance and allelism tests with 17 ramosus mutants. Pisum Genet. 1992;24:17–31.
Napoli CA, Beveridge CA, Snowden KC. Reevaluating concepts of apical dominance and the control of axillary bud outgrowth. Curr Top Dev Biol. 1999;44:127–69 [Current Topics in Developmental Biology].
Grbić V, Bleecker AB. Axillary meristem development in Arabidopsis thaliana. Plant J. 2000;21:215–23.
Stirnberg P, Van De Sande K, Leyser HMO. MAX1 and MAX2 control shoot lateral branching in Arabidopsis. Development. 2002;129:1131–41.
Kebrom TH, Burson BL, Finlayson SA. Phytochrome B represses teosinte branched1 expression and induces sorghum axillary bud outgrowth in response to light signals. Plant Physiol. 2006;140:1109–17.
Snowden KC, Napoli CA. A quantitative study of lateral branching in petunia. Funct Plant Biol. 2003;30:987–94.
Cline MG. Exogenous auxin effects on lateral bud outgrowth in decapitated shoots. Ann Bot. 1996;78:255–66.
Napoli CA, Ruehle J. New mutations affecting meristem growth and potential in Petunia hybrida Vilm. J Hered. 1996;87:371–7.
Cline MG. Apical dominance. Bot Rev. 1991;57:318–58.
Shimizu-Sato S, Mori H. Control of outgrowth and dormancy in axillary buds. Plant Physiol. 2001;127:1405–13.
Schumacher K, Schmitt T, Rossberg M, Schmitz G, Theres K. The Lateral suppressor (Ls) gene of tomato encodes a new member of the VHIID protein family. Proc Natl Acad Sci U S A. 1999;96:290–5.
Li X, Qian Q, Fu Z, Wang Y, Xiong G, Zeng D, et al. Control of tillering in rice. Nature. 2003;422:618–21.
Schmitz G, Tillmann E, Carriero F, Fiore C, Cellini F, Theres K. The tomato Blind gene encodes a MYB transcription factor that controls the formation of lateral meristems. Proc Natl Acad Sci U S A. 2002;99:1064–9.
Dechaine JM, Burger JC, Chapman MA, Seiler GJ, Brunick R, Knapp SJ, et al. Fitness effects and genetic architecture of plant-herbivore interactions in sunflower crop-wild hybrids. New Phytol. 2009;184:828–41.
Huang X, Wei X, Sang T, Zhao Q, Feng Q, Zhao Y, et al. Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet. 2010;42:961–7.
Zhao K, Tung CW, Eizenga GC, Wright MH, Ali ML, Price AH, et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat Commun. 2011;2:467.
Weng J, Xie C, Hao Z, Wang J, Liu C, Li M, et al. Genome-wide association study identifies candidate genes that affect plant height in Chinese elite maize (Zea mays L.) inbred lines. PLoS One. 2011;6:e29229. doi:10.1371/journal.pone.0029229.
Aranzana MJ, Kim S, Zhao K, Bakker E, Horton M, Jakob K, et al. Genome-wide association mapping in Arabidopsis identifies previously known flowering time and pathogen resistance genes. PLoS Genet. 2005;1:e60. 10.1371/journal.pgen.0010060.
DeRose-Wilson L, Gaut BS. Mapping salinity tolerance during Arabidopsis thaliana germination and seedling growth. PLoS One. 2011;6:e22832. 10.1371/journal.pone.0022832.
Dong G, Ma D, Li J. The histone methyltransferase SDG8 regulates shoot branching in Arabidopsis. Biochem Biophys Res Commun. 2008;373:659–64.
Doebley J, Stec A, Gustus C. teosinte branched1 and the origin of maize: Evidence for epistasis and the evolution of dominance. Genetics. 1995;141:333–46.
Doebley J, Stec A, Hubbard L. The evolution of apical dominance in maize. Nature. 1997;386:485–8.
Cubas P, Lauter N, Doebley J, Coen E. The TCP domain: a motif found in proteins regulating plant growth and development. Plant J. 1999;18:215–22.
Takeda T, Suwa Y, Suzuki M, Kitano H, Ueguchi-Tanaka M, Ashikari M, et al. The OsTB1 gene negatively regulates lateral branching in rice. Plant J. 2003;33:513–20.
Chapman MA, Leebens-Mack JH, Burke JM. Positive selection and expression divergence following gene duplication in the sunflower CYCLOIDEA gene family. Mol Biol Evol. 2008;25:1260–73.
We thank Eric Elsner, Jenny Leverett, Irvine Larsen, Tanya MacInnes, Anh Nguyen, and many undergraduate workers for their assistance in the field. This research was supported by USDA National Institute of Food and Agriculture (2008-35300-19263), the NSF Plant Genome Research Program (DBI-0820451), the USDA/DOE Plant Feedstocks Genomics Joint Program (ER 64664), Genome BC, and Genome Canada (http://www.genomecanada.ca).
The authors declare that they have no competing interests.
SUN, JRM, LHR, SJK, and JMB conceived the study and designed the experiments. JRM and JMB were involved in the development of the association mapping population. SUN, JRM, LFM, and DE performed phenotypic analyses in the field. Data analyses were performed by SUN, JRM, JEB, LFM, DE, and JC. The manuscript was prepared by SUN, JRM, and JMB. All the authors read and approved the final manuscript.
Candidate branching genes. Genes involved in axillary meristem initiation and branch outgrowth in Arabidopsis thaliana (At), Petunia x hybrida (Ph), Pisum sativum (Ps), and Oryza sativa (Os).
List of candidate genes with significant hits to contigs within the draft sunflower genome assembly. For each candidate gene, up to eight sunflower blast hits were considered, with the actual number being indicated as the hit number. The positions of the mapped genes within the contigs are also listed. See text for details. Additional file 1 about the candidate genes.
Phenotypic data for branching in the association mapping population. Also included are the classification (i.e., oil vs. nonoil, breeding group, etc.) and USDA designation of each line. For more specific descriptions of these lines, refer to Mandel et al. . Locations are G=Georgia, I=Iowa, and C=British Columbia. Q indicates quarter. The various branching types are illustrated in Figure 1.
Quantile-quantile (Q-Q) plots of branching associations based on models accounting for the effects of kinship and/or population structure vs. a naive model. Results for branching in each of the four quarters at all three locations were computed using three models, K (red), P + K (blue), and Q + K (grey) and compared against the results of a naive model (brown). The expected values are shown in green. Locations are indicated as G=Georgia, I=Iowa, and C=British Columbia. Q indicates quarter. The branching types are illustrated in Figure 1.
Sunflower genetic map depicting the positions of significant SNPs associated with branching traits and candidate branching genes. Shown are all 17 LGs and positions of all significant SNPs found in the study. Candidate branching genes are in black. Branching traits associated with the SNPs are color coded as in Figure 5 with apical (A; in green), mid-apical (MA; in dark blue), mid (M; in orange), mid-basal (MB; in magenta), basal (B; in brown), and secondary branching (2B; in light blue). Locations are indicated as follows: G=Georgia, I=Iowa, and C=British Columbia). For LG 10 only, SNPs associated with all branching traits other than 2B are indicated as whole plant branching (WPB) and the numbers denote differences among geographical locations for a given branching type.
Average r 2 values for SNPs along LGs that contained candidate branching genes that co-localized with significant branching associations. Average r 2 values were calculated as described in Figure 6. SNPs exhibiting significant branching associations and located in close proximity to a candidate branching gene are depicted in red. The positions of the putative candidate genes are indicated by blue arrows on the x-axis; letters refer to specific genes.