Genome-wide Expression Analysis Research Articles

Genome-wide identification and expression analysis of Dmrt gene family and their role in gonad development of Pacific oyster (Crassostrea gigas).

Doublesex and Mab-3-related transcription factor (Dmrt) is a type of transcription factor with a zinc-finger DM structural domain, which plays a significant role in sex determination and differentiation in animals. Although Dmrt has been reported in many vertebrates and invertebrates, it has rarely been studied in bivalves. In this study, a total of three members of the Dmrt gene family were identified and characterized in Crassostrea gigas, and all these CgDmrt genes contained a conserved DM domain. Analysis of the phylogenetic tree and gene structure revealed that Dmrt genes clustered on one branch may have similar functions in bivalves. Expression profiling of CgDmrt mRNA in different tissues and stages of gonad development indicated that CgDmrt3 exhibited sexually dimorphic expression and played an important role in the development of the male gonad in C. gigas. Furthermore, analysis of CgDmrt mRNA expression between fertile triploids and sterile triploids showed that CgDmrt3 may be involved in sperm production. Collectively, the systematic analysis of the CgDmrt genes will provide potential insights into the function of these genes in gonadal development.

Genome-wide characterization, phylogenetic and expression analysis of Galectin gene family in Golden pompano Trachinotus ovatus.

Galectins (Gals) are a type of S-type lectin that are widespread and evolutionarily conserved among metazoans, and can act as pattern recognition receptors (PRRs) to recognize pathogen-associated molecular patterns (PAMPs). In this study, 10 Gals (ToGals) were identified in the Golden pompano (Trachinotus ovatus), and their conserved domains, motifs, and collinearity relationships were analyzed. The expression of ToGals was regulated following infection to Cryptocaryon irritans and Streptococcus agalactiae, indicating that ToGals participate in immune responses against microbial pathogens. Further analysis was conducted on one important member, Galectin-3, subcellular localization showing that ToGal-3like protein is expressed both in the nucleus and cytoplasm. Recombinant protein obtained through prokaryotic expression showed that rToGal-3like can agglutinate red blood cells of rabbit, carp and golden pompano and also agglutinate and kill Staphylococcus aureus, Bacillus subtilis, Vibrio vulnificus, S. agalactiae, Pseudomonas aeruginosa, and Aeromonas hydrophila. This study lays the foundation for further research on the immune roles of Gals in teleosts.

Genome-wide identification, characterization and expression analysis of the bZIP transcription factors in garlic (Allium sativum L.).

The bZIP genes (bZIPs) are essential in numerous biological processes, including development and stress responses. Despite extensive research on bZIPs in many plants, a comprehensive genome-wide analysis of bZIPs in garlic has yet to be undertaken. In this study, we identified and classified 64 AsbZIP genes (AsbZIPs) into 10 subfamilies. A systematic analysis of the evolutionary characteristics of these AsbZIPs, including chromosome location, gene structure, conserved motifs, and gene duplication, was conducted. Furthermore, we also examined the nucleotide diversity, cis-acting elements, and expression profiles of AsbZIPs in various tissues and under different abiotic stresses and hormone treatments. Our findings revealed that gene replication plays a crucial role in the expansion of AsbZIPs, with a minor genetic bottleneck observed during domestication. Moreover, the identification of cis-acting elements suggested potential associations of AsbZIPs with garlic development, hormone, and stress responses. Several AsbZIPs exhibited tissue-preferential and stress/hormone-responsive expression patterns. Additionally, Asa7G01972 and Asa7G01379 were notably differentially expressed under various stresses and hormone treatments. Subsequent yeast two-hybridization and yeast induction experiments validated their interactions with Asa1G01577, a homologue of ABI5, reinforcing their importance in hormone and abiotic stress responses. This study unveiled the characteristics of the AsbZIP superfamily and lays a solid foundation for further functional analysis of AsbZIP in garlic.

Genome-wide identification and expression analysis of TCP gene family of Andrographis paniculata under abiotic stress

Andrographis paniculata is an important medicinal plant in the Lingnan region of China, which has the functions of clearing heat, removing toxins, and resisting bacteria and inflammation. The TCP gene family is a class of transcription factors that regulate plant growth, development, and stress response. In order to analysis the role of the TCP gene family under abiotic stress in A. paniculata, this study identified the TCP gene family of A. paniculata at the genome-wide level and analyzed its expression pattern in response to abiotic stress. The results showed that the A. paniculata TCP gene family had 23 members, with length of amino acid ranging from 136 to 508, the relative molecular mass between 14 854.71 and 55 944.90 kDa, and the isoelectric point between 5.67 and 10.39. All members were located in the nucleus and unevenly distributed on 13 chromosomes. Phylogenetic analysis classified them into three subfamilies: PCF, CIN and CYC/TB1. Gene structure and conserved motif analysis showed that most members of the TCP gene family contained motif 1, motif 2, motif 3 in the same order and 1-3 CDS. The analysis of promoter cis-acting elements showed that the transcriptional expression of the TCP gene family in A. paniculata might be induced by light, hormones, and adversity stress. In light of the expression pattern analysis and qRT-PCR verification, the expression of ApTCP4, ApTCP5, ApTCP6, and ApTCP11 involved in response by various abiotic stresses such as drought, high temperature, and MeJA. This study lays the foundation for in-depth exploration of the functions of A. paniculata TCP genes in response to abiotic stress.

Genome-Wide Identification and Expression Analysis of the <i>GSK3</i> Gene Family in Sunflower under Various Abiotic Stresses

Genome-Wide Identification and Expression Analysis of the <i>GSK3</i> Gene Family in Sunflower under Various Abiotic Stresses

Genome-wide identification and expression analysis of TPP gene family under salt stress in peanut (Arachis hypogaea L.).

Trehalose-6-phosphate phosphatase (TPP), a key enzyme for trehalose biosynthesis in plants, plays a pivotal role in the growth and development of higher plants, as well as their adaptations to various abiotic stresses. Employing bioinformatics techniques, 45 TPP genes distributed across 17 chromosomes were identified with conserved Trehalose-PPase domains in the peanut genome, aiming to screen those involved in salt tolerance. Collinearity analysis showed that 22 TPP genes from peanut formed homologous gene pairs with 9 TPP genes from Arabidopsis and 31 TPP genes from soybean, respectively. Analysis of cis-acting elements in the promoters revealed the presence of multiple hormone- and abiotic stress-responsive elements in the promoter regions of AhTPPs. Expression pattern analysis showed that members of the TPP gene family in peanut responded significantly to various abiotic stresses, including low temperature, drought, and nitrogen deficiency, and exhibited certain tissue specificity. Salt stress significantly upregulated AhTPPs, with a higher number of responsive genes observed at the seedling stage compared to the podding stage. The intuitive physiological effect was reflected in the significantly higher accumulation of trehalose content in the leaves of plants under salt stress compared to the control. These findings indicate that the TPP gene family plays a crucial role in peanut's response to abiotic stresses, laying the foundation for further functional studies and utilization of these genes.

Genome-Wide Identification and Expression Analysis of the Trehalose-6-phosphate Synthase and Trehalose-6-phosphate Phosphatase Gene Families in Rose (Rosa hybrida cv 'Carola') under Different Light Conditions.

Trehalose, trehalose-6-phosphate synthase (TPS),and trehalose-6-phosphatase (TPP) have been reported to play important roles in plant abiotic stress and growth development. However, their functions in the flowering process of Rosa hybrida have not been characterized. In this study we found that, under a short photoperiod or weak light intensity, the content of trehalose in the shoot apical meristem of Rosa hybrida cv 'Carola' significantly decreased, leading to delayed flowering time. A total of nine RhTPSs and seven RhTPPs genes were identified in the genome. Cis-element analysis suggested that RhTPS and RhTPP genes were involved in plant hormones and environmental stress responses. Transcriptome data analysis reveals significant differences in the expression levels of RhTPSs and RhTPPs family genes in different tissues and indicates that RhTPPF and RhTPPJ are potential key genes involved in rose flower bud development under different light environments. The results of quantitative real-time reverse transcription (qRT-PCR) further indicate that under short photoperiod and weak light intensity all RhTPP members were significantly down-regulated. Additionally, RhTPS1a, RhTPS10, and RhTPS11 were up-regulated under a short photoperiod and showed a negative correlation with flowering time and trehalose content decrease. Under weak light intensity, RhTPS11 was up-regulated and negatively regulated flowering, while RhTPS5, RhTPS6, RhTPS7b, RhTPS9, and RhTPS10 were down-regulated and positively regulated flowering. This work lays the foundation for revealing the functions of RhTPS and RhTPP gene families in the regulation of rose trehalose.

Genome-Wide Identification and Expression Analysis of Calmodulin and Calmodulin-like Genes, Revealing CaM3 and CML13 Participating in Drought Stress in Phoebe bournei.

Calmodulin (CaM) and calmodulin-like (CML) proteins are major Ca2+ sensors involved in the regulation of plant development and stress responses by converting Ca2+ signals into appropriate cellular responses. However, characterization and expression analyses of CaM/CML genes in the precious species, Phoebe bournei, remain limited. In this study, five PbCaM and sixty PbCML genes were identified that only had EF-hand motifs with no other functional domains. The phylogenetic tree was clustered into 11 subgroups, including a unique clade of PbCaMs. The PbCaMs were intron-rich with four EF-hand motifs, whereas PbCMLs had two to four EF-hands and were mostly intronless. PbCaMs/CMLs were unevenly distributed across the 12 chromosomes of P. bournei and underwent purifying selection. Fragment duplication was the main driving force for the evolution of the PbCaM/CML gene family. Cis-acting element analysis indicated that PbCaMs/CMLs might be related to hormones, growth and development, and stress response. Expression analysis showed that PbCaMs were generally highly expressed in five different tissues and under drought stress, whereas PbCMLs showed specific expression patterns. The expression levels of 11 candidate PbCaMs/CMLs were responsive to ABA and MeJA, suggesting that these genes might act through multiple signaling pathways. The overexpression of PbCaM3/CML13 genes significantly increased the tolerance of yeast cells to drought stress. The identification and characterization of the CaM/CML gene family in P. bournei laid the foundation for future functional studies of these genes.

Genome-wide identification and expression analysis of the chlorophyll a/b binding protein gene family in oilseed (Brassica napus L.) under salt stress conditions

To elucidate the functions and molecular mechanisms of BnLhc gene family members in Brassica napus under different exogenous hormones and abiotic stress, we identified 57 BnLhc gene family members from the entire genome of B. napus variety ZS11 using Arabidopsis Lhc protein as a seed sequence combined with phyll A-B binding protein PF00504. We performed comprehensive analyses including physicochemical property analysis, structural analysis, phylogenetic relationship analysis, homology relationship analysis, protein interaction network analysis, microRNA-target gene interaction analysis and GO enrichment analysis. Additionally, based on transcriptome data, we investigated response patterns of BnLhc genes to various exogenous stress conditions. Comparative evaluation between non-salt (NS:0) and salt-treated (TS: 1) rapeseed plants revealed significant inhibition in leaf chlorophyll content, biomass (W), leaf area (LA), root surface area (RA), as well as chloroplasts and roots under salt treatment. Furthermore, qPCR was employed to assess the relative expression levels of BnLhc genes associated with salt tolerance. Our findings indicate that while BnLhcs were upregulated in leaves upon exposure to salt signals suggesting their potential involvement in enhancing rape's resistance to salt stress through ion redistribution processes; they were downregulated in roots. This study provides valuable insights into the functionality and molecular genetic mechanisms underlying the role BnLhc gene family

Genome-Wide Identification and Expression Analysis of the SUT Family from Three Species of Sapindaceae Revealed Their Role in the Accumulation of Sugars in Fruits.

Sapindaceae is an economically important family of Sapindales and includes many fruit crops. The dominant transport and storage form of photoassimilates in higher plants is sucrose. Sucrose transporter proteins play an irreplaceable role in the loading, transportation, unloading, and distribution of sucrose. A few SUT (sugar transporter) family genes have been identified and characterized in various plant species. In this study, 15, 15, and 10 genes were identified in litchi, longan, and rambutan, respectively, via genome-wide screening. These genes were divided into four subgroups based on phylogenetics. Gene duplication analysis suggested these genes underwent potent purifying selection and tandem duplications during evolution. The expression levels of SlSut01 and SlSut08 were significantly increased in the fruits of Sapindaceae members. The homologs of these two genes in longan and rambutan were also highly expressed in the fruits. The expression pattern of SUTs in three organs of the two varieties was also explored. Subcellular colocalization experiments revealed that the proteins encoded by both genes were present in the plasma membrane. This report provides data for the functional study of SUTs in litchi and provides a basis for screening sugar accumulation-related genes in fruits of Sapindaceae.

Genome-Wide Identification and Expression Analysis of Auxin Response Factor (ARF) Gene Family in Betula pendula

As the key transcription factors regulating auxin responsive genes expression, auxin response factors (ARFs) play critical roles in diverse aspects of plant growth and development. Betula pendula is a valuable ornamental tree, and the information on ARF gene family of B. pendula is needed for better understanding. The publication of the genome sequence of B. pendula enable to analyze the bioinformatics information and expression pattern of BpeARF gene family on the genome-wide basis. In this study, physical and chemical properties, chromosome location, phylogenetic relationship, gene structure, conserved domain, motif composition, and cis-acting element of BpeARF gene family were analyzed, and expression patterns of BpeARF genes were investigated using completely random design in different tissues and under exogenous NAA and drought treatments. A total of 17 BpeARF genes was identified from B. pendula genome, which were unevenly distributed on 13 chromosomes and encoded adequate proteins ranging from 613 to 1135 amino acids in length. Three BpeARF gene pairs were formed by segmental duplication, and the Ka/Ks values of these BpeARF gene pairs were less than 1. According to the phylogenetic relationship among B. pendula, Betula platyphylla, Populus trichocarpa, and Arabidopsis thaliana, the BpeARF genes were divided into four classes, and the intron/exon structure, conserved domain, and motif composition showed high similarity among the BpeARF genes within the same class. The cis-acting elements in the promoter regions of BpeARF genes were related to tissue development, hormone response, and stress resistance. Quantitative real-time PCR exhibited diverse expression patterns of BpeARF genes in different tissues and in response to exogenous auxin treatment and drought stress. The expressions of one, ten, seven, and three BpeARF genes were the high levels in buds, young leaves, stems, and roots, respectively. Under exogenous NAA treatment, six BpeARF genes in stems and roots were upregulated expression at all timepoints. Under drought stress, BpeARF7 and BpeARF15 were upregulated in stems and roots, and BpeARF5 and BpeARF6 were downregulated in leaves, stems, and roots. Our results provided valuable information for the classification and putative functions of BpeARF gene family, which may be helpful for selecting candidate genes and verifying gene function in the genetic engineering of birch trees in further research.

Genome-Wide Identification, Expression and Interaction Analyses of PP2C Family Genes in Chenopodium quinoa.

Plant protein phosphatase 2Cs (PP2Cs) function as inhibitors in protein kinase cascades involved in various processes and are crucial participants in both plant development and signaling pathways activated by abiotic stress. In this study, a genome-wide study was conducted on the CqPP2C gene family. A total of putative 117 CqPP2C genes were identified. Comprehensive analyses of physicochemical properties, chromosome localization and subcellular localization were conducted. According to phylogenetic analysis, CqPP2Cs were divided into 13 subfamilies. CqPP2Cs in the same subfamily had similar gene structures, and conserved motifs and all the CqPP2C proteins had the type 2C phosphatase domains. The expansion of CqPP2Cs through gene duplication was primarily driven by segmental duplication, and all duplicated CqPP2Cs underwent evolutionary changes guided by purifying selection. The expression of CqPP2Cs in various tissues under different abiotic stresses was analyzed using RNA-seq data. The findings indicated that CqPP2C genes played a role in regulating both the developmental processes and stress responses of quinoa. Real-time quantitative reverse transcription PCR (qRT-PCR) analysis of six CqPP2C genes in subfamily A revealed that they were up-regulated or down-regulated under salt and drought treatments. Furthermore, the results of yeast two-hybrid assays revealed that subfamily A CqPP2Cs interacted not only with subclass III CqSnRK2s but also with subclass II CqSnRK2s. Subfamily A CqPP2Cs could interact with CqSnRK2s in different combinations and intensities in a variety of biological processes and biological threats. Overall, our results will be useful for understanding the functions of CqPP2C in regulating ABA signals and responding to abiotic stress.

Genome-Wide Identification and Expression Analysis of the SBP-Box Gene Family in Loquat Fruit Development.

The loquat (Eriobotrya japonica L.) is a special evergreen tree, and its fruit is of high medical and health value as well as having stable market demand around the world. In recent years, research on the accumulation of nutrients in loquat fruit, such as carotenoids, flavonoids, and terpenoids, has become a hotspot. The SBP-box gene family encodes transcription factors involved in plant growth and development. However, there has been no report on the SBP-box gene family in the loquat genome and their functions in carotenoid biosynthesis and fruit ripening. In this study, we identified 28 EjSBP genes in the loquat genome, which were unevenly distributed on 12 chromosomes. We also systematically investigated the phylogenetic relationship, collinearity, gene structure, conserved motifs, and cis-elements of EjSBP proteins. Most EjSBP genes showed high expression in the root, stem, leaf, and inflorescence, while only five EjSBP genes were highly expressed in the fruit. Gene expression analysis revealed eight differentially expressed EjSBP genes between yellow- and white-fleshed fruits, suggesting that the EjSBP genes play important roles in loquat fruit development at the breaker stage. Notably, EjSBP01 and EjSBP19 exhibited completely opposite expression patterns between white- and yellow-fleshed fruits during fruit development, and showed a close relationship with SlCnr involved in carotenoid biosynthesis and fruit ripening, indicating that these two genes may participate in the synthesis and accumulation of carotenoids in loquat fruit. In summary, this study provides comprehensive information about the SBP-box gene family in the loquat, and identified two EjSBP genes as candidates involved in carotenoid synthesis and accumulation during loquat fruit development.

Genome-Wide Identification, Characterization, and Expression Analysis of the Amino Acid Permease Gene Family in Soybean

Amino acid permeases (AAPs) play important roles in transporting amino acids in plant species, leading to increased low-nitrogen tolerance, grain yield, or protein content. However, very few AAPs have been characterized in soybean (Glycine max). In this study, we scanned the soybean reference genome and identified a total of 36 AAP genes (named GmAAP). The GmAAPs were phylogenetically divided into three evolutionary clades, with the genes in the same clades sharing similar gene structures and domain organization. We also showed that seventeen GmAAP genes on ten chromosomes were in collinearity, likely due to whole-genome duplication. Further analysis revealed a variety of cis-acting regulatory elements (such as hormone response elements (ABRE, ERE, GARE, P-box, and TGA-element), stress response elements (LTR, MBS, MYB-related components, TC-rich repeats, TCA-element, and WUN-motif), the tissue expression element (GCN4-motif), and the circadian regulatory element (circadian) present in the 2 kb region of the GmAAP promoter region, demonstrating functional diversity and expression specificity. RNA-Seq data and quantitative real-time PCR identified five GmAAPs showing differential expression under nitrogen limitation, including GmAAP3, GmAAP5, and GmAAP8 showing downregulation while GmAAP14, GmAAP29 showed upregulation, suggesting their involvement in low-nitrogen stress response. These results provide comprehensive information on soybean AAP genes in nitrogen stress, and provide putative candidates with possible roles in enhancing amino acid delivery to seeds for yield improvement.

The CDPK Gene Family in Mustard (Brassica juncea L.): Genome-Wide Identification and Expression Analysis under Cold Stress

The CDPK family genes play crucial roles in signal transduction pathways during plant development and stress response. In this study, we comprehensively analyzed the CDPK family genes in mustard (Brassica juncea L.), resulting in the identification of 101 genes (BjuCDPK1–101) located on chromosomes AA_Chr01 to BB_Chr08. RNA-seq analysis showed that most BjuCDPK genes were predominantly expressed in the root and flower bud, suggesting their organ-specific expression patterns. Furthermore, the expression levels of most BjuCDPK genes were significantly altered by cold stress. The promoter regions of BjuCDPK genes had abundant cis-acting elements related to hormones and stress. The interaction network analysis predicted that BjuCDPK proteins respond to abiotic stress probably via interaction with ABF1, ABF4, and/or DI19. Our findings provide valuable information for further the investigation of the cold stress adaption of mustard via the CDPK signaling pathway.

The SWEET gene family in watermelon: genome-wide identification, phylogeny, duplication and expression analyses in male sterile buds and under sucrose, fructose and glucose treatments

The SWEET gene family in watermelon: genome-wide identification, phylogeny, duplication and expression analyses in male sterile buds and under sucrose, fructose and glucose treatments

Genome-Wide Identification and Expression Analysis of the Stearoyl-Acyl Carrier Protein Δ9 Desaturase Gene Family under Abiotic Stress in Barley.

Stearoyl-acyl carrier protein (ACP) Δ9 desaturase (SAD) is a critical fatty acid dehydrogenase in plants, playing a prominent role in regulating the synthesis of unsaturated fatty acids (UFAs) and having a significant impact on plant growth and development. In this study, we conducted a comprehensive genomic analysis of the SAD family in barley (Hordeum vulgare L.), identifying 14 HvSADs with the FA_desaturase_2 domain, which were divided into four subgroups based on sequence composition and phylogenetic analysis, with members of the same subgroup possessing similar genes and motif structures. Gene replication analysis suggested that tandem and segmental duplication may be the major reasons for the expansion of the SAD family in barley. The promoters of HvSADs contained various cis-regulatory elements (CREs) related to light, abscisic acid (ABA), and methyl jasmonate (MeJA). In addition, expression analysis indicated that HvSADs exhibit multiple tissue expression patterns in barley as well as different response characteristics under three abiotic stresses: salt, drought, and cold. Briefly, this evolutionary and expression analysis of HvSADs provides insight into the biological functions of barley, supporting a comprehensive analysis of the regulatory mechanisms of oil biosynthesis and metabolism in plants under abiotic stress.

Genome-wide identification and expression analysis of the glutamate receptor gene family in sweet potato and its two diploid relatives.

Plant glutamate receptor (GLR) homologs are crucial calcium channels that play an important role in plant development, signal transduction, and response to biotic and abiotic stresses. However, the GLR gene family has not yet been thoroughly and systematically studied in sweet potato. In this study, a total of 37 GLR genes were identified in the cultivated hexaploid sweet potato (Ipomoea batatas), and 32 GLR genes were discovered in each of the two diploid relatives (Ipomoea trifida and Ipomoea triloba) for the first time. Based on their evolutionary relationships to those of Arabidopsis, these GLRs were split into five subgroups. We then conducted comprehensive analysis to explore their physiological properties, protein interaction networks, promoter cis-elements, chromosomal placement, gene structure, and expression patterns. The results indicate that the homologous GLRs of the cultivated hexaploid sweet potato and its two relatives are different. These variations are reflected in their functions related to plant growth, hormonal crosstalk, development of tuberous roots, resistance to root rot, and responses to abiotic stress factors, all of which are governed by specific individual GLR genes. This study offers a comprehensive analysis of GLR genes in sweet potato and its two diploid relatives. It also provides a theoretical basis for future research into their regulatory mechanisms, significantly influencing the field of molecular breeding in sweet potatoes.

Genome-Wide Identification and Expression Analysis Reveals the B3 Superfamily Involved in Embryogenesis and Hormone Responses in Dimocarpus longan Lour.

B3 family transcription factors play an essential regulatory role in plant growth and development processes. This study performed a comprehensive analysis of the B3 family transcription factor in longan (Dimocarpus longan Lour.), and a total of 75 DlB3 genes were identified. DlB3 genes were unevenly distributed on the 15 chromosomes of longan. Based on the protein domain similarities and functional diversities, the DlB3 family was further clustered into four subgroups (ARF, RAV, LAV, and REM). Bioinformatics and comparative analyses of B3 superfamily expression were conducted in different light and with different temperatures and tissues, and early somatic embryogenesis (SE) revealed its specific expression profile and potential biological functions during longan early SE. The qRT-PCR results indicated that DlB3 family members played a crucial role in longan SE and zygotic embryo development. Exogenous treatments of 2,4-D (2,4-dichlorophenoxyacetic acid), NPA (N-1-naphthylphthalamic acid), and PP333 (paclobutrazol) could significantly inhibit the expression of the DlB3 family. Supplementary ABA (abscisic acid), IAA (indole-3-acetic acid), and GA3 (gibberellin) suppressed the expressions of DlLEC2, DlARF16, DlTEM1, DlVAL2, and DlREM40, but DlFUS3, DlARF5, and DlREM9 showed an opposite trend. Furthermore, subcellular localization indicated that DlLEC2 and DlFUS3 were located in the nucleus, suggesting that they played a role in the nucleus. Therefore, DlB3s might be involved in complex plant hormone signal transduction pathways during longan SE and zygotic embryo development.

Comprehensive genome‑wide analysis of the chicken heat shock protein family: identification, genomic organization, and expression profiles in indigenous chicken with highly pathogenic avian influenza infection

BackgroundHeat shock proteins (HSPs) function as molecular chaperones with critical roles in chicken embryogenesis, immune response to infectious diseases, and response to various environmental stresses. However, little is known on HSP genes in chicken. In this study, to understand the roles of chicken HSPs, we performed genome-wide identification, expression, and functional analyses of the HSP family genes in chicken.ResultsA total of 76 HSP genes were identified in the chicken genome, which were further classified into eight distinct groups (I-VIII) based on phylogenetic tree analysis. The gene-structure analysis revealed that the members of each clade had the same or similar exon-intron structures. Chromosome mapping suggested that HSP genes were widely dispersed across the chicken genome, except in chromosomes 16, 18, 22, 25, 26, and 28–32, which lacked chicken HSP genes. On the other hand, the interactions among chicken HSPs were limited, indicating that the remaining functions of HSPs could be investigated in chicken. Moreover, KEGG pathway analysis showed that the HSP gene family was involved in the regulation of heat stress, apoptotic, intracellular signaling, and immune response pathways. Finally, RNA sequencing data revealed that, of the 76 chicken HSP genes, 46 were differentially expressed at 21 different growth stages in chicken embryos, and 72 were differentially expressed on post-infection day 3 in two indigenous Ri chicken lines infected with highly pathogenic avian influenza.ConclusionsThis study provides significant insights into the potential functions of HSPs in chicken, including the regulation of apoptosis, heat stress, chaperone activity, intracellular signaling, and immune response to infectious diseases.