IntroductionDuring the COVID-19 Delta variant surge, the CLAIRE cross-sectional study sampled saliva from 120 hospitalized patients, 116 of whom had a positive COVID-19 PCR test. Patients received antibiotics upon admission due to possible secondary bacterial infections, with patients at risk of sepsis receiving broad-spectrum antibiotics (BSA).MethodsThe saliva samples were analyzed with shotgun DNA metagenomics and respiratory RNA virome sequencing. Medical records for the period of hospitalization were obtained for all patients. Once hospitalization outcomes were known, patients were classified based on their COVID-19 disease severity and the antibiotics they received.ResultsOur study reveals that BSA regimens differentially impacted the human salivary microbiome and disease progression. 12 patients died and all of them received BSA. Significant associations were found between the composition of the COVID-19 saliva microbiome and BSA use, between SARS-CoV-2 genome coverage and severity of disease. We also found significant associations between the non-bacterial microbiome and severity of disease, with Candida albicans detected most frequently in critical patients. For patients who did not receive BSA before saliva sampling, our study suggests Staphylococcus aureus as a potential risk factor for sepsis.DiscussionOur results indicate that the course of the infection may be explained by both monitoring antibiotic treatment and profiling a patient’s salivary microbiome, establishing a compelling link between microbiome and the specific antibiotic type and timing of treatment. This approach can aid with emergency room triage and inpatient management but also requires a better understanding of and access to narrow-spectrum agents that target pathogenic bacteria.

LIM domain-binding 3 (LDB3) is a member of the Enigma family of PDZ–LIM proteins. LDB3 has been reported as a striated muscle-specific Z-band alternatively spliced protein that plays an important role in mechanosensory actin cytoskeleton remodeling. This study shows that LDB3 is broadly expressed in the central and peripheral nervous system of human and mouse. LDB3 is predominantly expressed in the adult stages compared to early development and at a significantly higher level in the spinal cord than in the brain. As in skeletal muscle and heart, LDB3 is extensively alternatively spliced in the neurons. Three novel splice isoforms were identified suggesting splicing-dependent regulation of LDB3 expression in the nervous system. Expression of LDB3 in the motor cortex, cerebellum, spinal motor neuron, peripheral nerve, and neuromuscular junction in addition to skeletal muscle indicates important roles for this PDZ–LIM family protein in motor planning and execution. Moreover, expression in the hippocampal neurons suggests roles for LDB3 in learning and memory. LDB3 interactors filamin C and myotilin are also expressed in the spinal motor neuron, nerve, and neuromuscular junction, thereby providing the basis for neurogenic manifestations in myopathies associated with mutations in these so-called muscle proteins.

Abstract During the COVID-19 Delta variant surge, the CLAIRE cross-sectional study sampled saliva from 120 hospitalized patients, 116 of them with a positive COVID-19 PCR test. Patients received antibiotics upon admission due to possible secondary bacterial infections, with patients at risk of sepsis receiving more potent antibiotics. The saliva samples were analyzed with shotgun DNA metagenomics and respiratory RNA virome sequencing. Medical records for the period of hospitalization were obtained for all patients. Once hospitalization outcome was known, patients were classified based on their COVID-19 disease severity and their sepsis-related antibiotic treatment. 12 patients passed away with all 12 having received sepsis treatment. Significant associations were found between the composition of the COVID-19 saliva microbiome and broad-spectrum antibiotic use, between SARS-CoV-2 genome coverage and severity of disease, and between the non-bacterial microbiome and severity of disease, with Candida albicans, the most significantly elevated fungus in critical patients. Our results suggest a use of the saliva microbiome as a biomarker that aids with emergency room triage and suggests the need to develop narrow-spectrum agents that target pathogenic bacteria.

Abstract LIM domain-binding 3 (LDB3) is a member of the Enigma family of PDZ-LIM proteins. LDB3 has been reported as a striated muscle-specific Z-band alternatively spliced protein that plays an important role in mechanosensory actin cytoskeleton remodeling. This study shows that LDB3 is broadly expressed in the central and peripheral nervous system of human and mouse. LDB3 is predominantly expressed in the adult stages compared to early development and at a significantly higher level in the spinal cord than in the brain. As in skeletal muscle and heart, LDB3 is extensively alternatively spliced in the neurons. Three novel splice isoforms were identified suggesting splicing-dependent regulation of LDB3 expression in the nervous system. Expression of LDB3 in the motor cortex, cerebellum, spinal motor neuron, peripheral nerve, and neuromuscular junction in addition to skeletal muscle indicates important roles for this PDZ-LIM family protein in motor planning and execution. Moreover, expression in the hippocampal neurons suggests roles for LDB3 in learning and memory. LDB3 interactors filamin C and myotilin are also expressed in the spinal motor neuron, nerve, and neuromuscular junction, thereby providing the basis for neurogenic manifestations in myopathies associated with mutations in these so-called muscle proteins.

Parasexuality contributes to diversity and adaptive evolution of haploid (monokaryotic) fungi. However, non-sexual genetic exchange mechanisms are not defined in dikaryotic fungi (containing two distinct haploid nuclei). Newly emerged strains of the wheat stem rust pathogen, Puccinia graminis f. sp. tritici (Pgt), such as Ug99, are a major threat to global food security. Here, we provide genomics-based evidence supporting that Ug99 arose by somatic hybridisation and nuclear exchange between dikaryons. Fully haplotype-resolved genome assembly and DNA proximity analysis reveal that Ug99 shares one haploid nucleus genotype with a much older African lineage of Pgt, with no recombination or chromosome reassortment. These findings indicate that nuclear exchange between dikaryotes can generate genetic diversity and facilitate the emergence of new lineages in asexual fungal populations.

Understanding HIV-1 persistence despite antiretroviral therapy (ART) is of paramount importance. Both single-genome sequencing (SGS) and integration site analysis (ISA) provide useful information regarding the structure of persistent HIV DNA populations; however, until recently, there was no way to link integration sites to their cognate proviral sequences. Here, we used multiple-displacement amplification (MDA) of cellular DNA diluted to a proviral endpoint to obtain full-length proviral sequences and their corresponding sites of integration. We applied this method to lymph node and peripheral blood mononuclear cells from 5 ART-treated donors to determine whether groups of identical subgenomic sequences in the 2 compartments are the result of clonal expansion of infected cells or a viral genetic bottleneck. We found that identical proviral sequences can result from both cellular expansion and viral genetic bottlenecks occurring prior to ART initiation and following ART failure. We identified an expanded T cell clone carrying an intact provirus that matched a variant previously detected by viral outgrowth assays and expanded clones with wild-type and drug-resistant defective proviruses. We also found 2 clones from 1 donor that carried identical proviruses except for nonoverlapping deletions, from which we could infer the sequence of the intact parental virus. Thus, MDA-SGS can be used for "viral reconstruction" to better understand intrapatient HIV-1 evolution and to determine the clonality and structure of proviruses within expanded clones, including those with drug-resistant mutations. Importantly, we demonstrate that identical sequences observed by standard SGS are not always sufficient to establish proviral clonality.

Background: Trichoplusia ni derived cell lines are commonly used to enable recombinant protein expression via baculovirus infection to generate materials approved for clinical use and in clinical trials. In order to develop systems biology and genome engineering tools to improve protein expression in this host, we performed de novo genome assembly of the Trichoplusia ni-derived cell line Tni-FNL. Methods: By integration of PacBio single-molecule sequencing, Bionano optical mapping, and 10X Genomics linked-reads data, we have produced a draft genome assembly of Tni-FNL. Results: Our assembly contains 280 scaffolds, with a N50 scaffold size of 2.3 Mb and a total length of 359 Mb. Annotation of the Tni-FNL genome resulted in 14,101 predicted genes and 93.2% of the predicted proteome contained recognizable protein domains. Ortholog searches within the superorder Holometabola provided further evidence of high accuracy and completeness of the Tni-FNL genome assembly. Conclusions: This first draft Tni-FNL genome assembly was enabled by complementary long-read technologies and represents a high-quality, well-annotated genome that provides novel insight into the complexity of this insect cell line and can serve as a reference for future large-scale genome engineering work in this and other similar recombinant protein production hosts.

Pacific Biosciences’ (PacBio) RS II sequencer, utilizing Single-Molecule, Real-Time (SMRT) technology, has revolutionized next-generation sequencing by providing an accurate long-read platform. PacBio single-molecule long reads have been used to delineate complex spliceoforms, detect mutations in highly homologous sequences, identify mRNA chimeras and chromosomal translocations, accurately haplotype phasing over multiple kilobase distances and aid in assembly of genomes with complex structural variation. The PacBio protocol for preparation of sequencing templates employs blunt-end hairpin adapter ligation, which enables a short turnaround time for sequence production. However, we have found a significant portion of sequencing yield contains chimeric reads resulting from blunt-end ligation of multiple template molecules to each other prior to adapter ligation. These artefactual fusion DNA sequences pose a major challenge to analysis and can lead to false-positive detection of fusion events. We assessed the frequency of artefactual fusion when using blunt-end adapter ligation and compared it to an alternative method using A/T overhang adapter ligation. The A/T overhang adapter ligation method showed a vast improvement in limiting artefactual fusion events and is now our recommended procedure for adapter ligation during PacBio library preparation.

ABSTRACT Oat crown rust, caused by the fungus Pucinnia coronata f. sp. avenae , is a devastating disease that impacts worldwide oat production. For much of its life cycle, P. coronata f. sp. avenae is dikaryotic, with two separate haploid nuclei that may vary in virulence genotype, highlighting the importance of understanding haplotype diversity in this species. We generated highly contiguous de novo genome assemblies of two P. coronata f. sp. avenae isolates, 12SD80 and 12NC29, from long-read sequences. In total, we assembled 603 primary contigs for 12SD80, for a total assembly length of 99.16 Mbp, and 777 primary contigs for 12NC29, for a total length of 105.25 Mbp; approximately 52% of each genome was assembled into alternate haplotypes. This revealed structural variation between haplotypes in each isolate equivalent to more than 2% of the genome size, in addition to about 260,000 and 380,000 heterozygous single-nucleotide polymorphisms in 12SD80 and 12NC29, respectively. Transcript-based annotation identified 26,796 and 28,801 coding sequences for isolates 12SD80 and 12NC29, respectively, including about 7,000 allele pairs in haplotype-phased regions. Furthermore, expression profiling revealed clusters of coexpressed secreted effector candidates, and the majority of orthologous effectors between isolates showed conservation of expression patterns. However, a small subset of orthologs showed divergence in expression, which may contribute to differences in virulence between 12SD80 and 12NC29. This study provides the first haplotype-phased reference genome for a dikaryotic rust fungus as a foundation for future studies into virulence mechanisms in P. coronata f. sp. avenae . IMPORTANCE Disease management strategies for oat crown rust are challenged by the rapid evolution of Puccinia coronata f. sp. avenae , which renders resistance genes in oat varieties ineffective. Despite the economic importance of understanding P. coronata f. sp. avenae , resources to study the molecular mechanisms underpinning pathogenicity and the emergence of new virulence traits are lacking. Such limitations are partly due to the obligate biotrophic lifestyle of P. coronata f. sp. avenae as well as the dikaryotic nature of the genome, features that are also shared with other important rust pathogens. This study reports the first release of a haplotype-phased genome assembly for a dikaryotic fungal species and demonstrates the amenability of using emerging technologies to investigate genetic diversity in populations of P. coronata f. sp. avenae .

Gene isoforms are commonly found in both prokaryotes and eukaryotes. Since each isoform may perform a specific function in response to changing environmental conditions, studying the dynamics of gene isoforms is important in understanding biological processes and disease conditions. However, genome-wide identification of gene isoforms is technically challenging due to the high degree of sequence identity among isoforms. Traditional targeted sequencing approach, involving Sanger sequencing of plasmid-cloned PCR products, has low throughput and is very tedious and time-consuming. Next-generation sequencing technologies such as Illumina and 454 achieve high throughput but their short read lengths are a critical barrier to accurate assembly of highly similar gene isoforms, and may result in ambiguities and false joining during sequence assembly. More recently, the third generation sequencer represented by the PacBio platform offers sufficient throughput and long reads covering the full length of typical genes, thus providing a potential to reliably profile gene isoforms. However, the PacBio long reads are error-prone and cannot be effectively analyzed by traditional assembly programs.

Targeted genetic correction of hematopoietic stem cells is applied to X-linked chronic granulomatous disease. Gene therapy with genetically modified human CD34+ hematopoietic stem and progenitor cells (HSPCs) may be safer using targeted integration (TI) of transgenes into a genomic 'safe harbor' site rather than random viral integration. We demonstrate that temporally optimized delivery of zinc finger nuclease mRNA via electroporation and adeno-associated virus (AAV) 6 delivery of donor constructs in human HSPCs approaches clinically relevant levels of TI into the AAVS1 safe harbor locus. Up to 58% Venus+ HSPCs with 6–16% human cell marking were observed following engraftment into mice. In HSPCs from patients with X-linked chronic granulomatous disease (X-CGD), caused by mutations in the gp91phox subunit of the NADPH oxidase, TI of a gp91phox transgene into AAVS1 resulted in ∼15% gp91phox expression and increased NADPH oxidase activity in ex vivo–derived neutrophils. In mice transplanted with corrected HSPCs, 4–11% of human cells in the bone marrow expressed gp91phox. This method for TI into AAVS1 may be broadly applicable to correction of other monogenic diseases.

Pneumocystis jirovecii is a major cause of life-threatening pneumonia in immunosuppressed patients including transplant recipients and those with HIV/AIDS, yet surprisingly little is known about the biology of this fungal pathogen. Here we report near complete genome assemblies for three Pneumocystis species that infect humans, rats and mice. Pneumocystis genomes are highly compact relative to other fungi, with substantial reductions of ribosomal RNA genes, transporters, transcription factors and many metabolic pathways, but contain expansions of surface proteins, especially a unique and complex surface glycoprotein superfamily, as well as proteases and RNA processing proteins. Unexpectedly, the key fungal cell wall components chitin and outer chain N-mannans are absent, based on genome content and experimental validation. Our findings suggest that Pneumocystis has developed unique mechanisms of adaptation to life exclusively in mammalian hosts, including dependence on the lungs for gas and nutrients and highly efficient strategies to escape both host innate and acquired immune defenses.

Abstract BACKGROUND: Helicobacter pylori, a gram-negative bacterium associated with a spectrum of benign and malignant gastric conditions, is one of the most genetically variable pathogens. Its genome encodes a large number of DNA methyltransferases targeting specific motif sequences of ∼4-15 bp. DNA base modifications epigenetically regulate gene expression and genome-wide methylation profiles (methylomes) have been hypothesized to be associated with virulence. A limited number of H. pylori methylomes have been published to date, precluding meaningful analyses of clinical correlates. METHODS: We sequenced 30 H. pylori clinical isolates from Mexico City using Pacific Biosciences’ Single Molecule Real-Time (SMRT) technology, a unique platform that simultaneously determines a wide range of DNA base modifications. The resulting methylome data were compared for the 15 strains from gastric cancer cases vs. the 15 from non-atrophic gastritis controls. RESULTS: All 30 strains exhibited high levels of DNA methylation throughout their genomes. Three types of base modification were detected: N6-methyladenine (m6A), N4-methylcytosine (m4C) and 5-methylcytosine (m5C). m6A was more common than m4C, and m5C was rare. Over 100 different motifs were detected, including many novel motifs not found in the three reference strains. For most methylated motifs, >95% of all occurrences in a given genome were methylated. Individual strains had between 14 and 28 different motifs. Some of the methylated motifs were conserved across all (e.g., Cm6ATG) or almost all (e.g., Gm6ACT) strains, but the majority were found in only one or two strains. Two conserved motifs, GTm6AC and GTNNm6AC, were significantly (P<0.05) more common in cancer cases than controls. CONCLUSIONS: The H. pylori methylome varies widely across strains, in concert with the high plasticity of its genome. Epigenetic modifications may modulate bacterial pathogenicity. These intriguing observations warrant replication with larger numbers and multiple populations in different geographic areas. The inferred novel methyltransferases remain to be identified and functionally characterized. This report greatly increases the number of H. pylori strains with defined methylomes, and represents the first effort to associate these features with clinical outcome. Availability of these genomic and epigenomic data will enhance understanding of H. pylori pathogenesis. Citation Format: M. Constanza Camargo, Javier Torres, Jay V. Solnick, Castle Raley, Dylan Storey, Marina Becker, Roberto Torres, Allison Weis, Carol Huang, Nguyet Kong, Lori M. Hanson, Yongmei Zhao, Xiongfong Chen, Bao Tran, Bart Weimer, Charles S. Rabkin. DNA methylation profiles of Helicobacter pylori strains from patients with gastric cancer and gastritis. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3438.

Development of HIV-1 drug resistance mutations (HDRMs) is one of the major reasons for the clinical failure of antiretroviral therapy. Treatment success rates can be improved by applying personalized anti-HIV regimens based on a patient's HDRM profile. However, the sensitivity and specificity of the HDRM profile is limited by the methods used for detection. Sanger-based sequencing technology has traditionally been used for determining HDRM profiles at the single nucleotide variant (SNV) level, but with a sensitivity of only ≥ 20% in the HIV population of a patient. Next Generation Sequencing (NGS) technologies offer greater detection sensitivity (~ 1%) and larger scope (hundreds of samples per run). However, NGS technologies produce reads that are too short to enable the detection of the physical linkages of individual SNVs across the haplotype of each HIV strain present. In this article, we demonstrate that the single-molecule long reads generated using the Third Generation Sequencer (TGS), PacBio RS II, along with the appropriate bioinformatics analysis method, can resolve the HDRM profile at a more advanced quasispecies level. The case studies on patients' HIV samples showed that the quasispecies view produced using the PacBio method offered greater detection sensitivity and was more comprehensive for understanding HDRM situations, which is complement to both Sanger and NGS technologies. In conclusion, the PacBio method, providing a promising new quasispecies level of HDRM profiling, may effect an important change in the field of HIV drug resistance research.

We have developed and validated an amplification-free method for generating DNA sequencing libraries from very low amounts of input DNA (500 picograms – 20 nanograms) for singlemolecule sequencing on the Pacific Biosciences (PacBio) RS II sequencer. The common challenge of high input requirements for single-molecule sequencing is overcome by using a carrier DNA in conjunction with optimized sequencing preparation conditions and re-use of the MagBead-bound complex. Here we describe how this method can be used to produce sequencing yields comparable to those generated from standard input amounts, but by using 1000-fold less starting material.

PacBio RS, a newly emerging third-generation DNA sequencing platform, is based on a real-time, single-molecule, nano-nitch sequencing technology that can generate very long reads (up to 20-kb) in contrast to the shorter reads produced by the first and second generation sequencing technologies.As a new platform, it is important to assess the sequencing error rate, as well as the quality control (QC) parameters associated with the PacBio sequence data.In this study, a mixture of 10 prior known, closely related DNA amplicons were sequenced using the PacBio RS sequencing platform.After aligning Circular Consensus Sequence (CCS) reads derived from the above sequencing experiment to the known reference sequences, we found that the median error rate was 2.5% without read QC, and improved to 1.3% with an SVM based multi-parameter QC method.In addition, a De Novo assembly was used as a downstream application to evaluate the effects of different QC approaches.This benchmark study indicates that even though CCS reads are post error-corrected it is still necessary to perform appropriate QC on CCS reads in order to produce successful downstream bioinformatics analytical results.