Profiling microbial communities in an extremely acidic environment influenced by a cold natural carbon dioxide spring: A study of the Mefite in Ansanto Valley, Southern Italy.

The Ansanto Valley's Mefite, one of the Earth's largest non-volcanic CO2 gas emissions, is distinguished by its cold natural carbon dioxide springs. These emissions originate from the intricate tectonics and geodynamics of the southern Apennines in Italy. Known for over two millennia for its lethal concentration of CO2 and other harmful gases, the Mefite has a reputation for being toxic and dangerous. Despite its historical significance and unique geological features, there is a lack of information on the microbial diversity associated with the Mefite's gas emissions. This study presents an integrated exploration of the microbial diversity in the mud soil, using high-throughput sequencing of 16S rRNA (Prokaryotes) and ITS2 (Fungi), alongside a geochemical site characterisation. Our findings reveal that the Mefite's unique environment imposes a significant bottleneck on microbial diversity, favouring a select few microbial groups such as Actinobacteria and Firmicutes for Prokaryotes, and Basidiomycota for Fungi.

Optimizing Fecal Occult Blood Test (FOBT) Colorectal Cancer Screening Using Gut Bacteriome as a Biomarker.

BACKGROUND: Colorectal cancer (CRC) is a major cause of cancer mortality in the world. One of the most widely used screening tests for CRC is the immunochemical fecal occult blood test (iFOBT), which detects human hemoglobin from patient's stool sample. Although it is highly efficient in detecting blood from patients with gastro-intestinal lesions, such as polyps and cancers, the iFOBT has a high rate of false positive discovery. Recent studies suggested gut bacteria as a promising noninvasive biomarker for improving the diagnosis of CRC. In this study, we examined the composition of gut bacteria using iFOBT leftover from patients undergoing screening test along with a colonoscopy. METHODS: After collecting data from more than 800 patients, we considered 4 groups for this study. The first and second groups were respectively "healthy" in which the patients had either no blood in their stool or had blood but no lesions. The third and fourth groups of patients had both blood in their stools with precancerous and cancerous lesions and considered either as low-grade and high-grade lesion groups, respectively. An amplification of 16S rRNA (V4 region) gene was performed, followed by sequencing along with various statistical and bioinformatic analysis. RESULTS: We analyzed the composition of the gut bacteriome at phylum, class, genus, and species levels. Although members of the Firmicute phylum increased in the 3 groups compared to healthy patients, the phylum Actinobacteriota was found to decrease. Moreover, Blautia obeum and Anaerostipes hadrus from the phylum Firmicutes were increased and Collinsella aerofaciens from phylum Actinobacteriota was found decreased when healthy group is compared to the patients with high-grade lesions. Finally, among the 5 machine learning algorithms used to perform our analysis, both elastic net (AUC > 0.7) and random forest (AUC > 0.8) performs well in differentiating healthy patients from 3 other patient groups having blood in their stool. CONCLUSION: Our study integrates the iFOBT screening tool with gut bacterial composition to improve the prediction of CRC lesions.

Ribosomal Protein uS5 and Friends: Protein-Protein Interactions Involved in Ribosome Assembly and Beyond.

Ribosomal proteins are fundamental components of the ribosomes in all living cells. The ribosomal protein uS5 (Rps2) is a stable component of the small ribosomal subunit within all three domains of life. In addition to its interactions with proximal ribosomal proteins and rRNA inside the ribosome, uS5 has a surprisingly complex network of evolutionarily conserved non-ribosome-associated proteins. In this review, we focus on a set of four conserved uS5-associated proteins: the protein arginine methyltransferase 3 (PRMT3), the programmed cell death 2 (PDCD2) and its PDCD2-like (PDCD2L) paralog, and the zinc finger protein, ZNF277. We discuss recent work that presents PDCD2 and homologs as a dedicated uS5 chaperone and PDCD2L as a potential adaptor protein for the nuclear export of pre-40S subunits. Although the functional significance of the PRMT3-uS5 and ZNF277-uS5 interactions remain elusive, we reflect on the potential roles of uS5 arginine methylation by PRMT3 and on data indicating that ZNF277 and PRMT3 compete for uS5 binding. Together, these discussions highlight the complex and conserved regulatory network responsible for monitoring the availability and the folding of uS5 for the formation of 40S ribosomal subunits and/or the role of uS5 in potential extra-ribosomal functions.

Distinct regulatory functions and biological roles of lncRNA splice variants.

Alternative splicing (AS) of RNA molecules is a key contributor to transcriptome diversity. In humans, 90%-95% of multi-exon genes produce alternatively spliced RNA transcripts. Therefore, every single gene has the opportunity of producing multiple splice variants, including long non-coding RNA (lncRNA) genes that undergo RNA maturation steps such as conventional and alternative splicing. Emerging evidence suggests significant roles for these lncRNA splice variants in many aspects of cell biology. Differential changes in expression of specific lncRNA splice variants have also been associated with many diseases including cancer. This review covers the current knowledge on this emerging topic of investigation. We provide exclusive insights on the AS landscape of lncRNAs and also describe at the molecular level the functional relevance of lncRNA splice variants, i.e., RNA-based differential functions, production of micropeptides, and generation of circular RNAs. Finally, we discuss exciting perspectives for this emerging field and outline the work required to further develop research endeavors in this field.

HNF4α Acts as Upstream Functional Regulator of Intestinal Wnt3 and Paneth Cell Fate.

BACKGROUND & AIMS: The intestinal epithelium intrinsically renews itself ex vivo via the proliferation of Lgr5+ intestinal stem cells, which is sustained by the establishment of an epithelial stem cell niche. Differentiated Paneth cells are the main source of epithelial-derived niche factor supplies and produce Wnt3 as an essential factor in supporting Lgr5+ stem cell activity in the absence of redundant mesenchymal Wnts. Maturation of Paneth cells depends on canonical Wnt signaling, but few transcriptional regulators have been identified to this end. The role of HNF4α in intestinal epithelial cell differentiation is considered redundant with its paralog HNF4γ. However, it is unclear whether HNF4α alone controls intrinsic intestinal epithelial cell growth and fate in the absence of a mesenchymal niche. METHODS: We used transcriptomic analyses to dissect the role of HNF4α in the maintenance of jejunal epithelial culture when cultured ex vivo as enteroids in the presence or absence of compensatory mesenchymal cells. RESULTS: HNF4α plays a crucial role in supporting the growth and survival of jejunal enteroids. Transcriptomic analyses revealed an autonomous function of HNF4α in Wnt3 transcriptional regulation and Paneth cell differentiation. We showed that Wnt3a supplementation or co-culture with intestinal subepithelial mesenchymal cells reversed cell death and transcriptional changes caused by the deletion of Hnf4a in jejunal enteroids. CONCLUSIONS: Our results support the intrinsic epithelial role of HNF4α in regulating Paneth cell homeostasis and intestinal epithelium renewal in the absence of compensatory Wnt signaling.

High-level resistance to bictegravir and cabotegravir in subtype A- and D-infected HIV-1 patients failing raltegravir with multiple resistance mutations.

OBJECTIVES: The second-generation integrase strand transfer inhibitor (INSTI) bictegravir is becoming accessible in low- and middle-income countries (LMICs), and another INSTI, cabotegravir, has recently been approved as a long-acting injectable. Data on bictegravir and cabotegravir susceptibility in raltegravir-experienced HIV-1 subtype A- and D-infected patients carrying drug resistance mutations (DRMs) remain very scarce in LMICs. PATIENTS AND METHODS: HIV-1 integrase (IN)-recombinant viruses from eight patients failing raltegravir-based third-line therapy in Uganda were genotypically and phenotypically tested for susceptibility to bictegravir and cabotegravir. Ability of these viruses to integrate into human genomes was assessed in MT-4 cells. RESULTS: HIV-1 IN-recombinant viruses harbouring single primary mutations (N155H or Y143R/S) or in combination with secondary INSTI mutations (T97A, M50I, L74IM, E157Q, G163R or V151I) were susceptible to both bictegravir and cabotegravir. However, combinations of primary INSTI-resistance mutations such as E138A/G140A/G163R/Q148R or E138K/G140A/S147G/Q148K led to decreased susceptibility to both cabotegravir (fold change in EC50 values from 429 to 1000×) and bictegravir (60 to 100×), exhibiting a high degree of cross-resistance. However, these same IN-recombinant viruses showed impaired integration capacity (14% to 48%) relative to the WT HIV-1 NL4-3 strain in the absence of drug. CONCLUSIONS: Though not currently widely accessible in most LMICs, bictegravir and cabotegravir offer a valid alternative to HIV-infected individuals harbouring subtype A and D HIV-1 variants with reduced susceptibility to first-generation INSTIs but previous exposure to raltegravir may reduce efficacy, more so with cabotegravir.

OpenProt 2021: deeper functional annotation of the coding potential of eukaryotic genomes.

OpenProt (www.openprot.org) is the first proteogenomic resource supporting a polycistronic annotation model for eukaryotic genomes. It provides a deeper annotation of open reading frames (ORFs) while mining experimental data for supporting evidence using cutting-edge algorithms. This update presents the major improvements since the initial release of OpenProt. All species support recent NCBI RefSeq and Ensembl annotations, with changes in annotations being reported in OpenProt. Using the 131 ribosome profiling datasets re-analysed by OpenProt to date, non-AUG initiation starts are reported alongside a confidence score of the initiating codon. From the 177 mass spectrometry datasets re-analysed by OpenProt to date, the unicity of the detected peptides is controlled at each implementation. Furthermore, to guide the users, detectability statistics and protein relationships (isoforms) are now reported for each protein. Finally, to foster access to deeper ORF annotation independently of one's bioinformatics skills or computational resources, OpenProt now offers a data analysis platform. Users can submit their dataset for analysis and receive the results from the analysis by OpenProt. All data on OpenProt are freely available and downloadable for each species, the release-based format ensuring a continuous access to the data. Thus, OpenProt enables a more comprehensive annotation of eukaryotic genomes and fosters functional proteomic discoveries.

Hemeprotein Tpx1 interacts with cell-surface heme transporter Str3 in Schizosaccharomyces pombe.

Str3 is a transmembrane protein that mediates low-affinity heme uptake in Schizosaccharomyces pombe. Under iron-limiting conditions, Str3 remains at the cell surface in the presence of increasing hemin concentrations. Using a proximity-dependent biotinylation approach coupled to mass spectrometry and coimmunoprecipitation assays, we report that the peroxiredoxin Tpx1 is a binding partner of Str3. Under microaerobic conditions, cells deficient in heme biosynthesis and lacking the heme receptor Shu1 exhibit poor hemin-dependent growth in the absence of Tpx1. Analysis of membrane protein preparations from iron-starved hem1Δ shu1Δ str3Δ tpx1Δ cells coexpressing Str3-GFP and TAP-Tpx1 showed that TAP-Tpx1 is enriched in membrane protein fractions in response to hemin. Bimolecular fluorescence complementation assays brought additional evidence that an interaction between Tpx1 and Str3 occurs at the plasma membrane. Results showed that Tpx1 exhibits an equilibrium constant value of 0.26 µM for hemin. The association of Tpx1 with hemin protects hemin from degradation by H2 O2 . The peroxidase activity of hemin is lowered when it is bound to Tpx1. Taken together, these results revealed that Tpx1 is a novel interacting partner of Str3. Our data are the first example of an interaction between a cytoplasmic heme-binding protein and a cell-surface heme transporter.

Accumulation of integrase strand transfer inhibitor resistance mutations confers high-level resistance to dolutegravir in non-B subtype HIV-1 strains from patients failing raltegravir in Uganda.

BACKGROUND: Increasing first-line treatment failures in low- and middle-income countries (LMICs) have led to increased use of integrase strand transfer inhibitors (INSTIs) such as dolutegravir. However, HIV-1 susceptibility to INSTIs in LMICs, especially with previous raltegravir exposure, is poorly understood due to infrequent reporting of INSTI failures and testing for INSTI drug resistance mutations (DRMs). METHODS: A total of 51 non-subtype B HIV-1 infected patients failing third-line (raltegravir-based) therapy in Uganda were initially selected for the study. DRMs were detected using Sanger and deep sequencing. HIV integrase genes of 13 patients were cloned and replication capacities (RCs) and phenotypic susceptibilities to dolutegravir, raltegravir and elvitegravir were determined with TZM-bl cells. Spearman's correlation coefficient was used to determine cross-resistance between INSTIs. RESULTS: INSTI DRMs were detected in 47% of patients. HIV integrase-recombinant virus carrying one primary INSTI DRM (N155H or Y143R/S) was susceptible to dolutegravir but highly resistant to raltegravir and elvitegravir (>50-fold change). Two patients, one with E138A/G140A/Q148R/G163R and one with E138K/G140A/S147G/Q148K, displayed the highest reported resistance to raltegravir, elvitegravir and even dolutegravir. The former multi-DRM virus had WT RC whereas the latter had lower RCs than WT. CONCLUSIONS: In HIV-1 subtype A- and D-infected patients failing raltegravir and harbouring INSTI DRMs, there is high-level resistance to elvitegravir and raltegravir. More routine monitoring of INSTI treatment may be advised in LMICs, considering that multiple INSTI DRMs may have accumulated during prolonged exposure to raltegravir during virological failure, leading to high-level INSTI resistance, including dolutegravir resistance.

Janus Kinase Mutations in Mice Lacking PU.1 and Spi-B Drive B Cell Leukemia through Reactive Oxygen Species-Induced DNA Damage.

Precursor B cell acute lymphoblastic leukemia (B-ALL) is caused by genetic lesions in developing B cells that function as drivers for the accumulation of additional mutations in an evolutionary selection process. We investigated secondary drivers of leukemogenesis in a mouse model of B-ALL driven by PU.1/Spi-B deletion (Mb1-CreΔPB). Whole-exome-sequencing analysis revealed recurrent mutations in Jak3 (encoding Janus kinase 3), Jak1, and Ikzf3 (encoding Aiolos). Mutations with a high variant-allele frequency (VAF) were dominated by C→T transition mutations that were compatible with activation-induced cytidine deaminase, whereas the majority of mutations, with a low VAF, were dominated by C→A transversions associated with 8-oxoguanine DNA damage caused by reactive oxygen species (ROS). The Janus kinase (JAK) inhibitor ruxolitinib delayed leukemia onset, reduced ROS and ROS-induced gene expression signatures, and altered ROS-induced mutational signatures. These results reveal that JAK mutations can alter the course of leukemia clonal evolution through ROS-induced DNA damage.

A targeted reactivation of latent HIV-1 using an activator vector in patient samples from acute infection.

BACKGROUND: During combined anti-retroviral treatment, a latent HIV reservoir persists within resting memory CD4 T cells that initiates viral recrudescence upon treatment interruption. Strategies for HIV-1 cure have largely focused on latency reversing agents (LRAs) capable of reactivating and eliminating this viral reservoir. Previously investigated LRAs have largely failed to achieve a robust latency reversal sufficient for reduction of latent HIV pool or the potential of virus-free remission in the absence of treatment. METHODS: We utilize a polyvalent virus-like particle (VLP) formulation called Activator Vector (ACT-VEC) to 'shock' provirus into transcriptional activity. Ex vivo co-culture experiments were used to evaluate the efficacy of ACT-VEC in relation to other LRAs in individuals diagnosed and treated during the acute stage of infection. IFN-γ ELISpot, qRT-PCR and Illumina MiSeq were used to evaluate antigenicity, latency reversal, and diversity of induced virus respectively. FINDINGS: Using samples from HIV+ patients diagnosed and treated at acute/early infection, we demonstrate that ACT-VEC can reverse latency in HIV infected CD4 T cells to a greater extent than other major recall antigens as stimuli or even mitogens such as PMA/Iono. Furthermore, ACT-VEC activates more latent HIV-1 than clinically tested HDAC inhibitors or protein kinase C agonists. INTERPRETATION: Taken together, these results show that ACT-VEC can induce HIV reactivation from latently infected CD4 T cells collected from participants on first line combined antiretroviral therapy for at least two years after being diagnosed and treated at acute/early stage of infection. These findings could provide guidance to possible targeted cure strategies and treatments. FUNDING: NIH and CIHR.

Human Hepatocyte Nuclear Factor 4-α Encodes Isoforms with Distinct Transcriptional Functions.

HNF4α is a nuclear receptor produced as 12 isoforms from two promoters by alternative splicing. To characterize the transcriptional capacities of all 12 HNF4α isoforms, stable lines expressing each isoform were generated. The entire transcriptome associated with each isoform was analyzed as well as their respective interacting proteome. Major differences were noted in the transcriptional function of these isoforms. The α1 and α2 isoforms were the strongest regulators of gene expression whereas the α3 isoform exhibited significantly reduced activity. The α4, α5, and α6 isoforms, which use an alternative first exon, were characterized for the first time, and showed a greatly reduced transcriptional potential with an inability to recognize the consensus response element of HNF4α. Several transcription factors and coregulators were identified as potential specific partners for certain HNF4α isoforms. An analysis integrating the vast amount of omics data enabled the identification of transcriptional regulatory mechanisms specific to certain HNF4α isoforms, hence demonstrating the importance of considering all isoforms given their seemingly diverse functions.

Tree shape-based approaches for the comparative study of cophylogeny.

Cophylogeny is the congruence of phylogenetic relationships between two different groups of organisms due to their long-term interaction. We investigated the use of tree shape distance measures to quantify the degree of cophylogeny. We implemented a reverse-time simulation model of pathogen phylogenies within a fixed host tree, given cospeciation probability, host switching, and pathogen speciation rates. We used this model to evaluate 18 distance measures between host and pathogen trees including two kernel distances that we developed for labeled and unlabeled trees, which use branch lengths and accommodate different size trees. Finally, we used these measures to revisit published cophylogenetic studies, where authors described the observed associations as representing a high or low degree of cophylogeny. Our simulations demonstrated that some measures are more informative than others with respect to specific coevolution parameters especially when these did not assume extreme values. For real datasets, trees' associations projection revealed clustering of high concordance studies suggesting that investigators are describing it in a consistent way. Our results support the hypothesis that measures can be useful for quantifying cophylogeny. This motivates their usage in the field of coevolution and supports the development of simulation-based methods, i.e., approximate Bayesian computation, to estimate the underlying coevolutionary parameters.

Evidence for a recombinant origin of HIV-1 Group M from genomic variation.

Reconstructing the early dynamics of the HIV-1 pandemic can provide crucial insights into the socioeconomic drivers of emerging infectious diseases in human populations, including the roles of urbanization and transportation networks. Current evidence indicates that the global pandemic comprising almost entirely of HIV-1/M originated around the 1920s in central Africa. However, these estimates are based on molecular clock estimates that are assumed to apply uniformly across the virus genome. There is growing evidence that recombination has played a significant role in the early history of the HIV-1 pandemic, such that different regions of the HIV-1 genome have different evolutionary histories. In this study, we have conducted a dated-tip analysis of all near full-length HIV-1/M genome sequences that were published in the GenBank database. We used a sliding window approach similar to the 'bootscanning' method for detecting breakpoints in inter-subtype recombinant sequences. We found evidence of substantial variation in estimated root dates among windows, with an estimated mean time to the most recent common ancestor of 1922. Estimates were significantly autocorrelated, which was more consistent with an early recombination event than with stochastic error variation in phylogenetic reconstruction and dating analyses. A piecewise regression analysis supported the existence of at least one recombination breakpoint in the HIV-1/M genome with interval-specific means around 1929 and 1913, respectively. This analysis demonstrates that a sliding window approach can accommodate early recombination events outside the established nomenclature of HIV-1/M subtypes, although it is difficult to incorporate the earliest available samples due to their limited genome coverage.

First-line HIV treatment failures in non-B subtypes and recombinants: a cross-sectional analysis of multiple populations in Uganda.

BACKGROUND: Our understanding of HIV-1 and antiretroviral treatment (ART) is strongly biased towards subtype B, the predominant subtype in North America and western Europe. Efforts to characterize the response to first-line treatments in other HIV-1 subtypes have been hindered by the availability of large study cohorts in resource-limited settings. To maximize our statistical power, we combined HIV-1 sequence and clinical data from every available study population associated with the Joint Clinical Research Centre (JCRC) in Uganda. These records were combined with contemporaneous ART-naive records from Uganda in the Stanford HIVdb database. METHODS: Treatment failures were defined by the presence of HIV genotype records with sample collection dates after the ART start dates in the JCRC database. Drug resistances were predicted by the Stanford HIVdb algorithm, and HIV subtype classification and recombination detection was performed with SCUEAL. We used Bayesian network analysis to evaluate associations between drug exposures and subtypes, and binomial regression for associations with recombination. RESULTS: This is the largest database of first-line treatment failures ([Formula: see text]) in Uganda to date, with a predicted statistical power of 80% to detect subtype associations at an odds ratio of [Formula: see text]. In the subset where drug regimen data were available, we observed that use of 3TC was associated with a higher rate of first line treatment failure, whereas regimens containing AZT and TDF were associated with reduced rates of failure. In the complete database, we found limited evidence of associations between HIV-1 subtypes and treatment failure, with the exception of a significantly lower frequency of failures among A/D recombinants that comprised about 7% of the population. First-line treatment failure was significantly associated with reduced numbers of recombination breakpoints across subtypes. CONCLUSIONS: Expanding access to first-line ART should confer the anticipated public health benefits in Uganda, despite known differences in the pathogenesis of HIV-1 subtypes. Furthermore, the impact of ART may actually be enhanced by frequent inter-subtype recombination in this region.

Detecting Amino Acid Coevolution with Bayesian Graphical Models.

The comparative study of homologous proteins can provide abundant information about the functional and structural constraints on protein evolution. For example, an amino acid substitution that is deleterious may become permissive in the presence of another substitution at a second site of the protein. A popular approach for detecting coevolving residues is by looking for correlated substitution events on branches of the molecular phylogeny relating the protein-coding sequences. Here we describe a machine learning method (Bayesian graphical models) implemented in the open-source phylogenetic software package HyPhy, http://hyphy.org , for extracting a network of coevolving residues from a sequence alignment.

An open-source k-mer based machine learning tool for fast and accurate subtyping of HIV-1 genomes.

For many disease-causing virus species, global diversity is clustered into a taxonomy of subtypes with clinical significance. In particular, the classification of infections among the subtypes of human immunodeficiency virus type 1 (HIV-1) is a routine component of clinical management, and there are now many classification algorithms available for this purpose. Although several of these algorithms are similar in accuracy and speed, the majority are proprietary and require laboratories to transmit HIV-1 sequence data over the network to remote servers. This potentially exposes sensitive patient data to unauthorized access, and makes it impossible to determine how classifications are made and to maintain the data provenance of clinical bioinformatic workflows. We propose an open-source supervised and alignment-free subtyping method (Kameris) that operates on k-mer frequencies in HIV-1 sequences. We performed a detailed study of the accuracy and performance of subtype classification in comparison to four state-of-the-art programs. Based on our testing data set of manually curated real-world HIV-1 sequences (n = 2, 784), Kameris obtained an overall accuracy of 97%, which matches or exceeds all other tested software, with a processing rate of over 1,500 sequences per second. Furthermore, our fully standalone general-purpose software provides key advantages in terms of data security and privacy, transparency and reproducibility. Finally, we show that our method is readily adaptable to subtype classification of other viruses including dengue, influenza A, and hepatitis B and C virus.

Sanger and next generation sequencing in the characterisation of arbuscular mycorrhizal fungi (AMF) in Pancratium maritimum L. (Amaryllidaceae), a representative plant species of Mediterranean sand dunes.

MAIN CONCLUSION: An interesting AMF colonization microcosm has been detected in the roots of Pancratium maritimum (sea daffodil). Both sequencing techniques (Sanger and NGS) have been used for AMF characterisation, showing a balanced trade-off between pros and cons. By Sanger and next generation sequencing of rRNA nuclear molecular markers (SSU-ITS-LSU and ITS2, respectively), the presence of AMF communities in the roots of P. maritimum was evaluated. Our results shed light on the presence of AMF in sea daffodil and the diversity of assemblages of AMF detected after Sanger sequencing of the SSU-ITS-LSU marker is much higher than that determined following NGS sequencing of ITS2 alone.

Absence of HIV-1 Drug Resistance Mutations Supports the Use of Dolutegravir in Uganda.

To screen for drug resistance and possible treatment with Dolutegravir (DTG) in treatment-naive patients and those experiencing virologic failure during first-, second-, and third-line combined antiretroviral therapy (cART) in Uganda. Samples from 417 patients in Uganda were analyzed for predicted drug resistance upon failing a first- (N = 158), second- (N = 121), or third-line [all 51 involving Raltegravir (RAL)] treatment regimen. HIV-1 pol gene was amplified and sequenced from plasma samples. Drug susceptibility was interpreted using the Stanford HIV database algorithm and SCUEAL was used for HIV-1 subtyping. Frequency of resistance to nucleoside reverse transcriptase inhibitors (NRTIs) (95%) and non-NRTI (NNRTI, 96%) was high in first-line treatment failures. Despite lack of NNRTI-based treatment for years, NNRTI resistance remained stable in 55% of patients failing second-line or third-line treatment, and was also at 10% in treatment-naive Ugandans. DTG resistance (n = 366) was not observed in treatment-naive individuals or individuals failing first- and second-line cART, and only found in two patients failing third-line cART, while 47% of the latter had RAL- and Elvitegravir-resistant HIV-1. Secondary mutations associated with DTG resistance were found in 2%-10% of patients failing third-line cART. Of 14 drugs currently available for cART in Uganda, resistance was readily observed to all antiretroviral drugs (except for DTG) in Ugandan patients failing first-, second-, or even third-line treatment regimens. The high NNRTI resistance in first-line treatment in Uganda even among treatment-naive patients calls for the use of DTG to reach the UNAIDS 90:90:90 goals.

Beta-Binomial Model for the Detection of Rare Mutations in Pooled Next-Generation Sequencing Experiments.

Against diminishing costs, next-generation sequencing (NGS) still remains expensive for studies with a large number of individuals. As cost saving, sequencing genome of pools containing multiple samples might be used. Currently, there are many software available for the detection of single-nucleotide polymorphisms (SNPs). Sensitivity and specificity depend on the model used and data analyzed, indicating that all software have space for improvement. We use beta-binomial model to detect rare mutations in untagged pooled NGS experiments. We propose a multireference framework for pooled data with ability being specific up to two patients affected by neuromuscular disorders (NMD). We assessed the results comparing with The Genome Analysis Toolkit (GATK), CRISP, SNVer, and FreeBayes. Our results show that the multireference approach applying beta-binomial model is accurate in predicting rare mutations at 0.01 fraction. Finally, we explored the concordance of mutations between the model and software, checking their involvement in any NMD-related gene. We detected seven novel SNPs, for which the functional analysis produced enriched terms related to locomotion and musculature.