Two Novel AP2/EREBP Transcription Factor Genes TaPARG Have Pleiotropic Functions on Plant Architecture and Yield-Related Traits in Common Wheat

AP2/EREBPs play significant roles in plant growth and development. A novel, pleiotropic TaPARG (PLANT ARCHITECTURE-RELATED GENE), a member of the AP2/EREBP transcription factor gene family, and its flanking sequences were isolated in wheat (Triticum aestivum L.). Two TaPARG genes were identified and named as TaPARG-2A and TaPARG-2D. Their amino acid sequences were highly similar especially in the functional domains. TaPARG-2A on chromosome 2A was flanked by markers Xwmc63 and Xgwm372. TaPARG-2D was mapped to chromosome 2D. Subcellular localization revealed that TaPARG-2D was localized in the nucleus. The results of tissue expression pattern, overexpression in rice, association analysis and distinct population verification jointly revealed that TaPARG functions during the entire growth cycle of wheat. Its functions include regulation of plant architecture-related and yield-related traits. Association analysis, geographic distribution and allelic frequencies suggested that favored haplotypes Hap-2A-2 and Hap-2A-3 were selected in Chinese wheat breeding programs. Both favored haplotypes might be caused by a single amino acid substitution (His/Tyr). These results suggest that TaPARG is a regulatory factor in plant growth and development, and that the favored alleles might be useful for improving plant architecture and grain yield of wheat.


INTRODUCTION
AP2/EREBPs (ethylene-responsive element-binding proteins) belong to a superfamily of plantspecific transcription factors characterized by presence of an AP2 DNA-binding domain of 60 amino acids (Weigel, 1995;Okamuro et al., 1997). The highly conserved amino acid residues of AP2 domains were also called the GCC-box binding domain (GBD), because of specific binding to the GCC-box (Ohme-Takagi and Shinshi, 1995;Allen et al., 1998). The AP2/EREBP superfamily is divided into two subfamilies based on the number of AP2 domains (Riechmann and Meyerowitz, 1998). EREBPs containing a single AP2 domain are involved in regulatory networks of response to hormones, pathogen attack, and environmental signals involving DREBs (dehydration responsive element binding proteins) and ERFs (ethylene responsive factors) (Xu et al., 2011;Licausi et al., 2013;Jisha et al., 2015). EREBPs with two tandom AP2 domains (known as repeat 1 and repeat 2) have roles in the regulation of plant development; the AP2 subfamily is a good example (Zhao et al., 2006;Wang et al., 2014;Kuluev et al., 2015).
Previous studies demonstrated that members of the AP2 subfamily have roles in controlling flower and seed developmental processes. APETALA2 (AP2), with an important role in flowering and seed development in Arabidopsis encodes a putative transcription factor distinguished by a novel DNA binding motif referred to as the AP2 domain (Okamuro et al., 1997). In Arabidopsis AP2 is also involved in determination of seed size, seed weight, and accumulation of seed oil and protein (Jofuku et al., 2005). The AINTEGUMENTA (ANT) gene regulates flower development and is essential for ovule formation (Klucher et al., 1996). BABY BOOM (BBM) triggers conversion from vegetative to embryonic growth during embryogenesis (Boutilier et al., 2002). The PLETHORA1 (PLT1) and PLETHORA2 (PLT2) genes are key effectors for establishment of the stem cell niche during embryonic pattern formation in Arabidopsis (Aida et al., 2004). The AINTEGUMENTA-like (AIL) genes are expressed in young tissues and may specify meristematic or divison-competent states (Nole-Wilson et al., 2005). AINTEGUMENTA-like gene NtANTL, a putative ortholog of AtANT, controls cell division and expansion in controlling organ growth in tobacco (Kuluev et al., 2015). Extensive research has been conducted on the AP2/EREBP family of transcription factors in wheat, but has mainly focused on the DREB and ERF subfamilies because their overexpression enhances tolerance to abiotic and biotic stresses (Xu et al., 2011); however, some AP2 genes have been studied and showed their functions in plant growth and development. For example, gene Q which has a high degree of similarity to members of AP2 family has pleiotropically influences on many domestication-related traits, such as glume shape and tenacity, rachis fragility, spike length, plant height (PH), and spike emergence time (Simons et al., 2006). The engineering of cleistogamous wheat requires the presence of a functional TaAP2 modification at each of the three homoeologs (Ning et al., 2013).
In this study, the novel, pleiotropic gene, designated TaPARG (PLANT ARCHITECTURE-RELATED GENE), was isolated and characterized in wheat. Our findings suggest that TaPARGs play a key role in growth and development, including regulation of plant architecture-related traits such as PH, tiller number, and yield-related traits, including fertility and 1,000 kernel weight (TKW). The coding region of TaPARG-2A was analyzed, and sequence polymorphisms among common wheat cultivars were identified. Finally, functional markers for various haplotypes were developed, and favored haplotypes were identified by association analysis of two germplasm populations. The favored haplotypes Hap-2A-2 and Hap-2A-3 were selected in Chinese wheat breeding programs. Both favored haplotypes probably arose by a single amino acid change (His/Tyr). Our overall results suggest that TaPARGs are pleiotropic genes involved in growth and development, and that two functional markers developed for TaPARG-2A might be helpful for marker-assisted selection of ideal plant architecture and high-yielding wheat genotypes.

Plant Materials and Measurement of Agronomic Traits
Hanxuan 10 is a drought-tolerant cultivar released in 1966 and still grown in arid and barren areas, and also one of the strongly drought-tolerant accessions selected from more than 20,000 wheat germplasm resources. Therefore, Hanxuan 10 is used as the plant material for cloning genes TaPARG-2A and TaPARG-2D, and to analysis the polymorphism of TaPARGs and its relation with the agronomic traits under diverse water conditions. Lumai 14 is a high-yielding cultivar adapted to abundant water and fertile conditions, and is widely grown during the 1990s in northern China. Using Hanxuan 10 and Lumai 14 as parents, we developed a doubled haploid (DH) population and mapped QTLs for a list of important agronomic traits (Yang et al., 2007;Wang et al., 2010;Wu et al., 2011Wu et al., , 2012. Hanxuan 10 was used for gene cloning and expression analysis. Thirtyfour diverse cultivars (Supplementary Table S1) were used to identify sequence polymorphism of TaPARG. A set of nullisomictetrasomic and ditelosomic lines of Chinese Spring were used for chromosome location of the gene. A doubled haploid (DH) population derived from the cross Hanxuan 10 × Lumai 14 (Yang et al., 2007) was employed for genetic mapping. Kitaake, a short life cycle rice cultivar, was used for genetic transformation.
Three common wheat germplasm populations were used in association studies of TaPARG haplotypes and phenotypic traits. The Population 1 (262 accessions) was firstly used; the accessions were mainly from the Northern Winter Wheat and Yellow and Huai River Valleys Facultative Wheat Zones . Population 2 (157 landraces) and Population 3 (348 modern cultivars) including genotypes from all 10 Chinese wheat production zones were used to verify the results of the initial association analysis, allelic frequencies and geographic distribution analysis. Population 2 entries mainly came from the Chinese wheat mini-core collection, which represents more than 70% of the total genetic diversity of the Chinese wheat germplasm collection, and Population 3 entries were from the Chinese wheat core collection (Hao et al., 2008(Hao et al., , 2011 The field trial was managed in randomized complete blocks with three replicates. Each plot contained three 1.5 m rows spaced 15 cm apart, with 15 plants in each row. At mature stage, 10 plants were randomly sampled from the middle position of each plot to measure tiller number and 1,000 kernel weight, the main panicles were used for measuring plant height and seed setting rate.

Subcelluar Localization of TaPARG Protein
TaPARG-2D was fused upstream of GFP in the pCAMBIA1300 vector under control of the CaMV 35S promoter. Primers 5 -CCC AAGCTTATGACCAACAACAACGGCA-3 (HindIII site underlined) and 5 -GGGGTACCTGCGTCGCTCCACGC-3 (KpnI site underlined) were used for sub-cloning. The constructs were transferred into wheat mesophyll protoplasts for subcellular localization in wheat protoplasts by the PEG-mediated method (Yoo et al., 2007). After incubation at 25 • C for 16 h, fluorescence signals were detected using a laser scanning confocal microscope (Leica TCS-NT, Germany). For observations of subcellular localization in tobacco leaf cells, constructs with the coding regions of TaPARG-2D and GFP as control, were transferred into tobacco leaves by Agrobacterium-mediated transformation (EHA105) (Liu et al., 2010), and cultured for 4 days at 22 • C in a 16 h light/8 h darkness photoperiod.

Analysis of Gene Expression
Total RNA was extracted with TRIzol, and cDNA was synthesized with a SuperScript R Double-Stranded cDNA Synthesis Kit (Invitrogen Com., USA). Semi-quantitative RT-PCR was performed to analyze the expression pattern of TaPARG in various tissues. Primers P-TaPARG (5 -GGC CCTTCCACCCATATCA-3 and 5 -GGTCACGCCACGGTA CATG-3 ) and P-tubulin (5 -TGAGGACTGGTGCTTACCGC-3 and 5 -GCACCATCAAACCTCAGGGA-3 ) were used to amplify TaPARGs and the tubulin gene, respectively. Because TaPARG-2A and TaPARG-2D were highly similar in sequence, it was hard to design specific primers to separately amplify these two copies, therefore, the common primer P-TaPARG was used to amplify both TaPARG-2A and TaPARG-2D in wheat tissue expression analysis. The target PCR product size of P-TaPARG was 150 bp.
Different wheat tissues including roots, root bases, leaf sheaths, leaf blades, nodes, internodes, spikes, and plumules were collected to examine the expression pattern of TaPARG at various developmental stages (germination, seedlings, and heading).

Functional Marker Development
Thirty-four accessions were initially chosen to detect nucleotide polymorphisms in TaPARG-2A coding and flanking regions. Three haplotypes were identified based on five SNPs (single nucleotide polymorphisms). Two specific functional markers (TaPARM1 and TaPARM2) were developed for the two SNPs in the exons. The primers of TaPARM1 were 5 -ACGCCACAGCCGAGAACGTCTTG-3 and 5 -GG GGTGTACCTGGTGACGCCACAGT-3 , and the primers of TaPARM2 were 5 -GTTTTCGCACTGTCCCAAATC-3 and 5 -GGACATGACGTCCTGGTTTC-3 . Genotyping was implemented by two rounds of PCR, the first was to amplify fragments from chromosome 2A in all cultivars with the genomic-specific primer pair F2/R2; the second was performed as follows: first round PCR products were diluted 100 times, and 1 µl of diluted solution was used as template for the second round using primers TaPARM1 and TaPARM2 (annealing at 55 • C for 30 s, and extension at 72 • C for 30 s). The PCR products were digested by the corresponding restriction enzymes and separated in 4% agarose gels.

Population Structure and Association Analysis
Using data from 209 whole-genome SSR markers , population structure was determined by STRUCTURE v2.3.4. Ten subpopulations (k = 1 to 10) were set with a burnin period of 50,000 iterations and a run of 500,000 replications of Markov Chain Monte Carlo after burn in. The k method was applied according to LnP (D) in the STRUCTURE, and the output and result were estimated (Li et al., 2015). The general linear model (GLM) in TASSEL v2.1 was used for association analysis, and P < 0.01 denotes significant association. Statistical analysis was conducted by SPSS v17.0, statistical significance was calculated using one-way ANOVA.

Generation of Transgenic Plants in Rice
The cDNA of TaPARG-2D containing the entire ORF was inserted into pCUbi1390 vectors using primers 5 -GGGGTACC ATGACCAACAACAACGGCA-3 (KpnI site underlined) and 5 -CCCAAGCTTTCATGCGTCGCTCCACG-3 (HindIII site underlined). The construct was transferred into rice. Positive plants were identified with hygromycin, and T3 generation transgenics overexpressing TaPARG-2D were reconfirmed by PCR.

Cloning and Sequence Analysis of TaPARG
Two types of cDNA sequence of TaPARG were obtained from wheat cv. Hanxuan 10, and designated TaPARG-2A (GenBank AC: KU593503) and TaPARG-2D (GenBank AC: KU593504) according to their genomic origins. TaPARG-2A and TaPARG-2D have the same structure, containing eight exons and encoding proteins containing 601 and 600 amino acids, respectively. The identity of the two putative TaPARG proteins was 97.67%. The identities of TaPARG-2A with Oryza sativa, Glycine soja, and Arabidopsis thaliana were 53.10, 50.47, and 55.85%, respectively. Both TaPARGs contain two AP2 domains ( Figure 1A). A neighbor-joining phylogenetic tree of AP2/EREBP2 was constructed for the full-length protein sequences of the TaPARGs with Clustal W (Figure 1B). The TaPARGs and the orthologous gene from A. tauschii were grouped into a separate clade. Protein identity and phylogeny analysis suggested that TaPARGs might have different functions compared to other AP2/EREBPs.
Since the two TaPARGs were highly similar in sequence, particularly in the AP2 domains, the D genomic member, TaPARG-2D, was chosen for functional analysis.

Subcellular Localization of TaPARG-2D Protein
AP2/EREBPs have roles as transcription factors (Shigyo et al., 2006). Prediction of subcellular localization using ProtComp Version 9.0 software suggested that TaPARG-2D was a typical nuclear localized protein. To further address this point, the recombinant construct of TaPARG-2D-GFP fusion protein was expressed in wheat protoplasts and tobacco leaves. TaPARG-2D was specifically located in the nucleus (Figure 2).

Tissue Expression Pattern of TaPARGs in Wheat
Semi-quantitative RT-PCR was performed to analyze the tissue expression pattern of TaPARGs at different developmental stages in wheat (Figure 3). TaPARGs were constitutively expressed in different tissues at the germination, seedling, and heading stages, and were highly expressed in nodes and root bases. Constitutive expression suggested that TaPARGs might have active roles during the entire growth cycle.

Overexpression of TaPARG-2D in Rice Influences Overall Plant Development
Since TaPARG was constitutively expressed, we investigated the effects of overexpression in transformed rice plants harboring TaPARG-2D constructs driven by the CaMV 35S promoter. Eighteen transgenic lines were obtained. Three transgenic lines with different expression quantity were used for phenotype identification (Figure 4A; Supplementary Figure S1). Effects on plant growth ( Figure 4B) and fertility ( Figure 4C) were detected. Several biological and agronomic traits were measured to characterize the effects, including plant height (Figure 4D), tiller number per plant (Figure 4E), seed setting rate (Figure 4F), and 1,000 kernel weight ( Figure 4G). With increasing relative expression levels of TaPARG-2D in different transgenic lines, plant height was reduced (from 75.2 to 62.3 cm) and tiller number per plant increased (from 23 to 55). Seed setting rate declined dramatically (from 78.9 to 1.3%) and 1,000 kernel weight decreased slightly (from 27.3 to 24.3 g). The overexpression lines exhibited significant changes in the growth and development of rice, which indicated TaPARG-2D performs functions in plant morphology formation.
The genomic sequences of TaPARG-2A and TaPARG-2D and flanking regions were obtained by searching the A and D genomic databases. Genome-specific primer pairs for the 2A (Primer F2/R2) and 2D (Primer F3/R3) genes were designed based on polymorphisms in the genome sequences. The TaPARG-2A and TaPARG-2D genomic fragments were 3,388 bp ( Figure 6A) and 3,745 bp (Supplementary Figure S2), respectively. There was no nucleotide difference in the genomic region of TaPARG-2D among 34 accessions, but five variations (T/C, C/T, C/T, G/A, and C/T) were detected in the genomic region of TaPARG-2A ( Figure 6B); three haplotypes were identified and named Hap-2A-1, Hap-2A-2, and Hap-2A-3. Only the SNP (C/T) in the second exon resulted in an amino acid change (His/Tyr; Figure 6C). Functional markers TaPARM1 and TaPARM2 were developed based on the SNPs at 900 bp (C/T) and 2637 bp (C/T), respectively. Amplified fragment length polymorphisms of TaPARM1 and TaPARM2 are shown in Figures 6D,E.

Association Analysis of TaPARG-2A Haplotypes and Agronomic Traits
Population 1 (262 accessions) with an abundance of genetic diversity consisted of two sub-populations, comprising 126 and 136 accessions, based on genotyping results of 209 SSR markers . Significant associations were identified between TaPARG-2A haplotypes and agronomic traits, including PH, PL, TSL, ETN, and TKW in all 10 environments (year × site × water regimes; Table 1) in any agronomic trait. The results were in good agreement in both DS and WW conditions (Figure 7). Thus, Hap-2A-2 and Hap-2A-3 might be favorable haplotypes for improvement of plant architecture and yield in wheat.

The Effects of TaPARG-2A Haplotypes Were Confirmed in Population 3
Population 3 (348 modern cultivars) accessions from the Chinese wheat core collection were used to further evaluate haplotype differences. Functional analysis of TaPARG

Geographic Distribution of TaPARG-2A Haplotypes in Chinese Wheat Production Zones
Common wheat in China is grouped into three types according to time of sowing and vernalization requirement, namely, spring, FIGURE 2 | Subcellular localization of TaPARG-2D in wheat protoplasts (A) and tobacco leaf cells (B). The vector control (35S::GFP) and fusion protein construct 35S::TaPARG-2D::GFP were introduced into wheat protoplasts and tobacco leaf cells, respectively. For wheat protoplast transformation, GFP was detected in transformed wheat lines after 24 h with a laser scanning confocal microscope; for tobacco, GFP was detected in leaf cells 4 days after transformation. Images are in dark field (1, 5, 9, and 12), bright field (2, 6, 10, and 13), and combined (4, 8, 11, and 14). Chloroplasts are indicated by red auto-fluorescence (3 and 7). Scale bars: 20 µm for wheat protoplasts; 50 µm for tobacco leaf cells.

FIGURE 3 | Tissue expression patterns of TaPARGs in wheat.
Expression patterns at different developmental stages were detected by semi-quantitative RT-PCR. R, root; S, leaf sheath; L, leaf blade; I, internode; N, node; SP, spike; P, plumule; RB, root base. facultative, and winter types. The national wheat production area is divided into 10 major agro-ecological production zones based on ecological conditions, variety type and cropping season (Jin, 1983). Accessions in Population 2 (157 landraces) and Population 3 (348 modern cultivars) covering all 10 zones were used to evaluate the geographic distributions of TaPARG-2A haplotypes. Hap-2A-1 is the dominant haplotype in landraces in all zones, except Zone VI, i.e., the Northeastern Spring Wheat Zone, where Hap-2A-3 is the dominant haplotype ( Figure 9A). However, among modern cultivars, the combined frequency of the two favored haplotypes (Hap-2A-2 and Hap-2A-3) were higher in Zones I-IV (70-80%) than in Zones V-X (50-75%; Figure 9B). The zones with higher proportions of favored haplotypes are the major production regions, where the frequency of release of new varieties is higher than other regions.

Favored Haplotypes of TaPARG-2A Were Selected by Wheat Breeders
To assess changes in haplotype frequency over time 335 accessions in Population 3 with known dates of release were divided into six subgroups according to decade of release (Pre-1950s, 1960s, 1970s, 1980s, 1990s, and Post-2000. The proportion of Hap-2A-1 genotypes decreased significantly, from 81.8% in the pre-1950s cohort to 20.6% in the 1980s, correspondingly the frequencies of Hap-2A-2 and Hap-2A-3 increased from 18.2 to 79.4% (Figure 10A), and the PH also significantly reduced and TKW significantly increased ( Figure 10B). These data strongly suggest that the two favored haplotypes, Hap-2A-2 and Hap-2A-3, were apparently selected    ETN, effective tiller number; TKW, 1,000 kernel weight; n.s., not significant; * P < 0.05, * * P < 0.01, and * * * P < 0.005, respectively. The 10 environments were at Changping (CP) and Shunyi (SY) under drought-stressed (DS) and well-watered (WW) conditions in 2010-2012. in Chinese wheat breeding programs before the 1980s. The proportions of the favored haplotypes have maintained a relatively high and stable level since the 1980s, but the PH and TKW continued earlier trends implying that other genes affecting these traits were selected after the 1980s.

DISCUSSION
Genes TaPARG-2A and TaPARG-2D, containing two highly conserved AP2 domains, belong to the AP2 subfamily of AP2/EREBP transcription factors, that are involved in regulation of diverse processes of plant development and stress response. Both AP2 and ANT have important roles in reproductive organ development in angiosperm species (Shigyo and Ito, 2004). A series of AP2-like transcription factors have been characterized. AINTEGUMENTA-like6 (AIL6) is an important regulator of floral growth and patterning (Krizek, 2009). Another AP2-like gene Indeterminate spikelet 1 (IDS1), suppresses indeterminate growth within the spikelet meristem (Chuck et al., 1998). The functions of a gene family are generally conserved, but in some cases functions may be divergent. For example, TaPARGs, members of the AP2/EREBP transcription factor family, differ from other proteins in the family and are therefore classified   , spike number per plant (Wu et al., 2011(Wu et al., , 2012 and grain weight (Li et al., 2012). Two conditional QTLs affecting PH were separately flanked by Xwmc51-Xgwm68.2 and Xgwm448-Xgwm558, and the genetic distances between QTLs for PH and TaPARG-2A were 2.4 and 1.0 cM, respectively. The nearest markers of QTLs for spike number per plant were Xwmc401 and Xgwm448, and the genetic distances between QTLs for spike number per plant and TaPARG-2A were 4.0 and 1.0 cM, respectively. Two QTLs for grain weight were flanked by Xgwm68.2-Xgwm71.7 and Xgwm372-Xgwm448, and the genetic distances between the QTLs and TaPARG-2A were 1.7 and 1.0 cM, respectively. TaPARG-2A and three QTLs flanked by Xgwm372-Xgwm558 were located in the centromere region of chromosome 2A 1 , controlling the traits of grain weight, plant height, and spike number per plant, respectively, and all of them with 1.0 cM away from TaPARG-2A. By the method proposed by Visscher and Goddard (2004), TaPARG-2A fell into the 95% confidence intervals of all of these three QTLs. Gene expression studies showed that TaPARGs were expressed in multiple tissues at three growth  -1950s, 1960s, 1970s, 1980s, 1990s, and post-2000, respectively. (B) PH and TKW changes in 335 modern cultivars (Population 3) over decades. Error bars denote 1 SE. stages, suggesting roles in multiple developmental processes as a pleiotropic gene. Overexpression of TaPARG-2D in rice caused a range of changes, including reduced plant height, seed setting rate and 1,000 kernel weight, but an increase in tiller number per plant. Association analysis of different germplasm populations indicated that TaPARG-2A affects several agronomic traits (PH,  PL, TSL, ETN, and TKW). All these data show that TaPARGs are novel, pleiotropic genes in the AP2/EREBP gene family, and are involved in regulation of plant growth and development.
The present research suggests that TaPARG is a regulatory factor in plant growth and development by rice transgenic lines and association analysis in wheat. In the future, the gene expression and function should be further verified in same plant species to minimize the influence of genetic backgrounds.
Wheat was domesticated about 10,000 years ago (Gill et al., 2004). In two rounds of polyploidization, common wheat evolved with only part of the genetic diversity of its progenitors and subsequently generated new diversity by mutation tolerated by genetic duplication and redundancy (Dubcovsky and Dvorak, 2007). We identified TaPARG-2D was highly conserved and involved in regulating several agronomic traits, including plant height, tiller number per plant, seed setting rate, and 1,000 kernel weight. Three haplotypes of TaPARG-2A in hexaploid wheat based on five SNPs. According to an association analysis of TaPARG-2A genotype and agronomic traits, the three haplotypes were divided into favored (Hap-2A-2 and Hap-2A-3) and unfavored (Hap-2A-1) haplotypes for wheat improvement.
The amino acid residues (Arg147 and Gly148) in the GBD have backbone-binding function (Allen et al., 1998), and the equivalent functional loci in TaPARG-2A were Arg265 and Gly266. Five SNPs in TaDGRG-2A were identified among the accessions examined, but only one nucleotide difference (C/T) at 900 bp led to an amino acid change, i.e., His264 to Tyr264 (Figure 11). Histidine is an alkaline amino acid, whereas Tyrosine is a polar neutral amino acid. As the SNP change was closely adjacent to the DNA binding sites Arg265 and Gly266, the change likely influences interaction between TaPARG-2A and its target cis-element, and thus could result in a functional difference. The same single amino acid difference was present in two agronomically favored TaPARG-2A haploytpes, Hap-2A-2 and Hap-2A-3.
Wheat domestication is a process of accumulating superior genotypes and phenotypes. Increasing applications of nitrogenous fertilizer in traditional wheat management did not initially increase grain yield because the crop became too tall and lodged. The adoption of lodging-resistant wheat permitted higher cereal crop yields from the 1960s, and the "Green Revolution" was fueled in part by the introduction of semidwarf wheat varieties carrying Rht mutations (Borner et al., 1996). Almost all current commercial wheat varieties carry one of the Rht mutant alleles (Silverstone and Sun, 2000). Thus rational development and utilization of height-reducing genes made a huge contribution to increasing wheat production. Our research showed that the favored Hap-2A-2 and Hap-2A-3 haplotypes of TaPARG-2A were selected in the process of wheat improvement in China. From landraces to recent modern cultivars, the proportion of the favored haplotypes increased in stepwise manner. This increase was more apparent in the major wheat-production regions where there was a greater turnover of cultivars. The frequencies of the two favored haplotypes rose steadily from the 1950 to 1980s, during which time Chinese wheat varieties underwent very significant reductions in plant height and increases in grain yield. Given that the favored haplotypes of TaPARG-2A were apparently selected in the Chinese wheat breeding history, further selection of favored alleles might be useful for continued improvement of plant architecture and grain yield.
Plant breeding through phenotypic selection is a timeconsuming and relatively inefficient process (Gedye et al., 2012). Fortunately, the wheat community has made remarkable progress in developing molecular resources for breeding (Paux et al., 2012). Functional markers also called perfect or diagnostic markers have been developed for many traits including plant architecture, disease resistance, and grain quality (Gedye et al., 2012). In this study, after functional verification of TaPARG, a functional marker (TaPARM1) was developed to identify the favored Hap-2A-2 and Hap-2A-3 haplotypes. Currently, the combined frequency of the two favored haplotypes is about 80%. Utilization of the functional marker should contribute to further improvement of plant architecture and grain yield in wheat.