<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2025.1652613</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Genome-wide identification of the <italic>COBRA-like</italic> gene family and expression pattern analysis under abiotic stresses of <italic>Sorghum bicolor</italic> (L.)</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Liu</surname>
<given-names>Shipeng</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3190186/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liang</surname>
<given-names>Shuang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3110711/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jing</surname>
<given-names>Tingrui</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3190186/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Xinyi</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3190185/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Hairuo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3190411/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Quan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Junshen</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Kai</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3109778/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>He</surname>
<given-names>Xiaolong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhao</surname>
<given-names>Haibing</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jiang</surname>
<given-names>Wenting</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Xiangqian</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Life Sciences, Yan&#x2019;an University</institution>, <addr-line>Yan&#x2019;an</addr-line>,&#xa0;<country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Engineering Research Center of Microbial Resources Development and Green Recycling, University of Shaanxi Province, Yan&#x2019;an University</institution>, <addr-line>Yan&#x2019;an</addr-line>,&#xa0;<country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/388702/overview">Rajesh Kumar Pathak</ext-link>, Chung-Ang University, Republic of Korea</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1335549/overview">Sujata Thakur</ext-link>, National Bureau of Plant Genetic Resources (ICAR), India</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1386382/overview">Parva Sharma</ext-link>, University of Maryland, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3116604/overview">Panrong Ren</ext-link>, Longdong University, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Shipeng Liu, <email xlink:href="mailto:Liushipeng2003@126.com">Liushipeng2003@126.com</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>09</day>
<month>09</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>16</volume>
<elocation-id>1652613</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>06</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>19</day>
<month>08</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Liu, Liang, Jing, Guo, Wang, Ma, Wang, Wang, He, Zhao, Jiang and Zhang.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Liu, Liang, Jing, Guo, Wang, Ma, Wang, Wang, He, Zhao, Jiang and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>
<italic>COBRA-Like</italic> (<italic>CBL</italic>) genes encode glycosylphosphatidylinositol (GPI) -anchored proteins specific to plants that play important roles in cellulose biosynthesis in primary and secondary cell walls.</p>
</sec>
<sec>
<title>Methods</title>
<p>This study used a bioinformatics approach to characterize the <italic>CBL</italic> family genes in <italic>Sorghum bicolor</italic> (<italic>S. bicolor</italic>) at the genome-wide level to investigate their potential functions in <italic>S. bicolor</italic> development.</p>
</sec>
<sec>
<title>Results</title>
<p>The results revealed the identification of 10 <italic>CBL</italic> genes in the BTx623 and E048 <italic>S. bicolor</italic> genomes, respectively. A comparative analysis of conserved Motifs revealed that all <italic>CBL</italic> family genes in <italic>S. bicolor</italic> possess CCVS conserved structural domains. Phylogenetic analysis revealed that the family can be divided into two subfamilies, with genes within each subfamily exhibiting similar gene structures and physicochemical properties. Whole Genome Duplication (WGD) played an important role in the expansion of <italic>SbCBL</italic> gene family. The tissue-specific expression patterns of <italic>SbCBL</italic> genes suggest varying expression levels across different organs and tissues in <italic>S. bicolor</italic>, with <italic>SbCBL1</italic>, <italic>SbCBL5</italic>, and <italic>SbCBL9</italic> showing significantly higher expression levels in roots. PEG and NaCl treatments significantly affected <italic>SbCBL</italic> expression levels. <italic>SbCBL4</italic> expression increased after PEG treatment, while <italic>SbCBL9</italic> expression decreased after NaCl treatment.</p>
</sec>
<sec>
<title>Conclusions</title>
<p>Overall, this study provides new insights into the role of the <italic>CBL</italic> gene family in <italic>S. bicolor</italic>.</p>
</sec>
</abstract>
<kwd-group>
<kwd>
<italic>SbCBL</italic> gene family</kwd>
<kwd>systematic evolution</kwd>
<kwd>whole genome duplication</kwd>
<kwd>abiotic stresses</kwd>
<kwd>
<italic>SbCBL4/9</italic>
</kwd>
</kwd-group>
<counts>
<fig-count count="9"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="37"/>
<page-count count="14"/>
<word-count count="6564"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Bioinformatics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The <italic>COBRA-Like</italic> (<italic>CBL</italic>) gene family encodes a class of glycosylphosphatidylinositol (GPI)-anchored proteins that act as cell surface receptors localized directly to the outer surface of plant plasma membranes and are involved in the perception and transduction of cell wall remodeling signals (<xref ref-type="bibr" rid="B29">Schindelman et&#xa0;al., 2001</xref>; <xref ref-type="bibr" rid="B27">Roudier et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B25">Ringli, 2010</xref>). <italic>CBL</italic> gene family was originally identified in mutants of <italic>Arabidopsis thaliana</italic> (<italic>A. thaliana</italic>) root cells that are abnormally expanded (<xref ref-type="bibr" rid="B4">Benfey et&#xa0;al., 1993</xref>). The number of members of the <italic>CBL</italic> gene family varies considerably among species: <italic>A. thaliana</italic>, <italic>Oryza sativa</italic> (<italic>O. sativa</italic>), <italic>Zea mays</italic> (<italic>Z. mays</italic>), and <italic>Gossypium hirsutum</italic> (<italic>G. hirsutum</italic>) contain 12, 11, 9, and 39 members, respectively (<xref ref-type="bibr" rid="B27">Roudier et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B21">Li et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B5">Brady et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B12">Fu et&#xa0;al., 2024</xref>). Typical structural features include an N-terminal signal peptide that mediates endoplasmic reticulum localization, an aromatic amino acid-rich region that constitutes the cellulose binding site, a central CCVS structural domain (containing cysteine clusters) that maintains protein conformation, and a C-terminal GPI-anchored signal sequence (containing a &#x3c9;-site) that mediates membrane localization (<xref ref-type="bibr" rid="B27">Roudier et&#xa0;al., 2002</xref>).</p>
<p>The <italic>cob</italic> mutant in <italic>A. thaliana</italic> exhibits abnormal root cell&#xa0;expansion and a dwarfing phenotype (<xref ref-type="bibr" rid="B26">Roudier et&#xa0;al., 2005</xref>). <italic>In situ</italic> hybridization results showed that the <italic>COB</italic> gene is highly expressed in the root elongation region, suggesting its involvement in regulating cellulose deposition during cell elongation (<xref ref-type="bibr" rid="B26">Roudier et&#xa0;al., 2005</xref>). Further studies revealed that <italic>COB</italic> gene mutations disrupt the orientation of cellulose microfilament arrangement, triggering a decrease in cellulose content (<xref ref-type="bibr" rid="B26">Roudier et&#xa0;al., 2005</xref>). The <italic>COBL4</italic> mutant, homologous to <italic>COB</italic>, exhibited a significant decrease in secondary wall cellulose content (<xref ref-type="bibr" rid="B6">Brown et&#xa0;al., 2005</xref>). In addition to CBL, five other cloned <italic>A. thaliana</italic> family members (e.g., <italic>AtCOBL4</italic>, <italic>AtCOBL9</italic> and <italic>AtCOBL2</italic>) are involved in cell wall synthesis. The <italic>AtCOBL4</italic> affects secondary wall cellulose synthesis (<xref ref-type="bibr" rid="B6">Brown et&#xa0;al., 2005</xref>); the <italic>AtCOBL9</italic> mutation results in defective polar root hair growth (<xref ref-type="bibr" rid="B16">Jones et&#xa0;al., 2006</xref>); and the <italic>AtCOBL2</italic> is involved in seed coat cellulose deposition (<xref ref-type="bibr" rid="B3">Ben et&#xa0;al., 2015</xref>). Notably, complete deletion of the <italic>COB</italic> gene triggers plant growth arrest and aberrant expression of defense-related genes (<xref ref-type="bibr" rid="B18">Ko et&#xa0;al., 2006</xref>). In monocotyledonous plants, a mutant of the <italic>O. sativa BC1</italic> gene (which encodes a CBL protein) exhibits a &#x201c;brittle rod&#x201d; phenotype and has 60.7% amino acid sequence homology with the <italic>A. thaliana</italic> COB protein (<xref ref-type="bibr" rid="B21">Li et&#xa0;al., 2003</xref>). This mutation results in reduced cell wall thickness, decreased cellulose content, abnormal lignin deposition, and decreased mechanical strength (<xref ref-type="bibr" rid="B21">Li et&#xa0;al., 2003</xref>). Wu Lab identified the <italic>OsBCL4</italic> gene, which encodes a CBL protein containing a typical GPI structural domain. Its T-DNA insertion mutant exhibits cell expansion, cellulose reduction, and pectin accumulation (<xref ref-type="bibr" rid="B10">Dai et&#xa0;al., 2011</xref>). In addition, the <italic>OsBCL5</italic> mutation affects male gametophyte transport (<xref ref-type="bibr" rid="B11">Dai et&#xa0;al., 2009</xref>), OsBC5 is involved in secondary wall formation in stem nodes (<xref ref-type="bibr" rid="B1">Aohara et&#xa0;al., 2009</xref>), and OsBC6 encodes a CESA-related protein that affects secondary wall synthesis (<xref ref-type="bibr" rid="B19">Kotake et&#xa0;al., 2011</xref>). In&#xa0;maize, ZmBK2L3, a member of the ZmBK2L family, is closely related to AtCOB, and its encoded protein retains conserved structural domains and is widely expressed at different developmental stages (<xref ref-type="bibr" rid="B5">Brady et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B7">Cao et&#xa0;al., 2012</xref>). Recent studies have confirmed that ZmBK2L1 is involved in the regulation of root hair development (<xref ref-type="bibr" rid="B15">Hochholdinger et&#xa0;al., 2008</xref>). Studies have shown that COBL genes play a role in how plants respond to abiotic stresses. For instance, DROT1 encodes a COBL protein. It was discovered that DROT1 increases cellulose content and maintains cellulose crystallinity in rice. This modulates the cell wall structure and enhances the plant&#x2019;s resistance to drought (<xref ref-type="bibr" rid="B33">Sun et&#xa0;al., 2022</xref>). In <italic>A. thaliana</italic>, overexpressing <italic>PtCOBL12</italic> promotes plant growth and increases cellulose content and relative crystallinity (<xref ref-type="bibr" rid="B13">Geng et&#xa0;al., 2023</xref>). It also improves growth under drought stress conditions. Additionally, the <italic>GhCOBL22</italic> gene plays a pivotal role in cotton&#x2019;s response to drought stress (<xref ref-type="bibr" rid="B12">Fu et&#xa0;al., 2024</xref>).</p>
<p>
<italic>Sorghum bicolor</italic> (<italic>S. bicolor</italic>), as the fifth largest cereal crop in the world, is an important energy and forage crop with excellent agronomic traits such as high photosynthetic efficiency, high nutritional value, high adaptability, and resistance to drought and salinity, and an important model crop for the study of other energy crops (<xref ref-type="bibr" rid="B31">Silva et&#xa0;al., 2022</xref>). BTx623 and E048 are two distinct sorghum varieties. BTx623 was the first variety for which a high-quality whole genome was sequenced and assembled (<xref ref-type="bibr" rid="B24">Paterson et&#xa0;al., 2009</xref>). Its genome has become a universal &#x201c;reference template&#x201d; for sorghum research, providing an important foundation for gene targeting and editing (e.g., CRISPR) as well as functional validation. E048, on the other hand, is derived from Sudanese tropical germplasm (Early Hegari) and differs significantly from BTx623 in terms of disease resistance gene clusters and metabolic pathways. These differences make E048 ideal for comparative genomics studies. For these reasons, BTx623 and E048 were selected for this study. In this study, we identified members of the <italic>S. bicolor CBL</italic> gene family at the genome-wide level using bioinformatics methods, analyzed their gene structures, evolutionary relationships, selective pressures, and expression patterns, and laid the groundwork for elucidating the functions of this gene family in the stress response of <italic>S. bicolor</italic>.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Identification of gene family members</title>
<p>In this study, HMMER 3.3.2 software was used to perform a homology search (E-value threshold of 1e<sup>-5</sup>) based on the Hidden Markov Model (HMM) of the CBL domain (PF04833) in the Pfam database (<ext-link ext-link-type="uri" xlink:href="https://pfam.xfam.org/">https://pfam.xfam.org/</ext-link>) (<xref ref-type="bibr" rid="B23">Mistry et&#xa0;al., 2021</xref>) for <italic>S. bicolor</italic> BTx623 and E048 protein sequences downloaded from the Phytozome v13 (<ext-link ext-link-type="uri" xlink:href="https://phytozome-next.jgi.doe.gov/">https://phytozome-next.jgi.doe.gov/</ext-link>) (<xref ref-type="bibr" rid="B14">Goodstein et&#xa0;al., 2012</xref>) and SGMD databases (<ext-link ext-link-type="uri" xlink:href="https://S.bicolor.genetics.ac.cn/SGMD">https://S.bicolor.genetics.ac.cn/SGMD</ext-link>, accessed on 5 May 2025) (<xref ref-type="bibr" rid="B8">Chen et&#xa0;al., 2025</xref>), respectively. The gene sequences obtained from the initial screening were further validated for conserved structural domains using NCBI Conserved Domain Database (CDD) (<ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi">https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi</ext-link>) to ensure that the identified sequences contained complete and typical CBL domains, thus accurately identifying members of the <italic>S. bicolor CBL</italic> gene family.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Analysis of gene structure and conserved motifs</title>
<p>The coding sequences (CDS) of the genes were aligned with the corresponding genomic sequences and visualized using TBtools software (<xref ref-type="bibr" rid="B9">Chen et&#xa0;al., 2023</xref>) to show the exon-intron structure of the <italic>SbCBL</italic> genes. Then, the conserved motifs of <italic>S. bicolor</italic> CBL proteins were predicted using MEME Suite 5.5.3 (<ext-link ext-link-type="uri" xlink:href="https://meme-suite.org/meme/tools/meme">https://meme-suite.org/meme/tools/meme</ext-link>) (<xref ref-type="bibr" rid="B2">Bailey et&#xa0;al., 2009</xref>) software, setting the maximum number of motifs to 10, the motif length range from 6-50 amino acids, and other parameters as default. In addition, the conserved structural domains of the <italic>S. bicolor</italic> CBL gene family were predicted by NCBI-CDD. Finally, the gene structures, conserved motifs, and conserved structural domains of <italic>S. bicolor CBL</italic> gene family members were comprehensively analyzed using TBtools software (<xref ref-type="bibr" rid="B9">Chen et&#xa0;al., 2023</xref>).</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Phylogenetic tree construction</title>
<p>The amino acid sequences of <italic>S. bicolor CBL</italic> gene family members and CBL protein sequences of <italic>A. thaliana</italic>, <italic>O. sativa</italic>, <italic>Z. mays</italic>, <italic>Solanum lycopersicum</italic> (<italic>S. lycopersicum</italic>) and <italic>Setaria italica</italic> (<italic>S. italica</italic>) were subjected to multiple sequence comparison using MEGA 7.0 software (<xref ref-type="bibr" rid="B20">Kumar et al., 2016</xref>). Neighbor-joining (NJ) was then used to construct the phylogenetic tree. In the parameter settings, the number of bootstrap tests was 1000, and the Poisson correction model was selected to calculate the genetic distance. Finally, the evolutionary tree was embellished by iTOL (<ext-link ext-link-type="uri" xlink:href="https://itol.embl.de/">https://itol.embl.de/</ext-link>) online website to show the evolutionary relationship between <italic>S. bicolor CBL</italic> gene family and members of this gene family in other plants, and to analyze the evolutionary pattern and classification of the gene family.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Analysis of replication events and selection pressure</title>
<p>Tandem and genome-wide replication events of the <italic>S. bicolor</italic> CBL gene family were analyzed using MCScan X (<xref ref-type="bibr" rid="B34">Wang et&#xa0;al., 2012</xref>). The Simple Ka/Ks Calculator module in TBtools software was used to input the coding sequences (CDS), protein sequences, and immediately homologous gene pairs of the genes, respectively, and the Ka (non-synonymous substitution rate)/Ks (synonymous substitution rate) values between homologous genes were calculated to estimate the selection pressure. In addition, MCScan X (<xref ref-type="bibr" rid="B34">Wang et&#xa0;al., 2012</xref>) was applied to analyze the covariation events between <italic>S. bicolor</italic> and <italic>A. thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, and <italic>S. italica CBL</italic> gene families.</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Gene expression analysis</title>
<p>Transcriptome data of different tissues of <italic>S. bicolor</italic> (including seedlings, leaves, roots, stems, inflorescences, and seeds, etc.) at different developmental stages were obtained from the <italic>S. bicolor</italic> Genome and Mutant Bank SGMD database (<xref ref-type="bibr" rid="B8">Chen et&#xa0;al., 2025</xref>). Then, we utilized Tbtools software to analyze gene expression data from various tissues, screening for CBL genes that exhibited high expression levels across different tissues and developmental stages. The data are shown in a heatmap with gene expression in different tissues with row-scaled transcriptome atlas (TPM values). Red and blue boxes indicate high and low expression levels of <italic>SbCBL</italic> genes.</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>RT-qPCR</title>
<p>
<italic>S. bicolor</italic> was grown in a growth chamber at Yan&#x2019;an University with 16.0 hours of light, temperature maintained at 25&#xb0;C and 70% humidity. To determine the expression level of <italic>CBL</italic> gene after NaCl and PEG treatments, <italic>S. bicolor</italic> seedlings at the three-leaf-one-heart stage were selected and treated with 150 mM NaCl and 15% PEG, respectively, followed by collection of <italic>S. bicolor</italic> root samples 7 days. All experiments were performed in three biological replicates (Three biological replicates and three technical replicates per biological sample were performed.). Total RNA was isolated using the Plant Total RNA Kit from Beijing Zhuangmeng International BioGenetics Co. Ltd. and reverse transcription was performed using HiScript IV All-in-One Ultra RT SuperMix for qPCR from Novozymes. The RT-qPCR amplification reaction system consisted of 5 &#x3bc;L 2&#xd7;SYBR, 3 &#x3bc;L ddH<sub>2</sub>O, 1 &#x3bc;L cDNA template, and 0.5 &#x3bc;M forward and reverse primers in a total volume of 10 &#x3bc;L. The expression level of <italic>SbCBL</italic> gene was analyzed by the 2<sup>-&#x394;&#x394;CT</sup> method in response to different stress treatments (<xref ref-type="bibr" rid="B35">Xue et&#xa0;al., 2024</xref>), and <italic>SbACTIN</italic> was used as an internal reference gene to analyze the expression level of <italic>SbCBL</italic> gene under different stress treatments.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Result</title>
<sec id="s3_1">
<label>3.1</label>
<title>Identification of <italic>S. bicolor CBL</italic> gene family members</title>
<p>Using HMMER 3.3.2 software based on the Hidden Markov Model of the CBL conservative domain (PF04833) in the Pfam database, we searched and validated <italic>S. bicolor</italic> protein sequences using NCBI-CDD. Finally, we identified ten <italic>CBL</italic> gene family members. These genes were named <italic>SbCBL1</italic>-<italic>SbCBL10</italic> and <italic>SbECBL1</italic>-<italic>SbECBL10</italic> based on their location on the chromosomes. The distribution of these 20 genes on the <italic>S. bicolor</italic> chromosomes was visualized using TBtools software, and the results are shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>. <italic>S. bicolor CBL</italic> gene family members were unevenly distributed across three chromosomes: five genes (<italic>SbCBL1</italic>-<italic>SbCBL5</italic>) were found on chromosome 1; four genes were found on chromosome 2; and one gene was found on chromosome 6. This uneven distribution pattern may be related to the evolution of gene families, chromosome structure, and gene function.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Chromosomal distribution of <italic>CBL</italic> genes in <italic>S. bicolor</italic>. <bold>(A)</bold> Chromosomal distribution of <italic>CBL</italic> genes in <italic>S. bicolor</italic> BTx623. <bold>(B)</bold> Chromosomal distribution of <italic>CBL</italic> genes in <italic>S. bicolor</italic> E048.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g001.tif">
<alt-text content-type="machine-generated">Diagram illustrating the genomic locations of CBL and ECBL genes on chromosomes. Panel A shows positions of SbCBL genes on chromosomes Chr01, Chr02, and Chr06. Panel B displays the SbECBL genes on the same chromosomes, with their respective locations labeled in red. Chromosome lengths are marked in megabases (Mb) along the left side.</alt-text>
</graphic>
</fig>
<p>This study systematically analyzed the amino acid length, molecular weight, isoelectric point, instability index, and hydrophilicity of these members. The results showed that the lengths of the amino acid sequences encoded by the <italic>S. bicolor</italic> CBL gene family members differed significantly. For instance, SbCBL1 encodes 187 amino acids, whereas SbCBL7 encodes a protein consisting of 673 amino acids (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2A</bold>
</xref>). Additionally, the amino acid lengths of CBL homologous genes differed between the two <italic>S. bicolor</italic> varieties. For instance, SbCBL9 corresponds to SbECBL9, which has 446 and 686 amino acids, respectively. This difference in length may reflect the functional complexity of different members. In terms of molecular weight, members of the <italic>S. bicolor</italic> CBL gene family exhibited similar diversity. The results showed that family proteins have a wide range of molecular weights. For example, SbCBL1 has a molecular weight of 46.832 kDa, while SbCBL7 has a molecular weight as high as 74.62 kDa. There are also members with a predicted molecular weight of about 36.72 kDa (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2B</bold>
</xref>). This difference in molecular weight may be closely related to the structural and functional diversity of the proteins. Notably, 70% of the <italic>SbCBL</italic> gene family members in Sorghum bicolor exhibit an isoelectric point (pI) exceeding 7. This prevalence of basic pI values underscores the potential functional adaptations of these proteins for operating in alkaline cellular environments or interacting with negatively charged macromolecules (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2C</bold>
</xref>). Additionally, 70% of SbCBL proteins were predicted to be unstable (PI &#x2265; 40), while 30% were stable (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2D</bold>
</xref>). This result establishes a foundation for the subsequent in-depth study of this family&#x2019;s protein functions. Through computational predictions of subcellular localization, all members of the SbCBL protein family exhibited exclusive targeting to the plasma membrane, indicating a highly conserved localization pattern. This striking uniformity in membrane association strongly suggests that SbCOBL proteins may play specialized roles in cell wall-plasma membrane interactions, potentially regulating cellulose deposition patterns or mediating mechanical stress responses at the cell surface.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Analysis of physicochemical properties of <italic>S. bicolor</italic> CBL protein. <bold>(A)</bold> Statistics of amino acid length of <italic>S. bicolor</italic> CBL protein. <bold>(B)</bold> Statistics of molecular weights of <italic>S. bicolor</italic> CBL protein. <bold>(C)</bold> Statistics of Isoelectric point of <italic>S. bicolor</italic> CBL protein. <bold>(D)</bold> Statistics of Instability Index of <italic>S. bicolor</italic> CBL protein.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g002.tif">
<alt-text content-type="machine-generated">Four pie charts, labeled A to D, show distributions of four characteristics. A: Protein lengths with 65% in 400-500 amino acids, 30% in &gt;600 amino acids, and 5% in 300-400 amino acids. B: Molecular weights with 45% in 40-50 kDa, 30% in &gt;60 kDa, and 25% in 50-60 kDa. C: Isoelectric points with 70% in 7-9.02 and 30% in 6-7. D: Size range with 70% in 29-40 II and 30% in 40-42 II.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Systematic evolutionary analysis</title>
<p>To resolve the evolutionary classification of the <italic>S. bicolor CBL</italic> gene family, this study used MEGA 7.0 software to conduct a phylogenetic analysis of the amino acid sequences of 10 SbCBL and 10 SbECBL members. The results showed that the CBL proteins of the two <italic>S. bicolor</italic> varieties can be categorized into three subfamilies: Group I, Group II, and Group III (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>). Group I contains 12 members: SbCBL1, SbCBL2, SbCBL3, SbCBL4, SbCBL7, SbCBL8, SbECBL1, SbECBL2, SbECBL3, SbECBL4, SbECBL7, and SbECBL8. Group II includes SbCBL10 and SbECBL10, and Group III includes SbCBL5, SbCBL6, SbCBL9, SbECBL5, SbECBL6, and SbECBL9.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Phylogenetic analysis of CBL proteins. <bold>(A)</bold> Phylogenetic analysis of CBL proteins in <italic>S. bicolor</italic>. The evolutionary tree was constructed using the neighbor-joining (NJ) method with 1,000 bootstrap replicates in MEGA 7. CBL proteins were classified into three subfamilies: Group I, Group II, and Group III, represented by purple, blue and pink branches, respectively. <bold>(B)</bold> Phylogenetic analysis of CBL proteins in <italic>S. bicolor</italic>, <italic>A thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>Z. mays</italic>. The evolutionary tree was constructed using the neighbor-joining (NJ) method with 1,000 bootstrap replicates in MEGA 7. CBL proteins were classified into three subfamilies: Group I, Group II, and Group III, represented by purple, blue and pink branches, respectively.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g003.tif">
<alt-text content-type="machine-generated">Phylogenetic analysis diagrams displaying two different representations. Diagram A shows a linear tree with three color-coded groups: Group I in purple, Group II in blue, and Group III in pink, each with labeled branches like SbCBL1 and SbCBL6. Diagram B is a circular tree with similar color coding and labels, organizing the data into three distinct radial groups.</alt-text>
</graphic>
</fig>
<p>To further elucidate the phylogenetic position of the <italic>S. bicolor CBL</italic> genes, this study integrated CBL protein sequences from <italic>S. bicolor</italic>, <italic>A. thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>Z. mays</italic> for a multiple comparison analysis (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3B</bold>
</xref>). Phylogenetic trees constructed using the neighbor-joining method showed that the <italic>CBL</italic> genes of all species could be divided into three significant branches (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4B</bold>
</xref>). Group I mainly included six SbCBLs, six SbECBLs, five AtCBLs, ten SlCBLs, seven OsCBLs, five ZmCBLs, and six SiCBLs. Group II mainly consisted of three SbCBLs, three SbECBLs, five AtCBLs, five SlCBLs, three OsCBLs, three ZmCBLs, and two SiCBLs. Group III mainly consisted of one SbCBL, one SbECBL, one AtCBL, two SlCBLs, one OsCBL, one ZmCBL, and one SiCBL. The results of the phylogenetic analysis showed that the <italic>S. bicolor CBL</italic> gene family has evolved to be related to <italic>CBL</italic> genes from other plants while maintaining its own specificity. ZmBk2L3, a COBRA family protein, functions in the regulation of cell wall dynamics and carbohydrate partitioning (<xref ref-type="bibr" rid="B17">Julius et&#xa0;al., 2021</xref>). Phylogenetic analysis revealed that its clustering with SbCBL1/2 and SbECBL1/2 suggests involvement in cell membrane-associated signaling or cell wall modification via similar mechanisms. ZmBk2 maintains the flexibility of plant organs by modulating the lignin-cellulose interaction pattern (<xref ref-type="bibr" rid="B32">Sindhu et&#xa0;al., 2007</xref>). Additionally, phylogenetic analyses revealed that ZmBk2 clusters with SbECBL4 and SbCBL4, indicating that the protein may perform similar functions via conserved molecular mechanisms.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Phylogenetic, conserved motif, domain, and gene structure analysis of <italic>S. bicolor</italic>, <italic>A thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>Z. mays</italic> CBL proteins. <bold>(A)</bold> Phylogenetic analysis of <italic>S. bicolor</italic>, <italic>A thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>Z. mays</italic> CBL proteins. The neighbor-joining (NJ) tree was constructed using MEGA 7 with 1,000 bootstrap replicates. <bold>(B)</bold> Intron-exon structure of <italic>S. bicolor</italic>, <italic>A. thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S.&#xa0;italica</italic> and <italic>Z. mays CBL</italic> genes. Visualization was performed using TBtools. <bold>(C)</bold> Conserved motif analysis of <italic>S. bicolor</italic>, <italic>A. thaliana</italic>, <italic>S. lycopersicum</italic>, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>Z. mays CBL</italic> proteins. Ten motifs were identified using the online tool MEME with default parameters.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g004.tif">
<alt-text content-type="machine-generated">Phylogenetic tree and gene structure visualization. Panel A shows a phylogenetic tree with genes grouped into three color-coded sections. Panel B displays gene structures with UTRs in green and CDS in yellow. Panel C shows motifs with various colors, each labeled from Motif 1 to Motif 10.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Gene structure and conserved motif characterization</title>
<p>A gene structure and conserved Motif analysis of <italic>S. bicolor CBL</italic> genes revealed significant structural and functional differences. Gene&#xa0;structure revealed notable variations in exon-intron structure among different subgroups. Specifically, members of subgroup I have three to seven exons and exhibit structural diversity. Members of subgroup II, such as SbCBL10 and SbECBL10, have six and four exons, respectively. Notably, subgroup III members have a more concise gene structure with only one or two exons. This diversity in gene structure may be closely related to functional differentiation and the evolutionary history of genes (<xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4A, B</bold>
</xref>). Different exon-intron structures may lead to variations in gene transcription and translation processes, which may affect gene function and expression regulation.</p>
<p>Conserved motif analysis performed by MEME Suite 5.5.3 identified a total of 10 characteristic motifs (Motif1-Motif10). Systematic analysis revealed that all motifs were intact in subgroup I, suggesting that these core elements may collectively maintain the basal biological functions of CBL proteins. Subgroup II exhibited a distinct motif deletion pattern: SbCOBL10 lacked Motif9, and SbECOBL10 lacked Motifs 1 and 9 (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>). Members of the third subgroup contained only Motifs 1, 2, 3, and 4, which were distributed in a manner that may be related to the functional differentiation of gene family members. The distribution of motifs 1, 2, 3, and 4 is specific and may be related to the functional differentiation of gene family members (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>). Analysis of gene structure and conserved motifs revealed that <italic>S. bicolor CBL</italic> gene family members are structurally conserved yet diverse. These conserved structural features may ensure the gene family&#x2019;s basic function, while diverse structures provide the basis for the genes&#x2019; functional differentiation and evolution, enabling different gene members to play unique roles in <italic>S. bicolor</italic> growth, development, and environmental adaptation.</p>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Analysis of duplication events</title>
<p>To reveal the expansion mechanism of the <italic>S. bicolor CBL</italic> gene family and its evolutionary constraints, this study systematically analyzed the types of replication events and selection pressures that characterize its members. As shown in the <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>, tandem duplication (TD) events drove the clustered distribution of CBL3 and CBL4. Meanwhile, Whole Genome Duplication (WGD) events contributed to the generation of two paralogous gene pairs: CBL3/CBL7 and CBL5/CBL9.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Analysis of gene duplication events of <italic>CBL</italic> genes in <italic>S. bicolor</italic>. <bold>(A)</bold> Gene duplication events of <italic>S. bicolor</italic> BTx623 <italic>CBL</italic> genes. Whole-genome duplication (WGD) events are indicated by purple lines, and tandem duplicated genes are labeled with purple gene IDs. <bold>(B)</bold> Gene duplication events of <italic>S. bicolor</italic> E048 <italic>CBL</italic> genes. Whole-genome duplication (WGD) events are indicated by purple lines, and tandem duplicated genes are labeled with purple gene IDs.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g005.tif">
<alt-text content-type="machine-generated">Circular diagrams labeled A and B display chromosomes numbered 1 to 10 in orange segments around the circle. Black lines indicate connections between genes labeled as SbCBL1 to SbCBL10 and SbECBL1 to SbECBL10. Purple curves highlight specific gene interactions.</alt-text>
</graphic>
</fig>
<p>Further analysis of the ratio of non-synonymous to synonymous substitution rates (Ka/Ks) revealed that the Ka/Ks values of these replication events were significantly less than one (CBL3/CBL7: Ka/Ks = 0.1; CBL5/CBL9: Ka/Ks = 0.27), suggesting that these genes underwent strong purifying selection (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;2</bold>
</xref>). This indicates that these genes experienced strong purifying selection during evolution.</p>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Synteny analysis of <italic>S. bicolor</italic> and other species</title>
<p>To elucidate the evolutionary trajectory of the <italic>CBL</italic> gene family in <italic>S. bicolor</italic>, this study performed whole-genome synteny analysis using MCScan X between S. bicolor and representative species, including the dicot model plant <italic>A. thaliana</italic>, monocot crops <italic>O. sativa</italic> and <italic>S. italica</italic>, as well as the solanaceous crop <italic>S. lycopersicum</italic>. A total of 3 (<italic>S. bicolor</italic>&#x2013;<italic>A. thaliana</italic>), 10 (<italic>S. bicolor</italic>&#x2013;<italic>S. lycopersicum</italic>), 10 (<italic>S. bicolor</italic>&#x2013;<italic>O. sativa</italic>), and 10 (<italic>S. bicolor</italic>&#x2013;<italic>S. italica</italic>) <italic>CBL</italic> homologous gene pairs were identified (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). Notably, <italic>O. sativa</italic>, <italic>S. italica</italic> and <italic>S. bicolor</italic>, as closely related species within the Gramineae family, exhibited significantly more syntenic gene pairs than <italic>A. thaliana</italic>, indicating that the <italic>CBL</italic> gene family retained higher genomic structural conservation after monocot&#x2013;dicot divergence. Further analysis revealed &#x201c;one-to-many&#x201d; homologous relationships between certain <italic>S. bicolor CBL</italic> genes and multiple species. For instance, <italic>SbCBL3</italic> showed synteny with both <italic>SlCBL1</italic> and <italic>SlCBL1</italic> in <italic>S. lycopersicum</italic>. Similarly, <italic>SbCBL3</italic> corresponded to three homologs in <italic>S. italica</italic> (<italic>SiCBL2</italic> and <italic>SiCBL7</italic>) and two in <italic>O. sativa</italic> (<italic>OsCBL1</italic> and <italic>OsCBL6</italic>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Synteny analysis of <italic>CBL</italic> genes with other species. Synteny between <italic>S. bicolor</italic> BTx623 and <italic>A. thaliana CBL</italic> genes. Synteny between <italic>S. bicolor</italic> and <italic>S. lycopersicum CBL</italic> genes. Synteny between <italic>S. bicolor</italic> and <italic>O. sativa CBL</italic> genes. Synteny between <italic>S. bicolor</italic> and <italic>S. italica CBL</italic> genes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g006.tif">
<alt-text content-type="machine-generated">Diagram illustrating synteny between the chromosomes of five plant species: Sorghum bicolor, Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa, and Setaria italica. Colored bars represent chromosomes, with connecting red lines showing syntenic relationships between the species' genomes.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3_6">
<label>3.6</label>
<title>Analysis of promoter cis-acting elements of the <italic>S. bicolor CBL</italic> gene family</title>
<p>To analyze the transcriptional and regulatory features of the <italic>S.&#xa0;bicolor CBL</italic> gene family, this study systematically analyzed the cis-acting elements in the upstream promoter region (2,000 bp prior to the transcriptional start site) of the gene using the PlantCARE database. The results showed that the <italic>CBL</italic> gene promoter region contained abundant regulatory elements, mainly categorized into four groups: hormone response, abiotic stress response, light signaling regulation, and growth and development (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;3</bold>
</xref>). All members carried abscisic acid (ABA) response elements (ABRE, ACGTG) and methyl jasmonate (MeJA) response elements (TGACG motif, TGACG), indicating that the <italic>S. bicolor CBL</italic> gene family may be involved in regulating adversity acclimatization through ABA and MeJA signaling pathways. Additionally, some members contained cis-regulatory elements related to growth hormones, salicylic acid, and gibberellin, implying that these members may be involved in multiple stress responses through hormone crosstalk. Abiotic stress response elements included the drought response element MBS (CAACTG), the low temperature response element LTR (CCGAAA), the anaerobic-induced element ARE (AAACCA), and the mechanical damage response element. Members of subgroup III (e.g., <italic>SbCBL7/9</italic>) showed a notably high frequency of the low-temperature response element LTR (CCGAAA) in the promoter region. This density was significantly higher than that of the other subgroups. This suggests that subgroup III may enhance cold hardiness by activating low-temperature acclimation pathways. Meanwhile, subgroups II and III specifically carried the mechanical damage response element, the WUN motif (AAATTACCT), which may respond to physical stress by regulating the cell wall.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Analysis of cis-regulatory elements in the promoter regions of <italic>S. bicolor CBL</italic> gene.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g007.tif">
<alt-text content-type="machine-generated">Phylogenetic tree and motif composition of sorghum SbCBL and SbEBCBL gene families. Gene names are listed alongside lines representing base pairs from 5 prime to 3 prime. Colored boxes indicate different motifs, with a legend showing colors for 24 motifs, including ABRE, AuxRR-core, CAT-box, and others.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3_7">
<label>3.7</label>
<title>Analysis of gene expression patterns</title>
<p>We analyzed the expression patterns of 10 <italic>S. bicolor CBL</italic> gene family members based on transcriptome data of different <italic>S. bicolor</italic> tissues (e.g., roots, stems, leaves, flowers, and seeds) at different developmental stages in the SGMD database (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). The data are presented in a heatmap illustrating gene expression across various tissues, utilizing a row-scaled transcriptome atlas based on TPM (Transcripts Per Million) values. In this visualization, red boxes signify high expression levels of <italic>SbCBL</italic> genes, while blue boxes indicate low expression levels. As illustrated in <xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>, various <italic>CBL</italic> genes exhibited significant expression variations across different <italic>S. bicolor</italic> tissues and developmental stages. In root tissues, the expression levels of <italic>SbCBL1</italic>, <italic>SbCBL5</italic> and <italic>SbCBL9</italic> were relatively high, suggesting that these genes play important roles in root growth and development. The normal development of roots is crucial for plant growth and survival because they are an important organ for water and nutrient uptake in plants. These genes may promote root growth and development by regulating cell wall synthesis and modification in root cells, as well as affecting cell elongation and differentiation. In stem tissues, <italic>SbCBL2</italic>, <italic>SbCBL4</italic> and <italic>SbCBL8</italic> exhibited high expression levels, suggesting their involvement in stem elongation and thickening processes. Stem growth and development play a key role in supporting the plant and transporting materials. These genes may regulate the arrangement and deposition of cellulose microfilaments in stems, enhancing their mechanical strength to support the plant during growth. <italic>SbCBL6</italic>, <italic>SbCBL7</italic> and <italic>SbCBL10</italic> are expressed at significantly higher levels in inflorescence tissues than in other tissues, suggesting that they play important roles in flower development and reproduction.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>The heat map shows the expression level of the <italic>S. bicolor CBL</italic> gene in different tissues. Red and blue boxes indicate high and low expression levels of <italic>SbCBL</italic> genes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g008.tif">
<alt-text content-type="machine-generated">Heatmap showing expression levels of SbCBL genes across different plant tissues, including root, seedling, leaves, stems, inflorescence, and seeds. Colors range from red (high expression) to blue (low expression). The legend on the right indicates expression range from -1 to 3.5.</alt-text>
</graphic>
</fig>
<p>The expression level of <italic>SbCBL3</italic> was significantly higher in the seed grain than in other tissues, suggesting that this gene plays an important role in seed grain development. Because material accumulation and maturation directly affect <italic>S. bicolor</italic> yield and quality, <italic>SbCBL3</italic> may influence seed grain size, shape, and starch accumulation by regulating cell wall synthesis and modification, ultimately affecting <italic>S. bicolor</italic> yield and quality. Gene expression pattern analysis revealed that members of the <italic>S. bicolor</italic> CBL gene family exhibit distinct expression patterns in various tissues and developmental stages. This differential expression may be closely related to the functional differentiation of these genes. These gene family members play unique roles in different stages of <italic>S. bicolor</italic> growth and development, as well as in different tissues. They are involved in plant growth, development, and reproduction by regulating cell wall synthesis and modification.</p>
</sec>
<sec id="s3_8">
<label>3.8</label>
<title>Expression profile of <italic>SbCBL</italic> in <italic>S. bicolor</italic> in response to NaCl and PEG</title>
<p>To explore the expression pattern of the <italic>CBL</italic> gene in <italic>S. bicolor</italic> under drought and salt stress, we selected a <italic>CBL</italic> gene with high expression levels in the roots and performed RT-qPCR analysis. As shown in <xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref>, the expression levels of several <italic>CBL</italic> genes changed significantly under NaCl and PEG treatment. Under NaCl treatment, the expression level of <italic>SbCBL4</italic> increased, while the expression levels of <italic>SbCBL3</italic>, <italic>SbCBL5</italic>, <italic>SbCBL7</italic>, <italic>SbCBL8</italic> and <italic>SbCBL9</italic> decreased. Conversely, PEG treatment up-regulated the expression of genes such as <italic>SbCBL4</italic> and <italic>SbCBL8</italic>, while down-regulating the expression of <italic>SbCBL1</italic>, <italic>SbCBL2</italic>, <italic>SbCBL3</italic>, <italic>SbCBL5</italic> and <italic>SbCBL9</italic>. Notably, the <italic>SbCBL3</italic>, <italic>SbCBL5</italic> and <italic>SbCBL9</italic> genes were down-regulated following both drought and salt stress treatments. These results suggest that the <italic>CBL</italic> gene in <italic>S. bicolor</italic> plays a role in the <italic>S. bicolor</italic> response to salt and drought stress.</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Analysis of the expression levels of <italic>CBL1</italic>, <italic>CBL2</italic>, <italic>CBL3</italic>, <italic>CBL4</italic>, <italic>CBL5</italic>, <italic>CBL6</italic>, <italic>CBL7</italic>, <italic>CBL8, CBL9</italic> and <italic>CBL10</italic> after treatment with NaCl and PEG. All data are means &#xb1; sd ( n &#x2265; 3). Letters a, b and c represent statistical significance, P &lt; 0.05.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-16-1652613-g009.tif">
<alt-text content-type="machine-generated">Bar graphs showing relative expression levels of CBL1 to CBL10 genes under Control, NaCl, and PEG conditions. Expression levels vary across treatments, indicated by letters a, b, and c, signifying statistical differences.</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<sec id="s4_1">
<label>4.1</label>
<title>Structural and functional speculation of the CBL gene family in <italic>S. bicolor</italic>
</title>
<p>In this study, 10 members of the <italic>CBL</italic> gene family were identified from the <italic>S. bicolor</italic> BTx623 and E048 genomes, respectively (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Gene structure analysis showed that the number of exons of these members ranged from 1-7, and there were differences in exon lengths and intron phases (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). Through conserved motif analysis, we identified 10 conserved motifs, of which Motif1, Motif2, Motif3 and Motif4 were distributed among all members, indicating that these motifs are important for maintaining the basic structure and function of CBL proteins (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). The structural features of genes are usually closely related to their functions (<xref ref-type="bibr" rid="B28">Sajjad et&#xa0;al., 2023</xref>). The diversity in exon-intron structure of <italic>S. bicolor CBL</italic> gene family members implies that they may be functionally differentiated. For example, <italic>SbCBL2</italic>, <italic>SbCBL3</italic>, <italic>SbCBL4</italic>, <italic>SbCBL7</italic> and <italic>SbCBL8</italic>, which have a higher number of exons, may encode proteins with more complex structures and functions, which are involved in the regulation of multiple processes in <italic>S. bicolor</italic> growth and development, whereas <italic>SbCBL5</italic>, <italic>SbCBL6</italic> and <italic>SbCBL9</italic>, which have a relatively low number of exons, may encode proteins that are simpler. This phenomenon has been observed in other species (<xref ref-type="bibr" rid="B28">Sajjad et&#xa0;al., 2023</xref>). Based on gene structure and conserved motif analyses combined with existing research reports, it is hypothesized that members of the <italic>S. bicolor CBL</italic> gene family are involved in several aspects of <italic>S. bicolor</italic> growth and development. However, these functional speculations require further experimental validation. Future studies could knock out or overexpress members of the <italic>S. bicolor CBL</italic> gene family using gene editing techniques, such as CRISPR-Cas9, to observe their effects on growth and development, cell wall structure, and related physiological processes. Simultaneously, proteomics and biochemical methods will be used to thoroughly study the interaction mechanism between CBL proteins and other molecules and clarify their specific functions and pathways of action within the cell.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Comparison with <italic>CBL</italic> gene families of other species</title>
<p>In this study, we identified a total of 10 <italic>CBL</italic> genes in <italic>S. bicolor</italic>, which aligns closely with the number of genes in this family reported in other monocotyledonous plants. For instance, <italic>O. sativa</italic> has 11 <italic>CBL</italic> genes (<xref ref-type="bibr" rid="B21">Li et&#xa0;al., 2003</xref>), <italic>Z. mays</italic> has 9 (<xref ref-type="bibr" rid="B5">Brady et&#xa0;al., 2007</xref>). In contrast, among dicotyledonous plants, <italic>A. thaliana</italic> possesses 12 <italic>CBL</italic> genes (<xref ref-type="bibr" rid="B27">Roudier et&#xa0;al., 2002</xref>), <italic>S. lycopersicum</italic> has 17 (<xref ref-type="bibr" rid="B7">Cao et&#xa0;al., 2012</xref>), <italic>P. trichocarpa</italic> 14 (<xref ref-type="bibr" rid="B28">Sajjad et&#xa0;al., 2023</xref>), and <italic>G. hirsutum</italic> has 39. This indicates a significant variation in the number of <italic>CBL</italic> family genes across different plant species.This suggests that the <italic>CBL</italic> gene family is highly conserved across monocotyledonous. Phylogenetic tree analysis revealed that the <italic>CBL</italic> genes of <italic>S. bicolor</italic> are interspersed with those of <italic>O. sativa</italic> and <italic>Z. mays</italic>, forming several small sub-branches (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). This indicates that <italic>S. bicolor</italic> is evolutionarily related to <italic>O. sativa</italic> and <italic>Z. mays</italic> in the <italic>CBL</italic> gene family and may share some functions with them. For instance, SbCOBL5, SbCOBL9, SbECOBL5, and SbECOBL9 clustered with OsBC1L1 and OsBC1L8, as well as ZmBK2L1 and ZmBK2L8, forming a sub-branch with some maize CBL genes (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<p>By comparing the conserved motifs of <italic>CBL</italic> genes in different species, it was found that several species have Motif1, Motif2, Motif3, and Motif4 (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). These motifs are found in CBL proteins from <italic>S. bicolor</italic>, <italic>A. thaliana</italic>, <italic>O. sativa</italic> and <italic>Z. mays</italic>. This suggests that they play important roles in the basic functions of the CBL gene family. They may be involved in the interaction of CBL proteins with other molecules or in maintaining the proteins&#x2019; structural stability. A comparative analysis of the CBL gene family with <italic>CBL</italic> gene families of other species reveals the evolutionary conservation and specificity of the <italic>S. bicolor CBL</italic> gene family. This conservation is reflected in similarities with other plants in terms of gene structure, conserved motifs, and evolutionary relationships. These similarities provide important clues for a deeper understanding of the CBL gene family&#x2019;s basic functions.</p>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Relationship between gene expression patterns and S. bicolor growth and development</title>
<p>Gene expression pattern analysis revealed that members of the <italic>S. bicolor CBL</italic> gene family exhibit distinct expression patterns in various tissues and developmental stages (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>). Root growth and development is a complex process involving cell division, elongation, and differentiation, and <italic>CBL</italic> genes may affect cell morphogenesis and physiological functions by regulating cell wall synthesis and modification in root cells (<xref ref-type="bibr" rid="B30">Shao et&#xa0;al., 2020</xref>). In <italic>A. thaliana</italic> and <italic>O. sativa</italic>, it has been shown that <italic>CBL</italic> genes regulate the initiation and growth of root hairs (<xref ref-type="bibr" rid="B26">Roudier et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B3">Ben et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B22">Li et&#xa0;al., 2022</xref>). In root tissues, the high expression of <italic>SbCBL1</italic>, <italic>SbCBL5</italic> and <italic>SbCBL9</italic> may regulate root growth and development (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>). For instance, these genes may regulate the orientation of cellulose microfilaments, affecting root cell elongation and differentiation. The highly expressed <italic>CBL</italic> genes in <italic>S. bicolor</italic> root tissues may have similar functions and play an important role in regulating root growth and development.</p>
<p>
<italic>SbCBL2</italic>, <italic>SbCBL4</italic>, and <italic>SbCBL8</italic> are highly expressed in stem tissues (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>). The mechanical strength of stems is crucial for upright growth and material transportation, as they are an important support structure in plants (<xref ref-type="bibr" rid="B37">Zhao et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B36">Zhang et&#xa0;al., 2022</xref>). <italic>SbCBL2</italic>, <italic>SbCBL4</italic>, and <italic>SbCBL8</italic> may enhance cell wall strength and toughness by regulating cellulose microfilament arrangement and deposition in the stem, thus promoting elongation and thickening. Phylogenetic analyses confirmed this hypothesis, showing that ZmBk2 maintains plant organ flexibility by regulating the interaction between lignin and cellulose (<xref ref-type="bibr" rid="B32">Sindhu et&#xa0;al., 2007</xref>). ZmBk2 clusters with SbECBL4 and SbCBL4, suggesting that it may perform similar functions through conserved molecular mechanisms. Different <italic>SbCBL</italic> genes showed differential expression patterns under drought and salt stress conditions (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref>). Specifically, NaCl treatment significantly increased the expression level of <italic>SbCBL4</italic> while decreasing the expression of <italic>SbCBL3</italic>, <italic>SbCBL5</italic>, <italic>SbCBL7</italic>, <italic>SbCBL8</italic> and <italic>SbCBL9</italic>. Conversely PEG treatment increased the expression of <italic>SbCBL4</italic> and <italic>SbCBL8</italic> while decreasing the expression of <italic>SbCBL1</italic>, <italic>SbCBL2</italic>, <italic>SbCBL3</italic>, <italic>SbCBL5</italic>, and <italic>SbCBL9</italic>. Analysis of the cis-acting elements of the <italic>SbCBL</italic> promoters revealed that <italic>SbCBL5</italic>, <italic>SbCBL7</italic>, <italic>SbCBL8</italic> and <italic>SbCBL9</italic> contain regulatory elements related to drought response. These results suggest that the <italic>SbCBL</italic> gene family plays a critical role in sorghum&#x2019;s response to drought and salt stress. This role has also been confirmed in other plants, such as rice and cotton (<xref ref-type="bibr" rid="B33">Sun et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B12">Fu et&#xa0;al., 2024</xref>). This finding is consistent with previous research indicating that <italic>CBL</italic> genes play important roles in drought and salt tolerance in other species. These findings provide important clues for understanding the molecular regulatory mechanisms of <italic>S. bicolor</italic> growth and development and potential gene targets for genetic improvement and molecular breeding of <italic>S. bicolor</italic>. Future studies can explore these genes&#x2019; specific mechanisms during <italic>S. bicolor</italic> growth and development through gene function validation experiments. This will provide stronger theoretical support for improving <italic>S. bicolor</italic> varieties and agricultural production.</p>
</sec>
<sec id="s4_4">
<label>4.4</label>
<title>Limitations and prospects of the study</title>
<p>In this study, a more comprehensive bioinformatics analysis&#xa0;of&#xa0;the <italic>CBL</italic> gene family of <italic>S. bicolor</italic> was conducted. However, certain&#xa0;limitations still exist. First, regarding gene function validation, this study only analyzed gene structure, evolution, and expression patterns using bioinformatics methods. Thus, experimental validation of gene functions has not yet been carried out. Although potential gene functions were hypothesized based on gene structure and expression patterns, these hypotheses must be verified using techniques such as gene editing and transgenesis. Second, due to limitations in experimental conditions and technical capabilities, this study only used transcriptomic data for gene expression analysis, lacking protein-level validation. It is important to note that gene expression ultimately reflects protein expression; however, transcriptome data only reflects changes in the transcription level of genes, which does not accurately reflect protein expression. Future research can be carried out in the following directions: 1. To verify gene function, use CRISPR-Cas9 and other&#xa0;gene-editing technologies to construct <italic>S. bicolor CBL</italic> gene family knockout mutants and overexpression plants. Through phenotyping and analyzing physiological and biochemical indexes, study the specific functions and mechanisms of the genes in <italic>S. bicolor</italic> growth, development, and response to adverse stress. 2. Combine proteome data with transcriptome data to study protein expression. Using proteomics technology, we will analyze protein expression in different <italic>S. bicolor</italic> tissues and developmental stages. This will allow us to verify the relationship between gene and protein expression and further clarify gene function.3. We will study the interactions between members of the <italic>CBL</italic> gene family and other genes. Then, we will construct a gene regulatory network to help us understand the molecular regulatory mechanisms of <italic>S. bicolor</italic> growth, development, and response to adversity.</p>
</sec>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions</title>
<p>A total of ten <italic>CBL</italic> genes were identified in the genomes of BTx623 and E048 S. bicolor. Phylogenetic analysis revealed that <italic>CBLs</italic> can be classified into three subfamilies: Group I, Group II, and Group III. Gene duplication events indicate that WGD was the primary driver of the expansion of the <italic>CBL</italic> gene family. The tissue-specific expression patterns of <italic>SbCBL</italic> genes suggest varying expression levels across different organs and tissues in <italic>S. bicolor</italic>, with <italic>SbCBL1</italic>, <italic>SbCBL5</italic>, and <italic>SbCBL9</italic> showing significantly higher expression levels in roots. Furthermore, treatments with PEG and NaCl markedly affected the expression levels of <italic>SbCBL</italic> genes; specifically, <italic>SbCBL4</italic> expression increased following PEG treatment, while <italic>SbCBL9</italic> expression decreased after NaCl treatment. Overall, this study provides valuable insights into the role of the <italic>CBL</italic> gene family in <italic>S. bicolor</italic>.</p>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material</bold>
</xref>.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>SPL: Funding acquisition, Supervision, Visualization, Writing &#x2013; original draft, Software. SL: Data curation, Formal Analysis, Validation, Writing &#x2013; review &amp; editing. TJ: Formal Analysis, Methodology, Software, Validation, Writing &#x2013; original draft. XG: Data curation, Methodology, Resources, Writing &#x2013; review &amp; editing. HW: Data curation, Resources, Software, Writing &#x2013; review &amp; editing.&#xa0;QM: Investigation, Resources, Validation, Writing &#x2013; review &amp; editing. JW: Data curation, Investigation, Validation, Writing &#x2013; review &amp; editing. KW: Conceptualization, Data curation, Investigation, Visualization, Writing &#x2013; review &amp; editing. XH: Conceptualization, Funding acquisition, Project administration, Supervision, Writing &#x2013; review &amp; editing. HZ: Funding acquisition, Supervision, Visualization, Writing &#x2013; review &amp; editing. WJ: Conceptualization, Funding acquisition, Supervision, Visualization, Writing &#x2013; review &amp; editing. XZ: Conceptualization, Funding acquisition, Project administration, Writing &#x2013; review &amp; editing.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research and/or publication of this article. Shaanxi Provincial Department of Science and Technology, General Project - Agricultural Fields 2023-YBNY-061; National Natural Science Foundation of China, Regional Fund, 32160761; Yan&#x2019;an University, Industry-University-Research Cooperation Cultivation Project, CXY202110; Shaanxi Provincial Department of Science and Technology, Key Industry Innovation Chain - Agriculture, 2023-ZDLNY-50.</p>
</sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="ai-statement">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2025.1652613/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2025.1652613/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table2.xls" id="SM2" mimetype="application/vnd.ms-excel"/>
<supplementary-material xlink:href="Table3.xls" id="SM3" mimetype="application/vnd.ms-excel"/>
<supplementary-material xlink:href="Table4.xls" id="SM4" mimetype="application/vnd.ms-excel"/>
<supplementary-material xlink:href="Table5.xls" id="SM5" mimetype="application/vnd.ms-excel"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aohara</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kotake</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Kaneko</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Takatsuji</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Tsumuraya</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Kawasaki</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Rice BRITTLE CULM 5 (BRITTLE NODE) is involved in secondary cell wall formation in the sclerenchyma tissue of nodes</article-title>. <source>Plant Cell Physiol.</source> <volume>50</volume>, <fpage>1886</fpage>&#x2013;<lpage>1897</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/pcp/pcp133</pub-id>, PMID: <pub-id pub-id-type="pmid">19812064</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bailey</surname> <given-names>T. L.</given-names>
</name>
<name>
<surname>Boden</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Buske</surname> <given-names>F. A.</given-names>
</name>
<name>
<surname>Frith</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Grant</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Clementi</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2009</year>). <article-title>MEME Suite: Tools for motif discovery and searching</article-title>. <source>Nucleic Acids Res.</source> <volume>37</volume>, <fpage>W202</fpage>&#x2013;<lpage>W208</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkp335</pub-id>, PMID: <pub-id pub-id-type="pmid">19458158</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ben</surname> <given-names>T. D.</given-names>
</name>
<name>
<surname>Abraham</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Stav</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Thompson</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Loraine</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Elbaum</surname> <given-names>R.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>COBRA-LIKE2, a member of the glycosylphosphatidylinositol-anchored COBRA-LIKE family, plays a role in cellulose deposition in arabidopsis seed coat mucilage secretory cells</article-title>. <source>Plant Physiol.</source> <volume>167</volume>, <fpage>711</fpage>&#x2013;<lpage>724</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.114.240671</pub-id>, PMID: <pub-id pub-id-type="pmid">25583925</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benfey</surname> <given-names>P. N.</given-names>
</name>
<name>
<surname>Linstead</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Roberts</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Schiefelbein</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Hauser</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Aeschbacher</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>Root development in Arabidopsis: four mutants with dramatically altered root morphogenesis</article-title>. <source>Development.</source> <volume>119</volume>, <fpage>57</fpage>&#x2013;<lpage>70</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1242/dev.119.1.57</pub-id>, PMID: <pub-id pub-id-type="pmid">8275864</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brady</surname> <given-names>S. M.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Dhugga</surname> <given-names>K. S.</given-names>
</name>
<name>
<surname>Rafalski</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Benfey</surname> <given-names>P. N.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Combining expression and comparative evolutionary analysis. The COBRA gene family</article-title>. <source>Plant Physiol.</source> <volume>143</volume>, <fpage>172</fpage>&#x2013;<lpage>187</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.106.087262</pub-id>, PMID: <pub-id pub-id-type="pmid">17098858</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Brown</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Zeef</surname> <given-names>L. A.</given-names>
</name>
<name>
<surname>Ellis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Goodacre</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Turner</surname> <given-names>S. R.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Identification of novel genes in Arabidopsis involved in secondary cell wall formation using expression profiling and reverse genetics</article-title>. <source>Plant Cell.</source> <volume>17</volume>, <fpage>2281</fpage>&#x2013;<lpage>2295</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1105/tpc.105.031542</pub-id>, PMID: <pub-id pub-id-type="pmid">15980264</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Giovannoni</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Functional characterization of a tomato COBRA-like gene functioning in fruit development and ripening</article-title>. <source>BMC Plant Biol.</source> <volume>12</volume>, <fpage>211</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/1471-2229-12-211</pub-id>, PMID: <pub-id pub-id-type="pmid">23140186</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Ge</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Du</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Sui</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2025</year>). <article-title>A comprehensive omics resource and genetic tools for functional genomics research and genetic improvement of sorghum</article-title>. <source>Mol. Plant</source> <volume>18</volume>, <fpage>703</fpage>&#x2013;<lpage>719</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2025.03.005</pub-id>, PMID: <pub-id pub-id-type="pmid">40055894</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Tbtools-II: A &#x201c;one for all, all for one&#x201d; bioinformatics platform for biological big-data mining</article-title>. <source>Mol. Plant</source> <volume>16</volume>, <fpage>1733</fpage>&#x2013;<lpage>1742</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2023.09.010</pub-id>, PMID: <pub-id pub-id-type="pmid">37740491</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname> <given-names>X.</given-names>
</name>
<name>
<surname>You</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>OsBC1L4 encodes a COBRA-like protein that affects cellulose synthesis in rice</article-title>. <source>Plant Mol. Biol.</source> <volume>75</volume>, <fpage>333</fpage>&#x2013;<lpage>345</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11103-011-9730-z</pub-id>, PMID: <pub-id pub-id-type="pmid">21264494</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname> <given-names>X.</given-names>
</name>
<name>
<surname>You</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Molecular characterization, expression pattern, and function analysis of the OsBC1L family in rice</article-title>. <source>Plant Mol. Biol.</source> <volume>71</volume>, <fpage>469</fpage>&#x2013;<lpage>481</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11103-009-9537-3</pub-id>, PMID: <pub-id pub-id-type="pmid">19688299</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname> <given-names>W. W.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Z. Y.</given-names>
</name>
<name>
<surname>Liusui</surname> <given-names>Y. H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>A. X.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>X. Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Genome-wide analysis of the cotton COBRA-like gene family and functional characterization of GhCOBL22 in relation to drought tolerance</article-title>. <source>BMC Plant Biol.</source> <volume>24</volume>, <fpage>1242</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-024-05965-x</pub-id>, PMID: <pub-id pub-id-type="pmid">39716062</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Geng</surname> <given-names>E. Y.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Genome-wide identification of the COBRA-Like gene family in Pinus tabuliformis and the role of PtCOBL12 in the regulation of cellulose biosynthesis</article-title>. <source>Ind. Crops Prod.</source> <volume>203</volume>, <fpage>117189</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.indcrop.2023.117189</pub-id>
</citation></ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goodstein</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Shu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Howson</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Neupane</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Hayes</surname> <given-names>R. D.</given-names>
</name>
<name>
<surname>Fazo</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Phytozome: A comparative platform for green plant genomics</article-title>. <source>Nucleic Acids Res.</source> <volume>40</volume>, <fpage>D1178</fpage>&#x2013;<lpage>D1186</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkr944</pub-id>, PMID: <pub-id pub-id-type="pmid">22110026</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hochholdinger</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Wen</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Zimmermann</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Chimot-Marolle</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Costa</surname> <given-names>S. O.</given-names>
</name>
<name>
<surname>Bruce</surname> <given-names>W.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>The maize (Zea mays L.) roothairless3 gene encodes a putative GPI-anchored, monocot-specific, COBRA-like protein that significantly affects grain yield</article-title>. <source>Plant J.</source> <volume>54</volume>, <fpage>888</fpage>&#x2013;<lpage>898</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-313X.2008.03459.x</pub-id>, PMID: <pub-id pub-id-type="pmid">18298667</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Raymond</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Smirnoff</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Analysis of the root-hair morphogenesis transcriptome reveals the molecular identity of six genes with roles in root-hair development in Arabidopsis</article-title>. <source>Plant J.</source> <volume>45</volume>, <fpage>83</fpage>&#x2013;<lpage>100</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-313X.2005.02609.x</pub-id>, PMID: <pub-id pub-id-type="pmid">16367956</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Julius</surname> <given-names>B. T.</given-names>
</name>
<name>
<surname>McCubbin</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Mertz</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Baert</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Knoblauch</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Grant</surname> <given-names>D. G.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Maize Brittle Stalk2-Like3, encoding a COBRA protein, functions in cell wall formation and carbohydrate partitioning</article-title>. <source>Plant Cell.</source> <volume>33</volume>, <fpage>3348</fpage>&#x2013;<lpage>3366</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/plcell/koab193</pub-id>, PMID: <pub-id pub-id-type="pmid">34323976</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ko</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Kim</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Jayanty</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Howe</surname> <given-names>G. A.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>K. H.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Loss of function of COBRA, a determinant of oriented cell expansion, invokes cellular defence responses in Arabidopsis thaliana</article-title>. <source>J. Exp. Bot.</source> <volume>57</volume>, <fpage>2923</fpage>&#x2013;<lpage>2936</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jxb/erl052</pub-id>, PMID: <pub-id pub-id-type="pmid">16873454</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotake</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Aohara</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Hirano</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Sato</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Kaneko</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tsumuraya</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Rice Brittle culm 6 encodes a dominant-negative form of CesA protein that perturbs cellulose synthesis in secondary cell walls</article-title>. <source>J. Exp. Bot.</source> <volume>62</volume>, <fpage>2053</fpage>&#x2013;<lpage>2062</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jxb/erq395</pub-id>, PMID: <pub-id pub-id-type="pmid">21209026</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumar</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Stecher</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Tamura</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets</article-title>. <source>Mol. Biol. Evol.</source> <volume>33</volume>, <fpage>1870</fpage>&#x2013;<lpage>1874</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/molbev/msw054</pub-id>, PMID: <pub-id pub-id-type="pmid">27004904</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Qian</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2003</year>). <article-title>BRITTLE CULM1, which encodes a COBRA-like protein, affects the mechanical properties of rice plants</article-title>. <source>Plant Cell.</source> <volume>15</volume>, <fpage>2020</fpage>&#x2013;<lpage>2031</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1105/tpc.011775</pub-id>, PMID: <pub-id pub-id-type="pmid">12953108</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Gong</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>P.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>COBL9 and COBL7 synergistically regulate root hair tip growth via controlling apical cellulose deposition</article-title>. <source>Biochem. Biophys. Res. Commun.</source> <volume>596</volume>, <fpage>6</fpage>&#x2013;<lpage>13</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.bbrc.2022.01.096</pub-id>, PMID: <pub-id pub-id-type="pmid">35104663</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mistry</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chuguransky</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Williams</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Qureshi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Salazar</surname> <given-names>G. A.</given-names>
</name>
<name>
<surname>Sonnhammer</surname> <given-names>E. L. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Pfam: The protein families database in 2021</article-title>. <source>Nucleic Acids Res.</source> <volume>49</volume>, <fpage>D412</fpage>&#x2013;<lpage>D419</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkaa913</pub-id>, PMID: <pub-id pub-id-type="pmid">33125078</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Paterson</surname> <given-names>A. H.</given-names>
</name>
<name>
<surname>Bowers</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Bruggmann</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Dubchak</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Grimwood</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Gundlach</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2009</year>). <article-title>The Sorghum bicolor genome and the diversification of grasses</article-title>. <source>Nature.</source> <volume>457</volume>, <fpage>551</fpage>&#x2013;<lpage>556</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature07723</pub-id>, PMID: <pub-id pub-id-type="pmid">19189423</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ringli</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Monitoring the outside: cell wall-sensing mechanisms</article-title>. <source>Plant Physiol.</source> <volume>153</volume>, <fpage>1445</fpage>&#x2013;<lpage>1452</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.110.154518</pub-id>, PMID: <pub-id pub-id-type="pmid">20508141</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roudier</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Fernandez</surname> <given-names>A. G.</given-names>
</name>
<name>
<surname>Fujita</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Himmelspach</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Borner</surname> <given-names>G. H.</given-names>
</name>
<name>
<surname>Schindelman</surname> <given-names>G.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>COBRA, an Arabidopsis extracellular glycosyl-phosphatidyl inositol-anchored protein, specifically controls highly anisotropic expansion through its involvement in cellulose microfibril orientation</article-title>. <source>Plant Cell.</source> <volume>17</volume>, <fpage>1749</fpage>&#x2013;<lpage>1763</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1105/tpc.105.031732</pub-id>, PMID: <pub-id pub-id-type="pmid">15849274</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roudier</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Schindelman</surname> <given-names>G.</given-names>
</name>
<name>
<surname>DeSalle</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Benfey</surname> <given-names>P. N.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>The COBRA family of putative GPI-anchored proteins in Arabidopsis. A new fellowship in expansion</article-title>. <source>Plant Physiol.</source> <volume>130</volume>, <fpage>538</fpage>&#x2013;<lpage>548</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.007468</pub-id>, PMID: <pub-id pub-id-type="pmid">12376623</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sajjad</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ahmad</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Riaz</surname> <given-names>M. W.</given-names>
</name>
<name>
<surname>Hussain</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Yasir</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>M. Z.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Recent genome resequencing paraded COBRA-Like gene family roles in abiotic stress and wood formation in Poplar</article-title>. <source>Front. Plant Sci.</source> <volume>14</volume>, <elocation-id>1242836</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2023.1242836</pub-id>, PMID: <pub-id pub-id-type="pmid">37780503</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schindelman</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Morikami</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Jung</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Baskin</surname> <given-names>T. I.</given-names>
</name>
<name>
<surname>Carpita</surname> <given-names>N. C.</given-names>
</name>
<name>
<surname>Derbyshire</surname> <given-names>P.</given-names>
</name>
<etal/>
</person-group>. (<year>2001</year>). <article-title>COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in Arabidopsis</article-title>. <source>Genes Dev.</source> <volume>15</volume>, <fpage>1115</fpage>&#x2013;<lpage>1127</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1101/gad.879101</pub-id>, PMID: <pub-id pub-id-type="pmid">11331607</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yuan</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>The YDA-MKK4/MKK5-MPK3/MPK6 cascade functions downstream of the RGF1-RGI ligand-receptor pair in regulating mitotic activity in root apical meristem</article-title>. <source>Mol. Plant</source> <volume>13</volume>, <fpage>1608</fpage>&#x2013;<lpage>1623</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2020.09.004</pub-id>, PMID: <pub-id pub-id-type="pmid">32916336</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Silva</surname> <given-names>T. N.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>J. B.</given-names>
</name>
<name>
<surname>Dahlberg</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Rhee</surname> <given-names>S. Y.</given-names>
</name>
<name>
<surname>Mortimer</surname> <given-names>J. C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Progress and challenges in sorghum biotechnology, a multipurpose feedstock for the bioeconomy</article-title>. <source>J. Exp. Bot.</source> <volume>73</volume>, <fpage>646</fpage>&#x2013;<lpage>664</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jxb/erab450</pub-id>, PMID: <pub-id pub-id-type="pmid">34644381</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sindhu</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Langewisch</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Olek</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Multani</surname> <given-names>D. S.</given-names>
</name>
<name>
<surname>McCann</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Vermerris</surname> <given-names>W.</given-names>
</name>
<etal/>
</person-group>. (<year>2007</year>). <article-title>Maize Brittle stalk2 encodes a COBRA-like protein expressed in early organ development but required for tissue flexibility at maturity</article-title>. <source>Plant Physiol.</source> <volume>145</volume>, <fpage>1444</fpage>&#x2013;<lpage>1459</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.107.102582</pub-id>, PMID: <pub-id pub-id-type="pmid">17932309</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xiong</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Natural variation of DROT1 confers drought adaptation in upland rice</article-title>. <source>Nat. Commun.</source> <volume>13</volume>, <fpage>4265</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-022-31844-w</pub-id>, PMID: <pub-id pub-id-type="pmid">35871266</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>DeBarry</surname> <given-names>J. D.</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity</article-title>. <source>Nucleic Acids Res.</source> <volume>40</volume>, <fpage>e49</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkr1293</pub-id>, PMID: <pub-id pub-id-type="pmid">22217600</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname> <given-names>B. P.</given-names>
</name>
<name>
<surname>Duan</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Gong</surname> <given-names>L. P.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>The OsDIR55 gene increases salt tolerance by altering the root diffusion barrier</article-title>. <source>Plant J.</source> <volume>118</volume>, <fpage>1550</fpage>&#x2013;<lpage>1568</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/tpj.16696</pub-id>, PMID: <pub-id pub-id-type="pmid">38412303</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Shan</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The MicroRNA397a-LACCASE17 module regulates lignin biosynthesis in Medicago ruthenica (L.)</article-title>. <source>Front. Plant Sci.</source> <volume>13</volume>, <elocation-id>978515</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2022.978515</pub-id>, PMID: <pub-id pub-id-type="pmid">36061772</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Luan</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xia</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Tang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tao</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Lignin provides mechanical support to herbaceous peony (Paeonia lactiflora Pall.) stems</article-title>. <source>Hortic. Res.</source> <volume>7</volume>, <fpage>213</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41438-020-00451-5</pub-id>, PMID: <pub-id pub-id-type="pmid">33372177</pub-id></citation></ref>
</ref-list>
</back>
</article>