Variety Of Sequences Research Articles

SARS-CoV-2 Omicron Specific Mutations Affecting Infectivity, Fusogenicity, and Partial TMPRSS2-Independency.

The COVID-19 pandemic resulted from the global spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Since its first appearance in 2019, new SARS-CoV-2 variants of concern (VOCs) have emerged frequently, changing the infection's dynamic. SARS-CoV-2 infects cells via two distinct entry routes; receptor-mediated endocytosis or membrane fusion, depending on the absence or presence of transmembrane serine protease 2 (TMPRSS2), respectively. In laboratory conditions, the Omicron SARS-CoV-2 strain inefficiently infects cells predominantly via endocytosis and is phenotypically characterized by decreased syncytia formation compared to the earlier Delta variant. Thus, it is important to characterize Omicron's unique mutations and their phenotypic manifestations. Here, by utilizing SARS-CoV-2 pseudovirions, we report that the specific Omicron Spike F375 residue decreases infectivity, and its conversion to the Delta S375 sequence significantly increases Omicron infectivity. Further, we identified that residue Y655 decreases Omicron's TMPRSS2 dependency and entry via membrane fusion. The Y655H, K764N, K856N and K969N Omicron revertant mutations, bearing the Delta variant sequence, increased the cytopathic effect of cell-cell fusion, suggesting these Omicron-specific residues reduced the severity of SARS-CoV-2. This study of the correlation of the mutational profile with the phenotypic outcome should sensitize our alertness towards emerging VOCs.

The N-terminus of CXCR4 splice variants determines expression and functional properties.

C-X-C motif chemokine ligand 12(CXCL12) is an essential chemokine for organ development and homeostasis in multiple tissues. Its receptor, C-X-C chemokine receptor type 4(CXCR4), is expressed on the surface of target cells. The chemokine and receptor are expressed almost ubiquitously in human tissues and cells throughout life, and abnormal expression of CXCL12 and CXCR4 is observed in pathological conditions, such as inflammation and cancer. CXCR4 is reportedly translated into five splicing variants of different lengths, which each have different amino acids in the N-terminus. As the N-terminus is the first recognition site for chemokines, CXCR4 variants may respond differently to CXCL12. Despite these differences, the molecular and functional properties of CXCR4 variants have not been thoroughly described or compared. Here, we explored the expression of CXCR4 variants in cell lines and analyzed their roles in cellular responses using biochemical approaches. RT-PCR revealed that most cell lines express more than one CXCR4 variant. When expressed in HEK293 cells, the CXCR4 variants differed in protein expression efficiency and cell surface localization. Although variant 2 demonstrated the strongest expression and cell surface localization, variants 1, 3, and 5 also mediated chemokine signaling and induced cellular responses. Our results demonstrate that the N-terminal sequences of each CXCR4 variant determine the expression of the receptor and affect ligand recognition. Functional analyses revealed that CXCR4 variants may also affect each other or interact during CXCL12-stimulated cellular responses. Altogether, our results suggest that CXCR4 variants may have distinct functional roles that warrant additional investigation and could contribute to future development of novel drug interventions.

Enzyme Cascade with Horseradish Peroxidase Readout for High-Throughput Screening and Engineering of Human Arginase-1

We report an enzyme cascade with horseradish peroxidase-based readout for screening human arginase-1 (hArg1) activity. We combined the four enzymes hArg1, ornithine decarboxylase, putrescine oxidase, and horseradish peroxidase in a reaction cascade that generated colorimetric or fluorescent signals in response to hArg1 activity and used this cascade to assay wild-type and variant hArg1 sequences as soluble enzymes and displayed on the surface of Escherichia coli. We screened a curated 13-member hArg1 library covering mutations that modified the electrostatic environment surrounding catalytic residues D128 and H141, and identified the R21E variant with a 13% enhanced catalytic turnover rate compared to wild type. Our scalable one-pot single-step arginase assay with continuous kinetic readout is amenable to high-throughput screening and directed evolution of arginase libraries and testing drug candidates for arginase inhibition.

Design and Rationale the SCAN-MP (Screening for Cardiac Amyloidosis With Nuclear Imaging in Minority Populations) Study.

Background Transthyretin amyloid cardiomyopathy (ATTR-CM) is an important cause of heart failure in older individuals. Misfolding and deposition of transthyretin or prealbumin protein causes ATTR-CM in the context of a normal (wild-type) or variant TTR sequence. Variant ATTR-CM is most commonly caused by the substitution of valine for isoleucine at position 122 in transthyretin (Val122Ile or pV142I, almost exclusively observed in individuals of West African ancestry), demonstrated in 3.4% of self-identified Black individuals in the United States with an estimated 1.5 million carriers. Despite the large number of known pV142I carriers, the proportion of older Black patients with heart failure attributable to ATTR-CM remains unknown. Methods To address this knowledge gap, the SCAN-MP (Screening for Cardiac Amyloidosis with Nuclear Imaging in Minority Populations) study was funded by the National Institutes of Health/National Heart, Lung, and Blood Institute (R01HL139671) to enroll a targeted population of self-identified, community-dwelling Black or Caribbean Hispanic patients (many of whom are of West African ancestry) >60 years of age with heart failure and identify ATTR-CM by noninvasive nuclear imaging. The principal objective of SCAN-MP is to determine the prevalence of ATTR-CM in this population. Secondary objectives will explore TTR genotype, demographics, progression of variant versus wild-type ATTR-CM, and biochemical mechanisms of transthyretin amyloid fibril formation. Conclusions The SCAN-MP study is the largest, prospective study of cardiac amyloidosis in Black and Hispanic individuals. Both wild-type and variant ATTR-CM are now treatable with the US Food and Drug-approved drug tafamidis. The insights gained from SCAN-MP are likely to improve those at risk for or afflicted with ATTR-CM. Registration URL: https://www.clinicaltrials.gov; Unique identifier: NCT03812172.

Date Palm Waste Compost Application Increases Soil Microbial Community Diversity in a Cropping Barley (Hordeum vulgare L.) Field.

Application of date palm waste compost is quite beneficial in improving soil properties and crop growth. However, the effect of its application on soil microbial communities is less understood. High-throughput sequencing and quantitative real-time PCR (qPCR) were used to evaluate the effect of compost application on the soil microbial composition in a barley field during the tillering, booting and ripening stages. The results showed that compost treatment had the highest bacterial and fungal abundance, and its application significantly altered the richness (Chao1 index) and α-diversity (Shannon index) of fungal and bacterial communities. The dominant bacterial phyla found in the samples were Proteobacteria and Actinobacteria while the dominant fungal orders were Ascomycota and Mortierellomycota. Interestingly, compost enriched the relative abundance of beneficial microorganisms such as Chaetomium, Actinobacteriota, Talaromyces and Mortierella and reduced those of harmful microorganisms such as Alternaria, Aspergillus and Neocosmospora. Functional prediction based on Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) showed that amplicon sequence variant (ASV) sequences related to energy metabolism, amino acid metabolism and carbohydrate metabolism were associated with compost-treated soil. Based on Fungi Functional Guild (FUNGuild), identified fungi community metabolic functions such as wood saprotroph, pathotroph, symbiotroph and endophyte were associated with compost-treated soil. Overall, compost addition could be considered as a sustainable practice for establishing a healthy soil microbiome and subsequently improving the soil quality and barley crop production.

A Novel Insertion in the Hepatitis B Virus Surface Protein Leading to Hyperglycosylation Causes Diagnostic and Immune Escape.

Chronic hepatitis B virus (HBV) infection is a global health threat. Mutations in the surface antigen of HBV (HBsAg) may alter its antigenicity, infectivity, and transmissibility. A patient positive for HBV DNA and detectable but low-level HBsAg in parallel with anti-HBs suggested the presence of immune and/or diagnostic escape variants. To support this hypothesis, serum-derived HBs gene sequences were amplified and cloned for sequencing, which revealed infection with exclusively non-wildtype HBV subgenotype (sgt) D3. Three distinct mutations in the antigenic loop of HBsAg that caused additional N-glycosylation were found in the variant sequences, including a previously undescribed six-nucleotide insertion. Cellular and secreted HBsAg was analyzed for N-glycosylation in Western blot after expression in human hepatoma cells. Secreted HBsAg was also subjected to four widely used, state-of-the-art diagnostic assays, which all failed to detect the hyperglycosylated insertion variant. Additionally, the recognition of mutant HBsAg by vaccine- and natural infection-induced anti-HBs antibodies was severely impaired. Taken together, these data suggest that the novel six-nucleotide insertion as well as two other previously described mutations causing hyperglycosylation in combination with immune escape mutations have a critical impact on in vitro diagnostics and likely increase the risk of breakthrough infection by evasion of vaccine-induced immunity.

The telomerase reverse transcriptase elongates reversed replication forks at telomeric repeats.

The telomerase reverse transcriptase elongates telomeres to prevent replicative senescence. This process requires exposure of the 3'-end, which is thought to occur when two sister telomeres are generated at replication completion. Using two-dimensional agarose gel electrophoresis (2D-gels) and electron microscopy, we found that telomeric repeats are hotspots for replication fork reversal. Fork reversal generates 3' telomeric ends before replication completion. To verify whether these ends are elongated by telomerase, we probed de novo telomeric synthesis in situ and at replication intermediates by reconstituting mutant telomerase that adds a variant telomere sequence. We found variant telomeric repeats overlapping with telomeric reversed forks in 2D-gels, but not with normal forks, nontelomeric reversed forks, or telomeric reversed forks with a C-rich 3'-end. Our results define reversed telomeric forks as a substrate of telomerase during replication.

Complete genome sequence of a novel variant of citrus tristeza virus infecting Chinese wild mandarin (Citrus daoxianensis S.W. He & G.F. Liu., syn. Citrus reticulata Blanco) in China

Complete genome sequence of a novel variant of citrus tristeza virus infecting Chinese wild mandarin (Citrus daoxianensis S.W. He & G.F. Liu., syn. Citrus reticulata Blanco) in China

Faecal Short-chain Fatty Acid and Early Introduction of Foods in the First 200 Days of Infant’s Life in the District of Abidjan (Ivory Coast)

Dosage of short chain fatty acids (SCFAs) according to food diet showed that the content of acetate was high in newborn feaces. Infants receiving food supplements have a complex and diverse gut microbiota. Moreover, the results show that infants from poor districts have an abundant concentration of SCFAs in their faeces compared to those living in places with a relatively high standard of living. Among infants receiving milk, the highest proportion of SCFA is acetate in breastfed infants (BF) at a rate of 15.025 ± 2.23 μmol/g, followed by propionate in infants receiving mixed feeding (BF+FF), at a rate of 13.58 ± 1.03 μmol/g and butyrate in infants taking mixed feeding at a rate of 0.32 ± 0.72 μmol/g. However, among infants starting early diet diversification, acetate is higher in infants receiving milk formula and diet diversification (FF+FD) with a concentration of 25.4 ± 0 μmol/g, followed by propionate (2.36 ± 0 μmol/g) in infants receiving mixed feeding (BF+FF) and butyrate in those fed with (BF+FD). Partial breastfeeding is associated with a higher proportion of acetate, butyrate and propionate. The study of the correlation between the different SCFAs produced and the ASV (Variants of Microbial Amplicon Sequences) of the intestinal community of the child, shows that acetate is positively correlated with Bifidobacterium and negatively with Streptococcus and Escherichia-Shigella. Propionate is positively correlated with Bifido bacterium and negatively with Escherichia-Shigella. Similarly, butyrate is positively correlated with Bifidobacterium and negatively with Escherichia-Shigella.

Vaccination using mutated receptor binding domains of SARS-CoV-2: Evidence for partial immune escape but not serotype formation.

SARS-CoV-2 has developed a number of Variants of Concern (VOC) with increased infectivity and/or reduced recognition by neutralizing antibodies specific for the receptor binding domain (RBD) of the spike protein. Extended studies of other viruses have shown that strong and broad viral escape from neutralizing serum antibodies is typically associated with the formation of serotypes. To address the question of serotype formation for SARS-CoV-2 in detail, we generated recombinant RBDs of VOCs and displayed them on virus-like particles (VLPs) for vaccination and specific antibody responses. As expected, mice immunized with wild type (wt) RBD generated antibodies that recognized wt RBD well but displayed reduced binding to VOC RBDs, in particular those with the E484K mutation. Unexpectedly, however, antibodies induced by the VOC vaccines typically recognized best the wt RBDs, often more than the homologous VOC RBDs used for immunization. Hence, these data do not reveal different serotypes but represent a newly observed viral evolution, suggesting a unique situation where inherent differences of RBDs are responsible for induction of neutralizing antibodies. Therefore, besides antibody (fine) specificity, other qualities of antibodies (e.g. their affinity) determine neutralizing capability. Immune escape of SARS-CoV-2 VOCs only affects a fraction of an individual's serum antibodies. Consequently, many neutralizing serum antibodies are cross-reactive and thus protective against multiple current and future VOCs. Besides considering variant sequences for next generation vaccines, broader protection will be achieved with vaccines that induce elevated titers of high-quality antibodies.

Variant-specific deleterious mutations in the SARS-CoV-2 genome reveal immune responses and potentials for prophylactic vaccine development.

Introduction: Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has had a disastrous effect worldwide during the previous three years due to widespread infections with SARS-CoV-2 and its emerging variations. More than 674 million confirmed cases and over 6.7 million deaths have been attributed to successive waves of SARS-CoV-2 infections as of 29th January 2023. Similar to other RNA viruses, SARS-CoV-2 is more susceptible to genetic evolution and spontaneous mutations over time, resulting in the continual emergence of variants with distinct characteristics. Spontaneous mutations of SARS-CoV-2 variants increase its transmissibility, virulence, and disease severity and diminish the efficacy of therapeutics and vaccines, resulting in vaccine-breakthrough infections and re-infection, leading to high mortality and morbidity rates. Materials and methods: In this study, we evaluated 10,531 whole genome sequences of all reported variants globally through a computational approach to assess the spread and emergence of the mutations in the SARS-CoV-2 genome. The available data sources of NextCladeCLI 2.3.0 (https://clades.nextstrain.org/) and NextStrain (https://nextstrain.org/) were searched for tracking SARS-CoV-2 mutations, analysed using the PROVEAN, Polyphen-2, and Predict SNP mutational analysis tools and validated by Machine Learning models. Result: Compared to the Wuhan-Hu-1 reference strain NC 045512.2, genome-wide annotations showed 16,954 mutations in the SARS-CoV-2 genome. We determined that the Omicron variant had 6,307 mutations (retrieved sequence:1947), including 67.8% unique mutations, more than any other variant evaluated in this study. The spike protein of the Omicron variant harboured 876 mutations, including 443 deleterious mutations. Among these deleterious mutations, 187 were common and 256 were unique non-synonymous mutations. In contrast, after analysing 1,884 sequences of the Delta variant, we discovered 4,468 mutations, of which 66% were unique, and not previously reported in other variants. Mutations affecting spike proteins are mostly found in RBD regions for Omicron, whereas most of the Delta variant mutations drawn to focus on amino acid regions ranging from 911 to 924 in the context of epitope prediction (B cell & T cell) and mutational stability impact analysis protruding that Omicron is more transmissible. Discussion: The pathogenesis of the Omicron variant could be prevented if the deleterious and persistent unique immunosuppressive mutations can be targeted for vaccination or small-molecule inhibitor designing. Thus, our findings will help researchers monitor and track the continuously evolving nature of SARS-CoV-2 strains, the associated genetic variants, and their implications for developing effective control and prophylaxis strategies.

Molecular Diagnosis of Duchenne Muscular Dystrophy Using Single NGS-Based Assay.

Duchenne Muscular Dystrophy (DMD) is an X-linked inherited neuromuscular disorder caused by pathogenic variants in the dystrophin gene (DMD; locus Xp21.2). The variant spectrum of DMD is unique in that 65% of causative mutations are intragenic deletions, with intragenic duplications and point mutations (along with other sequence variants) accounting for 6% to 10% and 30% to 35%, respectively. The traditional strategy for molecular diagnostic testing for DMD involves initial screening for deletions/duplications using microarray-based comparative genomic hybridization followed by a full-sequence analysis of DMD for sequence variants. This traditional strategy is expensive and time-consuming due to the involvement of two separate tests to detect all types of variants in the DMD gene. Recent advancements in next-generation sequencing (NGS) technology and improvements in analysis algorithms related to copy number variant detection ultimately resulted in the development of a single NGS-based assay to detect all variant types, including deletions/duplications and sequence variants. This article initially discusses the strategic algorithm for establishing a molecular diagnosis of DMD and later provides detailed molecular diagnostic protocols for DMD, including an NGS-based sequencing assay with sequence and copy number variant analysis. © 2023 Wiley Periodicals LLC. Basic Protocol: Next-generation sequencing of the entire genomic sequence of the DMD gene using IDT xGen Lockdown Probes.

Detection of a SARS-CoV-2 Beta-like variant in coastal Kenya after more than a year of disappearance

Continuous genomic surveillance is necessary and important to inform emergence of new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants and resurgence of previous circulating variants. In this study, we provide genomic characterisation of a Beta variant sequence identified through out genomic surveillance platform more than a year since the last reported case of Beta. This variant contained additional mutations associated with immune escape that have been observed in other newer variants such as Delta and Omicron implying ongoing convergent evolution of this variant in the community or in an immunocompromised patient.

High-throughput telomere length measurement at nucleotide resolution using the PacBio high fidelity sequencing platform

Telomeres are specialized nucleoprotein structures at the ends of linear chromosomes. The progressive shortening of steady-state telomere length in normal human somatic cells is a promising biomarker for age-associated diseases. However, there remain substantial challenges in quantifying telomere length due to the lack of high-throughput method with nucleotide resolution for individual telomere. Here, we describe a workflow to capture telomeres using newly designed telobaits in human culture cell lines as well as clinical patient samples and measure their length accurately at nucleotide resolution using single-molecule real-time (SMRT) sequencing. Our results also reveal the extreme heterogeneity of telomeric variant sequences (TVSs) that are dispersed throughout the telomere repeat region. The presence of TVSs disrupts the continuity of the canonical (5’-TTAGGG-3’)n telomere repeats, which affects the binding of shelterin complexes at the chromosomal ends and telomere protection. These findings may have profound implications in human aging and diseases.

Exploring broilers and native fowls of Andaman and Nicobar Islands as a source of β-lactamase-producing Enterobacteriaceae even with limited anthropogenic activities and docking-based identification of catalytic domains in novel β-lactamase variants.

The present study was conducted to detect the occurrence of β-lactamase and biofilm-producing Escherichia coli, Salmonella, and Klebsiella in broilers and native fowl reared in the Andaman and Nicobar Islands, India. The study also included molecular docking experiments to confirm the nature of the catalytic domains found in the β-lactamase variants obtained and to reveal the clonal relationship of the isolates with human clinical strains from the database. A total of 199 cloacal swabs were collected from five poultry breeds/varieties (broiler, Vanraja, Desi, Nicobari, and layer) in three districts of the Andaman and Nicobar Islands. E. coli, Salmonella enterica, and Klebsiella pneumoniae were isolated by standard techniques and confirmed by PCR. Phenotypical β-lactamase producers were identified by a double-disc test. The genes (bla CTX, bla SHV, bla TEM , and bla AmpC) were screened, and selected sequences of β-lactamase variants were submitted to DDBJ. Homology modeling, model validation, and active site identification of different β-lactamase variants were done by the SWISS-MODEL. Molecular docking was performed to identify the catalytic domains of the β-lactamase variants. The selected β-lactamase sequences were compared with the Indian ESBL sequences from human clinical strains in NCBI-GenBank. In total, 425 Enterobacteriaceae strains were isolated from the collected samples. Klebsiella pneumoniae (42.58%) was found to be the most prevalent, followed by Salmonella enterica (30.82%) and E. coli (26.58%). The phenotypical antibiogram of all 425 isolates showed the highest resistance against oxytetracycline (61-76%) and the lowest against gentamicin (15-20%). Phenotypical production of β-lactamase enzymes was observed in 141 (33.38%) isolates. The isolation rate of β-lactamase producing E. coli, Salmonella enterica, and Klebsiella pneumoniae was significantly higher (p < 0.05) in the birds reared in the South Andaman district (25.6, 17.5, and 18.7%, respectively) than in Nicobar (11.5, 7.6, 7.1%, respectively). Genotyping of the β-lactamase-producing isolates revealed the maximum possession of bla TEM, followed by bla SHV and bla CTX - M. The nucleotide sequences were found to be similar with bla CTX - M-15, bla SHV - 11, bla SHV - 27, bla SHV - 228, bla TEM - 1, and bla AmpC in BLAST search. Distribution of studied biofilm-associated genes in Enterobacteriaceae strains from different varieties of the birds revealed that the layer birds had the maximum possession, followed by Vanraja, Desi, broilers, and Nicobari fowls. The phylogenetic analysis of selected sequences revealed a partial clonal relationship with human clinical strains of the Indian subcontinent. Molecular docking depicted the Gibbs free energy release for 10 different macromolecules (proteins) and ligand (antibiotic) complexes, ranging from -8.1 (SHV-27 + cefotaxime) to -7 (TEM-1 + cefotaxime) kcal/mol. The study revealed β-lactamase variants circulating in the fowl population of the Andaman and Nicobar Islands (India), even in remote places with low anthropogenic activity. Most of the strains possessed bla TEM - 1, followed by bla CTX - M-15. Possession of bla SHV - 11, bla SHV - 27, and bla SHV - 228 in poultry Enterobacteriaceae strains was not reported earlier from any part of the world. The phylogenetic analysis revealed a partial clonal relationship of β-lactamase sequences with the human clinical strains isolated from the Indian subcontinent.

Factory Data Model for Planning Variant Sequences, Expansion Stages and Reconfigurations

A factory has to adapt to a future production program and increasing volumes. The planning of new factories must therefore already take expansion stages and reconfigurations into account. This leads to an increase in the complexity of planning, as different variants have to be planned and put into a chronologically reasonable and cost-optimized sequence. A factory data model that covers expansion stages manages the planning data amount and reduces complexity. Based on existing models from the literature, a factory data model was developed and modeled in the Unified Modeling Language standard. The factory data model shows which classes are interdependent and how planning results are stored reusable for other variants. An example is used to explain the developed factory data model and highlighting its advantages. The aim is to enable factory planners to document variant sequences and expansion stages in a comprehensible way, which supports the interdisciplinarity of factory planning projects.

VPrimer: A Method of Designing and Updating Primer and Probe With High Variant Coverage for RNA Virus Detection.

Fatal infectious diseases caused by RNA viruses, such as COVID-19, have emerged around the world. RT-PCR is widely employed for virus detection, and its accuracy depends on the primers and probes since RT-PCR can detect a virus only when the primers and probes bind to the target gene of the virus. Most of primer design methods are for a single host and so require a great deal of effort to design for RNA virus detection, including homology tests among the host and all the viruses for the host using BLAST-like tools. Furthermore, they do not consider variant sequences, which are very common in viruses. In this study, we describe VPrimer, a method of designing high-quality primer-probe sets for RNA viruses. VPrimer can find primer-probe sets that cover more than 95% of the variants of a target virus but do not cover any sequences of other viruses or the host. With VPrimer, we found 381,698,582 primer-probe sets for 3,104 RNA viruses. Multiplex PCR assays using the top 2 primer-probe sets suggested by VPrimer usually cover 100% of variants. To address the rapid changes in viral genomes, VPrimer finds the best and up-to-date primer-probe sets incrementally against the most recently reported variants.

Accurate prediction of functional effect of single amino acid variants with deep learning

The assessment of functional effect of amino acid variants is a critical biological problem in proteomics for clinical medicine and protein engineering. Although natively occurring variants offer insights into deleterious variants, high-throughput deep mutational experiments enable comprehensive investigation of amino acid variants for a given protein. However, these mutational experiments are too expensive to dissect millions of variants on thousands of proteins. Thus, computational approaches have been proposed, but they heavily rely on hand-crafted evolutionary conservation, limiting their accuracy. Recent advancement in transformers provides a promising solution to precisely estimate the functional effects of protein variants on high-throughput experimental data. Here, we introduce a novel deep learning model, namely Rep2Mut-V2, which leverages learned representation from transformer models. Rep2Mut-V2 significantly enhances the prediction accuracy for 27 types of measurements of functional effects of protein variants. In the evaluation of 38 protein datasets with 118,933 single amino acid variants, Rep2Mut-V2 achieved an average Spearman’s correlation coefficient of 0.7. This surpasses the performance of six state-of-the-art methods, including the recently released methods ESM, DeepSequence and EVE. Even with limited training data, Rep2Mut-V2 outperforms ESM and DeepSequence, showing its potential to extend high-throughput experimental analysis for more protein variants to reduce experimental cost. In conclusion, Rep2Mut-V2 provides accurate predictions of the functional effects of single amino acid variants of protein coding sequences. This tool can significantly aid in the interpretation of variants in human disease studies.

Difficulties Of Enterobacteriaceae Genome Annotation In Deciphering Gastrointestinal Microbiome Datasets Obtained By 16S RRNA Gene Amplicon Sequencing

Sequencing of the 16S rRNA gene amplicon is the cornerstone of the method for studying diverse bacteria in complex microbial communities. However, its use is complicated by an error rate of 10–17% when annotating 16S rRNA gene sequences. In our study, we examined the degree of accuracy of the taxonomic database of Enterobacteriaceae, compiled using the SILVA 132 reference database and a previously obtained dataset, viz. the microbiome of the gastrointestinal tract in adolescents with normal body weight and obesity. Material and Methods — In this study, previously obtained 16S rRNA gene amplicon sequencing data were used, and the deciphering was carried out using the QIIME2 2019.4 platform. Phylogenetic analysis was performed using MEGA X software. Results — Phylogenetic analysis of this family based on the studied V3–V4 fragment was hampered by polyphyly among some genera, and for half of the variants of the amplicon sequences it was not possible to clarify their genus. Statistical analysis did not reveal significant differences between the samples. Conclusion — Although the average values of bacterial genera in the studied groups intuitively differed from each other, statistical analysis did not reveal significant differences between the samples. However, it can be assumed that a more detailed study of taxonomic diversity, taking into account factors, such as enterotype, duration of breastfeeding and family history, may reveal differences in the frequency distribution.

Unsupervised outlier detection applied to SARS-CoV-2 nucleotide sequences can identify sequences of common variants and other variants of interest

As of June 2022, the GISAID database contains more than 11 million SARS-CoV-2 genomes, including several thousand nucleotide sequences for the most common variants such as delta or omicron. These SARS-CoV-2 strains have been collected from patients around the world since the beginning of the pandemic. We start by assessing the similarity of all pairs of nucleotide sequences using the Jaccard index and principal component analysis. As shown previously in the literature, an unsupervised cluster analysis applied to the SARS-CoV-2 genomes results in clusters of sequences according to certain characteristics such as their strain or their clade. Importantly, we observe that nucleotide sequences of common variants are often outliers in clusters of sequences stemming from variants identified earlier on during the pandemic. Motivated by this finding, we are interested in applying outlier detection to nucleotide sequences. We demonstrate that nucleotide sequences of common variants (such as alpha, delta, or omicron) can be identified solely based on a statistical outlier criterion. We argue that outlier detection might be a useful surveillance tool to identify emerging variants in real time as the pandemic progresses.