Multiple Sequence Alignment Analysis Research Articles

Emergence of allostery through reorganization of protein residue network architecture.

Despite more than a century of study, consensus on the molecular basis of allostery remains elusive. A comparison of allosteric and non-allosteric members of a protein family can shed light on this important regulatory mechanism, and the bacterial biotin protein ligases, which catalyze post-translational biotin addition, provide an ideal system for such comparison. While the Class I bacterial ligases only function as enzymes, the bifunctional Class II ligases use the same structural architecture for an additional transcription repression function. This additional function depends on allosterically activated homodimerization followed by DNA binding. In this work, we used experimental, computational network, and bioinformatics analyses to uncover distinguishing features that enable allostery in the Class II biotin protein ligases. Experimental studies of the Class II Escherichia coli protein indicate that catalytic site residues are critical for both catalysis and allostery. However, allostery also depends on amino acids that are more broadly distributed throughout the protein structure. Energy-based community network analysis of representative Class I and Class II proteins reveals distinct residue community architectures, interactions among the communities, and responses of the network to allosteric effector binding. Bioinformatics mutual information analyses of multiple sequence alignments indicate distinct networks of coevolving residues in the two protein families. The results support the role of divergent local residue community network structures both inside and outside of the conserved enzyme active site combined with distinct inter-community interactions as keys to the emergence of allostery in the Class II biotin protein ligases.

Identification and functional characterization of the SUMO system in sweet potato under salt and drought stress

Sumoylation is a crucial post-translation modification (PTM) that is the covalent attachment of SUMO molecules to the substrate catalyzed by enzyme cascade. Sumoylation is essential in almost every physiological process of plants, particularly in response to abiotic stress. However, little is known about sumoylation in sweet potato (Ipomoea batatas), the world’s seventh most important food crop. In this study, 17 sweet potato SUMO system genes have been cloned and functionally characterized. Multiple sequence alignment and phylogenetic analysis showed sweet potato SUMO system proteins had conserved domains and activity sites. IbSUMOs, IbSAE1, and IbSCE1 were localized in the cytoplasm and nucleus. E3 SUMO ligases showed nuclear or punctate localization. In vitro sumoylation assay confirmed the catalytic activity of sweet potato SUMO system components. Heterologous expression of IbSIZ1 genes in Arabidopsis atsiz1 mutant rescued the defective germination and growth phenotype. IbSCE1a/b and IbSIZ1a/b/c were salt and drought responsive genes. Heterologous expression of IbSCE1a/b/c improved the drought tolerance of Arabidopsis thaliana, while IbSIZ1a/b/c significantly enhanced the salt and drought tolerance. Our findings define that the SUMO system in sweet potato shared with conserved function but also possessed specific characterization. The resources presented here would facilitate uncovering the significance of sumoylation in sweet potato.

Transcriptome analysis of Populus × canadensis ‘Zhongliao1’ in response to low temperature stress

BackgroundLow temperatures are known to limit the growth and geographical distribution of poplars. Although some transcriptomic studies have been conducted to explore the response of poplar leaves to cold stress, only a few have comprehensively analyzed the effects of low temperature on the transcriptome of poplars and identified genes related to cold stress response and repair of freeze–thaw injury.ResultsWe exposed the Euramerican poplar Zhongliao1 to low temperatures; after stems were exposed to − 40℃, 4℃, and 20℃, the mixture of phloem and cambium was collected for transcriptome sequencing and bioinformatics analysis. A total of 29,060 genes were detected, including 28,739 known genes and 321 novel genes. Several differentially expressed genes (n = 36) were found to be involved in the Ca2+ signaling pathway, starch–sucrose metabolism pathway, abscisic acid signaling pathway, and DNA repair. They were functionally annotated; glucan endo-1,3-beta-glucosidase and UDP-glucuronosyltransferase genes, for instance, showed a close relationship with cold resistance. The expression of 11 differentially expressed genes was verified by qRT-PCR; RNA-Seq and qRT-PCR data were found to be consistent, which validated the robustness of our RNA-Seq findings. Finally, multiple sequence alignment and evolutionary analysis were performed, the results of which suggested a close association between several novel genes and cold resistance in Zhongliao1.ConclusionWe believe that the cold resistance and freeze–thaw injury repair genes identified in this study are of great significance for cold tolerance breeding.

Improvement of the stability and catalytic efficiency of heparan sulfate N-sulfotransferase for preparing N-sulfated heparosan.

The chemo-enzymatic and enzymatic synthesis of heparan sulfate and heparin are considered as an attractive alternative to the extraction of heparin from animal tissues. Sulfation of the hydroxyl group at position 2 of the deacetylated glucosamine is a prerequisite for subsequent enzymatic modifications. In this study, multiple strategies, including truncation mutagenesis based on B-factor values, site-directed mutagenesis guided by multiple sequence alignment, and structural analysis were performed to improve the stability and activity of human N-sulfotransferase. Eventually, a combined variant Mut02 (MBP-hNST-NΔ599-602/S637P/S741P/E839P/L842P/K779N/R782V) was successfully constructed, whose half-life at 37°C and catalytic activity were increased by 105-fold and 1.35-fold, respectively. After efficient overexpression using the Escherichia coli expression system, the variant Mut02 was applied to N-sulfation of the chemically deacetylated heparosan. The N-sulfation content reached around 82.87% which was nearly 1.88-fold higher than that of the wild-type. The variant Mut02 with high stability and catalytic efficiency has great potential for heparin biomanufacturing.

The burdock database: a multi-omic database for Arctium lappa, a food and medicinal plant

BackgroundBurdock is a biennial herb of Asteraceae found in Northern Europe, Eurasia, Siberia, and China. Its mature dry fruits, called Niu Bang Zi, are recorded in various traditional Chinese medicine books. With the development of sequencing technology, the mitochondrial, chloroplast, and nuclear genomes, transcriptome, and sequence-related amplified polymorphism (SRAP) fingerprints of burdock have all been reported. To make better use of this data for further research and analysis, a burdock database was constructed.ResultsThis burdock multi-omics database contains two burdock genome datasets, two transcriptome datasets, eight burdock chloroplast genomes, one burdock mitochondrial genome, one A. tomentosum chloroplast genome, one A. tomentosum mitochondrial genome, 26 phenotypes of burdock varieties, burdock rhizosphere-associated microorganisms, and chemical constituents of burdock fruit, pericarp, and kernel at different growth stages (using UPLC-Q-TOF–MS). The wild and cultivation distribution of burdock in China was summarized, and the main active components and pharmacological effects of burdock currently reported were concluded. The database contains ten central functional modules: Home, Genome, Transcriptome, Jbrowse, Search, Tools, SRAP fingerprints, Associated microorganisms, Chemical, and Publications. Among these, the “Tools” module can be used to perform sequence homology alignment (Blast), multiple sequence alignment analysis (Muscle), homologous protein prediction (Genewise), primer design (Primer), large-scale genome analysis (Lastz), and GO and KEGG enrichment analyses (GO Enrichment and KEGG Enrichment).ConclusionsThe database URL is http://210.22.121.250:41352/. This burdock database integrates molecular and chemical data to provide a comprehensive information and analysis platform for interested researchers and can be of immense help to the cultivation, breeding, and molecular pharmacognosy research of burdock.

A long-term retrospective analysis of the haemorrhagic fever with renal syndrome epidemic from 2005 to 2021 in Jiangxi Province, China

Jiangxi is one of the provinces in China most seriously affected by the haemorrhagic fever with renal syndrome (HFRS) epidemic. The aim of this paper was to systematically explore the HFRS epidemic in Jiangxi from the perspective of Hantavirus (HV) prevalence in rodents and humans and virus molecular characteristics. Individual information on all HFRS cases in Jiangxi from 2005 to 2021 was extracted from the China Information System for Disease Control and Prevention. All S and M fragment sequences of the Seoul virus and Hantan virus strains uploaded by Jiangxi and its neighbouring provinces and some representative sequences from provinces in China or some countries of Southeast Asia with the highest HV prevalence were retrieved and downloaded from NCBI GenBank. Periodogram and spatial autocorrelation were adopted for temporal periodicity and spatial clustering analysis of the HFRS epidemic. Joinpoint regression was utilized to explore the changing morbidity trend patterns of HFRS. Multiple sequence alignment and amino acid variation analysis were used to explore the homology and variation of strain prevalence in Jiangxi. Based on monthly morbidity time series, the periodogram analysis showed that the prevalence of HFRS had periodicities of 6 months and 12 months. Spatial autocorrelation analysis showed that HFRS distributed in Jiangxi was not random, with a “High-High” clustering area around Gaoan County. HFRS morbidity among the 0 ~ 15-year-old and ~ 61-year-old or older populations in Jiangxi increased significantly during the period of 2008–2015. Generally, HFRS morbidity was significantly positively correlated with the index of rat with virus (IRV) (r = 0.742) in the counties surrounding Gaoan from 2005 to 2019. HTNV strains in Jiangxi were in one independent branch, while the SEOV strains in Jiangxi were relatively more diverse. Both the YW89-15 and GAW30/2021 strains shared approximately 85% nucleotide homology and approximately 97% amino acid homology with their corresponding standard strains and vaccine strains. GAW30/2021 and YW89-15 had some amino acid site variations in nucleoprotein, glycoprotein precursor and RNA-dependent polymerase with their corresponding vaccine strains Z10 (HTNV) and Z37 (SEOV). The HFRS epidemic in Jiangxi has obvious temporal periodicity and spatial clustering, and the significant increase in the non-Immunization Expanded Program (EPI) targeted population (children and elderly) suggests that HFRS vaccination in this population needs to be considered. Although applying the EPI played a certain role in curbing the incidence of HFRS in Jiangxi from the perspective of ecological epidemiology, HTNV and SEOV strains prevalent in Jiangxi have some amino acid site variations compared to their corresponding vaccine strains, suggesting that HV variation needs to be continuously monitored in the future to observe vaccine protective efficiency.

An insight into the active sites of the catalytic basic protein subunit PB1 of the RNA polymerase of human influenza viruses

RNA polymerase from human influenza viruses is a heterotrimeric enzyme and performs the crucial functions of both transcription and replication for multiplication of the viruses in human cells. The heterotrimeric enzyme is made up of two basic protein subunits (PB1 and PB2) and an acidic protein subunit (PA). All the three subunits perform well-defined function(s) in the transcription and replication processes in the human cells. The basic protein subunit PB1 is shown to possess the polymerization activity, whereas the PB2 and the PA subunits are found to be involved in the cap-snatching and proofreading (PR) activities, respectively. The polymerase activity in the catalytic subunit, PB1, is found to be an RNA-dependent RNA polymerase (RdRp). Multiple sequence alignment (MSA) analysis of the PB1 subunits from all the three human influenza viruses, A, B and C shows large number of highly conserved peptides, amino acid motifs and invariant amino acids. Site-directed mutagenesis (SDM) analysis and X-ray crystallographic data have shown that two completely conserved motifs, viz. –GDN- and –SDD-, are involved in binding to the catalytic Mg2+ ion. These data are in close agreement with the MSA analysis data of the polymerases from all the three human influenza viruses. Furthermore, two highly conserved polymerase catalytic regions are identified in the PB1 subunits by sequence similarity to other DNA/RNA polymerases and hence, are proposed to function in the nucleotidyl transfer activities. Presence of the two catalytic regions suggest that the polymerase may function in a dual mode, i.e., in phase I, in association with the cap-snatching subunit PB2, it could be involved in the synthesis of mRNAs (transcription mode) and once enough proteins are made from the mRNAs, in the second phase, in association with PR exonuclease subunit PA, it could switch to the replication mode to synthesize error-free, exact copies of the viral genome. For both the activities, it could use the same invariant catalytic Mg2+-binding –GDN- and –SDD- motifs.

In Silico Phylogeny, Sequence and Structure Analyses of Fungal Thermoacidophilic GH11 Xylanases

Thermoacidophilic xylanase enzymes are mostly preferred for use as animal feed additives. In this study, we performed in silico phylogeny, sequence, structure, and enzyme-docked complex analyses of six thermoacidophilic GH11 xylanases belonging to various fungal species (Gymnopus androsaceus xylanase = GaXyl, Penicilliopsis zonata xylanase = PzXyl, Aspergillus neoniger xylanase = AnXyl, Calocera viscosa xylanase = CvXyl, Acidomyces richmondensis xylanase = ArXyl, Oidiodendron maius xylanase = OmXyl). To do this, amino acid sequences of six fungal thermoacidophilic GH11 xylanases, belonging to unreviewed protein entries in the UniProt/TrEMBL database, were investigated at molecular phylogeny and amino acid sequence levels. In addition, three-dimensional predicted enzyme models were built and then validated by using various bioinformatics programs computationally. The interactions between enzyme and the substrate were analyzed via docking program in the presence of two substrates (xylotetraose = X4 and xylopentaose = X5). According to molecular phylogeny analysis, three clusters of these enzymes occurred: the first group had PzXyl, AnXyl, and CvXyl, and the second group possessed GaXyl and OmXyl, and the third group included ArXyl. Multiple sequence alignment analysis demonstrated that the five xylanases (ArXyl, OmXyl, CvXyl, PzXyl, AnXyl) had longer N-terminal regions, indicating greater thermal stability, relative to the GaXyl. Homology modeling showed that all the predicted model structures were, to a great extent, conserved. Docking analysis results indicated that CvXyl, OmXyl, and AnXyl had higher binding efficiency to two substrates, compared to the GaXyl, PzXyl, and ArXyl xylanases, and CvXyl-X4 docked complex had the highest substrate affinity with a binding energy of -9.8 kCal/mol. CvXyl, OmXyl, and AnXyl enzymes commonly had arginine in B8 β-strand interacted with two substrates, different from the other enzymes having lower binding efficiency. As a result, it was concluded that the three thermoacidophilic xylanase enzymes might be better candidates as the animal feed additive.

Toll-like Receptor 3 in the Hybrid Yellow Catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂): Protein Structure, Evolution and Immune Response to Exogenous Aeromonas hydrophila and Poly (I:C) Stimuli.

As a major mediator of cellular response to viral infection in mammals, Toll-like receptor 3 (TLR3) was proved to respond to double-stranded RNA (dsRNA). However, the molecular mechanism by which TLR3 functions in the viral infection response in teleosts remains to be investigated. In this study, the Toll-like receptor 3 gene of the hybrid yellow catfish was identified and characterized by comparative genomics. Furthermore, multiple sequence alignment, genomic synteny and phylogenetic analysis suggested that the homologous TLR3 genes were unique to teleosts. Gene structure analysis showed that five exons and four introns were common components of TLR3s in the 12 examined species, and interestingly the third exon in teleosts was the same length of 194 bp. Genomic synteny analysis indicated that TLR3s were highly conserved in various teleosts, with similar organizations of gene arrangement. De novo predictions showed that TLR3s were horseshoe-shaped in multiple taxa except for avian (with a round-shaped structure). Phylogenetic topology showed that the evolution of TLR3 was consistent with the evolution of the studied species. Selection analysis showed that the evolution rates of TLR3 proteins were usually higher than those of TLR3-TIR domains, indicating that the latter were more conserved. Tissue distribution analysis showed that TLR3s were widely distributed in the 12 tested tissues, with the highest transcriptions in liver and intestine. In addition, the transcription levels of TLR3 were significantly increased in immune-related tissues after infection of exogenous Aeromonas hydrophila and poly (I:C). Molecular docking showed that TLR3 in teleosts forms a complex with poly (I:C). In summary, our present results suggest that TLR3 is a pattern recognition receptor (PRR) gene in the immune response to pathogen infections in hybrid yellow catfish.

Molecular Docking Studies of Rifampicin - rpoB complex: Repurposing Drug Design Implications for against Plasmodium falciparum Malaria through a Computational Approach.

Malaria is one of the world's most devastating diseases, infecting well over 300 million people annually and killing between 2 and 3 million worldwide. Increasing parasite resistance to many existing drugs is exacerbating disease. Resistance to commonly used malarial drugs is increasing the need to develop new drugs urgently. Due to the slow pace and substantial costs of new drug development, repurposing of old drugs which is recently increasingly becoming an attractive proposition of highly efficient and effective way of drug discovery led us to study the drug rifampicin for this purpose. The present paper aims to investigate the route of Plasmodium falciparum apicoplast-targeted proteins that putatively encode β subunits of RNA polymerase with an objective to develop an effective antimalarial drug. Homology searching for conserved binding site to the rifampicin drug and the functional analysis of rpoB gene were done. Multiple Sequence alignment analysis of rpoB was compared with that in E.coli - rpoB and M. tuberculosis - rpoB. Docking studies of Rifampicin - rpoB complex was also done for finding binding affinity. The results of computational studies showed that rifampicin is a potential drug for malaria.

Cloning and sequence analysis of 4-O-α-d-isomaltooligosaccharylmaltooligosaccharide 1,4-α-isomaltooligosaccharohydrolase from Sarocladium kiliense U4520.

A novel enzyme, 4-O-α-d-isomaltooligosaccharylmaltooligosaccharide 1,4-α-isomaltooligosaccharohydrolase (IMM-4IH), was previously discovered from Sarocladium kiliense U4520. In order to identify the factors underlying the unique substrate specificity of IMM-4IH, we endeavored to determine the amino acid sequence of the enzyme. By comparing the partial amino acid sequence of the enzyme to whole genome sequencing data of S. kiliense U4520, the IMM-4IH gene was estimated. The putative gene was expressed in Pichia pastoris, and its activity and properties were found to be consistent with those of the native enzyme. Comparing the amino acid sequence of IMM-4IH with those in the CAZy database led to classification in the glycoside hydrolase family 49 (GH49). Several amino acids important for catalysis (Asp406, Asp425, and Asp426) and substrate recognition at subsites+1 and -3 were estimated by multiple sequence alignment analysis. These results provide important information for characterizing IMM-4IH and other GH49 enzymes.

Sequence Conservation Analysis and Gene Relationships of Nucleocapsid (N) Gene in Orthocoronavirinae Subfamily

Coronavirus (CoV) is a virus that causes respiratory and gastrointestinal diseases in animals and humans. It belongs to the Orthocoronavirina. The nucleocapsid protein (N) plays multiple roles in virus assembly, RNA transcription, and interaction with host cells. This study aimed to analyse the N protein by identifying conserved residues and exploring the gene and protein relationships within the Orthocoronavirinae. Therefore the results of this study are expected to help identify conserved regions of N protein in SARS-CoV-2 which can be used as probes for the virus identification process and can be used as target areas in vaccine development. We used 159 N gene and protein sequences, including 64 from Alpha, 51 from Beta-, 11 from Delta-, and 20 from Gammacoronavirus genera of the Orthocoronavirinae. Three sequences from Tobaniviridae were used as outgroups. Multiple sequence alignment (MSA) and phylogenetic tree analysis were performed using the neighbour-joining and Maximum Likelihood. The MSA results revealed several conserved residues, ranging from 18 to 41, were located in the N-terminal and Cterminal domains, the linker region, Nuclear Localization Signal (NLS), Nuclear Export Signal (NES) motifs, and Packing Signal (PS) binding sites. The phylogenetic tree analysis indicated that Gammacoronavirus and Deltacoronavirus were closely related to Betacoronavirus, while Alfacoronavirus showed the most distant relationship. Furthermore, the study identified 23 conserved residues involved in RNA binding, including amino acids such as Ser89, Val111, Pro112, Gly124, Tyr125, Phe150, Tyr151, Gly154, Thr155, Gly156, Trp180, Val181, Gly409, Arg411, Asn419, Gly421, and Pro443. These residues interacted with phosphate groups, nitrogenous bases, and pentose sugars and exhibited non-specific interactions with RNA. In summary, this study investigated the N protein in the Orthocoronavirinae subfamily, providing insights into its function, structure, and evolutionary relationships.

Molecular cloning, expression analysis and immune-related functional identification of tumor necrosis factor alpha (TNFα) in Sepiella japonica under bacteria stress

Tumor necrosis factor α (TNFα), a cytokine mainly secreted by active macrophages and monocytes, causes hemorrhagic necrosis of tumor tissues, kills tumor cells, regulates inflammatory responses, and plays a crucial role in innate immunity. In this study, TNFα of Sepiella japonica (named as SjTNFα) was acquired, whose full-length cDNA was 1206 bp (GenBank accession no. ON357428), containing a 5′ UTR of 185 bp, a 3′ UTR of 137 bp and an open reading frame (ORF) of 1002bp to encode a putative peptide of 333 amino acids for constructing the transmembrane domain and the cytoplasmic TNF domain. Its predicted pI was 8.69 and the theoretical molecular weight was 44.72 KDa. Multiple sequence alignment and phylogenetic analysis showed that SjTNFα had the highest homology to Octopus sinensis, they fell into a unified branch and further clustered with other animals. Real-time PCR indicated that SjTNFα was widely expressed in all subject tissues, including spleen, pancreas, gill, heart, brain, optic lobe, liver and intestine, and exhibited the highest in the liver and the lowest in the brain. The relative expression of SjTNFα varied at the developmental period of juvenile stage, pre-spawning and oviposition in the squid, with the highest in the liver at the juvenile stage and oviposition, and in the optic lobe of pre-spawning. After being infected with Vibrio parahaemolyticus and Aeromonas hydrophila, the expression of SjTNFα in liver and gill were both upregulated with time, and the highest expression appeared at 24 h and 8 h in liver for different infection, and at 4 h in gill consistently. Cell localization showed that SjTNFα distributed on membrane of HEK293 cells because it was a type II soluble transmembrane protein. When HEK293 cells were stimulated with LPS of different concentrations, the NF-κB pathway was activated in the nucleus and the corresponding mRNA was transferred through the intracellular signal transduction pathway, resulting in the synthesis and release of TNFα, which made the expression of SjTNFα was up-regulated obviously. These findings showed that SjTNFα might play an essential role in the defense of S. japonica against bacteria challenge, which contributed to the understanding of the intrinsic immune signaling pathway of Cephalopoda and the further study of host-pathogen interactions.

Cloning, characterization and expression analysis of caspase-8 genes from the golden pompano (Trachinotus ovatus)

Bacterial infection triggers apoptosis and immunity by regulating the expression of caspases. In this article, we report two caspase-8 paralogs, termed ToCas8 and ToCas8-like, from an economically important teleost golden pompano (Trachinotus ovatus). Caspase-8 is a classic cysteine protease that is essential for the initiation and execution of apoptosis and is a critical regulator in initiating death receptor-induced programmed cell death. The two proteins encoded by ToCas8 and ToCas8-like were characterized by multiple sequence alignment and protein structure analysis showing similarity to their counterparts in vertebrates and exhibiting typical features of the Caspase-8 protein family. Gene synteny analysis further confirmed that two isoforms of caspase-8 genes exactly existed in golden pompano. Phylogenetic analysis showed that the Caspase-8 proteins of the selected species in the study were clustered together, and Cas8 and Cas8-like were clustered into two different subgroups. In addition, ToCas8 and ToCas8-like were widely distributed and both were expressed at high levels in gill, kidney and intestinal tissues. Transcript levels of ToCas8 and ToCas8-like were significantly upregulated in response to LPS and Poly(I:C) challenges, suggesting that caspase-8 genes are involved in the immune response to pathogens. In summary, we have identified two isoforms of caspase-8 in golden pompano that are involved in the immune response. These findings will expand our understanding of the evolutionary history and functional characteristics of caspase-8 genes in vertebrates.

Principal Component and Structural Element Analysis Provide Insights into the Evolutionary Divergence of Conotoxins.

Predatory cone snails (Conus) developed a sophisticated neuropharmacological mechanism to capture prey, escape against other predators, and deter competitors. Their venom's remarkable specificity for various ion channels and receptors is an evolutionary feat attributable to the venom's variety of peptide components (conotoxins). However, what caused conotoxin divergence remains unclear and may be related to the role of prey shift. Principal component analysis revealed clustering events within diet subgroups indicating peptide sequence similarity patterns based on the prey they subdue. Molecular analyses using multiple sequence alignment and structural element analysis were conducted to observe the events at the molecular level that caused the subgrouping. Three distinct subgroups were identified. Results showed homologous regions and conserved residues within diet subgroups but divergent between other groups. We specified that these structural elements caused subgrouping in alpha conotoxins that may play a role in function specificity. In each diet subgroup, amino acid character, length of intervening amino acids between cysteine residues, and polypeptide length influenced subgrouping. This study provides molecular insights into the role of prey shift, specifically diet preference, in conotoxin divergence.

In-silico Characterization of ICE1 Transcriptional Factor in Cold Signalling Network in Mungbean (Vigna radiata L.)

Aim: Mungbean (Vigna radiata L.) is a leguminous crop with high rich of protein and susceptible to cold stress. In silico characterization of ICE1, a cold tolerant transcriptional factor in Vigna radiata and transcriptional regulatory network with other species of Vigna and with special refrence to legume crops.
 Methodology: Genomic sequences of Vigna radiata, V. angularis and V. unguculata and other leguminous crops were retrieved from NCBI (https://pubmed.ncbi.nlm.nih.gov) and used for multiple sequence alignment, Phylogenetic and comparative analysis of ICE1 protein and motif analysis.
 Results: The present study showed that Vigna radiata ICE1 (VrICE1) gene was closely related to VaICE1 gene (Vigna angularis) and then VuICE1 gene (Vigna unguiculata). The study also revealed that the bHLH domain region (247-298aa) was the low complexity region and co-localization signals of V.radiataICE1. There were a similarities and dissimilarities of the ICE1 protein isoforms in V.radiata. It was noted that VrICE1 had a average of hydropathicity (GRAVY) of -0.593 in one and -0.6 in other. An instability index of 64.58 and 64.90. VrICE1 had a higher percentage of non-polar amino acid content and more number of random coils followed by alpha helix in secondary structure. The study also revealed that there was major involvement of VrICE1 proteins in biological and molecular functions under cold stress.
 Conclusion: The bioinformatics tools help researchers in getting information with regard to the functional aspect of the gene with respect to cold tolerance. A comprehensive analysis of the different physical and chemical properties of Vigna species could help us to identify their diversified usefulness.

Molecular characterization of Newcastle disease virus vaccines in Nigeria.

Newcastle disease (ND) caused by ND virus (NDV) is a serious impediment to effective poultry production in developing countries such as Nigeria. Despite employing vaccination and other control measures to curtail this disease, its severe forms still persist. This study aimed to confirm the virus strains in the NDV vaccine brands commonly used in Nigeria. We employed reverse-transcription polymerase chain reaction (RT-PCR), sequencing, and sequence analysis to characterize NDV strains in four NDV vaccines commonly used in Nigeria. Fragments of 300 bp from NDV fusion genes from the vaccines were amplified. Polymerase chain reaction products were sequenced and analyzed using multiple sequence alignment and phylogenetic analyses to characterize the vaccine viruses as pathotypes. All the vaccines gave positive results, confirming the presence of NDV. Multiple sequence alignment and phylogenetic analyses revealed that two of the vaccines had the lentogenic pathotype, while the other two had the mesogenic or velogenic pathotype. This study provides information to facilitate strategies for regular control of the quality of vaccines in Nigeria.

Schizosaccharomyces pombe Grx4, Fep1, and Php4: In silico analysis and expression response to different iron concentrations.

Due to iron's essential role in cellular metabolism, most organisms must maintain their homeostasis. In this regard, the fission yeast Schizosaccharomyces pombe (sp) uses two transcription factors to regulate intracellular iron levels: spFep1 under iron-rich conditions and spPhp4 under iron-deficient conditions, which are controlled by spGrx4. However, bioinformatics analysis to understand the role of the spGrx4/spFep1/spPhp4 axis in maintaining iron homeostasis in S. pombe is still lacking. Our study aimed to perform bioinformatics analysis on S. pombe proteins and their sequence homologs in Aspergillus flavus (af), Saccharomyces cerevisiae (sc), and Homo sapiens (hs) to understand the role of spGrx4, spFep1, and spPhp4 in maintaining iron homeostasis. The three genes' expression patterns were also examined at various iron concentrations. A multiple sequence alignment analysis of spGrx4 and its sequence homologs revealed a conserved cysteine residue in each PF00085 domain. Blast results showed that hsGLRX3 is most similar to spGrx4. In addition, spFep1 is most closely related in sequence to scDal80, whereas scHap4 is most similar to spFep1. We also found two highly conserved motifs in spFep1 and its sequence homologs that are significant for iron transport systems because they contain residues involved in iron homeostasis. The scHap4 is most similar to spPhp4. Using STRING to analyze protein-protein interactions, we found that spGrx4 interacts strongly with spPhp4 and spFep1. Furthermore, spGrx4, spPhp4, and spFep1 interact with spPhp2, spPhp3, and spPhp5, indicating that the three proteins play cooperative roles in iron homeostasis. At the highest level of Fe, spgrx4 had the highest expression, followed by spfep1, while spphp4 had the lowest expression; a contrast occurred at the lowest level of Fe, where spgrx4 expression remained constant. Our findings support the notion that organisms develop diverse strategies to maintain iron homeostasis.

Pdif-mediated antibiotic resistance genes transfer in bacteria identified by pdifFinder.

Modules consisting of antibiotic resistance genes (ARGs) flanked by inverted repeat Xer-specific recombination sites were thought to be mobile genetic elements that promote horizontal transmission. Less frequently, the presence of mobile modules in plasmids, which facilitate a pdif-mediated ARGs transfer, has been reported. Here, numerous ARGs and toxin-antitoxin genes have been found in pdif site pairs. However, the mechanisms underlying this apparent genetic mobility is currently not understood, and the studies relating to pdif-mediated ARGs transfer onto most bacterial genera are lacking. We developed the web server pdifFinder based on an algorithm called PdifSM that allows the prediction of diverse pdif-ARGs modules in bacterial genomes. Using test set consisting of almost 32 thousand plasmids from 717 species, PdifSM identified 481 plasmids from various bacteria containing pdif sites with ARGs. We found 28-bp-long elements from different genera with clear base preferences. The data we obtained indicate that XerCD-dif site-specific recombination mechanism may have evolutionary adapted to facilitate the pdif-mediated ARGs transfer. Through multiple sequence alignment and evolutionary analyses of duplicated pdif-ARGs modules, we discovered that pdif sites allow an interspecies transfer of ARGs but also across different genera. Mutations in pdif sites generate diverse arrays of modules which mediate multidrug-resistance, as these contain variable numbers of diverse ARGs, insertion sequences and other functional genes. The identification of pdif-ARGs modules and studies focused on the mechanism of ARGs co-transfer will help us to understand and possibly allow controlling the spread of MDR bacteria in clinical settings. The pdifFinder code, standalone software package and description with tutorials are available at https://github.com/mjshao06/pdifFinder.

Tracing the Genetic Evolution of Canine Parvovirus Type 2 (CPV-2) in Thailand.

Canine parvovirus type 2 (CPV-2) is responsible for hemorrhagic gastroenteritis in dogs worldwide. High genomic substitution rates in CPV-2 contribute to the progressive emergence of novel variants with increased ability to evade the host immune response. Three studies have analyzed the genomic mutations of CPV-2 variants in Thailand. These investigations were independently conducted at different timepoints. Thus, a retrospective integrated analysis of CPV-2 genomic mutations has not been fully performed. Our study aimed at evaluating the evolutionary changes in CPV-2 in Thailand from 2003 to 2019. Two hundred and sixty-eight Thai CPV-2 nucleotide sequences were used for multiple amino acid sequence alignment and phylogenetic analyses. From 2003 to 2010, CPV-2a and -2b were the only variants detected. CPV-2c, emerged in 2014, replacing CPV-2a and -2b, and has become a major variant in 2019. Phylogenetic analysis revealed that the proposed mutation pattern of VP2 amino acid residues could help distinguish Thai CPV-2 variants. This comprehensive examination provides insight into the genomic evolution of CPV-2 in Thailand since its first reporting in 2003, which may facilitate the surveillance of the potential genetic alteration of emergent CPV-2 variants.