Metagenomic analysis reveals differences in the co-occurrence and abundance of viral species in SARS-CoV-2 patients with different severity of disease

Background

SARS-CoV-2 infections have a wide spectrum of clinical manifestations whose causes are not completely understood. Some human conditions predispose to severe outcome, like old age or the presence of comorbidities, but many other facets, including coinfections with other viruses, remain poorly characterized.

Methods

In this study, the eukaryotic fraction of the respiratory virome of 120 COVID-19 patients was characterized through whole metagenomic sequencing.

Results

Genetic material from respiratory viruses was detected in 25% of all samples, whereas human viruses other than SARS-CoV-2 were found in 80% of them. Samples from hospitalized and deceased patients presented a higher prevalence of different viruses when compared to ambulatory individuals. Small circular DNA viruses from the Anneloviridae (Torque teno midi virus 8, TTV-like mini virus 19 and 26) and Cycloviridae families (Human associated cyclovirus 10), Human betaherpesvirus 6, were found to be significantly more abundant in samples from deceased and hospitalized patients compared to samples from ambulatory individuals. Similarly, Rotavirus A, Measles morbillivirus and Alphapapilomavirus 10 were significantly more prevalent in deceased patients compared to hospitalized and ambulatory individuals.

Conclusions

Results show the suitability of using metagenomics to characterize a broader peripheric virological landscape of the eukaryotic virome in SARS-CoV-2 infected patients with distinct disease outcomes. Identified prevalent viruses in hospitalized and deceased patients may prove important for the targeted exploration of coinfections that may impact prognosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07783-8.

Phylodynamics of HIV in the Mexico City Metropolitan Region

Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.

The Alpha Variant (B.1.1.7) of SARS-CoV-2 Failed to Become Dominant in Mexico

During the coronavirus disease 2019 (COVID-19) pandemic, the emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in the United Kingdom in September 2020, was well documented in different areas of the world and became a global public health concern because of its increased transmissibility. The B.1.1.7 lineage was first detected in Mexico during December 2020, showing a slow progressive increase in its circulation frequency, which reached its maximum in May 2021 but never became predominant. In this work, we analyzed the patterns of diversity and distribution of this lineage in Mexico using phylogenetic and haplotype network analyses. Despite the reported increase in transmissibility of the B.1.1.7 lineage, in most Mexican states, it did not displace cocirculating lineages, such as B.1.1.519, which dominated the country from February to May 2021. Our results show that the states with the highest prevalence of B.1.1.7 were those at the Mexico-U.S. border. An apparent pattern of dispersion of this lineage from the northern states of Mexico toward the center or the southeast was observed in the largest transmission chains, indicating possible independent introduction events from the United States. However, other entry points cannot be excluded, as shown by multiple introduction events. Local transmission led to a few successful haplotypes with a localized distribution and specific mutations indicating sustained community transmission. IMPORTANCE The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021. In this work, we analyzed the distribution of B.1.1.7 in Mexico using phylogenetic and haplotype network analyses. Our results show that the states with the highest prevalence of B.1.1.7 (around 30%) were those at the Mexico-U.S. border, which also exhibited the highest lineage diversity, indicating possible introduction events from the United States. Also, several haplotypes were identified with a localized distribution and specific mutations, indicating that sustained community transmission occurred in the country.

HIV Pretreatment Drug Resistance Trends in Mexico City, 2017-2020

In response to increasing pretreatment drug resistance (PDR), Mexico changed its national antiretroviral treatment (ART) policy, recommending and procuring second-generation integrase strand-transfer inhibitor (INSTI)-based regimens as preferred first-line options since 2019. We present a four-year observational study describing PDR trends across 2017-2020 at the largest HIV diagnosis and primary care center in Mexico City. A total of 6688 baseline protease-reverse transcriptase and 6709 integrase sequences were included. PDR to any drug class was 14.4% (95% CI, 13.6-15.3%). A significant increasing trend for efavirenz/nevirapine PDR was observed (10.3 to 13.6%, p = 0.02). No increase in PDR to second-generation INSTI was observed, remaining under 0.3% across the study period. PDR was strongly associated with prior exposure to ART (aOR: 2.9, 95% CI: 1.9-4.6, p < 0.0001). MSM had higher odds of PDR to efavirenz/nevirapine (aOR: 2.0, 95% CI: 1.0-3.7, p = 0.04), reflecting ongoing transmission of mutations such as K103NS and E138A. ART restarters showed higher representation of cisgender women and injectable drug users, higher age, and lower education level. PDR to dolutegravir/bictegravir remained low in Mexico City, although further surveillance is warranted given the short time of ART optimization. Our study identifies demographic characteristics of groups with higher risk of PDR and lost to follow-up, which may be useful to design differentiated interventions locally.

Genetic Analysis of SARS-CoV-2 Variants in Mexico during the First Year of the COVID-19 Pandemic

During the first year of the SARS-CoV-2 pandemic in Mexico, more than two million people were infected. In this study, we analyzed full genome sequences from 27 February 2020 to 28 February 2021 to characterize the geographical and temporal distribution of SARS-CoV-2 lineages and identify the most common circulating lineages during this period. We defined six different geographical regions with particular dynamics of lineage circulation. The Northeast and Northwest regions were the ones that exhibited the highest lineage diversity, while the Central south and South/Southeast regions presented less diversity with predominance of a certain lineage. Additionally, by late February 2021, lineage B.1.1.519 represented more than 89% of all circulating lineages in the country.

Emergence and spread of the potential variant of interest (VOI) B.1.1.519 of SARS-CoV-2 predominantly present in Mexico

SARS-CoV-2 variants emerged in late 2020, and at least three variants of concern (B.1.1.7, B.1.351, and P1) have been reported by WHO. These variants have several substitutions in the spike protein that affect receptor binding; they exhibit increased transmissibility and may be associated with reduced vaccine effectiveness. In the present work, we report the identification of a potential variant of interest, harboring the mutations T478K, P681H, and T732A in the spike protein, within the newly named lineage B.1.1.519, that rapidly outcompeted the preexisting variants in Mexico and has been the dominant virus in the country during the first trimester of 2021.

Pretreatment HIV drug resistance spread within transmission clusters in Mexico City

BACKGROUND

Pretreatment HIV drug resistance (HIVDR) to NNRTIs has consistently increased in Mexico City during the last decade.

OBJECTIVES

To infer the HIV genetic transmission network in Mexico City to describe the dynamics of the local HIV epidemic and spread of HIVDR.

PATIENTS AND METHODS

HIV pol sequences were obtained by next-generation sequencing from 2447 individuals before initiation of ART at the largest HIV clinic in Mexico City (April 2016 to June 2018). Pretreatment HIVDR was estimated using the Stanford algorithm at a Sanger-like threshold (≥20%). Genetic networks were inferred with HIV-TRACE, establishing putative transmission links with genetic distances <1.5%. We examined demographic associations among linked individuals with shared drug resistance mutations (DRMs) using a ≥ 2% threshold to include low-frequency variants.

RESULTS

Pretreatment HIVDR reached 14.8% (95% CI 13.4%-16.2%) in the cohort overall and 9.6% (8.5%-10.8%) to NNRTIs. Putative links with at least one other sequence were found for 963/2447 (39%) sequences, forming 326 clusters (2-20 individuals). The inferred network was assortative by age and municipality (P < 0.001). Clustering individuals were younger [adjusted OR (aOR) per year = 0.96, 95% CI 0.95-0.97, P < 0.001] and less likely to include women (aOR = 0.46, 95% CI 0.28-0.75, P = 0.002). Among clustering individuals, 175/963 (18%) shared DRMs (involving 66 clusters), of which 66/175 (38%) shared K103N/S (24 clusters). Eight municipalities (out of 75) harboured 65% of persons sharing DRMs. Among all persons sharing DRMs, those sharing K103N were younger (aOR = 0.93, 95% CI 0.88-0.98, P = 0.003).

CONCLUSIONS

Our analyses suggest age- and geographically associated transmission of DRMs within the HIV genetic network in Mexico City, warranting continuous monitoring and focused interventions.

Longitudinal variation in human immunodeficiency virus long terminal repeat methylation in individuals on suppressive antiretroviral therapy

BACKGROUND

Persistence of latent, replication-competent provirus in CD4+ T cells of human immunodeficiency virus (HIV)-infected individuals on antiretroviral treatment (ART) is the main obstacle for virus eradication. Methylation of the proviral 5' long terminal repeat (LTR) promoter region has been proposed as a possible mechanism contributing to HIV latency; however, conflicting observations exist regarding its relevance. We assessed 5'-LTR methylation profiles in total CD4+ T cells from blood of 12 participants on short-term ART (30 months) followed up for 2 years, and a cross-sectional group of participants with long-term ART (6-15 years), using next generation sequencing. We then looked for associations between specific 5'-LTR methylation patterns and baseline and follow-up clinical characteristics.

RESULTS

5'-LTR methylation was observed in all participants and behaved dynamically. The number of 5'-LTR variants found per sample ranged from 1 to 13, with median sequencing depth of 16270× (IQR 4107×-46760×). An overall significant 5'-LTR methylation increase was observed at month 42 compared to month 30 (median CpG Methylation Index: 74.7% vs. 0%, p = 0.025). This methylation increase was evident in a subset of participants (methylation increase group), while the rest maintained fairly high and constant methylation (constant methylation group). Persons in the methylation increase group were younger, had higher CD4+ T cell gain, larger CD8% decrease, and larger CD4/CD8 ratio change after 48 months on ART (all p < 0.001). Using principal component analysis, the constant methylation and methylation increase groups showed low evidence of separation along time (factor 2: p = 0.04). Variance was largely explained (21%) by age, CD4+/CD8+ T cell change, and CD4+ T cell subpopulation proportions. Persons with long-term ART showed overall high methylation (median CpG Methylation Index: 78%; IQR 71-87%). No differences were observed in residual plasma viral load or proviral load comparing individuals on short-term (both at 30 or 42 months) and long-term ART.

CONCLUSIONS

Our study shows evidence that HIV 5'-LTR methylation in total CD4+ T cells is dynamic along time and that it can follow different temporal patterns that are associated with a combination of baseline and follow-up clinical characteristics. These observations may account for differences observed between previous contrasting studies.

Pretreatment HIV-drug resistance in Mexico and its impact on the effectiveness of first-line antiretroviral therapy: a nationally representative 2015 WHO survey

BACKGROUND

WHO has developed a global HIV-drug resistance surveillance strategy, including assessment of pretreatment HIV-drug resistance. We aimed to do a nationally representative survey of pretreatment HIV-drug resistance in Mexico using WHO-recommended methods.

METHODS

Among 161 Ministry of Health antiretroviral therapy (ART) clinics in Mexico, the largest, including 90% of ART initiators within the Ministry of Health (66 in total), were eligible for the survey. We used a probability-proportional-to-size design method to sample 25 clinics throughout the country. Consecutive ART-naive patients with HIV about to initiate treatment were invited to participate in the survey; individuals with previous exposure to ART were excluded. We assessed pretreatment HIV-drug resistance by Sanger sequencing and next-generation sequencing of viruses from plasma specimens from eligible participants with Stanford University HIV Drug Resistance Database methods. We obtained follow-up data for a median of 9·4 months (range 6-12) after enrolment. We investigated possible relations between demographic variables and pretreatment drug resistance with univariate and multivariate logistic regression.

FINDINGS

Between Feb 3 and July 30, 2015, we screened 288 patients in 25 clinics, from whom 264 provided successfully sequenced viruses with no evidence of current exposure to antiretroviral drugs. With the Sanger method, of these 264 participants, 41 (15·5%, 95% CI 11·4-20·5) had pretreatment resistance to any antiretroviral drug and 28 (10·6%, 7·2-15·0) had pretreatment resistance to non-nucleoside reverse transcriptase inhibitors (NNRTIs). At least low-level pretreatment resistance (Stanford penalty score ≥15) was noted in 13 (4 · 9%) of participants to efavirenz and in 23 (8·7%) to the combination tenofovir plus emtricitabine plus efavirenz. With next-generation sequencing, of 264 participants, 38 (14·4%, 95% CI 10·4-19·2) had pretreatment resistance to any antiretroviral drug and 26 (9·8%, 6·5-14·1) had pretreatment resistance to NNRTIs. After median follow-up of 8 months (IQR 6·5-9·4, range 5-11) after ART initiation, 97 (72%) of 135 NNRTI initiators achieved viral suppression (<50 copies per mL) compared with ten (40%) of 25 individuals who started with protease inhibitor-based regimens (p=0·0045). After multivariate regression considering pretreatment resistance and initial ART regimen as composite variables, people starting NNRTIs with pretreatment drug resistance achieved significantly lower viral suppression (odds ratio 0·24, 95% CI 0·07-0·74; p=0·014) than patients without NNRTI resistance.

INTERPRETATION

High levels of pretreatment drug resistance were noted in Mexico, and NNRTI pretreatment drug resistance significantly reduced the effectiveness of first-line ART regimens based on these drugs. Baseline HIV-drug resistance testing for initial ART follow-up and decision making should be considered.

FUNDING

The Mexican Government and Consejo Nacional de Ciencia y Tecnología.

[Genetic characterization of adenovirus isolated from follicular conjunctivitis and epidemic keratoconjunctivitis in a group of Mexican patients]

OBJECTIVE

To determine the adenovirus serotype in Mexican patients with folicular conjunctivitis and keratoconjunctivitis.

METHODS

Adenovirus-specific PCR was used to analyze sample scrapings from the inferior fornix of patients with follicular conjunctivitis and clinical suspicion of adenovirus from January 2005 to December 2006. Identification of the serotype was made by automated sequencing. The nucleotide sequences obtained were compared with the reported sequences in GenBank. Descriptive statistical analyses were performed on the results.

RESULTS

Of the 77 samples with clinical data of follicular conjunctivitis that were analyzed, 43 (56%) presented adenovirus. The sequencing of each positive sample allowed the identification of Ad1, Ad2, Ad3 and Ad8; the sequences of the serotype were identical those reported in GenBank with accession numbers: AF 534906 and AY 224420 for a sequence of the gene coding for the filament of Ad1 and Ad2 respectively, and AY 854180 and DQ 149614 for a sequence of the gene that codes for the Hexon protein of Ad3 and Ad8 respectively. From the statistical analysis it was possible to determine that a preferential seasonality of the serotype does not exist.

CONCLUSION

In this work the Ad1, Ad2 and Ad3 serotypes were identified in patients with clinical diagnosis of follicular conjunctivitis in 2005. Ad2 was the predominant serotype. Ad8 was also detected in an outbreak of epidemic keratoconjunctivitis. From an epidemiological point of view, no serotype found seems to have a preferred seasonality.

[Identification of adenovirus associated with conjunctivitis by molecular methodology]

OBJECTIVE

To identify adenovirus in patients with conjunctivitis by molecular methods such as the Polymerase Chain Reaction (PCR) and DNA sequencing.

METHODS

Samples of scrapings from the inferior fornix of 51 patients (39 diagnosed with Follicular Conjunctivitis and 12 diagnosed with Vernal Conjunctivitis) were processed by generic PCR to identify adenovirus. All the samples that were PCR positive were cultured on VERO cells for virus isolation, with this being demonstrated by immunofluorescence. For the identification of the isolated serotype, the multiplex PCR was utilized and DNA automated sequencing was employed to identify the genetic variants.

RESULTS

Twenty-eight of the individuals diagnosed with Follicular Conjunctivitis and six of those diagnosed with Vernal Conjunctivitis, had positive results to adenovirus (67%). The cultures in VERO cells allowed the isolation of eight samples. Only three of the isolated viruses (one Ad1 and two Ad2) were identified by the multiplex PCR used to identify the subgenus C adenovirus. An Ad1 genetic variant was identified by automated sequencing while the Ad2 serotypes were identical to the ones reported by Genbank.

CONCLUSIONS

The Polymerase Chain Reaction and DNA sequencing are useful tools to identify and characterize microorganisms responsible for diseases such as conjunctivitis caused by adenoviruses.