| 
					
						|  |  
    					|  |  
    					| Genome-wide Analysis of the HD-ZIP_N Gene Family in Gossypium hirsutum |  
						| YE Si-Hong1,*, JI Mei-Jun2,*, CHEN Qi2, FENG Qun2, WANG Dong-Mei2, SUN Meng2, DENG Xiao-Nan1, WANG Rui1, HAN Wen-Bing1, WANG Bao-Hua2,** |  
						| 1 Cotton Research Institute, Anhui Academy of Agricultural Science, Hefei 230001, China; 2 School of Life Sciences, Nantong University, Nantong 226019, China
 |  
						|  |  
					
						| 
								
									| 
											
                        					 
												
													
													    |  |  
														| 
													
													    | Abstract  In previous research of our research group, a cotton (Gossypium) fiber strength candidate gene Gh_D09G1210 was identified by combining transcriptome sequencing and QTL mapping, which belongs to homeodomain-leucine zipper (HD-ZIP)__N gene family. The purpose of this study was to identify HD-ZIP_N gene family members in G. hirsutum and analyze their role in the formation of cotton fiber strength. Genome-wide analysis of HD-ZIP_N gene family in G. hirsutum was carried out, and the differential expression patterns of HD-ZIP_N genes in cotton fibers at different days post anthesis were analyzed by real-time quantitative PCR (qPCR). A total of 18 genes of the HD-ZIP_N gene family were identified in G. hirsutum genome conserved motif analysis and gene-promoter analysis showed that the gene family members had highly similar gene sequences and structural domains. The phylogenetic tree showed that all the gene family members could be divided into 3 subcategories. According to the gene location, HD-ZIP_N genes were located on 10 chromosomes in G. hirsutum genome. Most of the homologous genes had been purified and selected according to the Ka/Ks results. Collinearity analysis showed that the members of the HD-ZIP_N gene family distributed on homologous chromosomes had a strong connection. Cotton fibers of 17 and 21 d post-anthesis (dpa) were used to carry out qPCR verification experiments. The results showed that 3 genes Gh_D11G0809, Gh_D04G0859 and Gh_A11G0693 had a significant down-regulation trend in 'PD94042' with average fiber quality and 'IL9' with elite fiber quality (P<0.05), while Gh_A09G1204 had a significant up-regulation trend in 'IL9' (P<0.05), indicating that HD-ZIP_N gene family members had effects on regulation of fiber strength development. The findings of this study suggest that the HD-ZIP_N gene family members participate in regulating the development of G. hirsutum fiber strength, and provide genetic resources for improving cotton fiber strength through molecular breeding. |  
															| Received: 28 November 2022 |  
															|  |  
															| Corresponding Authors:
																**bhwang@ntu.edu.cn |  
															| About author:: *These authors contributed equally to this work |  |  |  |  
													
																												  
															| [1] 郭晓豪, 王寒涛, 魏鑫, 等. 2021. 基于两个陆地棉低世代群体定位纤维品质相关QTL[J]. 棉花学报, 33(1): 33-41. (Guo X H, Wang H T, Wei X, et al.2021. Mapping fiber quality-related QTLs based on two low-generation populations of upland cotton[J]. Cotton Science, 33(1): 33-41.)
 [2] 季美君, 曹孜怡, 王翌婷, 等. 2023. 陆地棉NCS1基因家族的全基因组成员鉴定及分析[J].浙江农业学报, 35(2): 249-258.
 (Ji M J, Cao Z Y, Wang Y T, et al.2023. Genome-wide identification and analysis of the NCS1 gene family in upland cotton[J]. Acta Agriculturae Zhejiangensis, 35(2): 249-258)
 [3] 李俐, 王义, 王康宇, 等. 2018. HD-Zip基因家族的结构、功能与表达模式[J]. 分子植物育种, 16(3): 781-790.
 (Li L, Wang Y, Wang K Y, et al.2018. Structure, function and expression pattern of HD-Zip gene family[J]. Molecular Plant Breeding, 16(3): 781-790.
 [4] 罗少, 唐翠明, 戴凡炜, 等. 2021. 植物转录因子HD-Zip基因家族参与逆境胁迫的研究进展[J]. 南方农业, 15(35): 173-175.
 (Luo S, Tang C M, Dai F W, et al.2021. Research progress of plant transcription factor HD-Zip gene family involved in stress[J]. South China Agriculture, 15(35): 173-175.)
 [5] 袁娜, 王彤, 刘廷利, 等. 2018. 棉花FAR1/FHY3基因家族的全基因组分析[J]. 棉花学报, 30(1): 1-11.
 (Yuan N, Wang T, Liu T L, et al.2018. Genome-wide analysis of the FAR1/FHY3 gene family in cotton[J]. Cotton Science, 30(1): 1-11)
 [6] Bailey T L, Williams N, Misleh C, et al.2006. MEME: discovering and analyzing DNA and protein sequence motifs[J]. Nucleic Acids Research, 34(2): 369-373.
 [7] Chen Q, Wang W, Khanal S, et al.2021. Transcriptome analysis reveals genes potentially related to high fiber strength in a Gossypium hirsutum line IL9 with Gossypium mustelinum introgression[J]. Genome, 64(11): 985-995.
 [8] Farshad R, Edwin G.2017. Regulatory function of homeodomain-leucine zipper (HD-ZIP) family proteins during embryogenesis[J]. The New Phytologist, 213(1): 95-104.
 [9] Ge H Q, Xu J J, Hua M Z, et al.2022. Genome-wide identification and analysis of ACP gene family in Sorghum bicolor (L.) Moench[J]. BMC Genomics, 23: 538.
 [10] Gong S H, Ding Y F, Hu S S, et al.2019. The role of HD-Zip class Ⅰ transcription factors in plant response to abiotic stresses[J]. Physiologia Plantarum, 167(4): 516-525.
 [11] Harris J C, Hrmova M, Lopato S, et al.2011. Modulation of plant growth by HD-Zip class Ⅰ and Ⅱ transcription factors in response to environmental stimuli[J]. The New Phytologist, 190(4): 823-837.
 [12] He G H, Liu P, Zhao H X, et al.2020. The HD-zip II transcription factors regulate plant architecture through the auxin pathway[J]. International Journal of Molecular Sciences, 21(9): 3250.
 [13] Hu B, Jin J P, Guo A Y, et al.2015. GSDS 2.0: An upgraded gene feature visualization serve[J]. Bioinformatics, 31(8): 1296-1297.
 [14] Islam S, Thyssen G N, Jenkins J N, et al.2016. A Magic population-based genome-wide association study reveals functional association of GhRBB1_A07 gene with superior fiber quality in cotton[J]. BMC genomics, 17(1): 903.
 [15] Ison J, Kalas M, Jonassen I, et al.2013. EDAM: An ontology of bioinformatics operations, types of data and identifiers, topics and formats[J]. Bioinformatics (Oxford, England), 29(10): 1325-1332.
 [16] Jackson S, Chen Z J.2010. Genomic and expression plasticity of polyploidy[J]. Current Opinion in Plant Biology, 13(2): 153-159.
 [17] Larkin M A, Blackshields G, Brown N P, et al.2007. Clustal W and Clustal X version 2.0[J]. Bioinformatics, 23(21): 2947-2948.
 [18] Lescot M, Déhais P, Thijs G, et al.2002. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences[J]. Nucleic Acids Research, 30(1): 325-327.
 [19] Li H T, Wang B, Zhang Q H, et al.2017. Genome-wide analysis of the auxin/indoleacetic acid (Aux/IAA) gene family in allotetraploid rapeseed (Brassica napus L.)[J]. BMC Plant Biology, 17(1): 204.
 [20] Li W, Dong J Y, Cao M X, et al.2019. Genome-wide identification and characterization of HD-ZIP genes in potato[J]. Gene, 697: 103-117.
 [21] Li Z, Gao Z Q, Li R H, et al.2020. Genome-wide identification and expression profiling of HD-ZIP gene family in Medicago truncatula[J]. Genomics, 112(5): 3624-3635.
 [22] May O L.1999. Registration of pd 94042 germplasm line of upland cotton with high yield and fiber maturity[J]. Crop Science, 39(2): 597-598.
 [23] Potter S C, Luciani A, Eddy S R, et al.2018. HMMER web server: 2018 update[J]. Nucleic Acids Research, 46(1): 200-204.
 [24] Raffaella C A, Andrea C, Samanta S, et al.2008. The Arabidopsis homeodomain-leucine zipper II gene family: Diversity and redundancy[J]. Plant Molecular Biology, 68(4-5): 465-478.
 [25] Sudhir K, Glen S, Koichiro T.2016. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology and Evolution, 33(7): 1870-1874.
 [26] Sun F D, Zhang J H, Wang S F, et al.2012. QTL mapping for fiber quality traits across multiple generations and environments in upland cotton[J]. Molecular Breeding, 30(1): 569-582.
 [27] Voorrips R E.2002. MapChart: Software for the graphical presentation of linkage maps and QTLs[J]. The Journal of Heredity, 93(1): 77-78.
 [28] William C J.2016. Auxin response factors[J]. Plant, Cell & Environment, 39(5): 1014-1028.
 [29] Zhang J.2020. Genome-wide analysis of PRR gene family uncovers their roles in circadian rhythmic changes and response to drought stress in Gossypium hirsutum L.[J]. Peer J, 8(4): 9936.
 [30] Zhang S J, Jiang Z X, Chen J, et al.2021. The cellulose synthase (CesA) gene family in four Gossypium species: phylogenetics, sequence variation and gene expression in relation to fiber quality in upland cotton[J]. Molecular Genetics and Genomics, 296(2): 355-368.
 [31] Zhang T Z, Hu Y, Jiang W K, et al.2015. Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement[J]. Nature Biotechnology, 33(5): 531-537.
 [32] Zhang Z, Ge Q, Liu A Y, et al.2017. Construction of a high-density genetic map and its application to QTL identification for fiber strength in upland cotton[J]. Crop Science, 57(2): 774-788.
 |  
											 
											 |  |  |