Assessment of phylogenetic sensitivity for reconstructing HIV-1 epidemiological relationships

The T350G Variation of Human Papillomavirus 16 E6 Gene Prevails in Oropharyngeal Cancer from a Small Cohort of Greek Patients

Recent trends have shown a dramatic rise in the incidence of oropharyngeal squamous cell carcinoma strongly associated with high-risk human papillomavirus (HPV) of type 16. The genetic variability of HPV16 has been extensively studied in cervical cancer but there are very limited published data concerning the genetic variations of this HPV type in oropharyngeal cancer. In the present study, the genetic variations of HPV16 E6 gene sequences originated from a small cohort of Greek patients diagnosed with oropharyngeal cancer were assessed. The vast majority of the sequences clustered within the European variant branch. The T350G variation was found to be the predominant one. This finding may indicate the need for further studies that could explain the possible impact of this variant in the pathomechanisms of oropharyngeal cancer.

Animal Coronaviruses Induced Apoptosis

Apoptosis is a form of programmed death that has also been observed in cells infected by several viruses. It is considered one of the most critical innate immune mechanisms that limits pathogen proliferation and propagation before the initiation of the adaptive immune response. Recent studies investigating the cellular responses to SARS-CoV and SARS-CoV-2 infection have revealed that coronaviruses can alter cellular homeostasis and promote cell death, providing evidence that the modulation of apoptotic pathways is important for viral replication and propagation. Despite the genetic diversity among different coronavirus clades and the infection of different cell types and several hosts, research studies in animal coronaviruses indicate that apoptosis in host cells is induced by common molecular mechanisms and apoptotic pathways. We summarize and critically review current knowledge on the molecular aspects of cell-death regulation during animal coronaviruses infection and the viral–host interactions to this process. Future research is expected to lead to a better understanding of the regulation of cell death during coronavirus infection. Moreover, investigating the role of viral proteins in this process will help us to identify novel antiviral targets related to apoptotic signaling pathways.

Spatiotemporal Distribution and Genetic Characterization of Measles Strains Circulating in Greece during the 2017-2018 Outbreak

Between May 2017 and November 2018, Greece has experienced a severe measles outbreak with a total of 3258 cases reported, after reaching its goal of eliminating measles since 2014-2015. In this study, we aimed to investigate the origin and the dispersal patterns of the measles strains that circulated in Greece during this outbreak and to identify possible transmission patterns of measles virus (MeV) in the country. Of the 832 measles suspect cases referred to the National Measles and Rubella Reference Laboratory for MeV RNA detection, 131 randomly selected positive samples, representative of the temporal and spatial distribution of the laboratory-confirmed measles cases in Greece, were processed for genotypic identification by an RT-PCR amplification of a 598 bp fragment containing the 450 bp hypervariable region of the measles virus N gene. Phylogenetic analysis was carried out by the approximate maximum likelihood method (ML) under the generalized time-reversible (GTR + cat) model. All samples analyzed were found to belong to genotype B3. Comparative analysis with other European and reference measles strains revealed three separate major clusters and other multiple viruses circulating simultaneously in Greece. They were all isolated from three main community groups, Greek-Roma children, non-minority Greek nationals and immigrants/refugees, a finding that is in accordance with what was also observed in the last two measles outbreaks in 2005-2006 and 2010-2011. Notably, for one of the three clusters, no similarity was detected with previously reported prototype strains. Our results indicate the need for a more intensive vaccination program against measles amongst minority populations and in refugee hot-spots as well as the importance of molecular surveillance as a tool for monitoring measles outbreaks.

SARS-CoV-2 exhibits intra-host genomic plasticity and low-frequency polymorphic quasispecies

In December 2019, an outbreak of atypical pneumonia (Coronavirus disease 2019 -COVID-19) associated with a novel coronavirus (SARS-CoV-2) was reported in Wuhan city, Hubei province, China. The outbreak was traced to a seafood wholesale market and human to human transmission was confirmed. The rapid spread and the death toll of the new epidemic warrants immediate intervention. The intra-host genomic variability of SARS-CoV-2 plays a pivotal role in the development of effective antiviral agents and vaccines, as well as in the design of accurate diagnostics. We analyzed NGS data derived from clinical samples of three Chinese patients infected with SARS-CoV-2, in order to identify small- and large-scale intra-host variations in the viral genome. We identified tens of low- or higher- frequency single nucleotide variations (SNVs) with variable density across the viral genome, affecting 7 out of 10 protein-coding viral genes. The majority of these SNVs (72/104) corresponded to missense changes. The annotation of the identified SNVs but also of all currently circulating strain variations revealed colocalization of intra-host as well as strain specific SNVs with primers and probes currently used in molecular diagnostics assays. Moreover, we de-novo assembled the viral genome, in order to isolate and validate intra-host structural variations and recombination breakpoints. The bioinformatics analysis disclosed genomic rearrangements over poly-A / poly-U regions located in ORF1ab and spike (S) gene, including a potential recombination hot-spot within S gene. Our results highlight the intra-host genomic diversity and plasticity of SARS-CoV-2, pointing out genomic regions that are prone to alterations. The isolated SNVs and genomic rearrangements reflect the intra-patient capacity of the polymorphic quasispecies, which may arise rapidly during the outbreak, allowing immunological escape of the virus, offering resistance to anti-viral drugs and affecting the sensitivity of the molecular diagnostics assays.

Risk of HIV transmission through condomless sex in serodifferent gay couples with the HIV-positive partner taking suppressive antiretroviral therapy (PARTNER): final results of a multicentre, prospective, observational study

BACKGROUND

The level of evidence for HIV transmission risk through condomless sex in serodifferent gay couples with the HIV-positive partner taking virally suppressive antiretroviral therapy (ART) is limited compared with the evidence available for transmission risk in heterosexual couples. The aim of the second phase of the PARTNER study (PARTNER2) was to provide precise estimates of transmission risk in gay serodifferent partnerships.

METHODS

The PARTNER study was a prospective observational study done at 75 sites in 14 European countries. The first phase of the study (PARTNER1; Sept 15, 2010, to May 31, 2014) recruited and followed up both heterosexual and gay serodifferent couples (HIV-positive partner taking suppressive ART) who reported condomless sex, whereas the PARTNER2 extension (to April 30, 2018) recruited and followed up gay couples only. At study visits, data collection included sexual behaviour questionnaires, HIV testing (HIV-negative partner), and HIV-1 viral load testing (HIV-positive partner). If a seroconversion occurred in the HIV-negative partner, anonymised phylogenetic analysis was done to compare HIV-1 pol and env sequences in both partners to identify linked transmissions. Couple-years of follow-up were eligible for inclusion if condomless sex was reported, use of pre-exposure prophylaxis or post-exposure prophylaxis was not reported by the HIV-negative partner, and the HIV-positive partner was virally suppressed (plasma HIV-1 RNA <200 copies per mL) at the most recent visit (within the past year). Incidence rate of HIV transmission was calculated as the number of phylogenetically linked HIV infections that occurred during eligible couple-years of follow-up divided by eligible couple-years of follow-up. Two-sided 95% CIs for the incidence rate of transmission were calculated using exact Poisson methods.

FINDINGS

Between Sept 15, 2010, and July 31, 2017, 972 gay couples were enrolled, of which 782 provided 1593 eligible couple-years of follow-up with a median follow-up of 2·0 years (IQR 1·1-3·5). At baseline, median age for HIV-positive partners was 40 years (IQR 33-46) and couples reported condomless sex for a median of 1·0 years (IQR 0·4-2·9). During eligible couple-years of follow-up, couples reported condomless anal sex a total of 76 088 times. 288 (37%) of 777 HIV-negative men reported condomless sex with other partners. 15 new HIV infections occurred during eligible couple-years of follow-up, but none were phylogenetically linked within-couple transmissions, resulting in an HIV transmission rate of zero (upper 95% CI 0·23 per 100 couple-years of follow-up).

INTERPRETATION

Our results provide a similar level of evidence on viral suppression and HIV transmission risk for gay men to that previously generated for heterosexual couples and suggest that the risk of HIV transmission in gay couples through condomless sex when HIV viral load is suppressed is effectively zero. Our findings support the message of the U=U (undetectable equals untransmittable) campaign, and the benefits of early testing and treatment for HIV.

FUNDING

National Institute for Health Research.

HIV-1 molecular transmission clusters in nine European countries and Canada: association with demographic and clinical factors

BACKGROUND

Knowledge of HIV-1 molecular transmission clusters (MTCs) is important, especially in large-scale datasets, for designing prevention programmes and public health intervention strategies. We used a large-scale HIV-1 sequence dataset from nine European HIV cohorts and one Canadian, to identify MTCs and investigate factors associated with the probability of belonging to MTCs.

METHODS

To identify MTCs, we applied maximum likelihood inferences on partial pol sequences from 8955 HIV-positive individuals linked to demographic and clinical data. MTCs were defined using two different criteria: clusters with bootstrap support >75% (phylogenetic confidence criterion) and clusters consisting of sequences from a specific region at a proportion of >75% (geographic criterion) compared to the total number of sequences within the network. Multivariable logistic regression analysis was used to assess factors associated with MTC clustering.

RESULTS

Although 3700 (41%) sequences belonged to MTCs, proportions differed substantially by country and subtype, ranging from 7% among UK subtype C sequences to 63% among German subtype B sequences. The probability of belonging to an MTC was independently less likely for women than men (OR = 0.66; P < 0.001), older individuals (OR = 0.79 per 10-year increase in age; P < 0.001) and people of non-white ethnicity (OR = 0.44; P < 0.001 and OR = 0.70; P = 0.002 for black and 'other' versus white, respectively). It was also more likely among men who have sex with men (MSM) than other risk groups (OR = 0.62; P < 0.001 and OR = 0.69; P = 0.002 for people who inject drugs, and sex between men and women, respectively), subtype B (ORs 0.36-0.70 for A, C, CRF01 and CRF02 versus B; all P < 0.05), having a well-estimated date of seroconversion (OR = 1.44; P < 0.001), a later calendar year of sampling (ORs 2.01-2.61 for all post-2002 periods versus pre-2002; all P < 0.01), and being naïve to antiretroviral therapy at sampling (OR = 1.19; P = 0.010).

CONCLUSIONS

A high proportion (>40%) of individuals belonged to MTCs. Notably, the HIV epidemic dispersal appears to be driven by subtype B viruses spread within MSM networks. Expansion of regional epidemics seems mainly associated with recent MTCs, rather than the growth of older, established ones. This information is important for designing prevention and public health intervention strategies.

Prevalence and Characteristics of Hepatitis B Virus (HBV) Coinfection among HIV-Positive Women in South Africa and Botswana

There is progressive concern about the evolving burden of morbidity and mortality caused by coinfection with HIV-1 and hepatitis B virus (HBV) in sub-Saharan Africa, but the epidemiology and impact of this problem are not well defined. We therefore set out to assimilate more information about the nature of HBV/HIV coinfection in this region by undertaking a retrospective observational study of southern African adult women. We used samples from previously recruited HIV-1 positive women attending antenatal clinics in three settings in South Africa and Botswana (n = 950) and added a small cohort of HIV-negative antenatal South African women for comparison (n = 72). We tested for HBsAg and followed up HBsAg-positive samples by testing for HBeAg, HBV DNA, HBV genotype, presence of drug-resistance associated mutations (RAMs) and HDV. We identified HBsAg in 72 individuals (7% of the whole cohort), of whom 27% were HBeAg-positive, and the majority HBV genotypes A1 and A2. We did not detect any HDV coinfection. HBV prevalence was significantly different between geographically distinct cohorts, but did not differ according to HIV status. Among adults from South Africa, HBV/HIV coinfected patients had lower CD4+ T cell counts compared to those with HIV-monoinfection (p = 0.02), but this finding was not replicated in the cohort from Botswana. Overall, these data provide a snapshot of the coinfection problem at the heart of the HIV/HBV co-epidemic, and are important to inform public health policy, resource allocation, education, surveillance and clinical care.

The genealogical population dynamics of HIV-1 in a large transmission chain: bridging within and among host evolutionary rates

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the ‘store and retrieve’ hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.

Since its discovery three decades ago, the HIV epidemic has unfolded into one of the most devastating pandemics in human history. When HIV replication cannot be completely inhibited, the fast-evolving retrovirus continuously evades intra-host immune and drug selective pressure, but diversifies according to more neutral epidemiological dynamics at the interhost level. Limited evidence suggests that the virus may evolve faster in a single host than in a population of hosts, and various hypotheses have been put forward to explain this phenomenon. Here, we develop a new computational approach aimed at integrating host transmission information with pathogen genealogical reconstructions. We apply this approach to comprehensive sequence data sets sampled from a large HIV-1 subtype C transmission chain, and in addition to providing several insights into the reconstruction of HIV-1 transmissions histories and its associated population dynamics, we find that transmission decreases the HIV-1 evolutionary rate. The fact that we also identify this decline for substitutions that do not alter amino acid substitutions provides evidence against hypotheses that invoke selection forces. Instead, our findings support earlier reports that new infections start preferentially with less evolved variants, which may be stored in latently infected cells, and this may vary among different HIV-1 subtypes.