Genome-wide characterization of aspartic protease (AP) gene family in Populus trichocarpa and identification of the potential PtAPs involved in wood formation

Comprehensive analysis of amino acid/auxin permease family genes reveal the positive role of GhAAAP128 in cotton tolerance to cold stress

Amino acid/auxin permeases (AAAPs) play crucial roles in plant development and response to environmental stimuli. They have been characterized at genome-wide levels in several plant species. However, little is known about the AAAP genes in Gossypium. Here, we identified 149, 141, 73, and 70 AAAPs from G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. All the AAAPs were categorized into eight subfamilies (Groups I-VIII). Moreover, we found that 182 and 179 AAAP paralogous gene pairs existed within G. hirsutum and G. barbadense genomes, respectively, and whole genome duplication (WGD) or segmental duplication contributed to the expansion of these AAAPs during evolution. Additionally, many cis-elements related to abiotic stress responses were detected in the promoter regions of GhAAAPs and GbAAAPs. Consistently, the expression of multiple AAAPs was significantly induced by NaCl, polyethylene glycol 6000, and especially cold stress. Among these GhAAAPs, GhAAAP128 had clearly positive roles in cotton and Arabidopsis seedling tolerance to cold stress. It may improve plant cold tolerance by up-regulating the expression of some anthocyanin synthesis genes rather than CBF (C-repeat binding factor) signaling genes. Our findings provide important information for further functional analysis of GhAAAPs in response to stressful cues, particularly cold stress in cotton plants.

Aspartic proteases gene family: Identification and expression profiles during stem vascular development in tobacco

Aspartic proteases (APs) constitute a large family in plants and are widely involved in diverse biological processes, like chloroplast metabolism, biotic and abiotic stress responses, and reproductive development. In this study, we focused on overall analysis of the APs genes in tobacco. Our analysis included the phylogeny and cis-elements in the cell wall-associated promoters of these genes. To characterize the expression patterns of APs genes in stem vascular development. The tissue expression analysis showed that NtAED3-like was preferentially expressed in the differentiating xylem and phloem cells of the vascular system. Based on histochemical staining analysis showed that the NtAED3-like gene was specifically expressed in stem vascular tissue, root vascular tissue, and petiole vascular tissue. The TdT-mediated dUTP nick-end labeling (TUNEL) assay illustrated a delayed progression of programmed cell death (PCD) within the xylem of the ko-ntaed3a-like mutant, relative to the wild type. The mutant ko-ntaed3a-like exhibited a phenotype of thinning stem circumference and changed in xylem structure and lignin content. In addition, the two-dimension heteronuclear single quantum coherent nuclear magnetic resonance (2D-HSQC) analysis of three milled wood lignins (MWLs) showed that the content of β-O-4 connection in ko-ntaed3a-like decreased slightly compared with wild type. In conclusion, this study provides our understanding of the regulation of vascular tissue development by the NtAED3-like gene in tobacco and provides a better basis for determining the molecular mechanism of the aspartic protease in secondary cell wall (SCW) development.

Integrated omics approach reveals the molecular pathways activated in tomato by Kocuria rhizophila, a soil plant growth-promoting bacterium

Plant microbial biostimulants application has become a promising and eco-friendly agricultural strategy to improve crop yields, reducing chemical inputs for more sustainable cropping systems. The soil dwelling bacterium Kocuria rhizophila was previously characterized as Plant Growth Promoting Bacteria (PGPB) for its multiple PGP traits, such as indole-3-acetic acid production, phosphate solubilization capability and salt and drought stress tolerance. Here, we evaluated by a multi-omics approach, the PGP activity of K. rhizophila on tomato, revealing the molecular pathways by which it promotes plant growth. Transcriptomic analysis showed several up-regulated genes mainly related to amino acid metabolism, cell wall organization, lipid and secondary metabolism, together with a modulation in the DNA methylation profile, after PGPB inoculation. In agreement, proteins involved in photosynthesis, cell division, and plant growth were highly accumulated by K. rhizophila. Furthermore, "amino acid and peptides", "monosaccharides", and "TCA" classes of metabolites resulted the most affected by PGPB treatment, as well as dopamine, a catecholamine neurotransmitter mediating plant growth through S-adenosylmethionine decarboxylase (SAMDC), a gene enhancing the vegetative growth, up-regulated in tomato by K. rhizophila treatment. Interestingly, eight gene modules well correlated with differentially accumulated proteins (DAPs) and metabolites (DAMs), among which two modules showed the highest correlation with nine proteins, including a nucleoside diphosphate kinase, and cytosolic ascorbate peroxidase, as well as with several amino acids and metabolites involved in TCA cycle. Overall, our findings highlighted that sugars and amino acids, energy regulators, involved in tomato plant growth, were strongly modulated by the K. rhizophila-plant interaction.

Genome and GWAS analysis identified genes significantly related to phenotypic state of Rhododendron bark

As an important horticultural plant, Rhododendron is often used in urban greening and landscape design. However, factors such as the high rate of genetic recombination, frequent outcrossing in the wild, weak linkage disequilibrium, and the susceptibility of gene expression to environmental factors limit further exploration of functional genes related to important horticultural traits, and make the breeding of new varieties require a longer time. Therefore, we choose bark as the target trait which is not easily affected by environmental factors, but also has ornamental properties. Genome-wide association study (GWAS) of Rhododendron delavayi (30 samples), R. irroratum (30 samples) and their F1 generation R. agastum (200 samples) was conducted on the roughness of bark phenotypes. Finally, we obtained 2416.31 Gbp of clean data and identified 5 328 800 high-quality SNPs. According to the P-value and the degree of linkage disequilibrium of SNPs, we further identified 4 out of 11 candidate genes that affect bark roughness. The results of gene differential expression analysis further indicated that the expression levels of Rhdel02G0243600 and Rhdel08G0220700 in different bark phenotypes were significantly different. Our study identified functional genes that influence important horticultural traits of Rhododendron, and illustrated the powerful utility and great potential of GWAS in understanding and exploiting wild germplasm genetic resources of Rhododendron.

Identifying genetic variation associated with environmental gradients and drought-tolerance phenotypes in ponderosa pine

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.

Multivariate analysis and genetic dissection of staygreen and stem reserve mobilisation under combined drought and heat stress in wheat (Triticum aestivum L.)

Introduction: Abiotic stresses significantly reduce crop yield by adversely affecting many physio-biochemical processes. Several physiological traits have been targeted and improved for yield enhancement in limiting environmental conditions. Amongst them, staygreen and stem reserve mobilisation are two important mutually exclusive traits contributing to grain filling under drought and heat stress in wheat. Henceforth, the present study was carried out to identify the QTLs governing these traits and to identify the superiors' lines through multi-trait genotype-ideotype distance index (MGIDI) Methods: A mapping population consisting of 166 recombinant inbred lines (RILs) developed from a cross between HD3086 and HI1500 was utilized in this study. The experiment was laid down in alpha lattice design in four environmental conditions viz. Control, drought, heat and combined stress (heat and drought). Genotyping of parents and RILs was carried out with 35 K Axiom® array (Wheat breeder array). Results and Discussion: Medium to high heritability with a moderate to high correlation between traits was observed. Principal component analysis (PCA) was performed to derive latent variables in the original set of traits and the relationship of these traits with latent variables.From this study, 14 QTLs were identified, out of which 11, 2, and 1 for soil plant analysis development (SPAD) value, leaf senescence rate (LSR), and stem reserve mobilisation efficiency (SRE) respectively. Quantitative trait loci (QTLs) for SPAD value harbored various genes like Dirigent protein 6-like, Protein FATTY ACID EXPORT 3, glucan synthase-3 and Ubiquitin carboxyl-terminal hydrolase, whereas QTLs for LSR were found to contain various genes like aspartyl protease family protein, potassium transporter, inositol-tetrakisphosphate 1-kinase, and DNA polymerase epsilon subunit D-like. Furthermore, the chromosomal region for SRE was found to be associated with serine-threonine protein kinase. Serine-threonine protein kinases are involved in many signaling networks such as ABA mediated ROS signaling and acclimation to environmental stimuli. After the validation of QTLs in multilocation trials, these QTLs can be used for marker-assisted selection (MAS) in breeding programs.

Genome-Wide Characterization of Calmodulin and Calmodulin-like Protein Gene Families in Paulownia fortunei and Identification of Their Potential Involvement in Paulownia Witches' Broom

As significant Ca2+ sensors, calmodulin (CaM) and calmodulin-like proteins (CML), have been associated with a variety of environmental conditions in plants. However, whether CaMs/CMLs are related to the stress of phytoplasma infection has not been reported in Paulownia fortunei. In the current study, 5 PfCaMs and 58 PfCMLs were detected through a genome-wide investigation. The number of EF-hand motifs in all PfCaMs/CMLs varied. Bioinformatics analyses, including protein characteristics, conserved domain, gene structure, cis-elements, evolutionary relationship, collinearity, chromosomal location, post-translation modification site, subcellular localization and expression pattern analyses, represented the conservation and divergence of PfCaMs/CMLs. Furthermore, some PfCaMs/CMLs might be involved in plants' reaction to phytoplasma infection and exogenous calcium therapy, indicating these genes may play a role in abiotic as well as biotic stress responses. In addition, subcellular localization analysis showed that PfCML10 was located in the cell membrane and nucleus. In summary, these findings establish a stronger platform for their subsequent functional investigation in trees and further characterize their roles in Paulownia witches' broom (PaWB) occurrence.

Phylogenetic and AlphaFold predicted structure analyses provide insights for A1 aspartic protease family classification in Arabidopsis

Aspartic proteases are widely distributed in animals, plants, fungi and other organisms. In land plants, A1 aspartic protease family members have been implicated to play important and varied roles in growth, development and defense. Thus a robust classification of this family is important for understanding their gene function and evolution. However, current A1 family members in Arabidopsis are less well classified and need to be re-evaluated. In this paper, 70 A1 aspartic proteases in Arabidopsis are divided into four groups (group I-IV) based on phylogenetic and gene structure analyses of 1200 A1 aspartic proteases which are obtained from 12 Embryophyta species. Group I-III members are further classified into 2, 4 and 7 subgroups based on the AlphaFold predicted structures. Furthermore, unique insights of A1 aspartic proteases have been unraveled by AlphaFold predicted structures. For example, subgroup II-C members have a unique II-C specific motif in the C-extend domain, and subgroup IV is a Spermatophyta conserved group without canonical DTGS/DSGT active sites. These results prove that AlphaFold combining phylogenetic analysis is a promising solution for complex gene family classification.

Genetic architecture of fresh-market tomato yield

BACKGROUND

The fresh-market tomato (Solanum lycopersicum) is bred for direct consumption and is selected for a high yield of large fruits. To understand the genetic variations (distinct types of DNA sequence polymorphism) that influence the yield, we collected the phenotypic variations in the yields of total fruit, extra-large-sized fruit, small-sized fruit, or red-colored fruit from 68 core inbred contemporary U.S. fresh-market tomatoes for three consecutive years and the genomic information in 8,289,741 single nucleotide polymorphism (SNP) positions from the whole-genome resequencing of these tomatoes.

RESULTS

Genome-wide association (GWA) mapping using the SNP data with or without SNP filtering steps using the regularization methods, validated with quantitative trait loci (QTL) linkage mapping, identified 18 significant association signals for traits evaluated. Among them, 10 of which were not located within genomic regions previously identified as being associated with fruit size/shape. When mapping-driven association signals [558 SNPs associated with 28 yield (component) traits] were used to calculate genomic estimated breeding values (GEBVs) of evaluated traits, the prediction accuracies of the extra-large-sized fruit and small-sized fruit yields were higher than those of the total and red-colored fruit yields, as we tested the generated breeding values in inbred tomatoes and F2 populations. Improved accuracy in GEBV calculation of evaluated traits was achieved by using 364 SNPs identified using the regularization methods.

CONCLUSIONS

Together, these results provide an understanding of the genetic variations underlying the heritable phenotypic variability in yield in contemporary tomato breeding and the information necessary for improving such economically important and complex quantitative trait through breeding.

Aspartic proteases modulate programmed cell death and secondary cell wall synthesis during wood formation in poplar

Programmed cell death (PCD) is essential for wood development in trees. However, the determination of crucial factors involved in xylem PCD of wood development is still lacking. Here, two Populus trichocarpa typical aspartic protease (AP) genes, AP17 and AP45, modulate xylem maturation, especially fibre PCD, during wood formation. AP17 and AP45 were dominantly expressed in the fibres of secondary xylem, as suggested by GUS expression in APpro::GUS transgenic plants. Cas9/gRNA-induced AP17 or AP45 mutants delayed secondary xylem fibre PCD, and ap17ap45 double mutants showed more serious defects. Conversely, AP17 overexpression caused premature PCD in secondary xylem fibres, indicating a positive modulation in wood fibre PCD. Loss of AP17 and AP45 did not alter wood fibre wall thickness, whereas the ap17ap45 mutants showed a low lignin content in wood. However, AP17 overexpression led to a significant decrease in wood fibre wall thickness and lignin content, revealing the involvement in secondary cell wall synthesis during wood formation. In addition, the ap17ap45 mutant and AP17 overexpression plants resulted in a significant increase in saccharification yield in wood. Overall, AP17 and AP45 are crucial modulators in xylem maturation during wood development, providing potential candidate genes for engineering lignocellulosic wood for biofuel utilization.

Lectin affinity-based glycoproteome analysis of the developing xylem in poplar

Glycosylation is a significant post-translational modification of proteins, and some glycoproteins serve as players in plant cell wall synthesis and modification. Wood is a highly developed cell wall organization, and protein glycosylation as a regulatory mechanism may be involved in wood formation. Here, a lectin affinity-based glycoproteome was performed in stem developing xylem of poplar. After enrichment, trypsin digestion, LC-MS/MS analysis and peptide identification, we identified 154 glycoproteins from poplar developing xylem, which were classified into nine functional groups mainly including protein acting on carbohydrates, oxido-reductase, proteases, and protein kinases. Further, N- and/or O-glycosylation sites of the identified proteins were analyzed using bioinformatic tools, and deglycosylation experiments in the selected PtSOD and PtHAD proteins verified the reliability of the identified glycoproteins. Analysis of protein subcellular localization showed that a total of 63% of the identified glycoproteins were extracellular proteins or located in the plasma membrane. Poplar eFP and RT-qPCR data showed that a number of the genes encoding these glycoproteins such as laccase, peroxidase and cysteine protease, have highly preferential expression profiles in the developing xylem. Together with previously published research, most identified glycoproteins could be involved in wood cell wall synthesis and modification in poplar. Thus, our study provides some potential wood formation-related glycoproteins to be determined during tree stem development.

Genome-wide identification, structural analysis and expression profiles of short internodes related sequence gene family in quinoa

Based on the whole genome data information of Chenopodium quinoa Willd, the CqSRS gene family members were systematically identified and analyzed by bioinformatics methods, and the responses of CqSRS genes to NaCl (100 mmol/L), salicylic acid (200 umol/L) and low temperature (4°C) were detected by qRT-PCR. The results showed that a total of 10 SHI related sequence genes were identified in quinoa, and they were distributed on 9 chromosomes, and there were four pairs of duplicated genes. The number of amino acids encoded ranged from 143 aa to 370 aa, and the isoelectric point ranged from 4.81 to 8.90. The secondary structure was mainly composed of random coil (Cc). Most of the SRS gene encoding proteins were located in the cytoplasm (5 CqSRS). Phylogenetic analysis showed that the CqSRS genes were divided into three groups, and the gene structure showed that the number of exons of CqSRS was between two-five. Promoter analysis revealed that there are a total of 44 elements related to plant hormone response elements, light response elements, stress response elements and tissue-specific expression in the upstream regin of the gene. Protein interaction showed that all 10 CqSRS proteins appeared in the known protein interaction network diagram in Arabidopsis. Expression profile analysis showed that CqSRS genes had different expression patterns, and some genes had tissue-specific expression. qRT-PCR showed that all SRS family genes responded to ABA、NaCl、drought and low-temperature treatments, but the expression levels of different CqSRS genes were significantly different under various stresses. This study lays a foundation for further analyzed the function of CqSRS genes.

The grapevine aspartic protease gene family: characterization and expression modulation in response to Plasmopara viticola

Grapevine aspartic proteases gene family is characterized and five VviAPs appear to be involved in grapevine defense against downy mildew. Grapevine (Vitis vinifera L.) is one of the most important crops worldwide. However, it is highly susceptible to the downy mildew disease caused by Plasmopara viticola (Berk. & Curt.) Berl. & De Toni. To minimize the use of fungicides used to control P. viticola, it is essential to gain a deeper comprehension on this pathosystem and proteases have gained particular interest in the past decade. Proteases were shown to actively participate in plant-pathogen interactions, not only in the processes that lead to plant cell death, stress responses and protein processing/degradation but also as components of the recognition and signalling pathways. The aim of this study was to identify and characterize the aspartic proteases (APs) involvement in grapevine defense against P. viticola. A genome-wide search and bioinformatics characterization of the V. vinifera AP gene family was conducted and a total of 81 APs proteins, coded by 65 genes, were found. VviAPs proteins can be divided into three categories, similar to those previously described for other plants. Twelve APs coding genes were selected, and expression analysis was conducted at several time-points after inoculation in both compatible and incompatible interactions. Five grapevine APs may be involved in grapevine tolerance against P. viticola. Our findings provide an overall understanding of the VviAPs gene family and establish better groundwork to further describe the roles of VviAPs in defense against P. viticola.

Genome-Wide Analyses of Aspartic Proteases on Potato Genome (Solanum tuberosum): Generating New Tools to Improve the Resistance of Plants to Abiotic Stress

Aspartic proteases are proteolytic enzymes widely distributed in living organisms and viruses. Although they have been extensively studied in many plant species, they are poorly described in potatoes. The present study aimed to identify and characterize S. tuberosum aspartic proteases. Gene structure, chromosome and protein domain organization, phylogeny, and subcellular predicted localization were analyzed and integrated with RNAseq data from different tissues, organs, and conditions focused on abiotic stress. Sixty-two aspartic protease genes were retrieved from the potato genome, distributed in 12 chromosomes. A high number of intronless genes and segmental and tandem duplications were detected. Phylogenetic analysis revealed eight StAP groups, named from StAPI to StAPVIII, that were differentiated into typical (StAPI), nucellin-like (StAPIIIa), and atypical aspartic proteases (StAPII, StAPIIIb to StAPVIII). RNAseq data analyses showed that gene expression was consistent with the presence of cis-acting regulatory elements on StAP promoter regions related to water deficit. The study presents the first identification and characterization of 62 aspartic protease genes and proteins on the potato genome and provides the baseline material for functional gene determinations and potato breeding programs, including gene editing mediated by CRISPR.

Genome-wide identification studies - A primer to explore new genes in plant species

Genome data have accumulated rapidly in recent years, doubling roughly after every 6 months due to the influx of next-generation sequencing technologies. A plethora of plant genomes are available in comprehensive public databases. This easy access to data provides an opportunity to explore genome datasets and recruit new genes in various plant species not possible a decade ago. In the past few years, many gene families have been published using these public datasets. These genome-wide studies identify and characterize gene members, gene structures, evolutionary relationships, expression patterns, protein interactions and gene ontologies, and predict putative gene functions using various computational tools. Such studies provide meaningful information and an initial framework for further functional elucidation. This review provides a concise layout of approaches used in these gene family studies and demonstrates an outline for employing various plant genome datasets in future studies.

Transcriptional landscape of highly lignified poplar stems at single-cell resolution

BACKGROUND

Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood.

RESULTS

Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees.

CONCLUSIONS

These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees.

Defense and Offense Strategies: The Role of Aspartic Proteases in Plant-Pathogen Interactions

Plant aspartic proteases (APs; E.C.3.4.23) are a group of proteolytic enzymes widely distributed among different species characterized by the conserved sequence Asp-Gly-Thr at the active site. With a broad spectrum of biological roles, plant APs are suggested to undergo functional specialization and to be crucial in developmental processes, such as in both biotic and abiotic stress responses. Over the last decade, an increasing number of publications highlighted the APs' involvement in plant defense responses against a diversity of stresses. In contrast, few studies regarding pathogen-secreted APs and AP inhibitors have been published so far. In this review, we provide a comprehensive picture of aspartic proteases from plant and pathogenic origins, focusing on their relevance and participation in defense and offense strategies in plant-pathogen interactions.

Genome-wide identification, evolution and expression analysis of the aspartic protease gene family during rapid growth of moso bamboo (Phyllostachys edulis) shoots

BACKGROUND

Aspartic proteases (APs) are a class of aspartic peptidases belonging to nine proteolytic enzyme families whose members are widely distributed in biological organisms. APs play essential functions during plant development and environmental adaptation. However, there are few reports about APs in fast-growing moso bamboo.

RESULT

In this study, we identified a total of 129 AP proteins (PhAPs) encoded by the moso bamboo genome. Phylogenetic and gene structure analyses showed that these 129 PhAPs could be divided into three categories (categories A, B and C). The PhAP gene family in moso bamboo may have undergone gene expansion, especially the members of categories A and B, although homologs of some members in category C have been lost. The chromosomal location of PhAPs suggested that segmental and tandem duplication events were critical for PhAP gene expansion. Promoter analysis revealed that PhAPs in moso bamboo may be involved in plant development and responses to environmental stress. Furthermore, PhAPs showed tissue-specific expression patterns and may play important roles in rapid growth, including programmed cell death, cell division and elongation, by integrating environmental signals such as light and gibberellin signals.

CONCLUSION

Comprehensive analysis of the AP gene family in moso bamboo suggests that PhAPs have experienced gene expansion that is distinct from that in rice and may play an important role in moso bamboo organ development and rapid growth. Our results provide a direction and lay a foundation for further analysis of plant AP genes to clarify their function during rapid growth.

Genome-wide identification of the aspartic protease gene family and their response under powdery mildew stress in wheat

Aspartic proteases (APs) are one of the four main protease super families. In plants, they are involved in many biological processes, such as biotic and abiotic stress resistance, protein processing and degradation, senescence, and programmed cell death. By performing a database (TGACv1) search and domain prediction, we identified 263 wheat AP (TaAP) proteins and observed 38 TaAP genes exhibiting alternative splicing. Moreover, by constructing a phylogenetic tree, we found that the TaAP proteins can be divided into three families and have a certain close evolutionary relationship to Arabidopsis thaliana and rice AP proteins. Transcriptome analysis showed that 29 genes in the TaAP family were up-regulated after being induced by powdery mildew. The expression of TaAP224 showed the most significant difference in transcriptome and qRT-PCR analyses. Subsequently, the promoters of these 29 genes were analysed, and we found that they contained multiple disease resistance and hormone elements, such as WRKY71OS, a common disease resistance element that is also involved in the GA signalling pathway and inhibits starch hydrolysis. The comprehensive annotation and expression profiling performed in this study increased our understanding of the TaAP family genes in wheat growth and development, and the results can be used as a basis for further study of candidate TaAP genes involved in powdery mildew resistance mechanisms.