Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory.

Invasive species offer insights into rapid adaptation to novel environments. The iconic cane toad (Rhinella marina) is an excellent model for studying rapid adaptation during invasion. Previous research using the mitochondrial NADH dehydrogenase 3 (ND3) gene in Hawai'ian and Australian invasive populations found a single haplotype, indicating an extreme genetic bottleneck following introduction. Nuclear genetic diversity also exhibited reductions across the genome in these two populations. Here, we investigated the mitochondrial genomics of cane toads across this invasion trajectory. We created the first reference mitochondrial genome for this species using long-read sequence data. We combined whole-genome resequencing data of 15 toads with published transcriptomic data of 125 individuals to construct nearly complete mitochondrial genomes from the native (French Guiana) and introduced (Hawai'i and Australia) ranges for population genomic analyses. In agreement with previous investigations of these populations, we identified genetic bottlenecks in both Hawai'ian and Australian introduced populations, alongside evidence of population expansion in the invasive ranges. Although mitochondrial genetic diversity in introduced populations was reduced, our results revealed that it had been underestimated: we identified 45 mitochondrial haplotypes in Hawai'ian and Australian samples, none of which were found in the native range. Additionally, we identified two distinct groups of haplotypes from the native range, separated by a minimum of 110 base pairs (0.6%). These findings enhance our understanding of how invasion has shaped the genetic landscape of this species.

CAR+ and CAR- T cells share a differentiation trajectory into an NK-like subset after CD19 CAR T cell infusion in patients with B cell malignancies.

Chimeric antigen receptor (CAR) T cell therapy is effective in treating B cell malignancies, but factors influencing the persistence of functional CAR+ T cells, such as product composition, patients' lymphodepletion, and immune reconstitution, are not well understood. To shed light on this issue, here we conduct a single-cell multi-omics analysis of transcriptional, clonal, and phenotypic profiles from pre- to 1-month post-infusion of CAR+ and CAR- T cells from patients from a CARTELL study (ACTRN12617001579381) who received a donor-derived 4-1BB CAR product targeting CD19. Following infusion, CAR+ T cells and CAR- T cells shows similar differentiation profiles with clonally expanded populations across heterogeneous phenotypes, demonstrating clonal lineages and phenotypic plasticity. We validate these findings in 31 patients with large B cell lymphoma treated with CD19 CAR T therapy. For these patients, we identify using longitudinal mass-cytometry data an association between NK-like subsets and clinical outcomes at 6 months with both CAR+ and CAR- T cells. These results suggest that non-CAR-derived signals can provide information about patients' immune recovery and be used as correlate of clinically relevant parameters.

Flexible and efficient handling of nanopore sequencing signal data with slow5tools.

Nanopore sequencing is being rapidly adopted in genomics. We recently developed SLOW5, a new file format with advantages for storage and analysis of raw signal data from nanopore experiments. Here we introduce slow5tools, an intuitive toolkit for handling nanopore data in SLOW5 format. Slow5tools enables lossless data conversion and a range of tools for interacting with SLOW5 files. Slow5tools uses multi-threading, multi-processing, and other engineering strategies to achieve fast data conversion and manipulation, including live FAST5-to-SLOW5 conversion during sequencing. We provide examples and benchmarking experiments to illustrate slow5tools usage, and describe the engineering principles underpinning its performance.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes.

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

Fast nanopore sequencing data analysis with SLOW5.

Nanopore sequencing depends on the FAST5 file format, which does not allow efficient parallel analysis. Here we introduce SLOW5, an alternative format engineered for efficient parallelization and acceleration of nanopore data analysis. Using the example of DNA methylation profiling of a human genome, analysis runtime is reduced from more than two weeks to approximately 10.5 h on a typical high-performance computer. SLOW5 is approximately 25% smaller than FAST5 and delivers consistent improvements on different computer architectures.

Chromosome-level de novo genome assembly of Telopea speciosissima (New South Wales waratah) using long-reads, linked-reads and Hi-C.

Telopea speciosissima, the New South Wales waratah, is an Australian endemic woody shrub in the family Proteaceae. Waratahs have great potential as a model clade to better understand processes of speciation, introgression and adaptation, and are significant from a horticultural perspective. Here, we report the first chromosome-level genome for T. speciosissima. Combining Oxford Nanopore long-reads, 10x Genomics Chromium linked-reads and Hi-C data, the assembly spans 823 Mb (scaffold N50 of 69.0 Mb) with 97.8% of Embryophyta BUSCOs "Complete". We present a new method in Diploidocus (https://github.com/slimsuite/diploidocus) for classifying, curating and QC-filtering scaffolds, which combines read depths, k-mer frequencies and BUSCO predictions. We also present a new tool, DepthSizer (https://github.com/slimsuite/depthsizer), for genome size estimation from the read depth of single-copy orthologues and estimate the genome size to be approximately 900 Mb. The largest 11 scaffolds contained 94.1% of the assembly, conforming to the expected number of chromosomes (2n = 22). Genome annotation predicted 40,158 protein-coding genes, 351 rRNAs and 728 tRNAs. We investigated CYCLOIDEA (CYC) genes, which have a role in determination of floral symmetry, and confirm the presence of two copies in the genome. Read depth analysis of 180 "Duplicated" BUSCO genes using a new tool, DepthKopy (https://github.com/slimsuite/depthkopy), suggests almost all are real duplications, increasing confidence in the annotation and highlighting a possible need to revise the BUSCO set for this lineage. The chromosome-level T. speciosissima reference genome (Tspe_v1) provides an important new genomic resource of Proteaceae to support the conservation of flora in Australia and further afield.

Oncogenic cooperation between TCF7-SPI1 and NRAS(G12D) requires ß-catenin activity to drive T-cell acute lymphoblastic leukemia.

Spi-1 Proto-Oncogene (SPI1) fusion genes are recurrently found in T-cell acute lymphoblastic leukemia (T-ALL) cases but are insufficient to drive leukemogenesis. Here we show that SPI1 fusions in combination with activating NRAS mutations drive an immature T-ALL in vivo using a conditional bone marrow transplant mouse model. Addition of the oncogenic fusion to the NRAS mutation also results in a higher leukemic stem cell frequency. Mechanistically, genetic deletion of the ß-catenin binding domain within Transcription factor 7 (TCF7)-SPI1 or use of a TCF/ß-catenin interaction antagonist abolishes the oncogenic activity of the fusion. Targeting the TCF7-SPI1 fusion in vivo with a doxycycline-inducible knockdown results in increased differentiation. Moreover, both pharmacological and genetic inhibition lead to down-regulation of SPI1 targets. Together, our results reveal an example where TCF7-SPI1 leukemia is vulnerable to pharmacological targeting of the TCF/ß-catenin interaction.

Antarctic desert soil bacteria exhibit high novel natural product potential, evaluated through long-read genome sequencing and comparative genomics.

Actinobacteria and Proteobacteria are important producers of bioactive natural products (NP), and these phyla dominate in the arid soils of Antarctica, where metabolic adaptations influence survival under harsh conditions. Biosynthetic gene clusters (BGCs) which encode NPs, are typically long and repetitious high G + C regions difficult to sequence with short-read technologies. We sequenced 17 Antarctic soil bacteria from multi-genome libraries, employing the long-read PacBio platform, to optimize capture of BGCs and to facilitate a comprehensive analysis of their NP capacity. We report 13 complete bacterial genomes of high quality and contiguity, representing 10 different cold-adapted genera including novel species. Antarctic BGCs exhibited low similarity to known compound BGCs (av. 31%), with an abundance of terpene, non-ribosomal peptide and polyketide-encoding clusters. Comparative genome analysis was used to map BGC variation between closely related strains from geographically distant environments. Results showed the greatest biosynthetic differences to be in a psychrotolerant Streptomyces strain, as well as a rare Actinobacteria genus, Kribbella, while two other Streptomyces spp. were surprisingly similar to known genomes. Streptomyces and Kribbella BGCs were predicted to encode antitumour, antifungal, antibacterial and biosurfactant-like compounds, and the synthesis of NPs with antibacterial, antifungal and surfactant properties was confirmed through bioactivity assays.

Draft genome assembly of the invasive cane toad, Rhinella marina.

Background: The cane toad (Rhinella marina formerly Bufo marinus) is a species native to Central and South America that has spread across many regions of the globe. Cane toads are known for their rapid adaptation and deleterious impacts on native fauna in invaded regions. However, despite an iconic status, there are major gaps in our understanding of cane toad genetics. The availability of a genome would help to close these gaps and accelerate cane toad research. Findings: We report a draft genome assembly for R. marina, the first of its kind for the Bufonidae family. We used a combination of long-read Pacific Biosciences RS II and short-read Illumina HiSeq X sequencing to generate 359.5 Gb of raw sequence data. The final hybrid assembly of 31,392 scaffolds was 2.55 Gb in length with a scaffold N50 of 168 kb. BUSCO analysis revealed that the assembly included full length or partial fragments of 90.6% of tetrapod universal single-copy orthologs (n = 3950), illustrating that the gene-containing regions have been well assembled. Annotation predicted 25,846 protein coding genes with similarity to known proteins in Swiss-Prot. Repeat sequences were estimated to account for 63.9% of the assembly. Conclusions: The R. marina draft genome assembly will be an invaluable resource that can be used to further probe the biology of this invasive species. Future analysis of the genome will provide insights into cane toad evolution and enrich our understanding of their interplay with the ecosystem at large.

High risk human papilloma viruses (HPVs) are present in benign prostate tissues before development of HPV associated prostate cancer.

BACKGROUND: Although high risk HPVs are associated with an increased risk of prostate cancer it is not known if they have a causal role. The purpose of this study is to investigate the potential role of human papilloma viruses (HPVs) in prostate cancer. The aims are (i) to investigate the presence and confirm the identity of high risk HPVs in benign prostate tissues prior to the development of HPV positive prostate cancer in the same patients, and (ii) to determine if HPVs are biologically active. METHODS: We used polymerase chain reaction (PCR) to identify HPVs in specimens from 52 Australian men with benign prostate biopsies who 1 to 10 years later developed prostate cancer. Immunohistochemistry (IHC) was used to assess the expression of HPV E7 oncoproteins, cytokeratin and prostate specific antigen (PSA). We used RNASeq data from The Cancer Genome Atlas (TCGA) to identify possible HPV RNA sequences in prostate cancer. RESULTS: HPV screening using standard PCR was conducted on 28 of the 52 sets of benign and later prostate cancers. HPV L1 genes were identified in 13 (46%) benign and 8 (29%) of 28 later prostate cancers in the same patients. HPV E7 genes were identified in 23 (82%) benign and 19 (68%) of 28 subsequent prostate cancers in the same patients. The same HPV types were present in both the benign and subsequent prostate cancers in 9 sets of specimens. HPV type 16 was identified in 15% of benign and 3% of prostate cancers. HPV type 18 was identified in 26% of benign and 16% of prostate cancers. Small numbers of HPV types 45, 47, 76 and 115 were also identified. High confidence RNA-Seq evidence for high risk HPV types 16 and 18 was identified in 12 (2%) of the 502 TCGA prostate cancer transcriptomes. High risk HPV E7 oncoprotein was positively expressed in 23 (82%) of 28 benign prostate specimens but only in 8 (29%) of 28 of the later prostate cancer specimens. This difference is statistically significant (p = 0.001). Prostate specific antigen (PSA) was more highly expressed in 26 (50%) of 52 prostate cancer specimens as compared to prior benign prostate specimens in the same patients. CONCLUSIONS: High risk HPVs are present in benign prostate tissues prior to the development of HPV positive prostate cancer. There is a significantly higher expression of HPV E7 oncoproteins in benign prostate tissues as compared to late prostate cancer that subsequently developed in the same patients. This observation suggests that HPV oncogenic activity is an early phenomenon in a majority of prostate oncogenesis. TCGA RNA-Seq data suggests that HPV is biologically active in some prostate tumour samples.

Small RNA interactome of pathogenic E. coli revealed through crosslinking of RNase E.

RNA sequencing studies have identified hundreds of non-coding RNAs in bacteria, including regulatory small RNA (sRNA). However, our understanding of sRNA function has lagged behind their identification due to a lack of tools for the high-throughput analysis of RNA-RNA interactions in bacteria. Here we demonstrate that in vivo sRNA-mRNA duplexes can be recovered using UV-crosslinking, ligation and sequencing of hybrids (CLASH). Many sRNAs recruit the endoribonuclease, RNase E, to facilitate processing of mRNAs. We were able to recover base-paired sRNA-mRNA duplexes in association with RNase E, allowing proximity-dependent ligation and sequencing of cognate sRNA-mRNA pairs as chimeric reads. We verified that this approach captures bona fide sRNA-mRNA interactions. Clustering analyses identified novel sRNA seed regions and sets of potentially co-regulated target mRNAs. We identified multiple mRNA targets for the pathotype-specific sRNA Esr41, which was shown to regulate colicin sensitivity and iron transport in E. coli Numerous sRNA interactions were also identified with non-coding RNAs, including sRNAs and tRNAs, demonstrating the high complexity of the sRNA interactome.

Transcriptome analysis of human brain tissue identifies reduced expression of complement complex C1Q Genes in Rett syndrome.

BACKGROUND: MECP2, the gene mutated in the majority of Rett syndrome cases, is a transcriptional regulator that can activate or repress transcription. Although the transcription regulatory function of MECP2 has been known for over a decade, it remains unclear how transcriptional dysregulation leads to the neurodevelopmental disorder. Notably, little convergence was previously observed between the genes abnormally expressed in the brain of Rett syndrome mouse models and those identified in human studies. METHODS: Here we carried out a comprehensive transcriptome analysis of human brain tissue from Rett syndrome brain using both RNA-seq and microarrays. RESULTS: We identified over two hundred differentially expressed genes, and identified the complement C1Q complex genes (C1QA, C1QB and C1QC) as a point of convergence between gene expression changes in human and mouse Rett syndrome brain. CONCLUSIONS: The results of our study support a role for alterations in the expression level of C1Q complex genes in RTT pathogenesis.

The Neonatal CD8+ T Cell Repertoire Rapidly Diversifies during Persistent Viral Infection.

CMV is the most common congenital infection in the United States. The major target of congenital CMV is the brain, with clinical manifestations including mental retardation, vision impairment, and sensorineural hearing loss. Previous reports have shown that CD8(+) T cells are required to control viral replication and significant numbers of CMV-specific CD8(+) T cells persist in the brain even after the initial infection has been cleared. However, the dynamics of CD8(+) T cells in the brain during latency remain largely undefined. In this report, we used TCR sequencing to track the development and maintenance of neonatal clonotypes in the brain and spleen of mice during chronic infection. Given the discontinuous nature of tissue-resident memory CD8(+) T cells, we hypothesized that neonatal TCR clonotypes would be locked in the brain and persist into adulthood. Surprisingly, we found that the Ag-specific T cell repertoire in neonatal-infected mice diversified during persistent infection in both the brain and spleen, while maintaining substantial similarity between the CD8(+) T cell populations in the brain and spleen in both early and late infection. However, despite the diversification of, and potential interchange between, the spleen and brain Ag-specific T cell repertoires, we observed that germline-encoded TCR clonotypes, characteristic of neonatal infection, persisted in the brain, albeit sometimes in low abundance. These results provide valuable insights into the evolution of CD8(+) T cell repertoires following neonatal CMV infection and thus have important implications for the development of therapeutic strategies to control CMV in early life.