Applying Next-Generation Sequencing to Track HIV-1 Drug Resistance Mutations Circulating in Portugal.

BACKGROUND: The global scale-up of antiretroviral treatment (ART) offers significant health benefits by suppressing HIV-1 replication and increasing CD4 cell counts. However, incomplete viral suppression poses a potential threat for the emergence of drug resistance mutations (DRMs), limiting ART options, and increasing HIV transmission. OBJECTIVE: We investigated the patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) among HIV-1 patients in Portugal. METHODS: Data were obtained from 1050 HIV-1 patient samples submitted for HIV drug resistance (HIVDR) testing from January 2022 to June 2023. Evaluation of DRM affecting viral susceptibility to nucleoside/tide reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), protease inhibitors (PIs), and integrase strand transfer inhibitors (INSTIs) was performed using an NGS technology, the Vela Diagnostics Sentosa SQ HIV-1 Genotyping Assay. RESULTS: About 71% of patients were ART naïve and 29% were experienced. Overall, 20% presented with any DRM. The prevalence of TDR and ADR was 12.6% and 41.1%, respectively. M184V, T215S, and M41L mutations for NRTI, K103N for NNRTI, and M46I/L for PIs were frequent in naïve and treated patients. E138K and R263K mutations against INSTIs were more frequent in naïve than treated patients. TDR and ADR to INSTIs were 0.3% and 7%, respectively. Patients aged 50 or over (OR: 1.81, p = 0.015), originating from Portuguese-speaking African countries (PALOPs) (OR: 1.55, p = 0.050), HIV-1 subtype G (OR: 1.78, p = 0.010), and with CD4 < 200 cells/mm3 (OR: 1.70, p = 0.043) were more likely to present with DRMs, while the males (OR: 0.63, p = 0.003) with a viral load between 4.1 to 5.0 Log10 (OR: 0.55, p = 0.003) or greater than 5.0 Log10 (OR: 0.52, p < 0.001), had lower chances of presenting with DRMs. CONCLUSIONS: We present the first evidence on TDR and ADR to INSTI regimens in followed up patients presenting for healthcare in Portugal. We observed low levels of TDR to INSTIs among ART-naïve and moderate levels in ART-exposed patients. Regimens containing PIs could be an alternative second line in patients with intermediate or high-level drug resistance, especially against second-generation INSTIs (dolutegravir, bictegravir, and cabotegravir).

Serological screening in a large-scale municipal survey in Cascais, Portugal, during the first waves of the COVID-19 pandemic: lessons for future pandemic preparedness efforts.

Background: Serological surveys for SARS-CoV-2 were used early in the COVID-19 pandemic to assess epidemiological scenarios. In the municipality of Cascais (Portugal), serological testing combined with a comprehensive socio-demographic, clinical and behavioral questionnaire was offered to residents between May 2020 and beginning of 2021. In this study, we analyze the factors associated with adherence to this municipal initiative, as well as the sociodemographic profile and chronic diseases clinical correlates associated to seropositivity. We aim to contribute with relevant information for future pandemic preparedness efforts. Methods: This was a cross-sectional study with non-probabilistic sampling. Citizens residing in Cascais Municipality went voluntarily to blood collection centers to participate in the serological survey. The proportion of participants, stratified by socio-demographic variables, was compared to the census proportions to identify the groups with lower levels of adherence to the survey. Univariate and multivariate logistic regression were used to identify socio-demographic, clinical and behavioral factors associated with seropositivity. Results: From May 2020 to February 2021, 19,608 participants (9.2% of the residents of Cascais) were included in the study. Based on the comparison to census data, groups with lower adherence to this survey were men, the youngest and the oldest age groups, individuals with lower levels of education and unemployed/inactive. Significant predictors of a reactive (positive) serological test were younger age, being employed or a student, and living in larger households. Individuals with chronic diseases generally showed lower seroprevalence. Conclusion: The groups with low adherence to this voluntary study, as well as the socio-economic contexts identified as more at risk of viral transmission, may be targeted in future pandemic situations. We also found that the individuals with chronic diseases, perceiving higher risk of serious illness, adopted protective behaviors that limited infection rates, revealing that health education on preventive measures was effective for these patients.

The Role of Late Presenters in HIV-1 Transmission Clusters in Europe.

BACKGROUND: Investigating the role of late presenters (LPs) in HIV-1 transmission is important, as they can contribute to the onward spread of HIV-1 virus before diagnosis, when they are not aware of their HIV status. OBJECTIVE: To characterize individuals living with HIV-1 followed up in Europe infected with subtypes A, B, and G and to compare transmission clusters (TC) in LP vs. non-late presenter (NLP) populations. METHODS: Information from a convenience sample of 2679 individuals living with HIV-1 was collected from the EuResist Integrated Database between 2008 and 2019. Maximum likelihood (ML) phylogenies were constructed using FastTree. Transmission clusters were identified using Cluster Picker. Statistical analyses were performed using R. RESULTS: 2437 (91.0%) sequences were from subtype B, 168 (6.3%) from subtype A, and 74 (2.8%) from subtype G. The median age was 39 y/o (IQR: 31.0-47.0) and 85.2% of individuals were males. The main transmission route was via homosexual (MSM) contact (60.1%) and 85.0% originated from Western Europe. In total, 54.7% of individuals were classified as LPs and 41.7% of individuals were inside TCs. In subtype A, individuals in TCs were more frequently males and natives with a recent infection. For subtype B, individuals in TCs were more frequently individuals with MSM transmission route and with a recent infection. For subtype G, individuals in TCs were those with a recent infection. When analyzing cluster size, we found that LPs more frequently belonged to small clusters (<8 individuals), particularly dual clusters (2 individuals). CONCLUSION: LP individuals are more present either outside or in small clusters, indicating a limited role of late presentation to HIV-1 transmission.

The International Virus Bioinformatics Meeting 2023.

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

HIV-1 subtype diversity and phylogenetic insight into non-B subtype transmission in Slovenia, 1989-2013.

INTRODUCTION: Disease progression, drug resistance mutations, and treatment strategies may vary by HIV-1 subtype. This study determined HIV-1 subtypes circulating in Slovenia, a Central European country with an HIV-1 epidemic driven by men who have sex with men, focusing on molecular epidemiology of non-B subtypes. METHODS: A total of 367 HIV-1 sequences were included. Subtype was assigned by employing eight different HIV subtyping tools coupled with maximum likelihood phylogenetic analyses. RESULTS: The subtyping tools COMET, jpHMM, and REGA 3.0 exhibited the best performance on the dataset studied. Phylogenetic analyses showed a 14.7% prevalence of non-B subtypes, with subtype A detected most frequently (4.9%), followed by CRF02_AG (2.4%), subtype C (1.1%), subtypes D, G, and CRF01_AE (0.8% each), and subtypes F and CRF22_01A1 (0.3% each). A subtype could not be assigned to 12 sequences (3.3%), indicating potential unique recombinant forms. Non-B subtypes were significantly associated with a heterosexual route of transmission and infection acquired in Eastern Europe, Africa, or Asia. CONCLUSION: In a country where subtype B is predominant, non-B subtypes were observed in one out of seven patients, a non-negligible proportion, which underlines the importance of systematic surveillance of HIV subtype diversity and the corresponding molecular epidemiology.

A Tale of Three Recent Pandemics: Influenza, HIV and SARS-CoV-2.

Emerging infectious diseases are one of the main threats to public health, with the potential to cause a pandemic when the infectious agent manages to spread globally. The first major pandemic to appear in the 20th century was the influenza pandemic of 1918, caused by the influenza A H1N1 strain that is characterized by a high fatality rate. Another major pandemic was caused by the human immunodeficiency virus (HIV), that started early in the 20th century and remained undetected until 1981. The ongoing HIV pandemic demonstrated a high mortality and morbidity rate, with discrepant impacts in different regions around the globe. The most recent major pandemic event, is the ongoing pandemic of COVID-19, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has caused over 5.7 million deaths since its emergence, 2 years ago. The aim of this work is to highlight the main determinants of the emergence, epidemic response and available countermeasures of these three pandemics, as we argue that such knowledge is paramount to prepare for the next pandemic. We analyse these pandemics' historical and epidemiological contexts and the determinants of their emergence. Furthermore, we compare pharmaceutical and non-pharmaceutical interventions that have been used to slow down these three pandemics and zoom in on the technological advances that were made in the progress. Finally, we discuss the evolution of epidemiological modelling, that has become an essential tool to support public health policy making and discuss it in the context of these three pandemics. While these pandemics are caused by distinct viruses, that ignited in different time periods and in different regions of the globe, our work shows that many of the determinants of their emergence and countermeasures used to halt transmission were common. Therefore, it is important to further improve and optimize such approaches and adapt it to future threatening emerging infectious diseases.

Trends of Transmitted and Acquired Drug Resistance in Europe From 1981 to 2019: A Comparison Between the Populations of Late Presenters and Non-late Presenters.

Background: The increased use of antiretroviral therapy (ART) has decreased mortality and morbidity of HIV-1 infected people but increasing levels of HIV drug resistance threatens the success of ART regimens. Conversely, late presentation can impact treatment outcomes, health costs, and potential transmission of HIV. Objective: To describe the patterns of transmitted drug resistance (TDR) and acquired drug resistance (ADR) in HIV-1 infected patients followed in Europe, to compare its patterns in late presenters (LP) vs non-late presenters (NLP), and to analyze the most prevalent drug resistance mutations among HIV-1 subtypes. Methods: Our study included clinical, socio-demographic, and genotypic information from 26,973 HIV-1 infected patients from the EuResist Integrated Database (EIDB) between 1981 and 2019. Results: Among the 26,973 HIV-1 infected patients in the analysis, 11,581 (42.9%) were ART-naïve patients and 15,392 (57.1%) were ART-experienced. The median age was 37 (IQR: 27.0-45.0) years old and 72.6% were males. The main transmission route was through heterosexual contact (34.9%) and 81.7% of patients originated from Western Europe. 71.9% of patients were infected by subtype B and 54.8% of patients were classified as LP. The overall prevalence of TDR was 12.8% and presented an overall decreasing trend (p for trend < 0.001), the ADR prevalence was 68.5% also with a decreasing trend (p for trend < 0.001). For LP and NLP, the TDR prevalence was 12.3 and 12.6%, respectively, while for ADR, 69.9 and 68.2%, respectively. The most prevalent TDR drug resistance mutations, in both LP and NLP, were K103N/S, T215rev, T215FY, M184I/V, M41I/L, M46I/L, and L90M. Conclusion: Our study showed that the overall TDR (12.8%) and ADR (68.5%) presented decreasing trends during the study time period. For LP, the overall TDR was slightly lower than for NLP (12.3 vs 12.6%, respectively); while this pattern was opposite for ADR (LP slightly higher than NLP). We suggest that these differences, in the case of TDR, can be related to the dynamics of fixation of drug resistance mutations; and in the case of ADR with the more frequent therapeutic failure in LPs.

Spectrum of Atazanavir-Selected Protease Inhibitor-Resistance Mutations.

Ritonavir-boosted atazanavir is an option for second-line therapy in low- and middle-income countries (LMICs). We analyzed publicly available HIV-1 protease sequences from previously PI-naïve patients with virological failure (VF) following treatment with atazanavir. Overall, 1497 patient sequences were identified, including 740 reported in 27 published studies and 757 from datasets assembled for this analysis. A total of 63% of patients received boosted atazanavir. A total of 38% had non-subtype B viruses. A total of 264 (18%) sequences had a PI drug-resistance mutation (DRM) defined as having a Stanford HIV Drug Resistance Database mutation penalty score. Among sequences with a DRM, nine major DRMs had a prevalence >5%: I50L (34%), M46I (33%), V82A (22%), L90M (19%), I54V (16%), N88S (10%), M46L (8%), V32I (6%), and I84V (6%). Common accessory DRMs were L33F (21%), Q58E (16%), K20T (14%), G73S (12%), L10F (10%), F53L (10%), K43T (9%), and L24I (6%). A novel nonpolymorphic mutation, L89T occurred in 8.4% of non-subtype B, but in only 0.4% of subtype B sequences. The 264 sequences included 3 (1.1%) interpreted as causing high-level, 14 (5.3%) as causing intermediate, and 27 (10.2%) as causing low-level darunavir resistance. Atazanavir selects for nine major and eight accessory DRMs, and one novel nonpolymorphic mutation occurring primarily in non-B sequences. Atazanavir-selected mutations confer low-levels of darunavir cross resistance. Clinical studies, however, are required to determine the optimal boosted PI to use for second-line and potentially later line therapy in LMICs.

Genome-wide diversity of Zika virus: Exploring spatio-temporal dynamics to guide a new nomenclature proposal.

The Zika virus (ZIKV) disease caused a public health emergency of international concern that started in February 2016. The overall number of ZIKV-related cases increased until November 2016, after which it declined sharply. While the evaluation of the potential risk and impact of future arbovirus epidemics remains challenging, intensified surveillance efforts along with a scale-up of ZIKV whole-genome sequencing provide an opportunity to understand the patterns of genetic diversity, evolution, and spread of ZIKV. However, a classification system that reflects the true extent of ZIKV genetic variation is lacking. Our objective was to characterize ZIKV genetic diversity and phylodynamics, identify genomic footprints of differentiation patterns, and propose a dynamic classification system that reflects its divergence levels. We analysed a curated dataset of 762 publicly available sequences spanning the full-length coding region of ZIKV from across its geographical span and collected between 1947 and 2021. The definition of genetic groups was based on comprehensive evolutionary dynamics analyses, which included recombination and phylogenetic analyses, within- and between-group pairwise genetic distances comparison, detection of selective pressure, and clustering analyses. Evidence for potential recombination events was detected in a few sequences. However, we argue that these events are likely due to sequencing errors as proposed in previous studies. There was evidence of strong purifying selection, widespread across the genome, as also detected for other arboviruses. A total of 50 sites showed evidence of positive selection, and for a few of these sites, there was amino acid (AA) differentiation between genetic clusters. Two main genetic clusters were defined, ZA and ZB, which correspond to the already characterized 'African' and 'Asian' genotypes, respectively. Within ZB, two subgroups, ZB.1 and ZB.2, represent the Asiatic and the American (and Oceania) lineages, respectively. ZB.1 is further subdivided into ZB.1.0 (a basal Malaysia sequence sampled in the 1960s and a recent Indian sequence), ZB.1.1 (South-Eastern Asia, Southern Asia, and Micronesia sequences), and ZB.1.2 (very similar sequences from the outbreak in Singapore). ZB.2 is subdivided into ZB.2.0 (basal American sequences and the sequences from French Polynesia, the putative origin of South America introduction), ZB.2.1 (Central America), and ZB.2.2 (Caribbean and North America). This classification system does not use geographical references and is flexible to accommodate potential future lineages. It will be a helpful tool for studies that involve analyses of ZIKV genomic variation and its association with pathogenicity and serve as a starting point for the public health surveillance and response to on-going and future epidemics and to outbreaks that lead to the emergence of new variants.

Readdressing the genetic diversity and taxonomy of the Mesoniviridae family, as well as its relationships with other nidoviruses and putative mesonivirus-like viral sequences.

Research on the recently established Mesoniviridae family (Order Nidovirales), RNA genome insect-specific viruses, has been steadily growing in the last decade. However, after the last detailed phylogenetic characterization of mesoniviruses in 2014, numerous new sequences, even in organisms other than mosquitos, have been identified and characterized. In this study, we analyzed nucleotide and protein sequences of mesoniviruses with a wide range of molecular tools including genetic distance, Shannon entropy, selective pressure analysis, polymorphism identification, principal coordinate analysis, likelihood mapping and phylodynamic reconstruction. We also sought to revaluate new mesoniviruses sequence positions within the family, proposing a taxonomic revision. The different sub-lineages of mosquito mesoniviruses sequences presented low sequence diversity and entropy, with incongruences to the existing taxonomy being found after an extensive phylogenetic characterization. High sequence discrepancy and differences in genome organization were found between mosquito mesoniviruses and other mesoniviruses, so their future classification, as other meso-like viruses that are found in other organisms, should be approached with caution. No evidence of frequent recombination was found, and mesonivirus genomes seem to evolve under strong purifying selection. Insufficient data by root-to-tip analysis did not yet allow for an adequate phylogeographic reconstruction.

Insect-specific viruses in the Parvoviridae family: Genetic lineage characterization and spatiotemporal dynamics of the recently established Brevihamaparvovirus genus.

The analysis of the viruses allocated to the recently established Brevihamaparvovirus genus (Parvoviridae family), which includes all previously known brevidensoviruses, has not yet been carried out on an extensive basis. As a result, no detailed genetic lineage characterization has ever been performed for this group of insect-specific viruses. Using a wide range of molecular tools, we have explored this taxon by calculating Shannon entropy values, intra- and inter-taxon genetic distances, analysed sequence polymorphisms, and evaluated selective pressures acting on the viral genome. While the calculated Brevihamaparvovirus mutation rates were within the range of those of other parvoviruses, their genomes look to be under strong purifying selection, and are also characterized by low diversity and entropy. Furthermore, even though recombination events are quite common among parvoviruses, no evidence of recombination (either intra or intergenic) was found in the Brevihamaparvoviruses sequences analyzed. An extended taxonomic analysis and reevaluation of existing Brevihamaparvoviruses sequences, many still unclassified, was performed using cut-off values defining NS1 identity between viral sequences from the Parvovirus family. Two existing genetic lineages, Dipteran Brevihamaparvovirus 1 and Dipteran Brevihamaparvovirus 2, were rearranged and the creation of a new one, Dipteran Brevihamaparvovirus 3, was suggested. Finally, despite the uncertainties associated with both the time estimates of the most recent common ancestors, which could span from twenty thousand years before the current era to way earlier (in the last century), and the dispersal routes proposed for Brevihamaparvoviruses sequences by phylodynamic reconstruction, the analyses here presented could help define how future studies should be conducted as more isolates continue to be identified in the future, and contribute to eliminating possible analytical biases.

Differential patterns of postmigration HIV-1 infection acquisition among Portuguese immigrants of different geographical origins.

OBJECTIVE: To investigate the dynamics of phylogenetic transmission clusters involving immigrants of Portuguese Speaking Countries living in Portugal. DESIGN/METHODS: We included genomic sequences, sociodemographic and clinical data from 772 HIV migrants followed in Portugal between 2001 and 2017. To reconstruct HIV-1 transmission clusters, we applied phylogenetic inference from 16 454 patients: 772 migrants, 2973 Portuguese and 12 709 global controls linked to demographic and clinical data. Transmission clusters were defined using: clusters with SH greater than 90% (phylogenetic support), genetic distance less than 3.5% and clusters that included greater than 66% of patients from one specific geographic origin compared with the total of sequences within the cluster. Logistic regression was performed to assess factors associated with clustering. RESULTS: Three hundred and six (39.6%) of migrants were included in transmission clusters. This proportion differed substantially by region of origin [Brazil 54% vs. Portuguese Speaking African Countries (PALOPs) 36%, P < 0.0001] and HIV-1 infecting subtype (B 52%, 43% subtype G and 32% CRF02_AG, P < 0.001). Belonging to a transmission cluster was independently associated with treatment-naive patients, CD4+ greater than 500, with recent calendar years of sampling, origin from PALOPs and with seroconversion. Among Brazilian migrants - mainly infected with subtype B - 40.6% were infected by Portuguese. Among migrants from PALOPs - mainly infected with subtypes G and CFR02_AG - the transmission occurred predominantly within the migrants' community (53 and 80%, respectively). CONCLUSION: The acquisition of infection among immigrants living in Portugal differs according to the country of origin. These results can contribute to monitor the HIV epidemic and prevent new HIV infections among migrants.

Genetic lineage characterization and spatiotemporal dynamics of classical insect-specific flaviviruses: outcomes and limitations.

The genus Flavivirus incorporates bona fide arboviruses, as well as others viruses with restricted replication in insect cells. Among the latter, a large monophyletic cluster of viruses, known as cISF (classical insect-specific flaviviruses), has been sampled in many species of mosquitoes collected over a large geographic range. In this study, we investigated nucleotide and protein sequences with a suite of molecular characterization approaches including genetic distance, Shannon entropy, selective pressure analysis, polymorphism identification, principal coordinate analysis, likelihood mapping, phylodynamic reconstruction, and spatiotemporal dispersal, to further characterize this diverse group of insect-viruses. The different lineages and sub-lineages of viral sequences presented low sequence diversity and entropy (though some displayed lineage-specific polymorphisms), did not show evidence of frequent recombination and evolved under strong purifying selection. Moreover, the reconstruction of the evolutionary history and spatiotemporal dispersal was highly impacted by overall low signals of sequence divergence throughout time but suggested that cISF distribution in space and time is dynamic and may be dependent on human activities, including commercial trading and traveling.

Short Communication: HIV and Tuberculosis Co-Infection Among Migrants in Portugal: A Brief Study on Their Sociodemographic, Clinical, and Genomic Characteristics.

HIV and tuberculosis (TB) are among the global deadliest diseases. Migrant populations are particularly vulnerable to these infections. Yet, literature is still scarce on the epidemiology of HIV-TB co-infection among migrants. In this study, we characterized native and migrant HIV patients followed in Portuguese hospitals, who were diagnosed with TB, regarding their sociodemographic, clinical, and genomic characteristics. Among 67 patients with HIV and TB diagnoses, there were 24 migrants, most from sub-Saharan Africa. Most patients had CD4+ T cell counts below 350 cells/µL, and were diagnosed simultaneously for HIV and TB. When compared to natives, migrants presented a higher proportion of non-B HIV-1 infections. Patients infected with these non-B HIV-1 strains presented higher viral loads, which can have an important impact for the transmissibility and pathogenicity of both diseases. Future studies should investigate different HIV strains and consequences for TB and HIV transmission and disease outcomes, especially among vulnerable populations.

An Evolutionary Model-Based Approach To Quantify the Genetic Barrier to Drug Resistance in Fast-Evolving Viruses and Its Application to HIV-1 Subtypes and Integrase Inhibitors.

Viral pathogens causing global disease burdens are often characterized by high rates of evolutionary changes. The extensive viral diversity at baseline can shorten the time to escape from therapeutic or immune selective pressure and alter mutational pathways. The impact of genotypic background on the barrier to resistance can be difficult to capture, particularly for agents in experimental stages or that are recently approved or expanded into new patient populations. We developed an evolutionary model-based counting method to quickly quantify the population genetic potential to resistance and assess population differences. We demonstrate its applicability to HIV-1 integrase inhibitors, as their increasing use globally contrasts with limited availability of non-B subtype resistant sequence data and corresponding knowledge gap. A large sequence data set encompassing most prevailing HIV-1 subtypes and resistance-associated mutations of currently approved integrase inhibitors was investigated. A complex interplay between codon predominance, polymorphisms, and associated evolutionary costs resulted in a subtype-dependent varied genetic potential for 15 resistance mutations against integrase inhibitors. While we confirm the lower genetic barrier of subtype B for G140S, we convincingly discard a similar effect previously suggested for G140C. A supplementary analysis for HIV-1 reverse transcriptase inhibitors identified a lower genetic barrier for K65R in subtype C through differential codon usage not reported before. To aid evolutionary interpretations of genomic differences for antiviral strategies, we advanced existing counting methods with increased sensitivity to identify subtype dependencies of resistance emergence. Future applications include novel HIV-1 drug classes or vaccines, as well as other viral pathogens.

Characterization of a large cluster of HIV-1 A1 infections detected in Portugal and connected to several Western European countries.

HIV-1 subtypes associate with differences in transmission and disease progression. Thus, the existence of geographic hotspots of subtype diversity deepens the complexity of HIV-1/AIDS control. The already high subtype diversity in Portugal seems to be increasing due to infections with sub-subtype A1 virus. We performed phylogenetic analysis of 65 A1 sequences newly obtained from 14 Portuguese hospitals and 425 closely related database sequences. 80% of the A1 Portuguese isolates gathered in a main phylogenetic clade (MA1). Six transmission clusters were identified in MA1, encompassing isolates from Portugal, Spain, France, and United Kingdom. The most common transmission route identified was men who have sex with men. The origin of the MA1 was linked to Greece, with the first introduction to Portugal dating back to 1996 (95% HPD: 1993.6-1999.2). Individuals infected with MA1 virus revealed lower viral loads and higher CD4+ T-cell counts in comparison with those infected by subtype B. The expanding A1 clusters in Portugal are connected to other European countries and share a recent common ancestor with the Greek A1 outbreak. The recent expansion of this HIV-1 subtype might be related to a slower disease progression leading to a population level delay in its diagnostic.

Geographic dispersal and genetic diversity of tick-borne phleboviruses (Phenuiviridae, Phlebovirus) as revealed by the analysis of L segment sequences.

The large diversity of new tick-borne phleboviruses, and the negative impacts of the virulent viruses on human/animal health have led to a growing interest in their analysis. In this report, new insights are brought out into the diversity of putative phleboviruses circulating in Portugal (both the continental territory and the islands of São Miguel, in the Azores, and Madeira), as well as in the Spanish western regions of Extremadura and Castilla and León. Phlebovirus sequences were frequently detected (L-segment) from both questing and feeding ticks, but especially in Rhipicephalus sanguineus sensu lato (s.l.) specimens. These sequences were detected in adult ticks, as well as nymphs and eggs, supporting the hypothesis of viral maintenance by vertical transmission. Though multiple genetic groups could be identified in phylogenetic trees (AnLuc, KarMa, RiPar virus 1, and Spanish group 1 and 2), all the sequences from Portugal and Spain shared common ancestry with other viral sequence obtained from samples collected over a large geographic coverage. Spatiotemporal analysis placed Middle-East as the geographic origin of the most recent common ancestor (MRCA) of all phleboviruses analysed in the present study. More recent viral transitions might include migrations from Spain to continental Portugal, and from there to the Portuguese Islands. Our findings suggest that the time of the MRCA of phleboviruses was dated around 225 years ago [95% HPD: 124-387 year before the last sampling date].

Tuberculosis care for migrant patients in Portugal: a mixed methods study with primary healthcare providers.

BACKGROUND: Tuberculosis (TB) is still a major global health problem. The increasing number of cases observed among foreign-born populations contrasts with the decreasing trends observed in later years in some high-income countries. Healthcare providers are key interveners in the control of TB and HIV-TB infections. In this study, we aimed to explore the perspectives of healthcare providers working in primary care in Portugal about the provision of TB care for migrant patients with TB or HIV-TB co-infection. METHODS: We applied a mixed-methods approach using an online survey and semi-structured interviews with primary healthcare providers. A total of 120 Portuguese healthcare providers participated in the survey, and 17 were interviewed. Survey and interview data were analysed applying descriptive statistics and thematic analysis, respectively. RESULTS: Migrants' lack of knowledge on TB disease and its symptoms was the main reason for advanced-stage presentation of cases. Their high mobility and social isolation affect adherence to treatment. The providers also listed several barriers to migrants' access and use of TB care. The most frequently referred were limited socioeconomic resources, complex bureaucracy at the point of access and registration for healthcare services, especially for undocumented migrants, and obstacles for social protection. Providers also advocated more training initiatives on migrants' health, social and cultural contexts, on HIV and TB integrated care, and on TB scientific update for general practitioners and nurses working at primary healthcare centres. CONCLUSIONS: Future efforts should provide measures to overcome social, economic and administrative obstacles to care for TB-infected migrants, and promote regular training initiatives for national healthcare providers in order to raise awareness and facilitate better care to culturally diverse populations with TB.

Tracing the Impact of Public Health Interventions on HIV-1 Transmission in Portugal Using Molecular Epidemiology.

BACKGROUND: Estimation of temporal changes in human immunodeficiency virus (HIV) transmission patterns can help to elucidate the impact of preventive strategies and public health policies. METHODS: Portuguese HIV-1 subtype B and G pol genetic sequences were appended to global reference data sets to identify country-specific transmission clades. Bayesian birth-death models were used to estimate subtype-specific effective reproductive numbers (Re). Discrete trait analysis (DTA) was used to quantify mixing among transmission groups. RESULTS: We identified 5 subtype B Portuguese clades (26-79 sequences) and a large monophyletic subtype G Portuguese clade (236 sequences). We estimated that major shifts in HIV-1 transmission occurred around 1999 (95% Bayesian credible interval [BCI], 1998-2000) and 2000 (95% BCI, 1998-2001) for subtypes B and G, respectively. For subtype B, Re dropped from 1.91 (95% BCI, 1.73-2.09) to 0.62 (95% BCI,.52-.72). For subtype G, Re decreased from 1.49 (95% BCI, 1.39-1.59) to 0.72 (95% BCI, .63-.8). The DTA suggests that people who inject drugs (PWID) and heterosexuals were the source of most (>80%) virus lineage transitions for subtypes G and B, respectively. CONCLUSIONS: The estimated declines in Re coincide with the introduction of highly active antiretroviral therapy and the scale-up of harm reduction for PWID. Inferred transmission events across transmission groups emphasize the importance of prevention efforts for bridging populations.

VIRULIGN: fast codon-correct alignment and annotation of viral genomes.

SUMMARY: Virus sequence data are an essential resource for reconstructing spatiotemporal dynamics of viral spread as well as to inform treatment and prevention strategies. However, the potential benefit of these applications critically depends on accurate and correctly annotated alignments of genetically heterogeneous data. VIRULIGN was built for fast codon-correct alignments of large datasets, with standardized and formalized genome annotation and various alignment export formats. AVAILABILITY AND IMPLEMENTATION: VIRULIGN is freely available at https://github.com/rega-cev/virulign as an open source software project. SUPPLEMENTARY INFORMATION: Supplementary data is available at Bioinformatics online.