Re-evaluation of the carcinogenic significance of hepatitis B virus integration in hepatocarcinogenesis

To examine the role of hepatitis B virus (HBV) integration in hepatocarcinogenesis, a systematic comparative study of both tumor and their corresponding non-tumor derived tissue has been conducted in a cohort of 60 HBV associated hepatocellular carcinoma (HCC) patients. By using Alu-polymerase chain reaction (PCR) and ligation-mediated PCR, 233 viral-host junctions mapped across all human chromosomes at random, no difference between tumor and non-tumor tissue was observed, with the exception of fragile sites (P = 0.0070). HBV insertions in close proximity to cancer related genes such as hTERT were found in this study, however overall they were rare events. No direct correlation between chromosome aberrations and the number of HBV integration events was found using a sensitive array-based comparative genomic hybridization (aCGH) assay. However, a positive correlation was observed between the status of several tumor suppressor genes (TP53, RB1, CDNK2A and TP73) and the number of chromosome aberrations (r = 0.6625, P = 0.0003). Examination of the viral genome revealed that 43% of inserts were in the preC/C region and 57% were in the HBV X gene. Strikingly, approximately 24% of the integrations examined had a breakpoint in a short 15 nt viral genome region (1820-1834 nt). As a consequence, all of the confirmed X gene insertions were C-terminal truncated, losing their growth-suppressive domain. However, the same pattern of X gene C-terminal truncation was found in both tumor and non-tumor derived samples. Furthermore, the integrated viral sequences in both groups had a similar low frequency of C1653T, T1753V and A1762T/G1764A mutations. The frequency and patterns of HBV insertions were similar between tumor and their adjacent non-tumor samples indicating that the majority of HBV DNA integration events are not associated with hepatocarcinogenesis

Random Forest as a tumour genetic marker extractor

Identifying tumour genetic markers is an essential task for biomedicine. In this thesis, we analyse a dataset of chromosomal rearrangements of cancer samples and present a methodology for extracting genetic markers from this dataset by using a Random Forest as a feature selection tool

Mutational signatures in tumours induced by high and low energy radiation in Trp53 deficient mice.

Ionising radiation (IR) is a recognised carcinogen responsible for cancer development in patients previously treated using radiotherapy, and in individuals exposed as a result of accidents at nuclear energy plants. However, the mutational signatures induced by distinct types and doses of radiation are unknown. Here, we analyse the genetic architecture of mammary tumours, lymphomas and sarcomas induced by high (56Fe-ions) or low (gamma) energy radiation in mice carrying Trp53 loss of function alleles. In mammary tumours, high-energy radiation is associated with induction of focal structural variants, leading to genomic instability and Met amplification. Gamma-radiation is linked to large-scale structural variants and a point mutation signature associated with oxidative stress. The genomic architecture of carcinomas, sarcomas and lymphomas arising in the same animals are significantly different. Our study illustrates the complex interactions between radiation quality, germline Trp53 deficiency and tissue/cell of origin in shaping the genomic landscape of IR-induced tumours

Progressive Mauve: Multiple alignment of genomes with gene flux and rearrangement

Multiple genome alignment remains a challenging problem. Effects of recombination including rearrangement, segmental duplication, gain, and loss can create a mosaic pattern of homology even among closely related organisms. We describe a method to align two or more genomes that have undergone large-scale recombination, particularly genomes that have undergone substantial amounts of gene gain and loss (gene flux). The method utilizes a novel alignment objective score, referred to as a sum-of-pairs breakpoint score. We also apply a probabilistic alignment filtering method to remove erroneous alignments of unrelated sequences, which are commonly observed in other genome alignment methods. We describe new metrics for quantifying genome alignment accuracy which measure the quality of rearrangement breakpoint predictions and indel predictions. The progressive genome alignment algorithm demonstrates markedly improved accuracy over previous approaches in situations where genomes have undergone realistic amounts of genome rearrangement, gene gain, loss, and duplication. We apply the progressive genome alignment algorithm to a set of 23 completely sequenced genomes from the genera Escherichia, Shigella, and Salmonella. The 23 enterobacteria have an estimated 2.46Mbp of genomic content conserved among all taxa and total unique content of 15.2Mbp. We document substantial population-level variability among these organisms driven by homologous recombination, gene gain, and gene loss. Free, open-source software implementing the described genome alignment approach is available from http://gel.ahabs.wisc.edu/mauve .Comment: Revision dated June 19, 200

Twisted Tales: Insights into Genome Diversity of Ciliates Using Single-Cell 'Omics.

The emergence of robust single-cell 'omics techniques enables studies of uncultivable species, allowing for the (re)discovery of diverse genomic features. In this study, we combine single-cell genomics and transcriptomics to explore genome evolution in ciliates (a > 1 Gy old clade). Analysis of the data resulting from these single-cell 'omics approaches show: 1) the description of the ciliates in the class Karyorelictea as "primitive" is inaccurate because their somatic macronuclei contain loci of varying copy number (i.e., they have been processed by genome rearrangements from the zygotic nucleus); 2) gene-sized somatic chromosomes exist in the class Litostomatea, consistent with Balbiani's (1890) observation of giant chromosomes in this lineage; and 3) gene scrambling exists in the underexplored Postciliodesmatophora (the classes Heterotrichea and Karyorelictea, abbreviated here as the Po-clade), one of two major clades of ciliates. Together these data highlight the complex evolutionary patterns underlying germline genome architectures in ciliates and provide a basis for further exploration of principles of genome evolution in diverse microbial lineages

Extensive Genomic Diversity among Bovine-Adapted Staphylococcus aureus: Evidence for a Genomic Rearrangement within CC97

peer-reviewedStaphylococcus aureus is an important pathogen associated with both human and veterinary disease and is a common cause of bovine mastitis. Genomic heterogeneity exists between S. aureus strains and has been implicated in the adaptation of specific strains to colonise particular mammalian hosts. Knowledge of the factors required for host specificity and virulence is important for understanding the pathogenesis and management of S. aureus mastitis. In this study, a panel of mastitis-associated S. aureus isolates (n = 126) was tested for resistance to antibiotics commonly used to treat mastitis. Over half of the isolates (52%) demonstrated resistance to penicillin and ampicillin but all were susceptible to the other antibiotics tested. S. aureus isolates were further examined for their clonal diversity by Multi-Locus Sequence Typing (MLST). In total, 18 different sequence types (STs) were identified and eBURST analysis demonstrated that the majority of isolates grouped into clonal complexes CC97, CC151 or sequence type (ST) 136. Analysis of the role of recombination events in determining S. aureus population structure determined that ST diversification through nucleotide substitutions were more likely to be due to recombination compared to point mutation, with regions of the genome possibly acting as recombination hotspots. DNA microarray analysis revealed a large number of differences amongst S. aureus STs in their variable genome content, including genes associated with capsule and biofilm formation and adhesion factors. Finally, evidence for a genomic arrangement was observed within isolates from CC97 with the ST71-like subgroup showing evidence of an IS431 insertion element having replaced approximately 30 kb of DNA including the ica operon and histidine biosynthesis genes, resulting in histidine auxotrophy. This genomic rearrangement may be responsible for the diversification of ST71 into an emerging bovine adapted subgroup

Concordance of copy number abnormality detection using SNP arrays and Multiplex Ligation-dependent Probe Amplification (MLPA) in acute lymphoblastic leukaemia

In acute lymphoblastic leukaemia, MLPA has been used in research studies to identify clinically relevant copy number abnormality (CNA) profiles. However, in diagnostic settings other techniques are often employed. We assess whether equivalent CNA profiles are called using SNP arrays, ensuring platform independence. We demonstrate concordance between SNP6.0 and MLPA CNA calling on 143 leukaemia samples from two UK trials; comparing 1,287 calls within eight genes and a region. The techniques are 99% concordant using manually augmented calling, and 98% concordant using an automated pipeline. We classify these discordant calls and examine reasons for discordance. In nine cases the circular binary segmentation (CBS) algorithm failed to detect focal abnormalities or those flanking gaps in IKZF1 probe coverage. Eight cases were discordant due to probe design differences, with focal abnormalities detectable using one technique not observable by the other. Risk classification using manually augmented array calling resulted in four out of 143 patients being assigned to a different CNA risk group and eight patients using the automated pipeline. We conclude that MLPA defined CNA profiles can be accurately mirrored by SNP6.0 or similar array platforms. Automated calling using the CBS algorithm proved successful, except for IKZF1 which should be manually inspected