SARS-CoV-2 variants with reduced infectivity and varied sensitivity to the BNT162b2 vaccine are developed during the course of infection

In-depth analysis of SARS-CoV-2 quasispecies is pivotal for a thorough understating of its evolution during infection. The recent deployment of COVID-19 vaccines, which elicit protective anti-spike neutralizing antibodies, has stressed the importance of uncovering and characterizing SARS-CoV-2 variants with mutated spike proteins. Sequencing databases have allowed to follow the spread of SARS-CoV-2 variants that are circulating in the human population, and several experimental platforms were developed to study these variants. However, less is known about the SARS-CoV-2 variants that are developed in the respiratory system of the infected individual. To gain further insight on SARS-CoV-2 mutagenesis during natural infection, we preformed single-genome sequencing of SARS-CoV-2 isolated from nose-throat swabs of infected individuals. Interestingly, intra-host SARS-CoV-2 variants with mutated S genes or N genes were detected in all individuals who were analyzed. These intra-host variants were present in low frequencies in the swab samples and were rarely documented in current sequencing databases. Further examination of representative spike variants identified by our analysis showed that these variants have impaired infectivity capacity and that the mutated variants showed varied sensitivity to neutralization by convalescent plasma and to plasma from vaccinated individuals. Notably, analysis of the plasma neutralization activity against these variants showed that the L1197I mutation at the S2 subunit of the spike can affect the plasma neutralization activity. Together, these results suggest that SARS-CoV-2 intra-host variants should be further analyzed for a more thorough characterization of potential circulating variants.

Simplified, Shear Induced Generation of Double Emulsions for Robust Compartmentalization during Single Genome Analysis

Drop microfluidics has driven innovations for high throughput, low input analysis techniques such as single-cell RNA-seq. However, the instability of single emulsion (SE) drops occasionally causes significant merging during drop processing, limiting most applications to single-step reactions in drops. Here, we show that double emulsion (DE) drops address this critical limitation and completely prevent content mixing, which is essential for single entity analysis. DEs show excellent stability during thermal cycling. More importantly, DEs undergo rupture into the continuous phase instead of merging, preventing content mixing and eliminating unstable drops from the downstream analysis. Due to the lack of drop merging, the monodispersity of drops is maintained throughout a workflow, enabling the deterministic manipulation of drops downstream. We also developed a simple, one-layer fabrication method for DE drop makers. This design is powerful as it allows robust production of single-core DEs at a wide range of flow rates and better control over the shell thickness, both of which have been significant limitations of conventional two-layer devices. This approach makes the fabrication of DE devices much more accessible, facilitating its broader adoption. Finally, we show that DE droplets effectively maintain the compartmentalization of single virus genomes during PCR-based amplification and barcoding, while SEs mixed contents due to merging. With their resistance to content mixing, DE drops have key advantages for multistep reactions in drops, which is limited in SEs due to merging and content mixing.

Increased Frequency of Inter-Subtype HIV-1 Recombinants Identified by Near Full-Length Virus Sequencing in Rwandan Acute Transmission Cohorts

Most studies of HIV-1 transmission have focused on subtypes B and C. In this study, we determined the genomic sequences of the transmitted founder (TF) viruses from acutely infected individuals enrolled between 2005 and 2011 into IAVI protocol C in Rwanda and have compared these isolates to viruses from more recent (2016–2019) acute/early infections in three at risk populations – MSM, high risk women (HRW), and discordant couples (DC). For the Protocol C samples, we utilized near full-length single genome (NFLG) amplification to generate 288 HIV-1 amplicons from 26 acutely infected seroconverters (SC), while for the 21 recent seroconverter samples (13 from HRW, two from DC, and six from MSM), we PCR amplified overlapping half-genomes. Using PacBio SMRT technology combined with the MDPseq workflow, we performed multiplex sequencing to obtain high accuracy sequences for each amplicon. Phylogenetic analyses indicated that the majority of recent transmitted viruses from DC and HRW clustered within those of the earlier Protocol C cohort. However, five of six sequences from the MSM cohort branched together and were greater than 97% identical. Recombination analyses revealed a high frequency (6/26; 23%) of unique inter-subtype recombination in Protocol C with 19% AC and 4% CD recombinant viruses, which contrasted with only 6.5% of recombinants defined by sequencing of the pol gene previously. The frequency of recombinants was significantly higher (12/21; 57%) in the more recent isolates, although, the five related viruses from the MSM cohort had identical recombination break points. While major drug resistance mutations were absent from Protocol C viruses, 4/21 of recent isolates exhibited transmitted nevirapine resistance. These results demonstrate the ongoing evolution and increased prevalence of recombinant and drug resistant transmitted viruses in Rwanda and highlight the importance of defining NFLG sequences to fully understand the nature of TF viruses and in particular the prevalence of unique recombinant forms (URFs) in transmission cohorts.

NanoHIV: A Bioinformatics Pipeline for Producing Accurate, Near Full-Length HIV Proviral Genomes Sequenced Using the Oxford Nanopore Technology

HIV-1 proviral single-genome sequencing by limiting-dilution polymerase chain reaction (PCR) amplification is important for differentiating the sequence-intact from defective proviruses that persist during antiretroviral therapy (ART). Intact proviruses may rebound if ART is interrupted and are the barrier to an HIV cure. Oxford Nanopore Technologies (ONT) sequencing offers a promising, cost-effective approach to the sequencing of long amplicons such as near full-length HIV-1 proviruses, but the high diversity of HIV-1 and the ONT sequencing error render analysis of the generated data difficult. NanoHIV is a new tool that uses an iterative consensus generation approach to construct accurate, near full-length HIV-1 proviral single-genome sequences from ONT data. To validate the approach, single-genome sequences generated using NanoHIV consensus building were compared to Illumina® consensus building of the same nine single-genome near full-length amplicons and an average agreement of 99.4% was found between the two sequencing approaches.

Characterizing the Relationship Between Neutralization Sensitivity and env Gene Diversity During ART Suppression

Although antiretroviral therapy (ART) successfully suppresses HIV-1 replication, ART-treated individuals must maintain therapy to avoid rebound from an integrated viral reservoir. Strategies to limit or clear this reservoir are urgently needed. Individuals infected for longer periods prior to ART appear to harbor more genetically diverse virus, but the roles of duration of infection and viral diversity in the humoral immune response remain to be studied. We aim to clarify a role, if any, for autologous and heterologous antibodies in multi-pronged approaches to clearing infection. To that end, we have characterized the breadths and potencies of antibody responses in individuals with varying durations of infection and HIV-1 envelope (env) gene diversity as well as the sensitivity of their inducible virus reservoir to broadly neutralizing antibodies (bNAbs). Plasma was collected from 8 well-characterized HIV-1+ males on ART with varied durations of active infection. HIV envs from reservoir-derived outgrowth viruses were amplified and single genome sequenced in order to measure genetic diversity in each participant. IgG from plasma was analyzed for binding titers against gp41 and gp120 proteins, and for neutralizing titers against a global HIV-1 reference panel as well as autologous outgrowth viruses. The sensitivity to bNAbs of these same autologous viruses was measured. Overall, we observed that greater env diversity was associated with higher neutralizing titers against the global panel and also increased resistance to certain bNAbs. Despite the presence of robust anti-HIV-1 antibody titers, we did not observe potent neutralization against autologous viruses. In fact, 3 of 8 participants harbored viruses that were completely resistant to the highest tested concentration of autologous IgG. That this lack of neutralization was observed regardless of ART duration or viral diversity suggests that the inducible reservoir harbors ‘escaped’ viruses (that co-evolved with autologous antibody responses), rather than proviruses archived from earlier in infection. Finally, we observed that viruses resistant to autologous neutralization remained sensitive to bNAbs, especially CD4bs and MPER bNAbs. Overall, our data suggest that the inducible reservoir is relatively resistant to autologous antibodies and that individuals with limited virus variation in the env gene, such as those who start ART early in infection, are more likely to be sensitive to bNAb treatment.

Patterns of selection in the evolution of a transposable element

Transposable elements (TEs) are a major component of most eukaryotic genomes. Here, we present a new approach which allows us to study patterns of natural selection in the evolution of TEs over short time scales. The method uses the alignment of all elements with intact gag/pol genes of a TE family from a single genome. We predict that the ratio of non-synonymous to synonymous variants (vN/vS) in the alignment should decrease as a function of the frequency of the variants, because elements with non-synonymous variants that reduce transposition will have fewer progeny. We apply our method to Sirevirus LTR retrotransposons that are abundant in maize and other plant species and show that vN/vS declines as variant frequency increases, indicating that negative selection is acting strongly on the Sirevirus genome. The asymptotic value of vN/vS suggests that at least 85% of all non-synonymous mutations in the TE reduce transposition. Crucially, these patterns in vN/vS are only predicted to occur if the gene products from a particular TE insertion preferentially promote the transposition of the same insertion. Overall, this study is the first to use large numbers of intact elements to shed new light on the selective processes that act on TEs.

Detection of F1 hybrids from single-genome data reveals frequent hybridization in Hymenoptera and particularly ants

Hybridization occupies a central role in many fundamental evolutionary processes, such as speciation or adaptation. Yet, despite its pivotal importance in evolution, little is known about the actual prevalence and distribution of hybridization across the tree of life. Here we develop and implement a new statistical method enabling the detection of F1 hybrids from single-individual genome sequencing data. Using simulations and sequencing data from known hybrid systems, we first demonstrate the specificity of the method, and identify its statistical limits. Next, we showcase the method by applying it to available sequencing data from more than 1500 species of Arthropods, including Hymenoptera, Hemiptera, Coleoptera, Diptera and Archnida. Among these taxa, we find Hymenoptera, and especially ants, to display the highest number of candidate F1 hybrids, suggesting higher rates of recent hybridization in these groups. The prevalence of F1 hybrids was heterogeneously distributed across ants, with taxa including many candidates tending to harbor specific ecological and life history traits. This work shows how large-scale genomic comparative studies of recent hybridization can be implemented, uncovering the determinants of hybridization frequency across whole taxa.