HIVepsilon-seq-scalable characterization of intact persistent proviral HIV reservoirs in women.

IMPORTANCE: The lack of a reliable method to accurately detect when replication-competent HIV has been cleared is a major challenge in developing a cure. This study introduces a new approach called the HIVepsilon-seq (HIVε-seq) assay, which uses long-read sequencing technology and bioinformatics to scrutinize the HIV genome at the nucleotide level, distinguishing between defective and intact HIV. This study included 30 participants on antiretroviral therapy, including 17 women, and was able to discriminate between defective and genetically intact viruses at the single DNA strand level. The HIVε-seq assay is an improvement over previous methods, as it requires minimal sample, less specialized lab equipment, and offers a shorter turnaround time. The HIVε-seq assay offers a promising new tool for researchers to measure the intact HIV reservoir, advancing efforts towards finding a cure for this devastating disease.

Molecular barcoding of native RNAs using nanopore sequencing and deep learning.

Nanopore sequencing enables direct measurement of RNA molecules without conversion to cDNA, thus opening the gates to a new era for RNA biology. However, the lack of molecular barcoding of direct RNA nanopore sequencing data sets severely affects the applicability of this technology to biological samples, where RNA availability is often limited. Here, we provide the first experimental protocol and associated algorithm to barcode and demultiplex direct RNA nanopore sequencing data sets. Specifically, we present a novel and robust approach to accurately classify raw nanopore signal data by transforming current intensities into images or arrays of pixels, followed by classification using a deep learning algorithm. We demonstrate the power of this strategy by developing the first experimental protocol for barcoding and demultiplexing direct RNA sequencing libraries. Our method, DeePlexiCon, can classify 93% of reads with 95.1% accuracy or 60% of reads with 99.9% accuracy. The availability of an efficient and simple multiplexing strategy for native RNA sequencing will improve the cost-effectiveness of this technology, as well as facilitate the analysis of lower-input biological samples. Overall, our work exemplifies the power, simplicity, and robustness of signal-to-image conversion for nanopore data analysis using deep learning.

Unexpected host dependency of Antarctic Nanohaloarchaeota.

In hypersaline environments, Nanohaloarchaeota (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaeota [DPANN] superphylum) are thought to be free-living microorganisms. We report cultivation of 2 strains of Antarctic Nanohaloarchaeota and show that they require the haloarchaeon Halorubrum lacusprofundi for growth. By performing growth using enrichments and fluorescence-activated cell sorting, we demonstrated successful cultivation of Candidatus Nanohaloarchaeum antarcticus, purification of Ca. Nha. antarcticus away from other species, and growth and verification of Ca. Nha. antarcticus with Hrr. lacusprofundi; these findings are analogous to those required for fulfilling Koch's postulates. We use fluorescent in situ hybridization and transmission electron microscopy to assess cell structures and interactions; metagenomics to characterize enrichment taxa, generate metagenome assembled genomes, and interrogate Antarctic communities; and proteomics to assess metabolic pathways and speculate about the roles of certain proteins. Metagenome analysis indicates the presence of a single species, which is endemic to Antarctic hypersaline systems that support the growth of haloarchaea. The presence of unusually large proteins predicted to function in attachment and invasion of hosts plus the absence of key biosynthetic pathways (e.g., lipids) in metagenome assembled genomes of globally distributed Nanohaloarchaeota indicate that all members of the lineage have evolved as symbionts. Our work expands the range of archaeal symbiotic lifestyles and provides a genetically tractable model system for advancing understanding of the factors controlling microbial symbiotic relationships.

Adaptation of Oxford Nanopore technology for hepatitis C whole genome sequencing and identification of within-host viral variants.

BACKGROUND: Hepatitis C (HCV) and many other RNA viruses exist as rapidly mutating quasi-species populations in a single infected host. High throughput characterization of full genome, within-host variants is still not possible despite advances in next generation sequencing. This limitation constrains viral genomic studies that depend on accurate identification of hemi-genome or whole genome, within-host variants, especially those occurring at low frequencies. With the advent of third generation long read sequencing technologies, including Oxford Nanopore Technology (ONT) and PacBio platforms, this problem is potentially surmountable. ONT is particularly attractive in this regard due to the portable nature of the MinION sequencer, which makes real-time sequencing in remote and resource-limited locations possible. However, this technology (termed here 'nanopore sequencing') has a comparatively high technical error rate. The present study aimed to assess the utility, accuracy and cost-effectiveness of nanopore sequencing for HCV genomes. We also introduce a new bioinformatics tool (Nano-Q) to differentiate within-host variants from nanopore sequencing. RESULTS: The Nanopore platform, when the coverage exceeded 300 reads, generated comparable consensus sequences to Illumina sequencing. Using HCV Envelope plasmids (~ 1800 nt) mixed in known proportions, the capacity of nanopore sequencing to reliably identify variants with an abundance as low as 0.1% was demonstrated, provided the autologous reference sequence was available to identify the matching reads. Successful pooling and nanopore sequencing of 52 samples from patients with HCV infection demonstrated its cost effectiveness (AUD$ 43 per sample with nanopore sequencing versus $100 with paired-end short read technology). The Nano-Q tool successfully separated between-host sequences, including those from the same subtype, by bulk sorting and phylogenetic clustering without an autologous reference sequence (using only a subtype-specific generic reference). The pipeline also identified within-host viral variants and their abundance when the parameters were appropriately adjusted. CONCLUSION: Cost effective HCV whole genome sequencing and within-host variant identification without haplotype reconstruction are potential advantages of nanopore sequencing.

How to Define the Latent Reservoir: Tools of the Trade.

HIV is a devastating worldwide epidemic that has had substantial social and economic impacts throughout the globe. Due to the presence of a small pool of latently infected cells that persists during antiretroviral therapy (ART), HIV is not curable. Because of the high cost of ART and the lack of reliable accessibility across the globe, life-long ART is unfortunately not a feasible solution for the epidemic. Therefore, new strategies need to be developed and implemented to address HIV-1 infection. Several approaches toward this end are currently under investigation (Ebina et al. in Sci Rep 3:2510, 2013; Archin et al. in Nature 487:482-5, 2012; Elliott et al. in PLoS Pathog 10:e1004473, 2014; Rasmussen et al. in Lancet HIV 1:e13-e21, 2014; Tebas et al. in N Engl J Med 370:901-10, 2014; Archin et al. in Nat Rev Microbiol 12:750-64, 2014; Barton et al. in PLoS One 9:e102684, 2014; Sogaard et al. in PLoS Pathog 11:e1005142, 2015). Initial studies have proven promising, but have highlighted the need for sensitive and accurate assays to detect changes in very low concentrations of virus to allow confident interpretation of the success of curative approaches. This review will focus on assays that are currently available and the advantages and limitations of each.

The genetic bottleneck in vertical transmission of subtype C HIV-1 is not driven by selection of especially neutralization-resistant virus from the maternal viral population.

Subtype C human immunodeficiency virus type 1 (HIV-1C) continues to cause the majority of new cases of mother-to-child transmission (MTCT), and yet there are limited data on HIV-1C transmission. We amplified env from plasma RNA for 19 HIV-1C MTCT pairs, 10 transmitting in utero (IU) and 9 transmitting intrapartum (IP). There was a strong genetic bottleneck between all mother-infant pairs, with a majority of transmission events involving the transmission of a single virus. env genes of viruses transmitted to infants IP, but not IU, encoded Env proteins that were shorter and had fewer putative N-linked glycosylation sites in the V1-V5 region than matched maternal sequences. Viruses pseudotyped with env clones representative of each maternal and infant population were tested for neutralization sensitivity. The 50% inhibitory concentration of autologous serum was similar against both transmitted (infant) and nontransmitted (maternal) viruses in a paired analysis. Mother and infant Env proteins were also similar in sensitivity to soluble CD4, to a panel of monoclonal antibodies, and to heterologous HIV-1C sera. In addition, there was no difference in the breadth or potency of neutralizing antibodies between sera from 50 nontransmitting and 23 IU and 23 IP transmitting HIV-1C-infected women against four Env proteins from heterologous viruses. Thus, while a strong genetic bottleneck was detected during MCTC, with viruses of shorter and fewer glycosylation sites in env present in IP transmission, our data do not support this bottleneck being driven by selective resistance to antibodies.

Barriers and potential solutions to international collaboration in neuro-oncology clinical trials: Challenges from the Australian perspective.

AIM: The neuro-oncology community in Australia is well positioned to collaborate internationally, with a motivated trials group, strong regulatory bodies and an attractive fiscal environment. We sought to identify gaps in the Australian neuro-oncology clinical trials landscape and describe strategies to increase international trial access in Australia. METHODS: We searched clinical trial registries to identify active adult primary brain cancer trials. We compared the participation rate and phase of these trials between tumour types and countries. A survey was distributed to the Cooperative Trials Group for Neuro-Oncology membership to identify barriers and solutions to effective international collaboration. RESULTS: Globally, 307 trials for adult primary brain cancers were identified. These included 50% pharmaceutical agents, 18% cellular therapies and 9% radiation therapy. Twelve adult primary brain cancer trials were actively recruiting in Australia at the time the survey was sent out. There were more early phase brain cancer trials (34%) compared with colorectal and breast cancer (21% and 24%, respectively). In Australia, 92% of brain cancer trials were involving pharmaceutical agents. The most commonly cited barrier was lack of funding for international trials (86%) and insufficient research time (75%). High ranking solutions included increasing the availability of funding for international trials and creating opportunities to develop personal relationships with collaborators. Accreditation of clinical research key performance indicators into practice (88%) and hospital accreditation (73%) also ranked highly. CONCLUSIONS: Participation in international research in Australia could be improved by embedding clinical research targets into institutional funding, provision of funding for early phase studies and streamlining mutual ethics schemes.

sangeranalyseR: Simple and Interactive Processing of Sanger Sequencing Data in R.

sangeranalyseR is feature-rich, free, and open-source R package for processing Sanger sequencing data. It allows users to go from loading reads to saving aligned contigs in a few lines of R code by using sensible defaults for most actions. It also provides complete flexibility for determining how individual reads and contigs are processed, both at the command-line in R and via interactive Shiny applications. sangeranalyseR provides a wide range of options for all steps in Sanger processing pipelines including trimming reads, detecting secondary peaks, viewing chromatograms, detecting indels and stop codons, aligning contigs, estimating phylogenetic trees, and more. Input data can be in either ABIF or FASTA format. sangeranalyseR comes with extensive online documentation and outputs aligned and unaligned reads and contigs in FASTA format, along with detailed interactive HTML reports. sangeranalyseR supports the use of colorblind-friendly palettes for viewing alignments and chromatograms. It is released under an MIT licence and available for all platforms on Bioconductor (https://bioconductor.org/packages/sangeranalyseR, last accessed February 22, 2021) and on Github (https://github.com/roblanf/sangeranalyseR, last accessed February 22, 2021).

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Single molecule, near full-length genome sequencing of dengue virus.

Current methods for dengue virus (DENV) genome amplification, amplify parts of the genome in at least 5 overlapping segments and then combine the output to characterize a full genome. This process is laborious, costly and requires at least 10 primers per serotype, thus increasing the likelihood of PCR bias. We introduce an assay to amplify near full-length dengue virus genomes as intact molecules, sequence these amplicons with third generation "nanopore" technology without fragmenting and use the sequence data to differentiate within-host viral variants with a bioinformatics tool (Nano-Q). The new assay successfully generated near full-length amplicons from DENV serotypes 1, 2 and 3 samples which were sequenced with nanopore technology. Consensus DENV sequences generated by nanopore sequencing had over 99.5% pairwise sequence similarity to Illumina generated counterparts provided the coverage was > 100 with both platforms. Maximum likelihood phylogenetic trees generated from nanopore consensus sequences were able to reproduce the exact trees made from Illumina sequencing with a conservative 99% bootstrapping threshold (after 1000 replicates and 10% burn-in). Pairwise genetic distances of within host variants identified from the Nano-Q tool were less than that of between host variants, thus enabling the phylogenetic segregation of variants from the same host.

Canfam_GSD: De novo chromosome-length genome assembly of the German Shepherd Dog (Canis lupus familiaris) using a combination of long reads, optical mapping, and Hi-C.

BACKGROUND: The German Shepherd Dog (GSD) is one of the most common breeds on earth and has been bred for its utility and intelligence. It is often first choice for police and military work, as well as protection, disability assistance, and search-and-rescue. Yet, GSDs are well known to be susceptible to a range of genetic diseases that can interfere with their training. Such diseases are of particular concern when they occur later in life, and fully trained animals are not able to continue their duties. FINDINGS: Here, we provide the draft genome sequence of a healthy German Shepherd female as a reference for future disease and evolutionary studies. We generated this improved canid reference genome (CanFam_GSD) utilizing a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. The GSD assembly is ∼80 times as contiguous as the current canid reference genome (20.9 vs 0.267 Mb contig N50), containing far fewer gaps (306 vs 23,876) and fewer scaffolds (429 vs 3,310) than the current canid reference genome CanFamv3.1. Two chromosomes (4 and 35) are assembled into single scaffolds with no gaps. BUSCO analyses of the genome assembly results show that 93.0% of the conserved single-copy genes are complete in the GSD assembly compared with 92.2% for CanFam v3.1. Homology-based gene annotation increases this value to ∼99%. Detailed examination of the evolutionarily important pancreatic amylase region reveals that there are most likely 7 copies of the gene, indicative of a duplication of 4 ancestral copies and the disruption of 1 copy. CONCLUSIONS: GSD genome assembly and annotation were produced with major improvement in completeness, continuity, and quality over the existing canid reference. This resource will enable further research related to canine diseases, the evolutionary relationships of canids, and other aspects of canid biology.

Chiral DNA sequences as commutable controls for clinical genomics.

Chirality is a property describing any object that is inequivalent to its mirror image. Due to its 5'-3' directionality, a DNA sequence is distinct from a mirrored sequence arranged in reverse nucleotide-order, and is therefore chiral. A given sequence and its opposing chiral partner sequence share many properties, such as nucleotide composition and sequence entropy. Here we demonstrate that chiral DNA sequence pairs also perform equivalently during molecular and bioinformatic techniques that underpin genetic analysis, including PCR amplification, hybridization, whole-genome, target-enriched and nanopore sequencing, sequence alignment and variant detection. Given these shared properties, synthetic DNA sequences mirroring clinically relevant or analytically challenging regions of the human genome are ideal controls for clinical genomics. The addition of synthetic chiral sequences (sequins) to patient tumor samples can prevent false-positive and false-negative mutation detection to improve diagnosis. Accordingly, we propose that sequins can fulfill the need for commutable internal controls in precision medicine.

Targeted, High-Resolution RNA Sequencing of Non-coding Genomic Regions Associated With Neuropsychiatric Functions.

The human brain is one of the last frontiers of biomedical research. Genome-wide association studies (GWAS) have succeeded in identifying thousands of haplotype blocks associated with a range of neuropsychiatric traits, including disorders such as schizophrenia, Alzheimer's and Parkinson's disease. However, the majority of single nucleotide polymorphisms (SNPs) that mark these haplotype blocks fall within non-coding regions of the genome, hindering their functional validation. While some of these GWAS loci may contain cis-acting regulatory DNA elements such as enhancers, we hypothesized that many are also transcribed into non-coding RNAs that are missing from publicly available transcriptome annotations. Here, we use targeted RNA capture ('RNA CaptureSeq') in combination with nanopore long-read cDNA sequencing to transcriptionally profile 1,023 haplotype blocks across the genome containing non-coding GWAS SNPs associated with neuropsychiatric traits, using post-mortem human brain tissue from three neurologically healthy donors. We find that the majority (62%) of targeted haplotype blocks, including 13% of intergenic blocks, are transcribed into novel, multi-exonic RNAs, most of which are not yet recorded in GENCODE annotations. We validated our findings with short-read RNA-seq, providing orthogonal confirmation of novel splice junctions and enabling a quantitative assessment of the long-read assemblies. Many novel transcripts are supported by independent evidence of transcription including cap analysis of gene expression (CAGE) data and epigenetic marks, and some show signs of potential functional roles. We present these transcriptomes as a preliminary atlas of non-coding transcription in human brain that can be used to connect neurological phenotypes with gene expression.

High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes.

High-throughput single-cell RNA sequencing is a powerful technique but only generates short reads from one end of a cDNA template, limiting the reconstruction of highly diverse sequences such as antigen receptors. To overcome this limitation, we combined targeted capture and long-read sequencing of T-cell-receptor (TCR) and B-cell-receptor (BCR) mRNA transcripts with short-read transcriptome profiling of barcoded single-cell libraries generated by droplet-based partitioning. We show that Repertoire and Gene Expression by Sequencing (RAGE-Seq) can generate accurate full-length antigen receptor sequences at nucleotide resolution, infer B-cell clonal evolution and identify alternatively spliced BCR transcripts. We apply RAGE-Seq to 7138 cells sampled from the primary tumor and draining lymph node of a breast cancer patient to track transcriptome profiles of expanded lymphocyte clones across tissues. Our results demonstrate that RAGE-Seq is a powerful method for tracking the clonal evolution from large numbers of lymphocytes applicable to the study of immunity, autoimmunity and cancer.

Intronic ATTTC repeat expansions in STARD7 in familial adult myoclonic epilepsy linked to chromosome 2.

Familial Adult Myoclonic Epilepsy (FAME) is characterised by cortical myoclonic tremor usually from the second decade of life and overt myoclonic or generalised tonic-clonic seizures. Four independent loci have been implicated in FAME on chromosomes (chr) 2, 3, 5 and 8. Using whole genome sequencing and repeat primed PCR, we provide evidence that chr2-linked FAME (FAME2) is caused by an expansion of an ATTTC pentamer within the first intron of STARD7. The ATTTC expansions segregate in 158/158 individuals typically affected by FAME from 22 pedigrees including 16 previously reported families recruited worldwide. RNA sequencing from patient derived fibroblasts shows no accumulation of the AUUUU or AUUUC repeat sequences and STARD7 gene expression is not affected. These data, in combination with other genes bearing similar mutations that have been implicated in FAME, suggest ATTTC expansions may cause this disorder, irrespective of the genomic locus involved.

Romidepsin-induced HIV-1 viremia during effective antiretroviral therapy contains identical viral sequences with few deleterious mutations.

OBJECTIVE: To investigate the origin of the HIV-1 viremia induced by the latency-reversing agent romidepsin. DESIGN: Six individuals on suppressive antiretroviral therapy received romidepsin administered intravenously once weekly for 3 consecutive weeks. CD4 T cells were obtained at baseline, following the second and third romidepsin infusion, and 10 weeks after the final romidepsin treatment. Plasma samples were collected 24 and 72 h after romidepsin infusions. METHODS: Single-genome sequencing of the env and p24-RT region was used to genetically characterize the virus from proviral DNA, the transcribed cell-associated RNA and the plasma RNA pool. RESULTS: In three of six participants with available plasma samples we identified plasma HIV-1 RNA sequences that were identical to DNA and/or cell-associated RNA sequences from peripheral blood CD4 T cells. In two participants, plasma RNA sequences contained expansions of identical sequences, corresponding to 62 and 100% of the total sequences, respectively. Plasma HIV-1 RNA had very low amounts of defective viruses compared to cell-associated RNA (odds ratio 20.85, P < 0.001) and to DNA (odds ratio 7.07, P = 0.011) during romidepsin therapy. CONCLUSIONS: Romidepsin induced transcription from proviruses in peripheral blood cells, which contributed to viremia in patients on suppressive therapy. The intermingling of these cell-associated HIV-1 RNA with DNA sequences indicates transcription from a diverse range of proviruses, but the expansions of identical viral plasma sequences with few defects indicate that the romidepsin-induced viremia arises from intact proviruses with highly similar or identical genetic backgrounds.

Identification of Genetically Intact HIV-1 Proviruses in Specific CD4+ T Cells from Effectively Treated Participants.

Latent replication-competent HIV-1 persists in individuals on long-term antiretroviral therapy (ART). We developed the Full-Length Individual Proviral Sequencing (FLIPS) assay to determine the distribution of latent replication-competent HIV-1 within memory CD4+ T cell subsets in six individuals on long-term ART. FLIPS is an efficient, high-throughput assay that amplifies and sequences near full-length (∼9 kb) HIV-1 proviral genomes and determines potential replication competency through genetic characterization. FLIPS provides a genome-scale perspective that addresses the limitations of other methods that also genetically characterize the latent reservoir. Using FLIPS, we identified 5% of proviruses as intact and potentially replication competent. Intact proviruses were unequally distributed between T cell subsets, with effector memory cells containing the largest proportion of genetically intact HIV-1 proviruses. We identified multiple identical intact proviruses, suggesting a role for cellular proliferation in the maintenance of the latent HIV-1 reservoir.

HIV-1 Reservoirs During Suppressive Therapy.

The introduction of antiretroviral therapy (ART) 20 years ago has dramatically reduced morbidity and mortality associated with HIV-1. Initially there was hope that ART would be curative, but it quickly became clear that even though ART was able to restore CD4(+) T cell counts and suppress viral loads below levels of detection, discontinuation of treatment resulted in a rapid rebound of infection. This is due to persistence of a small reservoir of latently infected cells with a long half-life, which necessitates life-long ART. Over the past few years, significant progress has been made in defining and characterizing the latent reservoir of HIV-1, and here we review how understanding the latent reservoir during suppressive therapy will lead to significant advances in curative approaches for HIV-1.

Broad activation of latent HIV-1 in vivo.

The 'shock and kill' approach to cure human immunodeficiency virus (HIV) includes transcriptional induction of latent HIV-1 proviruses using latency-reversing agents (LRAs) with targeted immunotherapy to purge infected cells. The administration of LRAs (panobinostat or vorinostat) to HIV-1-infected individuals on antiretroviral therapy induces a significant increase in cell-associated unspliced (CA-US) HIV-1 RNA from CD4(+) T cells. However, it is important to discern whether the increases in CA-US HIV-1 RNA are due to limited or broad activation of HIV-1 proviruses. Here we use single-genome sequencing to find that the RNA transcripts observed following LRA administration are genetically diverse, indicating activation of transcription from an extensive range of proviruses. Defective sequences are more frequently found in CA HIV-1 RNA than in HIV-1 DNA, which has implications for developing an accurate measure of HIV-1 reservoir size. Our findings provide insights into the effects of panobinostat and vorinostat as LRAs for latent HIV-1.