NCBI's Conserved Domain Database and Tools for Protein Domain Analysis

Characterization and genomic analysis of the lytic bacteriophage vB_EclM_HK6 as a potential approach to biocontrol the spread of Enterobacter cloacae contaminating food

BACKGROUND

Increased prevalence of Enterobacter cloacae within food products underscores food as an underexplored reservoir for antibiotic resistance, thus requiring particular intervention. Bacteriophages have been explored as a promising approach for controlling bacterial growth in different matrices. Moreover, their specific interaction and self-replication, put them apart from traditional methods for controlling bacteria in different matrices.

METHODS

Sixteen Enterobacter cloacae strains were recovered from raw chicken. These strains were used to isolate bacteriophages using enrichment protocol. The broad-spectrum bacteriophage was evaluated in terms of thermal, pH, shearing stress and storge. Moreover, its infection kinetics, in vitro antibacterial activity, cytotoxicity were also assessed. Genomic sequencing was performed to exclude any potential virulence or resistance genes. Finally, the capability of the isolated phages to control bacterial growth in different chicken samples was assessed alone and in combination with sodium nitrite.

RESULTS

The lytic bacteriophage vB_EclM_HK6 was isolated and showed the broadest spectrum being able to infect 8/16 E. cloacae strains with a lytic activity against its host strain, E. cloacae EC21, as low as MOI of 10-6. The phage displays a latent period of 10 min and burst size of 115 ± 44 and resistance frequency of 5.7 × 10-4 ± 3.0 × 10-4. Stability assessment revealed a thermal tolerance up to 60 ˚C, wide range pH stability (3-10) and the ability to withstand shearing stress up to 250 rpm. HK6 shows no cytotoxicity against oral epithelial cells up to 1012 PFU/ml. Genomic analysis revealed a Strabovirus with total size of 177,845 bp that is free from known resistance and virulence genes. Finally, HK6 pretreatment of raw chicken, chicken nuggets and ready-made cheese salad shows a reduced bacterial count up to 4.6, 2.96 and 2.81 log-units, respectively. Moreover, combing HK6 with sodium nitrite further improved the antibacterial activity in both raw chicken and chicken nuggets without significant enhancement in case of cheese salad.

CONCLUSION

Enterobacter bacteriophage vB_EclM_HK6 presents a safe and effective approach for controlling E. cloacae contaminating stored chicken food samples. Moreover, they could be combined with a reduced concentrations of sodium nitrite to improve the killing capacity.

Genome-wide identification of the endonuclease family genes implicates potential roles of TaENDO23 in drought-stressed response and grain development in wheat

BACKGROUND

Endonucleases play a crucial role in plant growth and stress response by breaking down nuclear DNA. However, the specific members and biological functions of the endonuclease encoding genes in wheat remain to be determined.

RESULTS

In this study, we identified a total of 26 TaENDO family genes at the wheat genome-wide level. These genes were located on chromosomes 2 A, 2B, 2D, 3 A, 3B, and 3D and classified into four groups, each sharing similar gene structures and conserved motifs. Furthermore, we identified diverse stress-response and growth-related cis-elements in the promoter of TaENDO genes, which were broadly expressed in different organs, and several TaENDO genes were significantly induced under drought and salt stresses. We further examined the biological function of TaENDO23 gene since it was rapidly induced under drought stress and exhibited high expression in spikes and grains. Subcellular localization analysis revealed that TaENDO23 was localized in the cytoplasm of wheat protoplasts. qRT-PCR results indicated that the expression of TaENDO23 increased under PEG6000 and abscisic acid treatments, but decreased under NaCl treatment. TaENDO23 mainly expressed in leaves and spikes. A kompetitive allele-specific PCR (KASP) marker was developed to identify single nucleotide polymorphisms in TaENDO23 gene in 256 wheat accessions. The alleles with TaENDO23-HapI haplotypes had higher grain weight and size compared to TaENDO23-HapII. The geographical and annual frequency distributions of the two TaENDO23 haplotypes revealed that the elite haplotype TaENDO23-HapI was positively selected in the wheat breeding process.

CONCLUSION

We systematically analyzed the evolutionary relationships, gene structure characteristics, and expression patterns of TaENDO genes in wheat. The expression of TaENDO23, in particular, was induced under drought stress, mainly expressed in the leaves and grains. The KASP marker of TaENDO23 gene successfully distinguished between the wheat accessions, revealing TaENDO23-HapI as the elite haplotype associated with improved grain weight and size. These findings provide insights into the evolution and characteristics of TaENDO genes at the genome-wide level in wheat, laying the foundation for further biological analysis of TaENDO23 gene, especially in response to drought stress and grain development.

MOSCA 2.0: A bioinformatics framework for metagenomics, metatranscriptomics and metaproteomics data analysis and visualization

The analysis of meta-omics data requires the utilization of several bioinformatics tools and proficiency in informatics. The integration of multiple meta-omics data is even more challenging, and the outputs of existing bioinformatics solutions are not always easy to interpret. Here, we present a meta-omics bioinformatics pipeline, Meta-Omics Software for Community Analysis (MOSCA), which aims to overcome these limitations. MOSCA was initially developed for analysing metagenomics (MG) and metatranscriptomics (MT) data. Now, it also performs MG and metaproteomics (MP) integrated analysis, and MG/MT analysis was upgraded with an additional iterative binning step, metabolic pathways mapping, and several improvements regarding functional annotation and data visualization. MOSCA handles raw sequencing data and mass spectra and performs pre-processing, assembly, annotation, binning and differential gene/protein expression analysis. MOSCA shows taxonomic and functional analysis in large tables, performs metabolic pathways mapping, generates Krona plots and shows gene/protein expression results in heatmaps, improving omics data visualization. MOSCA is easily run from a single command while also providing a web interface (MOSGUITO). Relevant features include an extensive set of customization options, allowing tailored analyses to suit specific research objectives, and the ability to restart the pipeline from intermediary checkpoints using alternative configurations. Two case studies showcased MOSCA results, giving a complete view of the anaerobic microbial communities from anaerobic digesters and insights on the role of specific microorganisms. MOSCA represents a pivotal advancement in meta-omics research, offering an intuitive, comprehensive, and versatile solution for researchers seeking to unravel the intricate tapestry of microbial communities.

Genomic, LC-MS, and FTIR Analysis of Plant Probiotic Potential of Bacillus albus for Managing Xanthomonas oryzae via Different Modes of Application in Rice (Oryza sativa L.)

Xanthomonas oryzae causes tremendous damage in rice plants (Oryza sativa L). Therefore, this study is focused on siderophore-producing Bacillus albus (CWTS 10) for managing BLB disease caused by X. oryzae. Both B. albus and its crude siderophore (methanolic and diethyl ether) extracts inhibited X. oryzae (10-12 mm). Fourier transform infrared spectroscopy (FTIR) analysis of the extracts indicated the presence of catecholate siderophore functional groups. Liquid chromatography-mass spectrometry (LC-MS) analysis revealed the presence of antimicrobial compounds such as 2-deoxystreptamine, miserotoxin, fumitremorgin C, pipercide, pipernonaline, gingerone A, and deoxyvasicinone. Complete genome sequencing revealed the gene clusters for antibiotic, siderophore, antibacterial, antifungal, and secondary metabolite production. An in vivo study revealed that bacteria (CWTS 10) and their siderophore extracts effectively inhibited X. oryzae. The mode of application of bacterial or siderophore extracts in terms of DI and DSI percentage was as follows: soak method > inoculation method > spray method. In addition to providing enhanced antagonistic activity, there was a significant increase in root and shoot length and weight (wet and dry) of treated plants compared to control plants challenged with X. oryzae. Thus, the results clearly indicate that siderophore-producing B. albus and its siderophore extracts strongly inhibited X. oryzae. However, further field experiments are required before being formulated to protect rice crops from X. oryzae.

Analysis of the Rice Raffinose Synthase (OsRS) Gene Family and Haplotype Diversity

Based on the genome information of rice (Nipponbare), this study screened and identified six raffinose synthase (RS) genes and analyzed their physical and chemical properties, phylogenetic relationship, conserved domains, promoter cis-acting elements, and the function and genetic diversity of the gene-CDS-haplotype (gcHap). The results showed that these genes play key roles in abiotic stress response, such as OsRS5, whose expression in leaves changed significantly under high salt, drought, ABA, and MeJA treatments. In addition, the OsRS genes showed significant genetic variations in different rice populations. The main gcHaps of most OsRS loci had significant effects on key agronomic traits, and the frequency of these alleles varied significantly among different rice populations and subspecies. These findings provide direction for studying the RS gene family in other crops.

Genome-wide investigation of the PIF gene family in alfalfa (Medicago sativa L.) expression profiles during development and stress

BACKGROUND

Phytochrome-interacting factors (PIFs) plays an important role in plants as hubs for intracellular signaling regulation. The PIF gene family has been identified and characterized in many plants, but alfalfa (Medicago sativa L.), an important perennial high-quality legume forage, has not been reported on the PIF gene family.

RESULTS

In this study, we presented the identification and characterization of five MsPIF genes in alfalfa (Medicago sativa L.). Phylogenetic analysis indicated that PIFs from alfalfa and other four plant species could be divided into three groups supported by similar motif analysis. The collinearity analysis of the MsPIF gene family showed the presence of two gene pairs, and the collinearity analysis with AtPIFs showed three gene pairs, indicating that the evolutionary process of this family is relatively conservative. Analysis of cis-acting elements in promoter regions of MsPIF genes indicated that various elements were related to light, abiotic stress, and plant hormone responsiveness. Gene expression analyses demonstrated that MsPIFs were primarily expressed in the leaves and were induced by various abiotic stresses.

CONCLUSION

This study conducted genome-wide identification, evolution, synteny analysis, and expression analysis of the PIFs in alfalfa. Our study lays a foundation for the study of the biological functions of the PIF gene family and provides a useful reference for improving abiotic stress resistance in alfalfa.

Birth of protein folds and functions in the virome

The rapid evolution of viruses generates proteins that are essential for infectivity and replication but with unknown functions, due to extreme sequence divergence1. Here, using a database of 67,715 newly predicted protein structures from 4,463 eukaryotic viral species, we found that 62% of viral proteins are structurally distinct and lack homologues in the AlphaFold database2,3. Among the remaining 38% of viral proteins, many have non-viral structural analogues that revealed surprising similarities between human pathogens and their eukaryotic hosts. Structural comparisons suggested putative functions for up to 25% of unannotated viral proteins, including those with roles in the evasion of innate immunity. In particular, RNA ligase T-like phosphodiesterases were found to resemble phage-encoded proteins that hydrolyse the host immune-activating cyclic dinucleotides 3',3'- and 2',3'-cyclic GMP-AMP (cGAMP). Experimental analysis showed that RNA ligase T homologues encoded by avian poxviruses similarly hydrolyse cGAMP, showing that RNA ligase T-mediated targeting of cGAMP is an evolutionarily conserved mechanism of immune evasion that is present in both bacteriophage and eukaryotic viruses. Together, the viral protein structural database and analyses presented here afford new opportunities to identify mechanisms of virus-host interactions that are common across the virome.

Freshwater Mussel Viromes Increase Rapidly in Diversity and Abundance When Hosts Are Released from Captivity into the Wild

Freshwater mussels (order: Unionida) are highly imperiled globally and are increasingly the focus of captive propagation efforts to protect and restore wild populations. The Upper Tennessee River Basin (UTRB) in Virginia is a freshwater biodiversity hotspot hosting at least 45 of North America's ~300 species of freshwater mussels, including 21 threatened and endangered species listed under the U.S. Endangered Species Act. Recent studies have documented that viruses and other microbes have contributed to freshwater mussel population declines in the UTRB. We conducted a multi-year longitudinal study of captive-reared hatchery mussels released to restoration sites throughout the UTRB to evaluate their viromes and compare them to captive hatchery environments. We documented 681 viruses from 27 families. The hatchery mussels had significantly less viruses than those deployed to wild sites, with only 20 viruses unique to the hatchery mussels. After the hatchery mussels were released into the wild, their number of viruses initially spiked and then increased steadily over time, with 451 viruses in total unique to the mussels in the wild. We found Clinch densovirus 1 (CDNV-1), a virus previously associated with mass mortality events in the Clinch River, in all samples, but the wild site mussels consistently had significantly higher CDNV-1 levels than those held in the hatchery. Our data document substantial differences between the viruses in the mussels in the hatchery and wild environments and rapid virome shifts after the mussels are released to the wild sites. These findings indicate that mussel release programs might benefit from acclimatization periods or other measures to mitigate the potential negative effects of rapid exposure to infectious agents found in natural environments.

Genome-Wide Analysis and Expression Profiling of Soybean RbcS Family in Response to Plant Hormones and Functional Identification of GmRbcS8 in Soybean Mosaic Virus

Rubisco small subunit (RbcS), a core component with crucial effects on the structure and kinetic properties of the Rubisco enzyme, plays an important role in response to plant growth, development, and various stresses. Although Rbcs genes have been characterized in many plants, their muti-functions in soybeans remain elusive. In this study, a total of 11 GmRbcS genes were identified and subsequently divided into three subgroups based on a phylogenetic relationship. The evolutionary analysis revealed that whole-genome duplication has a profound effect on GmRbcSs. The cis-acting elements responsive to plant hormones, development, and stress-related were widely found in the promoter region. Expression patterns based on the RT-qPCR assay exhibited that GmRbcS genes are expressed in multiple tissues, and notably Glyma.19G046600 (GmRbcS8) exhibited the highest expression level compared to other members, especially in leaves. Moreover, differential expressions of GmRbcS genes were found to be significantly regulated by exogenous plant hormones, demonstrating their potential functions in diverse biology processes. Finally, the function of GmRbcS8 in enhancing soybean resistance to soybean mosaic virus (SMV) was further determined through the virus-induced gene silencing (VIGS) assay. All these findings establish a strong basis for further elucidating the biological functions of RbcS genes in soybeans.

Truncation of a novel C-terminal domain of a β-glucanase improves its thermal stability and specific activity

Enzymes that degrade β-glucan play important roles in various industries, including those related to brewing, animal feed, and health care. Csph16A, an endo-β-1,3(4)-glucanase encoded by a gene from the halotolerant, xerotolerant, and radiotrophic black fungus Cladosporium sphaerospermum, was cloned and expressed in Pichia pastoris. Two isoforms (Csph16A.1 and Csph16A.2) are produced, arising from differential glycosylation. The proteins were predicted to contain a catalytic Lam16A domain, along with a C-terminal domain (CTD) of unknown function which exhibits minimal secondary structure. Employing PCR-mediated gene truncation, the CTD of Csph16A was excised to assess its functional impact on the enzyme and determine potential alterations in biotechnologically relevant characteristics. The truncated mutant, Csph16A-ΔC, exhibited significantly enhanced thermal stability at 50°C, with D-values 14.8 and 23.5 times greater than those of Csph16A.1 and Csph16A.2, respectively. Moreover, Csph16A-ΔC demonstrated a 20%-25% increase in halotolerance at 1.25 and 1.5 M NaCl, respectively, compared to the full-length enzymes. Notably, specific activity against cereal β-glucan, lichenan, and curdlan was increased by up to 238%. This study represents the first characterization of a glucanase from the stress-tolerant fungus C. sphaerospermum and the first report of a halotolerant and engineered endo-β-1,3(4)-glucanase. Additionally, it sheds light on a group of endo-β-1,3(4)-glucanases from Antarctic rock-inhabiting black fungi harboring a Lam16A catalytic domain and a novel CTD of unknown function.

Evolution of ion channels in cetaceans: a natural experiment in the tree of life

Cetaceans represent a natural experiment within the tree of life in which a lineage changed from terrestrial to aquatic habitats. This shift involved phenotypic modifications, representing an opportunity to explore the genetic bases of phenotypic diversity. Among the different molecular systems that maintain cellular homeostasis, ion channels are crucial for the proper physiological functioning of all living species. This study aims to explore the evolution of ion channels during the evolutionary history of cetaceans. To do so, we created a bioinformatic pipeline to annotate the repertoire of ion channels in the genome of the species included in our sampling. Our main results show that cetaceans have, on average, fewer protein-coding genes and a higher percentage of annotated ion channels than non-cetacean mammals. Signals of positive selection were detected in ion channels related to the heart, locomotion, visual and neurological phenotypes. Interestingly, we predict that the NaV1.5 ion channel of most toothed whales (odontocetes) is sensitive to tetrodotoxin, similar to NaV1.7, given the presence of tyrosine instead of cysteine, in a specific position of the ion channel. Finally, the gene turnover rate of the cetacean crown group is more than three times faster than that of non-cetacean mammals.

SpSIZ1 from hyperaccumulator Sedum plumbizincicola orchestrates SpABI5 to fine-tune cadmium tolerance

Sedum plumbizincicola is a renowned hyperaccumulator of cadmium (Cd), possesses significant potential for eco-friendly phytoremediation of soil contaminated with Cd. Nevertheless, comprehension of the mechanisms underpinning its Cd stress response remains constrained, primarily due to the absence of a comprehensive genome sequence and an established genetic transformation system. In this study, we successfully identified a novel protein that specifically responds to Cd stress through early comparative iTRAQ proteome and transcriptome analyses under Cd stress conditions. To further investigate its structure, we employed AlphaFold, a powerful tool for protein structure prediction, and found that this newly identified protein shares a similar structure with Arabidopsis AtSIZ1. Therefore, we named it Sedum plumbizincicola SIZ1 (SpSIZ1). Our study revealed that SpSIZ1 plays a crucial role in positively regulating Cd tolerance through its coordination with SpABI5. Overexpression of SpSIZ1 significantly enhanced plant resistance to Cd stress and reduced Cd accumulation. Expression pattern analysis revealed higher levels of SpSIZ1 expression in roots compared to stems and leaves, with up-regulation under Cd stress induction. Importantly, overexpressing SpSIZ1 resulted in lower Cd translocation factors (Tfs) but maintained relatively constant Cd levels in roots under Cd stress, leading to enhanced Cd stress resistance in plants. Protein interaction analysis revealed that SpSIZ1 interacts with SpABI5, and the expression of genes responsive to abscisic acid (ABA) through SpABI5-dependent signaling was significantly up-regulated in SpSIZ1-overexpressing plants with Cd stress treatment. Collectively, our results illustrate that SpSIZ1 interacts with SpABI5, enhancing the expression of ABA downstream stress-related genes through SpABI5, thereby increasing Cd tolerance in plants.

Genome-wide identification, characterization, and evolutionary analysis of the barley TALE gene family and its expression profiles in response to exogenous hormones

Three-amino-loop-extension (TALE) family belongs to the homeobox gene superfamily and occurs widely in plants, playing a crucial role in regulating their growth and development. Currently, genome-wide analysis of the TALE family has been completed in many plants. However, the systematic identification and hormone response analysis of the TALE gene family in barley are still lacking. In this study, 21 TALE candidate genes were identified in barley, which can be divided into KNOX and BELL subfamilies. Barley TALE members in the same subfamily of the phylogenetic tree have analogically conserved motifs and gene structures, and segmental duplications are largely responsible for the expansion of the HvTALE family. Analysis of TALE orthologous and homologous gene pairs indicated that the HvTALE family has mainly undergone purifying selective pressure. Through spatial structure simulation, HvKNOX5-HvKNOX6 and HvKNOX5-HvBELL11 complexes are all formed through hydrogen bonding sites on both the KNOX2 and homeodomain (HD) domains of HvKNOX5, which may be essential for protein interactions among the HvTALE family members. Expression pattern analyses reveal the potential involvement of most HvTALE genes in responses to exogenous hormones. These results will lay the foundation for regulation and function analyses of the barley TALE gene family in plant growth and development by hormone regulation.

Genome-Wide Identification of the Whirly Gene Family and Its Potential Function in Low Phosphate Stress in Soybean (Glycine max)

The Whirly (WHY) gene family, functioning as transcription factors, plays an essential role in the regulation of plant metabolic responses, which has been demonstrated across multiple species. However, the WHY gene family and its functions in soybean remains unclear. In this paper, we conducted genome-wide screening and identification to characterize the WHY gene family. Seven WHY members were identified and randomly distributed across six chromosomes. The phylogenetic evolutionary tree of WHY genes in soybean and other species was divided into five clades. An in-depth analysis revealed that segmental duplications significantly contributed to the expansion of GmWHYs, and the GmWHY gene members may have experienced evolutionary pressure for purifying selection in soybeans. The analysis of promoter Cis-elements in GmWHYs suggested their potential significance in addressing diverse stress conditions. The expression patterns of GmWHYs exhibited tissue-specific variations throughout the different stages of soybean development. Additionally, six GmWHY genes exhibited different responses to low phosphate stress. These findings will provide a theoretical basis and valuable reference for the future exploration of WHY gene function.

Unveiling the mysteries: Functional insights into hypothetical proteins from Bacteroides fragilis 638R

Humans benefit from a vast community of microorganisms in their gastrointestinal tract, known as the gut microbiota, numbering in the tens of trillions. An imbalance in the gut microbiota known as dysbiosis, can lead to changes in the metabolite profile, elevating the levels of toxins like Bacteroides fragilis toxin (BFT), colibactin, and cytolethal distending toxin. These toxins are implicated in the process of oncogenesis. However, a significant portion of the Bacteroides fragilis genome consists of functionally uncharacterized and hypothetical proteins. This study delves into the functional characterization of hypothetical proteins (HPs) encoded by the Bacteroides fragilis genome, employing a systematic in silico approach. A total of 379 HPs were subjected to a BlastP homology search against the NCBI non-redundant protein sequence database, resulting in 162 HPs devoid of identity to known proteins. CDD-Blast identified 106 HPs with functional domains, which were then annotated using Pfam, InterPro, SUPERFAMILY, SCANPROSITE, SMART, and CATH. Physicochemical properties, such as molecular weight, isoelectric point, and stability indices, were assessed for 60 HPs whose functional domains were identified by at least three of the aforementioned bioinformatic tools. Subsequently, subcellular localization analysis was examined and the gene ontology analysis revealed diverse biological processes, cellular components, and molecular functions. Remarkably, E1WPR3 was identified as a virulent and essential gene among the HPs. This study presents a comprehensive exploration of B. fragilis HPs, shedding light on their potential roles and contributing to a deeper understanding of this organism's functional landscape.

Genomic-wide analysis reveals seven in absentia genes regulating grain development in wheat (Triticum aestivum L.)

Seven in absentia proteins, which contain a conserved SINA domain, are involved in regulating various aspects of wheat (Triticum aestivum L.) growth and development, especially in response to environmental stresses. However, it is unclear whether TaSINA family members are involved in regulating grain development until now. In this study, the expression pattern, genomic polymorphism, and relationship with grain-related traits were analyzed for all TaSINA members. Most of the TaSINA genes identified showed higher expression levels in young wheat spikes or grains than other organs. The genomic polymorphism analysis revealed that at least 62 TaSINA genes had different haplotypes, where the haplotypes of five genes were significantly correlated with grain-related traits. Kompetitive allele-specific PCR markers were developed to confirm the single nucleotide polymorphisms in TaSINA101 and TaSINA109 among the five selected genes in a set of 292 wheat accessions. The TaSINA101-Hap II and TaSINA109-Hap II haplotypes had higher grain weight and width compared to TaSINA101-Hap I and TaSINA109-Hap I in at least three environments, respectively. The qRT-PCR assays revealed that TaSINA101 was highly expressed in the palea shell, seed coat, and embryo in young wheat grains. The TaSINA101 protein was unevenly distributed in the nucleus when transiently expressed in the protoplast of wheat. Three homozygous TaSINA101 transgenic lines in rice (Oryza sativa L.) showed higher grain weight and size compared to the wild type. These findings provide valuable insight into the biological function and elite haplotype of TaSINA family genes in wheat grain development at a genomic-wide level.

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

CsPAT1, a GRAS transcription factor, promotes lignin accumulation by antagonistic interacting with CsWRKY13 in tea plants

Lignin is an important component of plant cell walls and plays crucial roles in the essential agronomic traits of tea quality and tenderness. However, the molecular mechanisms underlying the regulation of lignin biosynthesis in tea plants remain unclear. CsWRKY13 acts as a negative regulator of lignin biosynthesis in tea plants. In this study, we identified a GRAS transcription factor, phytochrome A signal transduction 1 (CsPAT1), that interacts with CsWRKY13. Silencing CsPAT1 expression in tea plants and heterologous overexpression in Arabidopsis demonstrated that CsPAT1 positively regulates lignin accumulation. Further investigation revealed that CsWRKY13 directly binds to the promoters of CsPAL and CsC4H and suppresses transcription of CsPAL and CsC4H. CsPAT1 indirectly affects the promoter activities of CsPAL and CsC4H by interacting with CsWRKY13, thereby facilitating lignin biosynthesis in tea plants. Compared with the expression of CsWRKY13 alone, the co-expression of CsPAT1 and CsWRKY13 in Oryza sativa significantly increased lignin biosynthesis. Conversely, compared with the expression of CsPAT1 alone, the co-expression of CsPAT1 and CsWRKY13 in O. sativa significantly reduced lignin accumulation. These results demonstrated the antagonistic regulation of the lignin biosynthesis pathway by CsPAT1 and CsWRKY13. These findings improve our understanding of lignin biosynthesis mechanisms in tea plants and provide insights into the role of the GRAS transcription factor family in lignin accumulation.

Unveiling the genetic architecture and transmission dynamics of a novel multidrug-resistant plasmid harboring blaNDM-5 in E. Coli ST167: implications for antibiotic resistance management

BACKGROUND

The emergence of multidrug-resistant (MDR) Escherichia coli strains poses significant challenges in clinical settings, particularly when these strains harbor New Delhi metallo-ß-lactamase (NDM) gene, which confer resistance to carbapenems, a critical class of last-resort antibiotics. This study investigates the genetic characteristics and implications of a novel blaNDM-5-carrying plasmid pNDM-5-0083 isolated from an E. coli strain GZ04-0083 from clinical specimen in Zhongshan, China.

RESULTS

Phenotypic and genotypic evaluations confirmed that the E. coli ST167 strain GZ04-0083 is a multidrug-resistant organism, showing resistance to diverse classes of antibiotics including ß-lactams, carbapenems, fluoroquinolones, aminoglycosides, and sulfonamides, while maintaining susceptibility to monobactams. Investigations involving S1 pulsed-field gel electrophoresis, Southern blot analysis, and conjugation experiments, alongside genomic sequencing, confirmed the presence of the blaNDM-5 gene within a 146-kb IncFIB plasmid pNDM-5-0083. This evidence underscores a significant risk for the horizontal transfer of resistance genes among bacterial populations. Detailed annotations of genetic elements-such as resistance genes, transposons, and insertion sequences-and comparative BLAST analyses with other blaNDM-5-carrying plasmids, revealed a unique architectural configuration in the pNDM-5-0083. The MDR region of this plasmid shares a conserved gene arrangement (repA-IS15DIV-blaNDM-5-bleMBL-IS91-suI2-aadA2-dfrA12) with three previously reported plasmids, indicating a potential for dynamic genetic recombination and evolution within the MDR region. Additionally, the integration of virulence factors, including the iro and sit gene clusters and enolase, into its genetic architecture poses further therapeutic challenges by enhancing the strain's pathogenicity through improved host tissue colonization, immune evasion, and increased infection severity.

CONCLUSIONS

The detailed identification and characterization of pNDM-5-0083 enhance our understanding of the mechanisms facilitating the spread of carbapenem resistance. This study illuminates the intricate interplay among various genetic elements within the novel blaNDM-5-carrying plasmid, which are crucial for the stability and mobility of resistance genes across bacterial populations. These insights highlight the urgent need for ongoing surveillance and the development of effective strategies to curb the proliferation of antibiotic resistance.

Genome sequencing and comparative genomics reveal insights into pathogenicity and evolution of Fusarium zanthoxyli, the causal agent of stem canker in prickly ash

BACKGROUND

Fusarium zanthoxyli is a destructive pathogen causing stem canker in prickly ash, an ecologically and economically important forest tree. However, the genome lack of F. zanthoxyli has hindered research on its interaction with prickly ash and the development of precise control strategies for stem canker.

RESULTS

In this study, we sequenced and annotated a relatively high-quality genome of F. zanthoxyli with a size of 43.39 Mb, encoding 11,316 putative genes. Pathogenicity-related factors are predicted, comprising 495 CAZymes, 217 effectors, 156 CYP450s, and 202 enzymes associated with secondary metabolism. Besides, a comparative genomics analysis revealed Fusarium and Colletotrichum diverged from a shared ancestor approximately 141.1 ~ 88.4 million years ago (MYA). Additionally, a phylogenomic investigation of 12 different phytopathogens within Fusarium indicated that F. zanthoxyli originated approximately 34.6 ~ 26.9 MYA, and events of gene expansion and contraction within them were also unveiled. Finally, utilizing conserved domain prediction, the results revealed that among the 59 unique genes, the most enriched domains were PnbA and ULP1. Among the 783 expanded genes, the most enriched domains were PKc_like kinases and those belonging to the APH_ChoK_Like family.

CONCLUSION

This study sheds light on the genetic basis of F. zanthoxyli's pathogenicity and evolution which provides valuable information for future research on its molecular interactions with prickly ash and the development of effective strategies to combat stem canker.

Functional domains of Acinetobacter bacteriophage tail fibers

A rapid increase in antimicrobial resistant bacterial infections around the world is causing a global health crisis. The Gram-negative bacterium Acinetobacter baumannii is categorized as a Priority 1 pathogen for research and development of new antimicrobials by the World Health Organization due to its numerous intrinsic antibiotic resistance mechanisms and ability to quickly acquire new resistance determinants. Specialized phage enzymes, called depolymerases, degrade the bacterial capsule polysaccharide layer and show therapeutic potential by sensitizing the bacterium to phages, select antibiotics, and serum killing. The functional domains responsible for the capsule degradation activity are often found in the tail fibers of select A. baumannii phages. To further explore the functional domains associated with depolymerase activity, tail-associated proteins of 71 sequenced and fully characterized phages were identified from published literature and analyzed for functional domains using InterProScan. Multisequence alignments and phylogenetic analyses were conducted on the domain groups and assessed in the context of noted halo formation or depolymerase characterization. Proteins derived from phages noted to have halo formation or a functional depolymerase, but no functional domain hits, were modeled with AlphaFold2 Multimer, and compared to other protein models using the DALI server. The domains associated with depolymerase function were pectin lyase-like (SSF51126), tailspike binding (cd20481), (Trans)glycosidases (SSF51445), and potentially SGNH hydrolases. These findings expand our knowledge on phage depolymerases, enabling researchers to better exploit these enzymes for therapeutic use in combating the antimicrobial resistance crisis.

Trace gas oxidation sustains energy needs of a thermophilic archaeon at suboptimal temperatures

Diverse aerobic bacteria use atmospheric hydrogen (H2) and carbon monoxide (CO) as energy sources to support growth and survival. Such trace gas oxidation is recognised as a globally significant process that serves as the main sink in the biogeochemical H2 cycle and sustains microbial biodiversity in oligotrophic ecosystems. However, it is unclear whether archaea can also use atmospheric H2. Here we show that a thermoacidophilic archaeon, Acidianus brierleyi (Thermoproteota), constitutively consumes H2 and CO to sub-atmospheric levels. Oxidation occurs across a wide range of temperatures (10 to 70 °C) and enhances ATP production during starvation-induced persistence under temperate conditions. The genome of A. brierleyi encodes a canonical CO dehydrogenase and four distinct [NiFe]-hydrogenases, which are differentially produced in response to electron donor and acceptor availability. Another archaeon, Metallosphaera sedula, can also oxidize atmospheric H2. Our results suggest that trace gas oxidation is a common trait of Sulfolobales archaea and may play a role in their survival and niche expansion, including during dispersal through temperate environments.

Isolation, characterisation and description of the roseoflavin producer Streptomyces berlinensis sp. nov

The Gram-positive bacteria Streptomyces davaonensis and Streptomyces cinnabarinus have been the only organisms known to produce roseoflavin, a riboflavin (vitamin B2) derived red antibiotic. Using a selective growth medium and a phenotypic screening, we were able to isolate a novel roseoflavin producer from a German soil sample. The isolation procedure was repeated twice, that is, the same strain could be isolated from the same location in Berlin 6 months and 12 months after its first isolation. Whole genome sequencing of the novel roseoflavin producer revealed an unusual chromosomal arrangement and the deposited genome sequence of the new isolate (G + C content of 71.47%) contains 897 genes per inverted terminal repeat, 6190 genes in the core and 107 genes located on an illegitimate terminal end. We identified the roseoflavin biosynthetic genes rosA, rosB and rosC and an unusually high number of riboflavin biosynthetic genes. Overexpression of rosA, rosB and rosC in Escherichia coli and enzyme assays confirmed their predicted functions in roseoflavin biosynthesis. A full taxonomic analysis revealed that the isolate represents a previously unknown Streptomyces species and we propose the name Streptomyces berlinensis sp. nov. for this roseoflavin producer.

Genome-wide identification of GATA transcription factors in tetraploid potato and expression analysis in differently colored potato flesh

The GATA gene family belongs to a kind of transcriptional regulatory protein featuring a zinc finger motif, which is essential for plant growth and development. However, the identification of the GATA gene family in tetraploid potato is still not performed. In the present research, a total of 88 GATA genes in the tetraploid potato C88.v1 genome were identified by bioinformatics methods. These StGATA genes had an uneven distribution on 44 chromosomes, and the corresponding StGATA proteins were divided into four subfamilies (I-IV) based on phylogenetic analysis. The cis-elements of StGATA genes were identified, including multiple cis-elements related to light-responsive and hormone-responsive. The collinearity analysis indicates that segmental duplication is a key driving force for the expansion of GATA gene family in tetraploid potato, and that the GATA gene families of tetraploid potato and Arabidopsis share a closer evolutionary relationship than rice. The transcript profiling analysis showed that all 88 StGATA genes had tissue-specific expression, indicating that the StGATA gene family members participate in the development of multiple potato tissues. The RNA-seq analysis was also performed on the tuber flesh of two potato varieties with different color, and 18 differentially expressed GATA transcription factor genes were screened, of which eight genes were validated through qRT-PCR. In this study, we identified and characterized StGATA transcription factors in tetraploid potato for the first time, and screened differentially expressed genes in potato flesh with different color. It provides a theoretical basis for further understanding the StGATA gene family and its function in anthocyanin biosynthesis.

Genome-Wide Analysis of VILLIN Gene Family Associated with Stress Responses in Cotton (Gossypium spp.)

The VILLIN (VLN) protein plays a crucial role in regulating the actin cytoskeleton, which is involved in numerous developmental processes, and is crucial for plant responses to both biotic and abiotic factors. Although various plants have been studied to understand the VLN gene family and its potential functions, there has been limited exploration of VLN genes in Gossypium and fiber crops. In the present study, we characterized 94 VLNs from Gossypium species and 101 VLNs from related higher plants such as Oryza sativa and Zea mays and some fungal, algal, and animal species. By combining these VLN sequences with other Gossypium spp., we classified the VLN gene family into three distinct groups, based on their phylogenetic relationships. A more in-depth examination of Gossypium hirsutum VLNs revealed that 14 GhVLNs were distributed across 12 of the 26 chromosomes. These genes exhibit specific structures and protein motifs corresponding to their respective groups. GhVLN promoters are enriched with cis-elements related to abiotic stress responses, hormonal signals, and developmental processes. Notably, a significant number of cis-elements were associated with the light responses. Additionally, our analysis of gene-expression patterns indicated that most GhVLNs were expressed in various tissues, with certain members exhibiting particularly high expression levels in sepals, stems, and tori, as well as in stress responses. The present study potentially provides fundamental insights into the VLN gene family and could serve as a valuable reference for further elucidating the diverse functions of VLN genes in cotton.

Genome-wide analysis of Glutathione peroxidase (GPX) gene family in Chickpea (Cicer arietinum L.) under salinity stress

Chickpea is the second most widely grown legume in the world. Its cultivation is highly affected by saline soils. Salt stress damages its all growth stages from germination to maturity. It has a huge genetic diversity containing adaptation loci that can help produce salt-tolerant cultivars. The glutathione peroxidase (GPX) gene family plays an important role in regulating plant response to abiotic stimuli and protects cells from oxidative damage. In current research, the role of GPX genes is studied for inducing salt tolerance in chickpea. This study identifies the GPX gene family in Cicer arietinum. In response to the NaCl stress, the gene expression profiles of CaGPX3 were examined using real-time qRT-PCR. The results of phylogenetic analysis show that CaGPX genes have an evolutionary relationship with monocots, dicots, chlorophytes, and angiosperms. Gene structure analysis showed that CaGPX3, CaGPX4, and CaGPX5 have six, CaGPX2 has five, and CaGPX1 contains 9 exons. According to the Ka and Ks analysis chickpea has one pair of duplicated genes of GPX and the duplication was tandem with negative (purifying) selection Ka < Ks (<1). In-silico gene expression analysis revealed that CaGPX3 is a salt stress-responsive gene among all other five GPX members in chickpea. The qRT-PCR results showed that the CaGPX3 gene expression was co-ordinately regulated under salt stress conditions, confirming CaGPX3's key involvement in salt tolerance.

Transcriptome analysis in Aegilops tauschii unravels further insights into genetic control of stripe rust resistance

MAIN CONCLUSION

The Aegilops tauschii resistant accession prevented the pathogen colonization by controlling the sugar flow and triggering the hypersensitive reaction. This study suggested that NBS-LRRs probably induce resistance through bHLH by controlling JA- and SA-dependent pathways. Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst) is one of wheat's most destructive fungal diseases that causes a severe yield reduction worldwide. The most effective and economically-friendly strategy to manage this disease is genetic resistance which can be achieved through deploying new and effective resistance genes. Aegilops tauschii, due to its small genome and co-evolution with Pst, can provide detailed information about underlying resistance mechanisms. Hence, we used RNA-sequencing approach to identify the transcriptome variations of two contrasting resistant and susceptible Ae. tauschii accessions in interaction with Pst and differentially expressed genes (DEGs) for resistance to stripe rust. Gene ontology, pathway analysis, and search for functional domains, transcription regulators, resistance genes, and protein-protein interactions were used to interpret the results. The genes encoding NBS-LRR, CC-NBS-kinase, and phenylalanine ammonia-lyase, basic helix-loop-helix (bHLH)-, basic-leucine zipper (bZIP)-, APETALA2 (AP2)-, auxin response factor (ARF)-, GATA-, and LSD-like transcription factors were up-regulated exclusively in the resistant accession. The key genes involved in response to salicylic acid, amino sugar and nucleotide sugar metabolism, and hypersensitive response contributed to plant defense against stripe rust. The activation of jasmonic acid biosynthesis and starch and sucrose metabolism pathways under Pst infection in the susceptible accession explained the colonization of the host. Overall, this study can fill the gaps in the literature on host-pathogen interaction and enrich the Ae. tauschii transcriptome sequence information. It also suggests candidate genes that could guide future breeding programs attempting to develop rust-resistant cultivars.

Genomic Analyses of Major SARS-CoV-2 Variants Predicting Multiple Regions of Pathogenic and Transmissive Importance

The rapid evolution of SARS-CoV-2 has fueled its global proliferation since its discovery in 2019, with several notable variants having been responsible for increases in cases of coronavirus disease 2019 (COVID-19). Analyses of codon bias and usage in these variants between phylogenetic clades or lineages may grant insights into the evolution of SARS-CoV-2 and identify target codons indicative of evolutionary or mutative trends that may prove useful in tracking or defending oneself against emerging strains. We processed a cohort of 120 SARS-CoV-2 genome sequences through a statistical and bioinformatic pipeline to identify codons presenting evidence of selective pressure as well as codon coevolution. We report the identification of two codon sites in the orf8 and N genes demonstrating such evidence with real-world impacts on pathogenicity and transmissivity.

Computational investigations into structure and function impact of novel mutations identified in targeted exons from ovarian cancer cell lines

The lack of sensitive and specific biomarkers for ovarian cancer leads to late stage diagnosis of the disease in a majority of the cases. Mutation accumulation is the basis for cancer progression, thus identifying mutations is an important step in the disease diagnosis. In the present study, a comprehensive analysis of fifteen Next Generation Sequencing samples from thirteen ovarian cancer cell lines was carried out for the identification of new mutations. The study revealed eight clinically significant novel mutations in six ovarian cancer oncogenes, viz. SMARCA4, ARID1A, PPP2R1A, CTNNB1, DICER1 and PIK3CA. In-depth computational analysis revealed that the mutations affected the structure of the proteins in terms of stability, solvent accessible surface area and molecular dynamics. Moreover, the mutations were present in functionally significant domains of the proteins, thereby adversely affecting the protein functionality. PPI network for SMARCA4, CTNNB1, DICER1, PIK3CA, PPP2R1A and ARID1A showed that these genes were involved in certain significant pathways affecting various hallmarks of cancer. For further validation, in vitro studies were performed that revealed hypermutability of the CTNNB1 gene. Through this study we have identified some key mutations and have analysed their structural and functional impact. The study establishes some key mutations, which can be potentially explored as biomarker and drug target.Communicated by Ramaswamy H. Sarma.

Helicovermis profundi gen. nov., sp. nov., a novel mesophilic, asporogenous bacterium within the Clostridia isolated from a deep-sea hydrothermal vent chimney

A novel mesophilic bacterial strain, designated S502T, was isolated from a deep-sea hydrothermal vent at Suiyo Seamount, Japan. Cells were Gram-positive, asporogenous, motile, and curved rods, measuring 1.6-5.6 µm in length. The strain was an obligate anaerobe that grew fermentatively on complex substrates such as yeast extract and Bacto peptone. Elemental sulfur stimulated the growth of the strain, and was reduced to hydrogen sulfide. The strain grew within a temperature range of 10-23 °C (optimum at 20 °C), pH range of 4.8-8.3 (optimum at 7.4), and a NaCl concentration range of 1.0-4.0% (w/v) (optimum at 3.0%, w/v). Phylogenetic analysis based on the 16S rRNA gene sequence revealed that the isolate was a member of the class Clostridia, with Fusibacter paucivorans strain SEBR 4211T (91.1% sequence identity) being its closest relative. The total size of the genome of the strain was 3.12 Mbp, and a G + C content was 28.2 mol%. The highest values for average nucleotide identity (ANI), average amino acid identity (AAI), and digital DNA-DNA hybridization (dDDH) value of strain S502T with relatives were 67.5% (with Marinisporobacter balticus strain 59.4MT), 51.5% (with M. balticus strain 59.4MT), and 40.9% (with Alkaliphilus serpentinus strain LacTT), respectively. Based on a combination of phylogenetic, genomic, and phenotypic characteristics, we propose strain S502T to represent a novel genus and species, Helicovermis profundi gen. nov., sp. nov., with the type strain S502T (= DSM 112048T = JCM 39167T).

Genomic insights into symbiosis and host adaptation of sponge-associated novel bacterium, Rossellomorea orangium sp. nov

Sponge-associated microorganisms play vital roles in marine sponge ecology. This study presents a genomic investigation of Rossellomorea sp. MCCB 382, isolated from Stelletta sp., reveals insights into its adaptations and symbiotic roles. Phylogenomic study and Overall Genomic Relatedness Index (OGRI) classify MCCB 382 as a novel species, Rossellomorea orangium sp. nov. The genome encodes numerous carbohydrate metabolism enzymes (CAZymes), likely aiding nutrient cycling in the sponge host. Unique eukaryotic-like protein domains hint at potential mechanisms of symbiosis. Defence mechanisms include CRISPR, restriction-modification systems, DNA phosphorothioation, toxin-antitoxin systems, and heavy metal and multidrug resistance genes, indicating adaptation to challenging marine environments. Unlike obligate mutualists, MCCB 382 shows no genome reduction. Furthermore, the presence of mobile genetic elements, horizontal gene transfer, and prophages suggest genetic versatility, implying flexible metabolic potential and capacity for rapid adaptation and symbiosis shifts. MCCB 382 possesses six biosynthetic gene clusters for secondary metabolites, including both type II and III polyketide synthases (PKS), terpenes, (NRPS), NRPS-independent-siderophore, and lassopeptide. Further genome mining using BiGScape revealed four distinct gene cluster families, T2PKS, NRPS-independent-siderophore, lasso peptide, and terpene, presenting opportunities for novel compound elucidation. Our study reveals a symbiotic lifestyle of MCCB 382 with the host sponge, highlighting symbiont factors that aid in establishing and sustaining this relationship. This is the pioneering genomic characterization of a novel Rossellomorea sp. within the sponge Stelletta sp. holobiont.

GTax: improving de novo transcriptome assembly by removing foreign RNA contamination

The cost and complexity of generating a complete reference genome means that many organisms lack an annotated reference. An alternative is to use a de novo reference transcriptome. This technology is cost-effective but is susceptible to off-target RNA contamination. In this manuscript, we present GTax, a taxonomy-structured database of genomic sequences that can be used with BLAST to detect and remove foreign contamination in RNA sequencing samples before assembly. In addition, we use a de novo transcriptome assembly of Solanum lycopersicum (tomato) to demonstrate that removing foreign contamination in sequencing samples reduces the number of assembled chimeric transcripts.

SingPro: a knowledge base providing single-cell proteomic data

Single-cell proteomics (SCP) has emerged as a powerful tool for detecting cellular heterogeneity, offering unprecedented insights into biological mechanisms that are masked in bulk cell populations. With the rapid advancements in AI-based time trajectory analysis and cell subpopulation identification, there exists a pressing need for a database that not only provides SCP raw data but also explicitly describes experimental details and protein expression profiles. However, no such database has been available yet. In this study, a database, entitled 'SingPro', specializing in single-cell proteomics was thus developed. It was unique in (a) systematically providing the SCP raw data for both mass spectrometry-based and flow cytometry-based studies and (b) explicitly describing experimental detail for SCP study and expression profile of any studied protein. Anticipating a robust interest from the research community, this database is poised to become an invaluable repository for OMICs-based biomedical studies. Access to SingPro is unrestricted and does not mandate a login at: http://idrblab.org/singpro/.

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

BACKGROUND

γ-glutamylcyclotransferase (GGCT), an enzyme to maintain glutathione homeostasis, plays a vital role in the response to plant growth and development as well as the adaptation to various stresses. Although the GGCT gene family analysis has been conducted in Arabidopsis and rice, the family genes have not yet been well identified and analyzed at the genome-wide level in wheat (Triticum aestivum L.).

RESULTS

In the present study, 20 TaGGCT genes were identified in the wheat genome and widely distributed on chromosomes 2A, 2B, 2D, 3A, 4A, 5A, 5B, 5D, 6A, 6B, 6D, 7A, 7B, and 7D. Phylogenetic and structural analyses showed that these TaGGCT genes could be classified into three subfamilies: ChaC, GGGACT, and GGCT-PS. They exhibited similar motif compositions and distribution patterns in the same subgroup. Gene duplication analysis suggested that the expansion of TaGGCT family genes was facilitated by segmental duplications and tandem repeats in the wheat evolutionary events. Identification of diverse cis-acting response elements in TaGGCT promoters indicated their potential fundamental roles in response to plant development and abiotic stresses. The analysis of transcriptome data combined with RT-qPCR results revealed that the TaGGCTs genes exhibited ubiquitous expression across plant organs, with highly expressed in roots, stems, and developing grains. Most TaGGCT genes were up-regulated after 6 h under 20% PEG6000 and ABA treatments. Association analysis revealed that two haplotypes of TaGGCT20 gene displayed significantly different Thousand-kernel weight (TKW), Kernel length (KL), and Kernel width (KW) in wheat. The geographical and annual distribution of the two haplotypes of TaGGCT20 gene further revealed that the frequency of the favorable haplotype TaGGCT20-Hap-I was positively selected in the historical breeding process of wheat.

CONCLUSION

This study investigated the genome-wide identification, structure, evolution, and expression analysis of TaGGCT genes in wheat. The motifs of TaGGCTs were highly conserved throughout the evolutionary history of wheat. Most TaGGCT genes were highly expressed in roots, stems, and developing grains, and involved in the response to drought stresses. Two haplotypes were developed in the TaGGCT20 gene, where TaGGCT20-Hap-I, as a favorable haplotype, was significantly associated with higher TKW, KL, and KW in wheat, suggesting that the haplotype is used as a function marker for the selection in grain yield in wheat breeding.

Rational peptide design for targeting cancer cell invasion

Cell invasion is an important process in cancer progression and recurrence. Invasion and implantation of cancer cells from their original place to other tissues, by disabling vital organs, challenges the treatment of cancer patients. Given the importance of the matter, many molecular treatments have been developed to inhibit cancer cell invasion. Because of their low production cost and ease of production, peptides are valuable therapeutic molecules for inhibiting cancer cell invasion. In recent years, advances in the field of computational biology have facilitated the design of anti-cancer peptides. In our investigation, using computational biology approaches such as evolutionary analysis, residue scanning, protein-peptide interaction analysis, molecular dynamics, and free energy analysis, our team designed a peptide library with about 100 000 candidates based on A6 (acetyl-KPSSPPEE-amino) sequence which is an anti-invasion peptide. During computational studies, two of the designed peptides that give the highest scores and showed the greatest sequence similarity to A6 were entered into the experimental analysis workflow for further analysis. In experimental analysis steps, the anti-metastatic potency and other therapeutic effects of designed peptides were evaluated using MTT assay, RT-qPCR, zymography analysis, and invasion assay. Our study disclosed that the IK1 (acetyl-RPSFPPEE-amino) peptide, like A6, has great potency to inhibit the invasion of cancer cells.

Novel, pathogenic insertion variant of GSDME associates with autosomal dominant hearing loss in a large Chinese pedigree

Nonsyndromic hearing loss (NSHL) is a genetically diverse, highly heterogeneous condition characterised by deafness, and Gasdermin E (GSDME) variants have been identified as directly inducing autosomal dominant NSHL. While many NSHL cases associated with GSDME involve the skipping of exon 8, there is another, less understood pathogenic insertion variant specifically found in Chinese pedigrees that causes deafness, known as autosomal dominant 5 (DFNA5) hearing loss. In this study, we recruited a large Chinese pedigree, conducted whole-exome and Sanger sequencing to serve as a comprehensive clinical examination, and extracted genomic DNA samples for co-segregation analysis of the members. Conservation and expression analyses for GSDME were also conducted. Our clinical examinations revealed an autosomal dominant phenotype of hearing loss in the family. Genetic analysis identified a novel insertion variant in GSDME exon 8 (GSDME: NM_004403.3: c.1113_1114insGGGGTGCAGCTTACAGGGTGGGTGT: p. P372fs*36). This variant is segregated with the deafness phenotype of this pedigree. The GSDME gene was highly conserved in the different species we analysed, and its mRNA expression was ubiquitously low in different human tissues. In conclusion, we have successfully identified a novel pathogenic insertion variant of GSDME in a Chinese pedigree that causes deafness, shedding light on the genetic basis of hearing loss within this specific family. Our findings expand the spectrum of known variants associated with GSDME-related deafness and may further support both the underlying gain-of-function mechanism and functional associations of GSDME hearing loss variants.

Real-time expression and in silico characterization of pea genes involved in salt and water-deficit stress

BACKGROUND

To tolerate salt and water-deficit stress, the plant adapts to the adverse environment by regulating its metabolism and expressing certain stress-induced metabolic pathways. This research analyzed the relative expression of four pea genes (P5CR, PAL1, SOD, and POX) in three pea varieties (Climax, Green grass, and Meteor) under different levels of salt and water-deficit stress.

METHODS AND RESULTS

The experiments on salt stress and water-deficit stress were carried out within greenhouse settings under controlled environment. The saturation percentage was employed to create artificial salinity conditions: Control without NaCl treatment, Treatment 1: 50 mM NaCl treatment, Treatment 2: 75 mM NaCl treatment, and Treatment 3: 100 mM NaCl treatment. Field capacity (FC) was used for the development of artificial water-deficit treatments in the pots, i.e., Treatment 1 (Control; water application 100% of FC), Treatment 2 (water application 75% of FC), and Treatment 3 (water application 50% of FC). Pea genes involved in biosynthetic pathways of proline, flavonoids, and enzymatic antioxidant enzymes including P5CR, PAL1, SOD, and POX were selected based on literature. Quantitative real-time PCR using cDNA as a template was used to analyze the gene expression. Pea genes were analyzed for phylogenetic analysis in closely related crops having similarity percent identity 80% and above. In silico characterization of selected proteins including the family classification was done by the NCBI CDD and INTERPRO online servers. Results from RT-qPCR analysis showed increased expression of P5CR, PAL1, and POX genes, while SOD gene expression decreased under both stresses. Climax exhibited superior stress tolerance with elevated expression of P5CR and PAL1, while Meteor showed better tolerance through increased POX expression. Phylogenetic analysis revealed common ancestry with other species like chickpea, red clover, mung bean, and barrel clover, suggesting the cross relationship among these plant species. Conserved domain analysis of respective proteins revealed that these proteins contain PLNO 2688, PLN02457, Cu-Zn Superoxide dismutase, and secretory peroxidase conserved domains. Furthermore, protein family classification indicated that the oxidation-reduction process is the most common chemical process involved in these stresses given to pea plant which validates the relationship of these proteins.

CONCLUSIONS

Salt and water-deficit stresses trigger distinct metabolic pathways, leading to the up-regulation of specific genes and the synthesis of corresponding proteins. These findings further emphasize the conservation of stress-tolerance-related genes and proteins across various plant species. This knowledge enhances our understanding of plant adaptation to stress and offers opportunities for developing strategies to improve stress resilience in crops, thereby addressing global food security challenges.

Identification of a functional peptide of a probiotic bacterium-derived protein for the sustained effect on preventing colitis

Several probiotic-derived factors have been identified as effectors of probiotics for exerting beneficial effects on the host. However, there is a paucity of studies to elucidate mechanisms of their functions. p40, a secretory protein, is originally isolated from a probiotic bacterium, Lactobacillus rhamnosus GG. Thus, this study aimed to apply structure-functional analysis to define the functional peptide of p40 that modulates the epigenetic program in intestinal epithelial cells for sustained prevention of colitis. In silico analysis revealed that p40 is composed of a signal peptide (1-28 residues) followed by a coiled-coil domain with uncharacterized function on the N-terminus, a linker region, and a β-sheet domain with high homology to CHAP on the C-terminus. Based on the p40 three-dimensional structure model, two recombinant p40 peptides were generated, p40N120 (28-120 residues) and p40N180 (28-180 residues) that contain first two and first three coiled coils, respectively. Compared to full-length p40 (p40F) and p40N180, p40N120 showed similar or higher effects on up-regulating expression of Setd1b (encoding a methyltransferase), promoting mono- and trimethylation of histone 3 on lysine 4 (H3K4me1/3), and enhancing Tgfb gene expression and protein production that leads to SMAD2 phosphorylation in human colonoids and a mouse colonic epithelial cell line. Furthermore, supplementation with p40F and p40N120 in early life increased H3K4me1, Tgfb expression and differentiation of regulatory T cells (Tregs) in the colon, and mitigated disruption of epithelial barrier and inflammation induced by DSS in adult mice. This study reveals the structural feature of p40 and identifies a functional peptide of p40 that could maintain intestinal homeostasis.

gjSOX9 Cloning, Expression, and Comparison with gjSOXs Family Members in Gekko japonicus

SOX9 plays a crucial role in the male reproductive system, brain, and kidneys. In this study, we firstly analyzed the complete cDNA sequence and expression patterns for SOX9 from Gekko japonicus SOX9 (gjSOX9), carried out bioinformatic analyses of physiochemical properties, structure, and phylogenetic evolution, and compared these with other members of the gjSOX family. The results indicate that gjSOX9 cDNA comprises 1895 bp with a 1482 bp ORF encoding 494aa. gjSOX9 was not only expressed in various adult tissues but also exhibited a special spatiotemporal expression pattern in gonad tissues. gjSOX9 was predicted to be a hydrophilic nucleoprotein with a characteristic HMG-Box harboring a newly identified unique sequence, "YKYQPRRR", only present in SOXE members. Among the 20 SOX9 orthologs, gjSOX9 shares the closest genetic relationships with Eublepharis macularius SOX9, Sphacrodactylus townsendi SOX9, and Hemicordylus capensis SOX9. gjSOX9 and gjSOX10 possessed identical physicochemical properties and subcellular locations and were tightly clustered with gjSOX8 in the SOXE group. Sixteen gjSOX family members were divided into six groups: SOXB, C, D, E, F, and H with gjSOX8, 9, and 10 in SOXE among 150 SOX homologs. Collectively, the available data in this study not only facilitate a deep exploration of the functions and molecular regulation mechanisms of the gjSOX9 and gjSOX families in G. japonicus but also contribute to basic research regarding the origin and evolution of SOX9 homologs or even sex-determination mode in reptiles.

Near-complete de novo assembly of Tricholoma bakamatsutake chromosomes revealed the structural divergence and differentiation of Tricholoma genomes

Tricholoma bakamatsutake, which is an edible ectomycorrhizal fungus associated with Fagaceae trees, may have diverged before the other species in Tricholoma section Caligata. We generated a highly contiguous whole-genome sequence for T. bakamatsutake SF-Tf05 isolated in an Oak (Quercus salicina) forest in Japan. The assembly of high-fidelity long reads, with a median read length of 12.3 kb, resulted in 13 chromosome-sized contigs comprising 142,068,211 bases with an average guanine and cytosine (GC) content of 43.94%. The 13 chromosomes were predicted to encode 11,060 genes. A contig (122,566 bases) presumably containing the whole circular mitochondrial genome was also recovered. The chromosome-wide comparison of T. bakamatsutake and Tricholoma matsutake (TMA_r1.0) indicated that the basic number of chromosomes (13) was conserved, but the structures of the corresponding chromosomes diverged, with multiple inversions and translocations. Gene conservation and cluster analyses revealed at least 3 phylogenetic clades in Tricholoma section Caligata. Specifically, all T. bakamatsutake strains belonged to the "bakamatsutake" clade, which is most proximal to the "caligatum" clade consisting of Tricholoma caligatum and Tricholoma fulvocastaneum. The constructed highly contiguous nearly telomere-to-telomere genome sequence of a T. bakamatsutake isolate will serve as a fundamental resource for future research on the evolution and differentiation of Tricholoma species.

Overexpressing CsPSY1 Gene of Tea Plant, Encoding a Phytoene Synthase, Improves α-Carotene and β-Carotene Contents in Carrot

Tea plants (Camellia sinensis (L.) O. Kuntze) belong to Theaceae family, in the section Thea. Tea plants are widely distributed in subtropical and tropical regions in the word. α-carotene and β-carotene in the tea leaves belong to carotenoids, which are associated with the aroma and color of the tea. Phytoene synthase (PSY) is a rate-limiting enzyme in carotenoids biosynthesis. We identified three CsPSY genes in 'Shuchazao', named CsPSY1, CsPSY2, and CsPSY3. Structural analysis of three CsPSY genes showed that CsPSY1 had a longer intro structure. The cis-acting elements of CsPSYs promoter were mainly associated with light-responsiveness, abiotic stress-responsiveness, and hormone-responsiveness. CsPSY1 exhibited expression in all tissues of the tea plants, whereas CsPSY2 and CsPSY3 were trace expression levels in all tissues. The positive expression of CsPSY1 under hormonal and abiotic stresses suggested its role in plant development and defense responses. The amino acid sequence of CsPSY1 was highly conserved in eight tea cultivars. The recombinant vector pCAMBIA1301-CsPSY1 was constructed to stabilize the overexpression of CsPSY1 in carrot. The contents of α-carotene and β-carotene in transgenic carrot callus were significantly increased. This study provides a foundational basis for further research on the function of CsPSYs and carotenoids accumulation in tea plants.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

Pyrofollis japonicus gen. nov. sp. nov., a novel member of the family Pyrodictiaceae isolated from the Iheya North hydrothermal field

A novel hyperthermophilic, heterotrophic archaeon, strain YC29T, was isolated from a deep-sea hydrothermal vent in the Mid-Okinawa Trough, Japan. Cells of strain YC29T were non-motile, irregular cocci with diameters of 1.2-3.0 µm. The strain was an obligatory fermentative anaerobe capable of growth on complex proteinaceous substrates. Growth was observed between 85 and 100 °C (optimum 90-95 °C), pH 4.9-6.4 (optimum 5.1), and in the presence of 1.4-4.0% (w/v) NaCl (optimum 3.0%). Inorganic carbon was required as a carbon source. Phylogenetic analysis based on the 16S rRNA gene sequence revealed that the isolate was a member of the family Pyrodictiaceae. The genome size was 2.02 Mbp with a G+C content of 49.4%. The maximum values for average nucleotide identity (ANI), average amino acid identity (AAI), and in silico DNA-DNA hybridization (dDDH) value of strain YC29T with relatives were 67.9% (with Pyrodictium abyssi strain AV2T), 61.1% (with Pyrodictium occultum strain PL-19T), and 33.8% (with Pyrolobus fumarii strain 1AT), respectively. Based on the phylogenetic, genomic, and phenotypic characteristics, we propose that strain YC29T represents a novel genus and species, Pyrofollis japonicus gen. nov., sp. (type strain YC29T = DSM 113394T = JCM 39171T).

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

Genome-wide identification and expression analysis of Hsp70 family genes in Cassava (Manihot esculenta Crantz)

Hsp70 proteins function as molecular chaperones, regulating various cellular processes in plants. In this study, a genome-wide analysis led to the identification of 22 Hsp70 (MeHsp70) genes in cassava. Phylogenetic relationship studies with other Malpighiales genomes (Populus trichocarpa, Ricinus communis and Salix purpurea) classified MeHsp70 proteins into eight groups (Ia, Ib, Ic, Id, Ie, If, IIa and IIb). Promoter analysis of MeHsp70 genes revealed the presence of tissue-specific, light, biotic and abiotic stress-responsive cis-regulatory elements showing their functional importance in cassava. Meta-analysis of publically available RNA-seq transcriptome datasets showed constitutive, tissue-specific, biotic and abiotic stress-specific expression patterns among MeHsp70s in cassava. Among 22 Hsp70, six MeHsp70s viz., MecHsp70-3, MecHsp70-6, MeBiP-1, MeBiP-2, MeBiP-3 and MecpHsp70-2 displayed constitutive expression, while three MecHsp70s were induced under both drought and cold stress conditions. Five MeHsp70s, MecHsp70-7, MecHsp70-11, MecHsp70-12, MecHsp70-13, and MecHsp70-14 were induced under drought stress conditions. We predicted that 19 MeHsp70 genes are under the regulation of 24 miRNAs. This comprehensive genome-wide analysis of the Hsp70 gene family in cassava provided valuable insights into their functional roles and identified various potential Hsp70 genes associated with stress tolerance and adaptation to environmental stimuli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03760-3.

Identification and Characterization of Novel Malto-Oligosaccharide-Forming Amylase AmyCf from Cystobacter sp. Strain CF23

Malto-oligosaccharides (MOSs) from starch conversion is advantageous for food and pharmaceutical applications. In this study, an efficient malto-oligosaccharide-forming α-amylase AmyCf was identified from myxobacter Cystobacter sp. strain CF23. AmyCf is composed of 417 amino acids with N-terminal 41 amino acids as the signal peptide, and conserved glycoside hydrolase family 13 (GH13) catalytic module and predicted C-terminal domain with β-sheet structure are also identified. Phylogenetic and functional analysis demonstrated that AmyCf is a novel member of GH13_6 subfamily. The special activity of AmyCf toward soluble starch and raw wheat starch is 9249 U/mg and 11 U/mg, respectively. AmyCf has broad substrate specificity toward different types of starches without requiring Ca2+. Under ideal circumstances of 60 °C and pH 7.0, AmyCf hydrolyzes gelatinized starch into maltose and maltotriose and maltotetraose as the main hydrolytic products with more than 80% purity, while maltose and maltotriose are mainly produced from the hydrolysis of raw wheat starch with more than 95% purity. The potential applicability of AmyCf in starch processing is highlighted by its capacity to convert gelatinized starch and raw starch granules into MOSs. This enzymatic conversion technique shows promise for the low-temperature enzymatic conversion of raw starch.

Genome-wide identification, evolution and expression pattern analysis of the GATA gene family in Sorghum bicolor

The GATA family of transcription factors is zinc finger DNA binding proteins involved in a variety of biological processes, including plant growth and development and response to biotic/abiotic stresses, and thus play an essential role in plant response to environmental changes. However, the GATA gene family of Sorghum (SbGATA) has not been systematically analyzed and reported yet. Herein, we used a variety of bioinformatics methods and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) to explore the evolution and function of the 33 SbGATA genes identified. These SbGATA genes, distributed on 10 chromosomes, are classified into four subfamilies (I-IV) containing one pair of tandem duplications and nine pairs of segment duplications, which are more closely related to the monocot Brachypodium distachyon and Oryza sativa GATA genes. The physicochemical properties of the SbGATAs are significantly different among the subfamilies, while the protein structure and conserved protein motifs are highly conserved in the subfamilies. In addition, the transcription of SbGATAs is tissue-specific during Sorghum growth and development, which allows for functional diversity in response to stress and hormones. Collectively, our study lays a theoretical foundation for an in-depth analysis of the functions, mechanisms and evolutionary relationships of SbGATA during plant growth and development.

Analysis of Raffinose Synthase Gene Family in Bread Wheat and Identification of Drought Resistance and Salt Tolerance Function of TaRS15-3B

Raffinose synthase (RS) plays a crucial role in plant growth and development, as well as in responses to biotic stresses and abiotic stresses, yet few studies have been conducted on its role in bread wheat. Therefore, in this study we screened and identified a family of bread wheat raffinose synthase genes based on bread wheat genome information and analyzed their physicochemical properties, phylogenetic evolutionary relationships, conserved structural domains, promoter cis-acting elements, and expression patterns. The BSMV-induced silencing of TaRS15-3B resulted in the bread wheat seedlings being susceptible to drought and salt stress and reduced the expression levels of stress-related and ROS-scavenging genes in bread wheat plants. This further affected the ability of bread wheat to cope with drought and salt stress. In conclusion, this study revealed that the RS gene family in bread wheat plays an important role in plant response to abiotic stresses and that the TaRS15-3B gene can improve the tolerance of transgenic bread wheat to drought and salt stresses, provide directions for the study of other RS gene families in bread wheat, and supply candidate genes for use in molecular breeding of bread wheat for stress resistance.

Identification, characterization, and validation of NBS-encoding genes in grass pea

Grass pea is a promising crop with the potential to provide food and fodder, but its genomics has not been adequately explored. Identifying genes for desirable traits, such as drought tolerance and disease resistance, is critical for improving the plant. Grass pea currently lacks known R-genes, including the nucleotide-binding site-leucine-rich repeat (NBS-LRR) gene family, which plays a key role in protecting the plant from biotic and abiotic stresses. In our study, we used the recently published grass pea genome and available transcriptomic data to identify 274 NBS-LRR genes. The evolutionary relationships between the classified genes on the reported plants and LsNBS revealed that 124 genes have TNL domains, while 150 genes have CNL domains. All genes contained exons, ranging from 1 to 7. Ten conserved motifs with lengths ranging from 16 to 30 amino acids were identified. We found TIR-domain-containing genes in 132 LsNBSs, with 63 TIR-1 and 69 TIR-2, and RX-CCLike in 84 LsNBSs. We also identified several popular motifs, including P-loop, Uup, kinase-GTPase, ABC, ChvD, CDC6, Rnase_H, Smc, CDC48, and SpoVK. According to the gene enrichment analysis, the identified genes undergo several biological processes such as plant defense, innate immunity, hydrolase activity, and DNA binding. In the upstream regions, 103 transcription factors were identified that govern the transcription of nearby genes affecting the plant excretion of salicylic acid, methyl jasmonate, ethylene, and abscisic acid. According to RNA-Seq expression analysis, 85% of the encoded genes have high expression levels. Nine LsNBS genes were selected for qPCR under salt stress conditions. The majority of the genes showed upregulation at 50 and 200 μM NaCl. However, LsNBS-D18, LsNBS-D204, and LsNBS-D180 showed reduced or drastic downregulation compared to their respective expression levels, providing further insights into the potential functions of LsNBSs under salt stress conditions. They provide valuable insights into the potential functions of LsNBSs under salt stress conditions. Our findings also shed light on the evolution and classification of NBS-LRR genes in legumes, highlighting the potential of grass pea. Further research could focus on the functional analysis of these genes, and their potential use in breeding programs to improve the salinity, drought, and disease resistance of this important crop.

Genome-wide investigation of NLP gene family members in alfalfa (Medicago sativa L.): evolution and expression profiles during development and stress

BACKGROUND

NIN-like protein (NLP) transcription factors (TFs) compose a plant-specific gene family whose members play vital roles in plant physiological processes, especially in the regulation of plant growth and the response to nitrate-nitrogen. However, no systematic identification or analysis of the NLP gene family has been reported in alfalfa. The recently completed whole-genome sequence of alfalfa has allowed us to investigate genome-wide characteristics and expression profiles.

RESULTS

53 MsNLP genes were identified from alfalfa and renamed according to their respective chromosome distributions. Phylogenetic analysis demonstrated that these MsNLPs can be classified into three groups on the basis of their conserved domains. Gene structure and protein motif analyses showed that closely clustered MsNLP genes were relatively conserved within each subgroup. Synteny analysis revealed four fragment duplication events of MsNLPs in alfalfa. The ratios of nonsynonymous (Ka) and synonymous (Ks) substitution rates of gene pairs indicated that the MsNLP genes underwent purifying selection during evolution. Examination of the expression patterns of different tissues revealed specific expression patterns of the MsNLP genes in the leaves, indicating that these genes are involved in plant functional development. Prediction of cis-acting regulatory elements and expression profiles further demonstrated that the MsNLP genes might play important roles in the response to abiotic stress and in phytohormone signal transduction processes.

CONCLUSION

This study represents the first genome-wide characterization of MsNLP in alfalfa. Most MsNLPs are expressed mainly in leaves and respond positively to abiotic stresses and hormonal treatments. These results provide a valuable resource for an improved understanding of the characteristics and biological roles of the MsNLP genes in alfalfa.