Distribution Of Motifs Research Articles

Interplay of intracellular and trans-cellular DNA methylation in natural archaeal consortia.

DNA methylation serves a variety of functions across all life domains. In this study, we investigated archaeal methylomics within a tripartite xylanolytic halophilic consortium. This consortium includes Haloferax lucertense SVX82, Halorhabdus sp. SVX81, and an ectosymbiotic Candidatus Nanohalococcus occultus SVXNc, a nano-sized archaeon from the DPANN superphylum. We utilized PacBio SMRT and Illumina cDNA sequencing to analyse samples from consortia of different compositions for methylomics and transcriptomics. Endogenous cTAG methylation, typical of Haloferax, was accompanied in this strain by methylation at four other motifs, including GDGcHC methylation, which is specific to the ectosymbiont. Our analysis of the distribution of methylated and unmethylated motifs suggests that autochthonous cTAG methylation may influence gene regulation. The frequency of GRAGAaG methylation increased in highly expressed genes, while CcTTG and GTCGaGG methylation could be linked to restriction-modification (RM) activity. Generally, the RM activity might have been reduced during the evolution of this archaeon to balance the protection of cells from intruders, the reduction of DNA damage due to self-restriction in stressful environments, and the benefits of DNA exchange under extreme conditions. Our methylomics, transcriptomics and complementary electron cryotomography (cryo-ET) data suggest that the nanohaloarchaeon exports its methyltransferase to methylate the Haloferax genome, unveiling a new aspect of the interaction between the symbiont and its host.

Motif-Backdoor: Rethinking the Backdoor Attack on Graph Neural Networks via Motifs

Graph neural network (GNN) with a powerful representation capability has been widely applied to various areas. Recent works have exposed that GNN is vulnerable to the backdoor attack, i.e., models trained with maliciously crafted training samples are easily fooled by patched samples. Most of the proposed studies launch the backdoor attack using a trigger that is either the randomly generated subgraph e.g., erdős-rényi backdoor (ER-B) for less computational burden or the gradient-based generative subgraph e.g., graph trojaning attack (GTA) to enable a more effective attack. However, the interpretation of how is the trigger structure and the effect of the backdoor attack related has been overlooked in the current literature. Motifs, recurrent and statistically significant subgraphs in graphs, contain rich structure information. In this article, we are rethinking the trigger from the perspective of motifs and propose a motif-based backdoor attack, denoted as <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Motif-Backdoor</i> . It contributes from three aspects. First, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Interpretation</i> , it provides an in-depth explanation for backdoor effectiveness by the validity of the trigger structure from motifs, leading to some novel insights, e.g., using subgraphs that appear less frequently in the graph as the trigger can achieve better attack performance. Second, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Effectiveness</i> , Motif-Backdoor reaches the state-of-the-art (SOTA) attack performance in both black-box and defensive scenarios. Third, <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">Efficiency</i> , based on the graph motif distribution, Motif-Backdoor can quickly obtain an effective trigger structure without target model feedback or subgraph model generation. Extensive experimental results show that Motif-Backdoor realizes the SOTA performance on three popular models and four public datasets compared with five baselines, e.g., Motif-Backdoor improves the attack success rate (ASR) by 14.73% compared with baselines on average. In addition, under a possible defense, Motif-Backdoor still implements satisfying performance, highlighting the requirement of defenses against backdoor attacks on GNNs. The datasets and code are available at https://github.com/Seaocn/Motif-Backdoor.

Genome-wide identification and characterization of DIRIGENT gene family (CcDIR) in pigeonpea (Cajanus cajan L.) provide insights on their spatial expression pattern and relevance to stress response

This study is a thorough characterization of pigeonpea dirigent gene (CcDIR) family, an important component of the lignin biosynthesis pathway. Genome-wide analysis identified 25 CcDIR genes followed by a range of analytical approaches employed to unravel their structural and functional characteristics. Structural examination revealed a classic single exon and no intron arrangement in CcDIRs contributing to our understanding on evolutionary dynamics. Phylogenetic analysis elucidated evolutionary relationships among CcDIR genes with six DIR sub-families, while motif distribution analysis displayed and highlighted ten conserved protein motifs in CcDIRs. Promoter analyses of all the dirigent genes detected 18 stress responsive cis-acting elements offering insights into transcriptional regulation. While spatial expression analyses across six plant tissues showed preferential expression of CcDIR genes, exposure to salt (CcDIR2 and CcDIR9) and herbivory (CcDIR1, CcDIR2, CcDIR3 and CcDIR11), demonstrated potential roles of specific DIRs in plant defense. Interestingly, increased gene expression during herbivory, also correlated with increased lignin content authenticating the specific response. Furthermore, exogenous application of stress hormones, SA and MeJA on leaves significantly induced the expression of CcDIRs that responded to herbivory. Taken together, these findings contribute to a comprehensive understanding of CcDIR genes impacting development and stress response in the important legume pigeonpea.

Genome-wide identification and functional characterization of trans-IDS genes give new insights into terpene biosynthesis in Syringa oblata

Floral scents play a vital role in plant pollination and ornamental plant industry. Syringa oblata, a typical aromatic and medicinal plant, holds a high terpene concentration as the main fragrant volatiles and plays a crucial role in the breeding in the Syringa genus. The trans-isopentenyl diphosphate synthases (trans-IDSs) were critical enzymes in terpenoid diversity and accumulation. However, genome-wide characterization and systematic analysis of the lilac trans-IDS family are still under-explored. In this study, we identified sixteen trans-IDS genes belonging to five subfamilies on 10 chromosomes of the lilac genome. Each subfamily shared similar exon, intron, and motif distributions. Synteny analysis suggested that the lilac trans-IDS gene family expanded through segmental and tandem duplication events. We also analyzed their expression levels in various organs and flower development stages for terpene biosynthesis functional insights. Furthermore, spatiotemporal gene expression, virus-induced gene silencing (VIGS), transient overexpression, and transgenic Arabidopsis results collectively showed SoGPPS2 enhances terpene volatile synthesis. SoGPPS2 also impacted plant growth and flowering. And potential upstream transcription regulatory networks for SoGPPS2 were identified. In summary, this study enhances our comprehensive understanding of the lilac trans-IDS gene family for the first time, providing crucial insights into terpene synthesis and laying the groundwork for future floral aromatic breeding, medicinal and essential oil industries in lilac.

Abstract 1695: Genomic clustering tendency of transcription factors reflects phase-separated transcriptional condensates at cancer super-enhancers

Abstract Many transcription factors (TFs) have been shown to bind to super-enhancers, forming transcriptional condensates to activate transcription in various cellular systems including cancer. However, the genomic and epigenomic determinants of phase-separated transcriptional condensate formation remain poorly understood. Questions regarding which TFs tend to associate with transcriptional condensates and what factors influence their association are largely unanswered. Here we systematically analyzed 528 DNA sequence motifs across the human genome and 6,650 TF binding profiles across different cell types to identify the molecular features contributing to the formation of transcriptional condensates. We found that the genomic distributions of sequence motifs for different TFs exhibit distinct clustering tendencies. Notably, TF motifs with a high genomic clustering tendency are significantly associated with super-enhancers. TF binding profiles showing a high genomic clustering tendency are further enriched at cell-type-specific super-enhancers. TFs with a high binding clustering tendency also possess high liquid-liquid phase separation abilities. Compared to non-clustered TF binding, densely clustered TF binding sites are more enriched at cell-type-specific super-enhancers with higher chromatin accessibility, elevated chromatin interaction, and stronger association with cancer outcomes. Our results indicate that the clustered genomic binding patterns and the phase separation properties of TFs collectively contribute to the formation of transcriptional condensates at cancer super-enhancers. Citation Format: Zhenjia Wang, Shengyuan Wang, Chongzhi Zang. Genomic clustering tendency of transcription factors reflects phase-separated transcriptional condensates at cancer super-enhancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1695.

Genome-wide identification, characterization and evolutionary analysis of betaine aldehyde dehydrogenase (BADH), mitogen-activated protein kinase (MAPK) and sodium/hydrogen exchanger (NHX) genes in maize (Zea mays) under salt stress

AbstractBetaine aldehyde dehydrogenase (BADH), mitogen-activated protein kinase (MAPK) and sodium/hydrogen exchanger (NHX) play important roles in the response to salt stress. This is the first study to identify the BADH and NHX genes in maize (Zea mays) via genome-wide analysis. The qRT‒PCR results indicated that ZmNHX was upregulated by 4.38-fold, while a significant difference was not observed in ZmBADH or ZmMAPK, with fold changes of 0.96 and 1.06, respectively, under salinity stress. Genome-wide analysis revealed 8 ZmBADH, 19 ZmMAPK and 11 ZmNHX proteins in Z. mays. Domain analysis confirmed the presence of the aldehyde dehydrogenase superfamily (ALDH-SF), protein kinase and Na_H_Exchanger domains in the ZmBADH, ZmMAPK and ZmNHX proteins, respectively. Motif analysis indicated that the phylogenetic relationships were similar to the conserved motif distributions within the clade. The Ka/Ks ratio indicated that the ZmBADH, ZmMAPK and ZmNHX genes were influenced primarily by purifying selection. This study provides comprehensive identification, characterization, and evolutionary analysis for a better understanding of the ZmBADH, ZmMAPK and ZmNHX genes in maize.

Genome-wide identification and expression analysis of NF-Y gene family in tobacco (Nicotiana tabacum L.)

Nuclear factor Y (NF-Y) gene family is an important transcription factor composed of three subfamilies of NF-YA, NF-YB and NF-YC, which is involved in plant growth, development and stress response. In this study, 63 tobacco NF-Y genes (NtNF-Ys) were identified in Nicotiana tabacum L., including 17 NtNF-YAs, 30 NtNF-YBs and 16 NtNF-YCs. Phylogenetic analysis revealed ten pairs of orthologues from tomato and tobacco and 25 pairs of paralogues from tobacco. The gene structure of NtNF-YAs exhibited similarities, whereas the gene structure of NtNF-YBs and NtNF-YCs displayed significant differences. The NtNF-Ys of the same subfamily exhibited a consistent distribution of motifs and protein 3D structure. The protein interaction network revealed that NtNF-YC12 and NtNF-YC5 exhibited the highest connectivity. Many cis-acting elements related to light, stress and hormone response were found in the promoter of NtNF-Ys. Transcriptome analysis showed that more than half of the NtNF-Y genes were expressed in all tissues, and NtNF-YB9/B14/B15/B16/B17/B29 were specifically expressed in roots. A total of 15, 12, 5, and 6 NtNF-Y genes were found to respond to cold, drought, salt, and alkali stresses, respectively. The results of this study will lay a foundation for further study of NF-Y genes in tobacco and other Solanaceae plants.

Genome-Wide Analysis of the TCP Transcription Factor Gene Family in Pepper (Capsicum annuum L.).

TCP transcription factors play a key role in regulating various developmental processes, particularly in shoot branching, flower development, and leaf development, and these factors are exclusively found in plants. However, comprehensive studies investigating TCP transcription factors in pepper (Capsicum annuum L.) are lacking. In this study, we identified 27 CaTCP members in the pepper genome, which were classified into Class I and Class II through phylogenetic analysis. The motif analysis revealed that CaTCPs in the same class exhibit similar numbers and distributions of motifs. We predicted that 37 previously reported miRNAs target 19 CaTCPs. The expression levels of CaTCPs varied in various tissues and growth stages. Specifically, CaTCP16, a member of Class II (CIN), exhibited significantly high expression in flowers. Class I CaTCPs exhibited high expression levels in leaves, while Class II CaTCPs showed high expression in lateral branches, especially in the CYC/TB1 subclass. The expression profile suggests that CaTCPs play specific roles in the developmental processes of pepper. We provide a theoretical basis that will assist in further functional validation of the CaTCPs.

Genome-Wide Identification and Expression Analysis Unveil the Involvement of the Cold Shock Protein (CSP) Gene Family in Cotton Hypothermia Stress.

Cold shock proteins (CSPs) are DNA/RNA binding proteins with crucial regulatory roles in plant growth, development, and stress responses. In this study, we employed bioinformatics tools to identify and analyze the physicochemical properties, conserved domains, gene structure, phylogenetic relationships, cis-acting elements, subcellular localization, and expression patterns of the cotton CSP gene family. A total of 62 CSP proteins were identified across four cotton varieties (Gossypium arboreum, Gossypium raimondii, Gossypium barbadense, Gossypium hirsutum) and five plant varieties (Arabidopsis thaliana, Brassica chinensis, Camellia sinensis, Triticum aestivum, and Oryza sativa). Phylogenetic analysis categorized cotton CSP proteins into three evolutionary branches, revealing similar gene structures and motif distributions within each branch. Analysis of gene structural domains highlighted the conserved CSD and CCHC domains across all cotton CSP families. Subcellular localization predictions indicated predominant nuclear localization for CSPs. Examination of cis-elements in gene promoters revealed a variety of elements responsive to growth, development, light response, hormones, and abiotic stresses, suggesting the potential regulation of the cotton CSP family by different hormones and their involvement in diverse stress responses. RT-qPCR results suggested that GhCSP.A1, GhCSP.A2, GhCSP.A3, and GhCSP.A7 may play roles in cotton's response to low-temperature stress. In conclusion, our findings underscore the significant role of the CSP gene family in cotton's response to low-temperature stress, providing a foundational basis for further investigations into the functional aspects and molecular mechanisms of cotton's response to low temperatures.

Genome-wide survey of peptides containing tyrosine sulfation (PSY) gene family and potential PSY specific miRNA revealed their role in plant development and diverse stress conditions in rice (Oryza sativa L.)

Signalling peptides play a vital role in regulating a wide range of aspects, including plant growth, development, and responses to various biotic and abiotic conditions. In recent times, small post-translationally altered peptides have garnered significant attention from plant biologists due to their crucial role as signalling molecules in various biological processes. Until now, the biological function of peptides containing Tyrosine Sulfation (PSY) gene family members has only been elucidated in Arabidopsis. However, their functions in plant developmental processes and their responses to different biotic and abiotic stress conditions remain elusive in other crops. Therefore, in this study, we identified 109 PSY genes in the rice genome for the first time. However, our detailed in silico analysis focused specifically on 7 OsPSY genes of the O. sativa Japonica Group (Nippobare). Phylogenetic analysis of PSY genes from various plant species has classified these genes into eight distinct groups, designated as Group I to VIII. Interestingly, the 7 OsPSY genes identified in rice were found to be distributed across all four groups including I, V, VI and VIII. In addition, closely linked groups of PSY genes showed more similarities in terms of gene structures, conserved motif distributions, and corresponding subcellular locations. Further, the Ka/Ks score analysis revealed that OsPSY genes have undergone significant purifying selection throughout their evolutionary history. Furthermore, multiple CAREs were discovered in the OsPSY promoters, suggesting their potential importance in regulating plant growth, responses to plant hormones, and multifactorial stress conditions. Additionally, the potential miRNAs that specifically targeted to OsPSY genes were identified and also examined their expression kinetics in different tissues and developmental stage. Out of the seven OsPSY genes, six OsPSY genes were found to be targeted by twelve miRNAs. Moreover the expression profile of OsPSY genes showed a differential expression pattern in various tissues and multifactorial stress conditions. Therefore, these findings serve as a valuable reference for further molecular studies aimed at enhancing the yield and stress tolerance of rice. The understanding the functions and regulatory mechanisms of the OsPSY genes can potentially lead to the development of targeted approaches to improve crop performance under various environmental conditions, ultimately contributing to agricultural productivity and global food security.

Temperature-induced suppression of structural disproportionation in paramagnetic quantum materials

With the development of electronic structure theory, a new class of materials—quantum ones—has been recognized by the community. Traditionally, it has been believed that the properties of such compounds cannot be described within the framework of modern density functional theory, and indeed, more advanced post-mean-field theory methods are needed. Motivated by this, herein, we develop a fundamental understanding of such complex materials using the example of paramagnetic YNiO3, which is experimentally known to exhibit metal-to-insulator phase transition. We show that this material has a temperature-dependent distribution of local motifs. Thus, while at low temperatures, YNiO3 has distinct structural disproportionation with the formation of large and small octahedra, as the temperature increases, this disproportionation is suppressed. We also explain the paramagnetic monoclinic to paramagnetic orthorhombic phase transition within the double-well to single-well energy profile, predicting the variation in the corresponding energy profile as a function of octahedral size distribution. In this way, we demonstrate a fundamental understanding of structural phase transitions in quantum materials, giving insights into how they can be used for different applications and what minimum level of theory is needed to describe such types of complex materials at finite temperatures.

Genome-Wide Characterization and Phylogenetic and Stress Response Expression Analysis of the MADS-Box Gene Family in Litchi (Litchi chinensis Sonn.).

The MADS-box protein is an important transcription factor in plants and plays an important role in regulating the plant abiotic stress response. In this study, a total of 94 MADS-box genes were predicted in the litchi genome, and these genes were widely distributed on all the chromosomes. The LcMADS-box gene family was divided into six subgroups (Mα, Mβ, Mγ, Mδ, MIKC, and UN) based on their phylogenetical relationships with Arabidopsis, and the closely linked subgroups exhibited more similarity in terms of motif distribution and intron/exon numbers. Transcriptome analysis indicated that LcMADS-box gene expression varied in different tissues, which can be divided into universal expression and specific expression. Furthermore, we further validated that LcMADS-box genes can exhibit different responses to various stresses using quantitative real-time PCR (qRT-PCR). Moreover, physicochemical properties, subcellular localization, collinearity, and cis-acting elements were also analyzed. The findings of this study provide valuable insights into the MADS-box gene family in litchi, specifically in relation to stress response. The identification of hormone-related and stress-responsive cis-acting elements in the MADS-box gene promoters suggests their involvement in stress signaling pathways. This study contributes to the understanding of stress tolerance mechanisms in litchi and highlights potential regulatory mechanisms underlying stress responses.

Directed Network Comparison Using Motifs.

Analyzing and characterizing the differences between networks is a fundamental and challenging problem in network science. Most previous network comparison methods that rely on topological properties have been restricted to measuring differences between two undirected networks. However, many networks, such as biological networks, social networks, and transportation networks, exhibit inherent directionality and higher-order attributes that should not be ignored when comparing networks. Therefore, we propose a motif-based directed network comparison method that captures local, global, and higher-order differences between two directed networks. Specifically, we first construct a motif distribution vector for each node, which captures the information of a node's involvement in different directed motifs. Then, the dissimilarity between two directed networks is defined on the basis of a matrix, which is composed of the motif distribution vector of every node and the Jensen-Shannon divergence. The performance of our method is evaluated via the comparison of six real directed networks with their null models, as well as their perturbed networks based on edge perturbation. Our method is superior to the state-of-the-art baselines and is robust with different parameter settings.

Comprehensive analysis of WOX transcription factors provide insight into genes related to the regulation of unisexual flowers development in Akebia trifoliata

Akebia trifoliata is a fascinating economic and medicinal plant that produces functionally unisexual flowers due to stamen/pistil abortion during flower development, and the genetic regulation pathway of this process remain completely unknown. Here, 10 AktWOXs were identified for the first time, all contained a highly conserved homeodomain. AktWOXs were divided into three clades, each with the same or similar intron, exon, and motifs distribution. Many cis-elements related to stress response, growth and development, and hormone response were found in the AktWOXs promoter region. In addition, four candidate genes AktWOX8, AktWOX11, AktWOX13.2 and AktWUS that might be involved in unisexual flowers development were screened, all of which were located in the nucleus and showed transcriptional activation activity. Yeast one-hybrid showed that both AktKNU and AktAG1, the potential core transcription factors in the activity termination pathway of flower meristem stem cells, could bind to the promoter region of AktWUS. Dual-luciferase assay further confirmed that only AktKNU inhibited the expression of AktWUS. Collectively, this study revealed the mechanism of AktWUS that might affect the formation of unisexual flowers by regulating the timely termination of flower meristem in A. trifoliata.

Genome-wide identification, classification and expression analysis of MYB gene family in coconut (Cocos nucifera L.).

MYB transcription factors regulate the growth, development, and secondary metabolism of plant species. To investigate the origin of color variations in coconut pericarp, we identified and analyzed the MYB gene family present in coconut. According to the sequence of MYB genes in Arabidopsis thaliana, homologous MYB gene sequences were found in the whole genome database of coconut, the conserved sequence motifs within MYB proteins were analyzed by Motif Elicitation (MEME) tool, and the sequences without conservative structure were eliminated. Additionally, we employed RNA-seq technology to generate gene expression signatures of the R2R3-MYB genes across distinctive coconut parts exhibiting diverse colors. To validate these profiles, we conducted quantitative PCR (qPCR). Through comprehensive genome-wide screening, we successfully identified a collection of 179 MYB genes in coconut. Subsequent phylogenetic analysis categorized these 179 coconut MYB genes into 4-subfamilies: 124 R2R3-MYB, 4 3R-MYB types, 4 4R-MYB type, and 47 unknown types. Furthermore, these genes were further divided into 34 subgroups, with 28 of these subgroups successfully classified into known subfamilies found in Arabidopsis thaliana. By mapping the CnMYB genes onto the 16 chromosomes of the coconut genome, we unveiled a collinearity association between them. Moreover, a preservation of gene structure and motif distribution was observed across the CnMYB genes. Our research encompassed a thorough investigation of the R2R3-MYB genes present in the coconut genome, including the chromosomal localization, gene assembly, conserved regions, phylogenetic associations, and promoter cis-acting elements of the studied genes. Our findings revealed a collection of 12 R2R3-MYB candidate genes, namely CnMYB8, CnMYB15, CnMYB27, CnMYB28, CnMYB61, CnMYB63, CnMYB68, CnMYB94, CnMYB101, CnMYB150, CnMYB153, and CnMYB164. These genes showed differential expressions in diverse tissues and developmental stages of four coconut species, such as CnMYB68, CnMYB101, and CnMYB28 exhibited high expression in majority of tissues and coconut species, while CnMYB94 and CnMYB164 showed lower expression. These findings shed light on the crucial functional divergence of CnMYB genes across various coconut tissues, suggesting these genes as promising candidate genes for facilitating color development in this important crop.

Genome-wide identification and analysis of the GGCT gene family in wheat

Backgroundγ-glutamylcyclotransferase (GGCT), an enzyme to maintain glutathione homeostasis, plays a vital role in the response to plant growth and development as well as the adaptation to various stresses. Although the GGCT gene family analysis has been conducted in Arabidopsis and rice, the family genes have not yet been well identified and analyzed at the genome-wide level in wheat (Triticum aestivum L.).ResultsIn the present study, 20 TaGGCT genes were identified in the wheat genome and widely distributed on chromosomes 2A, 2B, 2D, 3A, 4A, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B, and 7D. Phylogenetic and structural analyses showed that these TaGGCT genes could be classified into three subfamilies: ChaC, GGGACT, and GGCT-PS. They exhibited similar motif compositions and distribution patterns in the same subgroup. Gene duplication analysis suggested that the expansion of TaGGCT family genes was facilitated by segmental duplications and tandem repeats in the wheat evolutionary events. Identification of diverse cis-acting response elements in TaGGCT promoters indicated their potential fundamental roles in response to plant development and abiotic stresses. The analysis of transcriptome data combined with RT-qPCR results revealed that the TaGGCTs genes exhibited ubiquitous expression across plant organs, with highly expressed in roots, stems, and developing grains. Most TaGGCT genes were up-regulated after 6 h under 20% PEG6000 and ABA treatments. Association analysis revealed that two haplotypes of TaGGCT20 gene displayed significantly different Thousand-kernel weight (TKW), Kernel length (KL), and Kernel width (KW) in wheat. The geographical and annual distribution of the two haplotypes of TaGGCT20 gene further revealed that the frequency of the favorable haplotype TaGGCT20-Hap-I was positively selected in the historical breeding process of wheat.ConclusionThis study investigated the genome-wide identification, structure, evolution, and expression analysis of TaGGCT genes in wheat. The motifs of TaGGCTs were highly conserved throughout the evolutionary history of wheat. Most TaGGCT genes were highly expressed in roots, stems, and developing grains, and involved in the response to drought stresses. Two haplotypes were developed in the TaGGCT20 gene, where TaGGCT20-Hap-I, as a favorable haplotype, was significantly associated with higher TKW, KL, and KW in wheat, suggesting that the haplotype is used as a function marker for the selection in grain yield in wheat breeding.

Systematic analysis of Heat Shock Protein 70 (HSP70) gene family in radish and potential roles in stress tolerance

The 70 kD heat shock proteins (HSP70s) represent a class of molecular chaperones that are widely distributed in all kingdoms of life, which play important biological roles in plant growth, development, and stress resistance. However, this family has not been systematically characterized in radish (Raphanus sativus L.). In this study, we identified 34 RsHSP70 genes unevenly distributed within nine chromosomes of R. sativus. Phylogenetic and multiple sequence alignment analyses classified the RsHSP70 proteins into six distinct groups (Group A–F). The characteristics of gene structures, motif distributions, and corresponding cellular compartments were more similar in closely linked groups. Duplication analysis revealed that segmental duplication was the major driving force for the expansion of RsHSP70s in radish, particularly in Group C. Synteny analysis identified eight paralogs (Rs-Rs) in the radish genome and 19 orthologs (Rs-At) between radish and Arabidopsis, and 23 orthologs (Rs-Br) between radish and Chinese cabbage. RNA-seq analysis showed that the expression change of some RsHSP70s were related to responses to heat, drought, cadmium, chilling, and salt stresses and Plasmodiophora brassicae infection, and the expression patterns of these RsHSP70s were significantly different among 14 tissues. Furthermore, we targeted a candidate gene, RsHSP70–23, the product of which is localized in the cytoplasm and involved in the responses to certain abiotic stresses and P. brassicae infection. These findings provide a reference for further molecular studies to improve yield and stress tolerance of radish.

Genome-wide identification, characterization, and expression analysis of m6A readers-YTH domain-containing genes in alfalfa

Eukaryotic messenger RNAs (mRNAs) are often modified with methyl groups at the N6 position of adenosine (m6A), and these changes are interpreted by YTH domain-containing proteins to regulate the metabolism of m6A-modified mRNAs. Although alfalfa (Medicago sativa) is an established model organism for forage development, the understanding of YTH proteins in alfalfa is still limited. In the present investigation, 53 putative YTH genes, each encoding a YT521 domain-containing protein, were identified within the alfalfa genome. These genes were categorized into two subfamilies: YTHDF (49 members) and YTHDC (four members). Each subfamily demonstrates analogous motif distributions and domain architectures. Specifically, proteins encoded by MsYTHDF genes incorporate a single domain structure, while those corresponding to MsYTH5, 8, 12, 16 who are identified as members of the MsYTHDC subfamily, exhibit CCCH-type zinc finger repeats at their N-termini. It is also observed that the predicted aromatic cage pocket that binds the m6A residue of MsYTHDC consists of a sequence of two tryptophan residues and one tyrosine residue (WWY). Conversely, in MsYTHDF, the binding pocket comprises two highly conserved tryptophan residues and either one tryptophan residue (WWW) or tyrosine residue (WWY) in MsYTHDF.Through comparative analysis of qRT-PCR data, we observed distinct expression patterns in specific genes under abiotic stress, indicating their potential regulatory roles. Notably, five genes (MsYTH2, 14, 26, 27, 48) consistently exhibit upregulation, and two genes (MsYTH33, 35) are downregulated in response to both cold and salt stress. This suggests a common mechanism among these YTH proteins in response to various abiotic stressors in alfalfa. Further, integrating qRT-PCR with RNA-seq data revealed that MsYTH2, MsYTH14, and MsYTH16 are highly expressed in leaves at various development stages, underscoring their potential roles in regulating the growth of these plant parts. The obtained findings shed further light on the biological functions of MsYTH genes and may aid in the selection of suitable candidate genes for future genetic enhancement endeavors aimed at improving salt and cold tolerance in alfalfa.

Genome-Wide Identification and Analysis of ZF-HD Gene Family in Moso Bamboo (Phyllostachys edulis)

Zinc finger-homeodomain (ZF-HD) proteins play essential roles in plant growth, development and stress responses. However, knowledge of the expression and evolutionary history of ZF-HD genes in moso bamboo remains limited. In this study, a total of 24 ZF-HD genes were found unevenly distributed on 12 chromosomes in moso bamboo (Phyllostachys edulis). Phylogenetic analysis indicated that PeZF-HDs were divided into two subfamilies: ZHD and MIF. The ZHD subfamily genes were further classified into seven groups according to their orthologous relationships among the rice and Arabidopsis ZF-HD gene family. The gene structures and conserved motifs of PeZF-HDs were analyzed. Whole-genome duplication (WGD) or segmental duplication promoted the evolution and expansion of the moso bamboo ZF-HD gene family. Ka/Ks ratios suggested that the twenty-four duplication pairs had undergone purifying selection. Promoter analysis showed that most PeZF-HDs contained cis-elements associated with stress responses and hormones. Expression analysis demonstrated that many PeZF-HDs were responsive to abiotic stress treatment. Overall, this work investigated PeZF-HD genes in moso bamboo using bioinformatic approaches. The evolutionary research on gene structure, motif distribution and cis-regulatory elements indicated that PeZF-HDs play distinct roles in biological processes, which provides a theoretical basis for exploring the physiological functions of ZF-HDs and selecting candidate stress-related genes in moso bamboo.

Identification and bioinformatics analysis of MADS-box family genes containing K-box domain in maize

The MADS-box family genes are involved in the development of plant roots, leaves, flowers, and fruits, and play a crucial role in plant growth and development. Studying MADS-box genes with K-box domain is crucial to distinguish different types of MADS-box genes. This study systematically analysed the genomic structural information of maize MADS-box family members containing the K-box Domain at the genome-wide level using the maize (Zea mays) B73 genome as the reference sequence, and provided insight into the biological functions of the maize MADS-box family containing the K-box domain. According to the findings, 52 MADS-box family genes with K-box domain were identified and divided into 4 subgroups. The distribution of motif in the same subgroup was found to be relatively conservative, and all of them had MADS-box conserved domain and K-box domain. Gene structure analysis showed that the introns and exons of the same subgroup genes have similar gene structure, and different types of genes containing the K-box domain showed different exon/intron structure characteristics. Chromosome mapping showed that 52 genes containing the K-box domain were unevenly distributed on the 10 chromosomes of maize, most of which were distributed at both ends of the chromosome and a small number of genes were distributed near the centromere. Based on the analysis of cis-acting elements of it up-stream promoter, it was found that MADS-box family genes may be involved in light response, IAA, GA, ABA, and LTR signal pathways, indicating that they play a certain role in stress response and hormone signal transduction. The expression analysis of genes with the K-box domain in maize leaves treated with auxin and gibberellin revealed that MADS-box genes may have a regulatory effect on certain plant hormones. Through the identification and bioinformatics analysis of MADS-box family genes containing the K-box domain, it is helpful to further study the function and pathway of MADS-box family genes, and provide a theoretical basis for further re-search for the molecular mechanism of maize growth and development.