Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes.

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Mining genic resources regulating nitrogen-use efficiency based on integrative biological analyses and their breeding applications in maize and other crops.

Nitrogen (N) is an essential factor for limiting crop yields, and cultivation of crops with low nitrogen-use efficiency (NUE) exhibits increasing environmental and ecological risks. Hence, it is crucial to mine valuable NUE improvement genes, which is very important to develop and breed new crop varieties with high NUE in sustainable agriculture system. Quantitative trait locus (QTL) and genome-wide association study (GWAS) analysis are the most common methods for dissecting genetic variations underlying complex traits. In addition, with the advancement of biotechnology, multi-omics technologies can be used to accelerate the process of exploring genetic variations. In this study, we integrate the substantial data of QTLs, quantitative trait nucleotides (QTNs) from GWAS, and multi-omics data including transcriptome, proteome, and metabolome and further analyze their interactions to predict some NUE-related candidate genes. We also provide the genic resources for NUE improvement among maize, rice, wheat, and sorghum by homologous alignment and collinearity analysis. Furthermore, we propose to utilize the knowledge gained from classical cases to provide the frameworks for improving NUE and breeding N-efficient varieties through integrated genomics, systems biology, and modern breeding technologies.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock.

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

Transcriptome analysis reveals that PbMYB61 and PbMYB308 are involved in the regulation of lignin biosynthesis in pear fruit stone cells.

Pear fruit stone cells have thick walls and are formed by the secondary deposition of lignin in the primary cell wall of thin-walled cells. Their content and size seriously affect fruit characteristics related to edibility. To reveal the regulatory mechanism underlying stone cell formation during pear fruit development and to identify hub genes, we examined the stone cell and lignin contents of 30 'Shannongsu' pear flesh samples and analyzed the transcriptomes of 15 pear flesh samples collected at five developmental stages. On the basis of the RNA-seq data, 35 874 differentially expressed genes were detected. Additionally, two stone cell-related modules were identified according to a WGCNA. A total of 42 lignin-related structural genes were subsequently obtained. Furthermore, nine hub structural genes were identified in the lignin regulatory network. We also identified PbMYB61 and PbMYB308 as candidate transcriptional regulators of stone cell formation after analyzing co-expression networks and phylogenetic relationships. Finally, we experimentally validated and characterized the candidate transcription factors and revealed that PbMYB61 regulates stone cell lignin formation by binding to the AC element in the PbLAC1 promoter to upregulate expression. However, PbMYB308 negatively regulates stone cell lignin synthesis by binding to PbMYB61 to form a dimer that cannot activate PbLAC1 expression. In this study, we explored the lignin synthesis-related functions of MYB family members. The results presented herein are useful for elucidating the complex mechanisms underlying lignin biosynthesis during pear fruit stone cell development.

Genome-wide identification and molecular evolution of Dof gene family in Camellia oleifera.

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera.

Genome-wide identification and molecular evolution of Dof transcription factors in Cyperus esculentus.

Dof transcription factor family in Cyperus esculentus genome was identified and analyzed using bioinformatics. The analysis results revealed that C.esculentus genome contains 29 Dof genes (CesDof), all of which are located in the nucleus according to subcellular localization prediction. CesDof proteinrs have a range of 124 to 512 amino acids, with most being basic proteins. Their secondary structure was mainly irregular curl. The promoter sequence of CesDof genes contains cis-acting elements that respond to light, drought, hormones, low temperature, and circadian rhythm. Codon preference analysis showed that CesDof genes' codon preference ends in T/A. Collinearity analysis revealed that C.esculentus had three pairs of collinear CesDof genes. Additionally, there were 15 pairs of collinear genes between C.esculentus and Arabidopsis thaliana. The genetic relationship between C.esculentus and Rhynchospora pubera was found to be the closest. Phylogenetic tree analysis revealed that 29 CesDof genes of C.esculentus can be classified into 4 subgroups. Additionally, 144 miRNAs were predicted to target these CesDof genes. Furthermore, protein interaction analysis indicated that 15 Dof proteins in C.esculentus had interactions. The qRT-PCR verification results of drought stress and salt stress treatment experiments showed that most CesDof genes were involved in drought stress and salt stress responses, and the gene expression trends under drought stress and salt stress conditions were consistent. These results lay a theoretical foundation for further studying the molecular functions of Dof gene family in C.esculentus and its molecular mechanisms in regulating the life activities of C.esculentus.

Genome-wide identification, phylogenetic, structural and functional evolution of the core components of ABA signaling in plant species: a focus on rice.

MAIN CONCLUSION: A genome-wide analysis had identified 642 ABA core component genes from 20 plant species, which were further categorized into three distinct subfamilies. The gene structures and evolutionary relationships of these genes had been characterized. PP2C_1, PP2C_2, and SnRK2_1 had emerged as key players in mediating the ABA signaling transduction pathway, specifically in rice, in response to abiotic stresses. The plant hormone abscisic acid (ABA) is essential for growth, development, and stress response, relying on its core components, pyrabactin resistance, pyrabactin resistance-like, and the regulatory component of ABA receptor (PYR/PYL/RCAR), 2C protein phosphatase (PP2C), sucrose non-fermenting-1-related protein kinase 2 (SnRK2). However, there's a lack of research on their structural evolution and functional differentiation across plants. Our study analyzed the phylogenetic, gene structure, homology, and duplication evolution of this complex in 20 plant species. We found conserved patterns in copy number and homology across subfamilies. Segmental and tandem duplications drove the evolution of these genes, while whole-genome duplication (WGD) expanded PYR/PYL/RCAR and PP2C subfamilies, enhancing environmental adaptation. In rice and Arabidopsis, the PYR/PYL/RCAR, PP2C, and SnRK2 genes showed distinct tissue-specific expression and responded to various stresses. Notably, PP2C_1 and PP2C_2 interacted with SnRK2_1 and were crucial for ABA signaling in rice. These findings offered new insights into ABA signaling evolution, interactions, and integration in green plants, benefiting future research in agriculture, evolutionary biology, ecology, and environmental science.

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera.

MAIN CONCLUSION: Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

A systematic genome-wide analysis and screening of the NAC family genes related to wood formation in Cinnamomum camphora.

Many processes, such as growth, aging, and adaptation to abiotic stress, are regulated in plants by NAC transcription factors. In woody plants, NAC transcription factors acts as a primary switch that regulates secondary xylem development by activating various downstream transcription factors and modulating expression levels of genes involved in the synthesis of the secondary cell wall. Our team had previously sequenced the whole genome of the camphor tree (Cinnamomum camphora). Here, we performed a detailed analysis of the NAC gene family of C. camphora and examined its evolutionary history. The genomic sequences of 121 NAC genes of C. camphora were identified and classified into 20 subfamilies in 2 major classes based on the phylogenetic analysis and structural features. Expansion of the CcNAC gene family occurred mainly by fragment replication and was influenced by the purifying selection. By analyzing predicted interactions of the homologous AtNAC proteins, we identified five CcNACs that potentially regulate xylem development in C. camphora. RNA sequencing revealed distinct expression profiles of CcNACs in seven different plant tissues. Subcellular localization prediction revealed that 120, 3, and 2 CcNACs have biological functions in the nucleus, cytoplasm, and chloroplast, respectively. Furthermore, we verified expression patterns of five CcNACs (CcNAC012, CcNAC028, CcNAC055, CcNAC080, and CcNAC119) in various tissue types using qRT-PCR. Our results will facilitate further in-depth studies of the molecular mechanisms by which CcNAC transcription factors regulate wood formation and other processes in C. camphora.

Understanding the Consequences of Collinearity for Multilevel Models: The Importance of Disaggregation Across Levels.

In multilevel models, disaggregating predictors into level-specific parts (typically accomplished via centering) benefits parameter estimates and their interpretations. However, the importance of level-specificity has been sparsely addressed in multilevel literature concerning collinearity. In this study, we develop novel insights into the interactivity of centering and collinearity in multilevel models. After integrating the broad literatures on centering and collinearity, we review level-specific and conflated correlations in multilevel data. Next, by deriving formal relationships between predictor collinearity and multilevel model estimates, we demonstrate how the consequences of collinearity change across different centering specifications and identify data characteristics that may exacerbate or mitigate those consequences. We show that when all or some level-1 predictors are uncentered, slope estimates can be greatly biased by collinearity. Disaggregation of all predictors eliminates the possibility that fixed effect estimates will be biased due to collinearity alone; however, under some data conditions, collinearity is associated with biased standard errors and random effect (co)variance estimates. Finally, we illustrate the importance of disaggregation for diagnosing collinearity in multilevel data and provide recommendations for the use of level-specific collinearity diagnostics. Overall, the necessity of disaggregation for identifying and managing collinearity's consequences in multilevel models is clarified in novel ways.

Phylogeny of the plant receptor-like kinase (RLK) gene family and expression analysis of wheat RLK genes in response to biotic and abiotic stresses.

BACKGROUND: The receptor-like kinase (RLK) gene families in plants contains a large number of members. They are membrane proteins with an extracellular receptor domain and participate in biotic and abiotic stress responses. RESULTS: In this study, we identified RLKs in 15 representative plant genomes, including wheat, and classified them into 64 subfamilies by using four types of phylogenetic trees and HMM models. Conserved exon‒intron structures with conserved exon phases in the kinase domain were found in many RLK subfamilies from Physcomitrella patens to Triticum aestivum. Domain distributions of RLKs were also diagrammed. Collinearity events and tandem gene clusters suggested that polyploidization and tandem duplication events contributed to the member expansions of T. aestivum RLKs. Global expression pattern analysis was performed by using public transcriptome data. These analyses were involved in T. aestivum, Aegilops tauschii and Brachypodium distachyon RLKs under biotic and abiotic stresses. We also selected 9 RLKs to validate the transcriptome prediction by using qRT‒PCR under drought treatment and with Fusarium graminearum infection. The expression trends of these 9 wheat RLKs from public transcriptome data were consistent with the results of qRT‒PCR, indicating that they might be stress response genes under drought or F. graminearum treatments. CONCLUSION: In this study, we identified, classified, evolved, and expressed RLKs in wheat and related plants. Thus, our results will provide insights into the evolutionary history and molecular mechanisms of wheat RLKs.

Genome wide characterization of R2R3 MYB transcription factor from Apocynum venetum revealed potential stress tolerance and flavonoid biosynthesis genes.

MYB transcription factors are crucial in regulating stress tolerance and expression of major genes involved in flavonoid biosynthesis. The functions of MYBs is well explored in a number of plants, yet no study is reported in Apocynum venetum. We identified a total of 163 MYB candidates, that comprised of 101 (61.96%) R2R3, 6 3R, 1 4R and 55 1R. Syntenic analysis of A. venetum R2R3 (AvMYBs) showed highest orthologous pairs with Vitis vinifera MYBs followed by Arabidopsis thaliana among the four species evaluated. Thirty segmental duplications and 6 tandem duplications were obtained among AvMYB gene pairs signifying their role in the MYB gene family expansion. Nucleotide substitution analysis (Ka/Ks) showed the AvMYBs to be under the influence of strong purifying selection. Expression analysis of selected AvMYBs under low temperature and cadmium stresses resulted in the identification of AvMYB48, AvMYB97, AvMYB8, AvMYB4 as potential stress responsive genes and AvMYB10 and AvMYB11 in addition, proanthocyanidin biosynthesis regulatory genes which is consistent with their annotated homologues in Arabidopsis. Tissue specific expression profile analysis of the AvMYBs further supported the qPCR analysis result. MYBs with higher transcript levels in root, stem and leaf like AvMYB4 for example, was downregulated under the stresses and such with low transcript level such as AvMYB48 which had low transcript in the leaf was upregulated under both stresses. Transcriptome and phylogenetic analyses suggested AvMYB42 as a potential regulator of anthocyanin biosynthesis. Thus, this study provided valuable information on AvR2R3-MYB gene family with respect to stress tolerance and flavonoid biosynthesis.

The BIDIAP index: a clinical, analytical and ultrasonographic score for the diagnosis of acute appendicitis in children.

BACKGROUND: Pediatric acute appendicitis (PAA) continues to be a diagnostic challenge today. The diagnostic performance of classical indices is only moderate, especially in pediatric population. This study aimed to define a clinical, radiological and analytical index for the diagnosis of PAA. MATERIALS AND METHODS: This prospective study included 151 patients divided into two groups: (1) 53 patients with non-surgical abdominal pain (NSAP) and (2) 98 patients with a confirmed PAA. Sociodemographic and clinical characteristics were compared between groups using the Mann-Whitney U test and the Fisher exact test. To identify the predictors of PAA, we performed a multivariable logistic regression using a forward stepwise analysis and we assigned multiples of integer values to the selected variables. The diagnostic performance of the index was assessed by calculating the area under the receiver operating characteristic curve. Intra-cohort calibration was assessed with the Hosmer-Lemeshow test. RESULTS: We developed the BIDIAP index (BIomarkers for the DIagnosis of Appendicitis in Pediatrics), which included three variables that independently predicted higher odds of PAA: appendiceal caliber (≥ 6.9 mm), systemic immune-inflammation index (≥ 890) and peritoneal irritation, which scored 4, 3 and 2 points, respectively. Mean (SD) score of the participants was 2.38 (2.06) in group 1 and 7.89 (1.50) in group 2. The area under the ROC was 0.97 (95% CI 0.95-0.99). The cut-off point was established at 4 points, resulting in a sensitivity of 98.98% and a specificity of 77.78%. CONCLUSIONS: The BIDIAP index has an exceptional diagnostic performance in PAA. The importance of these results lies in its novelty and in the simplicity of the index. Although external validation will be necessary, initial results look promising.

Collinearity Issues in Autoregressive Models with Time-Varying Serially Dependent Covariates.

First-order autoregressive models are popular to assess the temporal dynamics of a univariate process. Researchers often extend these models to include time-varying covariates, such as contextual factors, to investigate how they moderate processes' dynamics. We demonstrate that doing so has implications for how well one can estimate the autoregressive and covariate effects, as serial dependence in the variables can imply predictor collinearity. This is a noteworthy contribution, since in current practice serial dependence in a time-varying covariate is rarely considered important. We first recapitulate the role of predictor collinearity for estimation precision in an ordinary least squares context, by discussing how it affects estimator variances, covariances and correlations. We then derive a general formula detailing how predictor collinearity in first-order autoregressive models is impacted by serial dependence in the covariate. We provide a simulation study to illustrate the implications of the formula for different types of covariates. The simulation results highlight when the collinearity issue becomes severe enough to hamper interpretation of the effects. We also show that the effect estimates can be biased in small samples (i.e., 50 time points). Implications for study design, the use of time as a predictor, and related model variants are discussed.

Genome-Wide Identification of MADS-Box Genes in Taraxacum kok-saghyz and Taraxacum mongolicum: Evolutionary Mechanisms, Conserved Functions and New Functions Related to Natural Rubber Yield Formation.

MADS-box transcription regulators play important roles in plant growth and development. However, very few MADS-box genes have been isolated in the genus Taraxacum, which consists of more than 3000 species. To explore their functions in the promising natural rubber (NR)-producing plant Taraxacum kok-saghyz (TKS), MADS-box genes were identified in the genome of TKS and the related species Taraxacum mongolicum (TM; non-NR-producing) via genome-wide screening. In total, 66 TkMADSs and 59 TmMADSs were identified in the TKS and TM genomes, respectively. From diploid TKS to triploid TM, the total number of MADS-box genes did not increase, but expansion occurred in specific subfamilies. Between the two genomes, a total of 11 duplications, which promoted the expansion of MADS-box genes, were identified in the two species. TkMADS and TmMADS were highly conserved, and showed good collinearity. Furthermore, most TkMADS genes exhibiting tissue-specific expression patterns, especially genes associated with the ABCDE model, were preferentially expressed in the flowers, suggesting their conserved and dominant functions in flower development in TKS. Moreover, by comparing the transcriptomes of different TKS lines, we identified 25 TkMADSs related to biomass formation and 4 TkMADSs related to NR content, which represented new targets for improving the NR yield of TKS.

Genome-Wide Characterization and Analysis of bHLH Transcription Factors Related to Anthocyanin Biosynthesis in Cinnamomum camphora ('Gantong 1').

Cinnamomum camphora is one of the most commonly used tree species in landscaping. Improving its ornamental traits, particularly bark and leaf colors, is one of the key breeding goals. The basic helix-loop-helix (bHLH) transcription factors (TFs) are crucial in controlling anthocyanin biosynthesis in many plants. However, their role in C. camphora remains largely unknown. In this study, we identified 150 bHLH TFs (CcbHLHs) using natural mutant C. camphora 'Gantong 1', which has unusual bark and leaf colors. Phylogenetic analysis revealed that 150 CcbHLHs were divided into 26 subfamilies which shared similar gene structures and conserved motifs. According to the protein homology analysis, we identified four candidate CcbHLHs that were highly conserved compared to the TT8 protein in A. thaliana. These TFs are potentially involved in anthocyanin biosynthesis in C. camphora. RNA-seq analysis revealed specific expression patterns of CcbHLHs in different tissue types. Furthermore, we verified expression patterns of seven CcbHLHs (CcbHLH001, CcbHLH015, CcbHLH017, CcbHLH022, CcbHLH101, CcbHLH118, and CcbHLH134) in various tissue types at different growth stages using qRT-PCR. This study opens a new avenue for subsequent research on anthocyanin biosynthesis regulated by CcbHLH TFs in C. camphora.

Optical maps refine the bread wheat Triticum aestivum cv. Chinese Spring genome assembly.

Until recently, achieving a reference-quality genome sequence for bread wheat was long thought beyond the limits of genome sequencing and assembly technology, primarily due to the large genome size and > 80% repetitive sequence content. The release of the chromosome scale 14.5-Gb IWGSC RefSeq v1.0 genome sequence of bread wheat cv. Chinese Spring (CS) was, therefore, a milestone. Here, we used a direct label and stain (DLS) optical map of the CS genome together with a prior nick, label, repair and stain (NLRS) optical map, and sequence contigs assembled with Pacific Biosciences long reads, to refine the v1.0 assembly. Inconsistencies between the sequence and maps were reconciled and gaps were closed. Gap filling and anchoring of 279 unplaced scaffolds increased the total length of pseudomolecules by 168 Mb (excluding Ns). Positions and orientations were corrected for 233 and 354 scaffolds, respectively, representing 10% of the genome sequence. The accuracy of the remaining 90% of the assembly was validated. As a result of the increased contiguity, the numbers of transposable elements (TEs) and intact TEs have increased in IWGSC RefSeq v2.1 compared with v1.0. In total, 98% of the gene models identified in v1.0 were mapped onto this new assembly through development of a dedicated approach implemented in the MAGAAT pipeline. The numbers of high-confidence genes on pseudomolecules have increased from 105 319 to 105 534. The reconciled assembly enhances the utility of the sequence for genetic mapping, comparative genomics, gene annotation and isolation, and more general studies on the biology of wheat.

Genome-wide investigation of the ZF-HD gene family in two varieties of alfalfa (Medicago sativa L.) and its expression pattern under alkaline stress.

BACKGROUND: Zinc finger homeodomain (ZHD) protein is a plant-specific transcription factor and a potential regulator of phosphoenolpyruvate carboxylase (PEPCase)-coding genes, and it also participates in plant growth regulation and abiotic stress responses. To study the function of MsZF-HD genes in the alkaline stress response, this paper assessed biological information and performed transcriptome analysis of the MsZF-HD gene family by using the genomes of two different varieties of alfalfa (XinJiangDa Ye and Zhongmu No. 1). RESULTS: In total, 49 and 11 MsZF-HD genes were identified in the two different varieties respectively, including the alleles of XinJiangDa Ye. According to their phylogenetic relationships, the 60 MsZF-HD genes were divided into 5 ZHD subfamilies and 1 MIF subfamily. A total of 88.3% of MsZF-HD genes do not contain introns and are unevenly distributed among the 6 chromosomes of alfalfa. A collinearity analysis indicated that 26 genes of XinJiangDa Ye have no orthologous genes in Zhongmu No. 1, although these genes (such as ZHD-X1-2, ZHD-X3-2 and ZHD-X4-2) have homologous genes in Arabidopsis thaliana, Medicago truncatula and Glycine max. Through RNA-seq and qRT-PCR verification, it was found that MsZF-HD genes are downregulated to participate in the alkaline stress response. CONCLUSION: The results of this study may lay the foundation for the cloning and functional study of MsZF-HD genes and provide a theoretical basis for revealing the difference between XinJiangDa Ye and Zhongmu No. 1 at the genome level.

Genome-wide expression analysis of carboxylesterase (CXE) gene family implies GBCXE49 functional responding to alkaline stress in cotton.

BACKGROUND: Carboxylesterase (CXE) is a type of hydrolase with α/ß sheet hydrolase activity widely found in animals, plants and microorganisms, which plays an important role in plant growth, development and resistance to stress. RESULTS: A total of 72, 74, 39, 38 CXE genes were identified in Gossypium barbadense, Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively. The gene structure and expression pattern were analyzed. The GBCXE genes were divided into 6 subgroups, and the chromosome distribution of members of the family were mapped. Analysis of promoter cis-acting elements showed that most GBCXE genes contain cis-elements related to plant hormones (GA, IAA) or abiotic stress. These 6 genes we screened out were expressed in the root, stem and leaf tissues. Combined with the heat map, GBCXE49 gene was selected for subcellular locate and confirmed that the protein was expressed in the cytoplasm. CONCLUSIONS: The collinearity analysis of the CXE genes of the four cotton species in this family indicated that tandem replication played an indispensable role in the evolution of the CXE gene family. The expression patterns of GBCXE gene under different stress treatments indicated that GBCXE gene may significantly participate in the response to salt and alkaline stress through different mechanisms. Through the virus-induced gene silencing technology (VIGS), it was speculated that GBCXE49 gene was involved in the response to alkaline stress in G. barbadense.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern.

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.