Search

MOTIVATION: The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths. RESULTS: For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings. AVAILABILITY AND IMPLEMENTATION: AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts).

Network alignment aims to uncover topologically similar regions in the protein-protein interaction (PPI) networks of two or more species under the assumption that topologically similar regions tend to perform similar functions. Although there exist a plethora of both network alignment algorithms and measures of topological similarity, currently no gold standard exists for evaluating how well either is able to uncover functionally similar regions. Here we propose a formal, mathematically and statistically rigorous method for evaluating the statistical significance of shared GO terms in a global, 1-to-1 alignment between two PPI networks. Given an alignment in which k aligned protein pairs share a particular GO term g, we use a combinatorial argument to precisely quantify the p-value of that alignment with respect to g compared to a random alignment. The p-value of the alignment with respect to all GO terms, including their inter-relationships, is approximated using the Empirical Browns Method. We note that, just as with BLASTs p-values, this method is not designed to guide an alignment algorithm towards a solution; instead, just as with BLAST, an alignment is guided by a scoring matrix or function; the p-values herein are computed after the fact, providing independent feedback to the user on the biological quality of the alignment that was generated by optimizing the scoring function. Importantly, we demonstrate that among all GO-based measures of network alignments, ours is the only one that correlates with the precision of GO annotation predictions, paving the way for network alignment-based protein function prediction.

In this study, we utilize Protein Residue Networks (PRNs), constructed using Local Spatial Pattern (LSP) alignment, to explore the dynamic behavior of Catabolite Activator Protein (CAP) upon the sequential binding of cAMP. We employed the Degree Centrality of these PRNs to investigate protein dynamics on a sub-nanosecond time scale, hypothesizing that it would reflect changes in CAPs entropy related to its thermal motions. We show that the binding of the first cAMP led to an increase in stability in the Cyclic-Nucleotide Binding Domain A (CNBD-A) and destabilization in CNBD-B, agreeing with previous reports explaining the negative cooperativity of cAMP binding in terms of an entropy-driven allostery. LSP-based PRNs also allow for the study of Betweenness Centrality, another graph-theoretical characteristic of PRNs, providing insights into global residue connectivity within CAP. Using this approach, we were able to correctly identify amino acids that were shown to be critical in mediating allosteric interactions in CAP. The agreement between our studies and previous experimental reports validates our method, particularly with respect to the reliability of Degree Centrality as a proxy for entropy related to protein thermal dynamics. Because LSP-based PRNs can be easily extended to include dynamics of small organic molecules, polynucleotides, or other allosteric proteins, the methods presented here mark a significant advancement in the field, positioning them as vital tools for a fast, cost-effective, and accurate analysis of entropy-driven allostery and identification of allosteric hotspots.

Swine acute diarrhoea syndrome coronavirus (SADS-CoV) causes vomiting, severe diarrhoea and death in newborn piglets. The spike (S) protein plays a crucial role in promoting virus invasion and inducing neutralizing antibody production. In this study, the extracellular region of the S protein was used as an immunogen to immunize BALB/c mice. After immunization, B cells were collected, fused with SP2/0 myeloma cells, cultured and subcloned, and a cell line capable of secreting neutralizing antibodies was obtained and named as 5D6. Additionally, it was determined that the 5D6 mAb could be used as the primary antibody for western blotting and indirect immunofluorescence assay (IFA) to detect SADS-CoV. Further studies indicated that the 5D6 mAb binds to the 136STSHAAD142 motif, which located in the N-terminal domain (NTD) of the spike protein. This result suggested that the NTD of the S protein can induce the production of neutralizing antibodies. Amino acid sequence alignment revealed that the epitope of the 5D6 mAb was conserved among SADS-CoV strains. This study helps elucidate the S protein function of SADS-CoV, and the 5D6 mAb may be used to develop diagnostic and treatment tools for detecting SADS-CoV infection. [ABSTRACT FROM AUTHOR]

The ability to predict how amino acid changes affect proteins has a wide range of applications including in disease variant classification and protein engineering. Many existing methods focus on learning from patterns found in either protein sequences or protein structures. Here, we present a method for integrating information from sequence and structure in a single model that we term SSEmb (Sequence Structure Embedding). SSEmb combines a graph representation for the protein structure with a transformer model for processing multiple sequence alignments. We show that by integrating both types of information we obtain a variant effect prediction model that is robust when sequence information is scarce. We also show that SSEmb learns embeddings of the sequence and structure that are useful for other downstream tasks such as to predict protein-protein binding sites. We envisage that SSEmb may be useful both for variant effect predictions and as a representation for learning to predict protein properties that depend on sequence and structure. SSEmb is a multi-modal machine learning model that predicts how changes in a protein's amino acid sequence affect its function by combining information from a multiple sequence alignment and the three-dimensional structure. [ABSTRACT FROM AUTHOR]

Background: The construction of a pangenome graph is a fundamental task in pangenomics. A natural theoretical question is how to formalize the computational problem of building an optimal pangenome graph, making explicit the underlying optimization criterion and the set of feasible solutions. Current approaches build a pangenome graph with some heuristics, without assuming some explicit optimization criteria. Thus it is unclear how a specific optimization criterion affects the graph topology and downstream analysis, like read mapping and variant calling. Results: In this paper, by leveraging the notion of maximal block in a Multiple Sequence Alignment (MSA), we reframe the pangenome graph construction problem as an exact cover problem on blocks called Minimum Weighted Block Cover (MWBC). Then we propose an Integer Linear Programming (ILP) formulation for the MWBC problem that allows us to study the most natural objective functions for building a graph. We provide an implementation of the ILP approach for solving the MWBC and we evaluate it on SARS-CoV-2 complete genomes, showing how different objective functions lead to pangenome graphs that have different properties, hinting that the specific downstream task can drive the graph construction phase. Conclusion: We show that a customized construction of a pangenome graph based on selecting objective functions has a direct impact on the resulting graphs. In particular, our formalization of the MWBC problem, based on finding an optimal subset of blocks covering an MSA, paves the way to novel practical approaches to graph representations of an MSA where the user can guide the construction. [ABSTRACT FROM AUTHOR]

Laccase is an enzyme that plays an important role in fungi, including lignin degradation, stress defense, and formation of fruiting bodies. Auricularia auricula-judae is a white-rot fungus in the Basidiomycota phylum, capable of delignifying wood. In this study, seven genes belonging to the laccase family were identified through de novo sequencing, containing Cu-Oxidase, Cu-Oxidase_2, and Cu-Oxidase_3 domains. Subsequently, the physical characteristics, phylogenetic relationships, protein secondary structure, and tertiary structure of the laccase family (AaLac1–AaLac7) were analyzed. Prediction of N-glycosylation sites identified 2 to 10 sites in the laccase family, with AaLac7 having the highest number of sites at 10. Sequence alignment and analysis of the laccase family showed high consistency in signature sequences. Phylogenetic analysis confirmed the relationship among laccases within the family, with AaLac3–AaLac4 and AaLac5–AaLac6 being closely positioned on the tree, exhibiting high similarity in tertiary structure predictions. This study identified and analyzed laccase family genes in Auricularia auricula-judae using de novo sequencing, offering a simple method for identifying and analyzing the laccase family in organisms with unknown genetic information. [ABSTRACT FROM AUTHOR]

Heat shock factor 2 (HSF2) is a versatile transcription factor that regulates gene expression under stress conditions, during development, and in disease. Despite recent advances in characterizing HSF2‐dependent target genes, little is known about the protein networks associated with this transcription factor. In this study, we performed co‐immunoprecipitation coupled with mass spectrometry analysis to identify the HSF2 interactome in mouse testes, where HSF2 is required for normal sperm development. Endogenous HSF2 was discovered to form a complex with several adhesion‐associated proteins, a finding substantiated by mass spectrometry analysis conducted in human prostate carcinoma PC‐3 cells. Notably, this group of proteins included the focal adhesion adapter protein talin‐1 (TLN1). Through co‐immunoprecipitation and proximity ligation assays, we demonstrate the conservation of the HSF2‐TLN1 interaction from mouse to human. Additionally, employing sequence alignment analyses, we uncovered a TLN1‐binding motif in the HSF2 C terminus that binds directly to multiple regions of TLN1 in vitro. We provide evidence that the 25 C‐terminal amino acids of HSF2, fused to EGFP, are sufficient to establish a protein complex with TLN1 and modify cell–cell adhesion in human cells. Importantly, this TLN1‐binding motif is absent in the C‐terminus of a closely related HSF family member, HSF1, which does not form a complex with TLN1. These results highlight the unique molecular characteristics of HSF2 in comparison to HSF1. Taken together, our data unveil the protein partners associated with HSF2 in a physiologically relevant context and identifies TLN1 as the first adhesion‐related HSF2‐interacting partner. [ABSTRACT FROM AUTHOR]

WRKY transcription factor (TF) family is well known to govern essential physiological functioning as well as regulate plant response to biotic and abiotic stress. In this study, we have identified 108 OsWRKY genes in the genome of Oryza sativa subsp. japonica, using the updated genomic data from the Rice Annotation Project Database and Oryzabase, which were further used to conduct the phylogenetic study, motif analysis, gene structure analysis, chromosomal mapping, and prediction of sub‐cellular localization. The multiple sequence alignment OsWRKY proteins revealed the presence of nine different types of alterations in the conserved heptapeptide sequence WRKYGQK associated with 19 OsWRKY genes. Physiochemical analysis discloses the hydrophobic amino acid‐rich, thermally stable, and polar nature of OsWRKY proteins. These genes were noted as highly conserved between the two cultivated sub‐species of Oryza sativa, that is Indica and japonica type. Additionally, from motif analysis, we have found a new motif, which was categorized as hAT family C‐terminal dimerization region associated with four members of group IIc. We have identified 21 stress‐responsive OsWRKY genes, and their significance to the different biotic and abiotic stress‐mediated cascades was further evaluated by analysing 1500 kb upstream sequences and this disclosed the presence of important phytohormone‐responsive cis‐elements in the OsWRKY gene, suggesting its direct involvement in defence against a wide range of external stressors and these 21 OsWRKY genes are tentatively listed as possible candidates for more study. [ABSTRACT FROM AUTHOR]

In this study, we have demonstrated a complementary-determining region (CDR) grafting technology for the generation of rabbit scFvs with different antigen recognition and physicochemical properties. The antigen-binding affinity of the CDR-grafted anti-CRP scFv, C1R/B1R (V1), which was generated by the CDR/framework region (CDR/FR) definition based on the traditional numbering rule, was insufficient when compared to that of the original clone, C1R, suggesting that the amino acid residues outside the original CDRs might significantly contribute to antigen recognition in rabbit scFvs. We redefined new CDRs and FRs to maintain antigen-binding affinities through the extension of multiple amino acid residues for CDRH1 and CDRH2, based on the amino acid sequence alignments of rabbit scFvs isolated from phage libraries. The new version successfully maintained the antigen binding affinity. CDR-grafted scFvs possessing a common CDR sequence and different FR sequences were successfully generated based on this new CDR/FR definition, and their physicochemical properties were further investigated. The antigen-binding activities of rabbit scFvs on Maxisorp varied between the tested clones in sandwich ELISA, supporting the idea that the combination of CDR with different FRs might change the physicochemical properties of scFvs on a solid material. The CDR-grafted scFvs possessing a frame sequence of anti-CRP scFv C2R maintained the ability to bind to protein L and were successfully purified. Expression titers showed improved solubility by diminishing the amount of insoluble scFvs. Thus, the method developed in this study is promising for generating alternatives with strict antigen binding recognition and different physicochemical properties. [Display omitted] [ABSTRACT FROM AUTHOR]

We introduce the volume entropy semi-norm and the systolic volume semi-norm in real homology and show that they satisfy functorial properties similar to the ones of the simplicial volume. Along the way, we also establish a roughly optimal upper bound on the systolic volume of the multiples of any homology class. Finally, we prove that the volume entropy semi-norm, the systolic volume semi-norm and the simplicial volume semi-norm are equivalent in every dimension. [ABSTRACT FROM AUTHOR]

Male sterility is a common phenomenon in higher plants and often plays an important role in the selection of superior offspring. 'Xiang Yun' is a mutant of Lagerstroemia indica that does not bear fruit after flowering, and its flowering period is significantly longer than that of normal L. indica. To explore the timing and molecular mechanisms of sterility in 'Xiang Yun', this study determined the period of sterility through anatomical observation and compared the content of nutrients and the activity of antioxidative enzymes at different stages of flower development. Finally, sequence alignment and qPCR were used to analyze the differences in pollen development genes between 'Xiang Yun' and 'Hong Ye'. The results showed that the anthers of 'Xiang Yun' dispersed pollen normally, but the pollen grains could not germinate normally. Observations with scanning electron microscopy revealed that the pollen grains were uneven in size and shriveled in shape. Further observation of anther sections found that abnormal development of the microspores began at the S2 stage, with the callose wall between microspores of 'Xiang Yun' being thicker than that of 'Hong Ye'. In addition, during the flower development of 'Xiang Yun', the contents of soluble sugar, soluble protein, free proline, and triglycerides were deficient to varying degrees, and the activities of POD, SOD, and MDA were lower. Sequence alignment and qPCR showed that there were several mutations in EFD1, TPD1, and DEX1 of 'Xiang Yun' compared with 'Hong Ye', and the expression levels of these genes were abnormally elevated in the later stages of development. Our results clarified the timing and phenotype of male sterility in 'Xiang Yun'. This provides solid and valuable information for further research on the molecular mechanism of sterility in 'Xiang Yun' and the genetic breeding of crape myrtle. [ABSTRACT FROM AUTHOR]

A novel multifunctional isoamylase, MIsA from Myxococcus sp. strain V11, was expressed in Escherichia coli BL21(DE3). Sequence alignment revealed that MIsA is a typical isoamylase that belongs to glycoside hydrolase family 13 (GH 13). MIsA can hydrolyze the α-1,6-branches of amylopectin and pullulan, as well as the α-1,4-glucosidic bond in amylose. Additionally, MIsA demonstrates 4-α-D-glucan transferase activity, enabling the transfer of α-1,4-glucan oligosaccharides between molecules, particularly with linear maltooligosaccharides. The Km, Kcat, and Vmax values of the MIsA for amylopectin were 1.22 mM, 40.42 µmol·min–1·mg–1, and 4046.31 mM·min–1. The yields of amylopectin and amylose hydrolyzed into oligosaccharides were 10.16% and 11.70%, respectively. The hydrolysis efficiencies were 55%, 35%, and 30% for amylopectin, soluble starch, and amylose, respectively. In the composite enzyme hydrolysis of amylose, the yield of maltotetraose increased by 1.81-fold and 2.73-fold compared with that of MIsA and MTHase (MCK8499120) alone, respectively. [ABSTRACT FROM AUTHOR]

Proteins play a pivotal role in biological systems. The use of machine learning algorithms for protein classification can assist and even guide biological experiments, offering crucial insights for biotechnological applications. We introduce the support bio-sequence machine for proteins (SBSM-Pro), a model purpose-built for the classification of biological sequences. This model starts with raw sequences and groups amino acids based on their physicochemical properties. It incorporates sequence alignment to measure the similarities between proteins and uses a novel multiple kernel learning (MKL) approach to integrate various types of information, utilizing support vector machines for classification prediction. The results indicate that our model demonstrates commendable performance across ten datasets in terms of the identification of protein function and post translational modification. This research not only exemplifies state-of-the-art work in protein classification but also paves avenues for new directions in this domain, representing a beneficial endeavor in the development of platforms tailored for the classification of biological sequences. SBSM-Pro is available for access at . [ABSTRACT FROM AUTHOR]

Simple Summary: To obtain more female individuals in the breeding population of R. dybowskii, it is important to explore the sex determination mechanism based on whole genome analysis. In this study, we analyzed the genome survey of male individuals of R. dybowskii using next-generation sequencing technology. The estimated genome size of R. dybowskii was approximately 3585.05 M. The total size of contigs and scaffolds was 3,748,543,415 and 3,765,862,278 bp, respectively. The N50 length of contigs and scaffolds was 31,988 and 336,385,783, respectively. The Dmrt1 gene with an 18,893 bp length was obtained from the scaffolding genomes of R. dybowskii. Sequence alignment to Dmrt1 revealed that the two male-specific DNA markers (MSM-222 and MSM-261) were not located on Dmrt1. This indicated that Dmrt1 may not be the only key gene for sex determination in R. dybowskii. Rana dybowskii is one of the important aquaculture species in Northeast China. The fallopian tubes of female R. dybowskii are used to prepare oviductus ranae (an important traditional Chinese medicine). Therefore, R. dybowskii females have higher economical value than males. An increasing female R. dybowskii population can increase the benefits from R. dybowskii culture. However, the genome of amphibians is complex, making it difficult to investigate their sex determination mechanism. In this study, we analyzed the genome of male R. dybowskii using next-generation sequencing technology. A total of 200,046,452,400 bp of clean data were obtained, and the K-mer analysis indicated that the depth was 50×. The genome size of R. dybowskii was approximately 3585.05 M, with a heterozygosity rate, repeat sequence ratio, and genome GC content of 1.15%, 68.96%, and approximately 43.0%, respectively. In total, 270,785 contigs and 498 scaffolds were generated. The size of the contigs and scaffolds was 3,748,543,415 and 3,765,862,278 bp, respectively, with the N50 length of 31,988 and 336,385,783. The longest contig and scaffold were of the size 137,967,485 and 1,808,367,828 bp, respectively. The number of contigs and scaffolds > 10K nt was 99,620 and 451, respectively. Through annotation, 40,913 genes were obtained, including 156,609 CDS (i.e., 3.83 CDS per gene). Sequence alignment was performed with the assembled scaffolding genome in this study. Two and one fragment had high homology with two male-specific DNA molecular markers of R. dybowskii discovered previously (namely, MSM-222 and MSM-261, respectively). In addition, the Dmrt1 gene of R. dybowskii was obtained with a length of 18,893 bp by comparison and splicing. The forward primers amplifying MSM-222 and MSM-261 were located at 322–343 and 14,501–14,526 bp of Dmrt1, respectively. However, sequence alignment revealed that MSM-222 and MSM-261 were not located on Dmrt1, and only some homologous parts were observed. This indicated that in addition to Dmrt1, other important genes may play a crucial role in the sex determination mechanism of R. dybowskii. Our study provided a foundation for the subsequent high-quality genome construction and provided important genomic resources for future studies on R. dybowskii. [ABSTRACT FROM AUTHOR]

Background: Multiple sequence alignment (MSA) has proven extremely useful in computational biology, especially in inferring evolutionary relationships via phylogenetic analysis and providing insight into protein structure and function. An alternative to the standard MSA model is partial order alignment (POA), in which aligned sequences are represented as paths in a graph rather than rows in a matrix. While the POA model has proven useful in several applications (e.g. sequencing reads assembly and pangenome structure exploration), we lack efficient visualization tools that could highlight its advantages. Results: We propose Sequence Flow – a web application designed to address the above problem. Sequence Flow presents the POA as a Sankey diagram, a kind of graph visualisation typically used for graphs representing flowcharts. Sequence Flow enables interactive alignment exploration, including fragment selection, highlighting a selected group of sequences, modification of the position of graph nodes, structure simplification etc. After adjustment, the visualization can be saved as a high-quality graphic file. Thanks to the use of SanKEY.js – a JavaScript library for creating Sankey diagrams, designed specifically to visualize POAs, Sequence Flow provides satisfactory performance even with large alignments. Conclusions: We provide Sankey diagram-based POA visualization tools for both end users (Sequence Flow) and bioinformatic software developers (SanKEY.js). Sequence Flow webservice is available at https://sequenceflow.mimuw.edu.pl/. The source code for SanKEY.js is available at https://github.com/Krzysiekzd/SanKEY.js and for Sequence Flow at https://github.com/Krzysiekzd/SequenceFlow. [ABSTRACT FROM AUTHOR]

Rapid Alkalization Factor (RALF) is a signaling molecule in plants that plays a crucial role in growth and development, reproductive processes, and responses to both biotic and abiotic stresses. Although RALF peptides have been characterized in Arabidopsis and rice, a comprehensive bioinformatics analysis of the ZmRALF gene family in maize is still lacking. In this study, we identified 20 RALF genes in the maize genome. Sequence alignment revealed significant structural variation among the ZmRALF family genes. Phylogenetic analysis indicates that RALF proteins from Arabidopsis, rice, and maize can be classified into four distinct clades. Duplication events suggest that the expansion of the RALF gene family in maize primarily relies on whole-genome duplication. ZmRALF genes are widely expressed across various tissues; ZmRALF1/15/18/19 are highly expressed in roots, while ZmRALF6/11/14/16 are predominantly expressed in anthers. RNA-seq and RT-qPCR demonstrated that the expression levels of ZmRALF7, ZmRALF9, and ZmRALF13 were significantly up-regulated and down-regulated in response to PEG and NaCl stresses, respectively. Overall, our study provides new insights into the role of the RALF gene family in abiotic stress. [ABSTRACT FROM AUTHOR]

Tetrapyrroles such as heme, chlorophyll, and vitamin B12 are essential for various metabolic pathways. They derive from 5-aminolevulinic acid (5-ALA), which can be synthesized by a single enzyme (5-ALA synthase or AlaS, Shemin pathway) or by a two-enzyme pathway. The genomes of some bacteriophages from aquatic environments carry various tetrapyrrole biosynthesis genes. Here, we analyze available metagenomic datasets and identify alaS homologs (viral alaS, or valaS) in sequences corresponding to marine and freshwater phages. The genes are found individually or as part of complete or truncated three-gene loci encoding heme-catabolizing enzymes. Amino-acid sequence alignments and three-dimensional structure prediction support that the valaS sequences likely encode functional enzymes. Indeed, we demonstrate that is the case for a freshwater phage valaS sequence, as it can complement an Escherichia coli 5-ALA auxotroph, and an E. coli strain overexpressing the gene converts the typical AlaS substrates glycine and succinyl-CoA into 5-ALA. Thus, our work identifies valaS as an auxiliary metabolic gene in phage sequences from aquatic environments, further supporting the importance of tetrapyrrole metabolism in bacteriophage biology. Tetrapyrroles such as heme and chlorophyll, which are essential for various metabolic pathways, are synthesized from 5-aminolevulinic acid (5-ALA). Here, the authors identify functional 5-ALA synthase genes in bacteriophage sequences from aquatic environments, thus highlighting the importance of tetrapyrrole metabolism in phage biology. [ABSTRACT FROM AUTHOR]

Pyocins are high molecular weight bacteriocins produced by Pseudomonas aeruginosa that can be retargeted to new bacterial species by exchanging the pyocin tail fibers with bacteriophage receptor binding proteins (RBPs). Here, we develop retargeted pyocins called campycins as new antibacterials to precisely and effectively kill the major foodborne pathogen Campylobacter jejuni. We used two diverse RBPs (H-fibers) encoded by CJIE1 prophages found in the genomes of C. jejuni strains CAMSA2147 and RM1221 to construct campycin 1 and campycin 2, respectively. Campycins 1 and 2 could target all C. jejuni strains tested due to complementary antibacterial spectra. In addition, both campycins led to more than 3 log reductions in C. jejuni counts under microaerobic conditions at 42 °C, whereas the killing efficiency was less efficient under anaerobic conditions at 5 °C. Furthermore, we discovered that both H-fibers used to construct the campycins bind to the essential major outer membrane protein (MOMP) present in all C. jejuni in a strain-specific manner. Protein sequence alignment and structural modeling suggest that the highly variable extracellular loops of MOMP form the binding sites of the diverse H-fibers. Further in silico analyses of 5000 MOMP sequences indicated that the protein falls into three major clades predicted to be targeted by either campycin 1 or campycin 2. Thus, campycins are promising antibacterials against C. jejuni and are expected to broadly target numerous strains of this human pathogen in nature and agriculture. Key points: • Campycins are engineered R-type pyocins containing H-fibers from C. jejuni prophages • Campycins reduce C. jejuni counts by >3 logs at conditions promoting growth • Campycins bind to the essential outer membrane protein MOMP in a strain-dependent way [ABSTRACT FROM AUTHOR]

Potassium (K+) is essential to plant growth and development. The absorption and transport of K+ by plant roots occurs mainly through K+ channels or transporters. In this study, an AKT1-type channel named PbrAKT1 was successfully cloned and characterised by a series of detailed procedures, including protein sequence alignment, expression and subcellular localisation analysis, as well as functional verification by heterologous expression function validation in yeast mutant R5421, Arabidopsis akt1 mutant, and Xenopus oocytes. The results showed that PbrAKT1 is preferentially expressed in pear roots and exhibits plasma membrane localisation, but its expression is unaffected by low K+ treatment. PbrAKT1 exhibits the function of an inwardly-rectifying K+ channel, and its activity is autonomously regulated when expressed in Xenopus oocytes. In addition, a series of heterologous expression experiments confirmed that PbrAKT1 plays an indispensable role in the processes of K+ transport, as well as alleviating the low-K+ -sensitive phenotype of yeast mutant R5421 and Arabidopsis akt1 mutant. In summary, PbrAKT1 plays an indispensable role in the processes of K+ transport and provide new information on the molecular mechanism of K+ absorption in pear roots. [ABSTRACT FROM AUTHOR]

Tight junctions play a pivotal role in the functional integrity of the human body by forming barriers that compartmentalize tissues and protect the body from external threats. Essential components of tight junctions are the transmembrane claudin proteins, which can polymerize into tight junction strands and meshworks. This study delves into the structural determinants of claudin polymerization, using the close homology yet strong difference in polymerization capacity between claudin-3 and claudin-4. Through a combination of sequence alignment and structural modeling, critical residues in the second extracellular segment are pinpointed. Molecular dynamics simulations provide insights into the interactions of and the conformational changes induced by the identified extracellular segment 2 residues. Live-stimulated emission depletion imaging demonstrates that introduction of these residues from claudin-3 into claudin-4 significantly enhances polymerization in nonepithelial cells. In tight junction-deficient epithelial cells, mutated claudin-4 not only influences tight junction morphology but also partially restores barrier function. Understanding the structural basis of claudin polymerization is crucial, as it offers insights into the dynamic nature of tight junctions. This knowledge could be applied to targeted therapeutic interventions, offer insight to repair or prevent barrier defects associated with pathological conditions, or introduce temporary barrier openings during drug delivery. [ABSTRACT FROM AUTHOR]

The water‐lily clade represents the second earliest‐diverging branch of angiosperms. Most of its species belong to Nymphaeaceae, of which the "core Nymphaeaceae"—comprising the genera Euryale, Nymphaea and Victoria—is the most diverse clade. Despite previous molecular phylogenetic studies on the core Nymphaeaceae, various aspects of their evolutionary relationships have remained unresolved. The length‐variable introns and intergenic spacers are known to contain most of the sequence variability within the water‐lily plastomes. Despite the challenges with multiple sequence alignment, any new molecular phylogenetic investigation on the core Nymphaeaceae should focus on these noncoding plastome regions. For example, a new plastid phylogenomic study on the core Nymphaeaceae should generate DNA sequence alignments of all plastid introns and intergenic spacers based on the principle of conserved sequence motifs. In this investigation, we revisit the phylogenetic history of the core Nymphaeaceae by employing such an approach. Specifically, we use a plastid phylogenomic analysis strategy in which all coding and noncoding partitions are separated and then undergo software‐driven DNA sequence alignment, followed by a motif‐based alignment inspection and adjustment. This approach allows us to increase the reliability of the character base compared to the default practice of aligning complete plastomes through software algorithms alone. Our approach produces significantly different phylogenetic tree reconstructions for several of the plastome regions under study. The results of these reconstructions underscore that Nymphaea is paraphyletic in its current circumscription, that each of the five subgenera of Nymphaea is monophyletic, and that the subgenus Nymphaea is sister to all other subgenera of Nymphaea. Our results also clarify many evolutionary relationships within the Nymphaea subgenera Brachyceras, Hydrocallis and Nymphaea. In closing, we discuss whether the phylogenetic reconstructions obtained through our motif‐based alignment adjustments are in line with morphological evidence on water‐lily evolution. [ABSTRACT FROM AUTHOR]

Soybean [Glycine max (L.) Merr.], as a globally commercialized crop, is an important source of protein and oil for both humans and livestock. With more frequent extreme weather disasters, abiotic stress has become one of the critical factors restricting soybean production. Proline (Pro) is a well-known substance in plants that responds to abiotic stress. To identify potential effector genes involved in soybean resistance to abiotic stress, we focused on the pyrroline-5-carboxylate dehydrogenase (P5CDH) which is a key enzyme in the degradation process of Pro. Through homologous sequence alignment, phylogenetic tree, and predicted expression, we chose GmP5CDH1 (Glyma.05G029200) for further research. Tissue-specific expression assay showed that GmP5CDH1 had higher expression levels in soybean seed and cotyledon development. Subcellular localization assay revealed that GmP5CDH1 was a nuclear-membrane-localized protein. As the result of the predicted cis-acting regulatory element indicates, the expression level of GmP5CDH1 was induced by low temperature, drought, salt stress, and ABA in soybean. Next, we constructed transgenic Arabidopsis overexpressing GmP5CDH1. The results showed that GmP5CDH1 also strongly responded to exogenous Pro, and overcame the toxicity of abiotic stress on plants by regulating the endogenous concentration of Pro. The interaction between GmP5CDH1 and GmSAM1 was validated through yeast two-hybrid, LUC fluorescence complementary, and BIFC. In conclusion, overexpression of a soybean pyrroline-5-carboxylate dehydrogenase GmP5CDH1 regulates the development of Arabidopsis thaliana by altering proline content dynamically under salt stress, especially improving the growth of plants under exogenous Pro. [ABSTRACT FROM AUTHOR]

Soil salinity poses a serious threat to rice production. The salt tolerance of rice at the germination stage is one of the major determinants of stable stand establishment, which is very important for direct seeding in saline soil. The complexity and polygenic nature of salt tolerance have limited the efficiency of discovering and cloning key genes in rice. In this study, an RIL population with an ultra-high-density genetic map was employed to investigate the salt-tolerant genetic basis in rice, and a total of 20 QTLs were detected, including a major and stable QTL (qRCL3-1). Subsequently, salt-specific DEGs from a comparative transcriptome analysis were overlaid onto annotated genes located on a stable QTL interval, and eight putative candidate genes were further identified. Finally, from the sequence alignment and variant analysis, OsCam1-1 was confirmed to be the most promising candidate gene for regulating salinity tolerance in rice. This study provides important information for elucidating the genetic and molecular basis of rice salt tolerance at the germination stage, and the genes detected here will be useful for improvements in rice salt tolerance. [ABSTRACT FROM AUTHOR]

Dihydrofolate reductase (DHFR) is ubiquitously present in all living organisms and plays a crucial role in the growth of the fungal pathogen R.solani. Sequence alignment confirmed the evolutionary conservation of the essential lid domain, with the amino acid 'P' within the PEKN lid domain appearing with a frequency of 89.5% in higher organisms and 11.8% in lower organisms. Consequently, a K65P variant was introduced into R.solani DHFR (rDHFR). Subsequent enzymatic kinetics assays were conducted for human DHFR (hDHFR), rDHFR, E. coli DHFR (eDHFR), and the K65P variant. hDHFR exhibited the highest kcat of 0.95 s−1, followed by rDHFR with 0.14 s−1, while eDHFR displayed the lowest kcat of 0.09 s−1. Remarkably, the K65P variant induced a significant reduction in Km, resulting in a 1.8-fold enhancement in catalytic efficiency (kcat/Km) relative to the wild type. Differential scanning fluorimetry and binding free energy calculations confirmed the enhanced substrate affinity for both folate and NADPH in the K65P variant. These results suggest that the K65P mutation enhances substrate affinity and catalytic efficiency in DHFR, highlighting the evolutionary and functional importance of the K65 residue. [ABSTRACT FROM AUTHOR]

DNA barcoding is an effective identification tool for tree ferns with heteromorphic generations with morphologically simple gametophytic phase. A complete reference dataset which includes coverage of the target local flora is necessary for accurate identification. Three different Cyathea sp.pl. were investigated to test the utility of plastid DNA barcode regions rbcL with the intention of developing an identification system for native pteridophytes. In addition, the inter-specific and inter-generic variations among the Cyathea sp.pl. were also focused. BLAST analysis revealed the similar rbcL sequences related to C. nilgirensis Holttum, C. gigantea (Wallich ex Hook. f.) Holttum and C. crinita (Hook.) Copel and a sum of 143 sequences viz., Cyathea sp.pl. (121), Gymnosphaera sp.pl.(4), Sphaeropteris sp.pl. (1) and Alsophila sp.pl. (17). The sequences of 143 plants were retrieved from GenBank in FASTA format. The sequences were aligned using multiple sequence alignment tool MUSCLE. The inter-specific and inter-generic distance percentage were calculated using Maximum Likelihood, Minimum Evolution, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) and Neighbor-Joining method. The analyses were constructed using MEGA X software. The rbcL gene provided sequence variation with its strong resolving power. This study demonstrates the overall effectiveness of DNA barcodes for species identification of tree fern genus Cyathea in the pteridophyte flora of Southern Western Ghats. [ABSTRACT FROM AUTHOR]

SUMMARY: Structural prediction by artificial intelligence can be powerful new instruments to discover novel protein–protein interactions, but the community still grapples with the implementation, opportunities and limitations. Here, we discuss and re‐analyse our in silico screen for novel pathogen‐secreted inhibitors of immune hydrolases to illustrate the power and limitations of structural predictions. We discuss strategies of curating sequences, including controls, and reusing sequence alignments and highlight important limitations caused by different platforms, sequence depth and computing times. We hope these experiences will support similar interactomic screens by the research community. [ABSTRACT FROM AUTHOR]

Allantus luctifer is one of the most serious pests of buckwheat, with its larvae mainly damaging leaves during the seedling and flowering stages. Entomogenous fungi play a crucial role as biological regulators of arthropod populations in nature. In this paper, a newly isolated strain 19GZAl-1, was identified as Cordyceps. cateniannulata through the observetions of macroscopic and microscopic morphological features, and the results of rDNA- ITS sequence alignment and phylogenetic analysis. By comparing the efficacy of ten entomopathogenic fungal strains, including six strains of C. cateniannulata and four other Cordyceps species, against A. luctifer larvae, the new strain 19GZAl-1 exhibited the highest mortality rate. Mortality rates of A. luctifer larvae treated by spraying with the strain 19GZAl-1 increased as larval stage decreased, with first instar larvae showing the highest mortality rate of 85%. The appressorium from conidia invaded the larval body through areas with thin cuticle, such as larva pod base and internode folds, after which the hyphae grew rapidly, filling the haemocoel, and finally protruding from the integument upon melanization and decomposition of the intestinal wall cells. The results suggest that the strain 19GZAL-1 of C. cateniannulata has potential against A. luctifer larvae, which is significant for further study on the infection mechanism of C. cateniannulata on insects. [ABSTRACT FROM AUTHOR]

Automatic assessment of the operation ability of operators based on computers is an essential approach for improving the effectiveness of equipment operation training and enhancing equipment safety. Present methods primarily focus on the operation results but pay less attention to the operation procedure. One reason is that there is a lack of a model that has the ability to describe all probable paths to accomplish the same task. Therefore, an operations chain model is put forward for the first time to describe the standard operation procedure and relationships among operations based on the decomposition of operational tasks and the relationships among the various operations required to fulfill the task. A specific operation task corresponds to an operations chain, which will form one or multiple standard operation sequences that will allow trainees to complete the same task through different paths. The Needleman–Wunsch sequence alignment algorithm is introduced to match the trainees' operation sequence with all standard sequences. The maximum alignment result is the score of the trainees' operations. An example shows that the operations chain model can accurately describe the complex structure of the standard operating procedures. The Needleman–Wunsch sequence alignment algorithm can objectively evaluate the trainee's operation capabilities. The combination of the operations chain model and sequence alignment algorithm can form a complete operation procedure assessment method that is friendlier to trainees and has more objective evaluation results. The method will help to improve the effectiveness of the competency assessment of equipment operators. [ABSTRACT FROM AUTHOR]

Objectives: This study aimed to discuss the essential amino acid residues and catalytic mechanism of trans-epoxysuccinate hydrolase from Pseudomonas koreensis for the production of meso-tartaric acid. Results: The optimum conditions of the enzyme were 45 °C and pH 9.0, respectively. It was strongly inhibited by Zn2+, Mn2+ and SDS. Michaelis–Menten enzyme kinetics analysis gave a Km value of 3.50 mM and a kcat of 99.75 s−1, with an exceptional EE value exceeding 99.9%. Multiple sequence alignment and homology modeling revealed that the enzyme belonged to MhpC superfamily and possessed a typical α/β hydrolase folding structure. Site-directed mutagenesis indicated H34, D104, R105, R108, D128, Y147, H149, W150, Y211, and H272 were important catalytic residues. The 18O-labeling study suggested the enzyme acted via two-step catalytic mechanism. Conclusions: The structure and catalytic mechanism of trans-epoxysuccinate hydrolase were first reported. Ten residues were critical for its catalysis and a two-step mechanism by an Asp-His-Asp catalytic triad was proposed. [ABSTRACT FROM AUTHOR]

The rise of Candida auris as a multidrug-resistant pathogen presents considerable challenges in the clinical environment, requiring the identification of new therapeutic targets. This study aimed to annotate a hypothetical gene product, Set5 (Su(var)3-9, Enhancer-of-zeste, Trithorax) lysine methyltransferase, in Candida auris using sequence and structure analysis. Utilizing advanced bioinformatics tools, including sequence alignment, phylogenetic analysis and structural modeling, we sought to provide a thorough annotation of this gene product. Our methodology involved identifying conserved motifs characteristic of SET and MYND (myeloid translocation protein, Nervy, Deaf, SMYD) methyltransferases and confirming the key catalytic and substrate-binding domains through structural modeling. The results showed significant sequence and structural similarities with known SMYD family members, supporting the annotation of the hypothetical gene product as Set5 methyltransferase. This discovery not only enhances our understanding of Candida auris biology, but also identifies potential therapeutic targets for combating infections caused by this pathogen. The identification of SMYD in Candida auris could lead to new antifungal targets, contributing to global efforts to manage and treat infections caused by this challenging pathogen. [ABSTRACT FROM AUTHOR]

Statistical analysis of multiple sequence alignments of homologous proteins has revealed groups of coevolving amino acids called sectors. These groups of amino-acid sites feature collective correlations in their amino-acid usage, and they are associated to functional properties. Modeling showed that nonlinear selection on an additive functional trait of a protein is generically expected to give rise to a functional sector. These modeling results motivated a principled method, called ICOD, which is designed to identify functional sectors, as well as mutational effects, from sequence data. However, a challenge for all methods aiming to identify sectors from multiple sequence alignments is that correlations in amino-acid usage can also arise from the mere fact that homologous sequences share common ancestry, i.e. from phylogeny. Here, we generate controlled synthetic data from a minimal model comprising both phylogeny and functional sectors. We use this data to dissect the impact of phylogeny on sector identification and on mutational effect inference by different methods. We find that ICOD is most robust to phylogeny, but that conservation is also quite robust. Next, we consider natural multiple sequence alignments of protein families for which deep mutational scan experimental data is available. We show that in this natural data, conservation and ICOD best identify sites with strong functional roles, in agreement with our results on synthetic data. Importantly, these two methods have different premises, since they respectively focus on conservation and on correlations. Thus, their joint use can reveal complementary information. Author summary: Proteins perform crucial functions in the cell. The biological function of a protein is encoded in its amino-acid sequence. Natural selection acts at the level of function, while mutations arise randomly on sequences. In alignments of sequences of homologous proteins, which share common ancestry and common function, the amino acid usages at different sites can be correlated due to functional constraints. In particular, groups of collectively correlated amino acids, termed sectors, tend to emerge due to selection on functional traits. However, correlations can also arise from the shared evolutionary history of homologous proteins, even without functional constraints. This may obscure the inference of functional sectors. By analyzing controlled synthetic data as well as natural protein sequence data, we show that two very different methods allow to identify sectors and mutational effects in a way that is most robust to phylogeny. We suggest that considering both of these methods allows a better identification of functionally important sites from protein sequences. These results have potential impact on the design of new functional sequences. [ABSTRACT FROM AUTHOR]

[Objective] To identify samples of suspected tumor disease from a local chicken farm in Anhui Province. [Method] The chicken liver was collected and the virus DNA was extracted. The primers were designed according to the sequence of ALV and MDV that GenBank had registered. The ALV env and MDV meq gene were detected by PCR. The amplified sequence was sequenced and analyzed by evolutionary tree and nucleotide homology. [Result] Detected bands of ALV and MD, and sequence alignment results showed that ALV belongs to subgroup J, while MD is a strong wild virus strain. [Conclusion] This study provided reference data for understanding the mixed infection situation of ALV and MDV in the area. [ABSTRACT FROM AUTHOR]

Bioinformatic analysis revealed that the genomes of ubiquitous Penicillium spp. might carry dozens of biosynthetic gene clusters (BGCs), yet many clusters have remained uncharacterized. In this study, a detailed investigation of co-culture fermentation including the basidiomycete Armillaria mellea CPCC 400891 and the P. brasilianum CGMCC 3.4402 enabled the isolation of five new compounds including two bisabolene-type sesquiterpenes (arpenibisabolanes A and B), two carotane-type sesquiterpenes (arpenicarotanes A and B), and one polyketide (arpenichorismite A) along with seven known compounds. The assignments of their structures were deduced by the extensive analyses of detailed spectroscopic data, electronic circular dichroism spectra, together with delimitation of the biogenesis. Most new compounds were not detected in monocultures under the same fermentation conditions. Arpenibisabolane A represents the first example of a 6/5-fused bicyclic bisabolene. The bioassay of these five new compounds exhibited no cytotoxic activities in vitro against three human cancer cell lines (A549, MCF-7, and HepG2). Moreover, sequence alignments and bioinformatic analysis to other metabolic pathways, two BGCs including Pb-bis and Pb-car, responsible for generating sesquiterpenoids from co-culture were identified, respectively. Furthermore, based on the chemical structures and deduced gene functions of the two clusters, a hypothetic metabolic pathway for biosynthesizing induced sesquiterpenoids was proposed. These results demonstrated that the co-culture approach would facilitate bioprospecting for new metabolites even from the well-studied microbes. Our findings would provide opportunities for further understanding of the biosynthesis of intriguing sesquiterpenoids via metabolic engineering strategies. Key points: • Penicillium and Armillaria co-culture facilitates the production of diverse secondary metabolites • Arpenibisabolane A represents the first example of 6/5-fused bicyclic bisabolenes • A hypothetic metabolic pathway for biosynthesizing induced sesquiterpenoids was proposed [ABSTRACT FROM AUTHOR]

With the rapid development of biotechnology, gene sequencing methods are gradually improved. The structure of gene sequences is also more complex. However, the traditional sequence alignment method is difficult to deal with the complex gene sequence alignment work. In order to improve the efficiency of gene sequence analysis, D2 series method of k-mer statistics is selected to build the model of gene sequence alignment analysis. According to the structure of the foreground sequence, the sequence to be aligned can be cut by different lengths and divided into multiple subsequences. Finally, according to the selected subsequences, the maximum dissimilarity in the alignment results is determined as the statistical result. At the same time, the research also designed an application system for the sequence alignment analysis of the model. The experimental results showed that the statistical power of the sequence alignment analysis model was directly proportional to the sequence coverage and cutting length, and inversely proportional to the K value and module length. At the same time, the model was applied to the system designed in this paper. The maximum storage capacity of the system was 71 GB, the maximum disk capacity was 135 GB, and the running time was less than 2.0s. Therefore, the k-mer statistic sequence alignment model and system proposed in this study have considerable application value in gene alignment analysis. [ABSTRACT FROM AUTHOR]

Background: Novel protein-coding genes were considered to be born by re-organization of pre-existing genes, such as gene duplication and gene fusion. However, recent progress of genome research revealed that more protein-coding genes than expected were born de novo, that is, gene origination by accumulating mutations in non-genic DNA sequences. Nonetheless, the in-depth process (scenario) for de novo origination is not well understood. Results: We have conceived bioinformatic analysis for sketching a scenario for de novo origination of protein-coding genes. For each de novo protein-coding gene, we firstly identified an edge of a given phylogenetic tree where the gene was born based on parsimony. Then, from a multiple sequence alignment of the de novo gene and its orthologous regions, we constructed ancestral DNA sequences of the gene corresponding to both end nodes of the edge. We finally revealed statistical features observed in evolution between the two ancestral sequences. In the analysis of the Saccharomyces cerevisiae lineage, we have successfully sketched a putative scenario for de novo origination of protein-coding genes. (1) In the beginning was GC-rich genome regions. (2) Neutral mutations were accumulated in the regions. (3) ORFs were extended/combined, and then (4) translation signature (Kozak consensus sequence) was recruited. Interestingly, as the scenario progresses from (2) to (4), the specificity of mutations increases. Conclusion: To the best of our knowledge, this is the first report outlining a scenario of de novo origination of protein-coding genes. Our bioinformatic analysis can capture events that occur during a short evolutionary time by directly observing the evolution of the ancestral sequences from non-genic to genic. This property is suitable for the analysis of fast evolving de novo genes. [ABSTRACT FROM AUTHOR]

Background: Cannabis sativa cultivars can be classified as marijuana or hemp, depending on its amount of the psychoactive cannabinoid Δ9‐tetrahydrocannabinolic acid (THCA). Hemp Cheungsam is a non-drug type Cannabis sativa that is characterized by low THCA content. However, the transcripts and expression profile of cannabinoid biosynthesis pathway genes of hemp Cheungsam have not been investigated. Methods: RNA-sequencing (RNA-seq) was performed on three different tissue types (flower, leaf, and stem) of hemp Cheungsam to understand their transcriptomes. The expression of cannabinoid biosynthesis pathway genes was further analyzed in each tissue type. Multiple sequence alignment and conserved domain analyses were used to investigate the homologs of cannbinoid biosynthesis genes. Results: We found that the cannabinoid biosynthesis pathway was mainly expressed in the flowers of hemp Cheungsam, similar to other Cannabis cultivars. However, expression of cannabidiolic acid (CBDA) synthase was much higher than THCA synthase and cannabichromenic acid (CBCA) synthase, suggesting that the transcription profile favors CBDA biosynthesis. Sequence analysis of cannabinoid biosynthesis pathway genes suggested the identity of orthologs in hemp Cheungsam. Conclusions: Cannabinoid biosynthesis in hemp Cheungsam mostly occurs in the flowers, compared to other plant organs. While CBDA synthase expression is high, THCA and CBCA synthase expression is considerably low, indicating lesser THCA biosynthesis and thus low THCA content. Sequence analysis of key genes (CBDA, THCA, and CBCA synthases) of the cannabinoid biosynthetic pathway indicates that orthologs are present in hemp Cheungsam. [ABSTRACT FROM AUTHOR]

Models have always been central to inferring molecular evolution and to reconstructing phylogenetic trees. Their use typically involves the development of a mechanistic framework reflecting our understanding of the underlying biological processes, such as nucleotide substitutions, and the estimation of model parameters by maximum likelihood or Bayesian inference. However, deriving and optimizing the likelihood of the data is not always possible under complex evolutionary scenarios or even tractable for large datasets, often leading to unrealistic simplifying assumptions in the fitted models. To overcome this issue, we coupled stochastic simulations of genome evolution with a new supervised deep-learning model to infer key parameters of molecular evolution. Our model is designed to directly analyze multiple sequence alignments and estimate per-site evolutionary rates and divergence without requiring a known phylogenetic tree. The accuracy of our predictions matched that of likelihood-based phylogenetic inference when rate heterogeneity followed a simple gamma distribution, but it strongly exceeded it under more complex patterns of rate variation, such as codon models. Our approach is highly scalable and can be efficiently applied to genomic data, as we showed on a dataset of 26 million nucleotides from the clownfish clade. Our simulations also showed that the integration of per-site rates obtained by deep learning within a Bayesian framework led to significantly more accurate phylogenetic inference, particularly with respect to the estimated branch lengths. We thus propose that future advancements in phylogenetic analysis will benefit from a semi-supervised learning approach that combines deep-learning estimation of substitution rates, which allows for more flexible models of rate variation, and probabilistic inference of the phylogenetic tree, which guarantees interpretability and a rigorous assessment of statistical support. [ABSTRACT FROM AUTHOR]

Gelatin and collagen are two animal-derived ingredients that are widely used in various industries. Both have distinctive physico-chemical characteristic that made them ingredients of interest for many industrial players to be applied as there are vast arrays of usage in the food, cosmetic and biomedical fields. However, the origin of gelatin and collagen poses ethical and religious concerns, especially for Muslims and Jews who have restrictions on food consumption. Porcine by-products are of concern for religious and health reasons, and there is a demand for precise and reliable detection techniques. The limitation of DNA detection is due to extreme environment in food processing which results in low extractability of DNA. Therefore, peptide-based detection using mass spectrometry is required. However, identify the suitable marker is like searching needle in haystacks. Hence, combination of bioinformatics and mass spectrometry is proposed. This study aims to identify the specific peptide biomarkers by employing bioinformatics technique which can be applied to identify gelatin and collagen sources with the aid of mass spectrometry. In these approach, combination of Petunia Trans-Proteomic Pipeline (TPP, version 5.2.0) and sequence alignment ClustalW were applied to facilitate the MS data (LC-QTOF-MS) and peptide identification. As a result, 69 fasta file of protein sequence from both UniProtKB and NCBInr have been collected, 81 collagen peptides sequence and 118 gelatine peptides has been attainable that have the potential to distinguish different species. In conclusion, in silico protein sequence approaches helps to enable rapid screening of proteotypic peptides that can serve as species biomarkers proficiently. [ABSTRACT FROM AUTHOR]

Levan is a highly water-soluble biopolymer that has extensive applications in the pharmaceutical, personal care, food, and industrial fields. Lactobacillus fermentum SHN1 was isolated from milk and dairy products as an exopolysaccharide (levan) producer and identified according to morphological, biochemical tests and using the 16S rRNA gene sequencing. The sequence alignment in the gene bank indicated that the isolate has a high percentage of similarity (100%) to the recovered sequence of Lactobacillus fermentum isolate. The produced biopolymer was characterized as levan by FTIR, ¹H-NMR, and 13C-NMR spectroscopy. The effect of various nutritional and physical factors including pH, carbon source, sucrose concentration, nitrogen source, inoculum size, and incubation period, on the synthesis of levan by L. fermentum SHN1 was studied. The findings indicated that the optimal carbon source for the production of levan was 100 g/L sucrose, yielded of 25.95 g/L. Additionally, a combination of peptone + yeast extract at a ratio of 2:0.5 was identified as the most effective nitrogen source for levan production, giving in a yield of 24.21 g/L. L. fermentum SHN1 produces the greatest quantity of levan, 23.89 g/L, under a pH of 6.5 and an inoculum size of 1%, resulting in a yield of 25 g/L. The most favorable period for incubation was 48 hours, resulting in a production yield of 24.87 g/L. [ABSTRACT FROM AUTHOR]

There is a collective push to diversify human genetic studies by including underrepresented populations. However, analyzing DNA sequence reads involves the initial step of aligning the reads to the GRCh38/hg38 reference genome which is inadequate for non-European ancestries. In this study, using long-read sequencing technology, we constructed de novo genome assemblies from two indigenous Americans from Arizona (IAZ). Each assembly included ∼17 Mb of DNA sequence not present [nonreference sequence (NRS)] in hg38, which consists mostly of repeat elements. Forty NRSs totaling 240 kb were uniquely anchored to the hg38 primary assembly generating a modified hg38-NRS reference genome. DNA sequence alignment and variant calling were then conducted with whole-genome sequencing (WGS) sequencing data from 387 IAZ using both the hg38 and modified hg38-NRS reference maps. Variant calling with the hg38-NRS map identified ∼50,000 single-nucleotide variants present in at least 5% of the WGS samples which were not detected with the hg38 reference map. We also directly assessed the NRSs positioned within genes. Seventeen NRSs anchored to regions including an identical 187 bp NRS found in both de novo assemblies. The NRS is located in HCN2 79 bp downstream of Exon 3 and contains several putative transcriptional regulatory elements. Genotyping of the HCN2 -NRS revealed that the insertion is enriched in IAZ (minor allele frequency = 0.45) compared to other reference populations tested. This study shows that inclusion of population-specific NRSs can dramatically change the variant profile in an underrepresented ethnic groups and thereby lead to the discovery of previously missed common variations. [ABSTRACT FROM AUTHOR]

Computer-aided analysis of proteins or nucleic acids seems like a matter of course nowadays; however, the history of Bioinformatics and Computational Biology is quite recent. The advent of high-throughput sequencing has led to the production of "big data", which has also affected the field of virology. The collaboration between the communities of bioinformaticians and virologists already started a few decades ago and it was strongly enhanced by the recent SARS-CoV-2 pandemics. In this article, which is the first in a series on how bioinformatics can enhance virus research, we show that highly useful information is retrievable from selected general and dedicated databases. Indeed, an enormous amount of information—both in terms of nucleotide/protein sequences and their annotation—is deposited in the general databases of international organisations participating in the International Nucleotide Sequence Database Collaboration (INSDC). However, more and more virus-specific databases have been established and are progressively enriched with the contents and features reported in this article. Since viruses are intracellular obligate parasites, a special focus is given to host-pathogen protein-protein interaction databases. Finally, we illustrate several phylogenetic and phylodynamic tools, combining information on algorithms and features with practical information on how to use them and case studies that validate their usefulness. Databases and tools for functional inference will be covered in the next article of this series: Bioinformatics goes viral: II. Sequence-based and structure-based functional analyses for boosting virus research. [ABSTRACT FROM AUTHOR]

Aluminum-activated malate transporter (ALMT) genes play an important role in aluminum ion (Al3+) tolerance, fruit acidity, and stomatal movement. Although decades of research have been carried out in many plants, there is little knowledge about the roles of ALMT in Orchidaceae. In this study, 34 ALMT genes were identified in the genomes of four orchid species. Specifically, ten ALMT genes were found in Dendrobium chrysotoxum and D. catenatum, and seven were found in Apostasia shenzhenica and Phalaenopsis equestris. These ALMT genes were further categorized into four clades (clades 1–4) based on phylogenetic relationships. Sequence alignment and conserved motif analysis revealed that most orchid ALMT proteins contain conserved regions (TM1, GABA binding motif, and WEP motif). We also discovered a unique motif (19) belonging to clade 1, which can serve as a specifically identified characteristic. Comparison with the gene structure of AtALMT genes (Arabidopsis thaliana) showed that the gene structure of ALMT was conserved across species, but the introns were longer in orchids. The promoters of orchid ALMT genes contain many light-responsive and hormone-responsive elements, suggesting that their expression may be regulated by light and phytohormones. Chromosomal localization and collinear analysis of D. chrysotoxum indicated that tandem duplication (TD) is the main reason for the difference in the number of ALMT genes in these orchids. D. catenatum was chosen for the RT-qPCR experiment, and the results showed that the DcaALMT gene expression pattern varied in different tissues. The expression of DcaALMT1-9 was significantly changed after ABA treatment. Combining the circadian CO2 uptake rate, titratable total acid, and RT-qPCR data analysis, most DcaALMT genes were highly expressed at night and around dawn. The result revealed that DcaALMT genes might be involved in photosynthate accumulation. The above study provides more comprehensive information for the ALMT gene family in Orchidaceae and a basis for subsequent functional analysis. [ABSTRACT FROM AUTHOR]

Background Bacterial RNA polymerase (RNAP) is a promising target for antimicrobial chemotherapy due to its indispensable role in bacterial growth and survival. Among its components, only the rpoB gene encoding the β-subunit is known for its association with rifampicin resistance. We recently identified a variant of the RNAP α-subunit (RpoA) in Pseudomonas aeruginosa, conferring heightened bacterial susceptibility to antimicrobials. This susceptibility was attributed to the specific down-regulation of the MexEF-OprN efflux pump. Objectives We asked how to distinguish antimicrobial-susceptible variant strains from clinical isolates. Methods In this study, we identified various P. aeruginosa RpoA variants from clinical sources. Using the sequence alignment of different bacterial RpoA species, we computed the positional conservation of substitutions in RpoA variants using Shannon Entropy. Results Our findings revealed that selective RpoA variant strains exhibited distinct profiles of antimicrobial susceptibility. Notably, RpoA variant strains, containing single-substitutions in the C-terminal domain (α-CTD) but not the N-terminal domain (α-NTD), showed attenuated MexEF-OprN expression and increased susceptibility to MexEF-OprN-specific antibiotics. Furthermore, we observed a close correlation between the susceptibility of these α-CTD RpoA variant strains to antibiotics and the conservation degrees of positional substitutions. Conclusions Our findings demonstrate the prevalence of antimicrobial-susceptible RpoA variant strains among P. aeruginosa clinical isolates. The identified positional conservation pattern in our study facilitates the rapid classification of RpoA variant strains with distinct drug resistances. Given the high conservation of RNAP across bacterial species, our findings open a new therapeutic perspective for precisely and efficiently combating pathogenic RpoA variant strains with specific antimicrobials. [ABSTRACT FROM AUTHOR]

The optimal efficacy of xenogeneically generated proteins intended for application in humans requires that their own antigenicity be minimized. This necessary adaptation of antibodies to a humanized version poses challenges since modifications even distant from the binding sites can greatly influence antigen recognition and this is the primary feature that must be maintained during all modifications. Current strategies often rely on grafting and/or randomization/selection to arrive at a humanized variant retaining the binding properties of the original molecule. However, in terms of speed and efficiency, rationally directed approaches can be superior, provided the requisite structural information is available. We present here a humanization procedure based on the high-resolution X-ray structure of a chimaeric IgG against a marker for multiple myeloma. Based on in silico modelling of humanizing amino acid substitutions identified from sequence alignments, we devised a straightforward cloning procedure to rapidly evaluate the proposed sequence changes. Careful inspection of the structure allowed the identification of a potentially problematic amino acid change that indeed disrupted antigen binding. Subsequent optimization of the antigen binding loop sequences resulted in substantial recovery of binding affinity lost in the completely humanized antibody. X-ray structures of the humanized and optimized variants demonstrate that the antigen binding mode is preserved, with surprisingly few direct contacts to antibody atoms. These results underline the importance of structural information for the efficient optimization of protein therapeutics. Key messages: Structure-based humanization of an IgG against BCMA, a marker for Multiple Myeloma. Identification of problematic mutations and unexpected modification sites. Structures of the modified IgG-antigen complexes verified predictions. Provision of humanized high-affinity IgGs against BCMA for therapeutic applications. [ABSTRACT FROM AUTHOR]

Most third-generation sequencing (TGS) processing tools rely on multiple sequence alignment (MSA) methods to manage sequencing errors. Despite the broad range of MSA approaches available, a limited selection of implementations are commonly used in practice for this type of application, and no comprehensive comparative assessment of existing tools has been undertaken to date. In this context, we have developed an automatic pipeline, named MSA Limit, designed to facilitate the execution and evaluation of diverse MSA methods across a spectrum of conditions representative of TGS reads. MSA Limit offers insights into alignment accuracy, time efficiency, and memory utilization. It serves as a valuable resource for both users and developers, aiding in the assessment of algorithmic performance and assisting users in selecting the most appropriate tool for their specific experimental settings. Through a series of experiments using real and simulated data, we demonstrate the value of such exploration. Our findings reveal that in certain scenarios, popular methods may not consistently exhibit optimal efficiency and that the choice of the most effective method varies depending on factors such as sequencing depth, genome characteristics, and read error patterns. MSA Limit is an open source and freely available tool. All code and data pertaining to it and this manuscript are available at . [ABSTRACT FROM AUTHOR]

Low temperatures are the most significant abiotic stressor for the conservation and production of Phalaenopsis in non-tropical areas. CBF/DREB1 transcription factors play an important role in the plant abiotic stress response. In this study, 31 DREB family members were identified in the Phalaenopsis genome. Expression pattern analysis showed that the expression of different PaDREB members varied among tissue sites. PaDREB1D was isolated from Phalaenopsis aphrodite, and multiple sequence alignment showed that PaDREB1D belonged to the A1 subgroup of the DREB family and was localized in the nucleus. PaDREB1D overexpression in protocorm-like bodies of Phalaenopsis reduced cell damage during low-temperature stress, increased antioxidant enzyme activity, and enhanced the low-temperature tolerance of protocorm-like bodies. The results of this study provide a theoretical basis for breeding for cold resistance and investigating the molecular mechanisms related to low-temperature responses in Phalaenopsis. [ABSTRACT FROM AUTHOR]

The expansin genes are commonly expressed in plant cells, and the encoded proteins influence plant growth and stress resistance by loosening the structure and increasing the flexibility of the cell wall. The objective of this study was to characterize expansin gene promoters in Populus trichocarpa to clarify the regulatory mechanisms underlying gene expression and evolution. Sequence alignments revealed that the similarity among 36 poplar expansin genes was greater for the coding sequences than for the promoter sequences, which suggested these promoter sequences evolved asynchronously. The bases flanking the start codon exhibited a usage bias, with sites +3, +4, and +5 biased toward GC, whereas the other sites were biased toward AT. The flanking sites were significantly correlated with gene expression, especially sites −10 and −17, in which C and G are the bases positively associated with gene expression. A total of 435 regulatory elements (61 types) were identified on the promoters of the poplar expansin genes; Skn-1 was the most common element in 23 promoters. Some expansin genes had more regulatory elements on their promoters (e.g., PtrEXPA4, PtrEXPA3, PtrEXPB3, and PtrEXPB1), whereas some others had less (e.g., PtrEXLA2, PtrEXLB1, and PtrEXPA23). Furthermore, 26 types of elements were involved in expansin gene expression, 25 of which positively affected expression in all analyzed samples. The exception was the endosperm expression-related element Skn-1, which negatively regulated expression in four tissues or treatments. Expression analysis showed that the expansin genes in Populus trichocarpa performed much differently under regular and abiotic stress conditions, which well matched the diversity of their promoter sequences. The results show that expansin genes play an important role in plant growth and development and stress resistance through expression adjustment. [ABSTRACT FROM AUTHOR]

Copyright of Acta Agriculturae Zhejiangensis is the property of Acta Agriculturae Zhejiangensis Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

MOTIVATION: Methods for concept recognition (CR) in clinical texts have largely been tested on abstracts or articles from the medical literature. However, texts from electronic health records (EHRs) frequently contain spelling errors, abbreviations, and other nonstandard ways of representing clinical concepts. RESULTS: Here, we present a method inspired by the BLAST algorithm for biosequence alignment that screens texts for potential matches on the basis of matching k-mer counts and scores candidates based on conformance to typical patterns of spelling errors derived from 2.9 million clinical notes. Our method, the Term-BLAST-like alignment tool (TBLAT) leverages a gold standard corpus for typographical errors to implement a sequence alignment-inspired method for efficient entity linkage. We present a comprehensive experimental comparison of TBLAT with five widely used tools. Experimental results show an increase of 10% in recall on scientific publications and 20% increase in recall on EHR records (when compared against the next best method), hence supporting a significant enhancement of the entity linking task. The method can be used stand-alone or as a complement to existing approaches. AVAILABILITY AND IMPLEMENTATION: Fenominal is a Java library that implements TBLAT for named CR of Human Phenotype Ontology terms and is available at https://github.com/monarch-initiative/fenominal under the GNU General Public License v3.0.