Background: Flowering time has an important effect on regional adaptation and yields for crops. The tyrosine kinase-like (TKL) gene family is widely existed and participates in many biological processes in plants. Furthermore, only few TKLs have been characterized functions in controlling flowering time in wheat. Results: Here, we report that TaCTR1, a tyrosine kinase-like (TKL) gene, regulates flowering time in wheat. Based on identification and evolutionary analysis of TKL_CTR1-DRK-2 subfamily in 15 plants, we proposed an evolutionary model for TaCTR1, suggesting that occurrence of some exon fusion events during evolution. The overexpression of TaCTR1 caused early flowering time in transgenic lines. Transcriptomics analysis enabled identification of mass differential expression genes including plant hormone (ET, ABA, IAA, BR) signaling, flavonoid biosynthesis, phenolamides and antioxidant, and flowering-related genes in TaCTR1 overexpression transgenic lines compared with WT plants. qRT–PCR results showed that the expression levels of ethylene (ET) signal-related genes (ETR, EIN, ERF) and flowering-related genes (FT, PPD1, CO, PRR, PHY) were altered in TaCTR1-overexpressing wheat compared with WT plants. Metabonomics analysis showed that flavonoid contents were altered. Conclusions: Thus, the results show that TaCTR1 plays a positive role in controlling flowering time by activating various signaling pathways and regulating flowering-related genes, and will provide new insights on the mechanisms of wheat flowering regulation.
Keywords: Serine/threonine-protein kinases; TaCTR1; Flowering; Transcriptomics analysis
Peisen Su and Jingyu Li contributed equally to this work.
Vegetative and reproductive growth are the two major phases of the plant life cycle. The floral transition is one of the major transitions occurring between these two phases, and it determines plant environmental adaptation, survival, and grain yield. Plants have evolved complex regulatory networks in response to environmental signals, including vernalization, temperature, photoperiod, the phytohormone gibberellin (GA), and age, to tightly control flowering [[
Eukaryotic protein kinases belong to a large gene superfamily. They contain a protein kinase catalytic domain that phosphorylates the serine, threonine, or tyrosine residues of target proteins [[
TKLs belong to the serine/threonine kinase family in plants and are widespread among plants [[
Wheat is an important cereal crop worldwide. Flowering time affects regional wheat adaptation and yields. Here, we identified the function of one wheat TKL_CTR1-DRK-2 gene, TaCTR1, in controlling flowering time. Analysis of gene identification and evolution indicated that two exon fusion events occurred during the evolution of the TKL_CTR1-DRK-2 subfamilies I-IV. The overexpression of TaCTR1 in wheat triggered early flowering through the activation of the ET signalling pathway and flowering-related genes. Our findings provide new insights into the regulatory roles of serine/threonine-protein kinases in determining wheat heading date.
The TKL_CTR1-DRK-2 gene family belongs to the PK (protein kinase) gene superfamily studied in our previous article [[
Graph: Fig. 1 Classification and phylogenetic tree of the TKL_CTR1-DRK-2 in 15 plants
In order to obtain further insights into TKL_CTR1-DRK-2 evolution, we diagrammed the exon-intron structures within the kinase domain in the 15 investigated plants (Fig. S2). We found that TKL_CTR1-DRK-2 subfamilies I–IV contained a conserved exon-intron structure within the exon phases in the kinase domain (Fig. S3). For instance, the subfamily I contained conserved exon-intron structures within the "1002-0002-02" exon phases in the kinase domain. Subfamily II contained "1002-0000-2", and subfamilies III and IV contained "1002-0002-2". Interestingly, P. patens contained only one member (Pp3c12_3550V3.2) of TKL_CTR1-DRK-2, and its exon phases were "1002-0002-02," hinting that subfamily I was the ancient subfamily. Indeed, Pp3c12_3550V3.2 was located in the root of subfamilies I and II clades among the bayesian and ML phylogenetic trees (Fig. S1).
We also studied the evolution of the TKL_CTR1-DRK-2 genes in Triticum species (Fig. S4). First, there was a "1-2-3-fold" trend in the distributions of members of the TKL_CTR1-DRK-2 subfamilies. In Zea mays and Brachypodium distachyon, the distributions of TKL_CTR1-DRK-2 subfamily members were the same (11 in total: 3 each in I, II, and III and two in IV). This distribution of subfamily III and IV members was similar in Aegilops tauschii (D genome, 12 in total: 5 in I, 2 in II, 3 in III and 2 in IV). Compared with those in Ae. tauschii, the number of members of subfamilies I-IV were twice as high in Triticum dicoccoides (AB subgenome, 25 in total: 11 in I, 4 in II, 6 in III and 4 in IV) and Triticum turgidum (AB subgenome, 26 in total: 12 in I, 4 in II, 6 in III and 4 in IV). Similarly, compared with those in Ae. tauschii, the number of members of subfamilies I-IV were three times as high in Triticum aestivum (ABD subgenome, 40 in total: 18 in I, 6 in II, 9 in III, 6 in IV, with 1 member excluded). The four phylogenetic trees also supported that the TKL_CTR1-DRK-2 genes of Triticum species underwent a genome triploidization event during their evolution (Fig. 1 and Fig. S1). For example, a subclade of subfamily I contained a "1-2-3-fold" trend with a bootstrap value of 100 in Ae. tauschii (D genome, Ata AET6Gv20792500.7), Triticum urartu (A genome, Tur TRIUR3 17,777-P1), T. dicoccoides (AB subgenome, Tdi TRIDC6BG056450.9 and Tdi TRIDC6AG048310.2), T. turgidum (AB subgenome, Ttu TRITD6Bv1G188650.3 and Ttu TRITD6Av1G195350.3), T. spelta (ABD subgenome, Tsp TraesTSP6D01G320200.1, Tsp TraesTSP6B01G378500.1 and Tsp TraesTSP6A01G304500.1) and T. aestivum (ABD subgenome, Tae TraesCS6D02G301700.1, Tae TraesCS6B02G352700.1 and Tae TraesCS6A02G321900.1).
How did these four types of exon phase arrangements of TKL_CTR1-DRK-2 evolve to form the subfamilies I–IV? Based on the analysis above, we inferred an exon-intron evolutionary model with two exon fusion events of TKL_CTR1-DRK-2 subfamilies I–IV (Fig. 2). This model suggests that two rounds of exon fusion events occurred during the emergence of Pteridophyta and Angiosperms (Fig. 2A). The first exon fusion event—two exons fused into one exon—emerged in subfamilies I and II during the Pteridophyta emergence. The exon-intron structure of P. patens TKL_CTR1-DRK-2 subfamily I (Pp3c12_3550V3.2) is "1002-0002-02", while the exon-intron structure of the subfamily II of TKL_CTR1-DRK-2 genes from S. moellendorffii to Angiosperms evolved into "1002-0000-2". Two exons behind within the exon phase "20" in the subfamily I fused into one exon within the exon phase "0" and formed the subfamily II of TKL_CTR1-DRK-2. The evidence of alignment within the breakpoints of exon fusion events between some subfamilies I and II genes supported our proposed evolutionary model (Fig. 2B). Similarly, the second exon fusion event—two exons fused into one exon—emerged in the subfamilies I, III, and IV during the Angiosperm emergence, or even earlier. The last two exons of the subfamily I within the exon phase "02" fused into one exon within the exon phase "2" and formed the ancestor of the TKL_CTR1-DRK-2 subfamilies III and IV. The exon-intron structure "1002-0002-2" of subfamilies III and IV remained in Angiosperms. We also aligned nucleotide sequences within the breakpoints of exon fusion events (Fig. S5).
Graph: Fig. 2 Exon − intron evolutionary model of TKL_CTR1-DRK-2 subfamilies I-IV. (A) Two exon fusion events of subfamily II and III-IV TKL_CTR1-DRK-2 in evolution. (B) Alightment of exon fusion blocks with amino acid sequences. Filled boxes: red represents the PK_Tyr_Ser-Thr domain; black boxes: untranslated regions (UTRs); white boxes: other exon regions; lines: introns; numbers 0, 1, and 2: exon phases. The lengths of the boxes and lines are scaled based on the lengths of the genes. The long introns are shorted by "//". Exon fusion events are circled by blue boxes. Exon phases are circled by red boxes and arrows
We diagrammed the domain structures of the 15 investigated plants (Fig. S6). The evolution of domains in representative species is summarized in Fig. 3. C. reinhardtii TKL_CTR1-DRK-2 (PNW85595) contains two domains (EDR1 and PK_Tyr_Ser-Thr), and P. patens TKL_CTR1-DRK-2 (Pp3c12_3550V3.2) contains three domains (PAS, EDR1, and PK_Tyr_Ser-Thr). Subfamilies I and II had already appeared in S. moellendorffii, and they contain different domain compositions (I: PAS and PK_Tyr_Ser-Thr; II: EDR1 and PK_Tyr_Ser-Thr). Later, the PAS domain of one subfamily I member became the PAS_9 domain in the ancestor of eudicots and monocots. Surprisingly, the subfamily I member within the PAS_9 domain was lost in the eudicots V. vinifera and A. thaliana but was retained in the primitive eudicot (A) trichopoda and some monocots, such as (B) distachyon and T. aestivum. Subfamilies III and IV emerged in eudicots and monocots, and their domain compositions were the same as those of subfamily II (EDR1 and PK_Tyr_Ser-Thr).
Graph: Fig. 3 Domain evolutionary model of TKL_CTR1-DRK-2 subfamilies I-IV. Filled boxes: red represents PK_Tyr_Ser-Thr domain; green represents PAS domain; purple represents PAS_9 domain; orange represents EDR1 domain. PAS and PAS_9 domains were circled by blue box in subfamily I. The lengths of the boxes and lines are scaled based on the length of proteins (X axis)
Gene families are generated from gene duplication events on genome-wide or local chromosomal scales. Type I functional divergence indicates that the amino acid pattern is highly conserved in one duplicate cluster but highly variable in the other cluster [[
We mapped the locations of 40 TKL_CTR1-DRK-2 genes on 21 T. aestivum chromosomes (Fig. S9). About half of T. aestivum chromosomes contained only one TKL_CTR1-DRK-2. However, the other T. aestivum chromosomes, including 3 A (2 TKL_CTR1-DRK-2 genes), 3D (
Graph: Fig. 4 Synteny analysis of TKL_CTR1-DRK-2 genes. This graph displays syntenic maps of T. aestivum associated with two Graminaceae (B. distachyon and O. sativa). Red curves represent syntenic gene pairs between the TKL_CTR1-DRK-2 genes, and grey curves represent other genes. (A) Synteny of TKL_CTR1-DRK-2 genes in T. aestivum. (B) Synteny of TKL_CTR1-DRK-2 genes and other genes in T. aestivum. (C) Synteny of TKL_CTR1-DRK-2 genes in T. aestivum and B. distachyon. (D) Synteny of TKL_CTR1-DRK-2 genes and other genes in T. aestivum and B. distachyon. (E) Synteny of TKL_CTR1-DRK-2 genes in T. aestivum and O. sativa. (F) Synteny of TKL_CTR1-DRK-2 genes and other genes in T. aestivum and O. sativa
The allohexaploid bread wheat (T. aestivum) genome contains three subgenomes (A, B, and D), which are the products of polyploidization. The A and B genomes diverged 6.5 Mya, and they generated the D genome 1–2 million years later [[
We assessed the expression patterns of wheat TKL genes from four public RNA-seq datasets from the NCBI SRA during development (Fig. S11). Based on the results of FastQC analysis, we performed a quality control assessment and excluded one RUN file from the following analysis (Table S8). (
Some TaCTR information, including protein sequence, coding sequence (CDS), chromosome information, GFF information, exon − intron and kinase domain diagrams, and BLAST results, are shown in Fig. S14. Analysis of public RNA-Seq data indicated that the expression of TaCTR1 was upregulated at the heading stage, followed by a "maintained high expression level" (Fig. S11). To further elucidate the biological function of TaCTR1, we generated TaCTR1-overexpressing transgenic lines in the common wheat "JWI" background and successfully genotyped three independent homozygous overexpression transgenic lines. The transcript levels of TaCTR1 in TaCTR1-overexpressing transgenic lines were analysed by qRT‒PCR. The results indicated that the expression levels of TaCTR1 were greater in the transgenic lines than in the WT plants (Fig. 5A). Subsequently, the transgenic and WT plants were planted in a greenhouse and subjected to a detailed phenotypic analysis. The TaCTR1 transgenic lines flowered seven days earlier than did the WT plants (Fig. 5B and C). In summary, the overexpression of TaCTR1 promoted flowering time in wheat.
Graph: Fig. 5 TaCTR1 overexpression causes early flowering in wheat. (A) The expression levels of TaCTR1 in transgenic lines and WT plants. (B) Comparison of flowering phenotypes among TaCTR1 overexpression transgenic and wild-type (WT) plants grown in the greenhouse. (C) Flowering time of TaCTR1 overexpression transgenic and WT plants. Data are means ± SD of 15 plants. (*P < 0.05; **P < 0.01; t-test)
In order to better understand the regulatory mechanisms by which TaCTR1 affects wheat flowering, we performed a transcriptome analysis comparing TaCTR1 overexpression transgenic lines and JWI wild-type plants. Our transcriptomic data indicated that 13,312 DEGs (8610 upregulated and 4702 downregulated) were identified between transgenic lines and JWI plants (Fig. S15A, Table S9). GO classification analyses revealed that molecular functions (lipase activity, glucosidase activity, chlorophyll binding, and microtubule binding), biological processes (cell tip growth, developmental cell growth, and post-embryonic plant organ morphogenesis), and cell components (photosystem II, secretory vesicle, and photosystem I) were highly enriched terms (Fig. S15B, Table S10). KEGG enrichment analyses indicated that DEGs were mainly enriched in the biosynthesis of secondary metabolites, plant hormone signal transduction, MAPK signaling pathways, carbon metabolism, and starch and sucrose metabolism (Fig. S15C, Table S10).
To further study how TaCTR1 affects plant flowering, we generated a network to illustrate the relationships between flavonoid biosynthesis and TaCTR1 overexpression. We found that the transcript levels of 77 genes in the flavonoid biosynthesis pathway were altered, among which naringenin 3-dioxygenase (F3H) was upregulated compared to the control, while phenylalanine ammonia-lyase (PAL) and flavonoid 3'-monooxygenase (CYP75B1) were downregulated (Fig. 6, Table S11). In addition, the expression profiles of flavanone 4-reductase (DFR), flavonoid 3',5'-hydroxylase (CYP75A), 4-coumarate–CoA ligase (4CL), anthocyanidin reductase (ANR), and flavonol-3-O-L-rhamnoside-7-O-glucosyltransferase (UGT73C6) were combined; detailed information is displayed in Table S11.
Graph: Fig. 6 TaCTR1 overexpression alters flavonoid biosynthesis pathway. The network of flavonoid biosynthesis pathway in TaCTR1 overexpression transgenic lines, compared with WT plants. PAL, phenylalanine ammonia-lyase; 4CL, 4-coumarate–CoA ligase; CHS, chalcone synthase; DFR, flavanone 4-reductase; ANR, anthocyanidin reductase; FLS, flavonol synthase; UGT73C6, flavonol-3-O-L-rhamnoside − 7-O-glucosyltransferase; F3H, naringenin 3-dioxygenase; CYP75A, flavonoid 3',5'-hydroxylase; CYP75B1, flavonoid 3'-monooxygenase, FG2, flavonol-3-O-glucoside L-rhamnosyltransferase. The red and green indicate high and low expression level, respectively
To further explore which type of flavonoid metabolites were altered due to TaCTR1 overexpression, we measured the flavonoid content in TaCTR1-overexpressing lines and WT plants. The results showed that 8 flavonoid-related metabolites, including swertisin, apigenin-7-O-(2"-glucosyl)arabinoside, and apigenin-7-O-rutinoside-4'-O-rhamnoside, were upregulated. Eighty-nine flavonoid-related metabolites, including kaempferol-3,7-O-diglucoside, apigenin-7,4'-dimethyl ether, 6-hydroxy-2'-methoxyflavone, and abrusin, decreased in abundance (Table S12). TaCTR1 overexpression promoted early flowering in the transgenic lines by stimulating flavonoid biosynthesis.
Our transcriptome data also indicated that phytohormone pathways were activated in TaCTR1 overexpression transgenic lines. Therefore, we generated a network to explore the relationships between phytohormone pathways and TaCTR1 overexpression. A large number of genes were enriched in plant hormone pathways: 32 ET-related genes were altered, among which adenosylhomocysteinase (AHCY) and DNA (cytosine-5)-methyltransferase 1 (DCM) were upregulated compared to the control, while ethylene-responsive transcription factors (ERF) were downregulated (Fig. 7, Table S11). We also performed qPCR to analyze the expression levels of known ET-related genes in transgenic lines and WT plants. The results indicated that the transcript levels of several ET-related genes, such as ethylene receptor (ETR), ethylene insensitive (EIN), and ethylene-responsive transcription factor (ERF), were significantly increased in TaCTR1 overexpression lines compared with WT plants (Fig. 8).
Graph: Fig. 7 TaCTR1 overexpression alters ethylene-related genes expression. The network of ET biosynthesis and signalling, phenolamides and antioxidant pathways in TaCTR1 overexpression transgenic lines, compared with WT plants. BHMT, homocysteine S-methyltransferase; DCM, DNA (cytosine-5)-methyltransferase 1; AHCY, adenosylhomocysteinase; ACO, aminocyclopropanecarboxylate oxidase; MKK, mitogen-activated protein kinase kinase; MPK, mitogen-activated protein kinase; EIN2, ethylene-insensitive protein 2; GPX, glutathione peroxidase; GST, glutathione S-transferase; APX, ascorbate peroxidase; RRM1, ribonucleoside-diphosphate reductase subunit M1; ODC1, ornithine decarboxylase; SMOX, spermine oxidase; ERF, ethylene-responsive transcription factor. The red and green indicate high and low expression level, respectively
Graph: Fig. 8 TaCTR1 overexpression alters the relative expression of ET signal-related genes. (A) The expression profile of ETR. (B) The expression profile of EIN. (C) The expression profile of ERF. All experiments included three replicates. The wheat gene 18SrRNA was used as an endogenous control. The gene expression profile were calculated using the 2–ΔΔCT method
In addition, 64 ABA-related genes were found, among which abscisic-aldehyde oxidase (AAO3) and protein phosphatase 2 C (PP2C) were upregulated compared to the control, while 15-cis-phytoene synthase (ctrB), prolycopene isomerase (crtH), zeaxanthin epoxidase (ABA1), and 9-cis-epoxycarotenoid dioxygenase (NCED) were downregulated (Fig. S16, Table S11). 131 auxin-related genes—including anthranilate synthase (TRP3), L-tryptophan-pyruvate aminotransferase (TAA1), indole-3-pyruvate monooxygenase (YUCCA), aldehyde dehydrogenase (ALDH), and auxin-responsive GH3 gene family (GH3)—33 SA-related genes (NPR1, TGA, PR-1), 51 JA-related genes (hydroperoxide dehydratase (AOS), lipoxygenase (LOX), COI1, and MYC2), and 110 BR-related genes (steroid 5-alpha-reductase (DET2), brassinosteroid 6-oxygenase (CYP85A1), BR-signaling kinase (BSK), protein brassinosteroid insensitive 1 (BRI1)) were altered upon TaCTR1 overexpression. The detailed information is displayed in Fig. 9, Table S11.
Graph: Fig. 9 TaCTR1 overexpression alters auxin signaling pathways. The network of auxin signaling pathways in TaCTR1 overexpression transgenic lines, compared with WT plants. TPR3, anthranilate synthase; DDC, aromatic-L-amino-acid; TAA1, L-tryptophan—pyruvate aminotransferase;ALDH, aldehyde dehydrogenase (NAD+); YUCCA, indole-3-pyruvate monooxygenase; AMI, amidase; AUX/IAA, auxin-responsive protein. The red and green indicate high and low expression level, respectively
In plants, S-adenosyl-L-methionine (SAM) is a precursor of ET and polyamine biosynthesis [[
Flowering time is regulated by flowering-related genes, such as CO, FT, and PPD1 [[
Graph: Fig. 10 TaCTR1 overexpression alters circadian rhythm pathway. (A) The network of circadian rhythm pathway in TaCTR1 overexpression transgenic lines, compared with WT plants. The blue indicated significantly differences in gene expression levels, while red (up-regulated) and green (down-regulated) represented specific expression between TaCTR1 overexpression transgenic lines and WT plants. PIF3, phytochrome-interacting factor 3; PRR, pseudo-response regulator; LHY, MYB-related transcription factor; CHE, protein CCA1 HIKING EXPEDITION; TOC1, pseudo-response regulator 1; COP1, E3 ubiquitin-protein ligase RFWD2; CK2α, casein kinase II subunit alpha; CHS, chalcone synthase; FT, FLOWERING LOCUS T. (B) The heatmap of DEGs in circadian rhythm pathway. The red and green indicate high and low expression level, respectively
Graph: Fig. 11 TaCTR1 overexpression alters expression of flowering-related genes. (A) The expression profile of FT. (B) The expression profile of PPD1. (C-D) The expression profile of CO. (E-F) The expression profile of PRRs. (G-K) The expression profile of PHY genes. The wheat gene 18SrRNA was used as an endogenous control. The gene expression profile were calculated using the 2–ΔΔCT method. All experiments included three replicates
The CTR1-DRK-2 gene family belongs to the PK gene superfamily, which has expanded in the flowering plant lineage [[
First, CTR1-DRK-2 subfamily I might have already appeared in a green algae-resembling ancestor. We identified one CTR1-DRK-2 gene (PNW85595) in green algae (C. reinhardtii), suggesting that the ancestral sequence of the CTR1-DRK-2 gene family might have arisen during the Phycophyta emergence. This sequence contained only two domains (EDR1 and PK_Tyr_Ser-Thr) (Fig. 3), and its exon‒intron structure was algae specific. As displayed in Fig. 3, P. patens had only one CTR1-DRK-2 gene (Ppa_Pp3c12_3550V3.2), but its domain arrangement had turned into three domains (PAS, EDR1, and PK_Tyr_Ser-Thr), suggesting that Ppa_Pp3c12_3550V3.2 obtained a new PAS domain during evolution. Moreover, the exon phases from the PK_Tyr_Ser-Thr domain changed to a "1002-0002-02" pattern, which was mostly retained during the evolution of eudicots and monocots. Based on the analysis of phylogenetic trees, the Ppa_Pp3c12_3550V3.2 sequence belongs to subfamily I, suggesting that subfamily I was the ancestral subfamily. Second, subfamilies I and II might have appeared earlier in a fern-resembling ancestor. Subfamily II may be the second subfamily because of its sequences (Smo_EFJ10671 and Smo_EFJ07496) in S. moellendorffii. These genes all contained two domains (EDR1 and PK_Tyr_Ser-Thr) that were retained in subfamilies II–IV during evolution. Moreover, the subfamily I sequence (Smo_EFJ04327) in S. moellendorffii contained two domains (PAS and PK_Tyr_Ser-Thr). Considering the three domains (PAS, EDR1, and PK_Tyr_Ser-Thr) in the P. patens sequence (Ppa_Pp3c12_3550V3.2), the subfamily I sequence (Smo_EFJ04327) in S. moellendorffii lost its EDR1 domain during evolution, and subfamily II (Smo_EFJ10671 and Smo_EFJ07496) lost its PAS domain. One exon fusion event occurred in S. moellendorffii subfamily II (Smo_EFJ10671 and Smo_EFJ07496) and changed to a "1002-000-02" structure. These structures were conserved from S. moellendorffii to T. aestivum. Third, the common ancestral sequence of subfamilies III and IV might have appeared during the emergence of angiosperms. Their two domain combinations (EDR1 and PK_Tyr_Ser-Thr) and "1002-0002-2" exon‒intron structure are conserved in both eudicots and monocots.
The PKs are known to function in plant development and growth. For example, the knockdown of CPK32 has been shown to result in late flowering by altering FLOWERING CONTROL LOCUS A (FCA) alternative polyadenylation and FLOWERING LOCUS C (FLC) transcription [[
Many studies have demonstrated that plant flowering time is regulated by various pathways, such as circadian rhythm and flavonoid biosynthesis [[
Here, we identified and characterized the role of the tyrosine kinase-like (TKL) gene TaCTR1 in flowering time. Phylogenetic analysis revealed the occurrence of some exon fusion events during evolution. The overexpression of TaCTR1 caused early flowering in transgenic wheat through activation of the plant hormone signaling, flavonoid biosynthesis, phenolamides and antioxidant, and flowering-related genes. This work and its findings provide insights on wheat controlling flowering time and the regulatory mechanism while illustrating the potential of TaCTR1 overexpression in altering flowering time during crop improvement programs.
The wheat (Triticum aestivum L.) cultivar "Sumai3" which was sourced from the College of Agronomy (Shandong Agricultural University), was used in this study. The cultivar "JWI" was used for transformation to construct overexpression lines. The transgenic wheat lines and "JWI" were grown at 20–25 °C with a photoperiod of 16 h light/8 h dark in a greenhouse.
The genome and protein sequences of T. aestivum, T. spelta, T. turgidum, T. dicoccoides, T. urartu, Ae. tauschii, B. distachyon, Z. mays, Oryza sativa, Arabidopsis thaliana, Vitis vinifera, Amborella trichopoda, Selaginella moellendorffii, Physcomitrella patens, and Chlamydomonas reinhardtii were downloaded from Ensembl Plant release-51 (
The alignment of truncated TKL_CTR1-DRK-2 sequences with the PK_Tyr_Ser-Thr domain was performed by ClustalX v2.0 [[
The domain and intron–exon structures of TKL_CTR1-DRK-2 sequences in these fifteen plants were generated by our Perl and R scripts based on the corresponding GFF file information from Ensembl Plant release-51 (
Type I functional divergence analysis was performed by DIVERGE (version 2.0) software [[
Type I functional divergence which resulted in altered functional constraints between duplicated genes, leads to one of the genes being conserved and the other gene varying highly. Coefficient of functional divergence (θ) is an indicator for the level of type I functional divergence among two homologous gene clusters. The sites (k) with contribution to the functional divergence were predicted according to their posterior probabilities (Qk). The sites with Qk > 0.67 were meaningful for the functional divergence. We chose 91 sequences (including Pkinase or PK_Tyr_Ser-Thr domain) of eight typical plants (including C. reinhardtii, P. patens, S. moellendorffii, A. trichopoda, (A) thaliana, (B) distachyon, Ae. tauschii and T. aestivum) to calculate the coefficients of functional divergence.
Based on the extracted information of GFF files from Ensembl Plants release-51 (
Public wheat (T. aestivum) expression datasets were retrieved from NCBI's Sequence Read Archive (SRA). The four categories identified from the RNA-seq data of wheat were grain development (Accession: PRJNA525250, ID: 525,250), early wheat spike development (Accession: PRJNA325489, ID: 325,489), four growth stages of wheat flag leaves (Accession: PRJNA656068, ID: 656,068), and wheat early and late heading date (Accession: PRJNA668815, ID: 668,815). Quality control assessment of raw data was performed using FastQC. High-quality RNA-seq reads were aligned to reference wheat (T. aestivum) genome of Ensembl Plants release-51 by software Hisat2. Counting of mapped reads was performed using Samtools and HTseq software. The reads per kilobase of exon model per million mapped reads (RPKM) algorithm was used to normalize the data. Heat maps of wheat TKL expression levels were generated using Mev4.9 [[
Total RNA from the wheat cultivar "Sumai 3" leaves was extracted with a TRIzol reagent (Transgen). The full-length cDNA sequence of TaCTR1 (TraesCS4D01G010200) was obtained from Ensembl plants (
The transcriptome sequencing was conducted on the leaves and stems of JWI and TaCTR1 overexpression transgenic lines collected at elongation stages. The total RNA of samples was extracted using the RNAprep Pure Plant Kit (Transgen) and was sent to Wuhan Metware Biotechnology Co., Ltd. (Wuhan, China) for RNA sequencing. The transcriptome analysis was performed using an Illumina HiSeq™2000. The clean reads were mapped to the wheat reference genome (
Total RNA of all samples was extracted using RNAprep Pure Plant Kit (Transgen), and then reverse-transcribed into cDNA. The cDNA was used as template for expression analysis. qRT–PCR were performed using the Roche LightCycler ®480 system (Roche, Germany). The wheat gene 18SrRNA was used as an endogenous control. The gene expression profile were calculated using the 2
The bayesian phylogenetic tree was constructed by the amino acid sequences of the PK_Tyr_Ser-Thr domain using MrBayes v3.2.7. (A) Subfamilies I-IV. Detailed information is provided in Supplementary Fig. S1a. (B) Expanded subfamily II and other compressed subfamilies. Our two studied TKL_CTR1-DRK-2 sequences (TraesCS4D02G010200.1 and KQJ91573) in T. aestivum and B. distachyon were circled by red boxes.
Not applicable.
PSS, XWH, JY and WYG conceived and designed the experiments. PSS performed most experiments. JY and WYG performed the bioinformatics analysis of TaCRT1, and revised the manuscript. YJL, YFN, FTS, TYS, BJH and ZJL performed the screening of homozygous transgenic lines, RNA extraciton and expression analysis. PSS and XWH wrote and revised the manuscript. All authors have read and approved the final manuscript.
This work was funded by Natural Science Foundation of Shandong Province (ZR2022QC129); Doctoral research start-up funds, Liaocheng University (318052018); State Key Laboratory of Crop Biology, Shandong Agricultural University (2021KF03); Undergraduate innovation and entrepreneurship training program, Liaocheng University (CXCY2023260, CXCY2023287) and Shandong Province Higher Educational Program for Introduction and Cultivation of Young Innovative Talents in 2021; the Construction fund of Liaocheng University for the 14th Five-Year Top Characteristic Disciplines Landscape Architecture (31946221201). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Data availabilityThe genomes, proteomes and GFF files of the investigated plants are available in Ensembl Plants 51-release (
The plant materials used in this study were grown in the greenhouse of Liaocheng university for research use only. The seeds of plant materials was obtained from the College of Agronomy, Shandong Agricultural University, Tai'an, Shandong, PR China and are not listed in the Convention on International Trade in Endangered Species of Wild Fauna and Flora Appendices I, II and III (valid from 4 April 2017, https://cites.org/eng/app/appendices.php). Collection of plant materials complied with the institutional, national and international guidelines.
Not applicable.
The authors declare no competing interests.
Below is the link to the electronic supplementary material.
Graph: Supplementary Material 1
Graph: Supplementary Material 2
Graph: Supplementary Material 3
Graph: Supplementary Material 4
Graph: Supplementary Material 5
Graph: Supplementary Material 6
Graph: Supplementary Material 7
Graph: Supplementary Material 8
Graph: Supplementary Material 9
Graph: Supplementary Material 10
Graph: Supplementary Material 11
Graph: Supplementary Material 12
Graph: Supplementary Material 13
Graph: Supplementary Material 14
Graph: Supplementary Material 15
Graph: Supplementary Material 16
Graph: Supplementary Material 17
Graph: Supplementary Material 18
Graph: Supplementary Material 19
Graph: Supplementary Material 20
Graph: Supplementary Material 21
Graph: Supplementary Material 22
Graph: Supplementary Material 23
Graph: Supplementary Material 24
Graph: Supplementary Material 25
Graph: Supplementary Material 26
Graph: Supplementary Material 27
Graph: Supplementary Material 28
Graph: Supplementary Material 29
Graph: Supplementary Material 30
Graph: Supplementary Material 31
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
By Peisen S. Su; Jingyu Li; Dongtian Zang; Zhiyu Wang; Yangyang Wu; Shatong Chi; Fanting Sun; Yufei Niu; Xuewen Hua; Jun Yan and Wenyang Ge
Reported by Author; Author; Author; Author; Author; Author; Author; Author; Author; Author; Author