Protein identification and analysis tools in the ExPASy server

Exploring the roles of ZmARM gene family in maize development and abiotic stress response

Armadillo (ARM) was a gene family important to plants, with crucial roles in regulating plant growth, development, and stress responses. However, the properties and functions of ARM family members in maize had received limited attention. Therefore, this study employed bioinformatics methods to analyze the structure and evolution of ARM-repeat protein family members in maize. The maize (Zea mays L.) genome contains 56 ARM genes distributed over 10 chromosomes, and collinearity analysis indicated 12 pairs of linkage between them. Analysis of the physicochemical properties of ARM proteins showed that most of these proteins were acidic and hydrophilic. According to the number and evolutionary analysis of the ARM genes, the ARM genes in maize can be divided into eight subgroups, and the gene structure and conserved motifs showed similar compositions in each group. The findings shed light on the significant roles of 56 ZmARM domain genes in development and abiotic stress, particularly drought stress. RNA-Seq and qRT-PCR analysis revealed that drought stress exerts an influence on specific members of the ZmARM family, such as ZmARM4, ZmARM12, ZmARM34 and ZmARM36. The comprehensive profiling of these genes in the whole genome, combined with expression analysis, establishes a foundation for further exploration of plant gene function in the context of abiotic stress and reproductive development.

Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding

The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.

Immunoinformatic exploration of a multi-epitope-based peptide vaccine candidate targeting emerging variants of SARS-CoV-2

Many countries around the world are facing severe challenges due to the recently emerging variants of SARS-CoV-2. Over the last few months, scientists have been developing treatments, drugs, and vaccines to subdue the virus and prevent its transmission. In this context, a peptide-based vaccine construct containing pathogenic proteins of the virus known to elicit an immune response was constructed. An analysis of the spike protein-based epitopes allowed us to design an "epitope-based subunit vaccine" against coronavirus using the approaches of "reverse vaccinology" and "immunoinformatics." Computational experimentation and a systematic, comprehensive protocol were followed with an aim to develop and design a multi-epitope-based peptide (MEBP) vaccine candidate. Our study attempted to predict an MEBP vaccine by introducing mutations of SARS-CoV-2 (Delta, Lambda, Iota, Omicron, and Kappa) in Spike glycoprotein and predicting dual-purpose epitopes (B-cell and T-cell). This was followed by screening the selected epitopes based on antigenicity, allergenicity, and population coverage and constructing them into a vaccine by using linkers and adjuvants. The vaccine construct was analyzed for its physicochemical properties and secondary structure prediction, and a 3D structure was built, refined, and validated. Furthermore, the peptide-protein interaction of the vaccine construct with Toll-like receptor (TLR) molecules was performed. Immune profiling was performed to check the immune response. Codon optimization of the vaccine construct was performed to obtain the GC content before cloning it into the E. coli genome, facilitating its progression it into a vector. Finally, an in-silico simulation of the vaccine-protein complex was performed to comprehend its stability and conformational behavior.

Prediction of Antigenic Vaccine Peptide Candidates From BfmRS Associated With Biofilm Formation in Acinetobacter baumannii

INTRODUCTION

 A. baumannii is categorized as a priority pathogen due to its propensity for multi-drug resistance, exhibiting resistance against the last resort of antibiotics. It is also considered a potent nosocomial pathogen, so targeting the microbe using novel strategies would be the need of the hour. In this context, the in-silico computational approach would serve the best to design the possible epitope peptides, which may be further considered for the experimental trials for their immunological response.  Objective: To predict the immune-dominant epitope peptide candidates against the bfmR and bfmS proteins mediating the two-component system adaptation in the formation of biofilm in A. baumannii.

MATERIALS AND METHODS

11 different FASTA sequences of bfmR and bfmS from A. baumannii strains retrieved based on the blast-p similarity search tool were subjected to linear epitope B-cell epitope predictions under the IEDB B-cell epitope prediction server. Further analysis on antigenicity, allergenicity, and toxigenicity was achieved using the AntigenPro, Vaxijen, and AlgPred tools, with the physical and chemical properties evaluated using the Expasy Protparam server. Selection of the immunodominant peptides for T-cells was done through the databases under IEDB. The final assessment of protein-TLR2 interactions was done by MHC cluster servers.

RESULTS

Four peptide sequences (E1-E4) were predicted for B-cell dominance, with E1, E2, and E4 as probable antigens. All were soluble and non-toxigenic. E1 and E3 were considered non-allergens. GRAVY values were negative for all the peptides, indicating the protein to be hydrophilic in nature. Analysis of the T-cell epitopes was promising, with 100% conservancy for class-I HLA alleles, high interaction scores for similarity with TLR2, and more hydrogen bonds for E2, followed by other epitope peptides.

CONCLUSION

The promising four epitopes, as predicted for bfmR and bfmS in the present study, suggest their potent role as possible candidates for the design of vaccines targeting the TCS of A. baumannii, recommending further in vitro and in vivo experimental validation.

PepAnalyzer: predicting peptide properties using its sequence

Peptides are short linear molecules consisting of amino acids that play an essential role in most biological processes. They can treat diseases by working as a vaccine or antimicrobial agent and serves as a cancer molecule to deliver the drug to the target site for the treatment of cancer. They have the potential to solve the drawbacks of current medications and can be industrially produced in large quantities at low cost. However, poor chemical and physical stability, short circulating plasma half-life, and solubility are some issues that need solutions before they can be used as therapeutics. PepAnalyzer tool is a user-friendly tool that predicts 15 different properties such as binding potential, half-life, transmembrane patterns, test tube stability, charge, isoelectric point, molecular weights, and molar extinction coefficients only using the sequence. The tool is designed using BioPython utility and has even results with standard tools, such as Expasy, EBI, Genecorner, and Geneinfinity. The tool assists students, researchers, and the pharmaceutical sector. The PepAnalyzer tool's online platform is accessible at the link: http://www.iksmbrlabdu.in/peptool .

Design of a novel multi-epitope vaccine candidate against endometrial cancer using immunoinformatics and bioinformatics approaches

Endometrial cancer (EC) is one of the most common cancers of the female reproductive system. Multi-epitope vaccine may be a promising and effective strategy against EC. In this study, we designed a novel multi-epitope vaccine based on the antigenic proteins PRAME and TMPRSS4 using immunoinformatics and bioinformatics approaches. After a rigorous selection process, 14 cytotoxic T lymphocyte (CTL) epitopes, 6 helper T lymphocyte (HTL) epitopes, and 8 B cell epitopes (BCEs) were finally selected for vaccine construction. To enhance the immunogenicity of the vaccine candidate, the pan HLA DR-binding epitope was included in the vaccine design as an adjuvant. The final vaccine construct had 455 amino acids and a molecular weight of 49.8 kDa, and was predicted to cover 95.03% of the total world population. Docking analysis showed that there were 10 hydrogen bonds and 19 hydrogen bonds in the vaccine-HLA-A*02:01 and vaccine-HLA-DRB1*01:01 complexes, respectively, indicating that the vaccine has a good affinity to MHC molecules. This was further supported by molecular dynamics (MD) simulation. Immune simulation showed that the designed vaccine was able to induce higher levels of immune cell activity, with the secretion of numerous cytokines. The codon adaptation index (CAI) value and GC content of the optimised codon sequences of the vaccine were 0.986 and 54.43%, respectively, indicating that the vaccine has the potential to be highly expressed. The in silico analysis suggested that the designed vaccine may provide a novel therapeutic option for the individualised treatment of EC patients in the future.Communicated by Ramaswamy H. Sarma.

Vaccinomics-based next-generation multi-epitope chimeric vaccine models prediction against Leishmania tropica - a hierarchical subtractive proteomics and immunoinformatics approach

Leishmania tropica is a vector-borne parasitic protozoa that is the leading cause of leishmaniasis throughout the global tropics and subtropics. L. tropica is a multidrug-resistant parasite with a diverse set of serological, biochemical, and genomic features. There are currently no particular vaccines available to combat leishmaniasis. The present study prioritized potential vaccine candidate proteins of L. tropica using subtractive proteomics and vaccinomics approaches. These vaccine candidate proteins were downstream analyzed to predict B- and T-cell epitopes based on high antigenicity, non-allergenic, and non-toxic characteristics. The top-ranked overlapping MHC-I, MHC-II, and linear B-cell epitopes were prioritized for model vaccine designing. The lead epitopes were linked together by suitable linker sequences to design multi-epitope constructs. Immunogenic adjuvant sequences were incorporated at the N-terminus of the model vaccine constructs to enhance their immunological potential. Among different combinations of constructs, four vaccine designs were selected based on their physicochemical and immunological features. The tertiary structure models of the designed vaccine constructs were predicted and verified. The molecular docking and molecular dynamic (MD) simulation analyses indicated that the vaccine design V1 demonstrated robust and stable molecular interactions with toll-like receptor 4 (TLR4). The top-ranked vaccine construct model-IV demonstrated significant expressive capability in the E. coli expression system during in-silico restriction cloning analysis. The results of the present study are intriguing; nevertheless, experimental bioassays are required to validate the efficacy of the predicted model chimeric vaccine.

A previously unrecognized superfamily of macro-conotoxins includes an inhibitor of the sensory neuron calcium channel Cav2.3

Animal venom peptides represent valuable compounds for biomedical exploration. The venoms of marine cone snails constitute a particularly rich source of peptide toxins, known as conotoxins. Here, we identify the sequence of an unusually large conotoxin, Mu8.1, which defines a new class of conotoxins evolutionarily related to the well-known con-ikot-ikots and 2 additional conotoxin classes not previously described. The crystal structure of recombinant Mu8.1 displays a saposin-like fold and shows structural similarity with con-ikot-ikot. Functional studies demonstrate that Mu8.1 curtails calcium influx in defined classes of murine somatosensory dorsal root ganglion (DRG) neurons. When tested on a variety of recombinantly expressed voltage-gated ion channels, Mu8.1 displayed the highest potency against the R-type (Cav2.3) calcium channel. Ca2+ signals from Mu8.1-sensitive DRG neurons were also inhibited by SNX-482, a known spider peptide modulator of Cav2.3 and voltage-gated K+ (Kv4) channels. Our findings highlight the potential of Mu8.1 as a molecular tool to identify and study neuronal subclasses expressing Cav2.3. Importantly, this multidisciplinary study showcases the potential of uncovering novel structures and bioactivities within the largely unexplored group of macro-conotoxins.

Exploring and targeting potential druggable antimicrobial resistance targets ArgS, SecY, and MurA in Staphylococcus sciuri with TCM inhibitors through a subtractive genomics strategy

Staphylococcus sciuri (also currently Mammaliicoccus sciuri) are anaerobic facultative and non-motile bacteria that cause significant human pathogenesis such as endocarditis, wound infections, peritonitis, UTI, and septic shock. Methicillin-resistant S. sciuri (MRSS) strains also infects animals that include healthy broilers, cattle, dogs, and pigs. The emergence of MRSS strains thereby poses a serious health threat and thrives the scientific community towards novel treatment options. Herein, we investigated the druggable genome of S. sciuri by employing subtractive genomics that resulted in seven genes/proteins where only three of them were predicted as final targets. Further mining the literature showed that the ArgS (WP_058610923), SecY (WP_058611897), and MurA (WP_058612677) are involved in the multi-drug resistance phenomenon. After constructing and verifying the 3D protein homology models, a screening process was carried out using a library of Traditional Chinese Medicine compounds (consisting of 36,043 compounds). The molecular docking and simulation studies revealed the physicochemical stability parameters of the docked TCM inhibitors in the druggable cavities of each protein target by identifying their druggability potential and maximum hydrogen bonding interactions. The simulated receptor-ligand complexes showed the conformational changes and stability index of the secondary structure elements. The root mean square deviation (RMSD) graph showed fluctuations due to structural changes in the helix-coil-helix and beta-turn-beta changes at specific points where the pattern of the RMSD and root mean square fluctuation (RMSF) (< 1.0 Å) support any major domain shifts within the structural framework of the protein-ligand complex and placement of ligand was well complemented within the binding site. The β-factor values demonstrated instability at few points while the radius of gyration for structural compactness as a time function for the 100-ns simulation of protein-ligand complexes showed favorable average values and denoted the stability of all complexes. It is assumed that such findings might facilitate researchers to robustly discover and develop effective therapeutics against S. sciuri alongside other enteric infections.

Combined genomic and transcriptomic analysis reveals the contribution of tandem duplication genes to low-temperature adaptation in perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is an agronomically important cool-season grass species that is widely used as forage for ruminant animal production and cultivated in temperate regions for the establishment of lawns. However, the underlying genetic mechanism of the response of L. perenne to low temperature is still unclear. In the present study, we performed a comprehensive study and identified 3,770 tandem duplication genes (TDGs) in L. perenne, and evolutionary analysis revealed that L. perenne might have undergone a duplication event approximately 7.69 Mya. GO and KEGG pathway functional analyses revealed that these TDGs were mainly enriched in photosynthesis, hormone-mediated signaling pathways and responses to various stresses, suggesting that TDGs contribute to the environmental adaptability of L. perenne. In addition, the expression profile analysis revealed that the expression levels of TDGs were highly conserved and significantly lower than those of all genes in different tissues, while the frequency of differentially expressed genes (DEGs) from TDGs was much higher than that of DEGs from all genes in response to low-temperature stress. Finally, in-depth analysis of the important and expanded gene family indicated that the members of the ELIP subfamily could rapidly respond to low temperature and persistently maintain higher expression levels during all low temperature stress time points, suggesting that ELIPs most likely mediate low temperature responses and help to facilitate adaptation to low temperature in L. perenne. Our results provide evidence for the genetic underpinning of low-temperature adaptation and valuable resources for practical application and genetic improvement for stress resistance in L. perenne.

Comprehensive identification of glutathione peroxidase (GPX) gene family in response to abiotic stress in pepper (Capsicum annuum L.)

Plant glutathione peroxidase (GPX) plays an important role in the maintenance of cell homeostasis and in the antioxidant response in plants. In this study, the peroxidase (GPX) gene family was identified in the whole genome of pepper using bioinformatic method. As a result, a total of 5 CaGPX genes were identified, which were unevenly distributed on 3 of the 12 chromosomes of pepper genome. Based on phylogenetic analysis, 89 GPX genes in 17 species from lower plants to higher plants can be divided into 4 groups GroupⅠ, Group Ⅱ, Group Ⅲ, Group Ⅳ). The MEME Suite analysis of GPX proteins shows that all these proteins contain four highly conserved motifs, as well as other conserved sequences and amino acid residues. Gene structure analysis revealed the conservative exon-intron organization pattern of these genes. In the promoter region of CaGPX genes, many cis elements of plant hormone and abiotic stress response were identified in each of CaGPX proteins. In addition, expression patterns of CaGPX genes in different tissues, developmental stages and responses to abiotic stress were also performed. The results of qRT-PCR showed that the transcripts of CaGPX genes varied greatly under abiotic stress at different time points. There results suggest that the GPX gene family of pepper may play a role in plant development andstress response. In conclusion, our research provides new insights into the evolution of pepper GPX gene family, and understanding for functional of these genes in response to abiotic stresses.

Hsf transcription factor gene family in peanut (Arachis hypogaea L.): genome-wide characterization and expression analysis under drought and salt stresses

Heat shock transcription factors (Hsfs) play important roles in plant developmental regulations and various stress responses. In present study, 46 Hsf genes in peanut (AhHsf) were identified and analyzed. The 46 AhHsf genes were classed into three groups (A, B, and C) and 14 subgroups (A1-A9, B1-B4, and C1) together with their Arabidopsis homologs according to phylogenetic analyses, and 46 AhHsf genes unequally located on 17 chromosomes. Gene structure and protein motif analysis revealed that members from the same subgroup possessed similar exon/intron and motif organization, further supporting the results of phylogenetic analyses. Gene duplication events were found in peanut Hsf gene family via syntenic analysis, which were important in Hsf gene family expansion in peanut. The expression of AhHsf genes were detected in different tissues using published data, implying that AhHsf genes may differ in function. In addition, several AhHsf genes (AhHsf5, AhHsf11, AhHsf20, AhHsf24, AhHsf30, AhHsf35) were induced by drought and salt stresses. Furthermore, the stress-induced member AhHsf20 was found to be located in nucleus. Notably, overexpression of AhHsf20 was able to enhance salt tolerance. These results from this study may provide valuable information for further functional analysis of peanut Hsf genes.

There are three major Neisseria gonorrhoeae β-lactamase plasmid variants which are associated with specific lineages and carry distinct TEM alleles

Neisseria gonorrhoeae is a significant threat to global health with an estimated incidence of over 80 million cases each year and high levels of antimicrobial resistance. The gonococcal β-lactamase plasmid, pbla, carries the TEM β-lactamase, which requires only one or two amino acid changes to become an extended-spectrum β-lactamase (ESBL); this would render last resort treatments for gonorrhoea ineffective. Although pbla is not mobile, it can be transferred by the conjugative plasmid, pConj, found in N. gonorrhoeae. Seven variants of pbla have been described previously, but little is known about their frequency or distribution in the gonococcal population. We characterised sequences of pbla variants and devised a typing scheme, Ng_pblaST that allows their identification from whole genome short-read sequences. We implemented Ng_pblaST to assess the distribution of pbla variants in 15 532 gonococcal isolates. This demonstrated that only three pbla variants commonly circulate in gonococci, which together account for >99 % of sequences. The pbla variants carry different TEM alleles and are prevalent in distinct gonococcal lineages. Analysis of 2758 pbla-containing isolates revealed the co-occurrence of pbla with certain pConj types, indicating co-operativity between pbla and pConj variants in the spread of plasmid-mediated AMR in N. gonorrhoeae. Understanding the variation and distribution of pbla is essential for monitoring and predicting the spread of plasmid-mediated β-lactam resistance in N. gonorrhoeae.

Genome-wide identification and analysis of the WRKY gene family and low-temperature stress response in Prunus sibirica

BACKGROUND

WRKY transcription factors are a prominent gene family in plants, playing a crucial role in various biological processes including development, metabolism, defense, differentiation, and stress response. Although the WRKY gene family has been extensively studied and analysed in numerous plant species, research on Prunus sibirica's WRKY genes (PsWRKY) remains lacking.

RESULTS

This study analysed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, and Gene ontology (GO) annotation of PsWRKY gene family members using bioinformatics methods based on the whole-genome data of P. sibirica. In total, 55 WRKYs were identified in P. sibirica and were heterogeneously distributed on eight chromosomes. Based on the phylogenetic analysis, these WRKYs were classified into three major groups: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III. Members of different subfamilies have different cis-acting elements, conserved motifs, and intron-exon structures, indicating functional heterogeneity of the WRKY family. Prediction of subcellular localisation indicated that PsWRKYs were mainly located in the nucleus. Twenty pairs of duplicated genes were identified, and segmental duplication events may play an important role in PsWRKY gene family expansion. Analysis of the Ka/Ks ratio showed that the PsWRKY family's homologous genes were primarily purified by selection. Additionally, GO annotation analysis showed that the WRKY gene family was mainly involved in responses to stimuli, immune system processes, and reproductive processes. Furthermore, quantitative real-time PCR (qRT-PCR) analysis showed that 23 PsWRKYs were highly expressed in one or more tissues (pistils and roots) and PsWRKYs showed specific expression patterns under different low-temperature stress conditions.

CONCLUSIONS

Our results provide a scientific basis for the further exploration and functional validation of WRKYs in P. sibirica.

Isolation of three MiDi19-4 genes from mango, the ectopic expression of which confers early flowering and enhances stress tolerance in transgenic Arabidopsis

MAIN CONCLUSION

Three Di19-4 genes were identified in mango. Overexpression of MiDi19-4B in A. thaliana promoted earlier flowering and enhanced drought, salt, and ABA resistance. Drought-induced protein 19 (Di19) is a drought-induced protein that is mainly involved in multiple stress responses. Here, three Di19-4 genes (MiDi19-4A/B/C) in mango (Mangifera indica L.) were identified, and the coding sequences (CDS) had lengths of 684, 666, and 672 bp and encoded proteins with 228, 222, and 224 amino acids, respectively. The promoters of the MiDi19-4 genes contained phytohormone-, light-, and abiotic stress-responsive elements. The MiDi19-4 genes were expressed in every tissue and highly expressed in leaves. Moreover, MiDi19-4 genes were highly correlated with the vegetative growth period and induced by polyethylene glycol (PEG) or salt stress. MiDi19-4B displayed the highest expression during the vegetative growth period and then showed decreased expression, and MiDi19-4B was highly expressed at both the late stage of the vegetative growth period and the initial stage of the flowering induction period. The 35S::GFP-MiDi19-4B fusion protein was located in the cell nucleus. The transgenic plants ectopically expressing MiDi19-4B exhibited earlier flowering and increased expression patterns of FRUITFULL (AtFUL), APETALA1 (AtAP1), and FLOWERING LOCUS T (AtFT). The drought and salt tolerance of MiDi19-4B transgenic plants was significantly increased, and these plants showed decreased sensitivity to abscisic acid (ABA) and considerably increased expression levels of drought- and salt-related genes and ABA signalling pathway genes. Additionally, bimolecular fluorescence complementation (BiFC) experiments revealed that the MiDi19-4B protein interacted with CAULIFLOWER (MiCAL1), MiCAL2, MiAP1-1, and MiAP1-2. Taken together, these results highlighted the important regulatory roles of MiDi19-4B in tolerance to multiple abiotic stresses and in flowering.

Chimeric vaccine design against the epidemic Langya Henipavirus using immunoinformatics and validation via immune simulation approaches

In July 2022, a new virus called Langya virus (LayV) was discovered in China in patients who had a fever. This virus is a type of Henipavirus (HNV) and is considered a potential threat as it could spread from animals to humans. It causes respiratory disease with symptoms including fever, coughing, and fatigue and is closely linked to two other henipaviruses that are known to infect humans, namely Hendra and Nipah viruses. These viruses may cause fatal respiratory illnesses. Investigators believe that the LayV is spread by shrews, and may have infected humans directly or via an intermediary species. Thus, the use of vaccines or immunizations against LayV is an alternate strategy for disease prevention. In this study, we employed various immunoinformatics methods to predict B cell, HTL and T cell epitopes from the LayV proteome in order to find the most promising candidate for a LayV vaccine. The most potent epitopes that are immunogenic and non-allergenic were joined with each other through suitable linkers. Human β-defensin 2 was employed as an adjuvant to increase the immunogenicity of the vaccine construct. The final sequence of a multi-epitope vaccine construct was modelled for docking with TLRs. Concisely, our results suggest that the docked complexes of vaccine-TLRs seemed to be stable. Additionally, in silico cloning was done using E. coli as the host in order to validate the expression of our designed vaccine construct. The GC content of 54.39% and CAI value of 0.94 revealed that the vaccine component expresses efficiently in the host. This study presents the novel vaccine construct for LayV which will be essential for further experimental validations to confirm the immunogenicity and safety of the proposed vaccine structure, and eventually to treat HNV-related diseases.

1H, 15N, 13C backbone and sidechain resonance assignments and secondary structure of mouse NOTCH1 EGF27

NOTCH1 is a transmembrane receptor in metazoans that is linked to a variety of disorders. The receptor contains an extracellular domain (ECD) with 36 tandem epidermal growth factor-like (EGF) repeats. The ECD is responsible for intercellular signaling via protein-ligand interactions with neighboring cells. Each EGF repeat consists of approximately 40 amino acids and 3 conserved disulfide bonds. The Abruptex region (EGF24-29) is critical for NOTCH1 signaling and is known for its missense mutations. Certain EGF repeats are modified with the addition of O-linked glycans and many have calcium binding sites, which give each EGF repeat a unique function. It has been shown that the loss of the O-fucose site of EGF27 alters NOTCH1 activity. To investigate the role of glycosylation in the NOTCH1 signaling pathway, nuclear magnetic resonance spectroscopy has been employed to study the structures of EGF27 and its glycoforms. Here, we report the backbone and sidechain 1H, 15N, and 13C-resonance assignments of the unmodified EGF27 protein and the predicted secondary structure derived from the assigned chemical shifts.

The effects of free Cys residues on the structure, activity, and tetrameric stability of mammalian uricase

Mammalian uricases contain four conserved cysteine (Cys) residues, but little is known about their structures and functions. In this study, we first confirmed that all four Cys residues are free and not involved in disulfide bond formation, using canine uricase as a model protein. Cys residues had a greater effect on stability than on activity based on single Cys-to-Ser (serine) substitutions. Circular dichroism (CD) and homology modeling indicated that C188S reduces β-sheet contents and inter- and intra-subunit hydrophobic interaction, potentially impairing the core tetrameric β-barrel structure of the tunneling-fold protein, and ultimately decreased the tetrameric stability. Additionally, the inactivation of C188S during the stability tests may be a complex process involving depolymerization followed by irregular aggregation. Double mutations or thiol blockage of Cys188 and Cys195 significantly disrupted the formation and stability of tetrameric uricase, which may be mediated by the free thiols in Cys residues. The present results demonstrated that the free Cys residues are essential for tetrameric formation and stability in mammalian uricase. This implies that free cysteine residues, although not involved in disulfide bonding, may play important structural roles in certain proteins, underscoring the significance of the hydrophobic characteristics of the free thiols in Cys residues. KEY POINTS: • Four Cys residues are not involved in disulfide bonding in mammalian uricase. • The hydrophobicity of free thiols is critical for tetrameric stability in uricase. • Free Cys residues can serve structural roles without involving in disulfide bonds.

New robust subtilisins from halotolerant and halophilic Bacillaceae

The aim of the present study was the characterisation of three true subtilisins and one phylogenetically intermediate subtilisin from halotolerant and halophilic microorganisms. Considering the currently growing enzyme market for efficient and novel biocatalysts, data mining is a promising source for novel, as yet uncharacterised enzymes, especially from halophilic or halotolerant Bacillaceae, which offer great potential to meet industrial needs. Both halophilic bacteria Pontibacillus marinus DSM 16465^T and Alkalibacillus haloalkaliphilus DSM 5271^T and both halotolerant bacteria Metabacillus indicus DSM 16189 and Litchfieldia alkalitelluris DSM 16976^T served as a source for the four new subtilisins SPPM, SPAH, SPMI and SPLA. The protease genes were cloned and expressed in Bacillus subtilis DB104. Purification to apparent homogeneity was achieved by ethanol precipitation, desalting and ion-exchange chromatography. Enzyme activity could be observed between pH 5.0-12.0 with an optimum for SPPM, SPMI and SPLA around pH 9.0 and for SPAH at pH 10.0. The optimal temperature for SPMI and SPLA was 70 °C and for SPPM and SPAH 55 °C and 50 °C, respectively. All proteases showed high stability towards 5% (w/v) SDS and were active even at NaCl concentrations of 5 M. The four proteases demonstrate potential for future biotechnological applications. KEY POINTS: • Halophilic and halotolerant Bacillaceae are a valuable source of new subtilisins. • Four new subtilisins were biochemically characterised in detail. • The four proteases show potential for future biotechnological applications.

Genome-wide identification of RUB activating enzyme and conjugating enzyme gene families and their expression analysis under abiotic stresses in Capsicum annuum

NEDD8/RUB, as a ubiquitin-like protein, participates in the post-translational modification of protein and requires unique E1, E2, and E3 enzymes to bind to its substrate. The RUB E1 activating enzyme and E2 conjugating enzyme play a significant role in the neddylation. However, it is unknown whether RUB E1 and E2 exist in pepper and what its function is. In this study, a total of three putative RUB E1 and five RUB E2 genes have been identified in the pepper genome. Subsequently, their physical and chemical properties, gene structure, conserved domains and motifs, phylogenetic relationship, and cis-acting elements were analyzed. The structure and conserved domain of RUB E1 and E2 are similar to that of Arabidopsis and tomato. The RUB E1 and E2 genes were randomly distributed on seven chromosomes, and there were two pairs of collinearity between pepper and Arabidopsis and eight pairs of collinearity between pepper and tomato. Phylogenetic analysis reveals that RUB E1 and E2 genes of pepper have a closer relationship with that of tomato, potato, and Nicotiana attenuate. The cis-elements of RUB E1 and E2 genes contained hormone response and stress response. RUB E1 and E2 genes were expressed in at least one tissue and CaRCE1.3 and CaRCE2.1 were exclusively expressed in flowers and anthers. Moreover, the expression of RUB E1 genes (CaECR1, CaAXR1.1, and CaAXR1.2) and RUB E2 genes (CaRCE1.1, CaRCE1.2, and CaRCE2.1) was increased to varying degrees under low-temperature, drought, salt, ABA, and IAA treatments, while CaRCE1.3 and CaRCE2.2 were down-regulated under low-temperature treatment. In addition, these genes were hardly expressed under MeJA treatment. In summary, this study provides a theoretical foundation to explore the role of RUB E1 and E2 in the response of plants to stress.

Screening and identification of cyprinid herpesvirus 2 (CyHV-2) ORF55-interacting proteins by phage display

BACKGROUND

Cyprinid herpesvirus 2 (CyHV-2) is a pathogenic fish virus belonging to family Alloherpesviridae. The CyHV-2 gene encoding thymidine kinase (TK) is an important virulence-associated factor. Therefore, we aimed to investigate the biological function of open reading frame 55 (ORF55) in viral replication.

METHODS

Purified CyHV-2 ORF55 protein was obtained by prokaryotic expression, and the interacting peptide was screened out using phage display. Host interacting proteins were then predicted and validated.

RESULTS

ORF55 was efficiently expressed in the prokaryotic expression system. Protein and peptide interaction prediction and dot-blot overlay assay confirmed that peptides identified by phage display could interact with the ORF55 protein. Comparing the peptides to the National Center for Biotechnology Information database revealed four potential interacting proteins. Reverse transcription quantitative PCR results demonstrated high expression of an actin-binding Rho-activating protein in the latter stages of virus-infected cells, and molecular docking, cell transfection and coimmunoprecipitation experiments confirmed that it interacted with the ORF55 protein.

CONCLUSION

During viral infection, the ORF55 protein exerts its biological function through interactions with host proteins. The specific mechanisms remain to be further explored.

Functional characterization of the GhNRT2.1e gene reveals its significant role in improving nitrogen use efficiency in Gossypium hirsutum

Background

Nitrate is the primary type of nitrogen available to plants, which is absorbed and transported by nitrate transporter 2 (NRT2) at low nitrate conditions.

Methods

Genome-wide identification of NRT2 genes in G. hirsutum was performed. Gene expression patterns were revealed using RNA-seq and qRT-PCR. Gene functions were characterized using overexpression in A. thaliana and silencing in G. hirsutum. Protein interactions were verified by yeast two-hybrid and luciferase complementation imaging (LCI) assays.

Results

We identified 14, 14, seven, and seven NRT2 proteins in G. hirsutum, G. barbadense, G. raimondii, and G. arboreum. Most NRT2 proteins were predicted in the plasma membrane. The NRT2 genes were classified into four distinct groups through evolutionary relationships, with members of the same group similar in conserved motifs and gene structure. The promoter regions of NRT2 genes included many elements related to growth regulation, phytohormones, and abiotic stresses. Tissue expression pattern results revealed that most GhNRT2 genes were specifically expressed in roots. Under low nitrate conditions, GhNRT2 genes exhibited different expression levels, with GhNRT2.1e being the most up-regulated. Arabidopsis plants overexpressing GhNRT2.1e exhibited increased biomass, nitrogen and nitrate accumulation, nitrogen uptake and utilization efficiency, nitrogen-metabolizing enzyme activity, and amino acid content under low nitrate conditions. In addition, GhNRT2.1e-silenced plants exhibited suppressed nitrate uptake and accumulation, hampered plant growth, affected nitrogen metabolism processes, and reduced tolerance to low nitrate. The results showed that GhNRT2.1e could promote nitrate uptake and transport under low nitrate conditions, thus effectively increasing nitrogen use efficiency (NUE). We found that GhNRT2.1e interacts with GhNAR2.1 by yeast two-hybrid and LCI assays.

Discussion

Our research lays the foundation to increase NUE and cultivate new cotton varieties with efficient nitrogen use.

Development and Clinical Application of a Multilocus Sequence Typing Scheme for Bacteroides fragilis Based on Whole-Genome Sequencing Data

Bacteroides fragilis is among the most abundant and pathogenic bacterial species in the gut microbiota and is associated with diarrheal disease in children, inflammatory bowel disease, and the development of colorectal cancer. It is increasingly resistant to potent antimicrobial agents such as carbapenems and metronidazole, making it among the most resistant anaerobic bacteria. These factors combined call for increased monitoring of B. fragilis and its population structure on a worldwide scale. Here, we present a possible solution through the development of a multilocus sequence typing scheme (MLST). The scheme is based on seven core gene fragments: groL (hsp60), rpoB, recA, dnaJ, rprX, prfA, and fusA. These gene fragments possess high discriminatory power while retaining concordance with whole core genome-based phylogenetic analysis. The scheme proved capable of differentiating B. fragilis isolates at the strain level. It offers a standardized method for molecular typing and can be applied to isolates from various sampling backgrounds, such as patient isolates, environmental samples, and strains obtained from food and animal sources. In total, 567 B. fragilis genomes were sequence typed and their isolate data collected. The MLST scheme clearly divided the B. fragilis population into two divisions based on the presence of the cfiA and cepA resistance genes. However, no other specific subpopulations within the analyzed genomes were found to be associated with any specific diseases or geographical location. With this MLST scheme, we hope to provide a powerful tool for the study and monitoring of B. fragilis on an international scale. IMPORTANCE Here, we present the first MLST scheme for Bacteroides fragilis, one of the most abundant pathogenic bacteria in the human gut microbiota. The scheme enables standard classification and monitoring of B. fragilis on a worldwide scale and groups the majority of current isolate data in one place. A more unified approach to the collection and analysis of B. fragilis data could provide crucial insights into how the pathogen operates and develops as a species. Close monitoring of B. fragilis is especially relevant, as it is increasingly resistant to potent antimicrobial agents and engages in horizontal gene transfer with other bacteria. Hopefully, this approach will guide new discoveries into how B. fragilis evolves and interacts with its human host. Additionally, the scheme could potentially be applied to other species of the genus Bacteroides.

In silico designing of a novel polyvalent multi-subunit peptide vaccine leveraging cross-immunity against human visceral and cutaneous leishmaniasis: an immunoinformatics-based approach

CONTEXT

Leishmaniasis is a group of vector-borne infectious diseases caused by over 20 pathogenic Leishmania species that are endemic in many tropical and subtropical countries. The emergence of drug-resistant strains, the adverse side effects of anti-Leishmania drugs, and the absence of a preventative vaccination strategy threaten the sensitive population. Recently, many groups of researchers have exploited the field of reverse vaccinology to develop vaccines, focusing chiefly on inducing immunity against either visceral or cutaneous leishmaniasis.

METHODS

This present work involves retrieving twelve experimentally validated leishmanial antigenic protein sequences from the UniProt database, followed by their antigenicity profiling employing ANTIGENpro and Vaxijen 2.0 servers. MHC-binding epitopes for the same were predicted using both NetCTL 1.2 and SYFPEITHI servers, while epitopes for B cell were computed using ABCpred and BepiPred 2.0 servers. The screened epitopes with significantly higher scores were utilized for designing the vaccine construct with appropriate linkers and natural adjuvant. The secondary and tertiary structures of the synthetic peptide were determined by conditional random fields, shallow neural networks, and profile-profile threading alignment with iterative structure assembly simulations, respectively. The 3-D vaccine model was validated through CASP10-tested refinement and the MolProbity web server. Molecular docking and multi-scale normal mode analysis simulation were performed to analyze the best vaccine-TLR complex. Finally, computational immune simulation findings revealed promising cellular and humoral immune responses, suggesting that the engineered chimeric peptide is a potential broad-spectrum vaccine against visceral and cutaneous leishmaniasis.

Genome-wide characterization and identification of Trihelix transcription factors and expression profiling in response to abiotic stresses in Chinese Willow (Salix matsudana Koidz)

Trihelix transcription factors (TTF) are a class of light-responsive proteins with a typical triple-helix structure (helix-loop-helix-loop-helix). Members of this gene family play an important role in plant growth and development, especially in various abiotic stress responses. Salix matsudana Koidz is an allotetraploid ornamental forest tree that is widely planted for its excellent resistance to stress, but no studies on its Trihelix gene family have been reported. In this study, the Trihelix gene family was analyzed at the genome-wide level in S. matsudana. A total of 78 S. matsudana Trihelix transcription factors (SmTTFs) were identified, distributed on 29 chromosomes, and classified into four subfamilies (GT-1, GT-2, SH4, SIP1) based on their structural features. The gene structures and conserved functional domains of these Trihelix genes are similar in the same subfamily and differ between subfamilies. The presence of multiple stress-responsive cis-elements on the promoter of the S. matsudana Trihelix gene suggests that the S. matsudana Trihelix gene may respond to abiotic stresses. Expression pattern analysis revealed that Trihelix genes have different functions during flooding stress, salt stress, drought stress and low temperature stress in S. matsudana. Given that SmTTF30, as a differentially expressed gene, has a faster response to flooding stress, we selected SmTTF30 for functional studies. Overexpression of SmTTF30 in Arabidopsis thaliana (Arabidopsis) enhances its tolerance to flooding stress. Under flooding stress, the leaf cell activity and peroxidase activity (POD) of the overexpression strain were significantly higher than the leaf cell activity and POD of the wild type, and the malondialdehyde (MDA) content was significantly lower than the MDA content of the wild type. Thus, these results suggest that SmTTF30 enhances plant flooding tolerance and plays a positive regulatory role in plant flooding tolerance.

Limosilactobacillus fermentum 3872 That Produces Class III Bacteriocin Forms Co-Aggregates with the Antibiotic-Resistant Staphylococcus aureus Strains and Induces Their Lethal Damage

LF3872 was isolated from the milk of a healthy lactating and breastfeeding woman. Earlier, the genome of LF3872 was sequenced, and a gene encoding unique bacteriocin was discovered. We have shown here that the LF3872 strain produces a novel thermolabile class III bacteriolysin (BLF3872), exhibiting antimicrobial activity against antibiotic-resistant Staphylococcus aureus strains. Sequence analysis revealed the two-domain structural (lysozyme-like domain and peptidase M23 domain) organization of BLF3872. At least 25% residues of this protein are expected to be intrinsically disordered. Furthermore, BLF3872 is predicted to have a very high liquid-liquid phase separation. According to the electron microscopy data, the bacterial cells of LF3872 strain form co-aggregates with the S. aureus 8325-4 bacterial cells. LF3872 produced bacteriolysin BLF3872 that lyses the cells of the S. aureus 8325-4 mastitis-inducing strain. The sensitivity of the antibiotic-resistant S. aureus collection strains and freshly isolated antibiotic-resistant strains was tested using samples from women with lactation mastitis; the human nasopharynx and oral cavity; the oropharynx of pigs; and the cows with a diagnosis of clinical mastitis sensitive to the lytic action of the LF3872 strain producing BLF3872. The co-cultivation of LF3872 strain with various antibiotic-resistant S. aureus strains for 24 h reduced the level of living cells of these pathogens by six log. The LF3872 strain was found to be able to co-aggregate with all studied S. aureus strains. The cell-free culture supernatant of LF3872 (CSLF3872) induced S. aureus cell damage and ATP leakage. The effectiveness of the bacteriolytic action of LF3872 strain did not depend on the origin of the S. aureus strains. The results reported here are important for the creation of new effective drugs against antibiotic-resistant strains of S. aureus circulating in humans and animals.

The design and functional characterization of a novel hybrid antimicrobial peptide from Esculentin-1a and melittin

Antimicrobial agents are one of the most widely used drugs in medicine. In the last fifty years, the misuse of these agents caused the emergence of resistant strains of bacteria that led to an increase in life-threatening infections. The need to develop new agents has become a priority, and antimicrobial peptides attained high consideration. The antimicrobial activities of a novel In-house designed hybrid cationic peptide (BKR1) were studied against different strains of Gram-negative bacteria. This was done using the broth dilution method as outlined by the Clinical and Laboratory Institute (CLSI). Checkerboard assy was employed to investigate the synergistic activity of BKR1 peptide with four antibiotics (Levofloxacin, chloramphenicol, rifampicin, and ampicillin). Finally, the cytotoxicity of BKR1 was evaluated against human blood cells and mammalian kidney cells (Vero cells). BKR1 displayed bactericidal activity against tested strains of Gram-negative bacteria, with zero hemolytic effects. It also acts as a strong adjuvant with levofloxacin, chloramphenicol, and rifampicin against resistant strains of P. aeruginosa and E. coli. This study represents the design and elucidation of the antimicrobial activities of a novel hybrid antimicrobial peptide named (BKR1). Our results indicate thar BKR1 is a promising candidate to treat resistant infectious diseases individually or as an adjuvant with conventional antibiotics.

CYP6DW3 Metabolizes Imidacloprid to Imidacloprid-urea in Whitefly (Bemisia tabaci)

Bemisia tabaci has developed high resistance to many insecticides and causes substantial agricultural and economic losses annually. The insecticide resistance of whitefly has been widely reported in previous studies; however, the underlying mechanism remains little known. In this study, we cloned two P450 genes: CYP6DW3 and CYP6DW5v1; these genes were markedly overexpressed in imidacloprid-resistant whitefly populations compared with susceptible populations, and knockdown of these genes decreased the imidacloprid resistance of whitefly. Moreover, heterologous expression of whitefly P450 genes in SF9 cells and metabolic studies showed that the CYP6DW3 protein could metabolize 14.11% imidacloprid and produced imidacloprid-urea in vitro. Collectively, the expression levels of CYP6DW3 and CYP6DW5v1 are positively correlated with imidacloprid resistance in B. tabaci. Our study further reveals that cytochrome P450 enzymes affect the physiological activities related to resistance in insects, which helps scholars more deeply understand the resistance mechanism, and contributes to the development of integrated pest management framework.

Genome-Wide Re-Identification and Analysis of CrRLK1Ls in Tomato

The Catharanthus roseus receptor-like kinase 1-like (CrRLK1L), which is a vital member of the plant receptor-like kinase family, plays versatile roles in plant growth, development, and stress response. Although the primary screening of tomato CrRLK1Ls has been reported previously, our knowledge of these proteins is still scarce. Using the latest genomic data annotations, a genome-wide re-identification and analysis of the CrRLK1Ls in tomatoes were conducted. In this study, 24 CrRLK1L members were identified in tomatoes and researched further. Subsequent gene structures, protein domains, Western blot analyses, and subcellular localization analyses all confirmed the accuracy of the newly identified SlCrRLK1L members. Phylogenetic analyses showed that the identified SlCrRLK1L proteins had homologs in Arabidopsis. Evolutionary analysis indicated that two pairs of the SlCrRLK1L genes had predicted segmental duplication events. Expression profiling analyses demonstrated that the SlCrRLK1L genes were expressed in various tissues, and most of them were up- or down-regulated by bacteria and PAMP treatments. Together, these results will lay the foundation for elaborating the biological roles of SlCrRLK1Ls in tomato growth, development, and stress response.

ADH Gene Cloning and Identification of Flooding-Responsive Genes in Taxodium distichum (L.) Rich

As a flooding-tolerant tree species, Taxodium distichum has been utilized in afforestation projects and proven to have important value in flooding areas. Alcohol dehydrogenase (ADH), which participates in ethanol fermentation, is essential for tolerance to the anaerobic conditions caused by flooding. In a comprehensive analysis of the ADH gene family in T. distichum, TdADHs were cloned on the basis of whole-genome sequencing, and then bioinformatic analysis, subcellular localization, and gene expression level analysis under flooding were conducted. The results show that the putative protein sequences of 15 cloned genes contained seven TdADHs and eight TdADH-like genes (one Class III ADH included) that were divided into five clades. All the sequences had an ADH_N domain, and except for TdADH-likeE2, all the other genes had an ADH_zinc_N domain. Moreover, the TdADHs in clades A, B, C, and D had a similar motif composition. Additionally, the number of TdADH amino acids ranged from 277 to 403, with an average of 370.13. Subcellular localization showed that, except for TdADH-likeD3, which was not expressed in the nucleus, the other genes were predominantly expressed in both the nucleus and cytosol. TdADH-likeC2 was significantly upregulated in all three organs (roots, stems, and leaves), and TdADHA3 was also highly upregulated under 24 h flooding treatment; the two genes might play key roles in ethanol fermentation and flooding tolerance. These findings offer a comprehensive understanding of TdADHs and could provide a foundation for the molecular breeding of T. distichum and current research on the molecular mechanisms driving flooding tolerance.

Design of a multi-epitope vaccine against six Nocardia species based on reverse vaccinology combined with immunoinformatics

Background

Nocardia genus, a complex group of species classified to be aerobic actinomycete, can lead to severe concurrent infection as well as disseminated infection, typically in immunocompromised patients. With the expansion of the susceptible population, the incidence of Nocardia has been gradually growing, accompanied by increased resistance of the pathogen to existing therapeutics. However, there is no effective vaccine against this pathogen yet. In this study, a multi-epitope vaccine was designed against the Nocardia infection using reverse vaccinology combined with immunoinformatics approaches.

Methods

First, the proteomes of 6 Nocardia subspecies Nocardia subspecies (Nocardia farcinica, Nocardia cyriacigeorgica, Nocardia abscessus, Nocardia otitidiscaviarum, Nocardia brasiliensis and Nocardia nova) were download NCBI (National Center for Biotechnology Information) database on May 1st, 2022 for the target proteins selection. The essential, virulent-associated or resistant-associated, surface-exposed, antigenic, non-toxic, and non-homologous with the human proteome proteins were selected for epitope identification. The shortlisted T-cell and B-cell epitopes were fused with appropriate adjuvants and linkers to construct vaccines. The physicochemical properties of the designed vaccine were predicted using multiple online servers. The Molecular docking and molecular dynamics (MD) simulation were performed to understand the binding pattern and binding stability between the vaccine candidate and Toll-like receptors (TLRs). The immunogenicity of the designed vaccines was evaluated via immune simulation.

Results

3 proteins that are essential, virulent-associated or resistant-associated, surface-exposed, antigenic, non-toxic, and non-homologous with the human proteome were selected from 218 complete proteome sequences of the 6 Nocardia subspecies epitope identification. After screening, only 4 cytotoxic T lymphocyte (CTL) epitopes, 6 helper T lymphocyte (HTL) epitopes, and 8 B cell epitopes that were antigenic, non-allergenic, and non-toxic were included in the final vaccine construct. The results of molecular docking and MD simulation showed that the vaccine candidate has a strong affinity for TLR2 and TLR4 of the host and the vaccine-TLR complexes were dynamically stable in the natural environment. The results of the immune simulation indicated that the designed vaccine had the potential to induce strong protective immune responses in the host. The codon optimization and cloned analysis showed that the vaccine was available for mass production.

Conclusion

The designed vaccine has the potential to stimulate long-lasting immunity in the host, but further studies are required to validate its safety and efficacy.

Loosenin-Like Proteins from Phanerochaete carnosa Impact Both Cellulose and Chitin Fiber Networks

Microbial expansin-related proteins are ubiquitous across bacterial and fungal organisms and reportedly play a role in the modification and deconstruction of cell wall polysaccharides, including lignocellulose. So far, very few microbial expansin-related proteins, including loosenins and loosenin-like (LOOL) proteins, have been functionally characterized. Herein, four LOOLs encoded by Phanerochaete carnosa and belonging to different subfamilies (i.e., PcaLOOL7 and PcaLOOL9 from subfamily A and PcaLOOL2 and PcaLOOL12 from subfamily B) were recombinantly produced and the purified proteins were characterized using diverse cellulose and chitin substrates. The purified PcaLOOLs weakened cellulose filter paper and cellulose nanofibril networks (CNF); however, none significantly boosted cellulase activity on the selected cellulose substrates (Avicel and Whatman paper). Although fusing the family 63 carbohydrate-binding module (CBM63) of BsEXLX1 encoded by Bacillus subtilis to PcaLOOLs increased their binding to cellulose, the CBM63 fusion appeared to reduce the cellulose filter paper weakening observed using wild-type proteins. Binding of PcaLOOLs to alpha-chitin was considerably higher than that to cellulose (Avicel) and was pH dependent, with the highest binding at pH 5.0. Amendment of certain PcaLOOLs in fungal liquid cultivations also impacted the density of the cultivated mycelia. The present study reveals the potential of fungal expansin-related proteins to impact both cellulose and chitin networks and points to a possible biological role in fungal cell wall processing. IMPORTANCE The present study deepens investigations of microbial expansin-related proteins and their applied significance by (i) reporting a detailed comparison of diverse loosenins encoded by the same organism, (ii) considering both cellulosic and chitin-containing materials as targeted substrates, and (iii) investigating the impact of the C-terminal carbohydrate binding module (CBM) present in other expansin-related proteins on loosenin function. By revealing the potential of fungal loosenins to impact both cellulose and chitin-containing networks, our study reveals a possible biological and applied role of loosenins in fungal cell wall processing.

Genome-wide identification of the class III peroxidase gene family of sugarcane and its expression profiles under stresses

INTRODUCTION

Plant-specific Class III peroxidases (PRXs) play a crucial role in lignification, cell elongation, seed germination, and biotic and abiotic stresses.

METHODS

The class III peroxidase gene family in sugarcane were identified by bioinformatics methods and realtime fluorescence quantitative PCR.

RESULTS

Eighty-two PRX proteins were characterized with a conserved PRX domain as members of the class III PRX gene family in R570 STP. The ShPRX family genes were divided into six groups by the phylogenetic analysis of sugarcane, Saccharum spontaneum, sorghum, rice, and Arabidopsis thaliana. The analysis of promoter cis-acting elements revealed that most ShPRX family genes contained cis-acting regulatory elements involved in ABA, MeJA, light responsiveness, anaerobic induction, and drought inducibility. An evolutionary analysis indicated that ShPRXs was formed after Poaceae and Bromeliaceae diverged, and tandem duplication events played a critical role in the expansion of ShPRX genes of sugarcane. Purifying selection maintained the function of ShPRX proteins. SsPRX genes were differentially expressed in stems and leaves at different growth stages in S. spontaneum. However, ShPRX genes were differentially expressed in the SCMV-inoculated sugarcane plants. A qRT-PCR analysis showed that SCMV, Cd, and salt could specifically induce the expression of PRX genes of sugarcane.

DISCUSSION

These results help elucidate the structure, evolution, and functions of the class III PRX gene family in sugarcane and provide ideas for the phytoremediation of Cd-contaminated soil and breeding new sugarcane varieties resistant to sugarcane mosaic disease, salt, and Cd stresses.

Genome-wide characterization of the PP2C gene family in peanut (Arachis hypogaea L.) and the identification of candidate genes involved in salinity-stress response

Plant protein phosphatase 2C (PP2C) play important roles in response to salt stress by influencing metabolic processes, hormone levels, growth factors, etc. Members of the PP2C family have been identified in many plant species. However, they are rarely reported in peanut. In this study, 178 PP2C genes were identified in peanut, which were unevenly distributed across the 20 chromosomes, with segmental duplication in 78 gene pairs. AhPP2Cs could be divided into 10 clades (A-J) by phylogenetic analysis. AhPP2Cs had experienced segmental duplications and strong purifying selection pressure. 22 miRNAs from 14 different families were identified, targeting 57 AhPP2C genes. Gene structures and motifs analysis exhibited PP2Cs in subclades AI and AII had high structural and functional similarities. Phosphorylation sites of AhPP2C45/59/134/150/35/121 were predicted in motifs 2 and 4, which located within the catalytic site at the C-terminus. We discovered multiple MYB binding factors and ABA response elements in the promoter regions of the six genes (AhPP2C45/59/134/150/35/121) by cis-elements analysis. GO and KEGG enrichment analysis confirmed AhPP2C-A genes in protein binding, signal transduction, protein modification process response to abiotic stimulus through environmental information processing. Based on RNA-Seq data of 22 peanut tissues, clade A AhPP2Cs showed a varying degree of tissue specificity, of which, AhPP2C35 and AhPP2C121 specifically expressed in seeds, while AhPP2C45/59/134/150 expressed in leaves and roots. qRT-PCR indicated that AhPP2C45 and AhPP2C134 displayed significantly up-regulated expression in response to salt stress. These results indicated that AhPP2C45 and AhPP2C134 could be candidate PP2Cs conferring salt tolerance. These results provide further insights into the peanut PP2C gene family and indicate PP2Cs potentially involved in the response to salt stress, which can now be further investigated in peanut breeding efforts to obtain cultivars with improved salt tolerance.

Genome-wide investigation of SnRK2 gene family in two jute species: Corchorus olitorius and Corchorus capsularis

BACKGROUND

Sucrose non-fermenting-1 (SNF1)-related protein kinase 2 (SnRK2), a plant-specific serine/threonine kinase family, is associated with metabolic responses, including abscisic acid signaling under biotic and abiotic stresses. So far, no information on a genome-wide investigation and stress-mediated expression profiling of jute SnRK2 is available. Recent whole-genome sequencing of two Corchorus species prompted to identify and characterize this SnRK2 gene family.

RESULT

We identified seven SnRK2 genes of each of Corchorus olitorius (Co) and C. capsularis (Cc) genomes, with similar physico-molecular properties and sub-group patterns of other models and related crops. In both species, the SnRK2 gene family showed an evolutionarily distinct trend. Highly variable C-terminal and conserved N-terminal regions were observed. Co- and CcSnRK2.3, Co- and CcSnRk2.5, Co- and CcSnRk2.7, and Co- and CcSnRK2.8 were upregulated in response to drought and salinity stresses. In waterlogging conditions, Co- and CcSnRk2.6 and Co- and CcSnRK2.8 showed higher activity when exposed to hypoxic conditions. Expression analysis in different plant parts showed that SnRK2.5 in both Corchorus species is highly expressed in fiber cells providing evidence of the role of fiber formation.

CONCLUSION

This is the first comprehensive study of SnRK2 genes in both Corchorus species. All seven genes identified in this study showed an almost similar pattern of gene structures and molecular properties. Gene expression patterns of these genes varied depending on the plant parts and in response to abiotic stresses.

Genome-wide identification and characterization of LcCCR13 reveals its potential role in lignin biosynthesis in Liriodendron chinense

Introduction

Wood formation is closely related to lignin biosynthesis. Cinnamoyl-CoA reductase (CCR) catalyzes the conversion of cinnamoyl-CoA to cinnamaldehydes, which is the initiation of the lignin biosynthesis pathway and a crucial point in the manipulation of associated traits. Liriodendron chinense is an economically significant timber tree. Nevertheless, the underlying mechanism of wood formation in it remains unknown; even the number of LcCCR family members in this species is unclear.

Materials and Results

This study aimed to perform a genome-wide identification of genes(s) involved in lignin biosynthesis in L. chinense via RT-qPCR assays and functional verification. Altogether, 13 LcCCR genes were identified that were divided into four major groups based on structural and phylogenetic features. The gene structures and motif compositions were strongly conserved between members of the same groups. Subsequently, the expression patterns analysis based on RNA-seq data indicated that LcCCR5/7/10/12/13 had high expression in the developing xylem at the stem (DXS). Furthermore, the RT-qPCR assays showed that LcCCR13 had the highest expression in the stem as compared to other tissues. Moreover, the overexpression of the LcCCR13 in transgenic tobacco plants caused an improvement in the CCR activity and lignin content, indicating that it plays a key role in lignin biosynthesis in the stems.

Discussion

Our research lays a foundation for deeper investigation of the lignin synthesis and uncovers the genetic basis of wood formation in L. chinense.

Genome-wide identification of the B3 gene family in soybean and the response to melatonin under cold stress

INTRODUCTION

Melatonin is a multipotent molecule that exists widely in animals and plants and plays an active regulatory role in abiotic stresses. The B3 superfamily is a ubiquitous transcription factor with a B3 functional domain in plants, which can respond temporally to abiotic stresses by activating defense compounds and plant hormones. Despite the fact that the B3 genes have been studied in a variety of plants, their role in soybean is still unknown.

METHODS

The regulation of melatonin on cold resistance of soybean and the response of B3 genes to cold stress were investigated by measuring biochemical indexes of soybean. Meanwhile, the genome-wide identification of B3 gene family was conducted in soybean, and B3 genes were analyzed based on phylogeny, motifs, gene structure, collinearity, and cis-regulatory elements analysis.

RESULTS

We found that cold stress-induced oxidative stress in soybean by producing excessive reactive oxygen species. However, exogenous melatonin treatment could increase the content of endogenous melatonin and other hormones, including IAA and ABA, and enhance the antioxidative system, such as POD activity, CAT activity, and GSH/GSSG, to scavenge ROS. Furthermore, the present study first revealed that melatonin could alleviate the response of soybean to cold stress by inducing the expression of B3 genes. In addition, we first identified 145 B3 genes in soybean that were unevenly distributed on 20 chromosomes. The B3 gene family was divided into 4 subgroups based on the phylogeny tree constructed with protein sequence and a variety of plant hormones and stress response cis-elements were discovered in the promoter region of the B3 genes, indicating that the B3 genes were involved in several aspects of the soybean stress response. Transcriptome analysis and results of qRT-PCR revealed that most GmB3 genes could be induced by cold, the expression of which was also regulated by melatonin. We also found that B3 genes responded to cold stress in plants by interacting with other transcription factors.

DISCUSSION

We found that melatonin regulates the response of soybean to cold stress by regulating the expression of the transcription factor B3 gene, and we identified 145 B3 genes in soybean. These findings further elucidate the potential role of the B3 gene family in soybean to resist low-temperature stress and provide valuable information for soybean functional genomics study.

A multiepitope vaccine candidate against infectious bursal disease virus using immunoinformatics-based reverse vaccinology approach

Infectious bursal disease virus is the causative agent of infectious bursal disease (Gumboro disease), a highly contagious immunosuppressive disease of chicken with a substantial economic impact on small- and large-scale poultry industries worldwide. Currently, live attenuated vaccines are widely used to control the disease in chickens despite their issues with safety (immunosuppression and bursal atrophy) and efficiency (breaking through the maternally-derived antibody titer). To overcome the drawbacks, the current study has, for the first time, attempted to construct a computational model of a multiepitope based vaccine candidate against infectious bursal disease virus, which has the potential to overcome the safety and protection issues found in the existing live-attenuated vaccines. The current study used a reverse vaccinology based immunoinformatics approach to construct the vaccine candidate using major and minor capsid proteins of the virus, VP2 and VP3, respectively. The vaccine construct was composed of four CD8⁺ epitopes, seven CD4⁺ T-cell epitopes, 11 B-cell epitopes and a Cholera Toxin B adjuvant, connected using appropriate flexible peptide linkers. The vaccine construct was evaluated as antigenic with VaxiJen Score of 0.6781, immunogenic with IEDB score of 2.89887 and non-allergenic. The 55.64 kDa construct was further evaluated for its physicochemical characteristics, which revealed that it was stable with an instability index of 16.24, basic with theoretical pI of 9.24, thermostable with aliphatic index of 86.72 and hydrophilic with GRAVY score of -0.256. The docking and molecular dynamics simulation studies of the vaccine construct with Toll-like receptor-3 revealed fair structural interaction (binding affinity of -295.94 kcal/mol) and complex stability. Further, the predicted induction of antibodies and cytokines by the vaccine construct indicated the possible elicitation of the host's immune response against the virus. The work is a significant attempt to develop next-generation vaccines against the infectious bursal disease virus though further experimental studies are required to assess the efficacy and protectivity of the proposed vaccine candidate in vivo.

Genome-wide identification and expression analysis of two-component system genes in sweet potato (Ipomoea batatas L.)

Two-component system (TCS), which comprises histidine kinases (HKs), histidine phosphotransfer proteins (HPs), and response regulators (RRs), plays essential roles in regulating plant growth, development, and response to various environmental stimuli. TCS genes have been comprehensively identified in various plants, while studies on the genome-wide identification and analysis of TCS in sweet potato were still not reported. Therefore, in this study, a total of 90 TCS members consisting of 20 HK(L)s, 11 HPs, and 59 RRs were identified in the genome of Ipomoea batatas. Furthermore, their gene structures, conserved domains, and phylogenetic relationships were analyzed in detail. Additionally, the gene expression profiles in various organs were analyzed, and response patterns to adverse environmental stresses were investigated. The results showed that these 90 TCS genes were mapped on 15 chromosomes with a notably uneven distribution, and the expansion of TCS genes in sweet potato was attributed to both segmental and tandem duplications. The majority of the TCS genes showed distinct organ-specific expression profiles, especially in three types of roots (stem roots, fibrous roots, tuberous roots). Moreover, most of the TCS genes were either induced or suppressed upon treatment with abiotic stresses (drought, salinity, cold, heat) and exogenous phytohormone abscisic acid (ABA). In addition, the yeast-two hybrid system was used to reveal the HK-HP-RR protein-protein interactions. IbHP1, IbHP2, IbHP4, and IbHP5 could interact with three HKs (IbHK1a, IbHK1b, and IbHK5), and also interact with majority of the type-B RRs (IbRR20-IbRR28), while no interaction affinity was detected for IbHP3. Our systematic analyses could provide insights into the characterization of the TCS genes, and further the development of functional studies in sweet potato.

In silico peptide-based therapeutics against human colorectal cancer by the activation of TLR5 signaling pathways

OBJECTIVE

Colorectal cancer (CRC) is the third leading cause of cancer-related deaths in both men and women. Toll-like receptor 5 (TLR5), an autoimmune signaling receptor that plays a role in cancer, can be exploited for the suppression of human colon cancer. Salmonella flagellin protein, a novel agonist of TLR5 activating downstream signaling, could be a basis for designing anticancer peptides.

METHODS

The three-dimensional crystal structure of TLR5 (PDB ID: 3J0A, Resolution = 26.0 Å) was optimized using the AMBER force field in the YASARA suit. In silico enzymatic digestion tool, PeptideCutter, was used to identify peptides from Salmonella flagellin, an agonist against human TLR5. The 3D structure of the peptides was generated using PEP-FOLD3. These peptides were screened against human TLR5 using shape complementarity principles based on the binding affinity and interactions with the active residue of TLR5 monomer, and the selected peptides were further validated by molecular dynamic (MD) simulation.

RESULTS

In this study, we generated 42 peptides from Salmonella flagellin protein by in silico protein digestion. Then, based on a new hidden Markov model sub-optimal conformation sampling approach as well as the size of the fragments, we select 38 effective peptides from these 42 cleavages. These peptides were screened against the monomeric Xray structure of human TLR5 using shape complementarity principles. Based on the binding affinity and interactions with the active residue of TLR5 monomer (residues 294 and 366 of TLR5), nine top-scored peptides were selected for the initial molecular dynamic (MD) simulation. Among these peptides, Clv10, Clv17, and Clv28 showed high stability and less flexibility during MD simulation. A 1 μs MD simulation was performed on TLR5-Clv10, TLR-Clv17, and TLR5-Clv28 complexes to further analyze the stability, conformational changes, and binding mode (Clv10, Clv17, and Clv28). During this MD study, the peptides showed high salt bridges and ionic interactions with residue ASP294 and residue ASP366 throughout the simulation and remained in the concave of the human TLR5 monomer. The RMSD and Rg values showed that the peptide-protein complexes become stable after 200 ns of contraction and extraction.

CONCLUSION

These findings can facilitate the rational design of selected peptides as an agonist of TLR5, which have antitumor activity, suppress colorectal cancer tumors, and can be used as promising candidates and novel agonists of TLR5.

Current Update on Rotavirus in-Silico Multiepitope Vaccine Design

Rotavirus gastroenteritis is one of the leading causes of pediatric morbidity and mortality worldwide in infants and under-five populations. The World Health Organization (WHO) recommended global incorporation of the rotavirus vaccine in national immunization programs to alleviate the burden of the disease. Implementation of the rotavirus vaccination in certain regions of the world brought about a significant and consistent reduction of rotavirus-associated hospitalizations. However, the efficacy of licensed vaccines remains suboptimal in low-income countries where the incidences of rotavirus gastroenteritis continue to happen unabated. The problem of low efficacy of currently licensed oral rotavirus vaccines in low-income countries necessitates continuous exploration, design, and development of new rotavirus vaccines. Traditional vaccine development is a complex, expensive, labor-intensive, and time-consuming process. Reverse vaccinology essentially utilizes the genome and proteome information on pathogens and has opened new avenues for in-silico multiepitope vaccine design for a plethora of pathogens, promising time reduction in the complete vaccine development pipeline by complementing the traditional vaccinology approach. A substantial number of reviews on licensed rotavirus vaccines and those under evaluation are already available in the literature. However, a collective account of rotavirus in-silico vaccines is lacking in the literature, and such an account may further fuel the interest of researchers to use reverse vaccinology to expedite the vaccine development process. Therefore, the main focus of this review is to summarize the research endeavors undertaken for the design and development of rotavirus vaccines by the reverse vaccinology approach utilizing the tools of immunoinformatics.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Siraitia siamensis

WRKY transcription factors, as the largest gene family in higher plants, play an important role in various biological processes including growth and development, regulation of secondary metabolites, and stress response. In this study, we performed genome-wide identification and analysis of WRKY transcription factors in S. siamensis. A total of 59 SsWRKY genes were identified that were distributed on all 14 chromosomes, and these were classified into three major groups based on phylogenetic relationships. Each of these groups had similar conserved motifs and gene structures. We compared all the S. siamensis SsWRKY genes with WRKY genes identified from three diverse plant species, and the results implied that segmental duplication and tandem duplication play an important roles in the evolution processes of the WRKY gene family. Promoter region analysis revealed that SsWRKY genes included many cis-acting elements related to plant growth and development, phytohormone response, and both abiotic and biotic stress. Expression profiles originating from the transcriptome database showed expression patterns of these SsWRKY genes in four different tissues and revealed that most genes are expressed in plant roots. Fifteen SsWRKY genes with low-temperature response motifs were surveyed for their gene expression under cold stress, showing that most genes displayed continuous up-regulation during cold treatment. Our study provides a foundation for further study on the function and regulatory mechanism of the SsWRKY gene family.

Fasciola gigantica vaccine construct: an in silico approach towards identification and design of a multi-epitope subunit vaccine using calcium binding EF-hand proteins

Continuous attempts have been made to pinpoint candidate vaccine molecules and evaluate their effectiveness in order to commercialise such vaccines for the treatment of tropical fascioliasis in livestock. The pathophysiology of fascioliasis can be related to liver damage brought on by immature flukes that migrate and feed, as well as immunological reactions to chemicals produced by the parasites and alarm signals brought on by tissue damage. Future research should, in our opinion, concentrate on the biology of invasive parasites and the resulting immune responses, particularly in the early stages of infection. The goal of the current study was to use the calcium-binding proteins from F. gigantica to create a multi-epitope subunit vaccine. The adjuvant, B-cell epitopes, CTL epitopes, and HTL epitopes that make up the vaccine construct are all connected by certain linkers. The antigenicity, allergenicity, and physiochemical properties of the vaccine construct were examined. The vaccine construct was docked with toll-like receptor 2, and simulations of the molecular dynamics of the complex's stability, interaction, and dynamics were run. After performing in silico cloning and immunosimulation, it was discovered that the construct was suitable for further investigation. New vaccination technologies and adjuvant development are advancing our food safety procedures since vaccines are seen as safe and are accepted by the user community. This research is also applicable to the F. hepatica system.

The functional analysis of sugar transporter proteins in sugar accumulation and pollen tube growth in pummelo (Citrus grandis)

Sugar transporter proteins (STPs) play vital roles in sugar transport and allocation of carbon sources in plants. However, the evolutionary dynamics of this important gene family and their functions are still largely unknown in citrus, which is the largest fruit crop in the world. In this study, fourteen non-redundant CgSTP family members were identified in pummelo (Citrus grandis). A comprehensive analysis based on the biochemical characteristics, the chromosomal location, the exon-intron structures and the evolutionary relationships demonstrated the conservation and the divergence of CgSTPs. Moreover, CgSTP4, 11, 13, 14 were proofed to be localized in plasma membrane and have glucose transport activity in yeast. The hexose content were significantly increased with the transient overexpression of CgSTP11 and CgSTP14. In addition, antisense repression of CgSTP4 induced the shorter pollen tube length in vitro, implying the potential role of CgSTP4 in pummelo pollen tube growth. Taken together, this work explored a framework for understanding the physiological role of CgSTPs and laid a foundation for future functional studies of these members in citrus species.

Genome-wide identification and analyses of cotton high-affinity nitrate transporter 2 family genes and their responses to stress

Nitrate transporters (NRTs) are crucial for the uptake, use, and storage of nitrogen by plants. In this study, 42 members of the GhNRT2 (Nitrate Transporter 2 family) were found in the four different cotton species. The conserved domains, phylogenetic relationships, physicochemical properties, subcellular localization, conserved motifs, gene structure, cis-acting elements, and promoter region expression patterns of these 42 members were analyzed. The findings confirmed that members of the NRT2 family behaved typically, and subcellular localization tests confirmed that they were hydrophobic proteins that were mostly located on the cytoplasmic membrane. The NRT2 family of genes with A.thaliana and rice underwent phylogenetic analysis, and the results revealed that GhNRT2 could be divided into three groups. The same taxa also shared similar gene structure and motif distribution. The composition of cis-acting elements suggests that most of the expression of GhNRT2 may be related to plant hormones, abiotic stress, and photoreactions. The GhNRT2 gene was highly expressed, mainly in roots. Drought, salt, and extreme temperature stress showed that GhNRT2 gene expression was significantly up-regulated or down-regulated, indicating that it may be involved in the stress response of cotton. In general, the genes of the NRT2 family of cotton were comprehensively analyzed, and their potential nitrogen uptake and utilization functions in cotton were preliminarily predicted. Additionally, we provide an experimental basis for the adverse stress conditions in which they may function.

Whole Exome Sequencing Identified Two Single Nucleotide Polymorphisms of Human Leukocyte Antigen-DRB5 in Familial Sarcoidosis in China

BACKGROUND

Sarcoidosis is a multisystem granulomatous disorder whose etiology is related to genetic and immunological factors. Familial aggregation and ethnic prevalence suggest a genetic predisposition and inherited susceptibility to sarcoidosis.

OBJECTIVE

This study aimed to identify suspected risk loci for familial sarcoidosis patients.

METHODS

We conducted whole exome sequencing on two sarcoidosis patients and five healthy family members in a Chinese family for a case-control study. The two sarcoidosis patients were siblings who showed chronic disease.

RESULTS

The Gene Ontology results showed single nucleotide polymorphisms in three genes, including human leukocyte antigen (HLA)-DRB1, HLA-DRB5, and KIR2DL4, associated with both 'antigen processing and presentation' and 'regulation of immune response.' Sanger sequencing verified two nonsynonymous mutations in HLA-DRB5 (rs696318 and rs115817940) located on 6p21.3 in the major histocompatibility complex (MHC) class II beta 1 region. The structural model simulated on Prot- Param protein analysis by the Expert Protein Analysis System predicted that the hydropathy index changed at two mutation sites (rs696318: p.F96L, -1.844 to -1.656 and rs115817940: p.T106N, -0.322 to -0.633), which indicated the probability of changes in peptide-binding selectivity.

CONCLUSION

Our results indicated that two nonsynonymous mutations of HLA-DRB5 have been identified in two sarcoidosis siblings, while their healthy family members do not have the mutations. The two HLA-DRB5 alleles may influence genetic susceptibility and chronic disease progression through peptide mutations on the MHC class II molecule among the two affected family members.

Genome-wide characterization of major latex protein gene family in peanut and expression analyses under drought and waterlogging stress

Peanut is an important oilseed crop around the world which provides vegetable oil, protein and vitamins for humans. Major latex-like proteins (MLPs) play important roles in plant growth and development, as well as responses to biotic and abiotic stresses. However, their biological function in peanut is still unclear. In this study, a genome-wide identification of MLP genes in cultivated peanut and two diploid ancestor species was analyzed to determine their molecular evolutionary characteristics and the expression profile under drought and waterlogging stress conditions. Firstly, a total of 135 MLP genes were identified from the genome of tetraploid peanut (Arachis hypogaea) and two diploid species Arachis. duranensis and Arachis. ipaensis. Then, phylogenetic analysis revealed that MLP proteins were divided into five different evolutionary groups. These genes were distributed unevenly at the ends of chromosomes 3, 5, 7, 8, 9 and 10 in three Arachis species. The evolution of MLP gene family in peanut was conserved and led by tandem and segmental duplication. The prediction analysis of cis-acting elements showed that the promoter region of peanut MLP genes contained different proportions of transcription factors, plant hormones-responsive elements and so on. The expression pattern analysis showed that they were differentially expressed under waterlogging and drought stress. These results of this study provide a foundation for further research on the function of the important MLP genes in peanut.

Genome-wide analysis of R2R3-MYB transcription factors in Boehmeria nivea (L.) gaudich revealed potential cadmium tolerance and anthocyanin biosynthesis genes

Gene family, especially MYB as one of the largest transcription factor family in plants, the study of its subfunctional characteristics is a key step in the study of plant gene function. The sequencing of ramie genome provides a good opportunity to study the organization and evolutionary characters of the ramie MYB gene at the whole genome level. In this study, a total of 105 BnGR2R3-MYB genes were identified from ramie genome and subsequently grouped into 35 subfamilies according to phylogeny divergence and sequences similarity. Chromosomal localization, gene structure, synteny analysis, gene duplication, promoter analysis, molecular characteristics and subcellular localization were accomplished using several bioinformatics tools. Collinearity analysis showed that the segmental and tandem duplication events is the dominant form of the gene family expansion, and duplications prominent in distal telomeric regions. Highest syntenic relationship was obtained between BnGR2R3-MYB genes and that of Apocynum venetum (88). Furthermore, transcriptomic data and phylogenetic analysis revealed that BnGMYB60, BnGMYB79/80 and BnGMYB70 might inhibit the biosynthesis of anthocyanins, and UPLC-QTOF-MS data further supported the results. qPCR and phylogenetic analysis revealed that the six genes (BnGMYB9, BnGMYB10, BnGMYB12, BnGMYB28, BnGMYB41, and BnGMYB78) were cadmium stress responsive genes. Especially, the expression of BnGMYB10/12/41 in roots, stems and leaves all increased more than 10-fold after cadmium stress, and in addition they may interact with key genes regulating flavonoid biosynthesis. Thus, a potential link between cadmium stress response and flavonoid synthesis was identified through protein interaction network analysis. The study thus provided significant information into MYB regulatory genes in ramie and may serve as a foundation for genetic enhancement and increased productivity.

Genome-wide identification and expression analysis of DREB genes in alfalfa (Medicago sativa) in response to cold stress

Dehydration-responsive element-binding proteins (DREBs) belong to members of the AP2/ERF transcription factor superfamily, which has been reported to involve various abiotic-stress responses and tolerance in plants. However, research on the DREB-family is still limited in alfalfa (Medicago sativa L.), a forage legume cultivated worldwide. The recent genome-sequence release of the alfalfa cultivar "XinJiangDaYe" allowed us to identify 172 DREBs by a multi-step homolog search. The phylogenetic analysis indicated that such MsDREBs could be classified into 5 groups, namely A-1 (56 members), A-2 (39), A-3 (3), A-4 (61) and 13 (A-5 (13), thus adding substantial new members to the DREB-family in alfalfa. Furthermore, a comprehensive survey in silico of conserved motif, gene structure, molecular weight, and isoelectric point (pI) as well as gene expression was conducted. The resulting data showed that, for cold-stress response, 33 differentially expressed MsDREBs were identified with a threshold of Log2-fold > 1, and most of which were transcriptionally upregulated within 48 h during a cold treatment(s). Moreover, the expression profiling of MsDREBs from two ecotypes of alfalfa subspecies i.e. M. sativa ssp. falcata (F56, from a colder region of Central Asia) and M. sativa ssp. sativa (B47, from Near East) revealed that most of the cold-stress responsive MsDREBs exhibited a significantly lower expression in F56, leading to a proposal of the existence of a distinct mechanism(s) for cold tolerance regulated by DREB-related action, which would have been evolved in alfalfa with a genotypic specificity. Additionally, by examining the transcriptome of a freezing-tolerance species (M. sativa cv. Zhaodong), eight DREBs were found to be implicated in a long-term freezing-stress adaptation with a great potential. Taken together, the current genome-wide identification in alfalfa points to the importance of some MsDREBs in the cold-stress response, providing some promising molecular targets to be functionally characterized for the improvement of cold tolerance in crops including alfalfa.