Comprehensive characterization of viral integrations and genomic aberrations in HBV-infected intrahepatic cholangiocarcinomas.

BACKGROUND AND AIMS: Despite the epidemiological association between intrahepatic cholangiocarcinoma (iCCA) and HBV infection, little is known about the relevant oncogenic effects. We sought to identify the landscape and mechanism of HBV integration, along with the genomic architecture of HBV-infected iCCA (HBV-iCCA) tumors. APPROACH AND RESULTS: We profiled a cohort of 108 HBV-iCCAs using whole-genome sequencing, deep sequencing, and RNA sequencing, together with preconstructed data sets of HBV-infected HCC (HBV-HCC; n = 167) and combined hepatocellular cholangiocarcinoma (HBV-cHCC/CCA; n = 59), and conventional (n = 154) and fluke-related iCCAs (n = 16). Platforms based on primary iCCA cell lines to evaluate the functional effects of chimeric transcripts were also used. We found that HBV had inserted at multiple sites in the iCCA genomes in 45 (41.7%) of the tumors. Recurrent viral integration breakpoints were found at nine different sites. The most common insertional hotspot (7 tumors) was in the TERT (telomerase reverse transcriptase) promoter, where insertions and mutations (11 tumors) were mutually exclusive, and were accompanied by promoter hyperactivity. Recurrent HBV integration events (5 tumors) were also detected in FAT2 (FAT atypical cadherin 2), and were associated with enrichment of epithelial-mesenchymal transition-related genes. A distinctive intergenic insertion (chr9p21.3), between DMRTA1 (DMRT like family A1) and LINC01239 (long intergenic non-protein coding RNA 1239), had oncogenic effects through activation of the mammalian target of rapamycin (mTOR)/4EBP/S6K pathway. Regarding the mutational profiles of primary liver cancers, the overall landscape of HBV-iCCA was closer to that of nonviral conventional iCCA, than to HBV-HCC and HBV-cHCC/CCA. CONCLUSIONS: Our findings provide insight into the behavior of iCCAs driven by various pathogenic mechanisms involving HBV integration events and associated genomic aberrations. This knowledge should be of use in managing HBV carriers.

A targeted amplicon sequencing panel to simultaneously identify mosquito species and Plasmodium presence across the entire Anopheles genus.

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel's power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control.

Patterns of within-host genetic diversity in SARS-CoV-2.

Monitoring the spread of SARS-CoV-2 and reconstructing transmission chains has become a major public health focus for many governments around the world. The modest mutation rate and rapid transmission of SARS-CoV-2 prevents the reconstruction of transmission chains from consensus genome sequences, but within-host genetic diversity could theoretically help identify close contacts. Here we describe the patterns of within-host diversity in 1181 SARS-CoV-2 samples sequenced to high depth in duplicate. 95.1% of samples show within-host mutations at detectable allele frequencies. Analyses of the mutational spectra revealed strong strand asymmetries suggestive of damage or RNA editing of the plus strand, rather than replication errors, dominating the accumulation of mutations during the SARS-CoV-2 pandemic. Within- and between-host diversity show strong purifying selection, particularly against nonsense mutations. Recurrent within-host mutations, many of which coincide with known phylogenetic homoplasies, display a spectrum and patterns of purifying selection more suggestive of mutational hotspots than recombination or convergent evolution. While allele frequencies suggest that most samples result from infection by a single lineage, we identify multiple putative examples of co-infection. Integrating these results into an epidemiological inference framework, we find that while sharing of within-host variants between samples could help the reconstruction of transmission chains, mutational hotspots and rare cases of superinfection can confound these analyses.

The COVID-19 pandemic has had major health impacts across the globe. The scientific community has focused much attention on finding ways to monitor how the virus responsible for the pandemic, SARS-CoV-2, spreads. One option is to perform genetic tests, known as sequencing, on SARS-CoV-2 samples to determine the genetic code of the virus and to find any differences or mutations in the genes between the viral samples. Viruses mutate within their hosts and can develop into variants that are able to more easily transmit between hosts. Genetic sequencing can reveal how genetically similar two SARS-CoV-2 samples are. But tracking how SARS-CoV-2 moves from one person to the next through sequencing can be tricky. Even a sample of SARS-CoV-2 viruses from the same individual can display differences in their genetic material or within-host variants. Could genetic testing of within-host variants shed light on factors driving SARS-CoV-2 to evolve in humans? To get to the bottom of this, Tonkin-Hill, Martincorena et al. probed the genetics of SARS-CoV-2 within-host variants using 1,181 samples. The analyses revealed that 95.1% of samples contained within-host variants. A number of variants occurred frequently in many samples, which were consistent with mutational hotspots in the SARS-CoV-2 genome. In addition, within-host variants displayed mutation patterns that were similar to patterns found between infected individuals. The shared within-host variants between samples can help to reconstruct transmission chains. However, the observed mutational hotspots and the detection of multiple strains within an individual can make this challenging. These findings could be used to help predict how SARS-CoV-2 evolves in response to interventions such as vaccines. They also suggest that caution is needed when using information on within-host variants to determine transmission between individuals.

Genetic surveillance in the Greater Mekong subregion and South Asia to support malaria control and elimination.

Background: National Malaria Control Programmes (NMCPs) currently make limited use of parasite genetic data. We have developed GenRe-Mekong, a platform for genetic surveillance of malaria in the Greater Mekong Subregion (GMS) that enables NMCPs to implement large-scale surveillance projects by integrating simple sample collection procedures in routine public health procedures. Methods: Samples from symptomatic patients are processed by SpotMalaria, a high-throughput system that produces a comprehensive set of genotypes comprising several drug resistance markers, species markers and a genomic barcode. GenRe-Mekong delivers Genetic Report Cards, a compendium of genotypes and phenotype predictions used to map prevalence of resistance to multiple drugs. Results: GenRe-Mekong has worked with NMCPs and research projects in eight countries, processing 9623 samples from clinical cases. Monitoring resistance markers has been valuable for tracking the rapid spread of parasites resistant to the dihydroartemisinin-piperaquine combination therapy. In Vietnam and Laos, GenRe-Mekong data have provided novel knowledge about the spread of these resistant strains into previously unaffected provinces, informing decision-making by NMCPs. Conclusions: GenRe-Mekong provides detailed knowledge about drug resistance at a local level, and facilitates data sharing at a regional level, enabling cross-border resistance monitoring and providing the public health community with valuable insights. The project provides a rich open data resource to benefit the entire malaria community. Funding: The GenRe-Mekong project is funded by the Bill and Melinda Gates Foundation (OPP11188166, OPP1204268). Genotyping and sequencing were funded by the Wellcome Trust (098051, 206194, 203141, 090770, 204911, 106698/B/14/Z) and Medical Research Council (G0600718). A proportion of samples were collected with the support of the UK Department for International Development (201900, M006212), and Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

Evolution and lineage dynamics of a transmissible cancer in Tasmanian devils.

Devil facial tumour 1 (DFT1) is a transmissible cancer clone endangering the Tasmanian devil. The expansion of DFT1 across Tasmania has been documented, but little is known of its evolutionary history. We analysed genomes of 648 DFT1 tumours collected throughout the disease range between 2003 and 2018. DFT1 diverged early into five clades, three spreading widely and two failing to persist. One clade has replaced others at several sites, and rates of DFT1 coinfection are high. DFT1 gradually accumulates copy number variants (CNVs), and its telomere lengths are short but constant. Recurrent CNVs reveal genes under positive selection, sites of genome instability, and repeated loss of a small derived chromosome. Cultured DFT1 cell lines have increased CNV frequency and undergo highly reproducible convergent evolution. Overall, DFT1 is a remarkably stable lineage whose genome illustrates how cancer cells adapt to diverse environments and persist in a parasitic niche.

Concordance for clonal hematopoiesis is limited in elderly twins.

Although acquisition of leukemia-associated somatic mutations by 1 or more hematopoietic stem cells is inevitable with advancing age, its consequences are highly variable, ranging from clinically silent clonal hematopoiesis (CH) to leukemic progression. To investigate the influence of heritable factors on CH, we performed deep targeted sequencing of blood DNA from 52 monozygotic (MZ) and 27 dizygotic (DZ) twin pairs (aged 70-99 years). Using this highly sensitive approach, we identified CH (variant allele frequency ≥0.5%) in 62% of individuals. We did not observe higher concordance for CH within MZ twin pairs as compared with that within DZ twin pairs, or to that expected by chance. However, we did identify 2 MZ pairs in which both twins harbored identical rare somatic mutations, suggesting a shared cell of origin. Finally, in 3 MZ twin pairs harboring mutations in the same driver genes, serial blood samples taken 4 to 5 years apart showed substantial twin-to-twin variability in clonal trajectories. Our findings propose that the inherited genome does not exert a dominant influence on the behavior of adult CH and provide evidence that CH mutations may be acquired in utero.

Sixteen diverse laboratory mouse reference genomes define strain-specific haplotypes and novel functional loci.

We report full-length draft de novo genome assemblies for 16 widely used inbred mouse strains and find extensive strain-specific haplotype variation. We identify and characterize 2,567 regions on the current mouse reference genome exhibiting the greatest sequence diversity. These regions are enriched for genes involved in pathogen defence and immunity and exhibit enrichment of transposable elements and signatures of recent retrotransposition events. Combinations of alleles and genes unique to an individual strain are commonly observed at these loci, reflecting distinct strain phenotypes. We used these genomes to improve the mouse reference genome, resulting in the completion of 10 new gene structures. Also, 62 new coding loci were added to the reference genome annotation. These genomes identified a large, previously unannotated, gene (Efcab3-like) encoding 5,874 amino acids. Mutant Efcab3-like mice display anomalies in multiple brain regions, suggesting a possible role for this gene in the regulation of brain development.

UTX-mediated enhancer and chromatin remodeling suppresses myeloid leukemogenesis through noncatalytic inverse regulation of ETS and GATA programs.

The histone H3 Lys27-specific demethylase UTX (or KDM6A) is targeted by loss-of-function mutations in multiple cancers. Here, we demonstrate that UTX suppresses myeloid leukemogenesis through noncatalytic functions, a property shared with its catalytically inactive Y-chromosome paralog, UTY (or KDM6C). In keeping with this, we demonstrate concomitant loss/mutation of KDM6A (UTX) and UTY in multiple human cancers. Mechanistically, global genomic profiling showed only minor changes in H3K27me3 but significant and bidirectional alterations in H3K27ac and chromatin accessibility; a predominant loss of H3K4me1 modifications; alterations in ETS and GATA-factor binding; and altered gene expression after Utx loss. By integrating proteomic and genomic analyses, we link these changes to UTX regulation of ATP-dependent chromatin remodeling, coordination of the COMPASS complex and enhanced pioneering activity of ETS factors during evolution to AML. Collectively, our findings identify a dual role for UTX in suppressing acute myeloid leukemia via repression of oncogenic ETS and upregulation of tumor-suppressive GATA programs.

Alcohol and endogenous aldehydes damage chromosomes and mutate stem cells.

Haematopoietic stem cells renew blood. Accumulation of DNA damage in these cells promotes their decline, while misrepair of this damage initiates malignancies. Here we describe the features and mutational landscape of DNA damage caused by acetaldehyde, an endogenous and alcohol-derived metabolite. This damage results in DNA double-stranded breaks that, despite stimulating recombination repair, also cause chromosome rearrangements. We combined transplantation of single haematopoietic stem cells with whole-genome sequencing to show that this damage occurs in stem cells, leading to deletions and rearrangements that are indicative of microhomology-mediated end-joining repair. Moreover, deletion of p53 completely rescues the survival of aldehyde-stressed and mutated haematopoietic stem cells, but does not change the pattern or the intensity of genome instability within individual stem cells. These findings characterize the mutation of the stem-cell genome by an alcohol-derived and endogenous source of DNA damage. Furthermore, we identify how the choice of DNA-repair pathway and a stringent p53 response limit the transmission of aldehyde-induced mutations in stem cells.

Design and Application of Multiplex PCR Seq for the Detection of Somatic Mutations Associated with Myeloid Malignancies.

Targeted sequencing, in which only a selected set of genomic loci are sequenced, enables a much higher coverage of each target than what is obtained using whole genome or exome sequencing. Multiplex PCR offers a simple and affordable technique for specific capture of target regions and can be easily adapted to generate next-generation sequencing (NGS)-ready amplicons. Here we describe a multiplex PCR (MxPCR) approach for capturing 13 leukemia-associated mutation hotspots followed by MiSeq sequencing that enables robust detection of mutations with a variant allele fraction (VAF) as low as 0.8% (0.008) in blood DNA.

Development and validation of a comprehensive genomic diagnostic tool for myeloid malignancies.

The diagnosis of hematologic malignancies relies on multidisciplinary workflows involving morphology, flow cytometry, cytogenetic, and molecular genetic analyses. Advances in cancer genomics have identified numerous recurrent mutations with clear prognostic and/or therapeutic significance to different cancers. In myeloid malignancies, there is a clinical imperative to test for such mutations in mainstream diagnosis; however, progress toward this has been slow and piecemeal. Here we describe Karyogene, an integrated targeted resequencing/analytical platform that detects nucleotide substitutions, insertions/deletions, chromosomal translocations, copy number abnormalities, and zygosity changes in a single assay. We validate the approach against 62 acute myeloid leukemia, 50 myelodysplastic syndrome, and 40 blood DNA samples from individuals without evidence of clonal blood disorders. We demonstrate robust detection of sequence changes in 49 genes, including difficult-to-detect mutations such as FLT3 internal-tandem and mixed-lineage leukemia (MLL) partial-tandem duplications, and clinically significant chromosomal rearrangements including MLL translocations to known and unknown partners, identifying the novel fusion gene MLL-DIAPH2 in the process. Additionally, we identify most significant chromosomal gains and losses, and several copy neutral loss-of-heterozygosity mutations at a genome-wide level, including previously unreported changes such as homozygosity for DNMT3A R882 mutations. Karyogene represents a dependable genomic diagnosis platform for translational research and for the clinical management of myeloid malignancies, which can be readily adapted for use in other cancers.

Leukemia-associated somatic mutations drive distinct patterns of age-related clonal hemopoiesis.

Clonal hemopoiesis driven by leukemia-associated gene mutations can occur without evidence of a blood disorder. To investigate this phenomenon, we interrogated 15 mutation hot spots in blood DNA from 4,219 individuals using ultra-deep sequencing. Using only the hot spots studied, we identified clonal hemopoiesis in 0.8% of individuals under 60, rising to 19.5% of those ≥90 years, thus predicting that clonal hemopoiesis is much more prevalent than previously realized. DNMT3A-R882 mutations were most common and, although their prevalence increased with age, were found in individuals as young as 25 years. By contrast, mutations affecting spliceosome genes SF3B1 and SRSF2, closely associated with the myelodysplastic syndromes, were identified only in those aged >70 years, with several individuals harboring more than one such mutation. This indicates that spliceosome gene mutations drive clonal expansion under selection pressures particular to the aging hemopoietic system and explains the high incidence of clonal disorders associated with these mutations in advanced old age.

Characterization of gene mutations and copy number changes in acute myeloid leukemia using a rapid target enrichment protocol.

Prognostic stratification is critical for making therapeutic decisions and maximizing survival of patients with acute myeloid leukemia. Advances in the genomics of acute myeloid leukemia have identified several recurrent gene mutations whose prognostic impact is being deciphered. We used HaloPlex target enrichment and Illumina-based next generation sequencing to study 24 recurrently mutated genes in 42 samples of acute myeloid leukemia with a normal karyotype. Read depth varied between and within genes for the same sample, but was predictable and highly consistent across samples. Consequently, we were able to detect copy number changes, such as an interstitial deletion of BCOR, three MLL partial tandem duplications, and a novel KRAS amplification. With regards to coding mutations, we identified likely oncogenic variants in 41 of 42 samples. NPM1 mutations were the most frequent, followed by FLT3, DNMT3A and TET2. NPM1 and FLT3 indels were reported with good efficiency. We also showed that DNMT3A mutations can persist post-chemotherapy and in 2 cases studied at diagnosis and relapse, we were able to delineate the dynamics of tumor evolution and give insights into order of acquisition of variants. HaloPlex is a quick and reliable target enrichment method that can aid diagnosis and prognostic stratification of acute myeloid leukemia patients.

Purification and characterization of a novel laccase from Fomitopsis pinicola mycelia.

A novel laccase was isolated from the culture filtrate of the brown-rot fungus, Fomitopsis pinicola. Enzyme production reached its highest level after cultivation for 8 days at 25°C. The enzyme was purified by ultrafiltration, ion exchange chromatography, gelfiltration chromatography, and hydrophobic interaction chromatography. Zymography analysis of the purified enzyme showed a laccase band with a molecular mass of 92 kDa. The molecular weight of the enzyme was 92 kDa as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and gel filtration chromatography. The enzyme also had an isoelectric point of 3.8. The optimum temperature and pH for enzyme activity were 80°C and 3.0, respectively. Enzyme activity was relatively stable in the pH range from 1.5 to 11.0 and at temperatures below 40°C. The N-terminal amino acid sequence of the enzyme was DTHKAEIACRFKDLG. Enzyme activity was potently inhibited by NaN3 and SDS. The enzyme showed the highest specific activity for 2,2-azino-bis(3-ethylthiazoline-6-sulfonate) (ABTS) as a substrate. The Km value of the enzyme for ABTS substrate was 0.28 mM with a Vmax value of 4.5 U/min. The enzyme degraded several recalcitrant dyes at different time intervals during decolorization. Therefore, the novel laccase from F. pinicola may be potentially useful in industry.

Purification and characterization of a novel fibrinolytic enzyme from fruiting bodies of Korean Cordyceps militaris.

A fibrinolytic enzyme has been purified from the fruiting bodies of Korean Cordyceps militaris. The molecular mass of the enzyme was estimated to be 34 kDa by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), fibrin-zymography, and gel filtration chromatography. The 15 amino acid residues of the N-terminal sequence of the enzyme were APVEQCDAPVGLARL, which is dissimilar to those of fibrinolytic enzymes from other mushrooms. Optimal pH and temperature values of the enzyme were 7.0 and 40°C, respectively. The enzyme activity was completely inhibited by phenylmethylsulfonyl fluoride (PMSF), TPCK, 1,10-phenanthroline, Cu(2+), and Ba(2+). It was also significantly inhibited by aprotinin, EDTA, and EGTA. The enzyme showed a higher specificity for a synthetic substrate, N-succinyl-Ala-Ala-Pro-Phe-p-nitroanilide, exhibiting that it is a chymotrypsin-like serine metalloprotease. The enzyme preferentially hydrolyzed the fibrinogen Aα-, followed by the Bß-chains and the γ-chain. The α, ß, and γ-γ chains of fibrin were also degraded by the enzyme.