Scientific Publications

Abstract +

Microbial communities might include distinct lineages of closely related organisms that complicate metagenomic assembly and prevent the generation of complete metagenome-assembled genomes (MAGs). Here we show that deep sequencing using long (HiFi) reads combined with Hi-C binning can address this challenge even for complex microbial communities. Using existing methods, we sequenced the sheep fecal metagenome and identified 428 MAGs with more than 90% completeness, including 44 MAGs in single circular contigs. To resolve closely related strains (lineages), we developed MAGPhase, which separates lineages of related organisms by discriminating variant haplotypes across hundreds of kilobases of genomic sequence. MAGPhase identified 220 lineage-resolved MAGs in our dataset. The ability to resolve closely related microbes in complex microbial communities improves the identification of biosynthetic gene clusters and the precision of assigning mobile genetic elements to host genomes. We identified 1,400 complete and 350 partial biosynthetic gene clusters, most of which are novel, as well as 424 (298) potential host–viral (host–plasmid) associations using Hi-C data.

Abstract +

Currently, protein-coding de novo variants and large copy number variants have been identified as important for ∼30% of individuals with autism. One approach to identify relevant variation in individuals who lack these types of events is by utilizing newer genomic technologies. In this study, highly accurate PacBio HiFi long-read sequencing was applied to a family with autism, treatment-refractory epilepsy, cognitive impairment, and mild dysmorphic features (two affected female full siblings, parents, and one unaffected sibling) with no known clinical variant. From our long-read sequencing data, a de novo missense variant in the KCNC2 gene (encodes Kv3.2 protein) was identified in both affected children. This variant was phased to the paternal chromosome of origin and is likely a germline mosaic. In silico assessment of the variant revealed it was in the top 0.05% of all conserved bases in the genome, and was predicted damaging by Polyphen2, MutationTaster, and SIFT. It was not present in any controls from public genome databases nor in a joint-call set we generated across 49 individuals with publicly available PacBio HiFi data. This specific missense mutation (Val473Ala) has been shown in both an ortholog and paralog of Kv3.2 to accelerate current decay, shift the voltage dependence of activation, and prevent the channel from entering a long-lasting open state. Seven additional missense mutations have been identified in other individuals with neurodevelopmental disorders (p = 1.03 × 10−5). KCNC2 is most highly expressed in the brain; in particular, in the thalamus and is enriched in GABAergic neurons. Long-read sequencing was useful in discovering the relevant variant in this family with autism that had remained a mystery for several years and will potentially have great benefits in the clinic once it is widely available.

Abstract +

Synpolydactyly 1 (SPD; MIM# 186000), also called syndactyly type II (SDTY2), is a genetic limb malformation characterized by polydactyly with syndactyly involving the webbing of the third and fourth fingers, and the fourth and fifth toes. It is caused by heterozygous alterations in HOXD13 with incomplete penetrance and phenotypic variability. In our study, a five-generation family with an SPD phenotype was enrolled in our Rare Disease Genomics Protocol. A comprehensive examination of three generations using Illumina short-read whole-genome sequencing (WGS) did not identify any causative variants. Subsequent WGS using Pacific Biosciences (PacBio) long-read HiFi Circular Consensus Sequencing (CCS) revealed a heterozygous 27-bp duplication in the polyalanine tract of HOXD13. Sanger sequencing of all available family members confirmed that the variant segregates with affected individuals. Re-analysis of an unrelated family with a similar SPD phenotype uncovered a 21-bp (7-alanine) duplication in the same region of HOXD13. Although ExpansionHunter identified these events in most individuals in a retrospective analysis, low sequence coverage due to high GC content in the HOXD13 polyalanine tract makes detection of these events challenging. Our findings highlight the value of long-read WGS in elucidating the molecular etiology of congenital limb malformation disorders.

Abstract +

Endosymbiotic Wolbachia bacteria are known to influence the host physiology, microbiota composition, and dissemination of pathogens. We surveyed a population of Tabanus nigrovittatus, commonly referred to as "greenheads," from Crane Beach (Ipswich, MA, USA) for the presence of the alphaproteobacterial symbiont Wolbachia. We studied the COI (mitochondrial cytochrome oxidase) marker gene to evaluate the phylogenetic diversity of the studied specimens. The DNA sequences show strong similarity (between 99.9 and 98%) among the collected specimens but lower similarity to closely related entries in the NCBI database (only between 96.3 and 94.7%), suggesting a more distant relatedness. Low levels of Wolbachia presence necessitated a nested PCR approach, and using 5 markers (ftsZ, fbpA, dnaA, coxA, and gatB), we determined that two recognized "supergroups" of Wolbachia species were represented in the studied specimens, members of clades A and B. Using next-generation sequencing, we also surveyed the insect gut microbiomes of a subset of flies, using Illumina and PacBio 16S rRNA gene sequencing with barcoded primers. The composition of Proteobacteria also varied from fly to fly, with components belonging to Gammaproteobacteria making up the largest percentage of organisms (30 to 70%) among the microbiome samples. Most of the samples showed the presence of Spiroplasma, a member of the phylum Mollicutes, although the frequency of its presence was variable, ranging from 2 to 57%. Another noteworthy bacterial phylum consistently identified was Firmicutes, though the read abundances were typically below 10%. Of interest is an association between Wolbachia presence and higher Alphaproteobacteria representation in the microbiomes, suggesting that the presence of Wolbachia affects the host microbiome. IMPORTANCE Tabanus nigrovittatus greenhead populations contain two supergroups of Wolbachia endosymbionts, members of supergroups A and B. Analysis of the greenhead microbiome using next-generation sequencing revealed that the majority of bacterial species detected belonged to Gammaproteobacteria, with most of the samples also showing the presence of Spiroplasma, a member of the Mollicutes phylum also known to infect insects. An association between Wolbachia presence and higher Alphaproteobacteria representation in the microbiomes suggests that Wolbachia presence affects the host microbiome composition.

Abstract +

Tracking evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) within infected individuals will help elucidate coronavirus disease 2019 (COVID-19) pathogenesis and inform use of antiviral interventions. In this study, we developed an approach for sequencing the region encoding the SARS-CoV-2 virion surface proteins from large numbers of individual virus RNA genomes per sample. We applied this approach to the WA-1 reference clinical isolate of SARS-CoV-2 passaged in vitro and to upper respiratory samples from 7 study participants with COVID-19. SARS-CoV-2 genomes from cell culture were diverse, including 18 haplotypes with non-synonymous mutations clustered in the spike NH2-terminal domain (NTD) and furin cleavage site regions. By contrast, cross-sectional analysis of samples from participants with COVID-19 showed fewer virus variants, without structural clustering of mutations. However, longitudinal analysis in one individual revealed 4 virus haplotypes bearing 3 independent mutations in a spike NTD epitope targeted by autologous antibodies. These mutations arose coincident with a 6.2-fold rise in serum binding to spike and a transient increase in virus burden. We conclude that SARS-CoV-2 exhibits a capacity for rapid genetic adaptation that becomes detectable in vivo with the onset of humoral immunity, with the potential to contribute to delayed virologic clearance in the acute setting.

Abstract +

Recoding viral genomes by introducing numerous synonymous but suboptimal codon pairs—called codon-pair deoptimization (CPD)—provides new types of live-attenuated vaccine candidates. The large number of nucleotide changes resulting from CPD should provide genetic stability to the attenuating phenotype, but this has not been rigorously tested. Human respiratory syncytial virus in which the G and F surface glycoprotein ORFs were CPD (called Min B) was temperature-sensitive and highly restricted in vitro. When subjected to selective pressure by serial passage at increasing temperatures, Min B substantially regained expression of F and replication fitness. Whole-genome deep sequencing showed many point mutations scattered across the genome, including one combination of six linked point mutations. However, their reintroduction into Min B provided minimal rescue. Further analysis revealed viral genomes bearing very large internal deletions (LD genomes) that accumulated after only a few passages. The deletions relocated the CPD F gene to the first or second promoter-proximal gene position. LD genomes amplified de novo in Min B–infected cells were encapsidated, expressed high levels of F, and complemented Min B replication in trans. This study provides insight on a variation of the adaptability of a debilitated negative-strand RNA virus, namely the generation of defective minihelper viruses to overcome its restriction. This is in contrast to the common “defective interfering particles” that interfere with the replication of the virus from which they originated. To our knowledge, defective genomes that promote rather than inhibit replication have not been reported before in RNA viruses.

Abstract +

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent{\textendash}child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average contig N50: 26 Mbp) integrate all forms of genetic variation even across complex loci. We identify 107,590 structural variants (SVs), of which 68\% are not discovered by short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterize 130 of the most active mobile element source elements and find that 63\% of all SVs arise by homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1,526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Abstract +

Antibodies targeting the SARS-CoV-2 spike receptor-binding domain (RBD) are being developed as therapeutics and are a major contributor to neutralizing antibody responses elicited by infection. Here, we describe a deep mutational scanning method to map how all amino-acid mutations in the RBD affect antibody binding and apply this method to 10 human monoclonal antibodies. The escape mutations cluster on several surfaces of the RBD that broadly correspond to structurally defined antibody epitopes. However, even antibodies targeting the same surface often have distinct escape mutations. The complete escape maps predict which mutations are selected during viral growth in the presence of single antibodies. They further enable the design of escape-resistant antibody cocktails—including cocktails of antibodies that compete for binding to the same RBD surface but have different escape mutations. Therefore, complete escape-mutation maps enable rational design of antibody therapeutics and assessment of the antigenic consequences of viral evolution.

Abstract +

We report monozygotic twin girls with syndromic intellectual disability who underwent exome sequencing but with negative pathogenic variants. To search for variants that are unrecognized by exome sequencing, high-fidelity long-read genome sequencing (HiFi LR-GS) was applied. A 12-kb copy-neutral inversion was precisely identified by HiFi LR-GS after trio-based variant filtering. This inversion directly disrupted two genes, CPNE9 and BRPF1, the latter of which attracted our attention because pathogenic BRPF1 variants have been identified in autosomal dominant intellectual developmental disorder with dysmorphic facies and ptosis (IDDDFP), which later turned out to be clinically found in the twins. Trio-based HiFi LR-GS together with haplotype phasing revealed that the 12-kb inversion occurred de novo on the maternally transmitted chromosome. This study clearly indicates that submicroscopic copy-neutral inversions are important but often uncharacterized culprits in monogenic disorders and that long-read sequencing is highly advantageous for detecting such inversions involved in genetic diseases.

Abstract +

Long-read sequencing (LRS) has the potential to comprehensively identify all medically relevant genome variation, including variation commonly missed by short-read sequencing (SRS) approaches. To determine this potential, we performed LRS around 15×–40× genome coverage using the Pacific Biosciences Sequel I System for five trios. The respective probands were diagnosed with intellectual disability (ID) whose etiology remained unresolved after SRS exomes and genomes. Systematic assessment of LRS coverage showed that ~35 Mb of the human reference genome was only accessible by LRS and not SRS. Genome-wide structural variant (SV) calling yielded on average 28,292 SV calls per individual, totaling 12.9 Mb of sequence. Trio-based analyses which allowed to study segregation, showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming that these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, ~100,000 were located in the ~35 Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known Mendelian disease genes. Our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.