Next Article in Journal
Transcriptome Analysis of Roots from Wheat (Triticum aestivum L.) Varieties in Response to Drought Stress
Next Article in Special Issue
First Draft Genome Assembly of Root-Lesion Nematode Pratylenchus scribneri Generated Using Long-Read Sequencing
Previous Article in Journal
Developmental Stage-Dependent Changes in Mitochondrial Function in the Brain of Offspring Following Prenatal Maternal Immune Activation
Previous Article in Special Issue
Identification and Characterization of PRE Genes in Moso Bamboo (Phyllostachys edulis)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systematical Characterization of the AT-Hook Gene Family in Juglans regia L. and the Functional Analysis of the JrAHL2 in Flower Induction and Hypocotyl Elongation

College of Forestry, Hebei Agricultural University, Baoding 071000, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(8), 7244; https://doi.org/10.3390/ijms24087244
Submission received: 9 March 2023 / Revised: 31 March 2023 / Accepted: 3 April 2023 / Published: 14 April 2023

Abstract

:
AT-hook motif nuclear localization (AHL) proteins play essential roles in various plant biological processes. Yet, a comprehensive understanding of AHL transcription factors in walnut (Juglans regia L.) is missing. In this study, 37 AHL gene family members were first identified in the walnut genome. Based on the evolutionary analysis, JrAHL genes were grouped into two clades, and their expansion may occur due to segmental duplication. The stress-responsive nature and driving of developmental activities of JrAHL genes were revealed by cis-acting elements and transcriptomic data, respectively. Tissue-specific expression analysis showed that JrAHLs had a profound transcription in flower and shoot tip, JrAHL2 in particular. Subcellular localization showed that JrAHL2 is anchored to the nucleus. Overexpression of JrAHL2 in Arabidopsis adversely affected hypocotyl elongation and delayed flowering. Our study, for the first time, presented a detailed analysis of JrAHL genes in walnut and provided theoretical knowledge for future genetic breeding programs.

1. Introduction

The AT-hook motif nuclear localized (AHL) gene family encodes the conserved transcription factor that widely distributed in plants species [1]. The AHL proteins are generally distinguished by two characteristic domains: the AT-hook motif and the Plant and Prokaryote Conserved (PPC) domain. The AT-hook motif was first described as a high mobility group (HMG) non-histone chromosomal protein HMG-I/Y [2]. The AT-hook motif contains conserved core sequence ‘Arg-Gly-Arg’ that can bind to the minor groove of the AT-rich B form of DNA [3]. AT-hook motifs can be divided into two categories according to the conserved amino acid sequence: Type-I AT-hook motif contains ‘Gly-Ser-Lys-Asn-Lys’ consensus sequences at the end of the ‘Arg-Gly-Arg’ core, while the Type-II contains ‘Arg-Lys-Tyr.’ The PPC domain, also known as the Domain of Unknown Function #296 (DUF296), is also found in, archaea, prokaryotes, and higher plants. The PPC domain was involved in the physical interaction with other AHL (or themselves) or nuclear proteins [1,4]. According to sequence differences, the PPC can be divided into two categories: Type-A and Type-B. Type-B contains the conserved core ‘Phe-Thr-Pro-His’ in the upstream region, while the core sequence within the PPC domain in Type-A is more variable (‘Leu-Arg-Ser-His’ or ‘Leu-Arg-Ala-His’, etc.) [5]. The AHL genes in land plants have been grouped into Clades based on the phylogenetic relationship: Clades-A and Clades-B. The Clade-A AHLs are intron-free and encode proteins containing only one Type-I AT-hook motif and Type-A PPC domain. The Clade-B AHLs contain introns, and Clade-B AHL proteins can be further divided into two types: one type of AHLs contains a Type-II AT-hook motif and a Type-B PPC domain, and the other type of AHLs contains two AT-hooks (Type-I and Type-II) and a Type-B PPC domain [6].
The AHL gene family plays vital roles in plant growth, development, and stimulus response. In Arabidopsis (Arabidopsis thaliana), the biological functions of several AHLs have been extensively studied. A number of AHL genes redundantly inhibit the elongation of hypocotyl, such as AHL22 [7], SOB3/AHL29 [6,8], and ESC/AHL27 [9]. AHL18, which is involved in the regulation of the length of the proliferation domain and the division in the root apical meristem, and inhibiting primary root growth and lateral root development [10]. In M. truncatula, the MtAHL1 and MtAHL2 were involved in the formation of nitrogen-fixing nodules [11].
The AHL gene is related to regulating flowering time and reproductive development differently. The TEK/AHL16 directly binds to the floral suppressor gene FLC’s repressive regulatory region, leading to its transcriptional silencing. Knockout of the TEK/AHL16 activates a variety of transposons, and up-regulate the expression of floral suppressor gene FLC through epigenetic regulation (histone acetylation and DNA methylation), thereby delaying flowering [12]. TEK/AHL16 also contributes to nexine formation and pollen wall development by directly binding to the nuclear matrix attachment region (MAR) and the promoter of AGP6, an arabinogalactan protein coding gene [13,14]. In addition, overexpressing AHL delays flowering by inhibiting FT expression [7,15]. AHL22 binds to the AT-rich sequence in the FT locus and then alters the structural design of the chromatin region by regulating both histone acetylation and methylation [15]. AHL18, which is in the same phylogenetic clade as AHL22, conferred a similar late-flowering phenotype when overexpression in Arabidopsis to that of AHL22 [7]. AHLs were also reported to be involved in hormone-related pathways. SOB3/AHL29 inhibited PIF4 protein accumulation, thereby reducing the binding of PIF4 to target loci involved in hormone transport, homeostasis, or response such as auxin, brassinosteroid, and ethylene [16]. SOB3/AHL29 also targeted downstream genes associated with auxin-related genes, such as YUCCA8 and SAUR19 [8]. AGF1/AHL25 regulated gibberellins homeostasis in Arabidopsis by targeting the GA 3-oxidase gene AtGA3ox1 [17]. In addition, AHLs participated in the biotic and abiotic stress responses. AHL13, AHL15, AHL19, AHL20, and AHL27 were reported to implicate in PAMP-triggered immunity [18,19,20]. Overexpression of AHL20 and its closely related family members (AHL15, AHL19, and AHL27) in protoplasts inhibited the innate immune response [19]. Especially in Arabidopsis, AHL13 regulated key factors in the jasmonic acid (JA) pathway and affected immunity toward Botrytis cinerea and Pseudomonas syringae pathogens [20]. AHL10 played a critical in response to low-water potential stress, and in the regulation of auxin- and jasmonic acid -related gene expression [21]. AtAHL4 was involved in the regulation of lipid catabolism during seedling establishment [22] and of root xylem development [4]. Interactions between AHL proteins are also crucial for their functions. In Arabidopsis roots, the AHL3 physically interacted with AHL4 and facilitated intercellular trafficking [23]. In maize (Zea mays), the AT-hook protein BAF1 is capable of forming homodimers and heterodimers with other members of the AT-hook family, controlling the formation of maize ears [24].
In recent years, the AHL gene family has been identified and analyzed in several species, such as Arabidopsis, soybean (Glycine max) [25], cotton (Gossypium raimondii) [26], Brassica napus [27], grape (Vitis vinifera) [28], carrot (Daucus carota subsp. sativus) [29], and so on. Walnut (Juglans regia L.) is an economically vital nut trees in China and other temperate parts of the world, but the long period before flowering and bearing fruiting after seed germination limits the progress of walnut breeding work. In addition, functional characterization of AHL genes in walnut is still missing. In this study, the AHL gene family was characterized in walnut based on the high-quality genome data [30]. The chromosomal locations, phylogenetic relationships, gene expansions and structures, conserved domains, putative protein interaction, and cis-elements in the promoter region were analyzed. Furthermore, the JrAHL2, a Clade-A member highly expressed in flower, was isolated, its protein subcellular localization was analyzed, and its functions involved in floral transition and hypocotyl elongation were dissected. The current study will provide a theoretical basis for further investigation of the functions of the JrAHL gene family, and an improved understanding of the molecular mechanisms underlying its role in regulating growth and development in walnut.

2. Results

2.1. AHL Genes in Walnut

The HMMscan was executed against the protein databases using the PPC/DUF296 (PF03479) and AT-hook as queries. After manual checking, the 37 unique AHL genes were obtained from walnut (Juglans regia L.). The JrAHL genes were named according to their chromosomal locations (JrAHL1JrAHL37) (Table 1). The gene lengths of JrAHL genes ranged from 667 bp to 7791 bp, which encode polypeptides from 221 to 406 amino acids with predicted molecular weights ranging from 20.04 kDa to 42.63 kDa. The theoretical pI ranged from 5.08 to 10.61. Approximately 92% of JrAHL proteins (except for JrAHL10/11/33) were predicted to be located in the nucleus (Supplementary Table S2), consistent with the previous report [1].

2.2. Phylogenetic Analysis of JrAHL Genes

To infer the evolutionary relationship among the JrAHL proteins, phylogenetic analysis was performed on the full-length protein sequences. The result showed that the JrAHL proteins in walnut could be divided into two clades, as in Arabidopsis (Figure 1)-Clade-A and Clade-B. JrAHL were further classified into Type-I, Type-II, and Type-III (Figure 2 left panel), based on the presence and combination of the two characteristic functional units-the PPC domain and the AT-hook motif (Figure 2 right panel). Clade-A contained 20 JrAHLs, and it was also classified into Type-I, harboring the Clade-A PPC domain and the Type-I AT-hook motifs in their putative protein sequence (Figure 2 right panel). The higher abundance of Type I members (54%) in walnut is also consistent with observations in other land plants [25,27,31]. There were 17 JrAHL members in Clade-B, which all harbored the Clade-B PPC domain. Clade-B was further divided into two types: 11 Type-II members and 6 Type-III members. Type-II members contained Type-I and Type-II AT-hook motifs, whereas the Type-III JrAHLs only contained the Type-II AT-hook motif (Figure 2 right panel).

2.3. Gene Structure and Duplication of JrAHLs

The distribution of introns and exons was investigated to explore the diversity of gene structure. The structures of JrAHL genes could be divided into two types, one contained no or only one intron, and the other contained multiple introns (Figure 3 left panel). All Type-II and Type-III JrAHLs contained four to five introns, while the JrAHLs in Type-I were intronless except for JrAHL21 and JrAHL26. Furthermore, the number and arrangement of introns are relatively conserved within the same clade, implicating the evolutionary similarity between these members. A total of seven conserved motifs in which two motifs (motif 1 and motif 2) have highly conserved ‘R-G-R-P’ amino acid sequences (Figure 3 right panel), which are the typical characteristic sequences of the AT-hook motif family.
Gene duplication leads to the expansion of AHL family through segmental and tandem duplications [1]. We speculated the exact mechanism has driven the expansion of the AHL gene family in walnut. We analyzed gene duplication events in the JrAHL gene family. As shown in Figure 4, a total of 37 JrAHL genes were unevenly mapped onto 16 chromosomes of walnut. Chromosomes 00 and 14 did not contain duplicated genes, whereas Chromosome 12 had the highest number of duplications. More than twenty pairs of JrAHL genes, such as JrAHL1/6, JrAHL4/9, JrAHL24/34, and JrAHL31/25, were located in duplicated genomic regions. Duplicated pairs showed similar motifs (Figure 2), patterns of gene structure (Figure 3), and clustered tightly in the evolutionary tree (Figure 1), suggesting the JrAHL genes have duplicated segmentally during the evolution.

2.4. Cis-Acting Element Analysis of JrAHLs Promoter

The promoter is critical for transcriptional regulation because of its responsiveness to specific hormone levels and environmental cues and plays a vital role in the biological regulation of gene expression under stress. We investigated the possible regulatory mechanisms by searching the promoter region of 1500 bp upstream of the JrAHL gene family. Four categories of cis-elements were identified related to hormone response, stress response, tissue-specific expression, and transcription factor binding (Figure 5). Hormone response elements, such as gibberellin response element (GARE-motif), salicylic acid response element (TCA-element), abscisic acid-responsive element (ABRE), auxin-responsive (AUXRR-core), and MeJA-responsiveness (CGTCA-motif) were abundant in JrAHL promoter regions. Stress response elements, including defense responsive element (TC-rich repeat), low-temperature response (LTR), and drought response element (MBS), were discovered in most of the promoter regions of JrAHL genes. As for cis-acting regulatory elements related to tissue-specific expression, endosperm- and meristem-expression elements were found in 4 and 14 JrAHL promoter regions, respectively (Figure 5 and Figure S2). Four JrAHL genes contained the AT-rich element, which provided the binding site of ATBP-1, a class of AT-rich DNA binding proteins. These results imply that JrAHLs are involved in many biological processes.

2.5. Protein-Protein Interaction Analysis JrAHLs

Members of the AHL family are partially involved in various plant biological processes via protein-protein interaction [23]. We constructed the protein interaction networks between the JrAHL proteins. Consistent with the previous report [16,24], the JrAHLs are shown to interact with each other (Figure 6). Additionally, JrAHL3/JrAHL8/JrAHL11/JrAHL13/JrAHL27/JrAHL33 were predicted in the networks, which were clustered with the physically-interacting proteins AHL3/AHL4 from Arabidopsis [23]. In Poncirus trifoliata, both PtrAHL14 and PtrAHL17 could interact with other proteins to target the same downstream genes [32]. However, the interaction between JrAHL proteins must be further verified in vitro and in vivo.

2.6. Gene Expression Patterns of JrAHLs

The JrAHLs expression profiles of F26 (anthracnose-resistant) and F423 (anthracnose-susceptible) fruit bracts at five-time points after infection were analyzed based on the transcriptomic data. In the anthracnose-resistant F26, more than half of the JrAHLs gene members were up-regulated within 24 hours after inoculation. However, the up-regulation of these JrAHL genes was delayed in the susceptible F423 (Figure 7). This result suggested that JrAHLs are involved in the immune regulation of pathogens. Expression profiles of JrAHLs at different developmental stages of endopleura, one of the walnut’s most commercially essential organs, were analyzed. More than half of JrAHL genes showed high transcription levels in the early endopleura development stage (Figure 8). The tissue-specific expression profile of some JrAHL genes from different clades was determined by real-time quantitative PCR (Figure 9). JrAHL1, JrAHL2, JrAHL7, JrAHL19, and JrAHL28 were highly expressed in the shoot tip. JrAHL1, JrAHL2, and JrAHL7 also showed high expression levels in pistillate flowers. Two genes were preferentially expressed in the vegetative organs, such as JrAHL10 in root and JrAHL33 in leaf. The JrAHL11 is highly expressed in fruit. In addition, JrAHL14 and JrAHL31 showed no tissue-specific expression.

2.7. Gene Cloning and Subcellular Localization of JrAHL2

JrAHL2 was cloned from walnut and found to contain a 927-bp CDS (Figure 10A) that encodes a protein of 307 amino acids. Multiple sequence alignment analysis results showed that the JrAHL2 shared 59.75% sequence identity with AtAHL22. The amino acid sequence within the AT-hook motif and the large PPC/DUF296 domain are highly conserved (Figure 10B). JrAHL2 protein is predicted to be localized in the nucleus (Supplementary Table S2). JrAHL2 was fused to a green fluorescent protein (EGFP) (Figure 10C) and transiently expressed in tobacco leaf cells. The GFP signal was targeted to the nucleus (Figure 10D), suggesting that JrAHL2 functions as a transcription factor.

2.8. Functional Analysis of JrAHL2 in Transgenic Arabidopsis

To dissect the function of JrAHL2, we generated JrAHL2-overexpressing (OE) transgenic Arabidopsis plants. The homozygous JrAHL2-OE plants flowered significantly later than the WT controls (Figure 11A,B). Additionally, the flowering genes AtFT is expressed at a lower level in the JrAHL2-OE Arabidopsis plants (Figure 11C). Overexpression of JrAHL2 also leads to a dominant-negative long hypocotyl phenotype under both long day and dark conditions (Figure 12). Thus, JrAHL2 inhibits floral transition and hypocotyl elongation.

3. Discussion

3.1. Identification and Analysis of JrAHL Gene Family Members

The walnut (Juglans regia L.) is one of the most economically essential nut trees, and the research on walnut focused on flowering and molecular breeding in recent years. However, less attention was paid to walnut than other plants because of the problems in the walnut industry development, such as the long juvenile phase, difficulty in flowering, and hypocotyl elongation through the husk. The high-quality walnut genome has been published [30], providing a reference for systematic research on the genetic composition and potential function. To date, some gene families from walnut have been identified, such as the Basic leucine zipper (bZIP) [33] and General Regulatory Factor (GRF) [34] genes, which may be involved in the regulation of flowering. However, the role of the AT-Hook Motif Containing Nuclear Localized (AHL) gene family has not been explored in higher plants beyond herbaceous plants. In this work, 37 JrAHL members were identified in the walnut genome (Table 1). More member was identified in walnut compared to Arabidopsis (29 members) [6], probably because the walnut genome (540 Mb) is larger than that of Arabidopsis (125 Mb). The JrAHL members were divided into two Clades, based on the phylogenetic tree (Figure 1) and PPC/DUF domain (Figure 2) as reported in other species [1]. AHL proteins interact with each other through their PPC domain [35]. Thus, the variability increase in PPC domains may extend the range of biological functions. The Clade-B is further divided into two types based on the AT-hook motif, Type-II and Type-III. At the same time, the Clade-A is referred to as Type-I (Figure 2). Because the AT-hook motif is responsible for recognizing and binding a specific nucleotide sequence as well as changing the chromatin structure [36], it is, of course, possible that JrAHLs from the different branch can target different downstream genes. The sequences of conserved motifs (Figure 2), the gene structure (Figure 3), and the arrangement of distinct domains (Figure 3) and are highly consistent within each evolutionary branch. High consistency may be caused by segmental duplications (Figure 4) or the whole-genome duplication process [27], which are the major driving forces of AHL family evolution. Similarly, domain-swapping analyses showed that the N-terminus of AHL4 is required for its DNA-binding activity [4]. All JrAHL genes in Clade-B contain multiple exons, whereas all the Clade-A members are intronless (Figure 3). Introns are ubiquitous in eukaryotic genomes, which affect gene function by affecting gene expression regulation [37]. For example, the number of introns in soybeans certainly affects the expression of AHL genes [25]. In addition to regulating the rate of gene evolution [38], introns also provide specific functional elements for selective splicing and exon shuffling, and some functions of non-coding DNA [39], and in so doing allow the occurrence of diversity in JrAHLs protein accumulation, structures (Supplementary Figure S1), and physically interactions (Figure 6). Because the eukaryote genomes gained introns during or after prokaryote-eukaryote divergence [13], an early divergence occurred between the Type-I JrAHLs and the other two types (II and III) of JrAHL genes, and the JrAHLs of Clade-B in the walnut genome have diverged from Clade-A. Gene duplication, one of the most important features of plant genome structure, contributed to the evolution of novel functions [40] (Panchy et al., 2016). The ratio of non-synonymous to synonymous mutations in genes (Ka/Ks) of duplicated gene pairs was < 1, suggesting that purifying selection was the main source [40] of evolution for the JrAHL gene family. It is thus plausible that introns acquisition may cause the structural and functional changes of JrAHLs during duplication (Figure 3 and Figure 4).
In addition to introns, the promoter, another important non-coding region, plays a critical role in regulating gene transcription. The cis-acting elements of some duplicated gene pairs are not the same, such as JrAHL1/6 and JrAHL4/9. These differences would allow the regulation of different aspects of plant development and the responses to hormones and stresses. Studies of promoters, especially cis-acting elements (CREs) are crucial for improving our fundamental understanding of gene regulation and function towards different environmental cues and stress. The CREs that regulate plant development and confer resistance stresses are found in the JrAHL promoters (Figure 5 and Figure S2). Stress-related hormone response elements, including ABA, MeJA, and SA, are widely distributed on multiple JrAHL promoters. Fifteen JrAHLs contained GARE-motif, especially JrAHL6, contains up to four of these elements, suggesting a cross-talk between JrAHLs and GA hormones, such as negative feedback regulation as previously reported [17]. Aux-responsive CREs are predicted in ten JrAHL promoters. Specifically, JrAHL28 is homologous to AHL10 in Arabidopsis, which regulates genes related to the auxin signaling pathway, suggesting a cross-talk between JrAHLs and auxin [8,21]. AHL also negatively regulates jasmonic acid (JA) levels and affects the gene expression in the JA pathway [20]. Notably, two AT-Hook proteins, PtrAHL14/17 in Poncirus trifoliata, were reported to modulate cold tolerance [32]. Adversity stress-related CREs, including low temperature, drought, and defense, suggest that JrAHLs are relevant to controlling these processes. Four JrAHL genes contained the AT-rich element, which provided the binding site of AT-rich DNA binding protein (ATBP-1) [41]. It is plausible that there is reciprocal transcriptional regulation among the JrAHLs family members, in addition to the predicted protein-protein interactions (Figure 6). In addition, tissue-specific elements such as endosperm- and meristem-expression (Figure 5 and Figure S2) suggested that JrAHLs regulate embryogenesis and meristem maturation like previously reported [35,42]. Specifically, JrAHL17 contains both endosperm- and meristem-expression CREs and is homologous to AtAHL10 from Arabidopsis, which affects the expression of the transcription factor STM (Shootmeristemless) required for meristem maintenance [21]. qRT-PCR revealed that the JrAHL genes containing meristem-expression tissue-specific elements in the promoter were highly expressed in the shoot tip (Figure 5 and Figure 9). Transcriptome data facilitate a comprehensive understanding of expression levels of multiple genes. Multiple JrAHL genes showed differential expression by analyzing the transcriptome data of two walnut materials after inoculation with Colletotrichum gloeosporioides. In response to C. gloeosporioides infection, more than half of the JrAHL genes were rapidly activated in F26 at the initial stage of inoculation. At the same time, the activation of JrAHL in F234 was relatively delayed, suggesting that JrAHL is involved in the immune regulation of pathogens and that early expression may confer high resistance to anthracnose in F26 fruits (Figure 7). In fact, there is already established evidence that AHL is involved in plant immunity and positively regulates plant resistance to pathogen infection [20,43]. More than half of the JrAHL genes are highly expressed at the beginning of endocarp development (Figure 8), including JrAHL11, which showed a profound expression in walnut fruit (Figure 9). JrAHL11 was clustered into the same clade with the Arabidopsis AHL4, and the latter is involved in the regulation of lipid mobilization and fatty acid β-oxidation [22], suggesting that the JrAHL11 plays an important role in fruit development. The three detected Type I JrAHL genes, JrAHL1, JrAHL2, and JrAHL7, were highly expressed in pistillate flowers and shoot tip (Figure 9). JrAHL1 was homologous to the AHL20, JrAHL2 and JrAHL7 were homologous to the AHL22 (Figure 1) in Arabidopsis, respectively. Both AHL20 and AHL22 were involved in controlling flowering time in Arabidopsis [7,15,44], suggesting that these three genes are involved in floral regulation.

3.2. JrAHL2 Delayed Flowering and Inhibited Hypocotyl Elongation

Hypocotyl elongation is an essential event in the photomorphogenesis of walnut seeds, and flowering is an important trait for economic forest species. The functions and mechanisms of AHL genes from Arabidopsis have been extensively studied, whereas the function dissection information about the specific JrAHL gene in walnut is still unavailable. The JrAHL2 were highly expressed in the pistillate flower and was homologous to the Arabidopsis AHL22, which stimulates us to investigate the function of JrAHL2 in controlling flowering and hypocotyl elongation. The transgenic Arabidopsis with constitutive expression of the JrAHL2 gene was prepared. Overexpression of JrAHL2 caused a significant inhibition of hypocotyl elongation, whereas the Arabidopsis ahl22 (homologous to JrAHL2) mutant is longer in the hypocotyl length than wide type (Figure 10) under both long-day conditions and darkness. It reported that AHL22, SOB3/AHL29, ESC/AHL27 [6,7,8] are negative modulators of hypocotyl growth in Arabidopsis [9]. This result indicates that JrAHL2 plays a negative role in controlling hypocotyl elongation and functions in both photomorphogensis and skotomorphogenesis. In addition to repressing the hypocotyl growth, JrAHL2 also inhibits the transcription of FT and delays flowering (Figure 9). Nevertheless, no perturbations to flower phenotype and FT expression were recorded in the ahl22 mutant (Figure 9), possibly because of extensive functional redundancy between AtAHL22 and other AHL genes [15]. These results suggested that the JrAHL2 has a much stronger effect on hypocotyl elongation than its effect on flowering regulation (Xiao et al., 2009). Additionally, despite the report that AHL regulates the expression of target genes, including FT, by directly binding to a specific nuclear MAR containing A/T-rich sequences within their promoter [15,32,44], the ubiquity of A/T-rich sequences in the promoter region makes it difficult to discriminate the elements specifically bound by JrAHL2. Thus, a systems approach such as ChIP-seq/qPCR could be used further to dissect the network of JrAHL.

4. Materials and Methods

4.1. Identification of the AHL Gene Family

For the discovery of the AHL family members in walnut, the high-quality common walnut (Juglans regia L.) reference genome data was downloaded [30]. The Hidden Markov Model scan (HMMscan) was used on the complete walnut protein database using PPC/DUF296 (PF03479) as well as AT-hook (PF02178) Pfam IDs obtained from the Pfam database (https://pfam.xfam.org/ (accessed on 7 January 2022)) [45] as queries. All non-redundant putative protein sequences were manually checked with the NCBI Conserved Domains Database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 12 January 2022)). The unique sequences containing the AT-hook motif and the PPC domain were selected as putative JrAHLs. The online ExPASy program (http://www.expasy.org/tools/ (accessed on 12 February 2022)) was used to predict protein physicochemical characteristics. The online system Cell-PLoc (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2/ (accessed on 15 February 2022)) [46] was used to predict the subcellular localization of these JrAHL proteins.

4.2. Phylogenetic Analysis of AHL Proteins

A Neighbor-Joining tree was constructed to infer the phylogenetic relationship between the AHL proteins from Arabidopsis and walnut. Sequence alignment was carried out using ClustalW in MEGA7 [47] with default parameters. The 1000 Bootstrap replicates were applied to present the evolutionary history. Multiple alignments were also performed with DNAMAN software (LynnonBiosoft, San Ramon, CA, USA. version 8.0).

4.3. Genes and Protein Structure Analysis

The GSDS 2.0 (Gene Structure Display Server, http://gsds.cbi.pku.edu.cn) was utilized to draw the structure of JrAHL genes. MEME (http://meme-suite.org/ (accessed on 15 April 2022)) was employed for functional motif analysis. The three-dimensional structure of JrAHL proteins was predicted by the PHYRE server v2.0 (http://www.sbg.bio.ic.ac.uk/phyre2 (accessed on 22 April 2022)).

4.4. Collinearity Analysis

The OrthoMCL algorithm [48] was used to identify paralogous genes within the walnut genome, and then the MCScan [49] algorithm was applied to detect syntenic blocks containing walnut AHL. The Circos [50] was used to visualize the syntenic relationships between the walnut genomes. The Ka/Ks ratio was calculated to estimate the selection pressure among duplicated JrAHL genes as previously described [51]. All plugins are used in TB tool software (TBtools-II v1.108) [52].

4.5. Promoter Sequence Analysis

The 1.5-kb upstream region of the start codon of each JrAHL gene was retrieved. The online PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 16 July 2022)) was used to identify the potential cis-acting elements.

4.6. Interactive Protein Network Analysis

The STRING database predicted the interaction network of AHL proteins (https://string-db.org/ (accessed on 25 July 2022)). All JrAHL proteins were used as the query to search for the walnut data in the STRING database. Both the functional and physical protein associations were included in the network. The minimum required interaction score was selected as medium confidence (0.400).

4.7. Expression Profiles of JrAHL Genes

Gene expression data from the anthracnose-resistant F26 and the anthracnose-susceptible F423 fruits in response to anthracnose (Colletotrichum gloeosporioides) infection were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/ (accessed on 20 July 2022)) at accession GSE147083 [53]. Transcriptomic data of the walnut endopleura during its development was also downloaded from GEO database at accession GSE185230 [54]. The log2 transformed FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) values were used for preparing the heatmaps. Root, stem, leaf, flowers, fruits, and shoot tip were taken, frozen immediately in liquid nitrogen, and stored at −80 °C until use.

4.8. Quantitative RT-PCR (qRT-PCR)

Total RNA was extracted from different tissues using the RNA extraction kit (OMEGA, Doraville, GA, USA). The RNA concentration was determined by spectrophotometer, and RNA integrity was verified by electrophoresis. First-strand cDNA was synthesized using the EasyScript® One-Step cDNA Synthesis SuperMix (TransGen Biotech, China) following the manufacturer’s instructions. Primer pairs for qRT-PCR were designed using Primer Premier 6.0 (Premier Biosoft, San Francisco, CA, USA) (Supplementary Table S1). qRT-PCR was performed on a Bio-Rad CFX96 thermal cycler (Bio-Rad Laboratories, Hercules, CA, USA). The 18S rRNA gene was used as a standard internal control. The relative expression levels were calculated by the comparative 2−ΔΔCt method [55].

4.9. Gene Cloning, Vector Construction, and Subcellular Localization

The coding sequence (CDS) of the JrAHL2 was cloned from the walnut pistillate flower’s cDNA and then inserted into the modified pCAMBIA2300-EGFP vector. Subcellular localization was carried out according to our previously published article [56]. The leaves of 5-week-old tobacco plants were infiltrated with the Agrobacterium strain GV3103 cells containing the recombinant vector. After 3 d, the green fluorescent protein signals were detected with an Olympus fluorescence microscope (Japan). All primers used for gene cloning and vector construction were listed in Supplementary Table S1.

4.10. Plant Transformation and Phenotype Analysis

Arabidopsis was transformed with the Agrobacterium strain GV3103 mentioned above, according to the A. tumefaciens mediated floral dip method [57]. The transformants were identified based on kanamycin resistance, and the T3 generation seedlings were used for functional analyses. The flowering phenotype of different Arabidopsis lines, including JrAHL2-overexpressing (JrAHL2-OE), wild type (Col-0), and AtAHL22 mutant (ahl22-1, Salk_018866), were identified according to the previous description [58]. To analyze the hypocotyl growth phenotype, Arabidopsis seeds were surface-sterilized and sown on Murashige & Skoog (MS) medium [56]. The plate was held at 4 °C for 3 days, exposed to white light for 4 h to stimulate germination, and then cultured for 6 days at 22 °C with a 16-h light/8-h dark photoperiod (light intensity of 150 μmol m−2 s−1) (long day) or darkness.

5. Conclusions

Thirty-seven AHL genes were identified in the walnut (Juglans regia L.) genome and found to be unequally distributed on 16 chromosomes. The phylogenetic tree divided JrAHLs into two clades. Segmental duplications were the major driving forces of AHL gene family expansion and evolution. Promoter and expression analysis indicated that JrAHLs were involved in multiple biological processes. Functional analysis showed that the nuclear-localized JrAHL2 plays a negative role in floral transition and hypocotyl elongation. The first systematic exploration of the AHL gene family in walnut will provide helpful information for future work in walnut breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24087244/s1.

Author Contributions

Conceptualization: P.J. and G.Q.; Data curation, P.J.; Formal analysis, K.Y., Q.D., H.L. (Haoan Luan) and S.G.; Funding acquisition, G.Q.; Investigation, H.L. (Han Li); Methodology, J.L.; Project administration, G.Q.; Resources, H.L. (Haoan Luan); Software, R.Y.; Supervision, G.Q.; Validation, X.Z.; Writing—review & editing, P.J. and G.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The National Key Research and Development Program grant number 2022YFD1600402. And The APC was funded by The National Key Research and Development Program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data are available from the corresponding author on request ([email protected]).

Acknowledgments

The authors would like to thank Rahat Sharif for the technical help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, J.; Favero, D.S.; Qiu, J.; Roalson, E.H.; Neff, M.M. Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants. BMC Plant Biol. 2014, 14, 266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhao, K.; Kas, E.; Gonzalez, E.; Laemmli, U.K. SAR-dependent mobilization of histone H1 by HMG-I/Y in vitro: HMG-I/Y is enriched in H1-depleted chromatin. EMBO J. 1993, 12, 3237–3247. [Google Scholar] [CrossRef] [PubMed]
  3. Huth, J.R.; Bewley, C.A.; Nissen, M.S.; Evans, J.N.; Reeves, R.; Gronenborn, A.M.; Clore, G.M. The solution structure of an HMG-I(Y)-DNA complex defines a new architectural minor groove binding motif. Nat. Struct. Biol. 1997, 4, 657–665. [Google Scholar] [CrossRef] [PubMed]
  4. Seo, M.; Lee, J.Y. Dissection of Functional Modules of AT-HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 4 in the Development of the Root Xylem. Front. Plant Sci. 2021, 12, 632078. [Google Scholar] [CrossRef]
  5. Zhang, W.M.; Cheng, X.Z.; Fang, D.; Cao, J. AT-HOOK MOTIF NUCLEAR LOCALIZED (AHL) proteins of ancient origin radiate new functions. Int. J. Biol. Macromol 2022, 214, 290–300. [Google Scholar] [CrossRef]
  6. Zhao, J.; Favero, D.S.; Peng, H.; Neff, M.M. Arabidopsis thaliana AHL family modulates hypocotyl growth redundantly by interacting with each other via the PPC/DUF296 domain. Proc. Natl. Acad. Sci. USA 2013, 110, E4688–E4697. [Google Scholar] [CrossRef] [Green Version]
  7. Xiao, C.; Chen, F.; Yu, X.; Lin, C.; Fu, Y.F. Over-expression of an AT-hook gene, AHL22, delays flowering and inhibits the elongation of the hypocotyl in Arabidopsis thaliana. Plant Mol. Biol. 2009, 71, 39–50. [Google Scholar] [CrossRef]
  8. Favero, D.S.; Jacques, C.N.; Iwase, A.; Le, K.N.; Zhao, J.; Sugimoto, K.; Neff, M.M. SUPPRESSOR OF PHYTOCHROME B4-#3 Represses Genes Associated with Auxin Signaling to Modulate Hypocotyl Growth. Plant Physiol. 2016, 171, 2701–2716. [Google Scholar] [CrossRef] [Green Version]
  9. Street, I.H.; Shah, P.K.; Smith, A.M.; Avery, N.; Neff, M.M. The AT-hook-containing proteins SOB3/AHL29 and ESC/AHL27 are negative modulators of hypocotyl growth in Arabidopsis. Plant J. 2008, 54, 1–14. [Google Scholar] [CrossRef]
  10. Sirl, M.; Snajdrova, T.; Gutierrez-Alanis, D.; Dubrovsky, J.G.; Vielle-Calzada, J.P.; Kulich, I.; Soukup, A. At-Hook Motif Nuclear Localised Protein 18 as a Novel Modulator of Root System Architecture. Int. J. Mol. Sci. 2020, 21, 1886. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, S.; Wang, T.; Lima, R.M.; Pettko-Szandtner, A.; Kereszt, A.; Downie, J.A.; Kondorosi, E. Widely conserved AHL transcription factors are essential for NCR gene expression and nodule development in Medicago. Nat. Plants 2023, 9, 280–288. [Google Scholar] [CrossRef]
  12. Xu, Y.; Gan, E.S.; Ito, T. The AT-hook/PPC domain protein TEK negatively regulates floral repressors including MAF4 and MAF5. Plant Signal. Behav. 2013, 8, e25006. [Google Scholar] [CrossRef] [Green Version]
  13. Jo, B.S.; Choi, S.S. Introns: The Functional Benefits of Introns in Genomes. Genomics. Inform. 2015, 13, 112–118. [Google Scholar] [CrossRef] [Green Version]
  14. Jia, Q.S.; Zhu, J.; Xu, X.F.; Lou, Y.; Zhang, Z.L.; Zhang, Z.P.; Yang, Z.N. Arabidopsis AT-hook protein TEK positively regulates the expression of arabinogalactan proteins for Nexine formation. Mol. Plant 2015, 8, 251–260. [Google Scholar] [CrossRef] [Green Version]
  15. Yun, J.; Kim, Y.S.; Jung, J.H.; Seo, P.J.; Park, C.M. The AT-hook motif-containing protein AHL22 regulates flowering initiation by modifying FLOWERING LOCUS T chromatin in Arabidopsis. J. Biol. Chem. 2012, 287, 15307–15316. [Google Scholar] [CrossRef] [Green Version]
  16. Favero, D.S.; Kawamura, A.; Shibata, M.; Takebayashi, A.; Jung, J.H.; Suzuki, T.; Jaeger, K.E.; Ishida, T.; Iwase, A.; Wigge, P.A.; et al. AT-Hook Transcription Factors Restrict Petiole Growth by Antagonizing PIFs. Curr. Biol. 2020, 30, 1454–1466.e1456. [Google Scholar] [CrossRef]
  17. Matsushita, A.; Furumoto, T.; Ishida, S.; Takahashi, Y. AGF1, an AT-hook protein, is necessary for the negative feedback of AtGA3ox1 encoding GA 3-oxidase. Plant Physiol. 2007, 143, 1152–1162. [Google Scholar] [CrossRef] [Green Version]
  18. Kim, S.Y.; Kim, Y.C.; Seong, E.S.; Lee, Y.H.; Park, J.M.; Choi, D. The chili pepper CaATL1: An AT-hook motif-containing transcription factor implicated in defence responses against pathogens. Mol. Plant Pathol. 2007, 8, 761–771. [Google Scholar] [CrossRef]
  19. Lu, H.; Zou, Y.; Feng, N. Overexpression of AHL20 negatively regulates defenses in Arabidopsis. J. Integr. Plant Biol. 2010, 52, 801–808. [Google Scholar] [CrossRef]
  20. Rayapuram, N.; Jarad, M.; Alhoraibi, H.M.; Bigeard, J.; Abulfaraj, A.A.; Volz, R.; Mariappan, K.G.; Almeida-Trapp, M.; Schloffel, M.; Lastrucci, E.; et al. Chromatin phosphoproteomics unravels a function for AT-hook motif nuclear localized protein AHL13 in PAMP-triggered immunity. Proc. Natl. Acad. Sci. USA 2021, 118, e2004670118. [Google Scholar] [CrossRef]
  21. Wong, M.M.; Bhaskara, G.B.; Wen, T.N.; Lin, W.D.; Nguyen, T.T.; Chong, G.L.; Verslues, P.E. Phosphoproteomics of Arabidopsis Highly ABA-Induced1 identifies AT-Hook-Like10 phosphorylation required for stress growth regulation. Proc. Natl. Acad. Sci. USA 2019, 116, 2354–2363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Cai, G.; Kim, S.C.; Li, J.; Zhou, Y.; Wang, X. Transcriptional Regulation of Lipid Catabolism during Seedling Establishment. Mol. Plant 2020, 13, 984–1000. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, J.; Wang, X.; Lee, J.Y.; Lee, J.Y. Cell-to-cell movement of two interacting AT-hook factors in Arabidopsis root vascular tissue patterning. Plant Cell 2013, 25, 187–201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gallavotti, A.; Malcomber, S.; Gaines, C.; Stanfield, S.; Whipple, C.; Kellogg, E.; Schmidt, R.J. BARREN STALK FASTIGIATE1 is an AT-hook protein required for the formation of maize ears. Plant Cell 2011, 23, 1756–1771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, M.; Chen, B.; Zhou, W.; Xie, L.; Wang, L.; Zhang, Y.; Zhang, Q. Genome-wide identification and expression analysis of the AT-hook Motif Nuclear Localized gene family in soybean. BMC Genomics 2021, 22, 361. [Google Scholar] [CrossRef]
  26. Zhao, L.; Lu, Y.; Chen, W.; Yao, J.; Li, Y.; Li, Q.; Pan, J.; Fang, S.; Sun, J.; Zhang, Y. Genome-wide identification and analyses of the AHL gene family in cotton (Gossypium). BMC Genom. 2020, 21, 69. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, W.M.; Fang, D.; Cheng, X.Z.; Cao, J.; Tan, X.L. Insights Into the Molecular Evolution of AT-Hook Motif Nuclear Localization Genes in Brassica napus. Front. Plant Sci. 2021, 12, 714305. [Google Scholar] [CrossRef]
  28. Li, X.; He, H.; Wang, H.; Wu, X.; Wang, H.; Mao, J. Identification and expression analysis of the AHL gene family in grape (Vitis vinifera). Plant Gene 2021, 26, 100285. [Google Scholar] [CrossRef]
  29. Machaj, G.; Grzebelus, D. Characteristics of the AT-Hook Motif Containing Nuclear Localized (AHL) Genes in Carrot Provides Insight into Their Role in Plant Growth and Storage Root Development. Genes 2021, 12, 764. [Google Scholar] [CrossRef]
  30. Zhang, J.; Zhang, W.; Ji, F.; Qiu, J.; Song, X.; Bu, D.; Pan, G.; Ma, Q.; Chen, J.; Huang, R.; et al. A high-quality walnut genome assembly reveals extensive gene expression divergences after whole-genome duplication. Plant Biotechnol. J. 2020, 18, 1848–1850. [Google Scholar] [CrossRef]
  31. Bishop, E.H.; Kumar, R.; Luo, F.; Saski, C.; Sekhon, R.S. Genome-wide identification, expression profiling, and network analysis of AT-hook gene family in maize. Genomics 2020, 112, 1233–1244. [Google Scholar] [CrossRef]
  32. Dahro, B.; Wang, Y.; Khan, M.; Zhang, Y.; Fang, T.; Ming, R.; Li, C.; Liu, J.H. Two AT-Hook proteins regulate A/NINV7 expression to modulate sucrose catabolism for cold tolerance in Poncirus trifoliata. New Phytol. 2022, 235, 2331–2349. [Google Scholar] [CrossRef]
  33. Zhang, Z.; Quan, S.; Niu, J.; Guo, C.; Kang, C.; Liu, J.; Yuan, X. Genome-Wide Identification, Classification, Expression and Duplication Analysis of bZIP Family Genes in Juglans regia L. Int. J. Mol. Sci. 2022, 23, 5961. [Google Scholar] [CrossRef]
  34. Zhang, Z.; Quan, S.; Niu, J.; Guo, C.; Kang, C.; Liu, J.; Yuan, X. Comprehensive Identification and Analyses of the GRF Gene Family in the Whole-Genome of Four Juglandaceae Species. Int. J. Mol. Sci. 2022, 23, 12663. [Google Scholar] [CrossRef]
  35. Karami, O.; Rahimi, A.; Khan, M.; Bemer, M.; Hazarika, R.R.; Mak, P.; Compier, M.; van Noort, V.; Offringa, R. A suppressor of axillary meristem maturation promotes longevity in flowering plants. Nat. Plants 2020, 6, 368–376. [Google Scholar] [CrossRef]
  36. Aravind, L.; Landsman, D. AT-hook motifs identified in a wide variety of DNA-binding proteins. Nucleic Acids Res. 1998, 26, 4413–4421. [Google Scholar] [CrossRef] [Green Version]
  37. Castillo-Davis, C.I.; Mekhedov, S.L.; Hartl, D.L.; Koonin, E.V.; Kondrashov, F.A. Selection for short introns in highly expressed genes. Nat. Genet. 2002, 31, 415–418. [Google Scholar] [CrossRef]
  38. Roy, S.W.; Gilbert, W. The evolution of spliceosomal introns: Patterns, puzzles and progress. Nat. Rev. Genet. 2006, 7, 211–221. [Google Scholar] [CrossRef]
  39. North, K.; Benbarche, S.; Liu, B.; Pangallo, J.; Chen, S.; Stahl, M.; Bewersdorf, J.P.; Stanley, R.F.; Erickson, C.; Cho, H.; et al. Synthetic introns enable splicing factor mutation-dependent targeting of cancer cells. Nat. Biotechnol. 2022, 40, 1103–1113. [Google Scholar] [CrossRef]
  40. Panchy, N.; Lehti-Shiu, M.; Shiu, S.H. Evolution of Gene Duplication in Plants. Plant Physiol. 2016, 171, 2294–2316. [Google Scholar] [CrossRef] [Green Version]
  41. Tjaden, G.; Coruzzi, G.M. A novel AT-rich DNA binding protein that combines an HMG I-like DNA binding domain with a putative transcription domain. Plant Cell 1994, 6, 107–118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Karami, O.; Rahimi, A.; Mak, P.; Horstman, A.; Boutilier, K.; Compier, M.; van der Zaal, B.; Offringa, R. An Arabidopsis AT-hook motif nuclear protein mediates somatic embryogenesis and coinciding genome duplication. Nat. Commun. 2021, 12, 2508. [Google Scholar] [CrossRef] [PubMed]
  43. Liu, H.; Dong, S.; Li, M.; Gu, F.; Yang, G.; Guo, T.; Chen, Z.; Wang, J. The Class III peroxidase gene OsPrx30, transcriptionally modulated by the AT-hook protein OsATH1, mediates rice bacterial blight-induced ROS accumulation. J. Integr Plant Biol. 2021, 63, 393–408. [Google Scholar] [CrossRef] [PubMed]
  44. Tayengwa, R.; Sharma Koirala, P.; Pierce, C.F.; Werner, B.E.; Neff, M.M. Overexpression of AtAHL20 causes delayed flowering in Arabidopsis via repression of FT expression. BMC Plant Biol. 2020, 20, 559. [Google Scholar] [CrossRef]
  45. Mistry, J.; Chuguransky, S.; Williams, L.; Qureshi, M.; Salazar, G.A.; Sonnhammer, E.L.L.; Tosatto, S.C.E.; Paladin, L.; Raj, S.; Richardson, L.J.; et al. Pfam: The protein families database in 2021. Nucleic. Acids Res. 2021, 49, D412–D419. [Google Scholar] [CrossRef]
  46. Chou, K.C.; Shen, H.B. Cell-PLoc: A package of Web servers for predicting subcellular localization of proteins in various organisms. Nat. Protoc. 2008, 3, 153–162. [Google Scholar] [CrossRef]
  47. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  48. Li, L.; Stoeckert, C.J., Jr.; Roos, D.S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 2003, 13, 2178–2189. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  50. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome. Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, Z. KaKs_Calculator 3.0: Calculating Selective Pressure on Coding and Non-coding Sequences. Genom. Proteom. Bioinform. 2022, 20, 536–540. [Google Scholar] [CrossRef]
  52. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  53. Fang, H.; Liu, X.; Dong, Y.; Feng, S.; Zhou, R.; Wang, C.; Ma, X.; Liu, J.; Yang, K.Q. Transcriptome and proteome analysis of walnut (Juglans regia L.) fruit in response to infection by Colletotrichum gloeosporioides. BMC Plant Biol. 2021, 21, 249. [Google Scholar] [CrossRef]
  54. Ji, F.; Ma, Q.; Zhang, W.; Liu, J.; Feng, Y.; Zhao, P.; Song, X.; Chen, J.; Zhang, J.; Wei, X.; et al. A genome variation map provides insights into the genetics of walnut adaptation and agronomic traits. Genome. Biol. 2021, 22, 300. [Google Scholar] [CrossRef]
  55. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  56. Jia, P.; Xing, L.; Zhang, C.; Zhang, D.; Ma, J.; Zhao, C.; Han, M.; Ren, X.; An, N. MdKNOX19, a class II knotted-like transcription factor of apple, plays roles in ABA signalling/sensitivity by targeting ABI5 during organ development. Plant Sci. 2021, 302, 110701. [Google Scholar] [CrossRef]
  57. Zhang, X.; Henriques, R.; Lin, S.S.; Niu, Q.W.; Chua, N.H. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nat. Protoc. 2006, 1, 641–646. [Google Scholar] [CrossRef]
  58. Jia, P.; Xing, L.B.; Zhang, C.G.; Chen, H.; Li, Y.M.; Zhang, D.; Ma, J.J.; Zhao, C.P.; Han, M.Y.; Ren, X.L.; et al. MdKNOX15, a class I knotted-like transcription factor of apple, controls flowering and plant height by regulating GA levels through promoting the MdGA2ox7 transcription. Environ. Exp. Bot. 2021, 185, 104411. [Google Scholar] [CrossRef]
Figure 1. Neighbor-joining tree representing phylogenetic relationships among AHL genes from walnut and Arabidopsis. The evolutionary history was inferred using the Neighbor-Joining matrix-based method and are in the units of the number of amino acid substitutions per site. The rate variation among sites was modeled with a gamma distribution. Evolutionary analyses were conducted in MEGA7.
Figure 1. Neighbor-joining tree representing phylogenetic relationships among AHL genes from walnut and Arabidopsis. The evolutionary history was inferred using the Neighbor-Joining matrix-based method and are in the units of the number of amino acid substitutions per site. The rate variation among sites was modeled with a gamma distribution. Evolutionary analyses were conducted in MEGA7.
Ijms 24 07244 g001
Figure 2. Phylogenetic (left panel) and sequence alignment (right panel) analysis of the walnut AHL proteins. The obtained phylogenetic tree by MEGA7 is shown on the left. Multiple sequence alignment analysis was performed using DNAMAN software, and the conserved domain is displayed on the right.
Figure 2. Phylogenetic (left panel) and sequence alignment (right panel) analysis of the walnut AHL proteins. The obtained phylogenetic tree by MEGA7 is shown on the left. Multiple sequence alignment analysis was performed using DNAMAN software, and the conserved domain is displayed on the right.
Ijms 24 07244 g002
Figure 3. Gene structure analysis and conservative motif prediction of the AT-hook motif gene family in walnut. The gene structure map was obtained with the GSDS online tool (left). The x-axis shows the inferred length of the different genes (5′ to 3′) and their respective CDS (red). Domain structure and motif were elucidated with the MEME platform (right).
Figure 3. Gene structure analysis and conservative motif prediction of the AT-hook motif gene family in walnut. The gene structure map was obtained with the GSDS online tool (left). The x-axis shows the inferred length of the different genes (5′ to 3′) and their respective CDS (red). Domain structure and motif were elucidated with the MEME platform (right).
Ijms 24 07244 g003
Figure 4. Chromosomal distribution and collinearity of walnut AHL genes. The gray lines in the background indicate the collinear blocks in the genomes of J. regia, while the colorful lines emphasize the segmental duplication of the walnut AHLs.
Figure 4. Chromosomal distribution and collinearity of walnut AHL genes. The gray lines in the background indicate the collinear blocks in the genomes of J. regia, while the colorful lines emphasize the segmental duplication of the walnut AHLs.
Ijms 24 07244 g004
Figure 5. Predicted cis-elements in the JrAHL promoters. The 1.5-kb sequence upstream from the start codon of JrAHL genes was analyzed using the PlantCARE database. Different colored graphics represent different cis-acting elements.
Figure 5. Predicted cis-elements in the JrAHL promoters. The 1.5-kb sequence upstream from the start codon of JrAHL genes was analyzed using the PlantCARE database. Different colored graphics represent different cis-acting elements.
Ijms 24 07244 g005
Figure 6. Predicted protein-protein interaction network of JrAHL protein. The network nodes represent proteins. The line color indicates the type of interaction evidence.
Figure 6. Predicted protein-protein interaction network of JrAHL protein. The network nodes represent proteins. The line color indicates the type of interaction evidence.
Ijms 24 07244 g006
Figure 7. JrAHL gene expression profiles in F26 and the F423 fruits in response to anthracnose infection. The heat map was generated based on the transcriptomic data of the anthracnose-resistant F26 fruit bracts and anthracnose-susceptible F423 fruit bracts at five time points after infection with C. gloeosporioides.
Figure 7. JrAHL gene expression profiles in F26 and the F423 fruits in response to anthracnose infection. The heat map was generated based on the transcriptomic data of the anthracnose-resistant F26 fruit bracts and anthracnose-susceptible F423 fruit bracts at five time points after infection with C. gloeosporioides.
Ijms 24 07244 g007
Figure 8. JrAHL gene expression profiles at different developmental stages of endopleura in walnut. The heat map was generated based on the transcriptomic data of the endopleura at different (35, 63, 91, 119, and 147) days after pollination (DAP).
Figure 8. JrAHL gene expression profiles at different developmental stages of endopleura in walnut. The heat map was generated based on the transcriptomic data of the endopleura at different (35, 63, 91, 119, and 147) days after pollination (DAP).
Ijms 24 07244 g008
Figure 9. JrAHL gene expression profiles in different tissues. Each value represents the mean ± standard error of three replicates. Different letters (a, b, c, d) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Figure 9. JrAHL gene expression profiles in different tissues. Each value represents the mean ± standard error of three replicates. Different letters (a, b, c, d) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Ijms 24 07244 g009
Figure 10. JrAHL2 gene cloning and subcellular localization. (A) Gene cloning of JrAHL2 coding region. (B) Putative protein sequence analysis of the cloned JrAHL2 gene. (C) Schematic diagram of the JrAHL2 overexpression vector for subcellular localization. (D) Subcellular localization of JrAHL2 in tobacco leaves.
Figure 10. JrAHL2 gene cloning and subcellular localization. (A) Gene cloning of JrAHL2 coding region. (B) Putative protein sequence analysis of the cloned JrAHL2 gene. (C) Schematic diagram of the JrAHL2 overexpression vector for subcellular localization. (D) Subcellular localization of JrAHL2 in tobacco leaves.
Ijms 24 07244 g010
Figure 11. JrAHL2 affects flowering in Arabidopsis. (A) Flowering phenotype of JrAHL2-overexpressing (JrAHL2-OE), wide type (Col-0), and AtAHL mutant (ahl22-1) Arabidopsis lines, respectively. (B) The number of rosette leaves at flowering. (C) The expression level of flowering-related gene FT. The columns represent the mean values of three replicates ± standard deviations. The asterisk (*) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Figure 11. JrAHL2 affects flowering in Arabidopsis. (A) Flowering phenotype of JrAHL2-overexpressing (JrAHL2-OE), wide type (Col-0), and AtAHL mutant (ahl22-1) Arabidopsis lines, respectively. (B) The number of rosette leaves at flowering. (C) The expression level of flowering-related gene FT. The columns represent the mean values of three replicates ± standard deviations. The asterisk (*) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Ijms 24 07244 g011
Figure 12. JrAHL2 affects hypocotyl elongation. Samples were photographed (A), (upper panel) for darkness, and (lower panel) for normal photoperiod and the length of the hypocotyl was analyzed (B). Bar = 2 mm. The columns represent the mean values of three replicates ± standard deviations. Different letters (a, b, c) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Figure 12. JrAHL2 affects hypocotyl elongation. Samples were photographed (A), (upper panel) for darkness, and (lower panel) for normal photoperiod and the length of the hypocotyl was analyzed (B). Bar = 2 mm. The columns represent the mean values of three replicates ± standard deviations. Different letters (a, b, c) indicate significant differences, as determined by Student’s t-test at the 0.05 level.
Ijms 24 07244 g012
Table 1. Information of the AHL genes in J. regia L.
Table 1. Information of the AHL genes in J. regia L.
NameGene IDLocation (Start)Location (End)StrandLengthCDS (bp)Peptide (aa)pIMW (kDa)
JrAHL1JreChr01G1298439,611,176 39,611,994 +819819 2726.0828.02
JrAHL2JreChr01G1299247,390,686 47,391,612 927927 3086.5932.84
JrAHL3JreChr01G1333347,470,900 47,477,455 +65561005 3349.5934.79
JrAHL4JreChr01G1376150,056,338 50,057,327 990990 3295.0834.44
JrAHL5JreChr02G108987,160,919 7,167,725 68071035 3449.9135.36
JrAHL6JreChr02G1146721,601,290 21,602,144 +855855 2846.4529.38
JrAHL7JreChr02G1147426,318,938 26,319,843 906906 3016.8031.86
JrAHL8JreChr02G1168726,371,944 26,376,879 +49361068 3558.8336.83
JrAHL9JreChr02G1197827,859,904 27,860,884 981981 3265.6034.15
JrAHL10JreChr03G133504,426,812 4,430,966 4155795 26410.4627.04
JrAHL11JreChr04G107418,368,413 8,372,787 +43751014 3379.2135.15
JrAHL12JreChr04G1221020,642,675 20,643,481 +807807 2689.2428.25
JrAHL13JreChr06G109706,003,983 6,011,773 +7791930 3099.6132.44
JrAHL14JreChr07G1038723,875,413 23,876,234 822822 2739.7428.01
JrAHL15JreChr07G1123016,643,389 16,644,297 909909 3025.7731.30
JrAHL16JreChr07G1123927,147,611 27,148,277 667666 22110.0923.04
JrAHL17JreChr07G1161627,252,429 27,256,966 +45381107 36810.0936.49
JrAHL18JreChr08G1192730,383,998 30,384,978 +981981 3267.0734.60
JrAHL19JreChr08G122205,871,308 5,872,264 957957 3186.3633.92
JrAHL20JreChr08G122258,400,985 8,408,164 71801182 39310.3339.80
JrAHL21JreChr09G122398,457,851 8,459,672 +1822696 2319.8723.85
JrAHL22JreChr10G11635644,863 648,220 +33581089 3628.4937.48
JrAHL23JreChr11G10743358,990 361,864 +28751068 3559.1236.87
JrAHL24JreChr11G1135523,940,903 23,941,778 876876 2915.9630.42
JrAHL25JreChr11G1147730,455,540 30,459,915 +43761017 33810.6135.04
JrAHL26JreChr11G1148131,577,477 31,578,376 900900 2996.1632.00
JrAHL27JreChr12G1063631,622,448 31,628,198 +5751939 3129.6331.89
JrAHL28JreChr12G1063818,010,455 18,017,606 +71521134 3779.3738.78
JrAHL29JreChr12G1064018,021,521 18,027,959 64391140 3799.5938.92
JrAHL30JreChr12G1099818,039,195 18,040,085 891891 2966.1031.07
JrAHL31JreChr12G1109221,387,711 21,392,430 +47201020 33910.0035.22
JrAHL32JreChr12G1109522,118,819 22,119,700 882882 2936.3031.42
JrAHL33JreChr13G1046422,173,675 22,179,704 +60301221 4069.5942.63
JrAHL34JreChr13G108461,739,954 1,740,790 837837 2785.8028.67
JrAHL35JreChr13G1096323,001,456 23,006,865 5410963 3209.7333.66
JrAHL36JreChr14G1030824,429,650 24,430,551 902795 2646.3726.99
JrAHL37JreChr01G129841,416,877 1,417,719 843843 2806.0429.02
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jia, P.; Liu, J.; Yan, R.; Yang, K.; Dong, Q.; Luan, H.; Zhang, X.; Li, H.; Guo, S.; Qi, G. Systematical Characterization of the AT-Hook Gene Family in Juglans regia L. and the Functional Analysis of the JrAHL2 in Flower Induction and Hypocotyl Elongation. Int. J. Mol. Sci. 2023, 24, 7244. https://doi.org/10.3390/ijms24087244

AMA Style

Jia P, Liu J, Yan R, Yang K, Dong Q, Luan H, Zhang X, Li H, Guo S, Qi G. Systematical Characterization of the AT-Hook Gene Family in Juglans regia L. and the Functional Analysis of the JrAHL2 in Flower Induction and Hypocotyl Elongation. International Journal of Molecular Sciences. 2023; 24(8):7244. https://doi.org/10.3390/ijms24087244

Chicago/Turabian Style

Jia, Peng, Jiale Liu, Rui Yan, Kaiyu Yang, Qinglong Dong, Haoan Luan, Xuemei Zhang, Han Li, Suping Guo, and Guohui Qi. 2023. "Systematical Characterization of the AT-Hook Gene Family in Juglans regia L. and the Functional Analysis of the JrAHL2 in Flower Induction and Hypocotyl Elongation" International Journal of Molecular Sciences 24, no. 8: 7244. https://doi.org/10.3390/ijms24087244

APA Style

Jia, P., Liu, J., Yan, R., Yang, K., Dong, Q., Luan, H., Zhang, X., Li, H., Guo, S., & Qi, G. (2023). Systematical Characterization of the AT-Hook Gene Family in Juglans regia L. and the Functional Analysis of the JrAHL2 in Flower Induction and Hypocotyl Elongation. International Journal of Molecular Sciences, 24(8), 7244. https://doi.org/10.3390/ijms24087244

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop