[Genome-wide identification of bZIP family genes and screening of candidate AarbZIPs involved in terpenoid biosynthesis in Artemisia argyi].

Artemisia argyi is an important medicinal and economic plant in China, with the effects of warming channels, dispersing cold, and relieving pain, inflammation, and allergy. The essential oil of this plant is rich in volatile terpenoids and widely used in moxi-bustion and healthcare products, with huge market potential. The bZIP transcription factors compose a large family in plants and are involved in the regulation of plant growth and development, stress response, and biosynthesis of secondary metabolites such as terpenoids. However, little is known about the bZIPs and their roles in A. argyi. In this study, the bZIP transcription factors in the genome of A. argyi were systematically identified, and their physicochemical properties, phylogenetic relationship, conserved motifs, and promoter-binding elements were analyzed. Candidate AarbZIP genes involved in terpenoid biosynthesis were screened out. The results showed that a total of 156 AarbZIP transcription factors were identified at the genomic level, with the lengths of 99-618 aa, the molecular weights of 11.7-67.8 kDa, and the theoretical isoelectric points of 4.56-10.16. According to the classification of bZIPs in Arabidopsis thaliana, the 156 AarbZIPs were classified into 12 subfamilies, and the members in the same subfamily had similar conserved motifs. The cis-acting elements of promoters showed that AarbZIP genes were possibly involved in light and hormonal pathways. Five AarbZIP genes that may be involved in the regulation of terpenoid biosynthesis were screened out by homologous alignment and phylogenetic analysis. The qRT-PCR results showed that the expression levels of the five AarbZIP genes varied significantly in different tissues of A. argyi. Specifically, AarbZIP29 and AarbZIP55 were highly expressed in the leaves and AarbZIP81, AarbZIP130, and AarbZIP150 in the flower buds. This study lays a foundation for the functional study of bZIP genes and their regulatory roles in the terpenoid biosynthesis in A. argyi.

Approaches to increase the validity of gene family identification using manual homology search tools.

Identifying homologs is an important process in the analysis of genetic patterns underlying traits and evolutionary relationships among species. Analysis of gene families is often used to form and support hypotheses on genetic patterns such as gene presence, absence, or functional divergence which underlie traits examined in functional studies. These analyses often require precise identification of all members in a targeted gene family. Manual pipelines where homology search and orthology assignment tools are used separately are the most common approach for identifying small gene families where accurate identification of all members is important. The ability to curate sequences between steps in manual pipelines allows for simple and precise identification of all possible gene family members. However, the validity of such manual pipeline analyses is often decreased by inappropriate approaches to homology searches including too relaxed or stringent statistical thresholds, inappropriate query sequences, homology classification based on sequence similarity alone, and low-quality proteome or genome sequences. In this article, we propose several approaches to mitigate these issues and allow for precise identification of gene family members and support for hypotheses linking genetic patterns to functional traits.

[Genome-wide identification of the banana GLR gene family and its expression analysis in response to low temperature and abscisic acid/methyl jasmonate].

Glutamate receptor-like (GLR) is an important class of Ca2+ channel proteins, playing important roles in plant growth and development as well as in response to biotic and abiotic stresses. In this paper, we performed genome-wide identification of banana GLR gene family based on banana genomic data. Moreover, we analyzed the basic physicochemical properties, gene structure, conserved motifs, promoter cis-acting elements, evolutionary relationships, and used real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) to verify the expression patterns of some GLR family members under low temperature of 4 ℃ and different hormone treatments. The results showed that there were 19 MaGLR family members in Musa acuminata, 16 MbGLR family members in Musa balbisiana and 14 MiGLR family members in Musa itinerans. Most of the members were stable proteins and had signal peptides, all of them had 3-6 transmembrane structures. Prediction of subcellular localization indicated that all of them were localized on the plasma membrane and irregularly distributed on the chromosome. Phylogenetic analysis revealed that banana GLRs could be divided into 3 subclades. The results of promoter cis-acting elements and transcription factor binding site prediction showed that there were multiple hormone- and stress-related response elements and 18 TFBS in banana GLR. RT-qPCR analysis showed that MaGLR1.1 and MaGLR3.5 responded positively to low temperature stress and were significantly expressed in abscisic acid/methyl jasmonate treatments. In conclusion, the results of this study suggest that GLR, a highly conserved family of ion channels, may play an important role in the growth and development process and stress resistance of banana.

Genome-wide identification of late embryogenesis abundant protein family and their key regulatory network in Pinus tabuliformis cold acclimation.

Cold acclimation is a crucial biological process that enables conifers to overwinter safely. The late embryogenesis abundant (LEA) protein family plays a pivotal role in enhancing freezing tolerance during this process. Despite its importance, the identification, molecular functions and regulatory networks of the LEA protein family have not been extensively studied in conifers or gymnosperms. Pinus tabuliformis, a conifer with high ecological and economic values and with high-quality genome sequence, is an ideal candidate for such studies. Here, a total of 104 LEA genes were identified from P. tabuliformis, and we renamed them according to their subfamily group: PtLEA1-PtLEA92 (group LEA1-LEA6), PtSMP1-PtSMP6 (group seed maturation protein) and PtDHN1-PtDHN6 (group Dehydrin). While the sequence structure of P. tabuliformis  LEA genes are conserved, their physicochemical properties exhibit unique characteristics within different subfamily groupings. Notably, the abundance of low-temperature responsive elements in PtLEA genes was observed. Using annual rhythm and temperature gradient transcriptome data, PtLEA22 was identified as a key gene that responds to low-temperature induction while conforming to the annual cycle of cold acclimation. Overexpression of PtLEA22 enhanced Arabidopsis freezing tolerance. Furthermore, several transcription factors potentially co-expressed with PtLEA22 were validated using yeast one-hybrid and dual-luciferase assays, revealing that PtDREB1 could directly bind PtLEA22 promoter to positively regulate its expression. These findings reveal the genome-wide characterization of P. tabuliformis  LEA genes and their importance in the cold acclimation, while providing a theoretical basis for studying the molecular mechanisms of cold acclimation in conifers.

Genome-wide identification and analysis of MIKC-type MADS-box genes expression in Chimonanthus salicifolius.

BACKGROUND: MIKC type MADS-box transcription factors are one of the largest gene families and play a pivotal role in flowering time and flower development. Chimonanthus salicifolius belongs to the family Calycanthaceae and has a unique flowering time and flowering morphology compared to other Chimonanthus species, but the research on MIKC type MADS-box gene family of C. salicifolius has not been reported. OBJECTIVE: Identification, comprehensive bioinformatic analysis, the expression pattern of MIKC-type MADS-box gene family from different tissues of C. salicifolius. METHODS: Genome-wide investigation and expression pattern under different tissues of the MIKC-type MADS-box gene family in C. salicifolius, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, promoter cis-acting element were performed. RESULTS: A total of 29 MIKC-type MADS-box genes were identified from the whole genome sequencing. Interspecies synteny analysis revealed more significant collinearity between C. salicifolius and the magnoliids species compared to eudicots and monocots. MIKC-type MADS-box genes from the same subfamily share similar distribution patterns, gene structure, and expression patterns. Compared with Arabidopsis thaliana, Nymphaea colorata, and Chimonanthus praecox, the FLC genes were absent in C. salicifolius, while the AGL6 subfamily was expanded in C. salicifolius. The selectively expanded promoter (AGL6) and lack of repressor (FLC) genes may explain the earlier flowering in C. salicifolius. The loss of the AP3 homologous gene in C. salicifolius is probably the primary cause of the morphological distinction between C. salicifolius and C. praecox. The csAGL6a gene is specifically expressed in the flowering process and indicates the potential function of promoting flowering. CONCLUSION: This study offers a genome-wide identification and expression profiling of the MIKC-types MADS-box genes in the C. salicifolius, and establishes the foundation for screening flowering development genes and understanding the potential function of the MIKC-types MADS-box genes in the C. salicifolius.

Genome-wide in silico identification of glutathione S-transferase (GST) gene family members in fig (Ficus carica L.) and expression characteristics during fruit color development.

Glutathione S-transferase (GSTs), a large and diverse group of multi-functional enzymes (EC 2.5.1.18), are associated with cellular detoxification, various biotic and abiotic stress responses, as well as secondary metabolites transportation. Here, 53 members of the FcGST gene family were screened from the genome database of fig (Ficus carica), which were further classified into five subfamilies, and the tau and phi were the major subfamilies. These genes were unevenly distributed over all the 13 chromosomes, and 12 tandem and one segmental duplication may contribute to this family expansion. Syntenic analysis revealed that FcGST shared closer genetic evolutionary origin relationship with species from the Ficus genus of the Moraceae family, such as F. microcarpa and F. hispida. The FcGST members of the same subfamily shared similar gene structure and motif distribution. The α helices were the chief structure element in predicted secondary and tertiary structure of FcGSTs proteins. GO and KEGG indicated that FcGSTs play multiple roles in glutathione metabolism and stress reactions as well as flavonoid metabolism. Predictive promoter analysis indicated that FcGSTs gene may be responsive to light, hormone, stress stimulation, development signaling, and regulated by MYB or WRKY. RNA-seq analysis showed that several FcGSTs that mainly expressed in the female flower tissue and peel during 'Purple-Peel' fig fruit development. Compared with 'Green Peel', FcGSTF1, and FcGSTU5/6/7 exhibited high expression abundance in the mature fruit purple peel. Additionally, results of phylogenetic sequences analysis, multiple sequences alignment, and anthocyanin content together showed that the expression changes of FcGSTF1, and FcGSTU5/6/7 may play crucial roles in fruit peel color alteration during fruit ripening. Our study provides a comprehensive overview of the GST gene family in fig, thus facilitating the further clarification of the molecular function and breeding utilization.

Genome-wide identification of the banana GLR gene family and its expression analysis in response to low temperature and abscisic acid/methyl jasmonate / 生物工程学报

Glutamate receptor-like (GLR) is an important class of Ca2+ channel proteins, playing important roles in plant growth and development as well as in response to biotic and abiotic stresses. In this paper, we performed genome-wide identification of banana GLR gene family based on banana genomic data. Moreover, we analyzed the basic physicochemical properties, gene structure, conserved motifs, promoter cis-acting elements, evolutionary relationships, and used real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) to verify the expression patterns of some GLR family members under low temperature of 4 ℃ and different hormone treatments. The results showed that there were 19 MaGLR family members in Musa acuminata, 16 MbGLR family members in Musa balbisiana and 14 MiGLR family members in Musa itinerans. Most of the members were stable proteins and had signal peptides, all of them had 3-6 transmembrane structures. Prediction of subcellular localization indicated that all of them were localized on the plasma membrane and irregularly distributed on the chromosome. Phylogenetic analysis revealed that banana GLRs could be divided into 3 subclades. The results of promoter cis-acting elements and transcription factor binding site prediction showed that there were multiple hormone- and stress-related response elements and 18 TFBS in banana GLR. RT-qPCR analysis showed that MaGLR1.1 and MaGLR3.5 responded positively to low temperature stress and were significantly expressed in abscisic acid/methyl jasmonate treatments. In conclusion, the results of this study suggest that GLR, a highly conserved family of ion channels, may play an important role in the growth and development process and stress resistance of banana.

Genome-wide identification of bZIP family genes and screening of candidate AarbZIPs involved in terpenoid biosynthesis in Artemisia argyi / 中国中药杂志

Artemisia argyi is an important medicinal and economic plant in China, with the effects of warming channels, dispersing cold, and relieving pain, inflammation, and allergy. The essential oil of this plant is rich in volatile terpenoids and widely used in moxi-bustion and healthcare products, with huge market potential. The bZIP transcription factors compose a large family in plants and are involved in the regulation of plant growth and development, stress response, and biosynthesis of secondary metabolites such as terpenoids. However, little is known about the bZIPs and their roles in A. argyi. In this study, the bZIP transcription factors in the genome of A. argyi were systematically identified, and their physicochemical properties, phylogenetic relationship, conserved motifs, and promoter-binding elements were analyzed. Candidate AarbZIP genes involved in terpenoid biosynthesis were screened out. The results showed that a total of 156 AarbZIP transcription factors were identified at the genomic level, with the lengths of 99-618 aa, the molecular weights of 11.7-67.8 kDa, and the theoretical isoelectric points of 4.56-10.16. According to the classification of bZIPs in Arabidopsis thaliana, the 156 AarbZIPs were classified into 12 subfamilies, and the members in the same subfamily had similar conserved motifs. The cis-acting elements of promoters showed that AarbZIP genes were possibly involved in light and hormonal pathways. Five AarbZIP genes that may be involved in the regulation of terpenoid biosynthesis were screened out by homologous alignment and phylogenetic analysis. The qRT-PCR results showed that the expression levels of the five AarbZIP genes varied significantly in different tissues of A. argyi. Specifically, AarbZIP29 and AarbZIP55 were highly expressed in the leaves and AarbZIP81, AarbZIP130, and AarbZIP150 in the flower buds. This study lays a foundation for the functional study of bZIP genes and their regulatory roles in the terpenoid biosynthesis in A. argyi.

[Genome-wide identification and effect of MdPEPC family genes during axillary bud outgrowth in apple (Malus domestica Borkh.)].

The PEPC family proteins are ubiquitous in various plants and play an important role in the process of photosynthetic carbon assimilation and have many non-photosynthetic biological functions. However, PEPC genes have not been reported in apple. In this study, the members of apple MdPEPC family were identified based on the new apple genome data by bioinformatics analysis, and their expression patterns in different tissues and the apple axillary bud transcriptome treated by decapitation and TDZ (cytokinin) were analyzed in order to explore the role of MdPEPC genes in apple axillary bud outgrowth. The results showed that 6 MdPEPC family members were identified in apple, which distributed on 6 different chromosomes, and had similar physicochemical characteristics. Phylogenetic tree and sequence alignment analysis showed that the MdPEPC could be divided into two subgroups (Group Ⅰ and Group Ⅱ), in which four members in MdPEPC family were clustered into Group Ⅰ, belonging to plant-type PEPCs. However, MdPEPC4 and MdPEPC5 were clustered into Group Ⅱ with AtPPC4, belonging to bacterial-type PEPCs. There were 7 pairs of fragments repeats among MdPEPC members, but no tandem repeats existed. The promoter cis-acting element analysis showed that MdPEPC genes were not only affected by light and stress, but also regulated by multiple hormones. The expression profiles showed that all MdPEPCs except MdPEPC4 and MdPEPC5 were expressed in different apple tissues. Transcriptome data analysis showed that the expression levels of MdPEPC1 and MdPEPC3 were up-regulated after decapitation and TDZ treatment, whereas MdPEPC2 was significantly down-regulated at 48 h after treatments. In conclusion, MdPEPC1, MdPEPC2 and MdPEPC3 were selected as the candidate genes involved in axillary bud outgrowth regulation for further study.

Identification and expression analysis of chlorophyll a/b binding protein gene family in grape (Vitis vinifera).

In higher plants, light capture of chlorophyll a/b binding protein (Lhc) plays a crucial role in the plant's response to adverse environment. So far, the family has not been systematically identified in grapes. In this study, 20 VvLhcs were identified in the grape genome, which were distributed in 13 of 19 grape chromosomes and divided into 7 developing branches. The results of gene duplication analysis showed that 6 VvLhcs formed fragment duplication events, while there was no tandem duplication in VvLhcs. Exon-intron structure analysis showed that they had a wide number of exons. Protein conserved motif analysis showed that VvLhcs contained more similar motif structures in the same phylogenetic branch. The cis-acting elements in the VvLhcs promoter region mainly respond to light, plant hormones and abiotic stresses. In addition, qRT-PCR results showed that different proportions of salt stress and red-blue light affected the expression of VvLhcs and the expression patterns of genes in different grape varieties were different. The results for further study on different grape varieties in different combinations of red and blue light of the Lhc provide a theoretical basis. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01204-5.

A Blueprint of Microstructures and Stage-Specific Transcriptome Dynamics of Cuticle Formation in Bombyx mori.

Insect cuticle is critical for the environmental adaptability and insecticide resistance of insects. However, there is no clear understanding of the structure and protein components of the cuticle during each developmental stage of holometabolous insects, and knowledge about the protein components within each layer is vague. We conducted serial sectioning, cuticular structure analysis, and transcriptome sequencing of the larval, pupal, and adult cuticles of Bombyx mori. The deposition processes of epicuticle, exocuticle, and endocuticle during larval, pupal, and adult cuticle formation were similar. Transcriptome analysis showed that these cuticle formations share 74% of the expressed cuticular protein (CP) genes and 20 other structural protein genes, such as larval serum protein and prisilkin. There are seven, six, and eleven stage-specific expressed CP genes in larval, pupal, and adult cuticles, respectively. The types and levels of CP genes may be the key determinants of the properties of each cuticular layer. For example, the CPs of the RR-2 protein family with high contents of histidine (His) are more essential for the exocuticle. Functional analysis suggested that BmorCPAP1-H is involved in cuticle formation. This study not only offers an in-depth understanding of cuticle morphology and protein components but also facilitates the elucidation of molecular mechanisms underlying cuticle formation in future studies.

DUF1005 Family Identification, Evolution Analysis in Plants, and Primary Root Elongation Regulation of CiDUF1005 From Caragana intermedia.

Proteins with a domain of unknown function (DUF) represent a number of gene families that encode functionally uncharacterized proteins in eukaryotes. In particular, members of the DUF1005 family in plants have a 411-amino-acid conserved domain, and this family has not been described previously. In this study, a total of 302 high-confidence DUF1005 family members were identified from 58 plant species, and none were found in the four algae that were selected. Thus, this result showed that DUF1005s might belong to a kind of plant-specific gene family, and this family has not been evolutionarily expanded. Phylogenetic analysis showed that the DUF1005 family genes could be classified into four subgroups in 58 plant species. The earliest group to emerge was Group I, including a total of 100 gene sequences, and this group was present in almost all selected species spanning from mosses to seed plants. Group II and Group III, with 69 and 74 members, respectively, belong to angiosperms. Finally, with 59 members, Group IV was the last batch of genes to emerge, and this group is unique to dicotyledons. Expression pattern analysis of the CiDUF1005, a member of the DUF1005 family from Caragana intermedia, showed that CiDUF1005 genes were differentially regulated under various treatments. Compared to the wild type, transgenic lines with heterologous CiDUF1005 expression in Arabidopsis thaliana had longer primary roots and more lateral roots. These results expanded our knowledge of the evolution of the DUF1005 family in plants and will contribute to elucidating biological functions of the DUF1005 family in the future.

A trouble shared is a trouble halved: The role of family identification and identification with humankind in well-being during the COVID-19 pandemic.

The COVID-19 pandemic has triggered health-related anxiety in ways that undermine peoples' mental and physical health. Contextual factors such as living in a high-risk area might further increase the risk of health deterioration. Based on the Social Identity Approach, we argue that social identities can not only be local that are characterized by social interactions, but also be global that are characterized by a symbolic sense of togetherness and that both of these can be a basis for health. In line with these ideas, we tested how identification with one's family and with humankind relates to stress and physical symptoms while experiencing health-related anxiety and being exposed to contextual risk factors. We tested our assumptions in a representative sample (N = 974) two-wave survey study with a 4-week time lag. The results show that anxiety at Time 1 was positively related to stress and physical symptoms at Time 2. Feeling exposed to risk factors related to lower physical health, but was unrelated to stress. Family identification and identification with humankind were both negatively associated with subsequent stress and family identification was negatively associated with subsequent physical symptoms. These findings suggest that for social identities to be beneficial for mental health, they can be embodied as well as symbolic.

Genome-Wide Identification and Expression Analysis of GA2ox, GA3ox, and GA20ox Are Related to Gibberellin Oxidase Genes in Grape (VitisVinifera L.).

Gibberellin (GAs) plays the important role in the regulation of grape developmental and growth processes. The bioinformatics analysis confirmed the differential expression of GA2, GA3, and GA20 gibberellin oxidase genes (VvGA2oxs, VvGA3oxs, and VvGA20oxs) in the grape genome, and laid a theoretical basis for exploring its role in grape. Based on the Arabidopsis GA2oxs, GA3oxs, and GA20oxs genes already reported, the VvGA2oxs, VvGA3oxs, and VvGA20oxs genes in the grape genome were identified using the BLAST software in the grape genome database. Bioinformatics analysis was performed using software such as DNAMAN v.5.0, Clustalx, MapGene2Chrom, MEME, GSDS v.2.0, ExPASy, DNAsp v.5.0, and MEGA v.7.0. Chip expression profiles were generated using grape Affymetrix GeneChip 16K and Grape eFP Browser gene chip data in PLEXdb. The expression of VvGA2oxs, VvGA3oxs, and VvGA20oxs gene families in stress was examined by qRT-PCR (Quantitative real-time-PCR). There are 24 GAoxs genes identified with the grape genome that can be classified into seven subgroups based on a phylogenetic tree, gene structures, and conserved Motifs in our research. The gene family has higher codon preference, while selectivity is negative selection of codon bias and selective stress was analyzed. The expression profiles indicated that the most of VvGAox genes were highly expressed under different time lengths of ABA (Abscisic Acid) treatment, NaCl, PEG and 5 °C. Tissue expression analysis showed that the expression levels of VvGA2oxs and VvGA20oxs in different tissues at different developmental stages of grapes were relatively higher than that of VvGA3oxs. Last but not least, qRT-PCR (Real-time fluorescent quantitative PCR) was used to determine the relative expression of the GAoxs gene family under the treatment of GA3 (gibberellin 3) and uniconazole, which can find that some VvGA2oxs was upregulated under GA3 treatment. Simultaneously, some VvGA3oxs and VvGA20oxs were upregulated under uniconazole treatment. In a nutshell, the GA2ox gene mainly functions to inactivate biologically active GAs, while GA20ox mainly degrades C20 gibberellins, and GA3ox is mainly composed of biologically active GAs. The comprehensive analysis of the three classes of VvGAoxs would provide a basis for understanding the evolution and function of the VvGAox gene family in a grape plant.

HIPPI: highly accurate protein family classification with ensembles of HMMs.

BACKGROUND: Given a new biological sequence, detecting membership in a known family is a basic step in many bioinformatics analyses, with applications to protein structure and function prediction and metagenomic taxon identification and abundance profiling, among others. Yet family identification of sequences that are distantly related to sequences in public databases or that are fragmentary remains one of the more difficult analytical problems in bioinformatics. RESULTS: We present a new technique for family identification called HIPPI (Hierarchical Profile Hidden Markov Models for Protein family Identification). HIPPI uses a novel technique to represent a multiple sequence alignment for a given protein family or superfamily by an ensemble of profile hidden Markov models computed using HMMER. An evaluation of HIPPI on the Pfam database shows that HIPPI has better overall precision and recall than blastp, HMMER, and pipelines based on HHsearch, and maintains good accuracy even for fragmentary query sequences and for protein families with low average pairwise sequence identity, both conditions where other methods degrade in accuracy. CONCLUSION: HIPPI provides accurate protein family identification and is robust to difficult model conditions. Our results, combined with observations from previous studies, show that ensembles of profile Hidden Markov models can better represent multiple sequence alignments than a single profile Hidden Markov model, and thus can improve downstream analyses for various bioinformatic tasks. Further research is needed to determine the best practices for building the ensemble of profile Hidden Markov models. HIPPI is available on GitHub at https://github.com/smirarab/sepp .

Family identification: a beneficial process for young adults who grow up in homes affected by parental intimate partner violence.

Exposure to parental intimate partner violence (parental IPV) is a complex trauma. Research within social psychology establishes that identification with social groups impacts positively on how we appraise, respond to and recover from traumatic events. IPV is also a highly stigmatized social phenomenon and social isolation is a major factor for families affected by IPV, yet strong identification with the family group may act as a beneficial psychological resource to young people who grew up in homes affected by IPV. The current study, an online survey of 355 students (M age = 20, 70% female), investigated if a psychosocial process, specifically identification with the family, may influence the relationship between the predictor, exposure to parental IPV, and outcomes, global self-esteem and state anxiety. Mediation analysis suggests that identification with the family has a positive influence on the relationship between exposure to parental IPV and psychological outcomes; exposure to parental IPV results in reduced family identification, but when family identification is strong it results in both reduced anxiety and increased self-esteem for young people. The findings highlight the importance of having a strong sense of belonging to the extended family for young people who were exposed to parental IPV, thus has implications for prevention, intervention, and social policy.