Catherine A. Shang

Whole-genome landscapes of major melanoma subtypes

Melanoma of the skin is a common cancer only in Europeans, whereas it arises in internal body surfaces (mucosal sites) and on the hands and feet (acral sites) in people throughout the world. Here we report analysis of whole-genome sequences from cutaneous, acral and mucosal subtypes of melanoma. The heavily mutated landscape of coding and non-coding mutations in cutaneous melanoma resolved novel signatures of mutagenesis attributable to ultraviolet radiation. However, acral and mucosal melanomas were dominated by structural changes and mutation signatures of unknown aetiology, not previously identified in melanoma. The number of genes affected by recurrent mutations disrupting non-coding sequences was similar to that affected by recurrent mutations to coding sequences. Significantly mutated genes included BRAF, CDKN2A, NRAS and TP53 in cutaneous melanoma, BRAF, NRAS and NF1 in acral melanoma and SF3B1 in mucosal melanoma. Mutations affecting the TERT promoter were the most frequent of all; however, neither they nor ATRX mutations, which correlate with alternative telomere lengthening, were associated with greater telomere length. Most melanomas had potentially actionable mutations, most in components of the mitogen-activated protein kinase and phosphoinositol kinase pathways. The whole-genome mutation landscape of melanoma reveals diverse carcinogenic processes across its subtypes, some unrelated to sun exposure, and extends potential involvement of the non-coding genome in its pathogenesis.

Growth hormone as a cytokine

1. The growth hormone (GH) receptor was the first of the class 1 cytokine receptors to be cloned. It shares a number of structural characteristics with other family members and common signalling mechanisms based on common usage of the Janus kinase 2 (JAK2).

Bread matters: a national initiative to profile the genetic diversity of Australian wheat

The large and complex genome of wheat makes genetic and genomic analysis in this important species both expensive and resource intensive. The application of next-generation sequencing technologies is particularly resource intensive, with at least 17 Gbp of sequence data required to obtain minimal (1×) coverage of the genome. A similar volume of data would represent almost 40× coverage of the rice genome. Progress can be made through the establishment of consortia to produce shared genomic resources. Australian wheat genome researchers, working with Bioplatforms Australia, have collaborated in a national initiative to establish a genetic diversity dataset representing Australian wheat germplasm based on whole genome next-generation sequencing data. Here, we describe the establishment and validation of this resource which can provide a model for broader international initiatives for the analysis of large and complex genomes.

Sequential Binding of MEIS1 and NKX2-5 on the Popdc2 Gene: A Mechanism for Spatiotemporal Regulation of Enhancers during Cardiogenesis

Summary The homeobox transcription factors NKX2-5 and MEIS1 are essential for vertebrate heart development and normal physiology of the adult heart. We show that, during cardiac differentiation, the two transcription factors have partially overlapping expression patterns, with the result that as cardiac progenitors from the anterior heart field differentiate and migrate into the cardiac outflow tract, they sequentially experience high levels of MEIS1 and then increasing levels of NKX2-5. Using the Popdc2 gene as an example, we also show that a significant proportion of target genes for NKX2-5 contain a binding motif recognized by NKX2-5, which overlaps with a binding site for MEIS1. Binding of the two factors to such overlapping sites is mutually exclusive, and this provides a simple regulatory mechanism for spatial and temporal synchronization of a common pool of targets between NKX2-5 and MEIS1.

Genes induced by growth hormone in a model of adipogenic differentiation.

A substantial number of GH regulated genes have been reported in mature hepatocytes, but genes involved in GH-initiated cell differentiation have not yet been identified. Here we have studied a well-characterised model of GH-dependent differentiation, adipogenesis of 3T3-F442A preadipocytes, to identify genes rapidly induced by GH. Using the suppression subtractive hybridisation technique, we have identified eight genes induced within 60 min of GH treatment, and verified these by northern analysis. Six were identifiable as Stat 2, Stat 3, thrombospondin-1, oncostatin M receptor beta chain, a DEAD box RNA helicase, and muscleblind, a developmental transcription factor. Bioinformatic approaches assigned one of the two remaining unknown genes as a novel 436 residue serine/threonine kinase. As each of the identified genes have important developmental roles, they may be important in initiating GH-induced adipogenesis.

Next-generation sequencing: a challenge to meet the increasing demand for training workshops in Australia

The widespread adoption of high-throughput next-generation sequencing (NGS) technology among the Australian life science research community is highlighting an urgent need to up-skill biologists in tools required for handling and analysing their NGS data. There is currently a shortage of cutting-edge bioinformatics training courses in Australia as a consequence of a scarcity of skilled trainers with time and funding to develop and deliver training courses. To address this, a consortium of Australian research organizations, including Bioplatforms Australia, the Commonwealth Scientific and Industrial Research Organisation and the Australian Bioinformatics Network, have been collaborating with EMBL-EBI training team. A group of Australian bioinformaticians attended the train-the-trainer workshop to improve training skills in developing and delivering bioinformatics workshop curriculum. A 2-day NGS workshop was jointly developed to provide hands-on knowledge and understanding of typical NGS data analysis workflows. The road show–style workshop was successfully delivered at five geographically distant venues in Australia using the newly established Australian NeCTAR Research Cloud. We highlight the challenges we had to overcome at different stages from design to delivery, including the establishment of an Australian bioinformatics training network and the computing infrastructure and resource development. A virtual machine image, workshop materials and scripts for configuring a machine with workshop contents have all been made available under a Creative Commons Attribution 3.0 Unported License. This means participants continue to have convenient access to an environment they had become familiar and bioinformatics trainers are able to access and reuse these resources.

Tumour procurement, DNA extraction, coverage analysis and optimisation of mutation-detection algorithms for human melanoma genomes

Summary Whole genome sequencing (WGS) of cancer patients’ tumours offers the most comprehensive method of identifying both novel and known clinically-actionable genomic targets. However, the practicalities of performing WGS on clinical samples are poorly defined. This study was designed to test sample preparation, sequencing specifications and bioinformatic algorithms for their effect on accuracy and cost-efficiency in a large WGS analysis of human melanoma samples. WGS was performed on melanoma cell lines (n = 15) and melanoma fresh frozen tumours (n = 222). The appropriate level of coverage and the optimal mutation detection algorithm for the project pipeline were determined. An incremental increase in sequencing coverage from 36X to 132X in melanoma tissue samples and 30X to 103X for cell lines only resulted in a small increase (1–2%) in the number of mutations detected, and the quality scores of the additional mutations indicated a low probability that the mutations were real. The results suggest that 60X coverage for melanoma tissue and 40X for melanoma cell lines empower the detection of 98–99% of informative single nucleotide variants (SNVs), a sensitivity level at which clinical decision making or landscape research projects can be carried out with a high degree of confidence in the results. Likewise the bioinformatic mutation analysis methodology strongly influenced the number and quality of SNVs detected. Detecting mutations in the blood genomes separate to the tumour genomes generated 41% more SNVs than if the blood and melanoma tissue genomes were analysed simultaneously. Therefore, simultaneous analysis should be employed on matched melanoma tissue and blood genomes to reduce errors in mutation detection. This study provided valuable insights into the accuracy of SNV with WGS at various coverage levels in human clinical cancer specimens. Additionally, we investigated the accuracy of the publicly available mutation detection algorithms to detect cancer specific SNVs which will aid researchers and clinicians in study design and implementation of WGS for the identification of somatic mutations in other cancers.

Towards an open, collaborative, reusable framework for sharing hands-on bioinformatics training workshops

Abstract There is a clear demand for hands-on bioinformatics training. The development of bioinformatics workshop content is both time-consuming and expensive. Therefore, enabling trainers to develop bioinformatics workshops in a way that facilitates reuse is becoming increasingly important. The most widespread practice for sharing workshop content is through making PDF, PowerPoint and Word documents available online. While this effort is to be commended, such content is usually not so easy to reuse or repurpose and does not capture all the information required for a third party to rerun a workshop. We present an open, collaborative framework for developing and maintaining, reusable and shareable hands-on training workshop content.

Development of a cloud-based Bioinformatics Training Platform

Abstract The Bioinformatics Training Platform (BTP) has been developed to provide access to the computational infrastructure required to deliver sophisticated hands-on bioinformatics training courses. The BTP is a cloud-based solution that is in active use for delivering next-generation sequencing training to Australian researchers at geographically dispersed locations. The BTP was built to provide an easy, accessible, consistent and cost-effective approach to delivering workshops at host universities and organizations with a high demand for bioinformatics training but lacking the dedicated bioinformatics training suites required. To support broad uptake of the BTP, the platform has been made compatible with multiple cloud infrastructures. The BTP is an open-source and open-access resource. To date, 20 training workshops have been delivered to over 700 trainees at over 10 venues across Australia using the BTP.

The cardiac-restricted protein ADP-ribosylhydrolase-like 1 is essential for heart chamber outgrowth and acts on muscle actin filament assembly

Adprhl1, a member of the ADP-ribosylhydrolase protein family, is expressed exclusively in the developing heart of all vertebrates. In the amphibian Xenopus laevis, distribution of its mRNA is biased towards actively growing chamber myocardium. Morpholino oligonucleotide-mediated knockdown of all Adprhl1 variants inhibits striated myofibril assembly and prevents outgrowth of the ventricle. The resulting ventricles retain normal electrical conduction and express markers of chamber muscle differentiation but are functionally inert. Using a cardiac-specific Gal4 binary expression system, we show that the abundance of Adprhl1 protein in tadpole hearts is tightly controlled through a negative regulatory mechanism targeting the 5′-coding sequence of Xenopus adprhl1. Over-expression of full length (40 kDa) Adprhl1 variants modified to escape such repression, also disrupts cardiac myofibrillogenesis. Disarrayed myofibrils persist that show extensive branching, with sarcomere division occurring at the actin-Z-disc boundary. Ultimately, Adprhl1-positive cells contain thin actin threads, connected to numerous circular branch points. Recombinant Adprhl1 can localize to stripes adjacent to the Z-disc, suggesting a direct role for Adprhl1 in modifying Z-disc and actin dynamics as heart chambers grow. Modelling the structure of Adprhl1 suggests this cardiac-specific protein is a pseudoenzyme, lacking key residues necessary for ADP-ribosylhydrolase catalytic activity.