Diversity and Long-Term Dynamics of Human Blood Anelloviruses.

Anelloviruses (AVs) are commensal members of the human blood virome. Even though it was estimated that over 90% of the human population carries AVs, the dynamics of the AV virome ("anellome") are unknown. We investigated the dynamics of blood anellomes in two healthy people followed up for more than 30 years. Both subjects were positive for AVs in the majority of samples. Alphatorquevirus (torque teno virus [TTV]) was the most common genus in both subjects, followed by Betatorquevirus (torque teno minivirus [TTMV]) and Gammatorquevirus (torque teno midivirus [TTMDV]). Almost five times more lineages were found in subject 1 than in subject 2, and the anellomes differed phylogenetically. Both anellomes remained compositionally stable, and 9 out of 64 AV lineages were detected in over half of the time points. We confirmed the long-term and short-term persistence of 13 lineages by specific quantitative PCR (qPCR). AV lineages were detected in blood for over 30 years. Noticeable differences in anellome richness were found between the tested subjects, but both anellomes remained compositionally stable over time. These findings demonstrate that the human blood anellome is personal and that AV infection is chronic and potentially commensal. IMPORTANCE Knowledge of the persistence of AVs in humans is crucial to our understanding of the nature of AV infection (chronic or acute) and the role of AV in the host. We therefore investigated the dynamics of anellovirus infection in two healthy people followed up for 30 years. Our findings suggest that the human blood anellovirus virome (anellome) remains stable and personal for decades.

Similar Risk of Severe Acute Respiratory Syndrome Coronavirus 2 Infection and Similar Nucleocapsid Antibody Levels in People With Well-Controlled Human Immunodeficiency Virus (HIV) and a Comparable Cohort of People Without HIV.

BACKGROUND: Within the ongoing AGEhIV Cohort Study in Amsterdam, we prospectively compared the incidence of and risk factors for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection between human immunodeficiency virus (HIV)-positive and HIV-negative participants. Moreover, we compared SARS-CoV-2 nucleocapsid antibody levels between participants with incident infection from both groups. METHODS: Starting in September 2020, consenting HIV-positive and HIV-negative participants were assessed every 6 months for incident SARS-CoV-2 infection, using combined immunoglobulin (Ig) A/IgM/IgG SARS-CoV-2 nucleocapsid antibody assay. Cumulative incidence of SARS-CoV-2 infection and associated risk factors were assessed from 27 February 2020 through 30 April 2021, using complementary log-log regression. In those with incident SARS-CoV-2 infection, nucleocapsid (N) antibody levels were compared between groups using linear regression. RESULTS: The study included 241 HIV-positive (99.2% virally suppressed) and 326 HIV-negative AGEhIV participants. The cumulative SARS-CoV-2 incidence by April 2021 was 13.4% and 11.6% in HIV-positive and HIV-negative participants, respectively (P = .61). Younger age and African origin were independently associated with incident infection. In those with incident infection, only self-reported fever, but not HIV status, was associated with higher N antibody levels. CONCLUSIONS: HIV-positive individuals with suppressed viremia and adequate CD4 cell counts had similar risk of SARS-CoV-2 acquisition and similar SARS-CoV-2 N antibody levels after infection compared with a comparable HIV-negative cohort. CLINICAL TRIAL REGISTRATION: NCT01466582.

Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe.

HIV-1 set-point viral load-the approximately stable value of viraemia in the first years of chronic infection-is a strong predictor of clinical outcome and is highly variable across infected individuals. To better understand HIV-1 pathogenesis and the evolution of the viral population, we must quantify the heritability of set-point viral load, which is the fraction of variation in this phenotype attributable to viral genetic variation. However, current estimates of heritability vary widely, from 6% to 59%. Here we used a dataset of 2,028 seroconverters infected between 1985 and 2013 from 5 European countries (Belgium, Switzerland, France, the Netherlands and the United Kingdom) and estimated the heritability of set-point viral load at 31% (CI 15%-43%). Specifically, heritability was measured using models of character evolution describing how viral load evolves on the phylogeny of whole-genome viral sequences. In contrast to previous studies, (i) we measured viral loads using standardized assays on a sample collected in a strict time window of 6 to 24 months after infection, from which the viral genome was also sequenced; (ii) we compared 2 models of character evolution, the classical "Brownian motion" model and another model ("Ornstein-Uhlenbeck") that includes stabilising selection on viral load; (iii) we controlled for covariates, including age and sex, which may inflate estimates of heritability; and (iv) we developed a goodness of fit test based on the correlation of viral loads in cherries of the phylogenetic tree, showing that both models of character evolution fit the data well. An overall heritability of 31% (CI 15%-43%) is consistent with other studies based on regression of viral load in donor-recipient pairs. Thus, about a third of variation in HIV-1 virulence is attributable to viral genetic variation.

Correction: Viral genetic variation accounts for a third of variability in HIV-1 set-point viral load in Europe.

[This corrects the article DOI: 10.1371/journal.pbio.2001855.].

Widespread hepatitis B virus genotype G (HBV-G) infection during the early years of the HIV epidemic in the Netherlands among men who have sex with men.

BACKGROUND: Hepatitis B virus (HBV) variants belong to different genotypes, A-J, whose worldwide distribution is linked with geography, probably because viral spread was associated with ancient human migrations. HBV genotype G (HBV-G) is an aberrant genotype with little sequence divergence, suggesting a recent origin. HBV-G is strongly associated with certain risk groups such as intravenous drug users (IDUs) and men who have sex with men (MSM), but hardly with geography. The origin and epidemiology of HBV-G remain unresolved, as is the disease association. METHODS: To estimate the prevalence and possible time of introduction of HBV-G into the MSM community in Amsterdam, the Netherlands, we have retrospectively analysed 226 blood serum samples from HBsAg positive MSM enrolled in the Amsterdam Cohort Studies (ACS) on HIV infection and AIDS dating from 1984 to 1999 using genotype-specific PCR assays. RESULTS: Of the 226 HBsAg-positive samples, 149 were HBV DNA positive. Of those, 104 were positive for HBV genotype A (HBV-A) and five for HBV-G, and 40 showed a dual infection with both HBV-A and HBV-G. Being HIV-infected was significantly associated with a reduced HBV DNA viral load in blood, but not with the prevalence of HBV-G. Early virus already contained stop codons in the precore region and a 36 bp insertion in the core gene which are the characteristics of HBV-G. CONCLUSIONS: HBV-G was introduced before 1985 into the Amsterdam MSM community. Early isolates show very limited sequence variation, confirming a low evolutionary rate. HBV-G acquisition was independent of HIV infection, but being HIV-infected was significantly associated with a reduced HBV viral load in blood, indicating a beneficial effect of early HIV infection in controlling HBV replication.

A novel genus in the order Picornavirales detected in human stool.

Recently, four new viruses belonging to an unassigned family within the order Picornavirales were identified in excrements of healthy carp (fisavirus) and pigs (posavirus 1, 2 and 3). We report the detection and characterization of a fifth virus present in human faeces. The virus, named human stool-associated RNA virus (husavirus), contains a single ORF encoding a putative 2993 AA polyprotein, with a Hel-Pro-Pol replication block, typical for the Picornavirales. Phylogenetic analysis revealed that the closest relative to husavirus is posavirus 1, and together they cluster with fisavirus, posavirus 2 and 3 and a roundworm (Ascaris suum) derived virus. Husavirus was detected in eight human stool samples collected in 1984 (n52), 1985 (n54), 1995 (n51) and 2014 (n51). From three strains of husavirus from 1984 and 1985 the full genome sequence was determined, showing less than 5% intraspecies variation in the nucleotide composition. The host of this virus remains to be determined.

Influenza-like illness symptoms due to endemic human coronavirus reinfections are not influenced by the length of the interval separating reinfections.

After 3 years of its introduction to humans, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been declared as endemic. Little is known about the severity of the disease manifestation that future infections may cause, especially when reinfections occur after humoral immunity from a previous infection or vaccination has waned. Such knowledge could inform policymakers regarding the frequency of vaccination. Reinfections by endemic human coronaviruses (HCoVs) can serve as a model system for SARS-CoV-2 endemicity. We monitored 44 immunocompetent male adults with blood sampling every 6 months (for 17 years), for the frequency of HCoV (re-)infections, using rises in N-antibodies of HCoV-NL63, HCoV-29E, HCoV-OC43, and HCoV-HKU1 as markers of infection. Disease associations during (re-)infections were examined by comparison of self-reporting records of influenza-like illness (ILI) symptoms, every 6 months, by all participants. During 8,549 follow-up months, we found 364 infections by any HCoV with a median of eight infections per person. Symptoms more frequently reported during HCoV infection were cough, sore throat, and myalgia. Two hundred fifty-one of the 364 infections were species-specific HCoV-reinfections, with a median interval of 3.58 (interquartile range 1.92-5.67) years. The length of the interval between reinfections-being either short or long-had no influence on the frequency of reporting ILI symptoms. All HCoV-NL63, HCoV-229E, HCoV-OC43, and HCoV-HKU1 (re-)infections are associated with the reporting of ILIs. Importantly, in immunocompetent males, these symptoms are not influenced by the length of the interval between reinfections. IMPORTANCE: Little is known about the disease following human coronavirus (HCoV) reinfection occurring years after the previous infection, once humoral immunity has waned. We monitored endemic HCoV reinfection in immunocompetent male adults for up to 17 years. We found no influence of reinfection interval length in the disease manifestation, suggesting that immunocompetent male adults are adequately protected against future HCoV infections.

Determinants of epidemic size and the impacts of lulls in seasonal influenza virus circulation.

During the COVID-19 pandemic, levels of seasonal influenza virus circulation were unprecedentedly low, leading to concerns that a lack of exposure to influenza viruses, combined with waning antibody titres, could result in larger and/or more severe post-pandemic seasonal influenza epidemics. However, in most countries the first post-pandemic influenza season was not unusually large and/or severe. Here, based on an analysis of historical influenza virus epidemic patterns from 2002 to 2019, we show that historic lulls in influenza virus circulation had relatively minor impacts on subsequent epidemic size and that epidemic size was more substantially impacted by season-specific effects unrelated to the magnitude of circulation in prior seasons. From measurements of antibody levels from serum samples collected each year from 2017 to 2021, we show that the rate of waning of antibody titres against influenza virus during the pandemic was smaller than assumed in predictive models. Taken together, these results partially explain why the re-emergence of seasonal influenza virus epidemics was less dramatic than anticipated and suggest that influenza virus epidemic dynamics are not currently amenable to multi-season prediction.

Translational HIV-1 research: from routine diagnostics to new virology insights in Amsterdam, the Netherlands during 1983-2013.

An HIV-1 diagnostic laboratory was established in the Academic Medical Center (AMC) of the University of Amsterdam after the discovery of human immunodeficiency virus (HIV) as the cause of the acquired immunodeficiency syndrome (AIDS). The first AIDS patients were diagnosed here in 1981 and since 1983 we have tested the samples of 50992 patients using a variety of assays that greatly improved over the years. We will describe some of the basic results from this diagnostic laboratory and then focus on the spin-off in terms of the development of novel virus assays to detect super-infections and ultra-sensitive assays to measure the intracellular HIV-1 RNA load. We also review several original research findings in the field of HIV-1 virology that stem from initial observations made in the diagnostic unit. This includes the study of genetic defects in the HIV-1 genome and time trends of the replication fitness over 30 years of viral evolution, but also the description of novel HIV-1 variants in difficult-to-diagnose clinical specimen.

Human immunodeficiency virus type 1 gp120 envelope characteristics associated with disease progression differ in family members infected with genetically similar viruses.

The human immunodeficiency virus type 1 (HIV-1) envelope protein provides the primary contact between the virus and host, and is the main target of the adaptive humoral immune response. The length of gp120 variable loops and the number of N-linked glycosylation events are key determinants for virus infectivity and immune escape, while the V3 loop overall positive charge is known to affect co-receptor tropism. We selected two families in which both parents and two children had been infected with HIV-1 for nearly 10 years, but who demonstrated variable parameters of disease progression. We analysed the gp120 envelope sequence and compared individuals that progressed to those that did not in order to decipher evolutionary alterations that are associated with disease progression when individuals are infected with genetically related virus strains. The analysis of the V3-positive charge demonstrated an association between higher V3-positive charges with disease progression. The ratio between the amino acid length and the number of potential N-linked glycosylation sites was also shown to be associated with disease progression with the healthier family members having a lower ratio. In conclusion in individuals initially infected with genetically linked virus strains the V3-positive charges and N-linked glycosylation are associated with HIV-1 disease progression and follow varied evolutionary paths for individuals with varied disease progression.

Identification of a new genotype of Torque Teno Mini virus.

BACKGROUND: Although human torque teno viruses (TTVs) were first discovered in 1997, still many associated aspects are not clarified yet. The viruses reveal a remarkable heterogeneity and it is possible that some genotypes are more pathogenic than others. The identification of all genotypes is essential to confirm previous pathogenicity data, and an unbiased search for novel viruses is needed to identify TTVs that might be related to disease. METHOD: The virus discovery technique VIDISCA-454 was used to screen serum of 55 HIV-1 positive injecting drug users, from the Amsterdam Cohort Studies, in search for novel blood-blood transmittable viruses which are undetectable via normal diagnostics or panvirus-primer PCRs. RESULTS: A novel torque teno mini virus (TTMV) was identified in two patients and the sequence of the full genomes were determined. The virus is significantly different from the known TTMVs (< 40% amino acid identity in ORF1), yet it contains conserved characteristics that are also present in other TTMVs. The virus is chronically present in both patients, and these patients both suffered from a pneumococcal pneumonia during follow up and had extremely low B-cells counts. CONCLUSION: We describe a novel TTMV which we tentatively named TTMV-13. Further research is needed to address the epidemiology and pathogenicity of this novel virus.

High prevalence of hepatitis B virus dual infection with genotypes A and G in HIV-1 infected men in Amsterdam, the Netherlands, during 2000-2011.

BACKGROUND: Hepatitis B virus (HBV) is divided into 8 definite (A-H) and 2 putative (I, J) genotypes that show a geographical distribution. HBV genotype G, however, is an aberrant genotype of unknown origin that demonstrates severe replication deficiencies and very little genetic variation. It is often found in co-infections with another HBV genotype and infection has been associated with certain risk groups such as intravenous drug users and men having sex with men (MSM). We aimed to estimate the prevalence of HBV-G in the Netherlands by analysing samples from HBV-positive patients visiting the Academic Medical Center in Amsterdam. METHODS: Ninety-six HBV-infected patients, genotyped as HBV-A or HBV-G infected, were retrieved from the clinical database. Blood plasma samples were analysed with a newly-developed real-time PCR assay that detects HBV-A and HBV-G. For three patients, the HBV plasma viral load (pVL) of both genotypes was followed longitudinally. In addition, three complete genomes of HBV-G were sequenced to determine their relationship to global HBV-G strains. RESULTS: Ten HBV-G infections were found in the selected Dutch patients. All concerned HIV-1 infected males with HBV-A co-infection. Dutch HBV-G strains were phylogenetically closely related to reference HBV-G strains. CONCLUSIONS: In this study, HBV-G infection in the Netherlands is found exclusively in HIV-1 infected men as co-infection with HBV-A. A considerable percentage (37%) of men infected with HBV and HIV-1 are actually co- infected with two HBV genotypes.

Modest nonadherence to antiretroviral therapy promotes residual HIV-1 replication in the absence of virological rebound in plasma.

BACKGROUND: Modern antiretroviral therapy (ART) regimens are widely assumed to forgive modest nonadherence, because virological suppression in plasma is common at adherence levels of >70%. Yet, it is unknown whether human immunodeficiency virus type 1 (HIV-1) replication is completely suppressed at these levels of adherence. METHODS: We longitudinally quantified levels of cell-associated HIV-1 RNA and DNA in 40 patients (median duration of successful ART before study initiation, 46 months), whose 1-week adherence to therapy prior to the sampling moments was measured electronically. RESULTS: Patients were constantly 100% adherent (the optimal-adherence group), demonstrated improving adherence over time (the improving-adherence group), or neither of the above (the poor-adherence group). Adherence never decreased to <70% in any patient, and no rebound in plasma virological levels was observed. Nevertheless, poor adherence but not optimal or improving adherence caused a significant longitudinal increase in cell-associated HIV RNA levels (P = .006). Time-weighted changes and regression slopes of viral RNA load for the poor-adherence group were significantly higher than those for the optimal-adherence group (P < .01). CONCLUSIONS: Because ART only blocks infection of new cells but not viral RNA transcription in cells infected before therapy initiation, the observed effects strongly suggest that modest nonadherence can cause new cycles of HIV-1 replication that are undetectable by commercial plasma viral load assays.

Anellovirus evolution during long-term chronic infection.

Human anelloviruses (AVs) are extremely genetically diverse, are widespread in the human population, and cause chronic infections. However, the evolutionary dynamics of AVs within single hosts is currently unknown, and it is unclear whether these changes have an implication on the long-term persistence of AVs in the host. Here, we assessed the evolutionary dynamics of six AV lineages during 30 years of chronic infection at single host resolution. The total number of substitutions and the number of variable sites increased over time. However, not all substitutions reached population fixation, showing that AV lineages form heterogeneous swarms within the host. Most substitutions occurred within a hypervariable region (HVR) located between nucleotide positions 800 and 1,300 of ORF1, which is known to be located within the spike domain. Different regions of the ORF1 gene undergo either positive or negative selection pressure. Sites under strong diversifying selection pressure were detected in the HVR, while the majority of the sites under purifying selection were detected outside this region. The HVR may play the role of an immunological decoy that prevents antibodies from binding to more vulnerable parts of ORF1. Moreover, the frequent substitutions in this region may increase the chances of AV particles escaping immune recognition.

Robust Vaccine-Induced as Well as Hybrid B- and T-Cell Immunity across SARS-CoV-2 Vaccine Platforms in People with HIV.

Few studies have comprehensively compared severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine-induced and hybrid B- and T-cell responses in people with HIV (PWH) to those in comparable controls without HIV. We included 195 PWH and 246 comparable controls from the AGEhIV COVID-19 substudy. A positive nucleocapsid antibody (INgezim IgA/IgM/IgG) or self-reported PCR test defined prior SARS-CoV-2 infection. SARS-CoV-2 anti-spike (anti-S) IgG titers and anti-S IgG production by memory B cells were assessed. Neutralizing antibody titers were determined in a subset of participants. T-cell responses were assessed by gamma interferon (IFN-γ) release and activation-induced marker assay. We estimated mean differences in postvaccination immune responses (ß) between levels of determinants. Anti-S IgG titers and anti-S IgG production by memory B cells were not different between PWH and controls. Prior SARS-CoV-2 infection (ß = 0.77), receiving mRNA vaccine (ß = 0.56), female sex (ß = 0.24), fewer days between last vaccination and sampling (ß = 0.07), and a CD4/CD8 ratio of <1.0 (ß = -0.39) were independently associated with anti-S IgG titers, but HIV status was not. Neutralization titers against the ancestral and Delta and Omicron SARS-CoV-2 variants were not different between PWH and controls. IFN-γ release was higher in PWH. Prior SARS-CoV-2 infection (ß = 2.39), HIV-positive status (ß = 1.61), and fewer days between last vaccination and sampling (ß = 0.23) were independently associated with higher IFN-γ release. The percentages of SARS-CoV-2-reactive CD4+ and CD8+ T cells, however, were not different between PWH and controls. Individuals with well-controlled HIV generally mount robust vaccine-induced as well as hybrid B- and T-cell immunity across SARS-CoV-2 vaccine platforms similar to controls. Determinants of a reduced vaccine response were likewise largely similar in both groups and included a lower CD4/CD8 ratio. IMPORTANCE Some studies have suggested that people with HIV may respond less well to vaccines against SARS-CoV-2. We comprehensively compared B- and T-cell responses to different COVID-19 vaccines in middle-aged persons with well-treated HIV and individuals of the same age without HIV, who were also highly comparable in terms of demographics and lifestyle, including those with prior SARS-CoV-2 infection. Individuals with HIV generally mounted equally robust immunity to the different vaccines. Even stronger immunity was observed in both groups after prior SARS-CoV-2 infection. These findings are reassuring with respect to the efficacy of SARS-Cov-2 vaccines for the sizable and increasing global population of people with HIV with access and a good response to HIV treatment.

HIV-1 dual infection is associated with faster CD4+ T-cell decline in a cohort of men with primary HIV infection.

BACKGROUND: In vitro, animal, and mathematical models suggest that human immunodeficiency virus (HIV) co- or superinfection would result in increased fitness of the pathogen and, possibly, increased virulence. However, in patients, the impact of dual HIV type 1 (HIV-1) infection on disease progression is unclear, because parameters relevant for disease progression have not been strictly analyzed. The objective of the present study is to analyze the effect of dual HIV-1 infections on disease progression in a well-defined cohort of men who have sex with men. METHODS: Between 2000 and 2009, 37 men who had primary infection with HIV-1 subtype B, no indication for immediate need of combination antiretroviral therapy (cART), and sufficient follow-up were characterized with regard to dual infection or single infection and to coreceptor use. Patients were followed to estimate the effect of these parameters on clinical disease progression, as defined by the rate of CD4(+) T-cell decline and the time to initiation of cART. RESULTS: Four patients presented with HIV-1 coinfection; 6 patients acquired HIV-1 superinfection, on average 8.5 months from their primary infection; and 27 patients remained infected with a single strain. Slopes of longitudinal CD4(+) T-cell counts and time-weighted changes from baseline were significantly steeper for patients with dual infection compared with patients with single infection. Multivariate analysis showed that the most important parameter associated with CD4(+) T-cell decline over time was dual infection (P = .001). Additionally, patients with HIV-1 coinfection had a significantly earlier start of cART (P < .0001). CONCLUSIONS: Dual HIV-1 infection is the main factor associated with CD4(+) T-cell decline in men who have untreated primary infection with HIV-1 subtype B.

Control of Human Anelloviruses by Cytosine to Uracil Genome Editing.

Anelloviruses are the most common viruses infecting humans. Every human carries a nonpathogenic personal anellovirus virome (anellome), yet it is unknown which mechanisms contribute to its stability. Here, we assessed the dynamics and impact of a host antiviral defense mechanism-cytidine deaminase activity leading to C to U editing in anelloviruses-on the stability of the anellome. We investigated anellome sequence data obtained from serum samples collected every 6 months from two healthy subjects followed for more than 30 years. The subjects were infected by a total of 64 anellovirus lineages. Minus-stranded C to U editing was observed in lineages belonging to the Alpha-, Beta-, and Gammatorquevirus genera. The edited genomes were present within virus particles, therefore editing must have occurred at the late stages of the virus life cycle. Editing was favored by 5'-TC contexts in the virus genome, indicating that apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like, catalytic subunit 3 or A3 (APOBEC3) proteins are involved. Within a lineage, mutational dynamics varied over time and few fixations of mutations were detected, indicating that C to U editing is a dead end for a virus genome. We detected an editing coldspot in the GC-rich regions, suggesting that the GC-rich region is crucial for genome packaging, since only packaged virus particles were included in the analysis. Finally, we noticed a lineage-specific reduced concentration after an editing event, yet no clearance. In conclusion, cytidine deaminase activity does not clear anelloviruses, nor does it play a major role in virus evolution, but it does contribute to the stability of the anellome. IMPORTANCE Despite significant attention on anellovirus research, the interaction between the anellovirus virome and the human host remains unknown. We show the dynamics of APOBEC3-mediated cytidine deaminase activity on anelloviruses during a 30-year period of chronic infection and postulate that this antiviral mechanism controls anelloviruses. These results expand our knowledge of anellovirus-host interactions, which may be important for the design of gene therapies.

Transmission of anelloviruses to HIV-1 infected children.

Anelloviruses (AVs) are widespread in the population and infect humans at the early stage of life. The mode of transmission of AVs is still unknown, however, mother-to-child transmission, e.g., via breastfeeding, is one of the likely infection routes. To determine whether the mother-to-child transmission of AVs may still occur despite the absence of natural birth and breastfeeding, 29 serum samples from five HIV-1-positive mother and child pairs were Illumina-sequenced. The Illumina reads were mapped to an AV lineage database "Anellometrix" containing 502 distinct ORF1 sequences. Although the majority of lineages from the mother were not shared with the child, the mother and child anellomes did display a significant similarity. These findings suggest that AVs may be transmitted from mothers to their children via different routes than delivery or breastfeeding.

Early-Life Colonization by Anelloviruses in Infants.

Anelloviruses (AVs) are found in the vast majority of the human population and are most probably part of a healthy virome. These viruses infect humans in the early stage of life, however, the characteristics of the first colonizing AVs are still unknown. We screened a collection of 107 blood samples from children between 0.4 and 64.8 months of age for the presence of three AV genera: the Alpha-, Beta- and Gammatorquevirus. The youngest child that was positive for AV was 1.2 months old, and a peak in prevalence (100% of samples positive) was reached between the twelfth and eighteenth months of life. Intriguingly, the beta- and gammatorqueviruses were detected most at the early stage of life (up to 12 months), whereas alphatorqueviruses, the most common AVs in adults, increased in prevalence in children older than 12 months. To determine whether that order of colonization may be related to oral transmission and unequal presence of AV genera in breast milk, we examined 63 breast milk samples. Thirty-two percent of the breast milk samples were positive in a qPCR detecting beta- and gammatorqueviruses, while alphatorqueviruses were detected in 10% of the samples, and this difference was significant (p = 0.00654). In conclusion, we show that beta- and gammatorqueviruses colonize humans in the first months of life and that breastfeeding could play a role in AV transmission.

Phylogenetic estimation of the viral fitness landscape of HIV-1 set-point viral load.

Set-point viral load (SPVL), a common measure of human immunodeficiency virus (HIV)-1 virulence, is partially determined by viral genotype. Epidemiological evidence suggests that this viral property has been under stabilising selection, with a typical optimum for the virus between 104 and 105 copies of viral RNA per ml. Here we aimed to detect transmission fitness differences between viruses from individuals with different SPVLs directly from phylogenetic trees inferred from whole-genome sequences. We used the local branching index (LBI) as a proxy for transmission fitness. We found that LBI is more sensitive to differences in infectiousness than to differences in the duration of the infectious state. By analysing subtype-B samples from the Bridging the Evolution and Epidemiology of HIV in Europe project, we inferred a significant positive relationship between SPVL and LBI up to approximately 105 copies/ml, with some evidence for a peak around this value of SPVL. This is evidence of selection against low values of SPVL in HIV-1 subtype-B strains, likely related to lower infectiousness, and perhaps a peak in the transmission fitness in the expected range of SPVL. The less prominent signatures of selection against higher SPVL could be explained by an inherent limit of the method or the deployment of antiretroviral therapy.