HIV-1 subtypes maintain distinctive physicochemical signatures in Nef domains associated with immunoregulation.

BACKGROUND: HIV subtype is associated with varied rates of disease progression. The HIV accessory protein, Nef, continues to be present during antiretroviral therapy (ART) where it has numerous immunoregulatory effects. In this study, we analyzed Nef sequences from HIV subtypes A1, B, C, and D using a machine learning approach that integrates functional amino acid information to identify if unique physicochemical features are associated with Nef functional/structural domains in a subtype-specific manner. METHODS: 2253 sequences representing subtypes A1, B, C, and D were aligned and domains with known functional properties were scored based on amino acid physicochemical properties. Following feature generation, we used statistical pruning and evolved neural networks (ENNs) to determine if we could successfully classify subtypes. Next, we used ENNs to identify the top five key Nef physicochemical features applied to specific immunoregulatory domains that differentiated subtypes. A signature pattern analysis was performed to the assess amino acid diversity in sub-domains that differentiated each subtype. RESULTS: In validation studies, ENNs successfully differentiated each subtype at A1 (87.2%), subtype B (89.5%), subtype C (91.7%), and subtype D (85.1%). Our feature-based domain scoring, followed by t-tests, and a similar ENN identified subtype-specific domain-associated features. Subtype A1 was associated with alterations in Nef CD4 binding domain; subtype B was associated with alterations with the AP-2 Binding domain; subtype C was associated with alterations in a structural Alpha Helix domain; and, subtype D was associated with alterations in a Beta-Sheet domain. CONCLUSIONS: Recent studies have focused on HIV Nef as a driver of immunoregulatory disease in those HIV infected and on ART. Nef acts through a complex mixture of interactions that are directly linked to the key features of the subtype-specific domains we identified with the ENN. The study supports the hypothesis that varied Nef subtypes contribute to subtype-specific disease progression.

Whole-genome sequencing of carbapenem-resistant Enterobacterales isolates in southeast Louisiana reveals persistent genetic clusters spanning multiple locations.

BACKGROUND: We investigated 51 g-negative carbapenem-resistant Enterobacterales (CRE) isolates collected from 22 patients over a five-year period from six health care institutions in the Ochsner Health network in southeast Louisiana. METHODS: Short genomic reads were generated using Illumina sequencing and assembled for each isolate. Isolates were classified as Enterobacter spp. (n = 20), Klebsiella spp. (n = 30), and Escherichia coli (n = 1) and grouped into 19 different multi-locus sequence types (MLST). Species and patient-specific core genomes were constructed representing ∼50% of the chromosomal genome. RESULTS: We identified two sets of patients with genetically related infections; in both cases, the related isolates were collected > 6 months apart, and in one case, the isolates were collected in different locations. On the other hand, we identified four sets of patients with isolates of the same species collected within 21 days from the same location; however, none had genetically related infections. Genes associated with resistance to carbapenem drugs (blaKPC and/or blaCTX-M-15) were found in 76% of the isolates. We found three blaKPC variants (blaKPC-2, blaKPC-3, and blaKPC-4) associated with four different Enterobacter MLST variants, and two blaKPC variants (blaKPC-2, blaKPC-3) associated with seven different Klebsiella MLST variants. CONCLUSIONS: Molecular surveillance is increasingly becoming a powerful tool to understand bacterial spread in both community and clinical settings. This study provides evidence that genetically related infections in clinical settings do not necessarily reflect temporal associations, and vice versa. Our results also highlight the regional genomic and resistance diversity within related bacterial lineages.

Comparing antimicrobial resistant genes and phenotypes across multiple sequencing platforms and assays for Enterobacterales clinical isolates.

INTRODUCTION: Whole genome sequencing (WGS) of bacterial isolates can be used to identify antimicrobial resistance (AMR) genes. Previous studies have shown that genotype-based AMR has variable accuracy for predicting carbapenem resistance in carbapenem-resistant Enterobacterales (CRE); however, the majority of these studies used short-read platforms (e.g. Illumina) to generate sequence data. In this study, our objective was to determine whether Oxford Nanopore Technologies (ONT) long-read WGS would improve detection of carbapenem AMR genes with respect to short-read only WGS for nine clinical CRE samples. We measured the minimum inhibitory breakpoint (MIC) using two phenotype assays (MicroScan and ETEST) for six antibiotics, including two carbapenems (meropenem and ertapenem) and four non-carbapenems (gentamicin, ciprofloxacin, cefepime, and trimethoprim/sulfamethoxazole). We generated short-read data using the Illumina NextSeq and long-read data using the ONT MinION. Four assembly methods were compared: ONT-only assembly; ONT-only assembly plus short-read polish; ONT + short-read hybrid assembly plus short-read polish; short-read only assembly. RESULTS: Consistent with previous studies, our results suggest that the hybrid assembly produced the highest quality results as measured by gene completeness and contig circularization. However, ONT-only methods had minimal impact on the detection of AMR genes and plasmids compared to short-read methods, although, notably, differences in gene copy number differed between methods. All four assembly methods showed identical presence/absence of the blaKPC-2 carbapenemase gene for all samples. The two phenotype assays showed 100% concordant results for the non-carbapenems, but only 65% concordance for the two carbapenems. The presence/absence of AMR genes was 100% concordant with AMR phenotypes for all four non-carbapenem drugs, although only 22%-50% sensitivity for the carbapenems. CONCLUSIONS: Overall, these findings suggest that the lack of complete correspondence between CRE AMR genotype and phenotype for carbapenems, while concerning, is independent of sequencing platform/assembly method.

The Persistence of HIV Diversity, Transcription, and Nef Protein in Kaposi's Sarcoma Tumors during Antiretroviral Therapy.

Epidemic Kaposi's sarcoma (KS), defined by co-infection with Human Herpes Virus 8 (HHV-8) and the Human Immunodeficiency Virus (HIV), is a major cause of mortality in sub-Saharan Africa. Antiretroviral therapy (ART) significantly reduces the risk of developing KS, and for those with KS, tumors frequently resolve with ART alone. However, for unknown reasons, a significant number of KS cases do not resolve and can progress to death. To explore how HIV responds to ART in the KS tumor microenvironment, we sequenced HIV env-nef found in DNA and RNA isolated from plasma, peripheral blood mononuclear cells, and tumor biopsies, before and after ART, in four Ugandan study participants who had unresponsive or progressive KS after 180-250 days of ART. We performed immunohistochemistry experiments to detect viral proteins in matched formalin-fixed tumor biopsies. Our sequencing results showed that HIV diversity and RNA expression in KS tumors are maintained after ART, despite undetectable plasma viral loads. The presence of spliced HIV transcripts in KS tumors after ART was consistent with a transcriptionally active viral reservoir. Immunohistochemistry staining found colocalization of HIV Nef protein and tissue-resident macrophages in the KS tumors. Overall, our results demonstrated that even after ART reduced plasma HIV viral load to undetectable levels and restored immune function, HIV in KS tumors continues to be transcriptionally and translationally active, which could influence tumor maintenance and progression.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes.

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes. Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level. Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method. Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Horizontal transfer and evolution of wall teichoic acid gene cassettes in Bacillus subtilis.

Background: Wall teichoic acid (WTA) genes are essential for production of cell walls in gram-positive bacteria and necessary for survival and variability in the cassette has led to recent antibiotic resistance acquisition in pathogenic bacteria.  Methods: Using a pan-genome approach, we examined the evolutionary history of WTA genes in  Bacillus subtilis ssp.  subtilis.  Results: Our analysis reveals an interesting pattern of evolution from the type-strain WTA gene cassette possibly resulting from horizontal acquisition from organisms with similar gene sequences. The WTA cassettes have a high level of variation which may be due to one or more independent horizontal transfer events during the evolution of  Bacillus subtilis ssp.  subtilis. This swapping of entire WTA cassettes and smaller regions within the WTA cassettes is an unusual feature in the evolution of the  Bacillus subtilis genome and highlights the importance of horizontal transfer of gene cassettes through homologous recombination within  B. subtilis or other bacterial species.  Conclusions: Reduced sequence conservation of these WTA cassettes may indicate a modified function like the previously documented WTA ribitol/glycerol variation. An improved understanding of high-frequency recombination of gene cassettes has ramifications for synthetic biology and the use of  B. subtilis in industry.

Emerging Patterns in HIV-1 gp120 Variable Domains in Anatomical Tissues in the Absence of a Plasma Viral Load.

The HIV envelope protein contains five hypervariable domains (V1-V5) that are fundamental for cell entry. We contrasted modifications in the variable domains derived from a panel of 24 tissues from 7 subjects with no measurable plasma viral load (NPVL) to variable domains from 76 tissues from 15 subjects who had a detectable plasma viral load (PVL) at death. NPVL subject's V1 and V2 domains were usually highly length variable, whereas length variation in PVL sequences was more conserved. Longer V1s contained more charged residues, whereas longer V2s were more glycosylated. Structural analysis demonstrated V1/V2 charge, and N-site additions/subtractions were localized to the CD4 binding pocket. Diversified envelopes in tissues during therapy may represent a mechanism for HIV persistence in tissues, as binding pocket complexity is associated with HIV that may escape neutralization, whereas shorter envelopes are associated with increased infectivity. Further analysis of tissue-derived envelope sequences may enable better understanding of potential immunological approaches targeting the persistent HIV reservoir.

Brain-specific HIV Nef identified in multiple patients with neurological disease.

HIV-1 Nef is a flexible, multifunctional protein with several cellular targets that is required for pathogenicity of the virus. This protein maintains a high degree of genetic variation among intra- and inter-host isolates. HIV Nef is relevant to HIV-associated neurological diseases (HAND) in patients treated with combined antiretroviral therapy because of the protein's role in promoting survival and migration of infected brain macrophages. In this study, we analyzed 2020 HIV Nef sequences derived from 22 different tissues and 31 subjects using a novel computational approach. This approach combines statistical regression and evolved neural networks (ENNs) to classify brain sequences based on the physical and chemical characteristics of functional Nef domains. Based on training, testing, and validation data, the method successfully classified brain Nef sequences at 84.5% and provided informative features for further examination. These included physicochemical features associated with the Src-homology-3 binding domain, the Nef loop (including the AP-2 Binding region), and a cytokine-binding domain. Non-brain sequences from patients with HIV-associated neurological disease were frequently classified as brain, suggesting that the approach could indicate neurological risk using blood-derived virus or for the development of biomarkers for use in assay systems aimed at drug efficacy studies for the treatment of HIV-associated neurological diseases.

Structure-based analysis of Bacilli and plasmid dihydrofolate reductase evolution.

Dihydrofolate reductase (DHFR), a key enzyme in tetrahydrofolate-mediated biosynthetic pathways, has a structural motif known to be highly conserved over a wide range of organisms. Given its critical role in purine and amino acid synthesis, DHFR is a well established therapeutic target for treating a wide range of prokaryotic and eukaryotic infections as well as certain types of cancer. Here we present a structural-based computer analysis of bacterial (Bacilli) and plasmid DHFR evolution. We generated a structure-based sequence alignment using 7 wild-type DHFR x-ray crystal structures obtained from the RCSB Protein Data Bank and 350 chromosomal and plasmid homology models we generated from sequences obtained from the NCBI Protein Database. We used these alignments to compare active site and non-active site conservation in terms of amino acid residues, secondary structure and amino acid residue class. With respect to amino acid sequences and residue classes, active-site positions in both plasmid and chromosomal DHFR are significantly more conserved than non-active site positions. Secondary structure conservation was similar for active site and non-active site positions. Plasmid-encoded DHFR proteins have greater degree of sequence and residue class conservation, particularly in sequence positions associated with a network of concerted protein motions, than chromosomal-encoded DHFR proteins. These structure-based were used to build DHFR specific phylogenetic trees from which evidence for horizontal gene transfer was identified.

HIV-1 Evolutionary Patterns Associated with Metastatic Kaposi's Sarcoma during AIDS.

Kaposi's sarcoma (KS) in HIV-infected individuals can have a wide range of clinical outcomes, from indolent skin tumors to a life-threatening visceral cancer. KS tumors contain endothelial-related cells and inflammatory cells that may be HIV-infected. In this study we tested if HIV evolutionary patterns distinguish KS tumor relatedness and progression. Multisite autopsies from participants who died from HIV-AIDS with KS prior to the availability of antiretroviral therapy were identified at the AIDS and Cancer Specimen Resource (ACSR). Two patients (KS1 and KS2) died predominantly from non-KS-associated disease and KS3 died due to aggressive and metastatic KS within one month of diagnosis. Skin and visceral tumor and nontumor autopsy tissues were obtained (n = 12). Single genome sequencing was used to amplify HIV RNA and DNA, which was present in all tumors. Independent HIV tumor clades in phylogenies differentiated KS1 and KS2 from KS3, whose sequences were interrelated by both phylogeny and selection. HIV compartmentalization was confirmed in KS1 and KS2 tumors; however, in KS3, no compartmentalization was observed among sampled tissues. While the sample size is small, the HIV evolutionary patterns observed in all patients suggest an interplay between tumor cells and HIV-infected cells which provides a selective advantage and could promote KS progression.

HIV DNA Is Frequently Present within Pathologic Tissues Evaluated at Autopsy from Combined Antiretroviral Therapy-Treated Patients with Undetectable Viral Loads.

UNLABELLED: HIV infection treatment strategies have historically defined effectiveness through measuring patient plasma HIV RNA. While combined antiretroviral therapy (cART) can reduce plasma viral load (pVL) to undetectable levels, the degree that HIV is eliminated from other anatomical sites remains unclear. We investigated the HIV DNA levels in 229 varied autopsy tissues from 20 HIV-positive (HIV(+)) cART-treated study participants with low or undetectable plasma VL and cerebrospinal fluid (CSF) VL prior to death who were enrolled in the National Neurological AIDS Bank (NNAB) longitudinal study and autopsy cohort. Extensive medical histories were obtained for each participant. Autopsy specimens, including at least six brain and nonbrain tissues per participant, were reviewed by study pathologists. HIV DNA, measured in tissues by quantitative and droplet digital PCR, was identified in 48/87 brain tissues and 82/142 nonbrain tissues at levels >200 HIV copies/million cell equivalents. No participant was found to be completely free of tissue HIV. Parallel sequencing studies from some tissues recovered intact HIV DNA and RNA. Abnormal histological findings were identified in all participants, especially in brain, spleen, lung, lymph node, liver, aorta, and kidney. All brain tissues demonstrated some degree of pathology. Ninety-five percent of participants had some degree of atherosclerosis, and 75% of participants died with cancer. This study assists in characterizing the anatomical locations of HIV, in particular, macrophage-rich tissues, such as the central nervous system (CNS) and testis. Additional studies are needed to determine if the HIV recovered from tissues promotes the pathogenesis of inflammatory diseases, such as HIV-associated neurocognitive disorders, cancer, and atherosclerosis. IMPORTANCE: It is well-known that combined antiretroviral therapy (cART) can reduce plasma HIV to undetectable levels; however, cART cannot completely clear HIV infection. An ongoing question is, "Where is HIV hiding?" A well-studied HIV reservoir is "resting" T cells, which can be isolated from blood products and succumb to cART once activated. Less-studied reservoirs are anatomical tissue samples, which have unknown cART penetration, contain a comparably diverse spectrum of potentially HIV-infected immune cells, and are important since <2% of body lymphocytes actually reside in blood. We examined 229 varied autopsy specimens from 20 HIV(+) participants who died while on cART and identified that >50% of tissues were HIV infected. Additionally, we identified considerable pathology in participants' tissues, especially in brain, spleen, lung, lymph node, liver, aorta, and kidney. This study substantiates that tissue-associated HIV is present despite cART and can inform future studies into HIV persistence.

On the Physicochemical and Structural Modifications Associated with HIV-1 Subtype B Tropism Transition.

HIV-1 enters immune cells via binding the viral envelope to a host cell CD4 receptor, and then a secondary co-receptor, usually CCR5 (R5) or CXCR4 (X4), and some HIV can utilize both co-receptors (R5X4). Although a small set of amino-acid properties such as charge and sequence length applied to HIV-1 envelope V3 loop sequence data can be used to predict co-receptor usage, we sought to expand the fundamental understanding of the physiochemical basis of tropism by analyzing many, perhaps less obvious, amino-acid properties over a diverse array of HIV sequences. We examined 74 amino-acid physicochemical scales over 1,559 V3 loop sequences with biologically tested tropisms downloaded from the Los Alamos HIV sequence database. Linear regressions were then calculated for each feature relative to three tropism transitions (R5→X4; R5→R5X4; R5X4→X4). Independent correlations were rank ordered to determine informative features. A structural analysis of the V3 loop was performed to better interpret these findings relative to HIV tropism states. Similar structural changes are required for R5 and R5X4 to transition to X4, thus suggesting that R5 and R5X4 types are more similar than either phenotype is to X4. Overall, the analysis suggests a continuum of viral tropism that is only partially related to charge; in fact, the analysis suggests that charge modification may be primarily attributed to decreased R5 usage, and further structural changes, particularly those associated with ß-sheet structure, are likely required for full X4 usage.

Identification of dual-tropic HIV-1 using evolved neural networks.

Blocking the binding of the envelope HIV-1 protein to immune cells is a popular concept for development of anti-HIV therapeutics. R5 HIV-1 binds CCR5, X4 HIV-1 binds CXCR4, and dual-tropic HIV-1 can bind either coreceptor for cellular entry. R5 viruses are associated with early infection and over time can evolve to X4 viruses that are associated with immune failure. Dual-tropic HIV-1 is less studied; however, it represents functional antigenic intermediates during the transition of R5 to X4 viruses. Viral tropism is linked partly to the HIV-1 envelope V3 domain, where the amino acid sequence helps dictate the receptor a particular virus will target; however, using V3 sequence information to identify dual-tropic HIV-1 isolates has remained difficult. Our goal in this study was to elucidate features of dual-tropic HIV-1 isolates that assist in the biological understanding of dual-tropism and develop an approach for their detection. Over 1559 HIV-1 subtype B sequences with known tropisms were analyzed. Each sequence was represented by 73 structural, biochemical and regional features. These features were provided to an evolved neural network classifier and evaluated using balanced and unbalanced data sets. The study resolved R5X4 viruses from R5 with an accuracy of 81.8% and from X4 with an accuracy of 78.8%. The approach also identified a set of V3 features (hydrophobicity, structural and polarity) that are associated with tropism transitions. The ability to distinguish R5X4 isolates will improve computational tropism decisions for R5 vs. X4 and assist in HIV-1 research and drug development efforts.

MicroRNA dynamics during human embryonic stem cell differentiation to pancreatic endoderm.

MicroRNAs (miRNAs) are small non-coding RNAs that have emerged as critical regulators of human embryonic stem cell (hESC) pluripotency and differentiation. Despite the wealth of information about the role individual that miRNAs play in these two processes, there has yet to be a large-scale temporal analysis of the dynamics of miRNA expression as hESCs move from pluripotency into defined lineages. In this report, we used Next Generation Sequencing (NGS) to map temporal expression of miRNAs over ten 24-hour intervals as pluripotent cells were differentiated into pancreatic endoderm. Of the 2042 known human miRNAs, 694 had non-zero expression on all 11 days. Of these 694 miRNAs, 494 showed statistically significant changes in expression during differentiation. Clusters of miRNAs were identified, each displaying unique expression profiles distributed over multiple days. Selected miRNAs associated with pluripotency/differentiation (miR-302/367 and miR-371/372/373) and development/growth (miR-21, miR-25, miR-103, miR-9, and miR-92a) were found to have distinct expression profiles correlated with changes in media used to drive the differentiation process. Taken together, the clustering of miRNAs to identify expression dynamics that occur over longer periods of time (days vs. hours) provides unique insight into specific stages of differentiation. Major shifts in defined stages of hESC differentiation appear to be heavily dependent upon changes in external environmental factors, rather than intrinsic conditions in the cells.

Tracking the Emergence of Host-Specific Simian Immunodeficiency Virus env and nef Populations Reveals nef Early Adaptation and Convergent Evolution in Brain of Naturally Progressing Rhesus Macaques.

UNLABELLED: While a clear understanding of the events leading to successful establishment of host-specific viral populations and productive infection in the central nervous system (CNS) has not yet been reached, the simian immunodeficiency virus (SIV)-infected rhesus macaque provides a powerful model for the study of human immunodeficiency virus (HIV) intrahost evolution and neuropathogenesis. The evolution of the gp120 and nef genes, which encode two key proteins required for the establishment and maintenance of infection, was assessed in macaques that were intravenously inoculated with the same viral swarm and allowed to naturally progress to simian AIDS and potential SIV-associated encephalitis (SIVE). Longitudinal plasma samples and immune markers were monitored until terminal illness. Single-genome sequencing was employed to amplify full-length env through nef transcripts from plasma over time and from brain tissues at necropsy. nef sequences diverged from the founder virus faster than gp120 diverged. Host-specific sequence populations were detected in nef (~92 days) before they were detected in gp120 (~182 days). At necropsy, similar brain nef sequences were found in different macaques, indicating convergent evolution, while gp120 brain sequences remained largely host specific. Molecular clock and selection analyses showed weaker clock-like behavior and stronger selection pressure in nef than in gp120, with the strongest nef selection in the macaque with SIVE. Rapid nef diversification, occurring prior to gp120 diversification, indicates that early adaptation of nef in the new host is essential for successful infection. Moreover, the convergent evolution of nef sequences in the CNS suggests a significant role for nef in establishing neurotropic strains. IMPORTANCE: The SIV-infected rhesus macaque model closely resembles HIV-1 immunopathogenesis, neuropathogenesis, and disease progression in humans. Macaques were intravenously infected with identical viral swarms to investigate evolutionary patterns in the gp120 and nef genes leading to the emergence of host-specific viral populations and potentially linked to disease progression. Although each macaque exhibited unique immune profiles, macaque-specific nef sequences evolving under selection were consistently detected in plasma samples at 3 months postinfection, significantly earlier than in gp120 macaque-specific sequences. On the other hand, nef sequences in brain tissues, collected at necropsy of two animals with detectable infection in the central nervous system (CNS), revealed convergent evolution. The results not only indicate that early adaptation of nef in the new host may be essential for successful infection, but also suggest that specific nef variants may be required for SIV to efficiently invade CNS macrophages and/or enhance macrophage migration, resulting in HIV neuropathology.

Factors related to HIV-associated neurocognitive impairment differ with age.

Over 50% of HIV-infected (HIV+) persons are expected to be over age 50 by 2015. The pathogenic effects of HIV, particularly in cases of long-term infection, may intersect with those of age-related illnesses and prolonged exposure to combined antiretroviral therapy (cART). One potential outcome is an increased prevalence of neurocognitive impairment in older HIV+ individuals, as well as an altered presentation of HIV-associated neurocognitive disorders (HANDs). In this study, we employed stepwise regression to examine 24 features sometimes associated with HAND in 40 older (55-73 years of age) and 30 younger (32-50 years of age) HIV+, cART-treated participants without significant central nervous system confounds. The features most effective in generating a true assessment of the likelihood of HAND diagnosis differed between older and younger cohorts, with the younger cohort containing features associated with drug abuse that were correlated to HAND and the older cohort containing features that were associated with lipid disorders mildly associated with HAND. As the HIV-infected population grows and the demographics of the epidemic change, it is increasingly important to re-evaluate features associated with neurocognitive impairment. Here, we have identified features, routinely collected in primary care settings, that provide more accurate diagnostic value than a neurocognitive screening measure among younger and older HIV individuals.