Genome-wide analyses and expression patterns under abiotic stress of NAC transcription factors in white pear (Pyrus bretschneideri)

CcNAC6 Acts as a Positive Regulator of Secondary Cell Wall Synthesis in Sudan Grass (Sorghum sudanense S.)

The degree of forage lignification is a key factor affecting its digestibility by ruminants such as cattle and sheep. Sudan grass (Sorghum sudanense S.) is a high-quality sorghum forage, and its lignocellulose is mostly stored in the secondary cell wall. However, the secondary cell wall synthesis mechanism of Sudan grass has not yet been studied in depth. To further study the secondary cell wall synthesis mechanism of Sudan grass using established transcriptome data, this study found that CcNAC6, a homologous gene of Arabidopsis AtSND2, is related to the secondary cell wall synthesis of Sudan grass. Accordingly, we constructed a CcNAC6-overexpressing line of Arabidopsis to investigate the function of the CcNAC6 gene in secondary cell wall synthesis. The results showed that the overexpression of the CcNAC6 gene could significantly increase the lignin content of Arabidopsis. Based on subcellular localization analysis, CcNAC6 is found in the nucleus. In addition, yeast two-hybridization screening showed that CcCP1, associated with secondary cell wall synthesis, can interact with CcNAC6. Therefore, the above results indicate that CcNAC6 has a positive regulatory effect on the secondary cell wall synthesis of Sudan grass, and it is speculated that CcNAC6 may be the main regulator of the secondary cell wall synthesis of Sudan grass through its interaction with another regulatory protein, CcCP1. This study provides a theoretical basis and new genetic resources for the creation of new Sudan grass germplasm with a low lignin content.

Transcription factor PbMYB80 regulates lignification of stone cells and undergoes RING finger protein PbRHY1-mediated degradation in pear fruit

The Chinese white pear (Pyrus bretschneideri) fruit carries a high proportion of stone cells, adversely affecting fruit quality. Lignin is a main component of stone cells in pear fruit. In this study, we discovered that a pear MYB transcription factor, PbMYB80, binds to the promoters of key lignin biosynthesis genes and inhibits their expression. Stable overexpression of PbMYB80 in Arabidopsis showed that lignin deposition and secondary wall thickening were inhibited, and the expression of the lignin biosynthesis genes in transgenic Arabidopsis was decreased. Transient overexpression of PbMYB80 in pear fruit inhibited lignin metabolism and stone cell development, and the expression of some genes in the lignin metabolism pathway was reduced. In contrast, silencing PbMYB80 with VIGS increased the lignin and stone cell content in pear fruit, and increased expression of genes in the lignin metabolism pathway. By screening a pear fruit cDNA library in yeast, we found that PbMYB80 binds to a RING finger (PbRHY1) protein. We also showed that PbRHY1 exhibits E3 ubiquitin ligase activity and degrades ubiquitinated PbMYB80 in vivo and in vitro. This investigation contributes to a better understanding of the regulation of lignin biosynthesis in pear fruit, and provides a theoretical foundation for increasing pear fruit quality at the molecular level.

Comprehensive Analysis of NAC Transcription Factors Reveals Their Evolution in Malvales and Functional Characterization of AsNAC019 and AsNAC098 in Aquilaria sinensis

NAC is a class of plant-specific transcription factors that are widely involved in the growth, development and (a)biotic stress response of plants. However, their molecular evolution has not been extensively studied in Malvales, especially in Aquilaria sinensis, a commercial and horticultural crop that produces an aromatic resin named agarwood. In this study, 1502 members of the NAC gene family were identified from the genomes of nine species from Malvales and three model plants. The macroevolutionary analysis revealed that whole genome duplication (WGD) and dispersed duplication (DSD) have shaped the current architectural structure of NAC gene families in Malvales plants. Then, 111 NAC genes were systemically characterized in A. sinensis. The phylogenetic analysis suggests that NAC genes in A. sinensis can be classified into 16 known clusters and four new subfamilies, with each subfamily presenting similar gene structures and conserved motifs. RNA-seq analysis showed that AsNACs presents a broad transcriptional response to the agarwood inducer. The expression patterns of 15 AsNACs in A. sinensis after injury treatment indicated that AsNAC019 and AsNAC098 were positively correlated with the expression patterns of four polyketide synthase (PKS) genes. Additionally, AsNAC019 and AsNAC098 were also found to bind with the AsPKS07 promoter and activate its transcription. This comprehensive analysis provides valuable insights into the molecular evolution of the NAC gene family in Malvales plants and highlights the potential mechanisms of AsNACs for regulating secondary metabolite biosynthesis in A. sinensis, especially for the biosynthesis of 2-(2-phenyl) chromones in agarwood.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Gene Family in Litchi chinensis

NAC proteins play an essential role in the growth and development of litchi, especially during reproductive development. However, a comprehensive analysis of the litchi NAC gene family is currently absent. Based on information from the litchi genome, we found that the 112 NAC genes of litchi show an uneven distribution on the chromosomes. Phylogenetic and conserved structural domain analyses indicated that different types of variability were exhibited in the family of litchi NACs (LcNACs). Gene covariance analysis showed that the LcNACs showed better similarity in the same genus than with Arabidopsis. We further investigated the differential expression patterns of LcNACs in buds and rudimentary leaves of litchi. qRT-PCR results implied that they were involved in the process. Profiling of LcNAC promoter elements in litchi showed that they were extensively involved in light response, phytohormone regulation, abiotic stress response, and plant growth and development processes. This study provides new insights into the identification, structural characterization, tissue-specific expression analysis, and promoter response elements of LcNACs. It reveals the characteristics of the LcNACs and lays the foundation for the subsequent understanding of its biological functions and molecular regulatory mechanisms.

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Genome-wide identification of the NAC transcription factor family and regulation of metabolites under salt stress in Isatis indigotica

NAC (NAM, ATAF1/2 and CUC2) transcription factors (TFs) are a class of TFs families unique to plants, which not only play an important role in the growth and developmental stages of plants but also function in response to stress and regulation of secondary metabolite biosynthesis. However, there are few studies on NAC genes in the medicinal plant Isatis indigotica. In this study, 96 IiNAC genes were identified based on the whole-genome data of I. indigotica, distributed in seven chromosomes and three contigs. IiNAC genes were structurally conserved and divided into 15 subgroups. Cis-elements were identified in the promoter region of the IiNAC gene in response to plant growth and development, abiotic stresses and hormones. In addition, transcriptome and metabolome data of I. indigotica leaves under salt stress were analyzed to construct a network of IiNAC gene co-expression and metabolite association. Ten differentially expressed IiNAC genes were co-expressed with 109 TFs, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that most of these genes were associated with plant growth and development and abiotic stress responses. Eleven IiNAC genes were positively associated with 72 metabolites. Eleven IiNAC genes were positively or negatively associated with 47 metabolites through 37 TFs. Commonly associated secondary metabolites include two terpenoids, abscisic acid and bilobalide, two flavonoids, dihydrokaempferol and syringaldehyde, a coumarin, 7-methoxycoumarin, an alkaloid, lupinine, and quinone dihydrotanshinone I. This study provides important data to support the identification of the NAC gene family in I. indigotica and the regulatory functions of IiNAC genes in metabolites under salt stress.

The NAC transcription factor family in Eucommia ulmoides: Genome-wide identification, characterization, and network analysis in relation to the rubber biosynthetic genes

The NAC transcription factor family is a large plant gene family, participating in plant growth and development, secondary metabolite synthesis, biotic and abiotic stresses responses, and hormone signaling. Eucommia ulmoides is a widely planted economic tree species in China that can produce trans-polyisoprene: Eucommia rubber (Eu-rubber). However, genome-wide identification of the NAC gene family has not been reported in E. ulmoides. In this study, 71 NAC proteins were identified based on genomic database of E. ulmoides. Phylogenetic analysis showed that the EuNAC proteins were distributed in 17 subgroups based on homology with NAC proteins in Arabidopsis, including the E. ulmoides-specific subgroup Eu_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 7, and multitudinous EuNAC genes contained two or three exons. Chromosomal location analysis revealed that the EuNAC genes were unevenly distributed on 16 chromosomes. Three pairs of genes of tandem duplicates genes and 12 segmental duplications were detected, which indicated that segmental duplications may provide the primary driving force of expansion of EuNAC. Prediction of cis-regulatory elements indicated that the EuNAC genes were involved in development, light response, stress response and hormone response. For the gene expression analysis, the expression levels of EuNAC genes in various tissues were quite different. To explore the effect of EuNAC genes on Eu-rubber biosynthesis, a co-expression regulatory network between Eu-rubber biosynthesis genes and EuNAC genes was constructed, which indicated that six EuNAC genes may play an important role in the regulation of Eu-rubber biosynthesis. In addition, this six EuNAC genes expression profiles in E. ulmoides different tissues were consistent with the trend in Eu-rubber content. Quantitative real-time PCR analysis showed that EuNAC genes were responsive to different hormone treatment. These results will provide a useful reference for further studies addressing the functional characteristics of the NAC genes and its potential role in Eu-rubber biosynthesis.

Comprehensive analysis of the NAC transcription factor gene family in Kandelia obovata reveals potential members related to chilling tolerance

BACKGROUND

Kandelia obovata is an important mangrove species extensively distributed in Eastern Asia that is susceptible to low-temperature stress. NAC (NAM, ATAF1/2 and CUC2) domain proteins are transcription factors (TFs) that play various roles in plant growth and development and in the plant response to environmental stresses. Nevertheless, genome-wide analyses of K. obovata NAC genes (KoNACs) and their responses to chilling stress have rarely been studied.

METHODS

The KoNAC gene family was identified and characterized using bioinformatic analysis, the subcellular location of some NAC proteins was confirmed using confocal microscopy analysis, and the KoNACs that responded to chilling stress were screened using RNA-seq and qRT-PCR analysis.

RESULTS

A total of 79 KoNACs were identified, and they were unequally distributed across all 18 chromosomes of K. obovata. The KoNAC proteins could be divided into 16 subgroups according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana. The KoNACs exhibited greater synteny with A. thaliana sequences than with Oryza sativa sequences, indicating that KoNACs underwent extensive evolution after the divergence of dicotyledons and monocotyledons. Segmental duplication was the main driving force of the expansions of KoNAC genes. Confocal microscopy analysis verified that the four randomly selected KoNACs localized to the nucleus, indicating the accuracy of the bioinformatic predictions. Tissue expression pattern analysis demonstrated that some KoNAC genes showed tissue-specific expression, suggesting that these KoNACs might be important for plant development and growth. Additionally, the expression levels of 19 KoNACs were significantly (15 positively and 4 negatively) induced by cold treatment, demonstrating that these KoNACs might play important roles during cold stress responses and might be candidate genes for the genetic engineering of K. obovata with enhanced chilling stress tolerance. Coexpression network analysis revealed that 381 coexpressed pairs (between 13 KoNACs and 284 other genes) were significantly correlated.

CONCLUSIONS

Seventy-nine KoNACs were identified in K. obovata, nineteen of which displayed chilling-induced expression patterns. These genes may serve as candidates for functional analyses of KoNACs engaged in chilling stress. Our results lay the foundation for evolutionary analyses of KoNACs and their molecular mechanisms in response to environmental stress.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

Genome-Wide Analysis of the NAC Family Associated with Two Paleohexaploidization Events in the Tomato

NAC transcription factors play an important regulatory role in tomato fruit ripening. We chose a novel perspective to explore the traces left by two paleopolyploidizations in the NAC family using a bioinformatics approach. We found that 85 (S. lycopersicum) and 88 (S. pennellii) members of the NAC family were present in two tomatoes, and most of them were amplified from two paleohexaploidizations. We differentiated NAC family members from the different paleohexaploidizations and found that the SWGT-derived NAC genes had more rearrangement events, so it was different from the DWGT-derived NAC genes in terms of physicochemical properties, phylogeny, and gene location. The results of selection pressure show that DWGT-derived NAC genes tended to be positively selected in S. lycopersicum and negatively selected in S. pennellii. A comprehensive analysis of paleopolyploidization and expression reveals that DWGT-derived NAC genes tend to promote fruit ripening, and are expressed at the early and middle stages, whereas SWGT-derived NAC genes tend to terminate fruit growth and are expressed at the late stages of fruit ripening. This study obtained NAC genes from different sources that can be used as materials for tomato fruit development, and the method in the study can be extended to the study of other plants.

Genome-Wide Identification and Analysis of the NAC Transcription Factor Gene Family in Garden Asparagus (Asparagus officinalis)

As a large plant-specific gene family, the NAC (NAM, ATAF1/2, and CUC2) transcription factor is related to plant growth, development, and response to abiotic stresses. Although the draft genome of garden asparagus (Asparagus officinalis) has been released, the genome-wide investigation of the NAC gene family is still unavailable. In this study, a total of 85 A. officinalis NAC genes were identified, and a comprehensive analysis of the gene family was performed, including physicochemical properties, phylogenetic relationship, chromosome localization, gene structure, conserved motifs, intron/exon, cis-acting elements, gene duplication, syntenic analysis, and differential gene expression analysis. The phylogenetic analysis demonstrated that there were 14 subgroups in both A. officinalis and Arabidopsis thaliana, and the genes with a similar gene structure and motif distribution were clustered in the same group. The cis-acting regulatory analysis of AoNAC genes indicated four types of cis-acting elements were present in the promoter regions, including light-responsive, hormone-responsive, plant-growth-and-development-related, and stress-responsive elements. The chromosomal localization analysis found that 81 NAC genes in A. officinalis were unevenly distributed on nine chromosomes, and the gene duplication analysis showed three pairs of tandem duplicated genes and five pairs of segmental duplications, suggesting that gene duplication is possibly associated with the amplification of the A. officinalis NAC gene family. The differential gene expression analysis revealed one and three AoNAC genes that were upregulated and downregulated under different types of salinity stress, respectively. This study provides insight into the evolution, diversity, and characterization of NAC genes in garden asparagus and will be helpful for future understanding of their biological roles and molecular mechanisms in plants.

NtNAC053, A Novel NAC Transcription Factor, Confers Drought and Salt Tolerances in Tobacco

The NAC (NAM, ATAF1/2, and CUC2) family acts as one of the largest families of the transcription factor in the plant kingdom and was revealed to function as the important regulators in various environmental stresses. However, a few studies were reported about the biofunctions of the NAC transcription factor in tobacco. In the current study, we characterized a novel NAC transcription factor encoding the gene NtNAC053 in tobacco, which was significantly up-regulated when exposed to salt and drought treatments. The results of cis-acting elements analysis suggested that the promoter region of NtNAC053 gene possesses a number of stress-responsive elements, and this gene could be induced by exogenous abscisic acid (ABA) treatment. Moreover, the NtNAC053-GFP fusion protein was localized in the cell nucleus and possessed a transactivation domain in its C-terminal, implying that NtNAC053 may undertake as a transcriptional activator in tobacco. Notably, the overexpression of NtNAC053 in tobacco resulted in hypersensitivity to ABA treatment. Furthermore, these overexpression lines showed significantly enhanced tolerances to drought and salt stresses. Under salt and drought stresses, these overexpression lines possessed higher superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) activities. Interestingly, the expressions of putative stress-related genes, including NtCOR15A, NtRAB18, NtDREB1A, NtERF5, NtKAT2, and NtERD11, were up-regulated in these overexpression lines when subjected to salt and drought stresses. The clues provided in our study suggested that the NtNAC053 gene encodes a novel NAC transcription factor and could confer the drought and salt stress tolerances by inspiring the downstream stress-responsive genes and antioxidant system in tobacco.

Influence of Climate Change on Metabolism and Biological Characteristics in Perennial Woody Fruit Crops in the Mediterranean Environment

The changes in the state of the climate have a high impact on perennial fruit crops thus threatening food availability. Indeed, climatic factors affect several plant aspects, such as phenological stages, physiological processes, disease-pest frequency, yield, and qualitative composition of the plant tissues and derived products. To mitigate the effects of climatic parameters variability, plants implement several strategies of defense, by changing phenological trends, altering physiology, increasing carbon sequestration, and metabolites synthesis. This review was divided into two sections. The first provides data on climate change in the last years and a general consideration on their impact, mitigation, and resilience in the production of food crops. The second section reviews the consequences of climate change on the industry of two woody fruit crops models (evergreen and deciduous trees). The research focused on, citrus, olive, and loquat as evergreen trees examples; while grape, apple, pear, cherry, apricot, almond, peach, kiwi, fig, and persimmon as deciduous species. Perennial fruit crops originated by a complex of decisions valuable in a long period and involving economic and technical problems that farmers may quickly change in the case of annual crops. However, the low flexibility of woody crops is balanced by resilience in the long-life cycle.

HuNAC20 and HuNAC25, Two Novel NAC Genes from Pitaya, Confer Cold Tolerance in Transgenic Arabidopsis

NAC transcription factors are one of the largest families of transcriptional regulators in plants, and members of the gene family play vital roles in regulating plant growth and development processes including biotic/abiotic stress responses. However, little information is available about the NAC family in pitaya. In this study, we conducted a genome-wide analysis and a total of 64 NACs (named HuNAC1-HuNAC64) were identified in pitaya (Hylocereus). These genes were grouped into fifteen subgroups with diversities in gene proportions, exon–intron structures, and conserved motifs. Genome mapping analysis revealed that HuNAC genes were unevenly scattered on all eleven chromosomes. Synteny analysis indicated that the segmental duplication events played key roles in the expansion of the pitaya NAC gene family. Expression levels of these HuNAC genes were analyzed under cold treatments using qRT-PCR. Four HuNAC genes, i.e., HuNAC7, HuNAC20, HuNAC25, and HuNAC30, were highly induced by cold stress. HuNAC7, HuNAC20, HuNAC25, and HuNAC30 were localized exclusively in the nucleus. HuNAC20, HuNAC25, and HuNAC30 were transcriptional activators while HuNAC7 was a transcriptional repressor. Overexpression of HuNAC20 and HuNAC25 in Arabidopsis thaliana significantly enhanced tolerance to cold stress through decreasing ion leakage, malondialdehyde (MDA), and H₂O₂ and O₂⁻ accumulation, accompanied by upregulating the expression of cold-responsive genes (AtRD29A, AtCOR15A, AtCOR47, and AtKIN1). This study presents comprehensive information on the understanding of the NAC gene family and provides candidate genes to breed new pitaya cultivars with tolerance to cold conditions through genetic transformation.

Identification of bZIP transcription factors and their responses to brown spot in pear

Basic leucine zipper (bZIP) is a conserved transcription factor (TF) widely present in eukaryotes, and it plays an important role in regulating plant growth and stress responses. To better understand the white pear bZIP gene family, comprehensive bioinformatics analysis of the pear genome was performed. A total of 84 PbbZIP genes were identified, which were divided into 13 subfamilies by phylogenetic analysis. The 84 PbbZIP genes were all located in the nucleus, and 77 of those genes were unevenly distributed across the 17 chromosomes of white pear. The other 7 PbbZIP genes were located on the scaffold. Subsequent expression profile analysis showed that PbbZIP genes in exocarp were significantly upregulated or downregulated in ‘Huangguan’ pear with brown spot (BS) compared with healthy pear and in response to hormonal treatment with gibberellin A₃ (GA₃). These results provide helpful insights into the characteristics of PbbZIP genes and their responses to BS in ‘Huangguan’ pear.

Peduncle Necking in Rosa hybrida Induces Stress-Related Transcription Factors, Upregulates Galactose Metabolism, and Downregulates Phenylpropanoid Biosynthesis Genes

Roses are highly valued as cut flowers worldwide but have limited vase life. Peduncle bending "bent neck" or "necking" is a major cause of reduced vase life, especially in some cultivars. Necking is thought to be caused by either an air embolism or accumulation of microorganisms at or within the stem end, blocking the xylem vessels and preventing water uptake. However, the underlying mechanisms of necking are poorly understood. Here, RNAseq analysis was applied to compare gene expression across three stages of peduncle necking (straight, <90°, and >90°), in the necking-susceptible Rosa hybrida cultivar H30. Most gene expression change was later in bending and there was, overall, more downregulation than upregulation of gene expression during necking. Photosynthetic, starch, and lignin biosynthesis genes were all downregulated, while genes associated with galactose metabolism, producing raffinose and trehalose that are both related to osmoprotection, were upregulated. Genes associated with starch breakdown, autophagy, and senescence were also upregulated, as were most of the NAC and WRKY transcription factors, involved in stress and senescence regulation. Microscopy showed a cellular collapse in the peduncle. These data support a possible mechanism, whereby a reduction in water transport leads to a cellular collapse in the peduncle, accompanied by upregulation of senescence and drought responses.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Genome-wide identification and expression analysis of the bZIP transcription factors, and functional analysis in response to drought and cold stresses in pear (Pyrus breschneideri)

BACKGROUND

Transcription factors (TFs) are involved in many important biological processes, including cell stretching, histological differentiation, metabolic activity, seed storage, gene regulation, and response to abiotic and biotic stresses. Little is known about the functions, evolutionary history, and expression patterns of basic region-leucine zipper TF family genes in pear, despite the release of the genome of Chinese white pears ("Dangshansuli").

RESULTS

Overall, 92 bZIP genes were identified in the pear genome (Pyrus breschneideri). Of these, 83 were randomly distributed on all 17 chromosomes except chromosome 4, and the other 9 genes were located on loose scaffolding. The genes were divided into 14 subgroups. Whole-genome duplications, dispersed duplication, and purifying selection for whole-genome duplications are the main reasons for the expansion of the PbrbZIP gene family. The analysis of functional annotation enrichment indicated that most of the functions of PbrbZIP genes were enriched in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways involved in the abiotic stress response. Next, expression analysis and virus-induced gene silencing results indicated that PbrbZIP genes might play critical roles in response to drought and cold stresses, especially for the genes from subgroups A, C, G, I, and S.

CONCLUSIONS

Ninety-two PbrbZIP genes were identified from the pear genome and classified into 14 subgroups. PbrbZIP genes were mainly expanded from whole-genome duplications and dispersed duplications and retained by purifying selection. PbrbZIP genes were induced by cold and drought stresses and played important roles in drought and cold tolerance. These results provided useful information for further increasing the tolerance of pears to stresses and a foundation to study the cold and drought tolerance mechanism of PbrbZIP genes.

Genome-Wide Identification and Analysis of NAC Transcription Factor Family in Two Diploid Wild Relatives of Cultivated Sweet Potato Uncovers Potential NAC Genes Related to Drought Tolerance

NAC (NAM, ATAF1/2, and CUC2) proteins play a pivotal role in modulating plant development and offer protection against biotic and abiotic stresses. Until now, no systematic knowledge of NAC family genes is available for the food security crop, sweet potato. Here, a comprehensive genome-wide survey of NAC domain-containing proteins identified 130 ItbNAC and 144 ItfNAC genes with full length sequences in the genomes of two diploid wild relatives of cultivated sweet potato, Ipomoea triloba and Ipomoea trifida, respectively. These genes were physically mapped onto 15 I. triloba and 16 I. trifida chromosomes, respectively. Phylogenetic analysis divided all 274 NAC proteins into 20 subgroups together with NAC transcription factors (TFs) from Arabidopsis. There were 9 and 15 tandem duplication events in the I. triloba and I. trifida genomes, respectively, indicating an important role of tandem duplication in sweet potato gene expansion and evolution. Moreover, synteny analysis suggested that most NAC genes in the two diploid sweet potato species had a similar origin and evolutionary process. Gene expression patterns based on RNA-Seq data in different tissues and in response to various hormone, biotic or abiotic treatments revealed their possible involvement in organ development and response to various biotic/abiotic stresses. The expression of 36 NAC TFs, which were upregulated in the five tissues and in response to mannitol treatment, was also determined by real-time quantitative polymerase chain reaction (RT-qPCR) in hexaploid cultivated sweet potato exposed to drought stress. Those results largely corroborated the expression profile of mannitol treatment uncovered by the RNA-Seq data. Some significantly up-regulated genes related to drought stress, such as ItbNAC110, ItbNAC114, ItfNAC15, ItfNAC28, and especially ItfNAC62, which had a conservative spatial conformation with a closely related paralogous gene, ANAC019, may be potential candidate genes for a sweet potato drought tolerance breeding program. This analysis provides comprehensive and systematic information about NAC family genes in two diploid wild relatives of cultivated sweet potato, and will provide a blueprint for their functional characterization and exploitation to improve the tolerance of sweet potato to abiotic stresses.

Genome-Wide Identification of NAC Transcription Factor Family in Juglans mandshurica and Their Expression Analysis during the Fruit Development and Ripening

The NAC (NAM, ATAF and CUC) gene family plays a crucial role in the transcriptional regulation of various biological processes and has been identified and characterized in multiple plant species. However, genome-wide identification of this gene family has not been implemented in Juglans mandshurica, and specific functions of these genes in the development of fruits remain unknown. In this study, we performed genome-wide identification and functional analysis of the NAC gene family during fruit development and identified a total of 114 JmNAC genes in the J. mandshurica genome. Chromosomal location analysis revealed that JmNAC genes were unevenly distributed in 16 chromosomes; the highest numbers were found in chromosomes 2 and 4. Furthermore, according to the homologues of JmNAC genes in Arabidopsis thaliana, a phylogenetic tree was constructed, and the results demonstrated 114 JmNAC genes, which were divided into eight subgroups. Four JmNAC gene pairs were identified as the result of tandem duplicates. Tissue-specific analysis of JmNAC genes during different developmental stages revealed that 39 and 25 JmNAC genes exhibited upregulation during the mature stage in walnut exocarp and embryos, indicating that they may serve key functions in fruit development. Furthermore, 12 upregulated JmNAC genes were common in fruit ripening stage in walnut exocarp and embryos, which demonstrated that these genes were positively correlated with fruit development in J. mandshurica. This study provides new insights into the regulatory functions of JmNAC genes during fruit development in J. mandshurica, thereby improving the understanding of characteristics and evolution of the JmNAC gene family.

High-resolution transcription factor binding sites prediction improved performance and interpretability by deep learning method

Transcription factors (TFs) are essential proteins in regulating the spatiotemporal expression of genes. It is crucial to infer the potential transcription factor binding sites (TFBSs) with high resolution to promote biology and realize precision medicine. Recently, deep learning-based models have shown exemplary performance in the prediction of TFBSs at the base-pair level. However, the previous models fail to integrate nucleotide position information and semantic information without noisy responses. Thus, there is still room for improvement. Moreover, both the inner mechanism and prediction results of these models are challenging to interpret. To this end, the Deep Attentive Encoder-Decoder Neural Network (D-AEDNet) is developed to identify the location of TFs-DNA binding sites in DNA sequences. In particular, our model adopts Skip Architecture to leverage the nucleotide position information in the encoder and removes noisy responses in the information fusion process by Attention Gate. Simultaneously, the Transcription Factor Motif Discovery based on Sliding Window (TF-MoDSW), an approach to discover TFs-DNA binding motifs by utilizing the output of neural networks, is proposed to understand the biological meaning of the predicted result. On ChIP-exo datasets, experimental results show that D-AEDNet has better performance than competing methods. Besides, we authenticate that Attention Gate can improve the interpretability of our model by ways of visualization analysis. Furthermore, we confirm that ability of D-AEDNet to learn TFs-DNA binding motifs outperform the state-of-the-art methods and availability of TF-MoDSW to discover biological sequence motifs in TFs-DNA interaction by conducting experiment on ChIP-seq datasets.

Genome-wide identification and functional analysis of U-box E3 ubiquitin ligases gene family related to drought stress response in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The plant U-box (PUB) proteins are a family of ubiquitin ligases (E3) enzymes that involved in diverse biological processes, as well as in responses to plant stress response. However, the characteristics and functional divergence of the PUB gene family have not yet been previously studied in the Chinese white pear (Pyrus bretschneideri).

RESULTS

In the present study, we identified 62 PbrPUBs in Chinese white pear genome. Based on the phylogenetic relationship, 62 PUB genes were clustered into five groups. The results of conserved motif and gene structure analysis supported the classification phylogenetic tree. The PbrPUB genes were unevenly distribution on 17 pear chromosomes, chromosome 15 housed most member of PUB family, with eight PUB genes. Cis-acting element analysis indicated that PUB genes might participate in diverse biological processes, especially in the response to abiotic stresses. Based on RNA-data from 'Dangshansuli' at seven tissues, we found that PUB genes exhibited diverse of expression level in seven tissues, and qRT-PCR experiment further supported the reliable of RNA-Seq data. To identify candidate genes associated with resistance, we conducted qRT-PCR experiment the expression level of pear seed plant under four abiotic stresses, including: ABA, dehydration, salt and cold treatment. One candidate PUB gene associated with dehydration stress was selected to conduct further functional experiment. Subcellular localization revealed PbrPUB18 protein was located on cell nucleus. Furthermore, heterologous over-expression of PbrPUB18 in Arabidopsis indicated that the over-expression of PbrPUB18 could enhance resistance in drought treatment. In conclusions, we systematically identified the PUB genes in pear, and provided useful knowledge for functional identification of PUB genes in pear.

Transcriptome-wide and expression analysis of the NAC gene family in pepino (Solanum muricatum) during drought stress

Solanum muricatum (Pepino) is an increasingly popular solanaceous crop and is tolerant of drought conditions. In this study, 71 NAC transcription factor family genes of S. muricatum were selected to provide a theoretical basis for subsequent in-depth study of their regulatory roles in the response to biological and abiotic stresses, and were subjected to whole-genome analysis. The NAC sequences obtained by transcriptome sequencing were subjected to bioinformatics prediction and analysis. Three concentration gradient drought stresses were applied to the plants, and the target gene sequences were analyzed by qPCR to determine their expression under drought stress. The results showed that the S. muricatum NAC family contains 71 genes, 47 of which have conserved domains. The protein sequence length, molecular weight, hydrophilicity, aliphatic index and isoelectric point of these transcription factors were predicted and analyzed. Phylogenetic analysis showed that the S. muricatum NAC gene family is divided into seven subfamilies. Some NAC genes of S. muricatum are closely related to the NAC genes of Solanaceae crops such as tomato, pepper and potato. The seedlings of S. muricatum were grown under different gradients of drought stress conditions and qPCR was used to analyze the NAC expression in roots, stems, leaves and flowers. The results showed that 13 genes did not respond to drought stress while 58 NAC genes of S. muricatum that responded to drought stress had obvious tissue expression specificity. The overall expression levels in the root were found to be high. The number of genes at extremely significant expression levels was very large, with significant polarization. Seven NAC genes with significant responses were selected to analyze their expression trend in the different drought stress gradients. It was found that genes with the same expression trend also had the same or part of the same conserved domain. Seven SmNACs that may play an important role in drought stress were selected for NAC amino acid sequence alignment of Solanaceae crops. Four had strong similarity to other Solanaceae NAC amino acid sequences, and SmNAC has high homology with the Solanum pennellii. The NAC transcription factor family genes of S. muricatum showed strong structural conservation. Under drought stress, the expression of NAC transcription factor family genes of S. muricatum changed significantly, which actively responded to and participated in the regulation process of drought stress, thereby laying foundations for subsequent in-depth research of the specific functions of NAC transcription factor family genes of S. muricatum.

Genome-wide identification, characterization, and expression analysis of the NAC transcription factor family in orchardgrass (Dactylis glomerata L.)

Background

Orchardgrass (Dactylis glomerata L.) is one of the most important cool-season perennial forage grasses that is widely cultivated in the world and is highly tolerant to stressful conditions. However, little is known about the mechanisms underlying this tolerance. The NAC (NAM, ATAF1/2, and CUC2) transcription factor family is a large plant-specific gene family that actively participates in plant growth, development, and response to abiotic stress. At present, owing to the absence of genomic information, NAC genes have not been systematically studied in orchardgrass. The recent release of the complete genome sequence of orchardgrass provided a basic platform for the investigation of DgNAC proteins.

Results

Using the recently released orchardgrass genome database, a total of 108 NAC (DgNAC) genes were identified in the orchardgrass genome database and named based on their chromosomal location. Phylogenetic analysis showed that the DgNAC proteins were distributed in 14 subgroups based on homology with NAC proteins in Arabidopsis, including the orchardgrass-specific subgroup Dg_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 15, and multitudinous DgNAC genes contained three exons. Chromosomal mapping analysis found that the DgNAC genes were unevenly distributed on seven orchardgrass chromosomes. For the gene expression analysis, the expression levels of DgNAC genes in different tissues and floral bud developmental stages were quite different. Quantitative real-time PCR analysis showed distinct expression patterns of 12 DgNAC genes in response to different abiotic stresses. The results from the RNA-seq data revealed that orchardgrass-specific NAC exhibited expression preference or specificity in diverse abiotic stress responses, and the results indicated that these genes may play an important role in the adaptation of orchardgrass under different environments.

Conclusions

In the current study, a comprehensive and systematic genome-wide analysis of the NAC gene family in orchardgrass was first performed. A total of 108 NAC genes were identified in orchardgrass, and the expression of NAC genes during plant growth and floral bud development and response to various abiotic stresses were investigated. These results will be helpful for further functional characteristic descriptions of DgNAC genes and the improvement of orchardgrass in breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07485-6.

Genome-wide identification of PbrbHLH family genes, and expression analysis in response to drought and cold stresses in pear (Pyrus bretschneideri)

BACKGROUND

The basic helix-loop-helix (bHLH) transcription factors play important roles in many processes in plant growth, metabolism and responses to abiotic stresses. Although, the sequence of Chinese white pear genome (cv. 'Dangshansuli') has already been reported, there is still a lack of clarity regarding the bHLH family genes and their evolutionary history.

RESULTS

In this work, a genome-wide identification of the bHLH genes in Chinese white pear was performed, and we characterized the functional roles of these PbrbHLH genes in response to abiotic stresses. Based on the phylogenetic analysis and structural characteristics, 197 identified bHLH genes could be well classified into 21 groups. Expansion of PbrbHLH gene family was mainly driven by WGD and dispersed duplication with the purifying selection from the recent WGD. The functional annotation enrichment showed that the majority of PbrbHLHs were enriched in the GO terms and KEGG pathways involved in responds to stress conditions as TFs. Transcriptomic profiles and qRT-PCR revealed that PbrbHLH7, PbrbHLH8, PbrbHLH128, PbrbHLH160, PbrbHLH161 and PbrbHLH195 were significantly up-regulated under cold and drought treatments. In addition, PbrbHLH195-silenced pear seedlings display significant reduced cold tolerance, exhibiting reduced chlorophyll content, as well as increased electrolyte leakage and concentrations of malondialdehyde and H₂O₂.

CONCLUSION

For the first time, a comprehensive analysis identified the bHLH genes in Chinese white pear and demonstrated that PbrbHLH195 is involved in the production of ROS in response to cold stress, suggesting that members of the PbrbHLH family play an essential role in the stress tolerance of pear.

Genome-Wide Analysis of the Role of NAC Family in Flower Development and Abiotic Stress Responses in Cleistogenes songorica

Plant-specific NAC (NAM, ATAF, CUC) transcription factor (TF) family plays important roles in biological processes such as plant growth and response to stress. Nevertheless, no information is known about NAC TFs in Cleistogenes songorica, a prominent xerophyte desert grass in northwestern China. In this study, 162 NAC genes were found from the Cleistogenes songorica genome, among which 156 C. songoricaNAC (CsNAC) genes (96.3%) were mapped onto 20 chromosomes. The phylogenetic tree constructed by CsNAC and rice NAC TFs can be separated into 14 subfamilies. Syntenic and Ka/Ks analyses showed that CsNACs were primarily expanded by genomewide replication events, and purifying selection was the primary force driving the evolution of CsNAC family genes. The CsNAC gene expression profiles showed that 36 CsNAC genes showed differential expression between cleistogamous (CL) and chasmogamous (CH) flowers. One hundred and two CsNAC genes showed differential expression under heat, cold, drought, salt and ABA treatment. Twenty-three CsNAC genes were commonly differentially expressed both under stress responses and during dimorphic floret development. Gene Ontology (GO) annotation, coexpression network and qRT-PCR tests revealed that these CsNAC genes may simultaneously regulate dimorphic floret development and the response to stress. Our results may help to characterize the NAC transcription factors in C. songorica and provide new insights into the functional research and application of the NAC family in crop improvement, especially in dimorphic floret plants.

Genome-wide characterization of the abscisic acid-, stress- and ripening-induced (ASR) gene family in wheat (Triticum aestivum L.)

BACKGROUND

Abscisic acid-, stress-, and ripening-induced (ASR) genes are a class of plant specific transcription factors (TFs), which play important roles in plant development, growth and abiotic stress responses. The wheat ASRs have not been described in genome-wide yet.

METHODS

We predicted the transmembrane regions and subcellular localization using the TMHMM server, and Plant-mPLoc server and CELLO v2.5, respectively. Then the phylogeny tree was built by MEGA7. The exon-intron structures, conserved motifs and TFs binding sites were analyzed by GSDS, MEME program and PlantRegMap, respectively.

RESULTS

In wheat, 33ASR genes were identified through a genome-wide survey and classified into six groups. Phylogenetic analyses revealed that the TaASR proteins in the same group tightly clustered together, compared with those from other species. Duplication analysis indicated that the TaASR gene family has expanded mainly through tandem and segmental duplication events. Similar gene structures and conserved protein motifs of TaASRs in wheat were identified in the same groups. ASR genes contained various TF binding cites associated with the stress responses in the promoter region. Gene expression was generally associated with the expected group-specific expression pattern in five tissues, including grain, leaf, root, spike and stem, indicating the broad conservation of ASR genes function during wheat evolution. The qRT-PCR analysis revealed that several ASRs were up-regulated in response to NaCl and PEG stress.

CONCLUSION

We identified ASR genes in wheat and found that gene duplication events are the main driving force for ASR gene evolution in wheat. The expression of wheat ASR genes was modulated in responses to multiple abiotic stresses, including drought/osmotic and salt stress. The results provided important information for further identifications of the functions of wheat ASR genes and candidate genes for high abiotic stress tolerant wheat breeding.

Genome-wide analysis of Jatropha curcas MADS-box gene family and functional characterization of the JcMADS40 gene in transgenic rice

Background

Physic nut (Jatropha curcas), an inedible oilseed plant, is among the most promising alternative energy sources because of its high oil content, rapid growth and extensive adaptability. Proteins encoded by MADS-box family genes are important transcription factors participated in regulating plant growth, seed development and responses to abiotic stress. However, there has been no in-depth research on the MADS-box genes and their roles in physic nut.

Results

In our study, 63 MADS-box genes (JcMADSs) were identified in the physic nut genome, and classed into five groups (MIKC^C, Mα, Mβ, Mγ, MIKC*) according to phylogenetic comparison with Arabidopsis homologs. Expression profile analysis based on RNA-seq suggested that many JcMADS genes had the strongest expression in seeds, and seven of them responded in leaves to at least one abiotic stressor (drought and/or salinity) at one or more time points. Transient expression analysis and a transactivation assay indicated that JcMADS40 is a nucleus-localized transcriptional activator. Plants overexpressing JcMADS40 did not show altered plant growth, but the overexpressing plants did exhibit reductions in grain size, grain length, grain width, 1000-seed weight and yield per plant. Further data on the reduced grain size in JcMADS40-overexpressing plants supported the putative role of JcMADS genes in seed development.

Conclusions

This study will be useful in order to further understand the process of MADS-box genes involved in regulating growth and development in addition to their functions in abiotic stress resistance, and will eventually provide a theoretical basis for the functional investigation and the exploitation of candidate genes for the molecular improvement of physic nut.

Genome-wide analysis of NAC transcription factor family in maize under drought stress and rewatering

The plant-specific NAC transcription factor (TFs) plays crucial role in plant growth as well as in stress resistance. In the present study, 87 Zea mays NAC TFs were obtained from the transcriptome analysis using drought-resistant maize inbred line Y882 as experimental material under PEG stress and rewatering treatment. Comprehensive analyses were conducted including genes structure, chromosomal localization, phylogenetic tree and motif prediction, cis-elements and expression patterns. The results showed that the 87 ZmNAC genes distributed on 10 chromosomes and were categorized into 15 groups based on their conserved gene structure and motifs. Phylogenetic tree analysis was also constructed referencing to the counterparts of Arabidopsis and rice, and the stress-related cis-elements in the promoter region were also analyzed. 87 ZmNAC genes exhibited different expression levels at 3 treatment points, indicating different response to drought stress. This genome-wide analysis of 87 ZmNAC genes will provide basis for further gene function detection.

Identification and Expression Analysis of the NAC Gene Family in Coffea canephora

The NAC gene family is one of the largest families of transcriptional regulators in plants, and it plays important roles in the regulation of growth and development as well as in stress responses. Genome-wide analyses have been performed in diverse plant species, but there is still no systematic analysis of the NAC genes of Coffea canephora Pierre ex A. Froehner. In this study, we identified 63 NAC genes from the genome of C. canephora. The basic features and comparison analysis indicated that the NAC gene members increased via duplication events during the evolution of the plant. Phylogenetic analysis divided the NAC proteins from C. canephora, Arabidopsis and rice into 16 subgroups. Analysis of the expression patterns of CocNACs under cold stress and coffee bean development indicated that 38 CocNACs were differentially expressed under cold stress; six genes may play important roles in the process of cold acclimation, and four genes among 54 CocNACs showing a variety of expression patterns during different developmental stages of coffee beans may be positively related to the bean development. This study can expand our understanding of the functions of the CocNAC gene family in cold responses and bean development, thereby potentially intensifying the molecular breeding programs of Coffea spp. plants.

Genome-Wide Investigation of the NAC Gene Family and Its Potential Association with the Secondary Cell Wall in Moso Bamboo

NAC (NAM, ATAF, and CUC) transcription factors (TFs) are implicated in the transcriptional regulation of diverse processes and have been characterized in a number of plant species. However, NAC TFs are still not well understood in bamboo, especially their potential association with the secondary cell wall (SCW). Here, 94 PeNACs were identified and characterized in moso bamboo (Phyllostachys edulis). Based on their gene structures and conserved motifs, the PeNACs were divided into 11 groups according to their homologs in Arabidopsis. PeNACs were expressed variously in different tissues of moso bamboo, suggesting their functional diversity. Fifteen PeNACs associated with the SCW were selected for co-expression analysis and validation. It was predicted that 396 genes were co-expressed with the 15 PeNACs, in which 16 and 55 genes were involved in the lignin catabolic process and cellulose biosynthetic process respectively. As the degree of lignification in the growing bamboo shoots increased, all 15 PeNACs were upregulated with a trend of rising first and then decreasing except PeNAC37, which increased continuously. These results indicated that these PeNACs might play important roles in SCW biosynthesis and lignification in bamboo shoots. Seven of 15 PeNACs had been found positively co-expressed with seven PeMYBs, and they had similar expression patterns with those of the PeMYBs in bamboo shoots. The targeted sites of miR164 were found in 16 PeNACs, of which three PeNACs associated with SCW were validated to have an opposite expression trend to that of miR164 in growing bamboo shoots. In addition, three PeNACs were selected and verified to have self-activation activities. These results provide comprehensive information of the NAC gene family in moso bamboo, which will be helpful for further functional studies of PeNACs to reveal the molecular regulatory mechanisms of bamboo wood property.