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.

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.

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.

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet.

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a 'poor man's crop', finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

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.

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.

TCR-mimic bispecific antibodies to target the HIV-1 reservoir.

HIV-1 infection is incurable due to the persistence of the virus in a latent reservoir of resting memory CD4+ T cells. "Shock-and-kill" approaches that seek to induce HIV-1 gene expression, protein production, and subsequent targeting by the host immune system have been unsuccessful due to a lack of effective latency-reversing agents (LRAs) and kill strategies. In an effort to develop reagents that could be used to promote killing of infected cells, we constructed T cell receptor (TCR)-mimic antibodies to HIV-1 peptide-major histocompatibility complexes (pMHC). Using phage display, we panned for phages expressing antibody-like variable sequences that bound HIV-1 pMHC generated using the common HLA-A*02:01 allele. We targeted three epitopes in Gag and reverse transcriptase identified and quantified via Poisson detection mass spectrometry from cells infected in vitro with a pseudotyped HIV-1 reporter virus (NL4.3 dEnv). Sequences isolated from phages that bound these pMHC were cloned into a single-chain diabody backbone (scDb) sequence, such that one fragment is specific for an HIV-1 pMHC and the other fragment binds to CD3ε, an essential signal transduction subunit of the TCR. Thus, these antibodies utilize the sensitivity of T cell signaling as readouts for antigen processing and as agents to promote killing of infected cells. Notably, these scDbs are exquisitely sensitive and specific for the peptide portion of the pMHC. Most importantly, one scDb caused killing of infected cells presenting a naturally processed target pMHC. This work lays the foundation for a novel therapeutic killing strategy toward elimination of the HIV-1 reservoir.

The XPB Subunit of the TFIIH Complex Plays a Critical Role in HIV-1 Transcription and XPB Inhibition by Spironolactone Prevents HIV-1 Reactivation from Latency.

HIV transcription requires assembly of cellular transcription factors at the HIV-1promoter. The TFIIH general transcription factor facilitates transcription initiation by opening the DNA strands around the transcription start site and phosphorylating the C-terminal domain for RNA polymerase II (RNAPII) for activation. Spironolactone (SP), an FDA approved aldosterone antagonist, triggers the proteasomal degradation of the XPB subunit of TFIIH, and concurrently suppresses acute HIV infection in vitro Here we investigated SP as a possible block-and-lock agent for a functional cure aimed at the transcriptional silencing of the viral reservoir. The long-term activity of SP was investigated in primary and cell line models of HIV-1 latency and reactivation. We show that SP rapidly inhibits HIV-1 transcription by reducing RNAPII recruitment to the HIV-1 genome. shRNA knockdown of XPB confirmed XPB degradation as the mechanism of action. Unfortunately, long-term pre-treatment with SP does not result in epigenetic suppression of HIV upon SP treatment interruption, since virus rapidly rebounds when XPB reemerges; however, SP alone without ART maintains the transcriptional suppression. Importantly, SP inhibits HIV reactivation from latency in both cell line models and resting CD4+T cells isolated from aviremic infected individuals upon cell stimulation with latency reversing agents. Furthermore, long-term treatment with concentrations of SP that potently degrade XPB does not lead to global dysregulation of cellular mRNA expression. Overall, these results suggest that XPB plays a key role in HIV transcriptional regulation and XPB degradation by SP strengthens the potential of HIV transcriptional inhibitors in block-and-lock HIV cure approaches.IMPORTANCE Antiretroviral therapy (ART) effectively reduces an individual's HIV loads to below the detection limit, nevertheless rapid viral rebound immediately ensues upon treatment interruption. Furthermore, virally suppressed individuals experience chronic immune activation from ongoing low-level virus expression. Thus, the importance of identifying novel therapeutics to explore in block-and-lock HIV functional cure approaches, aimed at the transcriptional and epigenetic silencing of the viral reservoir to block reactivation from latency. We investigated the potential of repurposing the FDA-approved spironolactone (SP), as one such drug. SP treatment rapidly degrades a host transcription factor subunit, XPB, inhibiting HIV transcription and blocking reactivation from latency. Long-term SP treatment does not affect cellular viability, cell cycle progression or global cellular transcription. SP alone blocks HIV transcription in the absence of ART but does not delay rebound upon drug removal as XPB rapidly reemerges. This study highlights XPB as a novel drug target in block-and-lock therapeutic approaches.

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.

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.