Antigen pressure from two founder viruses induces multiple insertions at a single antibody position to generate broadly neutralizing HIV antibodies.

Vaccination strategies aimed at maturing broadly neutralizing antibodies (bnAbs) from naïve precursors are hindered by unusual features that characterize these Abs, including insertions and deletions (indels). Longitudinal studies of natural HIV infection cases shed light on the complex processes underlying bnAb development and have suggested a role for superinfection as a potential enhancer of neutralization breadth. Here we describe the development of a potent bnAb lineage that was elicited by two founder viruses to inform vaccine design. The V3-glycan targeting bnAb lineage (PC39-1) was isolated from subtype C-infected IAVI Protocol C elite neutralizer, donor PC39, and is defined by the presence of multiple independent insertions in CDRH1 that range from 1-11 amino acids in length. Memory B cell members of this lineage are predominantly atypical in phenotype yet also span the class-switched and antibody-secreting cell compartments. Development of neutralization breadth occurred concomitantly with extensive recombination between founder viruses before each virus separated into two distinct population "arms" that evolved independently to escape the PC39-1 lineage. Ab crystal structures show an extended CDRH1 that can help stabilize the CDRH3. Overall, these findings suggest that early exposure of the humoral system to multiple related Env molecules could promote the induction of bnAbs by focusing Ab responses to conserved epitopes.

Association of Counseling and Psychotherapy on Retention in Medication for Addiction Treatment Within a Large Medicaid Population.

OBJECTIVES: Buprenorphine/naloxone is an effective medication for the treatment of opioid use disorder. Unlike methadone, which can only be dispensed in federally waived clinics and which must be combined with specific psychosocial treatment, buprenorphine can be dispensed by individual prescribers who have completed an 8-hour training program, with no requirement that patients receive concomitant psychotherapy. The objective of this study is to quantify the association of counseling and psychotherapy on retention in treatment. We also examine the effect of buprenorphine dosage on retention. METHODS: We examined a cohort of 4987 members of a not-for-profit managed care organization serving Medicaid members in 41 counties in Pennsylvania. This cohort was selected from all members who had a full year without any medication for opioid use disorder followed by initiation of treatment with buprenorphine/naloxone in 2016 to 2017 and who remained Medicaid eligible for at least 80% of the following 2 years. Outcomes were estimated using inverse probability weighted propensity scores. RESULTS: The addition of counseling and psychotherapy within the first 8 weeks of treatment was associated with greater total retention in treatment and there was a dose-response relationship. A 16 mg/d or greater dose of buprenorphine was also associated with greater retention. CONCLUSIONS: These results provide support for an integrated approach to treating people with an opioid use disorder, through a combination of buprenorphine pharmacotherapy and targeted counseling and psychotherapy within the first 2 months of treatment.

Inter-class Concomitant Pharmacotherapy in Medicaid-Insured Youth Receiving Psychiatric Residential Treatment.

Background: Concomitant pharmacotherapy has become increasingly common in the treatment of youth, including in psychiatric residential treatment facilities (PRTF) despite limited efficacy and safety data. Research is reported on the prevalence of any class and interclass concomitant pharmacotherapy, specific class combinations of psychotropics, and changes in number of medications from admission to discharge for Medicaid insured youth treated in PRTFs in one mid-Atlantic state. Methods: Medicaid administrative claims data were examined for youth under age 18 years who were discharged from one of 21 PRTFs during calendar year 2019. Descriptive statistics were calculated to examine patterns of service utilization 90 days prior to admission. The rates of concomitant psychotropic use at admission were compared to the rates at discharge. Logistic regression models were used to examine covariates associated with discharging on 4 or more medications. Results: Fifty-four % of youth were admitted on either two or three psychotropics, while 25% were admitted on four or more psychotropics. The proportion of youth admitting and discharging on 2 or 3 medications was stable. There was a 27% increase in number of youth discharging on 4 medications with a 24% decrease in those on a 5- drug regimen. Only the number of medications prescribed at admission was found to be significant (p < 0.001), with more medications at admission contributing to probability of discharging on 4 or more medications. Conclusions: Concomitant pharmacotherapy is common in PRTFs. These findings support the practice of deprescribing and underscore the need for further research.

ilastik: interactive machine learning for (bio)image analysis.

We present ilastik, an easy-to-use interactive tool that brings machine-learning-based (bio)image analysis to end users without substantial computational expertise. It contains pre-defined workflows for image segmentation, object classification, counting and tracking. Users adapt the workflows to the problem at hand by interactively providing sparse training annotations for a nonlinear classifier. ilastik can process data in up to five dimensions (3D, time and number of channels). Its computational back end runs operations on-demand wherever possible, allowing for interactive prediction on data larger than RAM. Once the classifiers are trained, ilastik workflows can be applied to new data from the command line without further user interaction. We describe all ilastik workflows in detail, including three case studies and a discussion on the expected performance.

Long-read amplicon denoising.

Long-read next-generation amplicon sequencing shows promise for studying complete genes or genomes from complex and diverse populations. Current long-read sequencing technologies have challenging error profiles, hindering data processing and incorporation into downstream analyses. Here we consider the problem of how to reconstruct, free of sequencing error, the true sequence variants and their associated frequencies from PacBio reads. Called 'amplicon denoising', this problem has been extensively studied for short-read sequencing technologies, but current solutions do not always successfully generalize to long reads with high indel error rates. We introduce two methods: one that runs nearly instantly and is very accurate for medium length reads and high template coverage, and another, slower method that is more robust when reads are very long or coverage is lower. On two Mock Virus Community datasets with ground truth, each sequenced on a different PacBio instrument, and on a number of simulated datasets, we compare our two approaches to each other and to existing algorithms. We outperform all tested methods in accuracy, with competitive run times even for our slower method, successfully discriminating templates that differ by a just single nucleotide. Julia implementations of Fast Amplicon Denoising (FAD) and Robust Amplicon Denoising (RAD), and a webserver interface, are freely available.

Full-Length Envelope Analyzer (FLEA): A tool for longitudinal analysis of viral amplicons.

Next generation sequencing of viral populations has advanced our understanding of viral population dynamics, the development of drug resistance, and escape from host immune responses. Many applications require complete gene sequences, which can be impossible to reconstruct from short reads. HIV env, the protein of interest for HIV vaccine studies, is exceptionally challenging for long-read sequencing and analysis due to its length, high substitution rate, and extensive indel variation. While long-read sequencing is attractive in this setting, the analysis of such data is not well handled by existing methods. To address this, we introduce FLEA (Full-Length Envelope Analyzer), which performs end-to-end analysis and visualization of long-read sequencing data. FLEA consists of both a pipeline (optionally run on a high-performance cluster), and a client-side web application that provides interactive results. The pipeline transforms FASTQ reads into high-quality consensus sequences (HQCSs) and uses them to build a codon-aware multiple sequence alignment. The resulting alignment is then used to infer phylogenies, selection pressure, and evolutionary dynamics. The web application provides publication-quality plots and interactive visualizations, including an annotated viral alignment browser, time series plots of evolutionary dynamics, visualizations of gene-wide selective pressures (such as dN/dS) across time and across protein structure, and a phylogenetic tree browser. We demonstrate how FLEA may be used to process Pacific Biosciences HIV env data and describe recent examples of its use. Simulations show how FLEA dramatically reduces the error rate of this sequencing platform, providing an accurate portrait of complex and variable HIV env populations. A public instance of FLEA is hosted at http://flea.datamonkey.org. The Python source code for the FLEA pipeline can be found at https://github.com/veg/flea-pipeline. The client-side application is available at https://github.com/veg/flea-web-app. A live demo of the P018 results can be found at http://flea.murrell.group/view/P018.

HIV Superinfection Drives De Novo Antibody Responses and Not Neutralization Breadth.

Eliciting antibodies that neutralize a broad range of circulating HIV strains (broadly neutralizing antibodies [bnAbs]) represents a key priority for vaccine development. HIV superinfection (re-infection with a second strain following an established infection) has been associated with neutralization breadth, and can provide insights into how the immune system responds to sequential exposure to distinct HIV envelope glycoproteins (Env). Characterizing the neutralizing antibody (nAb) responses in four superinfected women revealed that superinfection does not boost memory nAb responses primed by the first infection or promote nAb responses to epitopes conserved in both infecting viruses. While one superinfected individual developed potent bnAbs, superinfection was likely not the driver as the nAb response did not target an epitope conserved in both viruses. Rather, sequential exposure led to nAbs specific to each Env but did not promote bnAb development. Thus, sequential immunization with heterologous Envs may not be sufficient to focus the immune response onto conserved epitopes.

RIFRAF: a frame-resolving consensus algorithm.

Motivation: Protein coding genes can be studied using long-read next generation sequencing. However, high rates of indel sequencing errors are problematic, corrupting the reading frame. Even the consensus of multiple independent sequence reads retains indel errors. To solve this problem, we introduce Reference-Informed Frame-Resolving multiple-Alignment Free template inference algorithm (RIFRAF), a sequence consensus algorithm that takes a set of error-prone reads and a reference sequence and infers an accurate in-frame consensus. RIFRAF uses a novel structure, analogous to a two-layer hidden Markov model: the consensus is optimized to maximize alignment scores with both the set of noisy reads and with a reference. The template-to-reads component of the model encodes the preponderance of indels, and is sensitive to the per-base quality scores, giving greater weight to more accurate bases. The reference-to-template component of the model penalizes frame-destroying indels. A local search algorithm proceeds in stages to find the best consensus sequence for both objectives. Results: Using Pacific Biosciences SMRT sequences from an HIV-1 env clone, NL4-3, we compare our approach to other consensus and frame correction methods. RIFRAF consistently finds a consensus sequence that is more accurate and in-frame, especially with small numbers of reads. It was able to perfectly reconstruct over 80% of consensus sequences from as few as three reads, whereas the best alternative required twice as many. RIFRAF is able to achieve these results and keep the consensus in-frame even with a distantly related reference sequence. Moreover, unlike other frame correction methods, RIFRAF can detect and keep true indels while removing erroneous ones. Availability and implementation: RIFRAF is implemented in Julia, and source code is publicly available at https://github.com/MurrellGroup/Rifraf.jl. Supplementary information: Supplementary data are available at Bioinformatics online.

HIV Envelope Glycoform Heterogeneity and Localized Diversity Govern the Initiation and Maturation of a V2 Apex Broadly Neutralizing Antibody Lineage.

Understanding how broadly neutralizing antibodies (bnAbs) to HIV envelope (Env) develop during natural infection can help guide the rational design of an HIV vaccine. Here, we described a bnAb lineage targeting the Env V2 apex and the Ab-Env co-evolution that led to development of neutralization breadth. The lineage Abs bore an anionic heavy chain complementarity-determining region 3 (CDRH3) of 25 amino acids, among the shortest known for this class of Abs, and achieved breadth with only 10% nucleotide somatic hypermutation and no insertions or deletions. The data suggested a role for Env glycoform heterogeneity in the activation of the lineage germline B cell. Finally, we showed that localized diversity at key V2 epitope residues drove bnAb maturation toward breadth, mirroring the Env evolution pattern described for another donor who developed V2-apex targeting bnAbs. Overall, these findings suggest potential strategies for vaccine approaches based on germline-targeting and serial immunogen design.

Antibody 10-1074 suppresses viremia in HIV-1-infected individuals.

Monoclonal antibody 10-1074 targets the V3 glycan supersite on the HIV-1 envelope (Env) protein. It is among the most potent anti-HIV-1 neutralizing antibodies isolated so far. Here we report on its safety and activity in 33 individuals who received a single intravenous infusion of the antibody. 10-1074 was well tolerated and had a half-life of 24.0 d in participants without HIV-1 infection and 12.8 d in individuals with HIV-1 infection. Thirteen individuals with viremia received the highest dose of 30 mg/kg 10-1074. Eleven of these participants were 10-1074-sensitive and showed a rapid decline in viremia by a mean of 1.52 log10 copies/ml. Virologic analysis revealed the emergence of multiple independent 10-1074-resistant viruses in the first weeks after infusion. Emerging escape variants were generally resistant to the related V3-specific antibody PGT121, but remained sensitive to antibodies targeting nonoverlapping epitopes, such as the anti-CD4-binding-site antibodies 3BNC117 and VRC01. The results demonstrate the safety and activity of 10-1074 in humans and support the idea that antibodies targeting the V3 glycan supersite might be useful for the treatment and prevention of HIV-1 infection.

Early Antibody Lineage Diversification and Independent Limb Maturation Lead to Broad HIV-1 Neutralization Targeting the Env High-Mannose Patch.

The high-mannose patch on HIV Env is a preferred target for broadly neutralizing antibodies (bnAbs), but to date, no vaccination regimen has elicited bnAbs against this region. Here, we present the development of a bnAb lineage targeting the high-mannose patch in an HIV-1 subtype-C-infected donor from sub-Saharan Africa. The Abs first acquired autologous neutralization, then gradually matured to achieve breadth. One Ab neutralized >47% of HIV-1 strains with only ∼11% somatic hypermutation and no insertions or deletions. By sequencing autologous env, we determined key residues that triggered the lineage and participated in Ab-Env coevolution. Next-generation sequencing of the Ab repertoire showed an early expansive diversification of the lineage followed by independent maturation of individual limbs, several of them developing notable breadth and potency. Overall, the findings are encouraging from a vaccine standpoint and suggest immunization strategies mimicking the evolution of the entire high-mannose patch and promoting maturation of multiple diverse Ab pathways.

Influence of Lachancea thermotolerans on cv. Emir wine fermentation.

The present paper describes the behaviour of Lachancea thermotolerans and Saccharomyces cerevisiae in pure, co-cultured and sequential fermentations in cv. Emir grape must. Faster fermentation rates were observed in wine made with a pure culture of S. cerevisiae and wine produced with simultaneously inoculated cultures of L. thermotolerans and S. cerevisiae. Both L. thermotolerans and S. cerevisiae gave high population numbers. The use of L. thermotolerans in mixed and sequential cultures led to an increase in final total acidity content in the wines, varying in the range 5.40-6.28 g/l (as tartaric acid), compared to pure culture S. cerevisiae, which gave the lowest level of total acidity (5 g/l). The increase was in the order of 1.18-2.06 g/l total acidity. Increase in final acidity by the use of L. thermotolerans might be useful to improve wines with low acidity due to global climate change. Volatile acidity levels (as acetic acid) were in the range 0.53-0.73 g/l, while the concentration of ethyl alcohol varied in the range 10.76-11.62% v/v. Sequential fermentations of wines and pure culture fermentation of L. thermotolerans resulted in reduction in the concentrations of acetaldehyde and higher alcohols, with exception of N-propanol and esters. According to the sensory analysis, wine obtained with sequential inoculation of L. thermotolerans followed by inoculation of S. cerevisiae after 24 h, and simultaneous inoculation of these yeasts, was the most preferred. Copyright © 2016 John Wiley & Sons, Ltd.

Querying Co-regulated Genes on Diverse Gene Expression Datasets Via Biclustering.

Rapid development and increasing popularity of gene expression microarrays have resulted in a number of studies on the discovery of co-regulated genes. One important way of discovering such co-regulations is the query-based search since gene co-expressions may indicate a shared role in a biological process. Although there exist promising query-driven search methods adapting clustering, they fail to capture many genes that function in the same biological pathway because microarray datasets are fraught with spurious samples or samples of diverse origin, or the pathways might be regulated under only a subset of samples. On the other hand, a class of clustering algorithms known as biclustering algorithms which simultaneously cluster both the items and their features are useful while analyzing gene expression data, or any data in which items are related in only a subset of their samples. This means that genes need not be related in all samples to be clustered together. Because many genes only interact under specific circumstances, biclustering may recover the relationships that traditional clustering algorithms can easily miss. In this chapter, we briefly summarize the literature using biclustering for querying co-regulated genes. Then we present a novel biclustering approach and evaluate its performance by a thorough experimental analysis.

Rapid Sequencing of Complete env Genes from Primary HIV-1 Samples.

The ability to study rapidly evolving viral populations has been constrained by the read length of next-generation sequencing approaches and the sampling depth of single-genome amplification methods. Here, we develop and characterize a method using Pacific Biosciences' Single Molecule, Real-Time (SMRT®) sequencing technology to sequence multiple, intact full-length human immunodeficiency virus-1 env genes amplified from viral RNA populations circulating in blood, and provide computational tools for analyzing and visualizing these data.

Gene-wide identification of episodic selection.

We present BUSTED, a new approach to identifying gene-wide evidence of episodic positive selection, where the non-synonymous substitution rate is transiently greater than the synonymous rate. BUSTED can be used either on an entire phylogeny (without requiring an a priori hypothesis regarding which branches are under positive selection) or on a pre-specified subset of foreground lineages (if a suitable a priori hypothesis is available). Selection is modeled as varying stochastically over branches and sites, and we propose a computationally inexpensive evidence metric for identifying sites subject to episodic positive selection on any foreground branches. We compare BUSTED with existing models on simulated and empirical data. An implementation is available on www.datamonkey.org/busted, with a widget allowing the interactive specification of foreground branches.

A comparative analysis of biclustering algorithms for gene expression data.

The need to analyze high-dimension biological data is driving the development of new data mining methods. Biclustering algorithms have been successfully applied to gene expression data to discover local patterns, in which a subset of genes exhibit similar expression levels over a subset of conditions. However, it is not clear which algorithms are best suited for this task. Many algorithms have been published in the past decade, most of which have been compared only to a small number of algorithms. Surveys and comparisons exist in the literature, but because of the large number and variety of biclustering algorithms, they are quickly outdated. In this article we partially address this problem of evaluating the strengths and weaknesses of existing biclustering methods. We used the BiBench package to compare 12 algorithms, many of which were recently published or have not been extensively studied. The algorithms were tested on a suite of synthetic data sets to measure their performance on data with varying conditions, such as different bicluster models, varying noise, varying numbers of biclusters and overlapping biclusters. The algorithms were also tested on eight large gene expression data sets obtained from the Gene Expression Omnibus. Gene Ontology enrichment analysis was performed on the resulting biclusters, and the best enrichment terms are reported. Our analyses show that the biclustering method and its parameters should be selected based on the desired model, whether that model allows overlapping biclusters, and its robustness to noise. In addition, we observe that the biclustering algorithms capable of finding more than one model are more successful at capturing biologically relevant clusters.