NAC Transcription Factor Research Articles

RD29A and RD29B rearrange genetic and epigenetic markers in priming systemic defense responses against drought and salinity

Drought has become the most important limiting factor to crop productions. Research thus far has been devoted to identifying drought-responsive genes (DRGs) via breeding and engineering approaches. Still, these efforts have not resulted in a solution to combat drought’s effects because the ectopic expression of most DRGs causes adverse effects that reduce plant growth and yields. Lately, we discovered that two DRGs, Response to Desiccation (RD)29A and RD29B, induced by Paenibacillus polymyxa CR1, a plant growth-promoting rhizobacterium capable of priming drought tolerance and concurrently stimulating plant growth, play pivotal roles in defense responses against drought. In this study, we employ the ChlP and qRT-PCR analyses and further clarify that P. polymyxa CR1 reformats the chromatin/transcriptional memory of RD29s, positioned as upstream controllers that fine-tune the temporal dynamic of stress-regulating transcription factors (TFs) in elaborating induced systemic drought tolerance without growth penalties. Two genes coordinate the upregulation of NAC TFs, while feedback inhibiting CBF TFs, which regulate downstream DRG expressions. This supports that RD29s are unique, feasible transgene candidates for improving plants’ survival capacity in both optimal and drought conditions. However, the mode of action of RD29A and RD29B are partly independent, exerting distinct roles in disparate ecological states. When subjected to increasing NaCl concentrations, the KO mutant of RD29A (rd29a) displayed enhanced tolerance compared to WT and rd29b plants, proposing that RD29B, but not RD29A, a key player in conferring WT-like tolerance to salinity stress; further studies will be needed to optimize/maximize their applications in engineering for-profit drought and/or broad-spectrum stress tolerant crops.

Extra- and intranuclear heat perception and triggering mechanisms in plants.

The escalating impact of global warming on crop yield and quality poses a significant threat to future food supplies. Breeding heat-resistant crop varieties holds promise, but necessitates a deeper understanding of the molecular mechanisms underlying plant heat tolerance. Recent studies have shed light on the initial events of heat perception in plants. In this review, we provide a comprehensive summary of the recent progress made in unraveling the mechanisms of heat perception and response in plants. Calcium ion (Ca2+), hydrogen peroxide (H2O2), and nitric oxide (NO) have emerged as key participants in heat perception. Furthermore, we discuss the potential roles of the NAC transcription factor NTL3, thermo-tolerance 3.1 (TT3.1), and Target of temperature 3 (TOT3) as thermosensors associated with the plasma membrane. Additionally, we explore the involvement of cytoplasmic HISTONE DEACETYLASE 9 (HDA9), mRNA encoding the phytochrome-interacting factor 7 (PIF7), and chloroplasts in mediating heat perception. This review also highlights the role of intranuclear transcriptional condensates formed by phytochrome B (phyB), EARLY FLOWERING 3 (ELF3), and guanylate-binding protein (GBP)-like GTPase 3 (GBPL3) in heat perception. Finally, we raise the unresolved questions in the field of heat perception that require further investigation in the future.

Functional characterization of two NAC transcription factors HfNAP1 and HfNAC090 associated with flower programmed cell death in daylily (Hemerocallis fulva)

Daylily (Hemerocallis fulva) is one of the most widely used perennial flowers, but its ornamental and economic value is greatly limited due to its ephemeral flowering period. In general, the flower senescence is regulated by the developmental signals and considered as an irreversible process of programmed cell death (PCD). However, the molecular mechanism of flower PCD in daylily still remains unclear. In this study, two NAC transcription factors, namely HfNAP1 and HfNAC090, are first identified and found to be upregulated significantly in both the age-induced and the ABA-induced flower PCD processes in daylily. Then, the functions of HfNAP1 and HfNAC090 in regulating the flower PCD are investigated through transgenic phenotypes analysis. The results demonstrate that the ectopic and transient overexpression of these two genes can effectively regulate the flower PCD in tobacco and daylily. While the overexpression of HfNAP1 accelerates the flower PCD process, the overexpression of HfNAC090 significantly delays that. Furthermore, the yeast two-hybrid assay is performed to discover potential interactions related to these two genes, and the results demonstrate that HfNAP1 and HfNAC090 can interact with each other, or interact with other flower aging-related genes. Additionally, the yeast one-hybrid assay suggests that HfNAP1 and HfNAC090 can bind directly to the promoters of downstream senescence-associated genes HfSAG39 and HfSAG15. Taken overall, this study provides sufficient evidences to confirm that HfNAP1 and HfNAC090 play dominant roles in regulating the flower PCD in daylily, supporting the development of new strategies to prolong the longevity of daylily flowers.

Genome-wide identification of the NAC gene family and its functional analysis in Liriodendron

As one of the largest plant specific transcription factor families, NAC family members play an important role in plant growth, development and stress resistance. To investigate the function of NAC transcription factors during abiotic stress, as well as during somatic embryogenesis, we identified and characterized the NAC gene family in Liriodendron chinense. We found that most LcNAC members contain more than three exons, with a relatively conserved gene and motif structure, especially at the N-terminus. Interspecies collinearity analysis revealed a closer relationship between the L. chinense NACs and the P. trichocarpa NACs. We analyzed the expression of LcNAC in different tissues and under three abiotic stresses. We found that 12 genes were highly expressed during the ES3 and ES4 stages of somatic embryos, suggesting that they are involved in the development of somatic embryos. 6 LcNAC genes are highly expressed in flower organs. The expression pattern analysis of LcNACs based on transcriptome data and RT-qPCR obtained from L. chinense leaves indicated differential expression responses to drought, cold, and heat stress. Genes in the NAM subfamily expressed differently during abiotic stress, and LcNAC6/18/41/65 might be the key genes in response to abiotic stress. LcNAC6/18/41/65 were cloned and transiently transformed into Liriodendron protoplasts, where LcNAC18/65 was localized in cytoplasm and nucleus, and LcNAC6/41 was localized only in nucleus. Overall, our findings suggest a role of the NAC gene family during environmental stresses in L. chinense. This research provides a basis for further study of NAC genes in Liriodendron chinense.

A pathogen-induced putative NAC transcription factor mediates leaf rust resistance in barley

Leaf rust, caused by Puccinia hordei, is one of the most widespread and damaging foliar diseases affecting barley. The barley leaf rust resistance locus Rph7 has been shown to have unusually high sequence and haplotype divergence. In this study, we isolate the Rph7 gene using a fine mapping and RNA-Seq approach that is confirmed by mutational analysis and transgenic complementation. Rph7 is a pathogen-induced, non-canonical resistance gene encoding a protein that is distinct from other known plant disease resistance proteins in the Triticeae. Structural analysis using an AlphaFold2 protein model suggests that Rph7 encodes a putative NAC transcription factor with a zinc-finger BED domain with structural similarity to the N-terminal DNA-binding domain of the NAC transcription factor (ANAC019) from Arabidopsis. A global gene expression analysis suggests Rph7 mediates the activation and strength of the basal defence response. The isolation of Rph7 highlights the diversification of resistance mechanisms available for engineering disease control in crops.

The activation of the NAC transcription factor by Poplar PtrGATA9 facilitates the secondary cell wall deposition of interfascicular fiber cell in Arabidopsis

The activation of the NAC transcription factor by Poplar PtrGATA9 facilitates the secondary cell wall deposition of interfascicular fiber cell in Arabidopsis

Endosperm cell death promoted by NAC transcription factors facilitates embryo invasion in Arabidopsis

In flowering plants, two fertilization products develop within the limited space of the seed: the embryo and the surrounding nutritive endosperm. The final size of the endosperm is modulated by the degree of embryo growth. In Arabidopsis thaliana, the endosperm expands rapidly after fertilization, but later gets invaded by the embryo that occupies most of the seed volume at maturity, surrounded by a single remaining aleurone-like endosperm layer.1,2,3,4 Embryo invasion is facilitated by the endosperm-expressed bHLH-type transcription factor ZHOUPI, which promotes weakening of endosperm cell walls.5,6 Endosperm elimination in zou mutants is delayed, and embryo growth is severely affected; the endosperm finally collapses around the dwarf embryo, causing the shriveled appearance of mature zou seeds.5,6,7 However, whether ZHOUPI facilitates mechanical endosperm destruction by the invading embryo or whether an active programmed cell death (PCD) process causes endosperm elimination has been subject to debate.2,8 Here we show that developmental PCD controlled by multiple NAC transcription factors in the embryo-adjacent endosperm promotes gradual endosperm elimination. Misexpressing the NAC transcription factor KIRA1 in the entire endosperm caused total endosperm elimination, generating aleurone-less mature seeds. Conversely, dominant and recessive higher-order NAC mutants led to delayed endosperm elimination and impaired cell corpse clearance. Promoting PCD in the zhoupi mutant partially rescued its embryo growth defects, while the endosperm in a zhoupi nac higher-order mutant persisted until seed desiccation. These data suggest that a combination of cell wall weakening and PCD jointly facilitates embryo invasion by an active auto-elimination of endosperm cells.

NACs, generalist in plant life.

Plant-specific NAC proteins constitute a major transcription factor family that is well-known for its roles in plant growth, development, and responses to abiotic and biotic stresses. In recent years, there has been significant progress in understanding the functions of NAC proteins. NAC proteins have a highly conserved DNA-binding domain; however, their functions are diverse. Previous understanding of the structure of NAC transcription factors can be used as the basis for their functional diversity. NAC transcription factors consist of a target-binding domain at the N-terminus and a highly versatile C-terminal domain that interacts with other proteins. A growing body of research on NAC transcription factors helps us comprehend the intricate signalling network and transcriptional reprogramming facilitated by NAC-mediated complexes. However, most studies of NAC proteins have been limited to a single function. Here, we discuss the upstream regulators, regulatory components and targets of NAC in the context of their prospective roles in plant improvement strategies via biotechnology intervention, highlighting the importance of the NAC transcription factor family in plants and the need for further research.

An insertion of transposon in DcNAP inverted its function in the ethylene pathway to delay petal senescence in carnation (Dianthus caryophyllus L.).

Petal senescence is the final stage of flower development. Transcriptional regulation plays key roles in this process. However, whether and how post-transcriptional regulation involved is still largely unknown. Here, we identified an ethylene-induced NAC family transcription factor DcNAP in carnation (Dianthus caryophyllus L.). One allele, DcNAP-dTdic1, has an insertion of a dTdic1 transposon in its second exon. The dTdic1 transposon disrupts the structure of DcNAP and causes alternative splicing, which transcribes multiple domain-deleted variants (DcNAP2 and others). Conversely, the wild type allele DcNAP transcribes DcNAP1 encoding an intact NAC domain. Silencing DcNAP1 delays and overexpressing DcNAP1 accelerates petal senescence in carnation, while silencing and overexpressing DcNAP2 have the opposite effects, respectively. Further, DcNAP2 could interact with DcNAP1 and interfere the binding and activation activity of DcNAP1 to the promoters of its downstream target ethylene biosynthesis genes DcACS1 and DcACO1. Lastly, ethylene signalling core transcriptional factor DcEIL3-1 can activate the expression of DcNAP1 and DcNAP2 in the same way by binding their promoters. In summary, we discovered a novel mechanism by which DcNAP regulates carnation petal senescence at the post-transcriptional level. It may also provide a useful strategy to manipulate the NAC domains of NAC transcription factors for crop genetic improvement.

Integrated metabolic profiling and transcriptome analysis of Lonicera japonica flowers for chlorogenic acid, luteolin and endogenous hormone syntheses

The active ingredients of many medicinal plants are the secondary metabolites associated with the growth period. Lonicera japonica Thunb. is an important traditional Chinese medicine, and the flower development stage is an important factor that influences the quality of medicinal ingredients. In this study, transcriptomics and metabolomics were performed to reveal the regulatory mechanism of secondary metabolites during flowering of L. japonica. The results showed that the content of chlorogenic acid (CGA) and luteolin gradually decreased from green bud stage (Sa) to white flower stage (Sc), especially from white flower bud stage (Sb) to Sc. Most of the genes encoding the crucial rate-limiting enzymes, including PAL, C4H, HCT, C3′H, F3′H and FNSII, were down-regulated in three comparisons. Correlation analysis identified some members of the MYB, AP2/ERF, bHLH and NAC transcription factor families that are closely related to CGA and luteolin biosynthesis. Furthermore, differentially expressed genes (DEGs) involved in hormone biosynthesis, signalling pathways and flowering process were analysed in three flower developmental stage.

The Interplay between Enucleated Sieve Elements and Companion Cells.

In order to adapt to sessile life and terrestrial environments, vascular plants have developed highly sophisticated cells to transport photosynthetic products and developmental signals. Of these, two distinct cell types (i.e., the sieve element (SE) and companion cell) are arranged in precise positions, thus ensuring effective transport. During SE differentiation, most of the cellular components are heavily modified or even eliminated. This peculiar differentiation implies the selective disintegration of the nucleus (i.e., enucleation) and the loss of cellular translational capacity. However, some cellular components necessary for transport (e.g., plasmalemma) are retained and specific phloem proteins (P-proteins) appear. Likewise, MYB (i.e., APL) and NAC (i.e., NAC45 and NAC86) transcription factors (TFs) and OCTOPUS proteins play a notable role in SE differentiation. The maturing SEs become heavily dependent on neighboring non-conducting companion cells, to which they are connected by plasmodesmata through which only 20-70 kDa compounds seem to be able to pass. The study of sieve tube proteins still has many gaps. However, the development of a protocol to isolate proteins that are free from any contaminating proteins has constituted an important advance. This review considers the very detailed current state of knowledge of both bound and soluble sap proteins, as well as the role played by the companion cells in their presence. Phloem proteins travel long distances by combining two modes: non-selective transport via bulk flow and selective regulated movement. One of the goals of this study is to discover how the protein content of the sieve tube is controlled. The majority of questions and approaches about the heterogeneity of phloem sap will be clarified once the morphology and physiology of the plasmodesmata have been investigated in depth. Finally, the retention of specific proteins inside an SE is an aspect that should not be forgotten.

Genome-wide identification and molecular evolution of NAC gene family in Dendrobium nobile.

NAC transcription factors are an important genes that regulate plant growth and development, and can regulate functions such as fruit ripening in plants. Based on genome data of Dendrobium nobile, the NAC gene family was identified and analyzed by bioinformatics methods. In this study, we identified 85 NAC genes in Dendrobium nobile genome, and systematically analyzed the NAC gene family. We found that they were distributed unevenly in the nineteen chromosomes. The amino acid length of D. nobile NAC gene family (DnoNACs) ranged from 80 to 1065, molecular weight ranged from 22.17 to 119.02 kD, and isoelectric point ranged from 4.61~9.26. Its promoter region contains multiple stress responsive elements, including light responsive, gibberellin-responsive, abscisic acid responsiveness, MeJA-responsiveness and drought-inducibility elements. Phylogenetic analysis indicates that the D. nobile NAC gene family is most closely related to Dendrobium catenatum and Dendrobium chrysotoxum. Analysis of SSR loci indicates that the fraction of mononucleotide repeats was the largest, as was the frequency of A/T. Non-coding RNA analysis showed that these 85 NAC genes contain 397 miRNAs. The collinearity analysis shows that 9 collinear locis were found on the chromosomes of D. nobile with Arabidopsis thaliana, and 75 collinear locis with D.chrysotoxum. QRT-PCR experiment under different salt concentration and temperature conditions verified the response mechanism of DnoNAC gene family under stress conditions. Most DnoNAC genes are sensitive to salt stress and temperature stress. The results of this study provide a reference for further understanding the function of NAC gene in D. nobile.

Molecular regulation of apple and grape ripening: exploring common and distinct transcriptional aspects of representative climacteric and non-climacteric fruits.

Fleshy fruits of angiosperms are organs specialized for promoting seed dispersal by attracting herbivores and enticing them to consume the organ and the contained seeds. Ripening can be broadly defined as the processes serving as a plant strategy to make the fleshy fruit appealing to animals and leads to a coordinated series of changes in color, texture, aroma and flavor, as result of an intricate interplay of genetically and epigenetically programmed events. The ripening of fruits can be categorized into two types: climacteric, which is characterized by a rapid increase in respiration rate typically accompanied by a burst of ethylene production, and non-climacteric, where this pronounced peak in respiration is absent. Here we review the most current knowledge on transcriptomic changes taking place in apple (climacteric) and grapevine (non-climacteric) fruit during ripening, with the aim to highlight specific and common hormonal and molecular events governing the process in both species. In this perspective, we found that specific NAC transcription factor members participate in the ripening initiation in grape and are involved in the attempt to restore the normal physiological ripening progression in impaired fruit ripening physiology in apple. These elements suggest the existence of a common regulatory mechanism operated by NAC transcription factors and auxin in the two species.

Overexpression of BnNAC19 in Brassica napus enhances resistance to Leptosphaeria maculans, the blackleg pathogen of canola

AbstractLeptosphaeria maculans is a fungal pathogen that causes blackleg disease in canola (Brassica napus), resulting in significant yield and economic losses in Canada and many parts of the world. Plant NAC transcription factors play critical roles in plant development and response to biotic or abiotic stress. In this study, we identified and characterized a BnNAC19 gene from Brassica napus. The overexpression of BnNAC19 in transgenic canola plants contributed to the improvement of seedling resistance against L. maculans. The mycelial growth of a green fluorescent protein‐tagged strain of L. maculans and production of pycnidiospores were shown to be inhibited in the transgenic canola plants overexpressing BnNAC19. In addition, the canola transgenic line overexpressing BnNAC19 showed increased disease resistance in the adult plant, which was determined by quantitative resistance. Both increased seedling and adult plant resistance in transgenic canola plants overexpressing BnNAC19 indicate that the BnNAC19 gene plays a positive role against L. maculans. The expression pattern of genes BnNAC19 upstream and downstream of BnNAC19 that participate in plant defence pathways were investigated to elucidate the B. napus resistance mechanisms to L. maculans infection, and hence to aid the long‐term blackleg disease‐resistant breeding programmes.

Transcriptome analysis during vernalization in wheat (Triticum aestivum L.)

BackgroundVernalization, as a vital process in the life cycle of winter cereal, has important effects on floral organ formation and flowering time. Many morphological changes together with molecular changes occur during the vernalization period. Here, we used transcriptome sequencing to analyze the transcriptomic changes in wheat leaves before, during and after vernalization using the winter wheat cultivar ‘Shiluan02-1’.ResultsA total of 16,370 differentially expressed genes were obtained across different vernalization periods. Gene Ontology enrichment analysis revealed that photoperiodism, photoprotection, photosynthesis, lipid transport and biosynthetic process, and chlorophyll metabolic process were closely related to vernalization. In addition, AP2/ERF, C2H2, bHLH, WRKY, MYB, MYB-related, and NAC transcription factors were significantly enriched during vernalization, and the transcription factor expression patterns suggested the intricate regulation of transcription factor modules in plant vernalization pathways. Analysis of gene expression patterns of the MADS-box transcription factor genes showed different expression patterns during vernalization phases, among which VERNALIZATION1 (VRN1) genes were found to gradually increase during vernalization periods from V0 to V35, while decline in the V42 phase, then increase after vernalization. The Tavrt-2 gene cooperated with Tavrn1 to regulate flowering induced by vernalization, and its expression level was rapidly increased by vernalization but declined in the V42 phase and then increased after vernalization. Some genes from the ICE-CBF-COR pathway were also identified, and additional analysis indicated that some key genes related to phytohormone biosynthesis and signal transduction were enriched during the vernalization period, such as gibberellic acid, ethylene, abscisic acid and jasmonic acid biosynthesis and signaling pathway genes.ConclusionsOur study provides valuable molecular information for future studies on wheat vernalization regulation and also serves as an excellent reference for future wheat breeding.

The NAC transcription factors SNAP1/2/3/4 are central regulators mediating high nitrogen responses in mature nodules of soybean

Legumes can utilize atmospheric nitrogen via symbiotic nitrogen fixation, but this process is inhibited by high soil inorganic nitrogen. So far, how high nitrogen inhibits N2 fixation in mature nodules is still poorly understood. Here we construct a co-expression network in soybean nodule and find that a dynamic and reversible transcriptional network underlies the high N inhibition of N2 fixation. Intriguingly, several NAC transcription factors (TFs), designated as Soybean Nitrogen Associated NAPs (SNAPs), are amongst the most connected hub TFs. The nodules of snap1/2/3/4 quadruple mutants show less sensitivity to the high nitrogen inhibition of nitrogenase activity and acceleration of senescence. Integrative analysis shows that these SNAP TFs largely influence the high nitrogen transcriptional response through direct regulation of a subnetwork of senescence-associated genes and transcriptional regulators. We propose that the SNAP-mediated transcriptional network may trigger nodule senescence in response to high nitrogen.

Deciphering the regulatory network of the NAC transcription factor FvRIF, a key regulator of strawberry (Fragaria vesca) fruit ripening.

The NAC transcription factor ripening inducing factor (RIF) was previously reported to be necessary for the ripening of octoploid strawberry (Fragaria × ananassa) fruit, but the mechanistic basis of RIF-mediated transcriptional regulation and how RIF activity is modulated remains elusive. Here, we show that FvRIF in diploid strawberry, Fragaria vesca, is a key regulator in the control of fruit ripening and that knockout mutations of FvRIF result in a complete block of fruit ripening. DNA affinity purification sequencing coupled with transcriptome deep sequencing suggests that 2,080 genes are direct targets of FvRIF-mediated regulation, including those related to various aspects of fruit ripening. We provide evidence that FvRIF modulates anthocyanin biosynthesis and fruit softening by directly regulating the related core genes. Moreover, we demonstrate that FvRIF interacts with and serves as a substrate of MAP kinase 6 (FvMAPK6), which regulates the transcriptional activation function of FvRIF by phosphorylating FvRIF at Thr-310. Our findings uncover the FvRIF-mediated transcriptional regulatory network in controlling strawberry fruit ripening and highlight the physiological significance of phosphorylation modification on FvRIF activity in ripening.

Exogenous IAA application affects the specific characteristics of fluoranthene distribution in Arabidopsis

Indole-3-acetic acid (IAA) is a crucial growth regulator involved in the accumulation of polycyclic aromatic hydrocarbons (PAHs). However, the precise physiological and molecular mechanisms underlying IAA-mediated plant growth and PAH accumulation are not yet fully understood. In this study, two distinct IAA-sensitive genotypes of Arabidopsis thaliana (wild type and Axr5 mutant) were chosen to investigate the mechanisms of fluoranthene (Flu) uptake and accumulation in plant tissues (roots and leaves) through physiological and molecular analyses. The results revealed that the Flu concentration in Axr5 leaves was significantly higher than that in wild-type (WT) leaves. In roots, the Flu content decreased significantly with increasing IAA treatment, while no significant changes were observed with lower IAA treatment. Principal component analysis demonstrated that Flu accumulation in Arabidopsis roots was associated with IAA concentrations, whereas Flu accumulation in leaves was dependent on the genotype. Moreover, Flu accumulation showed a positive correlation with the activity of glutathione S-transferase (GST) and root length and a positive correlation with catalase (CAT) and peroxidase (POD) activity in the leaves. Transcriptome analysis confirmed that the expression of the ethylene-related gene ATERF6 and GST-related genes ATGSTF14 and ATGSTU27 in roots, as well as the POD-related genes AtPRX9 and AtPRX25 and CAT-related gene AtCAT3 in leaves, played a role in Flu accumulation. Furthermore, WRKY transcription factors (TFs) in roots and NAC TFs in leaves were identified as important regulators of Flu accumulation. Understanding the mechanisms of Flu uptake and accumulation in A. thaliana provides valuable insights for regulating PAH accumulation in plants.

The NAC Transcription Factor PgNAC41-2 Gene Involved in the Regulation of Ginsenoside Biosynthesis in Panax ginseng.

Ginseng (Panax ginseng C.A. Meyer) is a perennial herb of the Araliaceae family, a traditional and valuable Chinese herb in China. The main active component of ginseng is ginsenoside. The NAC transcription factors belong to a large family of plant-specific transcription factors, which are involved in growth and development, stress response and secondary metabolism. In this study, we mapped the NAC gene family on 24 pairs of ginseng chromosomes and found numerous gene replications in the genome. The NAC gene PgNAC41-2, found to be highly related to ginsenoside synthesis, was specifically screened. The phylogeny and expression pattern of the PgNAC41-2 gene were analyzed, along with the derived protein sequence, and a structure model was generated. Furthermore, the PgNAC41-2 gene was cloned and overexpressed by a Rhizobium rhizogenes mediated method, using ginseng petioles as receptor material. The saponin content of the transformed material was analyzed to verify the function of the NAC transcription factor in ginseng. Our results indicate that the PgNAC41-2 gene positively regulates the biosynthesis of saponins.

Genome-Wide Analysis of WRKY and NAC Transcription Factors in Carica papaya L. and Their Possible Role in the Loss of Drought Tolerance by Recent Cultivars through the Domestication of Their Wild Ancestors.

A genome-wide analysis for two families of key transcription factors (TF; WRKY and NAC) involved in drought response revealed 46 WRKY and 66 NAC members of the Carica papaya genome. A phylogenetic analysis grouped the CpWRKY proteins into three groups (I, II a, b, c, d, e and III), while the CpNAC proteins were clustered into 15 groups. The conserved domains, chromosomal localization and promoter cis-acting elements were also analyzed. In addition, from a previous transcriptome study of two contrasting genotypes in response to 14 days of water deficit stress (WDS), we found that 29 of the 46 CpWRKYs genes and 25 of the 66 CpNACs genes were differentially expressed in response to the WDS. In the present paper, the native wild genotype (WG) (collected in its center of origin) consistently showed a higher expression (transcripts per million; TPM and fold change; FC) than the commercial genotype (CG) in almost all the members of the CpWRKY and CpNAC gene families. To corroborate this, we selected CpWRKY50 and CpNAC83.1 for further evaluation by RT-qPCR. Consistently, the WG showed higher relative expression levels (REL) after 14 days of WDS than the CG, in both the leaves and roots. The results suggest that the CpWRKY and CpNAC TF families are important for drought tolerance in this species. The results may also suggest that, during the domestication process, the ability of the native (wild) C. papaya genotypes to respond to drought (including the overexpression of the CpWRKY and CpNAC genes) was somehow reduced in the current commercial genotypes.