Hierarchical transcription factor and regulatory network for drought response in Betula platyphylla

Abstract

Although many genes and biological processes involved in abiotic stress responses have been identified, how they are regulated remains largely unclear. Here, to study the regulatory mechanism of birch (Betula platyphylla) responding to drought induced by polyethylene glycol 6000 (20%, w/v), a partial correlation coefficient-based algorithm for constructing a gene regulatory network (GRN) was proposed, and a three-layer hierarchical GRN was constructed, including 68 transcription factors and 252 structural genes. A total of 1448 predicted regulatory relationships are included, and most of them are novel. The reliability of the GRN was verified by chromatin immunoprecipitation (ChIP)–PCR and qRT–PCR based on transient transformation. About 55% of genes in the bottom layer of the GRN could confer drought tolerance. We selected two TFs, BpMADS11 and BpNAC090, from the top layer and characterized their function in drought tolerance. Overexpression of BpMADS11 and BpNAC090 reduces electrolyte leakage, reactive oxygen species (ROS) and malondialdehyde (MDA) contents, giving greater drought tolerance than wild-type birch. According to this GRN, the important biological processes involved in drought were identified, including ‘signaling hormone pathways’, ‘water transport’, ‘regulation of stomatal movement’, and ‘response to oxidative stress’. This work indicated that BpERF017, BpAGL61, and BpNAC090 are the key upstream regulators of birch drought tolerance. Our data clearly revealed that upstream regulators and transcription factor–DNA interaction regulate different biological processes to adapt to drought stress.

References

[1] Huang J, Yu H, Guan X et al. Accelerated dryland expansion under climate change. Nat Clim Chang. 2016; 6: 166-71.

[2] Yin H, Li M, Li D et al. Transcriptome analysis reveals regulatory framework for salt and osmotic tolerance in a succulent xerophyte. BMC Plant Biol. 2019; 19: 1-15.

[3] Thirumalaikumar VP, Devkar V, Mehterov N et al. NAC transcription factor JUNGBRUNNEN1 enhances drought tolerance in tomato. Plant Biotechnol J. 2018; 16: 354-66.

[4] Li D, Yang J, Pak S et al. PuC3H35 confers drought tolerance by enhancing lignin and proanthocyanidin biosynthesis in the roots of Populus ussuriensis. New Phytol. 2022; 233: 390-408.

[5] Chen P, Zhi F, Li X et al. Zinc-finger protein MdBBX7/MdCOL9, a target of MdMIEL1 E3 ligase, confers drought tolerance in apple. Plant Physiol. 2022; 188: 540-59.

[6] Dong Q, Duan D, Zheng W et al. Overexpression of MdVQ37 reduces drought tolerance by altering leaf anatomy and SA homeostasis in transgenic apple. Tree Physiol 2022; 42: 160-74.

[7] Li J, Zhao S, Yu X et al. Role of Xanthoceras sorbifolium MYB44 in tolerance to combined drought and heat stress via modulation of stomatal closure and ROS homeostasis. Plant Physiol Biochem. 2021; 162: 410-20.

[8] Yu H, Gerstein M . Genomic analysis of the hierarchical structure of regulatory networks. Proc Natl Acad Sci USA 2006; 103: 14724-31.

[9] Wu J, Zhao X, Lin Z, Shao Z . Large scale gene regulatory network inference with a multi-level strategy. Mol BioSyst. 2015; 12: 588.

[10] Gunasekara C, Zhang K, Deng W et al. TGMI: an efficient algorithm for identifying pathway regulators through evaluation of triple-gene mutual interaction. Nucleic Acids Res. 2018; 46: e67.

[11] Kumari S, Deng W, Gunasekara C et al. Bottom-up GGM algorithm for constructing multilayered hierarchical gene regulatory networks that govern biological pathways or processes. BMC Bioinformatics. 2016; 17: 132.

[12] Margolin AA, Nemenman I, Basso K et al. ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinform. 2006; 7: S7.

[13] Deng W, Zhang K, Busov V, Wei H . Recursive random forest algorithm for constructing multilayered hierarchical gene regulatory networks that govern biological pathways. PLoS One. 2017; 12: e0171532.

[14] Shi J, Shi J, Zhao B et al. A phosphate starvation response-centered network regulates mycorrhizal symbiosis. Cell. 2021; 184: 5527-5540.e18.

[15] Vermeirssen V, De Clercq I, Van Parys T et al. Arabidopsis ensemble reverse-engineered gene regulatory network discloses interconnected transcription factors in oxidative stress. Plant Cell. 2015; 26: 4656-79.

[16] Guo L, Zhao G, Xu J-R et al. Compartmentalized gene regulatory network of the pathogenic fungus Fusarium graminearum. New Phytol 2016; 211: 527-41.

[17] Chen Y, Chen Y, Shi Z et al. Biosynthesis and signal transduction of ABA, JA, and BRs in response to drought stress of Kentucky bluegrass. Int J Mol Sci. 2019; 20: 1289.

[18] Gao Z, Li J, Li L et al. Structural and functional analyses of hub MicroRNAs in an integrated gene regulatory network of arabidopsis. Genomics Proteomics Bioinformatics. 2021

[19] Chen S, Wang Y, Yu L . Genome sequence and evolution of Betula platyphylla. Hortic Res. 2021; 8: 37.

[20] Sun N, Zhao H . Reconstructing transcriptional regulatory networks through genomics data. Stat Methods Med Res. 2009; 18: 595-617.

[21] Zang D, Wang C, Ji X, Wang Y . Tamarix hispida zinc finger protein ThZFP1 participates in salt and osmotic stress tolerance by increasing proline content and SOD and POD activities. Plant Sci. 2015; 235: 111-21.

[22] Ruan J, Zhou Y, Zhou M et al. Jasmonic acid signaling pathway in plants. Int J Mol Sci. 2019; 20: 2479.

[23] de Ollas C, Hernando B, Arbona V, Gómez-Cadenas A . Jasmonic acid transient accumulation is needed for abscisic acid increase in citrus roots under drought stress conditions. Physiol Plant 2013; 147: 296-306.

[24] Park HY, Seok H-Y, Woo D-H et al. AtERF71/HRE2 transcription factor mediates osmotic stress response as well as hypoxia response in Arabidopsis. Biochem Biophys Res Commun. 2011; 414: 135-41.

[25] Zhao PX, Zhang J, Chen S-Y et al. Arabidopsis MADS-box factor AGL16 is a negative regulator of plant response to salt stress by downregulating salt-responsive genes. New Phytol. 2021; 232: 2418-39.

[26] Yoshida T, Mogami J, Yamaguchi-Shinozaki K . ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr Opin Plant Biol. 2014; 21: 133-9.

[27] Lee HG, Seo PJ . The MYB96-HHP module integrates cold and abscisic acid signaling to activate the CBF-COR pathway in Arabidopsis. Plant J. 2015; 82: 962-77.

[28] Zang D, Wang J, Zhang X et al. Arabidopsis heat shock transcription factor HSFA7b positively mediates salt stress tolerance by binding to an E-box-like motif to regulate gene expression. J Exp Bot. 2019; 70: 5355-74.

[29] Singh R, Parihar P, Singh S et al. Reactive oxygen species signaling and stomatal movement: current updates and future perspectives. Redox Biol. 2017; 11: 213-8.

[30] Zhao PX, Miao Z-Q, Zhang J et al. Arabidopsis MADS-box factor AGL16 negatively regulates drought resistance via stomatal density and stomatal movement. J Exp Bot. 2020; 71: 6092-106.

[31] Pan Y, Wu LJ, Yu ZL . Effect of salt and drought stress on antioxidant enzymes activities and SOD isoenzymes of liquorice (Glycyrrhiza uralensis Fisch). Plant Growth Regul. 2006; 49: 157-65.

[32] Bouvier FCAU, Backhaus RA, Camara B . Induction and control of chromoplast-specific carotenoid genes by oxidative stress. J Biol Chem. 1998; 273: 30651-9.

[33] Michel B, Kaufmann M . The osmotic potential of polyethylene glycol 6000. Plant Physiol. 1973; 51: 914-6.

[34] Li B, Dewey CN . RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011; 12: 323.

[35] Langmead B, Salzberg SL . Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9: 357-9.

[36] Robinson MD, McCarthy DJ, Smyth GK . EdgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009; 26: 139-40.

[37] Wang L, Mo Q, Wang J . MIrExpress: a database for gene coexpression correlation in immune cells based on mutual information and Pearson correlation. J Immunol Res. 2015; 2015.

[38] Obayashi T, Kinoshita K . Rank of correlation coefficient as a comparable measure for biological significance of gene coexpression. DNA Res. 2009; 16: 249-60.

[39] Zang D, Wang L, Zhang Y et al. ThDof1.4 and ThZFP1 constitute a transcriptional regulatory cascade involved in salt or osmotic stress in Tamarix hispida. Plant Mol Biol. 2017; 94: 495-507.

[40] Haring M, Offermann S, Danker T et al. Chromatin immunoprecipitation: optimization, quantitative analysis and data normalization. Plant Methods. 2007; 3: 11.

[41] Dionisio-Sese ML, Tobita S . Antioxidant responses of rice seedlings to salinity stress. Plant Sci. 1998; 135: 1-9.

[42] Bates LS, Waldren RP, Teare ID . Rapid determination of free proline for water-stress studies. Plant Soil 1973; 39: 205-7.

[43] Madhava Rao KV, Sresty TVS . Antioxidative parameters in the seedlings of pigeonpea (Cajanus cajan (L.) Millspaugh) in response to Zn and Ni stresses. Plant Sci. 2000; 157: 113-28.

[44] Hsieh E, Cheng M, Lin T . Functional characterization of an abiotic stress-inducible transcription factor AtERF53 in Arabidopsis thaliana. Plant Mol Biol. 2013; 82: 223-37.

[45] Livak KJ, Schmittgen TD . Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001; 25: 402-8.