Karin S. Kassahn

Whole genomes redefine the mutational landscape of pancreatic cancer

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Whole–genome characterization of chemoresistant ovarian cancer

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.

MicroRNAs and their isomiRs function cooperatively to target common biological pathways

BackgroundVariants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules.ResultsTo assess the biological relevance of isomiRs, we have performed ultra-deep miRNA-seq on ten adult human tissues, and created an analysis pipeline called miRNA-MATE to align, annotate, and analyze miRNAs and their isomiRs. We find that isomiRs share sequence and expression characteristics with canonical miRNAs, and are generally strongly correlated with canonical miRNA expression. A large proportion of isomiRs potentially derive from AGO2 cleavage independent of Dicer. We isolated polyribosome-associated mRNA, captured the mRNA-bound miRNAs, and found that isomiRs and canonical miRNAs are equally associated with translational machinery. Finally, we transfected cells with biotinylated RNA duplexes encoding isomiRs or their canonical counterparts and directly assayed their mRNA targets. These studies allow us to experimentally determine genome-wide mRNA targets, and these experiments showed substantial overlap in functional mRNA networks suppressed by both canonical miRNAs and their isomiRs.ConclusionsTogether, these results find isomiRs to be biologically relevant and functionally cooperative partners of canonical miRNAs that act coordinately to target pathways of functionally related genes. This work exposes the complexity of the miRNA-transcriptome, and helps explain a major miRNA paradox: how specific regulation of biological processes can occur when the specificity of miRNA targeting is mediated by only 6 to 11 nucleotides.

PINA v2.0: mining interactome modules

The Protein Interaction Network Analysis (PINA) platform is a comprehensive web resource, which includes a database of unified protein–protein interaction data integrated from six manually curated public databases, and a set of built-in tools for network construction, filtering, analysis and visualization. The second version of PINA enhances its utility for studies of protein interactions at a network level, by including multiple collections of interaction modules identified by different clustering approaches from the whole network of protein interactions (‘interactome’) for six model organisms. All identified modules are fully annotated by enriched Gene Ontology terms, KEGG pathways, Pfam domains and the chemical and genetic perturbations collection from MSigDB. Moreover, a new tool is provided for module enrichment analysis in addition to simple query function. The interactome data are also available on the web site for further bioinformatics analysis. PINA is freely accessible at http://cbg.garvan.unsw.edu.au/pina/.

Animal performance and stress: responses and tolerance limits at different levels of biological organisation

Recent advances in molecular biology and the use of DNA microarrays for gene expression profiling are providing new insights into the animal stress response, particularly the effects of stress on gene regulation. However, interpretation of the complex transcriptional changes that occur during stress still poses many challenges because the relationship between changes at the transcriptional level and other levels of biological organisation is not well understood. To confront these challenges, a conceptual model linking physiological and transcriptional responses to stress would be helpful. Here, we provide the basis for one such model by synthesising data from organismal, endocrine, cellular, molecular, and genomic studies. We show using available examples from ectothermic vertebrates that reduced oxygen levels and oxidative stress are common to many stress conditions and that the responses to different types of stress, such as environmental, handling and confinement stress, often converge at the challenge of dealing with oxygen imbalance and oxidative stress. As a result, a common set of stress responses exists that is largely independent of the type of stressor applied. These common responses include the repair of DNA and protein damage, cell cycle arrest or apoptosis, changes in cellular metabolism that reflect the transition from a state of cellular growth to one of cellular repair, the release of stress hormones, changes in mitochondrial densities and properties, changes in oxygen transport capacities and changes in cardio‐respiratory function. Changes at the transcriptional level recapitulate these common responses, with many stress‐responsive genes functioning in cell cycle control, regulation of transcription, protein turnover, metabolism, and cellular repair. These common transcriptional responses to stress appear coordinated by only a limited number of stress‐inducible and redox‐sensitive transcription factors and signal transduction pathways, such as the immediate early genes c‐fos and c‐jun, the transcription factors NFκB and HIF‐1α, and the JNK and p38 kinase signalling pathways. As an example of environmental stress responses, we present temperature response curves at organismal, cellular and molecular levels. Acclimation and physiological adjustments that can shift the threshold temperatures for the onset of these responses are discussed and include, for example, adjustments of the oxygen delivery system, the heat shock response, cellular repair system, and transcriptome. Ultimately, however, an organism’s ability to cope with environmental change is largely determined by its ability to maintain aerobic scope and to prevent loss in performance. These systemic constraints can determine an organism’s long‐term survival well before cellular and molecular functions are disturbed. The conceptual model we propose here discusses some of the crosslinks between responses at different levels of biological organisation and the central role of oxygen balance and oxidative stress in eliciting these responses with the aim to help the interpretation of environmental genomic data in the context of organismal function and performance.

Evolution of gene function and regulatory control after whole-genome duplication: Comparative analyses in vertebrates

The significance of whole-genome duplications (WGD) for vertebrate evolution remains controversial, in part because the mechanisms by which WGD contributed to functional evolution or speciation are still incompletely characterized. Fish genomes provide an ideal context in which to examine the consequences of WGD, because the teleost lineage experienced an additional WGD soon after divergence from tetrapods and because five teleost genomes are available for comparative analysis. Here we present an integrated approach to characterize these post-duplication genomes based on genome-scale synteny, phylogenetic, temporal, and spatial gene expression and on protein sequence data. A minimum of 3%-4% of protein-coding loci have been retained in two copies in each of the five fish genomes, and many of these duplicates are key developmental genes that function as transcription factors or signaling molecules. Almost all duplicate gene pairs we examined have diverged in spatial and/or temporal expression during embryogenesis. A quarter of duplicate pairs have diverged in function via the acquisition of novel protein domains or via changes in the subcellular localization of their encoded proteins. We compared the spatial expression and protein domain architecture of zebrafish WGD-duplicates to those of their single mouse ortholog and found many examples supporting a model of neofunctionalization. WGD-duplicates have acquired novel protein domains more often than have single-copy genes. Post-WGD changes at the gene regulatory level were more common than changes at the protein level. We conclude that the most significant consequence of WGD for vertebrate evolution has been to enable more-specialized regulatory control of development via the acquisition of novel spatiotemporal expression domains. We find limited evidence that reciprocal gene loss led to reproductive isolation and speciation in this lineage.

Sleeping Beauty mutagenesis reveals cooperating mutations and pathways in pancreatic adenocarcinoma

Pancreatic cancer is one of the most deadly cancers affecting the Western world. Because the disease is highly metastatic and difficult to diagnosis until late stages, the 5-y survival rate is around 5%. The identification of molecular cancer drivers is critical for furthering our understanding of the disease and development of improved diagnostic tools and therapeutics. We have conducted a mutagenic screen using Sleeping Beauty (SB) in mice to identify new candidate cancer genes in pancreatic cancer. By combining SB with an oncogenic Kras allele, we observed highly metastatic pancreatic adenocarcinomas. Using two independent statistical methods to identify loci commonly mutated by SB in these tumors, we identified 681 loci that comprise 543 candidate cancer genes (CCGs); 75 of these CCGs, including Mll3 and Ptk2, have known mutations in human pancreatic cancer. We identified point mutations in human pancreatic patient samples for another 11 CCGs, including Acvr2a and Map2k4. Importantly, 10% of the CCGs are involved in chromatin remodeling, including Arid4b, Kdm6a, and Nsd3, and all SB tumors have at least one mutated gene involved in this process; 20 CCGs, including Ctnnd1, Fbxo11, and Vgll4, are also significantly associated with poor patient survival. SB mutagenesis provides a rich resource of mutations in potential cancer drivers for cross-comparative analyses with ongoing sequencing efforts in human pancreatic adenocarcinoma.

Computational approaches to identify functional genetic variants in cancer genomes

The International Cancer Genome Consortium (ICGC) aims to catalog genomic abnormalities in tumors from 50 different cancer types. Genome sequencing reveals hundreds to thousands of somatic mutations in each tumor but only a minority of these drive tumor progression. We present the result of discussions within the ICGC on how to address the challenge of identifying mutations that contribute to oncogenesis, tumor maintenance or response to therapy, and recommend computational techniques to annotate somatic variants and predict their impact on cancer phenotype.

The mammalian PYHIN gene family: Phylogeny, evolution and expression

BackgroundProteins of the mammalian PYHIN (IFI200/HIN-200) family are involved in defence against infection through recognition of foreign DNA. The family member absent in melanoma 2 (AIM2) binds cytosolic DNA via its HIN domain and initiates inflammasome formation via its pyrin domain. AIM2 lies within a cluster of related genes, many of which are uncharacterised in mouse. To better understand the evolution, orthology and function of these genes, we have documented the range of PYHIN genes present in representative mammalian species, and undertaken phylogenetic and expression analyses.ResultsNo PYHIN genes are evident in non-mammals or monotremes, with a single member found in each of three marsupial genomes. Placental mammals show variable family expansions, from one gene in cow to four in human and 14 in mouse. A single HIN domain appears to have evolved in the common ancestor of marsupials and placental mammals, and duplicated to give rise to three distinct forms (HIN-A, -B and -C) in the placental mammal ancestor. Phylogenetic analyses showed that AIM2 HIN-C and pyrin domains clearly diverge from the rest of the family, and it is the only PYHIN protein with orthology across many species. Interestingly, although AIM2 is important in defence against some bacteria and viruses in mice, AIM2 is a pseudogene in cow, sheep, llama, dolphin, dog and elephant. The other 13 mouse genes have arisen by duplication and rearrangement within the lineage, which has allowed some diversification in expression patterns.ConclusionsThe role of AIM2 in forming the inflammasome is relatively well understood, but molecular interactions of other PYHIN proteins involved in defence against foreign DNA remain to be defined. The non-AIM2 PYHIN protein sequences are very distinct from AIM2, suggesting they vary in effector mechanism in response to foreign DNA, and may bind different DNA structures. The PYHIN family has highly varied gene composition between mammalian species due to lineage-specific duplication and loss, which probably indicates different adaptations for fighting infectious disease. Non-genomic DNA can indicate infection, or a mutagenic threat. We hypothesise that defence of the genome against endogenous retroelements has been an additional evolutionary driver for PYHIN proteins.

Heterologous microarray experiments used to identify the early gene response to heat stress in a coral reef fish

Coral reef fishes are expected to experience rising sea surface temperatures due to climate change. How well tropical reef fishes will respond to these increased temperatures and which genes are important in the response to elevated temperatures is not known. Microarray technology provides a powerful tool for gene discovery studies, but the development of microarrays for individual species can be expensive and time‐consuming. In this study, we tested the suitability of a Danio rerio oligonucleotide microarray for application in a species with few genomic resources, the coral reef fish Pomacentrus moluccensis. Results from a comparative genomic hybridization experiment and direct sequence comparisons indicate that for most genes there is considerable sequence similarity between the two species, suggesting that the D. rerio array is useful for genomic studies of P. moluccensis. We employed this heterologous microarray approach to characterize the early transcriptional response to heat stress in P. moluccensis. A total of 111 gene loci, many of which are involved in protein processing, transcription, and cell growth, showed significant changes in transcript abundance following exposure to elevated temperatures. Changes in transcript abundance were validated for a selection of candidate genes using quantitative real‐time polymerase chain reaction. This study demonstrates that heterologous microarrays can be successfully employed to study species for which specific microarrays have not yet been developed, and so have the potential to greatly enhance the utility of microarray technology to the field of environmental and functional genomics.