Genome-Wide Identification and Expression Profiling of the TCP Family Genes in Spike and Grain Development of Wheat (Triticum aestivum L.).

The TCP family genes are plant-specific transcription factors and play important roles in plant development. TCPs have been evolutionarily and functionally studied in several plants. Although common wheat (Triticum aestivum L.) is a major staple crop worldwide, no systematic analysis of TCPs in this important crop has been conducted. Here, we performed a genome-wide survey in wheat and found 66 TCP genes that belonged to 22 homoeologous groups. We then mapped these genes on wheat chromosomes and found that several TCP genes were duplicated in wheat including the ortholog of the maize TEOSINTE BRANCHED 1. Expression study using both RT-PCR and in situ hybridization assay showed that most wheat TCP genes were expressed throughout development of young spike and immature seed. Cis-acting element survey along promoter regions suggests that subfunctionalization may have occurred for homoeologous genes. Moreover, protein-protein interaction experiments of three TCP proteins showed that they can form either homodimers or heterodimers. Finally, we characterized two TaTCP9 mutants from tetraploid wheat. Each of these two mutant lines contained a premature stop codon in the A subgenome homoeolog that was dominantly expressed over the B subgenome homoeolog. We observed that mutation caused increased spike and grain lengths. Together, our analysis of the wheat TCP gene family provides a start point for further functional study of these important transcription factors in wheat.

Introduction

TCP genes as plant-specific transcription factors (TFs) are widely present in plants. Previous study showed that TCP genes are involved in a number of critical biological processes, including plant growth, development, and stress responses (Palatnik et al., 2003; Ori et al., 2007; Nag et al., 2009; Danisman, 2016). The name of TCP is derived from three important genes, i.e., the maize TEOSINTE BRANCHED 1 (TB1) that is a major determinant of strong apical dominance in domesticated maize (Doebley et al., 1997), the snapdragon CYCLOIDEA (CYC) which is involved in the control of floral bilateral symmetry (Luo et al., 1996), and the rice PROLIFERATING CELL FACTORS 1/2 (PCF1/2) that regulate PROLIFERATING CELL NUCLEAR ANTIGEN (PCNA) for DNA replication and repair (Kosugi and Ohashi, 1997). The TCP domain is composed of 59 amino acid residues, forming a basic helix-loop-helix (bHLH) type of DNA-binding domain non-canonical to regular bHLH TFs (Murre et al., 1989; Cubas et al., 1999). Based on its TCP domain, the members of the TCP family can be grouped into two subfamilies: class I (PCF or TCP-P class) and class II (TCP-C class) (Kosugi and Ohashi, 2002; Navaud et al., 2007; Martín-Trillo and Cubas, 2010). The difference between the two is a four-amino-acid deletion in the TCP domain in class I compared with those of class II. The members of class II are quite heterogeneous and can be further divided into two subclades: CIN and CYC/TB1 subclades (Martín-Trillo and Cubas, 2010). Outside the TCP domain, several class II members process an 18–20-residue arginine-rich motif called the R domain with an unknown function (Cubas et al., 1999).

Homology searches in complete genomes identified 31 TCPs in maize (Zea mays) and 24 in sorghum (Sorghum bicolor L.) and rice (Oryza sativa) (Mondragon-Palomino and Trontin, 2011). The maize TCP genes are distributed unevenly on 10 chromosomes. Based on phylogeny, maize TCP genes are categorized into nine subclasses and purifying selection is assumed to be responsible for maintaining their functions (Chai et al., 2017). Further analysis of maize TCP genes show that they express in stem and ears, suggesting their roles in the development of these two organs. On the other hand, a systematic analysis of sorghum TCP genes (SbTCPs) show that, except for SbTCP8, all are either intronless or contain introns in the untranslated regions. Seven pairs of paralogous TCP genes are identified from sorghum, five of which seem to predate rice-sorghum divergence, with diverged expression patterns. Five sorghum TCPs are considered to regulate plant morphology, whereas three genes have been identified as candidates for engineering abiotic stress tolerance (Francis et al., 2016).

In rice, a transcriptome profiling revealed differential accumulation of TCP genes during panicle initiation, organ development, and early anther development. Meanwhile, rice TCP genes are also differentially expressed in response to abiotic stress treatments such as cold, drought, and salt (Sharma et al., 2010). The rice ortholog of the maize TB1 gene is also expressed in axillary buds and appears to play similar functions as a negative regulator for lateral branching (Takeda et al., 2003). Additional functions were found for other rice TCP genes. The rice REP1 (RETARDED PALEA1) encodes a homolog of the CYC-type gene and during early flower development only expresses in palea primordium. At later floral development stages, REP1 is radially dispersed in stamens, the vascular bundles of lemma, and palea (Yuan et al., 2009). Another rice TCP gene, OsTCP19, appears to be an important node in cell signaling which crosslinks stress and developmental pathways (Mukhopadhyay and Tyagi, 2015). OsPCF2 may activate OsNHX1 gene expression, which responds to salt and PEG-induced drought stress, while OsPCF2 may be associated with the salt and drought stress tolerance (Almeida et al., 2017). Like in other species, some rice TCP genes are targets of microRNA319 (miR319). Transgenic rice plants overexpressing miR319 or down regulating TCP21 exhibited disease-like phenotypes and showed significantly higher susceptibility to RRSV virus in comparison with the wild-type plants. In rice, the induction of miR319 by RRSV infection suppresses Jasmonic acid (JA)-mediated defense and facilitates virus infection and symptom development by down regulating TCP21 expression level (Zhang et al., 2016).

Common wheat (Triticum aestivum L.) is a staple crop worldwide. Recent availability of wheat genomes allows detailed analysis of gene families in the wheat genome[1]. The wheat TB1 ortholog, TaTB1, has been confirmed to be involved in regulating wheat plant architecture (Liu et al., 2017) and to coordinate axillary spikelet formation during the vegetative to floral transition (Dixon et al., 2018). Except for TaTB1, no further analysis of TCP genes and the gene family analysis have been conducted so far. Here, we performed a comprehensive study of TCP genes in wheat. A total of 66 TCP genes belonging to 22 homoeologous groups were found in common wheat and were characterized in detail, including gene/protein architectures, domain conservation, physical properties, chromosomal location, phylogenetic relationship, and cis-elements of promoters. We also studied their tissue-specific expression patterns at the wheat heading date, their patterns of expression during young spike and immature seed development, and protein–protein interaction capabilities of some TCP proteins. Finally, two TaTCP9 mutants from durum wheat were characterized for its functions in spike and grain lengths. Together, our data provide valuable information for further investigation of the molecular functions of TCP genes in wheat which may be useful for wheat genetic improvement.

Materials and Methods

Identification of TCP Genes in the Wheat Genome

Wheat protein dataset was downloaded from https://urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Assemblies/v1.0/ and were searched using Arabidopsis and rice TCP protein sequences as queries using the BLASTP program (p-value <= 1e-5). After removing the redundant hits, the presence of TCP domains was verified by searching Pfam and SMART[2], with ambiguous sequences being manually confirmed using the InterProScan program. Finally, the rice TCP genes were followed to name the wheat TCP genes. Biochemical properties, such as the molecular weight (kDa) and isoelectric point (pI) of each protein, were determined using the Compute pI/Mw tool on the ExPASy website[3].

Phylogenetic Analysis and Gene Localition

Multiple sequence alignments were generated using Cluster X (v2.0) with default settings (Larkin et al., 2007). An unrooted phylogenetic tree was constructed using MEGA7.0 software using the neighbor joining (NJ) method, with a bootstrap number of 1000. Chromosomal locations of TCP genes were determined by searching the wheat genomic sequences. The software CIRCOS was used to draw the diagram showing TCP locations and homology relationships on wheat chromosomes that was downloaded from http://mapinspect.software.informer.com/.

Gene Structure Analysis and Identification of Conserved Motifs

The genomic and coding sequences of TCP genes, together with their exon/intron structures, were extracted from the general feature format (GFF3) file of wheat genome sequences. Additional conserved motifs were checked using the online Multiple Expectation Maximization for Motif Elicitation (MEME) program. The repetition was set as any number with an optimal width of 6–200 residues and the maximum number of motifs as 10.

Cis-Acting Element Analysis and miR159/miR319 Target Site Prediction

The putative promoter sequence, 1.5 kb upstream the transcription start site of each TaTCP gene (Zhao et al., 2014; Wang and Liu, 2015), was extracted from wheat genome sequences according to the general feature format (GFF3) file. Cis-acting elements were predicted at PlantCARE and those for plant growth, plant development, and phytohormone responses were counted. Using the Analysis of Motif Enrichment (AME) function in the MEME program[4], enrichment analysis was performed to identify regulatory elements within a collection of promoter sequences from all genes. A set of randomly interrupted promoters were used as a control. The motif with an adjusted Fisher’s test p-value less than 0.05 was considered to be a significantly enriched one.

To predict miR159 and miR319 target sites, full-length TaTCPs nucleotide sequences were analyzed using the psRNATarget online application. The cutoff of the maximum expectation is four.

Plant Growth and Tissue Collection

The wheat landrace Chinese Spring (CS) was grown in a farm in Beijing (39.97°N, 116.34°E) for two consecutive springs (2016 and 2017). The developmental stages of the young inflorescences were checked under a stereomicroscope (S8 APO, Leica Microsystems). According to Waddington et al. (1983), spikes at seven developing stages were collected. These stages included W2, W2.5, W3.5, W5.5, W6.5, W7.5, and W8.5. Various tissues including glume, lemma, palea, anther, pistil, and rachis of spike from the middle section of a spike were collected and employed to confirm tissue-specific expression genes at the W9.5 stage. We determined W10 stage as 0 day after pollination (DAP), and collected seeds (4–5 spikelets in the middle section of the spike) for 0, 2, 4, 6, 8, and 10 DAP. Three biological replicates were collected for all samples which were hand-dissected and immediately submerged in liquid nitrogen. Then they were stored at -80°C before RNA extraction. Simultaneously, samples of W3.5, W5.5, and W6.5 were prepared for in situ hybridization assays. Seeds for in situ hybridization were collected from CS plants that were grown in a glass house with day/night temperatures as 22°C/20°C and light illumination of 16 h/8 h light/dark. The relative humidity was 50%.

Mutant Identification

Two EMS mutant lines L2431 and L3090 derived from Kronos (Uauy et al., 2009), a desert durum wheat cultivar, were obtained from Drs. Daolin Fu and Jiajie Wu of Shandong Agriculture University, Shandong, China that were originally generated by Dr. Jorge Dubcovsky at UC Davis after Blast searches using wheat TCP sequences (Krasileva et al., 2017). DNA was extracted using the CTAB method as described elsewhere. DNA of mutation sites were PCR amplified and confirmed before further analysis of the mutants.

RNA Isolation and Real Time Quantitative RT-PCR Analysis

Total seed RNA was extracted using RNAprep Pure Plant Kit (Polysaccharides and Polyphenolics-rich, TianGend) and other total RNA was prepared using TRIzol reagent (Invitrogen). Quantitative real time (RT) PCR was carried out using an ABI PRISM 7300 RT-PCR system. The thermal cycling conditions were one cycle at 95°C for 5 min and 40 cycles of 95°C for 10 s, 55°C for 10 s, and 72°C for 30 s. The relative mRNA level of a gene was calculated as 2-DDCT value using the wheat GAPDH (glyceraldehyde-3-phosphate dehydrogenase) gene as an internal control for normalization. Each cDNA sample was tested with three replications. The primer 5 software used for primer design from conserved regions between A, B, and D subgenome with amplified fragments 200–350 bp (Supplementary Table ).

Yeast Two-Hybrid Assay

The coding sequences of wheat TCP genes were cloned into the bait vector pGBKT7 and/or the prey vector pGADT7 at different restriction sites for yeast two-hybrid assay. The two vectors were then co-transformed into Y2HGold yeast strain (Clontech). After being incubated on double dropout (DDO) medium at 30°C for 2–3 days, yeast cells were selected on DDO medium lacking Leu and Trp. Colonies grown up were further selected on quadruple dropouts (QDO) media lacking Leu, Trp, His and Ade, or on triple dropouts (TDO) media lacking Leu, Trp, and His with 1 mM 3-amino-1,2,4-triazole (3-AT) to identify colonies with positive protein–protein interactions. Empty vectors were used as controls.

In situ Hybridization

RNA in situ hybridization was performed according to Liu et al. (2013). Young spikes and immature grains were fixed overnight in formalin-acetic acid-alcohol at 4°C. The samples were then dehydrated through a standard series of ethanol solutions and were embedded in Paraplast Plus (Sigma-Aldrich). A microtome (RM2235, Leica Microsystems) was used to cut tissues into 8 mm sections. Gene-specific regions were amplified and used to synthesizing Digoxigenin-labeled sense and antisense RNA probes using a DIG northern Starter Kit (Roche), according to the manufacturer’s instructions. The sequences of primers are listed in Supplementary Table .

Results

Identification, Phylogeny, and Classification of TCP Genes in Wheat

To identify wheat TCP genes, we used rice and Arabidopsis TCP protein sequences as queries and searched the wheat protein dataset at IWGSC[5] using Blastp. A total of 66 TCP genes were retrieved according to their similarity to query TCP genes and the presence of TCP domains and conserved motifs. The serial numbers of wheat TCP genes were named according to their best rice homologs (Table ). All wheat TCP genes had three homoeologs.

A phylogenetic tree was constructed using neighboring joining (NJ) method so as to distinguish the evolutionary relationship of wheat TCPs. Similar to that in other plants, wheat TCPs can be categorized into two classes: class I (PCF) and class II (CIN and CYC/TB1) (Figure ). Among them, 33 members (TaPCF1, 2, 3; TaTCP6, 7, 9, 17, 19, 25, 28, and 29) were in class I and the remaining fell into class II, including 12 CYC/TB1-type (TaTB1-1, TaTB1-2, TaTCP22, and 24) and 21 CIN-type (TaPCF5, 6, 7, 8; TaTCP5, 18, 21). This is comparable to those in rice where there are 13 class I TCPs and 12 class II TCPs, including three CYC/TB1 type TCPs and nine CIN type TCPs. It is interesting that the wheat TB1 locus was duplicated, one on chromosome 4 and the other on chromosome 5 and each having three homoeologs (Table and Figure ). The one on chromosome 4 has been shown to regulate inflorescence architecture and development in wheat (Dixon et al., 2018), while there is no functional report on the one on chromosome 5. To confirm this observation, we obtained all TB1 orthologs and paralogs from wheat donor species including diploid (Aegilops tauschii and Triticum urartu), tetraploid (Triticum turgidum), as well as some other monocots including barley, rice, maize, Brachypodium, and sorghum (Figure and Supplementary Table ). Phylogenetic analysis showed that wheat related species all had two TB1 homologs. Maize also had two TB1 homologs, but seemed to be arising by segmental duplication. In addition, all remaining species contained one TB1 gene including rice, Brachypodium, and sorghum. Then we calculated the Ks and Ka/Ks ratio for each duplicated TB1 gene pairs (Table ). The results showed that Ka/Ks of all five pairs were <1, indicating the effect of strong purifying selection and the slow evolution rate in TB1 protein sequences in wheat. The duplication of the five paralogous gene pairs was estimated to have occurred 0.62–0.84 million years ago (mya). We also found that there was no wheat ortholog for rice CIN-type TCPs OsTCP10 and 27 and for rice PCF-type TCPs OsTCP11, 12, and 20 (Figure ). Conversely, wheat TaTCP29 (TaTCP29-A, TaTCP29-B, and TaTCP29-D) had no ortholog in rice. To get a better understanding of the diversification of the TaTCP genes, we analyzed the exon/intron organization of TaTCPs and found that nearly all genes had their introns located at the 5′ or 3′ UTR regions, except for TaPCF2-A and TaTCP7-B which had their introns located in the middle of coding regions (Figure ).

Phylogenetic and gene structure analyses of wheat TCPs. (A) A Neighbor Joining tree of wheat TCP proteins. (B) Wheat TCP gene structures. UTR, CDS, and intron are represented with green boxes, yellow boxes, and black lines respectively.

Phylogenetic tree of TB1 orthologs and paralogs from wheat related species Aegilops tauschii (Aet, DD), Triticum urartu (Tu, AA), tetraploid Triticum turgidum (Tt, AABB), barley (Hv), and other monocot species Brachypodium (Bd), Zea mays (Zm), Oryza sativa (Os), and Sorghum (Sb). A, B, and D indicate wheat three subgenomes.

Genomic Organization and Duplication of Wheat TCP Proteins

Wheat TCP genes were found to be located on all chromosomes, with a maximum of seven genes on chromosome 5 and only one on chromosome 1 and 2 (Figure ). The A homoeolog of TaPCF2 was translocated from 5AL to 4AL which was consistent to previous study (Devos et al., 1995). The A homoeolog of TaTB1-1 and TaPCF3 were translocated from 4AS to 4AL, which was confirmed by check locations 100 Mb up and downstream regions of these two genes. Our pipeline also detected the 7BS and 4AL translocation that has been reported previously (Berkman et al., 2012). Interestingly, as mentioned earlier, TaTB1 was found on two loci, one on chromosome 4 and the other on chromosome 5. This is in contrast with those in rice and maize where only one TB1 gene was reported (Takeda et al., 2003; Clark et al., 2006).

Chromosomal localizations of TCP genes in wheat. TCP genes were mapped on wheat chromosomes by using genomic information available in (https://urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Annotations/v1.0/). Red lines represent the synteny of triplets. Different colors represent different linkage groups. Blue, the A subgenome; light purple, the B subgenome; and Purple, the D subgenome.

Conserved Motifs and Domains on TaTCP Proteins

To identify conserved motif and evaluate structural divergence, TaTCP proteins were further analyzed at the MEME web server. A total of 10 conserved motifs (motif 1–10) were identified (Supplementary Figure ). All but five (TaPCF1-A, D; TaTCP29-A, B, D) had conserved motif 1 and motif 2 that constituted the typical TCP domain. Similar to those reported by Chai et al. (2017), TCP proteins in the same subgroup share similar conserved motifs that were clearly distinguishable from those in other subgroups. CIN-type TCPs formed three subgroups: subgroup 1 had eight motifs (1, 2, 3, 4, 5, 6, 7, and 8); subgroup 2 had four motifs (1, 2, 4, and 8); and subgroup 3 had eight motifs (1, 2, 4, 5, 8, 9, and 10) (Figure and Supplementary Figure ). CYC/TB1-type TCPs had two subgroups: subgroup 1 conferred five motifs (1, 2, 3, 4, and 8) and subgroup 2 had four motifs (1, 2, 8, and 9), while PCF-type had five subgroups with subgroup 1 having three motifs (1, 2, and 3), subgroup 2 having seven motifs (1, 2, 3, 5, 6, 9, and 10), subgroup 3 having six motifs (1, 2, 3, 5, 6, and 8), subgroup 4 having seven motifs (1, 2, 3, 4, 5, 6, and 8), and subgroup 5 having seven motifs (1, 2, 3, 5, 7, 8, and 9).

To identify conserved domains of TCP proteins and to further classify TCP proteins, we searched the TCP domain using Scanprosite and found a total of four conserved motifs – Basic, Helix I, Loop, and Helix II (Figure ). In most cases, TCP proteins in the same subclade shared similar motif compositions. A four-amino-acid deletion was found to distinguish the class I TCPs from those of class II in their basic domain (Figure ).

Multiple sequence alignment of wheat TCP transcription factors. (A) Alignment of the TCP domain of wheat TCP proteins. Black boxes highlight residues conserved in all three TCP classes; Red, residues conserved in the PCF class; Green, residues conserved in the CIN class; and Purple, residues conserved in the CYC/TB1 class. Domains of the Basic, Helix I, Loop, and Helix II regions are underlined on the top of the alignment. (B) Alignment of the R domains.

A wheat TCP phylogenetic tree was built using full length protein sequences. The structure was similar to those of Arabidopsis and rice where all the proteins can be classified into two major classes, with those having a four-amino-acid deletion in the basic region of the TCP domain as class I TCP and those with TCP domains as class II which were further divided into CIN- and CYC/TB1-subclades (Yao et al., 2007). The wheat CYC/TB1 subclade contained four TaTCP genes (TaTB1-1, TaTB1-2, TaTCP22, and TaTCP24), and the CIN subclade contained seven (TaPCF5, TaPCF6, TaPCF7, TaPCF8, TaTCP5, TaTCP18, and TaTCP21; Figure ). The R domain was found in a subset of class II proteins (TaTB1-1, TaTB1-2, and TaTCP21) (Figure ).

Wheat TCP Genes With microRNA Target Sites

In Arabidopsis, CIN-TCP genes are post-transcriptionally regulated by miR319a (Palatnik et al., 2003). In case of the CIN-type TCPs, increased miR319a activity (e.g., jaw-D mutants and miR319a-overexpressing lines) causes simultaneous down regulation of five TCP genes (Nicolas and Cubas, 2016). We found in wheat that TaPCF5, 6, 8, and TaTCP21 contained sequences well matched with miR319a and might be the targets of microRNAs (Figure ). This is consistent with those in rice where five genes OsPCF5, 6, 8, OsTCP20, 21, and 27 are targeted by miR319. Moreover, four wheat TCPs: TaPCF5, 6, 8, and TaTCP21, may be the targets of the second microRNA tae-miR159a (Figure ).

Alignment of putative target areas for miR319 and miR159. Mismatches and G-U wobbles were represented by blue and purple, respectively.

Identification of Cis-Acting Elements in the Promoter of TaTCP Genes

The distribution of different cis-acting elements in the promoter of a gene may indicate the difference in its function and regulation. We extracted 1.5 kb of genomic sequence upstream the translation start site of each TCP gene and searched for cis-acting elements in the PlantCARE database (Lescot et al., 2002). Cis-elements responsible for plant growth and development and phytohormone responses were identified (Figure ). Two motifs, Skn-1 and GCN4, are involved in endosperm expression (Wang et al., 1999). Skn-1-motif was found in most TaTCP genes (at least in one homoeolog) except for four PCF-type genes (TaTCP9, 19, 25, and 28). The meristematic expression and specific activation elements CAT-box and CCGTCC-box were also found in most TaTCP genes. The PCF-type gene TaTCP7 and the CIN-type genes TaTCP5, TaTCP18, and TaPCF7 contained Skn-1 motifs, but no CCGTCC-box element was found in any of its homoeolog. All three homoeologs of CIN-type genes (TaTCP21, TaPCF6, 7, and 8), the PCF-type TaTCP29, and the CYC-type TaTB1-1 and TaTCP24 were found to bear circadian control element (Anderson et al., 1994). Wheat TCP gene promoters also possessed other cis-elements, such as those specific to seed and shoot development (e.g., RY and as-2-box) (Bobb et al., 1997) and those for zein metabolism regulation element (e.g., O2 site).

Cis-acting elements on promoters of wheat TCP genes.

For hormone-related cis-acting elements, the MeJA-responsive elements CGTCA and TGACG (Rouster et al., 1997) were most frequently identified at the wheat TCP gene promoters. The ABA-responsive element (ABRE; Shen and Ho, 1995) was also found in most of TaTCP genes. For TaTCP7, six and five cis-elements were identified on its A and B homoeolog promoters, respectively, but none on that of the D homoeolog. A number of other hormone related cis-elements, such as gibberellin responsive elements GARE and P-box and auxin-responsive elements TGA, AUXRR core, and TGA box (Kim et al., 1992; Ulmasov et al., 1997; Washida et al., 1999) were also present in promoters of some wheat TCP genes.

Since cis-acting elements are usually very short and not very well defined, to avoid the likelihood of detecting them by chance, we performed cis-acting element enrich analysis and confirmed that cis-elements responsible for plant growth and development were indeed enriched (Table ).

Spatial and Temporal Expression Patterns of Wheat TCP Genes

Previously, we generated RNA-seq data of wheat young inflorescence at double ridge (DR, Waddington scale W2), floret meristem (FM, W2.5 ∼ 3), anther primordium (AM, W3.5 ∼ 4), and tetrad (TS, W7.5) stages (Feng et al., 2017). We found that 31 homoeologs of 15 TCP genes were present in this dataset (Supplementary Table ). As shown in Figure and Supplementary Table , TCP genes with the highest expression level in TS stage were five class II genes, namely TaPCF5, 6, and TaTCP5, 18, and 21, while the class II gene TaPCF8 was most highly expressed at the stage W3.5 (Figure ), suggesting that TaTCP genes may be involved in floral organ development. For class II CYC/TB1 genes, TaTCP24 and TaTB1-1 showed the lowest expression level at the W7.5 stage, while TaTCP22 showed similar expression levels from W2.5 to W7.5. For class I PCF genes, TaPCF3, TaTCP9, and TaTCP17 were expressed with the highest level at W7.5, while TaTCP6 was expressed with the highest expression level at DR. Interestingly, TaTCP7 and TaTCP28 were constitutively expressed in all the four stages. We studied the expression patterns of wheat TCP genes at more time points using RT-PCR including W5.5, W6.5, and W8.5. The W3.5–W7.5 stages cover the key developmental stages of floral organ initiation and differentiation (Sreenivasulu and Schnurbusch, 2012). The six class II CIN-type TCPs (TaPCF5, 6, 8, TaTCP5, 18, 21) were expressed with high levels from W3.5 to W7.5 (Figure ), suggesting that these genes may play important roles in floral organ development. For class II CYC/TB1-type genes, TaTCP22 showed the highest expression level at W6.5 while TaTB1-1 was highly expressed at W2 (DR) which was consistent with that of RNA-seq. For six class I PCF-type genes, four (TaTCP6, 7, 17, and 28) were consistent with their RNA-seq patterns and two of them (TaPCF3 and TaTCP9) were different, probably caused by the changing of reference genomes for these two detecting methods. For instance, in the old version (Ensemble Triticum_aestivum.IWGSC2.25), only two homoeologs, A and D, were identified for TaTCP9, while in the new version (IWGSC RefSeq v1.0), all three homoeologs, A, B, and D, were identified. The change in homoeolog copies may cause inaccurate expression patterns in RNA-seq quantification, different to that of RT-PCR analysis which used primers consensus to all three homoeologs. From W2 to W7.5, the expression level of TaTCP17 increased significantly (Figure ). Such a pattern may indicate its importance during wheat inflorescence development and should be further studied.

Expression patterns wheat TCP genes during young spike development by RNA-seq. Developmental stages are indicated using Waddington (W) scales (Waddington et al., 1983). W2, double ridge (DR); W2.5, floral meristem stage (FM), W3.5, anther meristem stage (AM); and W7.5, tetrad stage (TS). Y axis is the RPKM (Reads per kilobase per million mapped reads) value; gene names are shown in Supplementary Table . Each point of the RPKM represents the average of two biological replicates. Error bars indicate SD.

TaTCP gene expression patterns during young spike development by RT-PCR. Developing young spike stages were determined using Waddington (W) scales (Waddington et al., 1983). Three biological duplicates were performed for each gene; Y axis represents the relative expression value to GAPDH gene, the internal control. Asterisks indicate significant difference relative to W2, as determined by Student’s t-test (∗p < 0.05, ∗∗p < 0.01).

Since most TaTCP genes were expressed at higher levels toward the later stages of floral organ development, we studied their organ-specific expression patterns at the heading date stage (W9.5). A total of six tissues were used including rachis, glume, lemma, palea, pistil, and stamen (Figure ). Two class I genes (TaTCP7, 19) were highly expressed in rachis. Seven TCPs including four class I TCPs (TaPCF3, TaTCP9, and 17, 29) and three class II CIN-type ones (TaPCF5, 6, and 8) were highly expressed in pistil. To determine whether wheat TCP genes also expressed in immature grain, we detected the expression patterns of the above seven pistil-specific genes in different grain developing stages. Four of them (TaPCF5, 6, 8, and 3) were down regulated after pollination which were further down regulated from 2 to 4 DAP (Figure ). Conversely, TaTCP9 and TaTCP17 were most highly expressed at 2 DAP, while TaTCP29 at 2-4 DAP, suggesting that these three genes may play a role in early grain development.

Tissue-specific expressions of wheat TCP genes. Tissues collected from plants at W9.5 and expression levels were detected by qRT-PCR. Three biological duplicates; Y axis represents the relative expression value to GAPDH gene, the internal expression control.

TaTCP gene expression patterns during grain development. Expression levels were detected by qRT-PCR with immature grains at six stages: 0, 2, 4, 6, 8, and 10 days after pollination (DAP). Three biological duplicates were performed; Y axis represents the relative expression value to GAPDH gene, the internal control. Asterisks indicate significant difference relative to W2, as determined by Student’s t-test (∗p < 0.05, ∗∗p < 0.01).

In situ Hybridization Assay of Wheat TCP Genes

To investigate whether TaTCP genes were expressed in specific tissues at early stages of spike development, we conducted in situ hybridization of four selected genes with relatively high expression level at W3.5, W5.5, and W6.5. We found that the four CIN-type genes (TaTCP18, 21 and TaPCF5, TaPCF6) were expressed from floret initiation to floral organ maturation and did not display specific expression domains along the spike. As shown in Figure , the four genes were expressed in multiple floral organs including lemma, palea, pistil, and stamen. The ubiquitous expression of these genes indicates that they may play multiple roles in maintaining cell specificity throughout early inflorescence development.

In situ hybridization assays of four TaTCP genes. (A–D) TaPCF5; (E–H) TaPCF6; (I–L) TaTCP21; and (M–P) TaTCP18. Spike developmental stages were W3.5: A, E, I, M; W5.5: B, F, G, N; W6.5: C, G, K, O. Sense probes were used for D, H, L, and P. Glu, glume; Lem, lemma; Pal, palea; Pis, pistil; Sta, stamen; and bar = 500 μm.

Protein–Protein Interaction Patterns Among Wheat TaTCPs

TCP proteins function by forming protein complexes (Kosugi and Ohashi, 2002; Parapunova et al., 2014). We analyzed the interaction capabilities of several TaTCP proteins not only with relative high expression level but also with the same expression pattern, specifically TaTCP18-D, TaTCP21-D and TaPCF3-A, using yeast two-hybrid assays. These TCP genes were expressed with similar patterns during grain development (Figure ). As shown in Figure and Table , class II CIN-type TCPs TaTCP18-D and TaTCP21-D formed both homodimers and heterodimers, while the class I PCF-type TaPCF3-A formed heterodimers with TaTCP18-D or TaTCP21-D but no homodimers were formed. These data indicate that wheat TCP proteins may carry out their functions by interacting with each other to form protein complexes, similar to that in rice and Arabidopsis.

Pair-wise protein-protein interactions among TaTCP18D, TaTCP21D, and TaPCF3A. Coding sequences of TaTCP genes were cloned to pGADT7 (AD) and pGBKT7 (BK) vectors. Interactions among the TaTCP proteins were analyzed by yeast two-hybrid assay. Transformants were assayed for growth on QDO nutritional selection medium.

Characterization of Tetraploid Mutants Carrying Loss-of-Function Alleles of TaTCP9

In rice, OsTCP19 appears to be an important node in cell signaling which crosslinks stress and developmental pathways (Mukhopadhyay and Tyagi, 2015). OsTCP9 and OsTCP19 belong to the same phylogeny branch (Figure ), but no functional study of OsTCP9 and its homologs in any other species is available. Here, we firstly studied the tissue specific expression patterns of TaTCP9 and found that it was highly expressed in pistil at the W9.5 stage (Figure ). A time course study in immature grains showed that its expression peaked at 2 DAP and decreased at later stages of grain development (Figure ), a pattern similar to the rice TGW6 gene (an IAA-glucose hydrolase) which enhances rice grain weight and increases yield (Ishimaru et al., 2013). To identify the precise expression domain of TaTCP9, in situ hybridization assays were performed. The results showed that TaCP9 began to express at early stages of grain development (2, 4, and 6 DAP) and was mainly expressed in endosperm transfer cells (ETC) and nucellar projection transfer cells (NPTCs) (Figures ). These cells are important because the first layer of ETC confers aleurone cell features and the second and third layers of ETC accumulate starch granules and protein bodies (Zheng and Wang, 2011) and actively transport sucrose from photosynthetic tissues to endosperm and embryo (Pfeifer et al., 2014).

Functional characterization of the wheat TaTCP9 gene. (A–D) In situ hybridization of TaTCP9 in immature seeds at 2 DAP (A), 4 DAP (B), and 6 DAP (C) in Chinese Spring. (D) sense control. Pe, percarp; VB, vascular bundle; NP, nucellar projection; EC, endosperm cavity; ETC, endosperm transfer cells; and NPTC, nucellar projection transfer cells. (E,F) Expression pattern of TaTCP9 during grain development in hexaploid wheat (E) and tetraploid wheat (F) Each point of the FPKM represents the average of three biological replicates. (G) the mutation sites on the TaTCP9-A gene that gave premature stop codons in two Kronos mutant lines L2431 and L3090. (H) grain morphology of the mutants L2431 and L3090 relative to the wild type Kronos. (I) Increased grain length in L2431 and L3090. Error bars indicate SD. A total of 23 spikes from six L2341 plants and 12 spikes from four L3090 plants were measured. Student’s t-test, ∗p < 0.05, ∗∗p < 0.01. Bars in A–D indicate 500 μm; Bar in H indicates 1 cm.

We also determined the expression patterns of TaTCP9 in tetraploid wheat (Triticum turgidum ssp. durum) because there are mutant lines that are generated using durum variety “Kronos” and can be used to genetically study the function of TaTCP9 (Uauy et al., 2009; Henry et al., 2014; Krasileva et al., 2017). TaTCP9 was highly expressed in Kronos grains of 2, 4, and 6 DAP (Figures ). Such a pattern is similar to that in common wheat, indicating that its functions in tetraploid and hexaploid wheat may be conserved. For the two homoeologs in Kronos, TaTCP9-A was expressed higher than TaTCP9-B (Figure ). In CS, the A genome homoeolog TaCTP9-A also expressed higher than the other two homoeologs TaTCP9-B and TaTCP9-D (Figure ). We searched the mutant sequence dataset using genomic sequence of TaTCP9 and two mutant lines L2341 and L3090 were found to contain C811T and C742T mutations, respectively, on their TaTCP9-A homoeolog leading to premature stop codons (Figure and Supplementary Figure ). Although these tetraploid mutants have been selected by their phenotypes for several generations (Krasileva et al., 2017), there are still additional mutation sites in them.

Morphologically, these two mutant lines showed significantly longer grains than those of the wild type (Figures ). In addition, grain width was wider in L2431 (Supplementary Figures ), resulting in higher 1,000-grain weight (Supplementary Figure ). Length of spike was also longer for the two mutant lines (Figures ). We found that the increased spike length was caused by both increased spikelet number and increased rachis internode length (Figures ).

Spike phenotypes of TaTCP9 mutants in tetraploid Kronos lines. (A) spike morphology of the mutants L2431 and L3090 relative to the wild type Kronos. (B) Increased internode numbers in L2431 and L3090. (C,D) Statistic analysis of spike length and spikelet numbers in L2431 and L3090. (E) Comparison of rachis internode lengths in the wild type and the two putative mutant lines L2431 and L3090. Error bars indicate standard deviation (SD). For spike length and spikelet number, 20 spikes from 20 L2341 to L3090 plants were measured. Student’s t-test, ∗p < 0.05, ∗∗p < 0.01. Bars in (A,B) indicate 1 cm.

Previous work showed that cell size and number increase most significantly in spikes from W2.5 to W3.5 (Feng et al., 2017). We then looked for all mutated genes (those with stop gain and altered splicing) for their expression patterns between these two stages (Supplementary Datasets –), particularly those that may be relevent to this process, such as those fell in the MapMan bins of cell cycle and devision, organization (RNA processing), recognition, and transcription, hormone, development (Supplementary Tables –). As shown in Supplementary Figures –, only two genes TaTCP9 (Traes_6AL_DA27ABCA61) and Traes_1BS_8A19C460B (an Auxin-responsive Aux/IAA gene) displayed a significant expression increase between W2.5 to W3.5, with the former having stop gain mutations in both L2431 and L3090 and the latter having such a mutation only in L2431. These data suggest that wheat TaTCP9 may indeed play a role in wheat grain development, although additional genetic study of these mutant lines regarding the function of TaTCP9 should be further performed.

Discussion

TCPs govern key plant developmental processes and have profound effects on meristem patterning and organ differentiation, often by mediating hormone mechanisms (Nicolas and Cubas, 2016). The polyploid nature of common wheat renders expansion of the TCP family in this important crop. Subsequent subfunctionalization of these homoeologs and more versatile gene dosages may contribute to the enhanced adaptation of wheat. In wheat, the general organization of the TCP family remained well conserved with significantly more members in the class I subfamily than those in the class II (Navaud et al., 2007).

The Duplication of Wheat TaTB1, an Important Tillering Gene

The importance of TCP genes were firstly demonstrated by the maize TB1 gene for its role in maize domestication as a negative regulator for the growth of axillary buds and hence branching (Doebley et al., 1995). The process may represent a common mechanism for post-embryonic development of the primary shoot architecture. Phytohormones, especially auxin and cytokinin, play pivotal roles during lateral branch development. The rice TB1 ortholog OsTB1 (or FINE CULM1, FC1) and the Arabidopsis BRANCHED 1 (BRC1) appear to play similar roles by negatively regulating lateral branching (Doebley et al., 1997; Takeda et al., 2003; Aguilar-Martinez et al., 2007). Moreover, OsTB1 is directly regulated by IPA1/OsSPL14 to suppress tillering in rice (Lu et al., 2013).

For wheat, in addition to the increased copy number of TB1 homoeologs due to polyploidization, it is duplicated within each subgenome, with a second copy found on chromosome 5, in addition to the orthologous one on chromosome 4. This is in contrast to rice, Brachypodium, and Arabidopsis in which only one TB1 gene is present. The orthologous chromosome 4 TaTB1-1 was shown recently as a key regulator for wheat inflorescence and plant architecture. It interacts with the Flowering Locus T1 (FT1) to regulate spikelet development (Dixon et al., 2018). The second copy of TaTB1 (TaTB1-2) had a different expression pattern from TaTB1-1 with the former having highest expression level at W6.5 when most floral organs were under development and the latter having its maximum expression level at the DR stage (W2.0). The subfunctionalization of TaTB1-2 is evident and should be further studied.

Nearly All Wheat TCP Genes Are Triplet With Three Homoeologs

Common wheat is a hexaploid (AABBDD, 2n = 6x = 42) that arose 7000–12,000 years ago following the hybridization of the early-domesticated allotetraploid T. turgidum ssp. dicoccon (AABB) and the diploid goat grass Aegilops tauschii (DD, 2n = 2× = 14) (Kihara, 1944; McFadden and Sears, 1946; Huang et al., 2002; Matsuoka, 2011). Thus, most genes are expected to have three homoeoalleles although in many cases the three homoeoalleles are not equally transcribed (Shitsukawa et al., 2007; Geng et al., 2013; Hu et al., 2013). Despite this, among the chromosome survey sequence (CSS) less than 23.16% genes (6576 triplets totaling 19,728 genes out of 85,173 high-confidence gene models) were identified as triplets, loci with three homoeologs, using a reciprocal best-hits approach with BLASTP e-value <= 1e-5 and identity > = 90% (Pfeifer et al., 2014).

Interestingly, all wheat TCP genes conferred three homoeologs indicating that they were preferentially retained after two rounds of polyploidization. Cis-element analysis showed clear differences between many homoeologs, consistent with their differential expression patterns in RNA-seq data. The maintenance of TCP homoeologs indicate they are under purifying selection while the divergence in cis-element patterns suggests on-going subfunctionalization of wheat homoeologs.

TCP Genes are Ubiquitously Expressed During Early Spike Development in Wheat

The profound functions of TCP genes from meristem maintenance to organ development require the expression of these genes in various tissues and organs. In this study, we focused on the expression patterns of wheat TCPs in spike and grain development. The relative low expression levels of wheat TCP genes are consistent with those identified in the orchid (De Paolo et al., 2015). For some wheat TCP genes, no signal can be detected by in situ hybridization. However, the expression range of wheat TCPs covers the whole processes of young spike and immature grain development. Such patterns indicate wide functions of this gene family in multiple aspects of wheat development.

Despite their versatile functions in plant development, to our knowledge, TCPs are rarely reported to be involved in grain development. We present preliminary data that wheat TaTCP9 gene might be involved in spike and grain development. TaTCP9 was expressed highly in pistil at the W9.5 stage and in immature grains at 2 DAP. In situ hybridization showed that TaTCP9 expressed in ETC and NPTCs. In Arabidopsis, as the key step to achieve final size and function for cells, endoreduplication is prevalent during plant development. AtTCP15, for instance, plays an important role in regulating endoreduplication during Arabidopsis development (Li et al., 2012). Moreover, mutants in two related class I TCP genes display a range of growth-related phenotypes, consistent with their dynamic expression patterns. The two genes influence plant stature by promoting cell division in young internodes (Kieffer et al., 2011). In our study, the increased spikelet number as well as increased rachis internode length may be accounted for in mutant plants. It is probable that the mutation in TaTCP9 might increase both cell number and cell size. Together with the phenotypes of the tetraploid mutants, we propose that TaTCP9 is highly possible to be involved in grain development and hence final grain yield, although further study is needed to reach a final conclusion.

Like TCPs in other plants, most wheat TCP genes expressed in multiple tissues and developmental stages, indicating their important roles in wheat development. Further molecular study of these genes should reveal more functional mechanisms for these genes. Since wheat is a polyploid, different to other major crops, study of wheat TCPs may help further understanding dosage-dependent working modes in this polyploid plant and may contribute to genetic engineering for wheat yield improvement.

Author Contributions

AL and LM planned and designed the research. JZ, YL, SG, GaS, MJ, FW, GuS, NF, XK, and LC performed the experiments. ZZ, JZ, YL, and JG analyzed the data. AL, LM, and JZ wrote the article, with contributions from all the authors.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewer NP and the handling Editor declared their shared affiliation.

Table 1

Wheat TCP family genes.

Gene name	Gene ID	ORF length	No. intron	pI	Mw	Chromosome location	Subclade
TaPCF5-A	TraesCS3A01G140100.1	1373	2	6.44	47915.75	chr3A:119137149-119140425	CIN
TaPCF5-B	TraesCS3B01G164500.1	1352	2	6.37	47007.34	chr3B:162222339-162224970	CIN
TaPCF5-D	TraesCS3D01G146900.1	1364	2	6.41	47521.58	chr3D:111206159-111208623	CIN
TaPCF6-A	TraesCS5A01G423900.1	1217	1	9.59	38470.11	chr5A:609619180-609621983	CIN
TaPCF6-B	TraesCS5B01G426200.1	1214	1	9.39	41923.59	chr5B:601864073-601867389	CIN
TaPCF6-D	TraesCS5D01G432500.1	1205	1	9.39	41382.26	chr5D:488677442-488679253	CIN
TaPCF7-A	TraesCS3A01G422500LC.1	1037	intronless	5.62	36063.32	chr3A:524756895-524757932	CIN
TaPCF7-B	TraesCS3B01G503200LC.1	1037	intronless	5.72	36269.33	chr3B:529404878-529405915	CIN
TaPCF7-D	TraesCS3D01G384600LC.1	1043	intronless	5.52	36205.32	chr3D:403111095-403112318	CIN
TaPCF8-A	TraesCS5A01G031900.1	893	2	6.76	32293.85	chr5A:29512544-29516131	CIN
TaPCF8-B	TraesCS5B01G032000.1	902	intronless	6.75	32434.76	chr5B:35188411-35189313	CIN
TaPCF8-D	TraesCS5D01G040100.1	887	2	6.7	31938.64	chr5D:41117072-41120665	CIN
TaTCP5-A	TraesCS3A01G302900.1	836	1	6.95	30140.75	chr3A:536895105-536897646	CIN
TaTCP5-B	TraesCS3B01G333400.1	842	1	7.3	30335.92	chr3B:539942261-539944819	CIN
TaTCP5-D	TraesCS3D01G298900.1	827	1	7.3	29822.66	chr3D:413676020-413678583	CIN
TaTCP18-A	TraesCS1A01G332600.1	920	1	6.72	32919.34	chr1A:520938291-520940709	CIN
TaTCP18-B	TraesCS1B01G346300.1	926	1	7.1	33085.47	chr1B:575440078-575442514	CIN
TaTCP18-D	TraesCS1D01G335300.1	941	intronless	6.41	33678.57	chr1D:425299621-425300562	CIN
TaTCP21-A	TraesCS2A01G251100.1	1646	1	10.02	59214.07	chr2A:379150934-379159267	CIN
TaTCP21-B	TraesCS2B01G265800.1	1274	1	9.18	45596.65	chr2B:357918017-357924451	CIN
TaTCP21-D	TraesCS2D01G252100.1	1250	1	9.42	44631.32	chr2D:302247032-302253739	CIN
TaTB1-1-A	TraesCS4A01G271300.1	1058	intronless	7.15	37911.58	chr4A:582839670-582841244	CYC/TB1
TaTB1-1-B	TraesCS4B01G042700.1	1064	intronless	7.72	38223.75	chr4B:30362277-30363341	CYC/TB1
TaTB1-1-D	TraesCS4D01G040100.1	1058	1	7.71	37835.57	chr4D:18463838-18472387	CYC/TB1
TaTB1-2-A	TraesCS5A01G001900.1	821	intronless	9.46	30056	chr5A:1346522-1347511	CYC/TB1
TaTB1-2-B	TraesCS5B01G002300.1	887	intronless	9.42	32627.27	chr5B:3474782-3475968	CYC/TB1
TaTB1-2-D	TraesCS5D01G002300.1	794	intronless	9.83	29191.66	chr5D:2511000-2512022	CYC/TB1
TaTCP22-A	TraesCS2A01G018000.1	614	intronless	5.42	22126.96	chr2A:8543341-8543955	CYC/TB1
TaTCP22-B	TraesCS2B01G025700.1	614	intronless	5.09	21865.79	chr2B:11836978-11837592	CYC/TB1
TaTCP22-D	TraesCS2D01G019300.1	710	intronless	5.46	25633.9	chr2D:9079476-9080529	CYC/TB1
TaTCP24-A	TraesCS5A01G207300.1	818	intronless	5.82	29889.42	chr5A:419207299-419208117	CYC/TB1
TaTCP24-B	TraesCS5B01G205600.1	824	intronless	5.94	30149.55	chr5B:374190914-374192086	CYC/TB1
TaTCP24-D	TraesCS5D01G213400.1	827	intronless	5.94	30230.53	chr5D:322577220-322578047	CYC/TB1
TaPCF1-A	TraesCSU01G021700LC.1	863	intronless	6.45	13470.32	chrUn:12584619-12585482	PCF
TaPCF1-B	TraesCS2B01G062900.1	482	intronless	9.79	16586.69	chr2B:30145219-30145701	PCF
TaPCF1-D	TraesCS2D01G049400.1	266	intronless	4.56	9326.79	chr2D:18210500-18210766	PCF
TaPCF2-A	TraesCS4A01G355900.1	889	2	5.28	30860.2	chr4A:629486061-629487401	PCF
TaPCF2-B	TraesCS5B01G516200.1	1049	intronless	5.43	36115.83	chr5B:680580273-680583936	PCF
TaPCF2-D	TraesCS5D01G516300.1	1034	1	5.54	35912.84	chr5D:539014262-539015611	PCF
TaPCF3-A	TraesCS4A01G222700.1	1157	intronless	6.98	40576.52	chr4A:530121678-530123374	PCF
TaPCF3-B	TraesCS4B01G093700.1	1145	intronless	6.81	40160.3	chr4B:96270932-96272626	PCF
TaPCF3-D	TraesCS4D01G090500.1	1157	intronless	6.98	40601.52	chr4D:65576568-65578302	PCF
TaTCP6-A	TraesCS3A01G436100.1	1025	1	6.27	35905.49	chr3A:679694378-679696479	PCF
TaTCP6-B	TraesCS3B01G470700.1	1049	2	6.35	36864.9	chr3B:718849770-718851635	PCF
TaTCP6-D	TraesCS3D01G428700.1	1010	2	6.23	35492.25	chr3D:543008179-543010159	PCF
TaTCP7-A	TraesCS6A01G233800.1	611	intronless	9.51	21371.71	chr6A:441486019-441486880	PCF
TaTCP7-B	TraesCS6B01G262600.1	598	intronless	9.51	20937.53	chr6B:473055822-473056442	PCF
TaTCP7-D	TraesCS6D01G216100.1	620	1	9.51	21813.98	chr6D:306129040-306130108	PCF
TaTCP9-A	TraesCS6A01G306500.1	1214	intronless	9.17	40755.88	chr6A:540144123-540146045	PCF
TaTCP9-B	TraesCS6B01G334900.1	1208	intronless	9.19	40408.73	chr6B:589127612-589129644	PCF
TaTCP9-D	TraesCS6D01G285600.1	1205	intronless	9.19	40460.8	chr6D:393934793-393936464	PCF
TaTCP17-A	TraesCS2A01G376000.1	617	intronless	9.61	21774.02	chr2A:618266685-618267625	PCF
TaTCP17-B	TraesCS2B01G392900.1	755	intronless	9.95	26968.79	chr2B:557226759-557227975	PCF
TaTCP17-D	TraesCS2D01G372200.1	617	intronless	9.61	21744.01	chr2D:476344966-476346144	PCF
TaTCP19-A	TraesCS7A01G244300LC.1	680	intronless	9.81	23772.15	chr7A:140536199-140537919	PCF
TaTCP19-B	TraesCS7B01G091900.1	1118	1	9.59	38470.11	chr7B:105322710-105325639	PCF
TaTCP19-D	TraesCS7D01G187600.1	1202	intronless	8.98	41495.36	chr7D:140232150-140234052	PCF
TaTCP25-A	TraesCS5A01G309300.1	1016	Intronless	5.06	34474.11	chr5A:521743536-521744868	PCF
TaTCP25-B	TraesCS5B01G310000.1	1019	intronless	4.95	34568.02	chr5B:491993823-491995127	PCF
TaTCP25-D	TraesCS5D01G316300.1	1004	intronless	5.03	33989.76	chr5D:410240816-410242165	PCF
TaTCP28-A	TraesCS5A01G121800.1	1259	intronless	5.53	43434.48	chr5A:264421099-264422849	PCF
TaTCP28-B	TraesCS5B01G117100.1	1253	intronless	5.66	43278.42	chr5B:206101060-206102857	PCF
TaTCP28-D	TraesCS5D01G199700LC.1	1655	intronless	6.47	57569.72	chr5D:205995024-205997202	PCF
TaTCP29-A	TraesCS6A01G412800.1	596	1	6.34	19973.82	chr6A:614381484-614383233	PCF
TaTCP29-B	TraesCS6B01G462100.1	596	intronless	6.58	20230.98	chr6B:715448937-715450031	PCF
TaTCP29-D	TraesCS6D01G397000.1	593	1	6.54	20051.93	chr6D:469205255-469207520	PCF

Table 2

Ka/Ks analysis and estimated divergence time for TaTB1 gene pairs in wheat and its donor species.

Paralogous pairs	Origin	Ka	Ks	Ka/Ks	Divergence time (mya)
TaTB1-1-A-TaTB1-2-A	AABBDD	0.422331	1.76758	0.238932	0.718717
TaTB1-1-B-TaTB1-3-B	AABBDD	0.39819	1.42517	0.279399	0.61836
TaTB1-1-D-TaTB1-4-D	AABBDD	0.391273	1.89824	0.206124	0.710809
TtTB1-1-A-TtTB1-2-A	AABB	0.389092	1.96989	0.19752	0.71316
AetTB1-1-AetTB1-2	DD	0.366917	2.7072	0.135534	0.839375
TuTB1-1-TuTB1-2	AA	0.366908	1.54729	0.23713	0.615564

Table 3

Enriched cis-acting elements in TaTCP gene promoter regions.

Motif altmetric ID	Consensus	Motif name	p-value	Adjusted p-value
Plant development
CCGTCC	CCGTCC	CCGTCC-box	4.12E-10	1.52E-07
GATAATGATG	GATAATGATG	as-2-box	1.14E-04	8.61E-03
CATGCATG	CATGCATG	RY-element	1.12E-06	3.34E-04
Hormone response
AAATGGAGA	AAATGGAGA	P-box	9.38E-04	0.204∗
AACGAC	AACGAC	TGA-element	8.02E-04	0.274∗
Biotic and abiotic response
AAACAGA	AAACAGA	ARE	2.80E-06	5.24E-04
AAAAAATTTC	AAAAAATTTC	HSE	9.13E-08	1.57E-05

Table 4

Protein-protein interaction of selected wheat TCP genes.

AD	TaTCP18-D	TaTCP21-D	TaPCF3-A
BD
TaTCP18-D	✓	✓	✓
TaTCP21-D	✓	✓	✓
TaPCF3-A	X	X	X