A long-read sequencing strategy with overlapping linkers on adjacent fragments (OLAF-Seq) for targeted resequencing and enrichment.

In this report, we present OLAF-Seq, a novel strategy to construct a long-read sequencing library such that adjacent fragments are linked with end-terminal duplications. We use the CRISPR-Cas9 nickase enzyme and a pool of multiple sgRNAs to perform non-random fragmentation of targeted long DNA molecules (> 300kb) into smaller library-sized fragments (about 20 kbp) in a manner so as to retain physical linkage information (up to 1000 bp) between adjacent fragments. DNA molecules targeted for fragmentation are preferentially ligated with adaptors for sequencing, so this method can enrich targeted regions while taking advantage of the long-read sequencing platforms. This enables the sequencing of target regions with significantly lower total coverage, and the genome sequence within linker regions provides information for assembly and phasing. We demonstrated the validity and efficacy of the method first using phage and then by sequencing a panel of 100 full-length cancer-related genes (including both exons and introns) in the human genome. When the designed linkers contained heterozygous genetic variants, long haplotypes could be established. This sequencing strategy can be readily applied in both PacBio and Oxford Nanopore platforms for both long and short genes with an easy protocol. This economically viable approach is useful for targeted enrichment of hundreds of target genomic regions and where long no-gap contigs need deep sequencing.

Long-term Sudan virus Ebola survivors maintain multiple antiviral defense mechanisms.

BACKGROUND: The critical issues of sustained memory immunity following ebolavirus disease among long-term survivors (EVD) are still unclear. METHODS: Here, we examine virus-specific immune and inflammatory responses in 12 Sudan virus (SUDV) long-term survivors from Uganda's 2000-1 Gulu outbreak, 15 years after recovery following in vitro challenge. Total RNA from isolated SUDV-stimulated and unstimulated PBMCs was extracted and analyzed. Matched serum samples were also collected to determine SUDV IgG levels and functionality. RESULTS: We detected persistent humoral (58%, 7 of 12) and cellular (33%, 4 of 12) immune responses in SUDV long-term survivors and identified critical molecular mechanisms of innate and adaptive immunity. Gene expression in immune pathways, the IFN signaling system, antiviral defense response, and activation and regulation of T- and B-cell responses were observed. SUDV long-term survivors also maintained robust virus-specific IgG antibodies capable of polyfunctional responses, including neutralizing and innate Fc effector functions. CONCLUSIONS: Data integration identified significant correlations among humoral and cellular immune responses and pinpointed a specific innate and adaptive gene expression signature associated with long-lasting immunity. This could help identify natural and vaccine correlates of protection against ebolavirus disease.

Phylogenomics of nontuberculous mycobacteria respiratory infections in people with cystic fibrosis.

Nontuberculous mycobacteria (NTM) can cause severe pulmonary disease in people with cystic fibrosis (pwCF). These infections present unique challenges for diagnosis and treatment, prompting a recent interest in understanding NTM transmission and pathogenesis during chronic infection. Major gaps remain in our knowledge regarding basic pathogenesis, immune evasion strategies, population dynamics, recombination potential, and the evolutionary implications of host and antibiotic pressures of long-term NTM infections in pwCF. Phylogenomic techniques have emerged as an important tool for tracking global patterns of transmission and are beginning to be used to ask fundamental biological questions about adaptation to the host during pathogenesis. In this review, we discuss the burden of NTM lung disease (NTM-LD), highlight the use of phylogenomics in NTM research, and address the clinical implications associated with these studies.

Genetic Adaptation and Acquisition of Macrolide Resistance in Haemophilus spp. during Persistent Respiratory Tract Colonization in Chronic Obstructive Pulmonary Disease (COPD) Patients Receiving Long-Term Azithromycin Treatment.

Patients with chronic obstructive pulmonary disease (COPD) benefit from the immunomodulatory effect of azithromycin, but long-term administration may alter colonizing bacteria. Our goal was to identify changes in Haemophilus influenzae and Haemophilus parainfluenzae during azithromycin treatment. Fifteen patients were followed while receiving prolonged azithromycin treatment (Hospital Universitari de Bellvitge, Spain). Four patients (P02, P08, P11, and P13) were persistently colonized by H. influenzae for at least 3 months and two (P04 and P11) by H. parainfluenzae. Isolates from these patients (53 H. influenzae and 18 H. parainfluenzae) were included to identify, by whole-genome sequencing, antimicrobial resistance changes and genetic variation accumulated during persistent colonization. All persistent lineages isolated before treatment were azithromycin-susceptible but developed resistance within the first months, apart from those belonging to P02, who discontinued the treatment. H. influenzae isolates from P08-ST107 acquired mutations in 23S rRNA, and those from P11-ST2480 and P13-ST165 had changes in L4 and L22. In H. parainfluenzae, P04 persistent isolates acquired changes in rlmC, and P11 carried genes encoding MefE/MsrD efflux pumps in an integrative conjugative element, which was also identified in H. influenzae P11-ST147. Other genetic variation occurred in genes associated with cell wall and inorganic ion metabolism. Persistent H. influenzae strains all showed changes in licA and hgpB genes. Other genes (lex1, lic3A, hgpC, and fadL) had variation in multiple lineages. Furthermore, persistent strains showed loss, acquisition, or genetic changes in prophage-associated regions. Long-term azithromycin therapy results in macrolide resistance, as well as genetic changes that likely favor bacterial adaptation during persistent respiratory colonization. IMPORTANCE The immunomodulatory properties of azithromycin reduce the frequency of exacerbations and improve the quality of life of COPD patients. However, long-term administration may alter the respiratory microbiota, such as Haemophilus influenzae, an opportunistic respiratory colonizing bacteria that play an important role in exacerbations. This study contributes to a better understanding of COPD progression by characterizing the clinical evolution of H. influenzae in a cohort of patients with prolonged azithromycin treatment. The emergence of macrolide resistance during the first months, combined with the role of Haemophilus parainfluenzae as a reservoir and source of resistance dissemination, is a cause for concern that may lead to therapeutic failure. Furthermore, genetic variations in cell wall and inorganic ion metabolism coding genes likely favor bacterial adaptation to host selective pressures. Therefore, the bacterial pathoadaptive evolution in these severe COPD patients raise our awareness of the possible spread of macrolide resistance and selection of host-adapted clones.

Transition of Serotype 35B Pneumococci From Commensal to Prevalent Virulent Strain in Children.

In our community-based prospective cohort study in young children, we observed a significant increase in pneumococcal serotype 35B nasopharyngeal (NP) commensal colonization during the 2011-2014 timeframe, but these strains were not associated with disease. Beginning in 2015 and continuing through to the present, the serotype 35B virulence changed, and it became the dominant bacteria isolated and associated with pneumococcal acute otitis-media (AOM) in our cohort. We performed comparative analyses of 250 35B isolates obtained from 140 children collected between 2006 and 2019. Changes in prevalence, clonal-complex composition, and antibiotic resistance were analyzed. Seventy-two (29%) of 35B isolates underwent whole-genome sequencing to investigate genomic changes associated with the shift in virulence that resulted in increased rates of 35B-associated AOM disease. 35B strains that were commensals and AOM disease-causing were mainly associated with sequence type (ST) 558. Antibiotic concentrations of ß-lactams and ofloxacin necessary to inhibit growth of 35B strains rose significantly (2006-2019) (p<0.005). However, only isolates from the 35B/ST558 showed significant increases in MIC50 of penicillin and ofloxacin between the years 2006-2014 and 2015-2019 (p=0.007 and p<0.0001). One hundred thirty-eight SNPs located in 34 different genes were significantly associated with post-2015 strains. SNPs were found in nrdG (metal binding, 10%); metP and metN (ABC transporter, 9%); corA (Mg2+ transporter, 6%); priA (DNA replication, 5%); and on the enzymic gene ldcB (LD-carboxypeptidase, 3%). Pneumococcal serotype 35B strains was a common NP commensal during 2010-2014. In 2015, a shift in increasing number of AOM cases occurred in young children caused by 35B, that was associated with changes in genetic composition and antibiotic susceptibility.

Learning, visualizing and exploring 16S rRNA structure using an attention-based deep neural network.

Recurrent neural networks with memory and attention mechanisms are widely used in natural language processing because they can capture short and long term sequential information for diverse tasks. We propose an integrated deep learning model for microbial DNA sequence data, which exploits convolutional neural networks, recurrent neural networks, and attention mechanisms to predict taxonomic classifications and sample-associated attributes, such as the relationship between the microbiome and host phenotype, on the read/sequence level. In this paper, we develop this novel deep learning approach and evaluate its application to amplicon sequences. We apply our approach to short DNA reads and full sequences of 16S ribosomal RNA (rRNA) marker genes, which identify the heterogeneity of a microbial community sample. We demonstrate that our implementation of a novel attention-based deep network architecture, Read2Pheno, achieves read-level phenotypic prediction. Training Read2Pheno models will encode sequences (reads) into dense, meaningful representations: learned embedded vectors output from the intermediate layer of the network model, which can provide biological insight when visualized. The attention layer of Read2Pheno models can also automatically identify nucleotide regions in reads/sequences which are particularly informative for classification. As such, this novel approach can avoid pre/post-processing and manual interpretation required with conventional approaches to microbiome sequence classification. We further show, as proof-of-concept, that aggregating read-level information can robustly predict microbial community properties, host phenotype, and taxonomic classification, with performance at least comparable to conventional approaches. An implementation of the attention-based deep learning network is available at https://github.com/EESI/sequence_attention (a python package) and https://github.com/EESI/seq2att (a command line tool).

Species-Level Profiling of Ixodes pacificus Bacterial Microbiomes Reveals High Variability Across Short Spatial Scales at Different Taxonomic Resolutions.

Background and Aims: Outbreaks of severe and chronic tick-borne diseases (TBDs) are on the rise. This is through the transmission of infectious disease agents to humans during tick feeding. The transmission rate and extent of microbial exchange, however, vary based on the tick microbiome composition. While select microbes are determined to be members of the normal tick microbiome and others are clearly recognized mammalian and/or avian pathogens, the status of many other members of the tick microbiota with respect to human and alternate host pathogenesis remains unclear. Moreover, the species-level 16S microbiome of prominent TBD vectors, including Ixodes pacificus, have not been extensively studied. To elucidate the I. pacificus microbiome composition, we performed a pan-domain species-specific characterization of the bacterial microbiome on adult I. pacificus ticks collected from two regional parks within Western California. Our methods provide for characterizing nuances within cohort microbiomes and their relationships to geo-locale of origin, surrounding fauna, and prevalences of known and suspected pathogens in relation to current TBD epidemiological zones. Methods: Ninety-two adult I. pacificus bacterial microbiomes were characterized using a high-fidelity, pan-domain, species-specific, full-length 16S rRNA amplification method using circular consensus sequencing performed on the Pacific Biosciences Sequel platform. Data analyses were performed with the MCSMRT data analysis package and database. Results: The species-specific I. pacificus microbiome composition illustrates a complex assortment of microflora, including over 900 eubacterial species with high taxonomic diversity, which was revealed to vary by sex and geo-locale, though the use of full-length 16S gene sequencing. The TBD-associated pathogens, such as Borrelia burgdorferi, Anaplasma phagocytophilum, and Rickettsia monacensis, were identified along with a host of bacteria previously unassociated with ticks. Conclusion: Species-level taxonomic classification of the I. pacificus microbiome revealed that full-length bacterial 16S gene sequencing is required for the granularity to elucidate the microbial diversity within and among ticks based on geo-locale.

Genome-wide analysis of DNA uptake across the outer membrane of naturally competent Haemophilus influenzae.

The genomes of naturally competent Pasteurellaceae and Neisseriaceae have many short uptake sequences (USS), which allow them to distinguish self-DNA from foreign DNA. To fully characterize this preference we developed genome-wide maps of DNA uptake using both a sequence-based computational model and genomic DNA that had been sequenced after uptake by and recovery from competent Haemophilus influenzae cells. When DNA fragments were shorter than the average USS spacing of ∼1,000 bp, sharp peaks of uptake were centered at USS and separated by valleys with 1000-fold lower uptake. Long DNA fragments (1.5-17 kb) gave much less variation, with 90% of positions having uptake within 2-fold of the mean. All detectable uptake biases arose from sequences that fit the USS uptake motif. Simulated competition predicted that, in its respiratory tract environment, H. influenzae will efficiently take up its own DNA even when human DNA is present in 100-fold excess.

Maf1 limits RNA polymerase III-directed transcription to preserve genomic integrity and extend lifespan.

A key to longevity assurance is the nutrient-sensing mTOR pathway. Inhibition of mTOR extends lifespan in a variety of organisms. However, the downstream effectors of the mTOR pathway for lifespan regulation are elusive. In a recent report, we described the role of Maf1 as a critical lifespan regulator downstream of the mTOR pathway in fission yeast. Maf1 is the master negative regulator of RNA polymerase III-directed transcription (e.g. tRNAs and 5S rRNAs) and is regulated by mTOR-mediated phosphorylation. We demonstrated that Maf1 is required for lifespan extension under calorie restriction or when mTOR is inhibited. We also showed that Maf1 prevents DNA damage at tRNA genes, which appears to contribute to lifespan maintenance by Maf1. Here we highlight these observations and present additional results to discuss the role of the mTOR-Maf1-Pol III axis in promoting genomic integrity in the face of DNA replication-transcription conflicts in order to maintain normal lifespan.

Maf1-dependent transcriptional regulation of tRNAs prevents genomic instability and is associated with extended lifespan.

Maf1 is the master repressor of RNA polymerase III responsible for transcription of tRNAs and 5S rRNAs. Maf1 is negatively regulated via phosphorylation by the mTOR pathway, which governs protein synthesis, growth control, and lifespan regulation in response to nutrient availability. Inhibiting the mTOR pathway extends lifespan in various organisms. However, the downstream effectors for the regulation of cell homeostasis that are critical to lifespan extension remain elusive. Here we show that fission yeast Maf1 is required for lifespan extension. Maf1's function in tRNA repression is inhibited by mTOR-dependent phosphorylation, whereas Maf1 is activated via dephosphorylation by protein phosphatase complexes, PP4 and PP2A. Mutational analysis reveals that Maf1 phosphorylation status influences lifespan, which is correlated with elevated tRNA and protein synthesis levels in maf1∆ cells. However, mTOR downregulation, which negates protein synthesis, fails to rescue the short lifespan of maf1∆ cells, suggesting that elevated protein synthesis is not a cause of lifespan shortening in maf1∆ cells. Interestingly, maf1∆ cells accumulate DNA damage represented by formation of Rad52 DNA damage foci and Rad52 recruitment at tRNA genes. Loss of the Rad52 DNA repair protein further exacerbates the shortened lifespan of maf1∆ cells. Strikingly, PP4 deletion alleviates DNA damage and rescues the short lifespan of maf1∆ cells even though tRNA synthesis is increased in this condition, suggesting that elevated DNA damage is the major cause of lifespan shortening in maf1∆ cells. We propose that Maf1-dependent inhibition of tRNA synthesis controls fission yeast lifespan by preventing genomic instability that arises at tRNA genes.

Exploring thematic structure and predicted functionality of 16S rRNA amplicon data.

Analysis of microbiome data involves identifying co-occurring groups of taxa associated with sample features of interest (e.g., disease state). Elucidating such relations is often difficult as microbiome data are compositional, sparse, and have high dimensionality. Also, the configuration of co-occurring taxa may represent overlapping subcommunities that contribute to sample characteristics such as host status. Preserving the configuration of co-occurring microbes rather than detecting specific indicator species is more likely to facilitate biologically meaningful interpretations. Additionally, analyses that use taxonomic relative abundances to predict the abundances of different gene functions aggregate predicted functional profiles across taxa. This precludes straightforward identification of predicted functional components associated with subsets of co-occurring taxa. We provide an approach to explore co-occurring taxa using "topics" generated via a topic model and link these topics to specific sample features (e.g., disease state). Rather than inferring predicted functional content based on overall taxonomic relative abundances, we instead focus on inference of functional content within topics, which we parse by estimating interactions between topics and pathways through a multilevel, fully Bayesian regression model. We apply our methods to three publicly available 16S amplicon sequencing datasets: an inflammatory bowel disease dataset, an oral cancer dataset, and a time-series dataset. Using our topic model approach to uncover latent structure in 16S rRNA amplicon surveys, investigators can (1) capture groups of co-occurring taxa termed topics; (2) uncover within-topic functional potential; (3) link taxa co-occurrence, gene function, and environmental/host features; and (4) explore the way in which sets of co-occurring taxa behave and evolve over time. These methods have been implemented in a freely available R package: https://cran.r-project.org/package=themetagenomics, https://github.com/EESI/themetagenomics.

Novel gRNA design pipeline to develop broad-spectrum CRISPR/Cas9 gRNAs for safe targeting of the HIV-1 quasispecies in patients.

The CRISPR/Cas9 system has been proposed as a cure strategy for HIV. However, few published guide RNAs (gRNAs) are predicted to cleave the majority of HIV-1 viral quasispecies (vQS) observed within and among patients. We report the design of a novel pipeline to identify gRNAs that target HIV across a large number of infected individuals. Next generation sequencing (NGS) of LTRs from 269 HIV-1-infected samples in the Drexel CARES Cohort was used to select gRNAs with predicted broad-spectrum activity. In silico, D-LTR-P4-227913 (package of the top 4 gRNAs) accounted for all detectable genetic variation within the vQS of the 269 samples and the Los Alamos National Laboratory HIV database. In silico secondary structure analyses from NGS indicated extensive TAR stem-loop malformations predicted to inactivate proviral transcription, which was confirmed by reduced viral gene expression in TZM-bl or P4R5 cells. Similarly, a high sensitivity in vitro CRISPR/Cas9 cleavage assay showed that the top-ranked gRNA was the most effective at cleaving patient-derived HIV-1 LTRs from five patients. Furthermore, the D-LTR-P4-227913 was predicted to cleave a median of 96.1% of patient-derived sequences from other HIV subtypes. These results demonstrate that the gRNAs possess broad-spectrum cutting activity and could contribute to an HIV cure.

Lung transcriptional unresponsiveness and loss of early influenza virus control in infected neonates is prevented by intranasal Lactobacillus rhamnosus GG.

Respiratory viral infections contribute substantially to global infant losses and disproportionately affect preterm neonates. Using our previously established neonatal murine model of influenza infection, we demonstrate that three-day old mice are exceptionally sensitive to influenza virus infection and exhibit high mortality and viral load. Intranasal pre- and post-treatment of neonatal mice with Lactobacillus rhamnosus GG (LGG), an immune modulator in respiratory viral infection of adult mice and human preterm neonates, considerably improves neonatal mice survival after influenza virus infection. We determine that both live and heat-killed intranasal LGG are equally efficacious in protection of neonates. Early in influenza infection, neonatal transcriptional responses in the lung are delayed compared to adults. These responses increase by 24 hours post-infection, demonstrating a delay in the kinetics of the neonatal anti-viral response. LGG pretreatment improves immune gene transcriptional responses during early infection and specifically upregulates type I IFN pathways. This is critical for protection, as neonatal mice intranasally pre-treated with IFNß before influenza virus infection are also protected. Using transgenic mice, we demonstrate that the protective effect of LGG is mediated through a MyD88-dependent mechanism, specifically via TLR4. LGG can improve both early control of virus and transcriptional responsiveness and could serve as a simple and safe intervention to protect neonates.

Species-level bacterial community profiling of the healthy sinonasal microbiome using Pacific Biosciences sequencing of full-length 16S rRNA genes.

BACKGROUND: Pan-bacterial 16S rRNA microbiome surveys performed with massively parallel DNA sequencing technologies have transformed community microbiological studies. Current 16S profiling methods, however, fail to provide sufficient taxonomic resolution and accuracy to adequately perform species-level associative studies for specific conditions. This is due to the amplification and sequencing of only short 16S rRNA gene regions, typically providing for only family- or genus-level taxonomy. Moreover, sequencing errors often inflate the number of taxa present. Pacific Biosciences' (PacBio's) long-read technology in particular suffers from high error rates per base. Herein, we present a microbiome analysis pipeline that takes advantage of PacBio circular consensus sequencing (CCS) technology to sequence and error correct full-length bacterial 16S rRNA genes, which provides high-fidelity species-level microbiome data. RESULTS: Analysis of a mock community with 20 bacterial species demonstrated 100% specificity and sensitivity with regard to taxonomic classification. Examination of a 250-plus species mock community demonstrated correct species-level classification of > 90% of taxa, and relative abundances were accurately captured. The majority of the remaining taxa were demonstrated to be multiply, incorrectly, or incompletely classified. Using this methodology, we examined the microgeographic variation present among the microbiomes of six sinonasal sites, by both swab and biopsy, from the anterior nasal cavity to the sphenoid sinus from 12 subjects undergoing trans-sphenoidal hypophysectomy. We found greater variation among subjects than among sites within a subject, although significant within-individual differences were also observed. Propiniobacterium acnes (recently renamed Cutibacterium acnes) was the predominant species throughout, but was found at distinct relative abundances by site. CONCLUSIONS: Our microbial composition analysis pipeline for single-molecule real-time 16S rRNA gene sequencing (MCSMRT, https://github.com/jpearl01/mcsmrt ) overcomes deficits of standard marker gene-based microbiome analyses by using CCS of entire 16S rRNA genes to provide increased taxonomic and phylogenetic resolution. Extensions of this approach to other marker genes could help refine taxonomic assignments of microbial species and improve reference databases, as well as strengthen the specificity of associations between microbial communities and dysbiotic states.

Antagonistic Pleiotropy in the Bifunctional Surface Protein FadL (OmpP1) during Adaptation of Haemophilus influenzae to Chronic Lung Infection Associated with Chronic Obstructive Pulmonary Disease.

Tracking bacterial evolution during chronic infection provides insights into how host selection pressures shape bacterial genomes. The human-restricted opportunistic pathogen nontypeable Haemophilus influenzae (NTHi) infects the lower airways of patients suffering chronic obstructive pulmonary disease (COPD) and contributes to disease progression. To identify bacterial genetic variation associated with bacterial adaptation to the COPD lung, we sequenced the genomes of 92 isolates collected from the sputum of 13 COPD patients over 1 to 9 years. Individuals were colonized by distinct clonal types (CTs) over time, but the same CT was often reisolated at a later time or found in different patients. Although genomes from the same CT were nearly identical, intra-CT variation due to mutation and recombination occurred. Recurrent mutations in several genes were likely involved in COPD lung adaptation. Notably, nearly a third of CTs were polymorphic for null alleles of ompP1 (also called fadL), which encodes a bifunctional membrane protein that both binds the human carcinoembryonic antigen-related cell adhesion molecule 1 (hCEACAM1) receptor and imports long-chain fatty acids (LCFAs). Our computational studies provide plausible three-dimensional models for FadL's interaction with hCEACAM1 and LCFA binding. We show that recurrent fadL mutations are likely a case of antagonistic pleiotropy, since loss of FadL reduces NTHi's ability to infect epithelia but also increases its resistance to bactericidal LCFAs enriched within the COPD lung. Supporting this interpretation, truncated fadL alleles are common in publicly available NTHi genomes isolated from the lower airway tract but rare in others. These results shed light on molecular mechanisms of bacterial pathoadaptation and guide future research toward developing novel COPD therapeutics.IMPORTANCE Nontypeable Haemophilus influenzae is an important pathogen in patients with chronic obstructive pulmonary disease (COPD). To elucidate the bacterial pathways undergoing in vivo evolutionary adaptation, we compared bacterial genomes collected over time from 13 COPD patients and identified recurrent genetic changes arising in independent bacterial lineages colonizing different patients. Besides finding changes in phase-variable genes, we found recurrent loss-of-function mutations in the ompP1 (fadL) gene. We show that loss of OmpP1/FadL function reduces this bacterium's ability to infect cells via the hCEACAM1 epithelial receptor but also increases its resistance to bactericidal fatty acids enriched within the COPD lung, suggesting a case of antagonistic pleiotropy that restricts ΔfadL strains' niche. These results show how H. influenzae adapts to host-generated inflammatory mediators in the COPD airways.

Broad-Spectrum and Personalized Guide RNAs for CRISPR/Cas9 HIV-1 Therapeutics.

The clustered regularly interspaced short palindromic repeats (CRISPR)-associated Cas9 system has been used to excise the HIV-1 proviral genome from latently infected cells, potentially offering a cure for HIV-infected patients. Recent studies have shown that most published HIV-1 guide RNAs (gRNAs) do not account for the diverse viral quasispecies within or among patients, which continue to diversify with time even in long-term antiretroviral therapy (ART)-suppressed patients. Given this observation, proviral genomes were deep sequenced from 23 HIV-1-infected patients in the Drexel Medicine CNS AIDS Research and Eradication Study cohort at two different visits. Based on the spectrum of integrated proviral DNA polymorphisms observed, three gRNA design strategies were explored: based on the patient's own HIV-1 sequences (personalized), based on consensus sequences from a large sample of patients [broad-spectrum (BS)], or a combination of both approaches. Using a bioinformatic algorithm, the personalized gRNA design was predicted to cut 46 of 48 patient samples at 90% efficiency, whereas the top 4 BS gRNAs (BS4) were predicted to excise provirus from 44 of 48 patient samples with 90% efficiency. Using a mixed design with the top three BS gRNAs plus one personalized gRNA (BS3 + PS1) resulted in predicted excision of provirus from 45 of 48 patient samples with 90% efficiency. In summary, these studies used an algorithmic design strategy to identify potential BS gRNAs to target a spectrum of HIV-1 long teriminal repeat (LTR) quasispecies for use with a small HIV-1-infected population. This approach should advance CRISPR/Cas9 excision technology taking into account the extensive molecular heterogeneity of HIV-1 that persists in situ after prolonged ART.

Designing broad-spectrum anti-HIV-1 gRNAs to target patient-derived variants.

Clustered regularly interspaced short palindromic repeats (CRISPR) CRISPR-associated protein 9 (Cas9), including specific guide RNAs (gRNAs), can excise integrated human immunodeficiency virus type 1 (HIV-1) provirus from host chromosomes. To date, anti-HIV-1 gRNAs have been designed to account for off-target activity, however, they seldom account for genetic variation in the HIV-1 genome within and between patients, which will be crucial for therapeutic application of this technology. This analysis tests the ability of published anti-HIV-1 gRNAs to cleave publicly available patient-derived HIV-1 sequences to inform gRNA design and provides basic computational tools to researchers in the field.

Phenotypic diversity and genotypic flexibility of Burkholderia cenocepacia during long-term chronic infection of cystic fibrosis lungs.

Chronic bacterial infections of the lung are the leading cause of morbidity and mortality in cystic fibrosis patients. Tracking bacterial evolution during chronic infections can provide insights into how host selection pressures-including immune responses and therapeutic interventions-shape bacterial genomes. We carried out genomic and phenotypic analyses of 215 serially collected Burkholderia cenocepacia isolates from 16 cystic fibrosis patients, spanning a period of 2-20 yr and a broad range of epidemic lineages. Systematic phenotypic tests identified longitudinal bacterial series that manifested progressive changes in liquid media growth, motility, biofilm formation, and acute insect virulence, but not in mucoidy. The results suggest that distinct lineages follow distinct evolutionary trajectories during lung infection. Pan-genome analysis identified 10,110 homologous gene clusters present only in a subset of strains, including genes restricted to different molecular types. Our phylogenetic analysis based on 2148 orthologous gene clusters from all isolates is consistent with patient-specific clades. This suggests that initial colonization of patients was likely by individual strains, followed by subsequent diversification. Evidence of clonal lineages shared by some patients was observed, suggesting inter-patient transmission. We observed recurrent gene losses in multiple independent longitudinal series, including complete loss of Chromosome III and deletions on other chromosomes. Recurrently observed loss-of-function mutations were associated with decreases in motility and biofilm formation. Together, our study provides the first comprehensive genome-phenome analyses of B. cenocepacia infection in cystic fibrosis lungs and serves as a valuable resource for understanding the genomic and phenotypic underpinnings of bacterial evolution.

Complete Genome Sequence of Aggregatibacter actinomycetemcomitans Strain IDH781.

We report here the complete genomic sequence and methylome of Aggregatibacter actinomycetemcomitans strain IDH781. This rough strain is used extensively as a model organism to characterize localized aggressive periodontitis pathogenesis, the basic biology and oral cavity colonization of A. actinomycetemcomitans, and its interactions with other members of the oral microbiome.

Mitochondrial Haplogroup Influences Motor Function in Long-Term HIV-1-Infected Individuals.

Evolutionary divergence of the mitochondrial genome has given rise to distinct haplogroups. These haplogroups have arisen in specific geographical locations and are responsible for subtle functional changes in the mitochondria that may provide an evolutionary advantage in a given environment. Based on these functional differences, haplogroups could define disease susceptibility in chronic settings. In this study, we undertook a detailed neuropsychological analysis of a cohort of long-term HIV-1-infected individuals in conjunction with sequencing of their mitochondrial genomes. Stepwise regression analysis showed that the best model for predicting both working memory and declarative memory were age and years since diagnosis. In contrast, years since diagnosis and sub-haplogroup were significantly predictive of psychomotor speed. Consistent with this, patients with haplogroup L3e obtained better scores on psychomotor speed and dexterity tasks when compared to the remainder of the cohort, suggesting that this haplogroup provides a protective advantage when faced with the combined stress of HIV-1 infection and long-term antiretroviral therapies. Differential performance on declarative memory tasks was noted for individuals with other sub-L haplogroups, but these differences were not as robust as the association between L3e and psychomotor speed and dexterity tasks. This work provides evidence that mitochondrial haplogroup is related to neuropsychological test performance among patients in chronic disease settings such as HIV-1 infection.