Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References.

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

A quantitative approach for measuring the reservoir of latent HIV-1 proviruses.

A stable latent reservoir for HIV-1 in resting CD4+ T cells is the principal barrier to a cure1-3. Curative strategies that target the reservoir are being tested4,5 and require accurate, scalable reservoir assays. The reservoir was defined with quantitative viral outgrowth assays for cells that release infectious virus after one round of T cell activation1. However, these quantitative outgrowth assays and newer assays for cells that produce viral RNA after activation6 may underestimate the reservoir size because one round of activation does not induce all proviruses7. Many studies rely on simple assays based on polymerase chain reaction to detect proviral DNA regardless of transcriptional status, but the clinical relevance of these assays is unclear, as the vast majority of proviruses are defective7-9. Here we describe a more accurate method of measuring the HIV-1 reservoir that separately quantifies intact and defective proviruses. We show that the dynamics of cells that carry intact and defective proviruses are different in vitro and in vivo. These findings have implications for targeting the intact proviruses that are a barrier to curing HIV infection.

Expression Microdissection for the Analysis of miRNA in a Single-Cell Type.

Cell-specific microRNA (miRNA) expression estimates are important in characterizing the localization of miRNA signaling within tissues. Much of these data are obtained from cultured cells, a process known to significantly alter miRNA expression levels. Thus, our knowledge of in vivo cell miRNA expression estimates is poor. We previously demonstrated expression microdissection-miRNA-sequencing (xMD-miRNA-seq) to acquire in vivo estimates, directly from formalin-fixed tissues, albeit with a limited yield. In this study, we optimized each step of the xMD process, including tissue retrieval, tissue transfer, film preparation, and RNA isolation, to increase RNA yields and ultimately show strong enrichment for in vivo miRNA expression by qPCR array. These method improvements, such as the development of a noncrosslinked ethylene vinyl acetate membrane, resulted in a 23- to 45-fold increase in miRNA yield, depending on the cell type. By qPCR, miR-200a increased by 14-fold in xMD-derived small intestine epithelial cells, with a concurrent 336-fold reduction in miR-143 relative to the matched nondissected duodenal tissue. xMD is now an optimized method to obtain robust in vivo miRNA expression estimates from cells. xMD will allow formalin-fixed tissues from surgical pathology archives to make theragnostic biomarker discoveries.

Direct observation of the evolution of cell-type-specific microRNA expression signatures supports the hematopoietic origin model of endothelial cells.

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate-specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative. We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, which was a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization, we localized Mir-126 in the tunicate Ciona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates. This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.

MEM-based pangenome indexing for k-mer queries.

Pangenomes are growing in number and size, thanks to the prevalence of high-quality long-read assemblies. However, current methods for studying sequence composition and conservation within pangenomes have limitations. Methods based on graph pangenomes require a computationally expensive multiple-alignment step, which can leave out some variation. Indexes based on k-mers and de Bruijn graphs are limited to answering questions at a specific substring length k. We present Maximal Exact Match Ordered (MEMO), a pangenome indexing method based on maximal exact matches (MEMs) between sequences. A single MEMO index can handle arbitrary-length queries over pangenomic windows. MEMO enables both queries that test k-mer presence/absence (membership queries) and that count the number of genomes containing k-mers in a window (conservation queries). MEMO's index for a pangenome of 89 human autosomal haplotypes fits in 2.04 GB, 8.8× smaller than a comparable KMC3 index and 11.4× smaller than a PanKmer index. MEMO indexes can be made smaller by sacrificing some counting resolution, with our decile-resolution HPRC index reaching 0.67 GB. MEMO can conduct a conservation query for 31-mers over the human leukocyte antigen locus in 13.89 seconds, 2.5x faster than other approaches. MEMO's small index size, lack of k-mer length dependence, and efficient queries make it a flexible tool for studying and visualizing substring conservation in pangenomes.

Uncalled4 improves nanopore DNA and RNA modification detection via fast and accurate signal alignment.

Nanopore signal analysis enables detection of nucleotide modifications from native DNA and RNA sequencing, providing both accurate genetic/transcriptomic and epigenetic information without additional library preparation. Presently, only a limited set of modifications can be directly basecalled (e.g. 5-methylcytosine), while most others require exploratory methods that often begin with alignment of nanopore signal to a nucleotide reference. We present Uncalled4, a toolkit for nanopore signal alignment, analysis, and visualization. Uncalled4 features an efficient banded signal alignment algorithm, BAM signal alignment file format, statistics for comparing signal alignment methods, and a reproducible de novo training method for k-mer-based pore models, revealing potential errors in ONT's state-of-the-art DNA model. We apply Uncalled4 to RNA 6-methyladenine (m6A) detection in seven human cell lines, identifying 26% more modifications than Nanopolish using m6Anet, including in several genes where m6A has known implications in cancer. Uncalled4 is available open-source at github.com/skovaka/uncalled4.

Gapless assembly of complete human and plant chromosomes using only nanopore sequencing.

The combination of ultra-long Oxford Nanopore (ONT) sequencing reads with long, accurate PacBio HiFi reads has enabled the completion of a human genome and spurred similar efforts to complete the genomes of many other species. However, this approach for complete, "telomere-to-telomere" genome assembly relies on multiple sequencing platforms, limiting its accessibility. ONT "Duplex" sequencing reads, where both strands of the DNA are read to improve quality, promise high per-base accuracy. To evaluate this new data type, we generated ONT Duplex data for three widely-studied genomes: human HG002, Solanum lycopersicum Heinz 1706 (tomato), and Zea mays B73 (maize). For the diploid, heterozygous HG002 genome, we also used "Pore-C" chromatin contact mapping to completely phase the haplotypes. We found the accuracy of Duplex data to be similar to HiFi sequencing, but with read lengths tens of kilobases longer, and the Pore-C data to be compatible with existing diploid assembly algorithms. This combination of read length and accuracy enables the construction of a high-quality initial assembly, which can then be further resolved using the ultra-long reads, and finally phased into chromosome-scale haplotypes with Pore-C. The resulting assemblies have a base accuracy exceeding 99.999% (Q50) and near-perfect continuity, with most chromosomes assembled as single contigs. We conclude that ONT sequencing is a viable alternative to HiFi sequencing for de novo genome assembly, and has the potential to provide a single-instrument solution for the reconstruction of complete genomes.

Convergent evolution of plant prickles by repeated gene co-option over deep time.

An enduring question in evolutionary biology concerns the degree to which episodes of convergent trait evolution depend on the same genetic programs, particularly over long timescales. In this work, we genetically dissected repeated origins and losses of prickles-sharp epidermal projections-that convergently evolved in numerous plant lineages. Mutations in a cytokinin hormone biosynthetic gene caused at least 16 independent losses of prickles in eggplants and wild relatives in the genus Solanum. Homologs underlie prickle formation across angiosperms that collectively diverged more than 150 million years ago, including rice and roses. By developing new Solanum genetic systems, we leveraged this discovery to eliminate prickles in a wild species and an indigenously foraged berry. Our findings implicate a shared hormone activation genetic program underlying evolutionarily widespread and recurrent instances of plant morphological innovation.

Filgotinib suppresses HIV-1-driven gene transcription by inhibiting HIV-1 splicing and T cell activation.

Despite effective antiretroviral therapy, HIV-1-infected cells continue to produce viral antigens and induce chronic immune exhaustion. We propose to identify HIV-1-suppressing agents that can inhibit HIV-1 reactivation and reduce HIV-1-induced immune activation. Using a newly developed dual-reporter system and a high-throughput drug screen, we identified FDA-approved drugs that can suppress HIV-1 reactivation in both cell line models and CD4+ T cells from virally suppressed HIV-1-infected individuals. We identified 11 cellular pathways required for HIV-1 reactivation as druggable targets. Using differential expression analysis, gene set enrichment analysis, and exon-intron landscape analysis, we examined the impact of drug treatment on the cellular environment at a genome-wide level. We identified what we believe to be a new function of a JAK inhibitor, filgotinib, that suppresses HIV-1 splicing. First, filgotinib preferentially suppresses spliced HIV-1 RNA transcription. Second, filgotinib suppresses HIV-1-driven aberrant cancer-related gene expression at the integration site. Third, we found that filgotinib suppresses HIV-1 transcription by inhibiting T cell activation and by modulating RNA splicing. Finally, we found that filgotinib treatment reduces the proliferation of HIV-1-infected cells. Overall, the combination of a drug screen and transcriptome analysis provides systematic understanding of cellular targets required for HIV-1 reactivation and drug candidates that may reduce HIV-1-related immune activation.

Longitudinal study reveals HIV-1-infected CD4+ T cell dynamics during long-term antiretroviral therapy.

Proliferation of CD4+ T cells harboring HIV-1 proviruses is a major contributor to viral persistence in people on antiretroviral therapy (ART). To determine whether differential rates of clonal proliferation or HIV-1-specific cytotoxic T lymphocyte (CTL) pressure shape the provirus landscape, we performed an intact proviral DNA assay (IPDA) and obtained 661 near-full-length provirus sequences from 8 individuals with suppressed viral loads on ART at time points 7 years apart. We observed slow decay of intact proviruses but no changes in the proportions of various types of defective proviruses. The proportion of intact proviruses in expanded clones was similar to that of defective proviruses in clones. Intact proviruses observed in clones did not have more escaped CTL epitopes than intact proviruses observed as singlets. Concordantly, total proviruses at later time points or observed in clones were not enriched in escaped or unrecognized epitopes. Three individuals with natural control of HIV-1 infection (controllers) on ART, included because controllers have strong HIV-1-specific CTL responses, had a smaller proportion of intact proviruses but a distribution of defective provirus types and escaped or unrecognized epitopes similar to that of the other individuals. This work suggests that CTL selection does not significantly check clonal proliferation of infected cells or greatly alter the provirus landscape in people on ART.

The Landscape of Persistent Viral Genomes in ART-Treated SIV, SHIV, and HIV-2 Infections.

Evaluation of HIV cure strategies is complicated by defective proviruses that persist in ART-treated patients but are irrelevant to cure. Non-human primates (NHP) are essential for testing cure strategies. However, the persisting proviral landscape in ART-treated NHPs is uncharacterized. Here, we describe viral genomes persisting in ART-treated, simian immunodeficiency virus (SIV)-infected NHPs, simian-human immunodeficiency virus (SHIV)-infected NHPs, and humans infected with HIV-2, an SIV-related virus. The landscapes of persisting SIV, SHIV, and HIV-2 genomes are also dominated by defective sequences. However, there was a significantly higher fraction of intact SIV proviral genomes compared to ART-treated HIV-1 or HIV-2 infected humans. Compared to humans with HIV-1, SIV-infected NHPs had more hypermutated genomes, a relative paucity of clonal SIV sequences, and a lower frequency of deleted genomes. Finally, we report an assay for measuring intact SIV genomes which may have value in cure research.