First reporting of BA.1* and BA.2* recombinant SARS-CoV-2 lineage XAP from Indian wastewaters

Tracking new variants of SARS-CoV-2 is vital for managing COVID-19 spread and allocating resources. Domestic antigen testing has created surveillance gaps that make it hard to identify new viral variants. We conducted whole genome sequencing of wastewater viral genes from major and minor treatment facilities in Dehradun from March 2022 onwards. Based on our analysis, the samples that achieved higher sequencing depth and covered >90 % of the viral genome uncovered a major variant pattern resembling the XAP recombinant lineage that is reported for the first time in the City of Dehradun, Uttrakhand and is the first ever records in India as on date. This novel XAP recombinant lineage had 9, 2, 30, 1, 2, 5, 1, 1, 1 aminoacid changes (total 54 mutations) in Orf1a, Orf1b, S, E, M, N, Orf3a, Orf6 and Orf8 regions of the gene respectively that shares 49 mutations common to the ancestral lineages BA.1* and BA.2*, with 6 unique mutations. Subsequent comparison and analysis of the clinical sequence data from the region post-detection of this rare and unusual variant showed no causalities infected with the newly detected XAP lineage. These findings are indicative of future alarming situation with plausible threats of fresh spur of Omicron variant led infections in the urban community.

The evolution and diversity of porcine reproductive and respiratory syndrome virus in China

The porcine reproductive and respiratory syndrome virus (PRRSV) has emerged as a significant pathogen in the global pork industry since the late 1980s, causing substantial economic losses due to its high contagiousness and genetic variability. China, with its complex epidemiological landscape, has witnessed the emergence of four distinct lineages of PRRSV-2 (Lineages 1, 3, 5, and 8) and occasional occurrences of PRRSV-1. This review summarizes the historical context and epidemiological trends that have led to the diversification of PRRSV in China, discusses the evolutionary dynamics behind the establishment of diverse genetic variants, as well as the impact of recombination and modified live vaccines (MLVs) on the virus's rapid evolution. The implications for disease management, including strategies to reduce the complexity of PRRSV epidemics and improve prevention and control measures, are also suggested. Understanding the evolutionary pattern and factors contributing to PRRSV diversity is crucial for enhancing our knowledge, control capabilities, and prevention strategies, which could be integrated into swine health management practices.

Drug resistance and genomic variations among Mycobacterium tuberculosis isolates from The Nile Delta, Egypt

Tuberculosis is a global public health concern. Earlier reports suggested the emergence of high rates of drug resistant tuberculosis in Egypt. This study included 102 isolates of Mycobacterium tuberculosis collected from two reference laboratories in Cairo and Alexandria. All clinical isolates were sub-cultured on Löwenstein-Jensen medium and analyzed using both BD BACTEC MGIT 960 SIRE Kit and standard diffusion disk assays to identify the antibiotic sensitivity profile. Extracted genomic DNA was subjected to whole genome sequencing (WGS) using Illumina platform. Isolates that belong to lineage 4 represented > 80%, while lineage 3 represented only 11% of the isolates. The percentage of drug resistance for the streptomycin, isoniazid, rifampicin and ethambutol were 31.0, 17.2, 19.5 and 20.7, respectively. Nearly 47.1% of the isolates were sensitive to the four anti-tuberculous drugs, while only one isolate was resistant to all four drugs. In addition, several new and known mutations were identified by WGS. High rates of drug resistance and new mutations were identified in our isolates. Tuberculosis control measures should focus on the spread of mono (S, I, R, E)- and double (S, E)-drug resistant strains present at higher rates throughout the whole Nile Delta, Egypt.

Genomic surveillance of SARS-CoV-2 in North Africa: 4 years of GISAID data sharing

Objectives

This study aimed to construct geographically, temporally, and epidemiologically representative data sets for SARS-CoV-2 in North Africa, focusing on Variants of Concern (VOCs), Variants of Interest (VOIs), and Variants Under Monitoring (VUMs).

Methods

SARS-CoV-2 genomic sequences and metadata from the EpiCoV database via the Global Initiative on Sharing All Influenza Data platform were analyzed. Data analysis included cases, deaths, demographics, patient status, sequencing technologies, and variant analysis.

Results

A comprehensive analysis of 10,783 viral genomic sequences from six North African countries revealed notable insights. SARS-CoV-2 sampling methods lack standardization, with a majority of countries lacking clear strategies. Over 59% of analyzed genomes lack essential clinical and demographic metadata, including patient age, sex, underlying health conditions, and clinical outcomes, which are essential for comprehensive genomic analysis and epidemiological studies, as submitted to the Global Initiative on Sharing All Influenza Data. Morocco reported the highest number of confirmed COVID-19 cases (1,272,490), whereas Tunisia leads in reported deaths (29,341), emphasizing regional variations in the pandemic's impact. The GRA clade emerged as predominant in North African countries. The lineage analysis showcased a diversity of 190 lineages in Egypt, 26 in Libya, 121 in Tunisia, 90 in Algeria, 146 in Morocco, and 10 in Mauritania. The temporal dynamics of SARS-CoV-2 variants revealed distinct waves driven by different variants.

Conclusions

This study contributes valuable insights into the genomic landscape of SARS-CoV-2 in North Africa, highlighting the importance of genomic surveillance in understanding viral dynamics and informing public health strategies.

Genetic diversity of Mycobacterium tuberculosis strains isolated from spiritual holy water site attendees in Northwest Ethiopia. A cross-sectional study

Background

The genetic diversity of Mycobacterium tuberculosis complex (MTBC) strains was characterized among isolates from individuals with pulmonary tuberculosis (PTB) symptoms attended holy water sites (HWSs) in the Amhara region, Ethiopia.

Methods

A cross-sectional study was done from June 2019 to March 2020 to describe the genetic diversity and drug-resistance profiles of MTBC isolates. Sputum specimens were collected and cultured in the Löwenstein-Jensen culture medium. Line Probe Assay, MTBDRplus VER 2.0, and MTBDRsl VER 2.0 were used to detect first-and second-line anti-TB drug-resistance patterns. A spoligotyping technique was utilized to characterize the genetic diversity. Statistical analysis was performed using STATA 15.

Results

Of 560 PTB-symptomatic participants, 122 (21.8%) were culture-positive cases. Spoligotyping of 116 isolates revealed diverse MTBC sublineages, with four major lineages: Euro-American (EA) (Lineage 4), East-African-Indian (EAI) (Lineage 3), Ethiopian (ETH) (Lineage 7), East Asian (EA) (Lineage 2). The majority (96.6%) of the isolates were EA (lineage 4) and EAI, with proportions of 54.3% and 42.2%, respectively. A total of 31 spoligotype patterns were identified, 26 of which were documented in the SITVIT2 database. Of these, there were 15 unique spoligotypes, while eleven were grouped with 2-17 isolates. SIT149/T3-ETH (n = 17), SIT26/CAS1-DELHI (n = 16), SIT25/CAS1-DELHI (n = 12), and SIT52/T2 (n = 11) spoligotypes were predominant. A rare spoligotype pattern: SIT41/Turkey and SIT1/Beijing, has also been identified in North Shewa. The overall clustering rate of sub-lineages with known SIT was 76.4%.Of the 122 culture-positive isolates tested, 16.4% were resistant to rifampicin (RIF) and/or isoniazid (INH). Multidrug-resistant TB (MDR-TB) was detected in 12.3% of isolates, five of which were fluoroquinolones (FLQs) resistant. SIT149/T3-ETH and SIT21/CAS1-KILI sublineages showed a higher proportion of drug resistance.

Conclusions

Diverse MTBC spoligotypes were identified, with the T and CAS families and EA (lineage 4) predominating. A high prevalence of drug-resistant TB, with SIT149/T3-ETH and CAS1-KILI sublineages comprising a greater share, was observed. A study with large sample size and a sequencing method with stronger discriminatory power is warranted to understand better the genetic diversity of circulating MTBC in this cohort of study, which would help to adopt targeted interventions.

Occurrence of Anaplasma spp. in wild lagomorphs from Southern Spain: Molecular detection of new Anaplasma bovis lineages

Wild lagomorphs can act as reservoirs of several pathogens of public and animal health concern. However, the number of studies assessing the presence of Anaplasma spp. in these species is scarce. The aim of the present study was to molecularly identify Anaplasma spp. in wild rabbits (Oryctolagus cuniculus) and Iberian hares (Lepus granatensis) from Southern Spain and assess their epidemiological role in the maintenance of the bacterium. During 2017-2021, spleen samples of 394 wild rabbits and 145 Iberian hares were collected. Anaplasma DNA was detected using different PCR assays (16S rRNA and groEL) and phylogenetic analyses were carried out by Bayesian approach. The possible influence of lagomorph species, age and sex on the prevalence of Anaplasma spp. was evaluated by a multiple logistic regression model. The 9.4% of the rabbits were positive to Anaplasma bovis, but all the hares were negative. No significant differences were found in Anaplasma spp. prevalence regarding to age or sex. This is the first report of A. bovis in lagomorphs from Europe. The phylogenetic analysis of A. bovis confirms the existence of different clusters suggesting the existence of several lineages. In addition, a high divergence of nucleotide identity was observed within the lineage 4, which could result in the under-detection of some strains when using A. bovis-specific PCR, hindering its detection and characterization. Since this analysis is based on a limited number of nucleotide bases and sequences, more studies are needed for further characterize A. bovis, as well as its relationship with other Anaplasma spp.

Molecular insights into the SARS-CoV-2 Omicron variant from Bangladesh suggest diverse and continuous evolution

The study analyzed the molecular dynamics of the circulating SARS-CoV-2 Omicron variant from its identification in November 2021 to January 2023. The SARS-CoV-2 sequences from Bangladesh revealed three distinct waves of the Omicron variant. More than 50 sub-lineages of Omicron variant were introduced into the country, with the majority belonging to the major lineages of BA.1-like (24.91%), BA.2-like (43.35%), BA.5-like (5.76%), XBB (10.47%), and "Others and Unassigned" (18.64%). Furthermore, the relative frequencies over time revealed that Omicron lineages existed for a short period of time before being replaced by other sub-lineages. Many potential mutations were found in the receptor binding domain of the Spike protein including G339D/H, S371 L/F, K417 N, T478K, E484A, Q493R, Q498R, and N501Y. In conclusion, the SARS-CoV-2 Omicron variant from Bangladesh showed diverse genetic features and continuous evolution. Therefore, the choice of vaccine and monitoring of hospitalized patients is important alongside genetic characterization of the circulating SARS-CoV-2.

Molecular investigations of Mycobacterium tuberculosis genotypes among baseline and follow-up strains circulating in four regions of Eswatini

BACKGROUND

The tuberculosis (TB) epidemic remains a major global health problem and Eswatini is not excluded. Our study investigated the circulating genotypes in Eswatini and compared them at baseline (start of treatment) and follow-up during TB treatment.

METHODS

Three hundred and ninety (n = 390) participants were prospectively enrolled from referral clinics and patients who met the inclusion criteria, were included in the study. A total of 103 participants provided specimens at baseline and follow-up within six months. Mycobacterium tuberculosis (M.tb) strains were detected by GeneXpert® MTB/RIF assay (Cephied, USA) and Ziehl -Neelsen (ZN) microscopy respectively at baseline and follow-up time-points respectively. The 206 collected specimens were decontaminated and cultured on BACTEC™ MGIT™ 960 Mycobacteria Culture System (Becton Dickinson, USA). Drug sensitivity testing was performed at both baseline and follow-up time points. Spoligotyping was performed on both baseline and follow-up strains after DNA extraction.

RESULTS

Resistance to at least one first line drug was detected higher at baseline compared to follow-up specimens with most of them developing into multidrug-resistant (MDR)-TB. A total of four lineages and twenty genotypes were detected. The distribution of the lineages varied among the different regions in Eswatini. The Euro-American lineage was the most prevalent with 46.12% (95/206) followed by the East Asian with 24.27% (50/206); Indo-Oceanic at 9.71% (20/206) and Central Asian at 1.94% (4/206). Furthermore, a high proportion of the Beijing genotype at 24.27% (50/206) and S genotype at 16.50% (34/206) were detected. The Beijing genotype was predominant in follow-up specimens collected from the Manzini region with 48.9% (23/47) (p = 0.001). A significant proportion of follow-up specimens developed MDR-TB (p = 0.001) with Beijing being the major genotype in most follow-up specimens (p < 0.000).

CONCLUSION

Eswatini has a high M.tb genotypic diversity. A significant proportion of the TB infected participants had the Beijing genotype associated with MDR-TB in follow-up specimens and thus indicate community wide transmission.

Analysis of self-renewing and differentiation-related microRNAs and transcription factors in multilineage mouse hematopoietic stem/progenitor cells induced by 1,4-benzoquinone

Background

HSPCs are targets for benzene-induced hematotoxicity and leukemogenesis. However, benzene toxicity targeting microRNAs (miRNAs) and transcription factors (TF) that are involve in regulating self-renewing and differentiation of HSPCs comprising of different hematopoietic lineages remains poorly understood. In this study, the effect of a benzene metabolite, 1,4-benzoquinone (1,4-BQ) exposure, in HSPCs focusing on the self-renewing (miRNAs: miR-196b and miR-29a; TF: HoxB4, Bmi-1) and differentiation (miRNAs: miR-181a, TF: GATA3) pathways were investigated.

Methods

Freshly isolated mouse BM cells were initially exposed to 1,4-BQ at 1.25 to 5 µM for 24 h, followed by miRNAs and TF studies in BM cells. Then, the miRNAs expression was further evaluated in HSPCs of different lineages comprised of myeloid, erythroid and pre-B lymphoid progenitors following 7-14 days of colony forming unit (CFU) assay.

Results

Exposure to 1,4-BQ in BM cells significantly (p < 0.05) reduced the miR-196b (2.5 and 5 µM), miR-181a (1.25, 2.5 and 5 µM) and miR-29a (1.25 µM) along with upregulation of miR-29a at 2.5 µM. Meanwhile, 1,4-BQ exposure in HSPCs significantly increased the miR-196b expression level (p < 0.05) only in myeloid and pre-B lymphoid progenitors at 2.5 and 5 µM. Significant (p < 0.05) reduction in expression of miR-181a in myeloid (1.25 µM), erythroid (5 µM) progenitors along with miR-29a in myeloid (1.25 µM) and pre-B lymphoid (5 µM) progenitors were noted following exposure to 1,4-BQ. Meanwhile, increased expression of miR-181a was observed in pre-B lymphoid progenitor upon exposure to 1,4-BQ, but only at 5 µM. As for TF studies, expression of HoxB4 protein was significantly increased (p < 0.05) at all 1,4-BQ concentrations as compared to Bmi-1 and GATA3, which were significantly (p < 0.05) elevated starting at 2.5 µM of 1,4-BQ.

Conclusion

1,4-BQ induces aberration of miRNAs and transcription factors protein expression that are involved in regulating self-renewing and differentiation pathways of HSPCs. Moreover, epigenetic toxicity as evidenced from the miRNAs expression was found to be mediated by a lineage-driven mechanism. The role of cell lineage in governing the toxicity of 1,4-BQ in HSPCs lineages deserves further investigation.

Monitoring SARS-CoV-2 variants in wastewater of Dhaka City, Bangladesh: approach to complement public health surveillance systems

BACKGROUND

Wastewater-based epidemiological surveillance has been considered a powerful tool for early detection and monitoring of the dynamics of SARS-CoV-2 and its lineages circulating in a community. This study is aimed to investigate the complexity of SARS-CoV-2 infection dynamics in Dhaka city by examining its genetic variants in wastewater. Also, the study seeks to determine a connection between the SARS-CoV-2 variations detected in clinical testing and those found in wastewater samples.

RESULTS

Out of 504 samples tested in RT-qPCR, 185 (36.7%) tested positive for SARS-CoV-2 viral RNA. The median log10 concentration of SARS-CoV-2 N gene copies/Liter of wastewater (gc/L) was 5.2, and the median log10 concentration of ORF1ab was 4.9. To further reveal the genetic diversity of SARS-CoV-2, ten samples with ORF1ab real-time RT-PCR cycle threshold (Ct) values ranging from 28.78 to 32.13 were subjected to whole genome sequencing using nanopore technology. According to clade classification, sequences from wastewater samples were grouped into 4 clades: 20A, 20B, 21A, 21J, and the Pango lineage, B.1, B.1.1, B.1.1.25, and B.1.617.2, with coverage ranging from 94.2 to 99.8%. Of them, 70% belonged to clade 20B, followed by 10% to clade 20A, 21A, and 21J. Lineage B.1.1.25 was predominant in Bangladesh and phylogenetically related to the sequences from India, the USA, Canada, the UK, and Italy. The Delta variant (B.1.617.2) was first identified in clinical samples at the beginning of May 2021. In contrast, we found that it was circulating in the community and was detected in wastewater in September 2020.

CONCLUSION

Environmental surveillance is useful for monitoring temporal and spatial trends of existing and emerging infectious diseases and supports evidence-based public health measures. The findings of this study supported the use of wastewater-based epidemiology and provided the baseline data for the dynamics of SARS-CoV-2 variants in the wastewater environment in Dhaka, Bangladesh.

Lineages to circuits: the developmental and evolutionary architecture of information channels into the central complex

The representation and integration of internal and external cues is crucial for any organism to execute appropriate behaviors. In insects, a highly conserved region of the brain, the central complex (CX), functions in the representation of spatial information and behavioral states, as well as the transformation of this information into desired navigational commands. How does this relatively invariant structure enable the incorporation of information from the diversity of anatomical, behavioral, and ecological niches occupied by insects? Here, we examine the input channels to the CX in the context of their development and evolution. Insect brains develop from ~ 100 neuroblasts per hemisphere that divide systematically to form "lineages" of sister neurons, that project to their target neuropils along anatomically characteristic tracts. Overlaying this developmental tract information onto the recently generated Drosophila "hemibrain" connectome and integrating this information with the anatomical and physiological recording of neurons in other species, we observe neuropil and lineage-specific innervation, connectivity, and activity profiles in CX input channels. We posit that the proliferative potential of neuroblasts and the lineage-based architecture of information channels enable the modification of neural networks across existing, novel, and deprecated modalities in a species-specific manner, thus forming the substrate for the evolution and diversification of insect navigational circuits.

Vaccine breakthrough and rebound infections modeling: Analysis for the United States and the ten U.S. HHS regions

A vaccine breakthrough infection and a rebound infection cases of COVID-19 are studied and analyzed for the ten U.S. Department of Health and Human Services (HHS) regions and the United States as a nation in this work. An innovative multi-strain susceptible-vaccinated-exposed-asymptomatic-symptomatic-recovered (SVEAIR) epidemic model is developed for this purpose for a population assumed to be susceptible to n-different variants of the disease, and those who are vaccinated and recovered from a specific strain k(k ≤ n) of the disease are immune to present strain and its predecessors j = 1, 2, …, k, but can still be infected by newer emerging strains j = k + 1, k + 2, …, n. The model is used to estimate epidemiological parameters, namely, the latent and infectious periods, the transmission rates, vaccination rates, recovery rates for each of the Delta B.1.617.2, Omicron B.1.1.529, and lineages BA.2, BA.2.12.1, BA.4, BA.5, BA.1.1, BA.4.6, and BA.5.2.6 for the United States and for each of the ten HHS regions. The transmission rate is estimated for both the asymptomatic and symptomatic cases. The effect of vaccines on each strain is analyzed. Condition that guarantees existence of an endemic with certain number of strains is derived and used to describe the endemic state of the population.

Identification of a new polymorphism on the wild-type canine distemper virus genome: could this contribute to vaccine failures?

The canine distemper virus (CDV) is responsible for a multisystem infectious disease with high prevalence in dogs and wild carnivores and has vaccination as the main control measure. However, recent studies show an increase in cases including vaccinated dogs in different parts of the world. There are several reasons for vaccine failures, including differences between vaccine strains and wild-type strains. In this study, a phylogenetic analysis of CDV strains from samples of naturally infected, vaccinated, and symptomatic dogs in Goiânia, Goiás, Brazil was performed with partial sequencing of the hemagglutinin (H) gene of CDV. Different sites of amino acid substitutions were found, and one strain had the Y549H mutation, typically present in samples from wild animals. Substitutions in epitopes (residues 367, 376, 379, 381, 386, and 388) that may interfere with the vaccine's ability to provide adequate protection against infection for CDV were observed. The identified strains were grouped in the South America 1/Europe lineage, with a significant difference from other lineages and vaccine strains. Twelve subgenotypes were characterized, considering a nucleotide identity of at least 98% among the strains. These findings highlight the relevance of canine distemper infection and support the need better monitoring of the circulating strains that contribute to elucidate if there is a need for vaccine update.

Exploration of the diversity of multi-drug resistant Mycobacterium tuberculosis complex in Lagos, Nigeria using WGS: Distribution of lineages, drug resistance patterns and genetic mutations

Multidrug-resistant (MDR) or extensively drug-resistant (XDR) Tuberculosis (TB) is a major challenge to global TB control. Therefore, accurate tracing of in-country MDR-TB transmission are crucial for the development of optimal TB management strategies. This study aimed to investigate the diversity of MTBC in Nigeria. The lineage and drug-resistance patterns of the clinical MTBC isolates of TB patients in Southwestern region of Nigeria were determined using the WGS approach. The phenotypic DST of the isolates was determined for nine anti-TB drugs. The sequencing achieved average genome coverage of 65.99X. The most represented lineages were L4 (n = 52, 83%), L1 (n = 8, 12%), L2 (n = 2, 3%) and L5 (n = 1, 2%), suggesting a diversified MTB population. In term of detection of M/XDR-TB, while mutations in katG and rpoB genes are the strong predictors for the presence of M/XDR-TB, the current study also found the lack of good genetic markers for drug resistance amongst the MTBC in Nigeria which may pose greater problems on local tuberculosis management efforts. This high-resolution molecular epidemiological data provides valuable insights into the mechanistic for M/XDR TB in Lagos, Nigeria.

An appraisal of gene targets for phylogenetic classification of canine distemper virus: Is the hemagglutinin the best candidate?

Sequence analysis of the canine distemper virus (CDV) hemagglutinin (H) gene may provide important insights on virus-host interactions and has also been frequently used for CDV phylogenetic classification. Herein, we performed an in silico analysis of CDV complete genomes (CGs) available in GenBank in order to investigate the suitability of H for CDV classification into lineages/genotypes. In addition, we analyzed the other viral genes for their potential use in CDV classification. Initially, we collected 116 CDV CGs from GenBank and compared their phylogenetic classification with that of their respective H nucleotide (nt) and amino acid (aa) sequences. Subsequently, we calculated the geodesic distance between the CG and H phylogenetic trees. These analyses were later performed with other CDV genes. All CDV CGs were also evaluated for possible recombination events. Nucleotide and aa analyses of H misclassified some Vaccine/America 1/Asia 3 lineage sequences compared to CG analysis, finding supported by both Maximum Likelihood (ML) and Bayesian Markov Chain Monte Carlo (B-MCMC) methods. Moreover, aa-based H analysis showed additional disagreements with the classification obtained by CG. The geodesic distance between the H and CG trees was 0.0680. Strong recombination signals were identified in the H gene, including Vaccine/America 1/Asia 3 lineage sequences. In contrast, C and P were the only genes that fully reproduced the CG classification (by ML and/or B-MCMC) and that did not show strong recombination signals. Furthermore, the P phylogenetic tree showed the lowest geodesic distance from the CG tree (0.0369). These findings suggest C and P as potential targets for CDV phylogenetic classification, especially when full genome sequencing is not possible. Finally, since our results were obtained considering the CDV CGs available to date, future analyses performed as more CDV sequences become available will be useful to assess probable issues of H-based phylogeny and to consolidate the suitability of the C and P genes for CDV classification.

Comparative genomics of drug-resistant strains of Mycobacterium tuberculosis in Ecuador

BACKGROUND

Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador.

RESULTS

Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected.

CONCLUSION

This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

Three highly variable genome regions of the four dengue virus serotypes can accurately recapitulate the CDS phylogeny

The circulation of the four-dengue virus (DENV) serotypes has significantly increased in recent years, accompanied by an increase in viral genetic diversity. In order to conduct disease surveillance and understand DENV evolution and its effects on virus transmission and disease, efficient and accurate methods for phylogenetic classification are required. Phylogenetic analysis of different viral genes sequences is the most used method, the envelope gene (E) being the most frequently selected target. We explored the genetic variability of the four DENV serotypes throughout their complete coding sequence (CDS) of sequences available in GenBank and used genomic regions of different variability rate to recapitulate the phylogeny obtained with the DENV CDS. Our results indicate that the use of high or low variable regions accurately recapitulate the phylogeny obtained with CDS of sequences from different DENV genotypes. However, when analyzing the phylogeny of a single genotype, highly variable regions performed better in recapitulating the distance branch length, topology, and support of the CDS phylogeny. The use of three concatenated highly variable regions was not statistically different in distance branch length and support to that obtained in CDS phylogeny.•This study demonstrated the ability of highly variable regions of the DENV genome to recapitulate the phylogeny obtained with the full coding sequence (CDS).•The use of genomic regions of high or low variability did not affect the performance in recapitulating the phylogeny obtained with CDS from different genotypes. However, when phylogeny was analyzed for sequences from a single genotype, highly variable regions performed better in recapitulating the distance branch length, topology, and support of the CDS phylogeny.•The use of concatenated highly variable genome regions represent a useful option for recapitulating genome-wide phylogenies in analyses of sequences belonging to the same DENV genotype.

Fading family lines- women and men without children, grandchildren and great-grandchildren in 19th, 20th and 21st Century Northern Sweden

We studied to what extent family lines die out over the course of 122 years based on Swedish population-level data. Our data included demographic and socioeconomic information for four generations in the Skellefteå region of northern Sweden from 1885 to 2007. The first generation in our sample consisted of men and women born between 1885 and 1899 (N = 5850), and we observed their children, grandchildren, and great-grandchildren. We found that 48% of the first generation did not have any living descendants (great-grandchildren) by 2007. The risk of a family line dying out within the four-generational framework was highest among those who had relatively low fertility in the first generation. Mortality during reproductive years was also a leading reason why individuals in the first generation ended up with a greater risk of not leaving descendants. We identified socioeconomic differences: both the highest-status and the lowest-status occupational groups saw an increased risk of not leaving any descendants. Almost all lineages that made it to the third generation also made it to the fourth generation.

Phylogenetic lineages of tuberculosis isolates and their association with patient demographics in Tanzania

BACKGROUND

Mycobacterium tuberculosis presents several lineages each with distinct characteristics of evolutionary status, transmissibility, drug resistance, host interaction, latency, and vaccine efficacy. Whole genome sequencing (WGS) has emerged as a new diagnostic tool to reliably inform the occurrence of phylogenetic lineages of Mycobacterium tuberculosis and examine their relationship with patient demographic characteristics and multidrug-resistance development.

METHODS

191 Mycobacterium tuberculosis isolates obtained from a 2017/2018 Tanzanian drug resistance survey were sequenced on the Illumina Miseq platform at Supranational Tuberculosis Reference Laboratory in Uganda. Obtained fast-q files were imported into tools for resistance profiling and lineage inference (Kvarq v0.12.2, Mykrobe v0.8.1 and TBprofiler v3.0.5). Additionally for phylogenetic tree construction, RaxML-NG v1.0.3(25) was used to generate a maximum likelihood phylogeny with 800 bootstrap replicates. The resulting trees were plotted, annotated and visualized using ggtree v2.0.4 RESULTS: Most [172(90.0%)] of the isolates were from newly treated Pulmonary TB patients. Coinfection with HIV was observed in 33(17.3%) TB patients. Of the 191 isolates, 22(11.5%) were resistant to one or more commonly used first line anti-TB drugs (FLD), 9(4.7%) isolates were MDR-TB while 3(1.6%) were resistant to all the drugs. Of the 24 isolates with any resistance conferring mutations, 13(54.2%) and 10(41.6%) had mutations in genes associated with resistance to INH and RIF respectively. The findings also show four major lineages i.e. Lineage 3[81 (42.4%)], followed by Lineage 4 [74 (38.7%)], the Lineage 1 [23 (12.0%)] and Lineages 2 [13 (6.8%)] circulaing in Tanzania.

CONCLUSION

The findings in this study show that Lineage 3 is the most prevalent lineage in Tanzania whereas drug resistant mutations were more frequent among isolates that belonged to Lineage 4.

Is Haemoproteus gabaldoni a valid species? An approach from morphology and molecular tools applied to parasites of Anseriformes

Currently, there are three recognized species of haemoproteids infecting Anseriformes: Haemoproteus nettionis, H. macrovacuolatus, and H. greineri. Unfortunately, genetic information associated with a morphotype is available only for H. macrovacuolatus. We recently found a parasite morphologically compatible with Haemoproteus gabaldoni, a species Bennet (1993) described in a Cairina moschata (Muscovy duck) from Venezuela. This species was synonymized to H. nettionis by Valkiūnas (2005), arguing not enough morphological differentiation between them; it was said that H. greineri could be as well a synonym of H. nettionis. In this study, we aimed to provide evidence to determine if Haemoproteus gabaldoni is a different species of H. nettionis and help to clarify other species status. We first performed morphological and morphometrical analyses and compared this information against the parahapantotypes of H. greineri, H. gabaldoni and material diagnosed as H. nettionis provided by the International Reference center for Avian Haematozoa (IRCAH), and H. macrovacuolatus from the Host-Parasite Relationship Study Group (GERPH, in Spanish Grupo de Estudio Relación Parásito Hospedero) biological collection. We used Principal Component Analysis (PCA) of dimensionless standard morphometrical variables from gametocytes. Furthermore, we amplified a small fragment of cytochrome b (cyt b) to compare the sequence with information in GenBank and Malavi through phylogenetic analyses and haplotype networks. PCA analyses revealed the presence of three distinct groups in the samples studied, supported in the morphological traits of each parasite species analyzed; phylogenetic analyses grouped parasite lineages separately according to the host and continent of provenance. Such results indicate that, H. gabaldoni, is a different species from H. nettionis. One more time, it is demonstrated the importance of linking barcode surveys to morphological studies. Finally, it is highlighted the importance of biological collections as repositories of worldwide biodiversity.

In silico analysis of mutant epitopes in new SARS-CoV-2 lineages suggest global enhanced CD8+ T cell reactivity and also signs of immune response escape

SARS-CoV-2 variants of concern have emerged since the COVID-19 outburst, notably the lineages detected in the UK, South Africa, and Brazil. Their increased transmissibility and higher viral load put them in the spotlight. Much has been investigated on the ability of those new variants to evade antibody recognition. However, little attention has been given to pre-existing and induced SARS-CoV-2-specific CD8+ T cell responses by new lineages. In this work, we predicted SARS-CoV-2-specific CD8+ T cell epitopes from the main variants of concern and their potential to trigger or hinder CD8+ T cell response by using HLA binding and TCR reactivity in silico predictions. Also, we estimated the population's coverage for different lineages, which accounts for the ability to present a set of peptides based on the most frequent HLA alleles of a given population. We considered binding predictions to 110 ccClass I HLA alleles from 29 countries to investigate differences in the fraction of individuals expected to respond to a given epitope set from new and previous lineages. We observed a higher population coverage for the variant detected in the UK (B.1.1.7), and South Africa (B.1.351), as well as for the Brazilian P.1 lineage, but not P.2, compared to the reference lineage. Moreover, individual mutations such as Spike N501Y and Nucleocapsid D138Y were predicted to have an overall stronger affinity through HLA-I than the reference sequence while Spike E484K shows signs of evasion. In summary, we provided evidence for the existence of potentially immunogenic and conserved epitopes across new SARS-CoV-2 variants, but also mutant peptides exhibiting diminished or abolished HLA-I binding. It also highlights the augmented population coverage for three new lineages. Whether these changes imply more T cell reactivity or potential to evade from CD8+ T cell responses requires experimental verification.

Era of TB elimination: Growing need to understand diversities of Mycobacterium tuberculosis lineages!

INTRODUCTION

The mycobacterium tuberculosis complex (MTBC) has highly clonal population structure which made the organism spread globally mirroring human migration out of Africa and resulted in the formation of seven lineages. We conducted this study to determine the proportion of spoligotype lineages and drug susceptibility profile of Mycobacterium tuberculosis isolates among smear positive TB patients attending a tertiary care hospital in Mysore, Karnataka, India.

METHODS

It is a descriptive study conducted at JSS Hospital a tertiary care centre at Mysore, India during 2018-19. The sputum smear positive samples were subjected to solid culture and drug susceptibility testing and spoligotyping for identification of lineages.

RESULTS

Of the 100 samples which were culture positive, 94 isolates were clustered into five spoligotype international types with SIT-126 (EAI-5) being the largest cluster of 46 (46%) isolates, followed by SIT-62 (H1) with 24 (24%), SIT -26 (CAS 1-DELHI) with 20 (20%), SIT-53 (T1) with 03 (3%) and SIT-482 (BOV-1) with 01 (1%). Among the remaining six isolates, two had unique Cameroon spoligotypes and four were orphans CONCLUSION: The study finding reveals that a diverse pattern of genotypes is circulating in the region of which EAI-5, Harleem (H1) and CAS-DELHI pattern forms the majority (88%). It is evident that there is a wide range of MTB genetic lineages in circulation and further research is needed to understand the diversity across the country.

Striking lineage diversity of severe acute respiratory syndrome coronavirus 2 from non-human sources

Due to the necessity to control human-to-human spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the overwhelming majority of the generated data on this virus was solely related to the genomic characteristics of strains infecting humans; conversely, this work aimed to recover and analyze the diversity of viral genomes from non-human sources. From a set of 3595 publicly available SARS-CoV-2 genome sequences, 128 lineages were identified in non-human hosts, the majority represented by the variants of concern Delta (n = 1105, 30.7%) and Alpha (n = 466, 12.9%), followed by B.1.1.298 lineage (n = 458, 12.7%). Environment, Neovison vison, Odocoileus virginianus and Felis catus were the non-human sources with the highest number of lineages (14, 12 and 10, respectively). Phylogenomic analyses showed viral clusters from environmental sources, N. vison, O. virginianus, Panthera tigris, and Panthera leo. These clusters were collectively related to human viruses as well as all other non-human sources that were heterogeneously distributed in the phylogenetic tree. Further, the genetic details of viral genomes from bats and pangolins were independently investigated owing to their high divergence, revealing five distinct clusters. Cluster 4 exclusively included bat-sourced genomes and the SARS-CoV-2 reference strain Wuhan-01. In summary, this study identified new genetic landmarks of SARS-CoV-2 evolution. We propose potential interspecies transmission routes of SARS-CoV-2 between animals and humans, which should be considered in order to establish better pathogen surveillance and containment strategies.

Structural and functional insights into the spike protein mutations of emerging SARS-CoV-2 variants

Since the emergence of the first case of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus (SARS-CoV-2), the viral genome has constantly undergone rapid mutations for better adaptation in the host system. These newer mutations have given rise to several lineages/ variants of the virus that have resulted in high transmission and virulence rates compared to the previously circulating variants. Owing to this, the overall caseload and related mortality have tremendously increased globally to > 233 million infections and > 4.7 million deaths as of Sept. 28th, 2021. SARS-CoV-2, Spike (S) protein binds to host cells by recognizing human angiotensin-converting enzyme 2 (hACE2) receptor. The viral S protein contains S1 and S2 domains that constitute the binding and fusion machinery, respectively. Structural analysis of viral S protein reveals that the virus undergoes conformational flexibility and dynamicity to interact with the hACE2 receptor. The SARS-CoV-2 variants and mutations might be associated with affecting the conformational plasticity of S protein, potentially linked to its altered affinity, infectivity, and immunogenicity. This review focuses on the current circulating variants of SARS-CoV-2 and the structure-function analysis of key S protein mutations linked with increased affinity, higher infectivity, enhanced transmission rates, and immune escape against this infection.

A rigorous framework for detecting SARS-CoV-2 spike protein mutational ensemble from genomic and structural features

The recent release of SARS-CoV-2 genomic data from several countries has provided clues into the potential antigenic drift of the coronavirus population. In particular, the genomic instability observed in the spike protein necessitates immediate action and further exploration in the context of viral-host interactions. By temporally tracking 527,988 SARS-CoV-2 genomes, we identified invariant and hypervariable regions within the spike protein. We evaluated combination of mutations from SARS-CoV-2 lineages and found that maximum number of lineage-defining mutations were present in the N-terminal domain (NTD). Based on mutant 3D-structural models of known Variants of Concern (VOCs), we found that structural properties such as accessibility, secondary structural type, and intra-protein interactions at local mutation sites are greatly altered. Further, we observed significant differences between intra-protein networks of wild-type and Delta mutant, with the latter showing dense intra-protein contacts. Extensive molecular dynamics simulations of D614G mutant spike structure with hACE2 further revealed dynamic features with 47.7% of mutations mapping on flexible regions of spike protein. Thus, we propose that significant changes within spike protein structure have occurred that may impact SARS-CoV-2 pathogenesis, and repositioning of vaccine candidates is required to contain the spread of COVID-19 pathogen.