2014 AGTA Conference Handbook
2014 AGTA Conference Handbook
2014 AGTA Conference Handbook
2014 AGTA
Conference
2
Handbook
2014 AGTA
Conference
CONTENTS
WELCOME 5
GENERAL INFORMATION 8
AGTA14 CONFERENCE PROGRAM 11
CONFERENCE SOCIAL PROGRAM 20
Abstracts and Biographies 21
POSTER PRESENTATIONS – Monday 59
POSTER PRESENTATIONS – Tuesday 62
SPONSORS 119
EXHIBITORS 121
EXHIBITION FLOOR PLAN 129
AGTA 2014 DELEGATE LIST 130
3
Handbook
2014 AGTA
Conference
AGTA
Executive Team
Professor Erik (Rik) Thompson, School of Ms Vivien Vasic, MHTP Medical Genomics
Biomedical Sciences and Queensland Institute Facility, Monash Health Translation Precinct,
of Technology, Brisbane Melbourne
4
Handbook
2014 AGTA
Conference
conference logo.
The AGTA conference enables the
opportunity for interaction between biologists,
bioinformaticians and technologists. This
unique mix is one of the reasons that the
Australian genomics community has a dynamic
cross-disciplinary and innovative approach
It is our great pleasure to welcome you to the
to genomic analysis, and is at the forefront of
14th annual conference of the Australasian
analysis tools for new types of ‘omics data. We
Genomic Technologies Association (AGTA)
have invited an outstanding list of international
at the Crown Promenade Hotel, Melbourne,
and national speakers plus have engaged
Australia, Sunday 12th to Wednesday 15th
industry to showcase the best in new genomic
October 2014.
technologies. We hope you make the most of
The AGTA conference (formerly known as the opportunity for some great networking and
AMATA) is Australia’s foremost genomic cutting edge science in the heart of the “most
technology conference. This year we have livable city” in the world.
worked hard to bring you an exciting program
We are excited to be holding the conference
that spans a diversity of research topics utilising
in Melbourne, one of the centres of Australian
cutting edge genomic technologies. Major
genomics research. We hope you enjoy the
themes include clinical genomics, functional
conference!
genomics, bioinformatics, metagenomics,
cancer genomics, epigenomics and plant
genomics. We have also scheduled a Alicia Oshlack
workshop focusing on the application of Murdoch Children’s Research Institute,
clinical genomics, an exciting area currently Melbourne
undergoing a major expansion in Australia.
Richard Tothill
Further features of our program include new
Peter MacCallum Cancer Centre, Melbourne
student social functions and our regular
conference dinner to be held at the Arts
Centre of Melbourne under the iconic
spire that has featured as part of this year’s
CONFERENCE MANAGERS
Leishman Associates
113 Harrington Street, Hobart TAS 7000
170 Elgin Street, Carlton VIC 3053
P. 03 6234 7844 F. 03 6234 5958
E. brigitte@leishman-associates.com.au
W. www.leishman-associates.com.au
6
Handbook
THE NEW CyTOF®2 MASS CYTOMETER
Discover more. Imagine more.
Mass Cytometry. The Future of Cytometry Today
Mass cytometry resolves multiple metal probes per cell with minimal signal
overlap, maximizing the information obtained from a single sample. With over
100 detector channels and over 30 isotopic probes, mass cytometry provides
researchers with an unparalleled ability to phenotypically and functionally
profile cells from normal and diseased states.
The DVS line of reagents, instrumentation, and data analysis tools provide all
the tools necessary for novel discovery applications using mass cytometry,
including:
®
www.mscience.com.au
sales@mscience.com.au • 1800 678 242
2014 AGTA
Conference
GENERAL
INFORMATION
8
Handbook
2014 AGTA
Conference
All delegates who are registered to attend For information on the chairperson attending
the dinner will receive a named sticker at your session, please see the Registration Desk.
registration. You MUST place your sticker on A technician will be present in the speaker’s
a table located on poster boards next to the Preparation Room during registration hours.
Registration Desk. You must allocate yourself There will be facility to test and modify your
to a table no later than 11.00 am Monday 13 presentation as required.
October 2014.
Oral Presentations
Student Functions
Please refer to the program for the time
All conference students and early career allocated for each presentation, as these do
researchers are invited to the meet-a-mentor vary. The chairperson for your session will give
session and casual function on Monday 13 you a 3 minute warning, however you are asked
October, however registration is required as to adhere to your time allocation so that the
spaces are limited. Please forward enquiries program remains on schedule.
to the Registration Desk where bookings
can be made subject to availability. Further
information about these events can be found
Poster Presentations
on pages 20. Posters will be displayed in the Trade Exhibition
Proudly supported by: Area for the duration of the conference. Poster
Session One will be held on Monday 13 October
2014 from 1.30pm to 3.15pm. Poster Session Two
will be held on Tuesday 14 October 2014 from
1.00pm to 3.00pm. Further information on the
Poster Sessions can be found on pages 58.
9
Handbook
Historical Breakthroughs
Breathtaking Progress
Radical Reinvention
Introducing NextSeq™
A Whole Human Genome on Your Desktop
SUNDAY
12 October 2014
Don’t forget
to get your Exhibitor
Passport stamped
to go into the draw
and win some great
prizes – including
an ipad mini!
MONDAY
13 October 2014
MONDAY
13 October 2014
TUESDAY
14 October 2014
Promenade 1
SESSION 4: CANCER GENOMICS
Chairs: Dr Mark Waltham and Dr Ian Majewski
0830 – 0915 Interrogating the architecture of cancer genomes
Dr Peter Campbell
Wellcome Trust Sanger Institute, Cambridge, UK
0915 – 0945 THE GENOMIC LANDSCAPE OF PRIMARY AND ACQUIRED RESISTANT
HIGH GRADE SEROUS OVARIAN CANCER
Professor David Bowtell
Peter MacCallum Cancer Centre, VIC
0945 – 1015 CIRCULATING TUMOUR DNA AS A LIQUID BIOPSY IN CANCER
Dr Sarah-Jane Dawson
Peter MacCallum Cancer Centre, VIC
1015 – 1030 A novel long noncoding RNA, lncUSMycN, promotes tumourigenesis
by binding to the RNA-binding protein NonO and up-regulating MYCN
oncogene expression
Dr Tao Liu
Histone Modification Group, Children’s Cancer Institute Australia,
NSW
1030 – 1045 Patterns of clonal evolution involved in treatment resistance in
Diffuse Large B-cell Lymphoma using tumour and plasma sequencing
Dr Ryan Morin
Simon Fraser University, Canada
1045 – 1115 Morning Refreshments & Trade Exhibition Promenade Foyer &
Promenade 2 & 3
TUESDAY
14 October 2014
Promenade 1
TUESDAY
14 October 2014
WEDNESDAY
15 October 2014
Promenade 1
Promenade 1
WEDNESDAY
15 October 2014
• High impact – ranked 6th among research journals in Genetics and Heredity
• Open access for all research
• Fast peer review and rapid publication after acceptance
• High visibility
• Inclusion in PubMed, PubMed Central, Medline and Web of Science
biomedcentral.com
A09435
2014 AGTA
Conference
CONFERENCE
SOCIAL PROGRAM
20
Handbook
2014 AGTA
Conference
Abstracts and
Biographies
21
Handbook
2014 AGTA
Conference
22
Handbook
2014 AGTA
Conference
Abstract
Genetic analysis in large structured tree
MONDAY populations using next-generation genomics
13 October 2014 technologies provides a powerful approach
to dissect genetic variation affecting
regulatory networks and biosynthetic pathways
underlying wood development in Eucalyptus
trees. We have performed comprehensive
genome mapping, transcriptome sequencing,
SESSION 1 metabolic profiling and wood cell wall analyses
in an interspecific F2 backcross population of
E. grandis x E. urophylla. Transcriptome profiling
(Illumina RNA-seq) of developing wood of 283
PLANT GENOMICS F2 progeny allowed genome-wide mapping of
expression QTLs (eQTLs) for more than 20,000
CHAIRS: CARSTEN KULHEIM xylem expressed genes. Several trans-eQTL
AND NOEL COGAN hotspots in both backcrosses mark the location
of polymorphisms affecting secondary cell
wall (SCW) biosynthetic genes and provide
0845 – 0930 evidence for segregating components of
the transcriptional network regulating SCW
formation. We further analysed the genetic
POPULATION GENOMICS UNRAVELS GENETIC architecture of transcript level variation for
DIVERSITY AND REGULATION OF GROWTH 353 xylem expressed transcription factors (TFs)
AND DEVELOPMENT IN EUCALYPTUS and show that shared eQTLs and expression
Professor Zander Myburg correlation of TF genes can be used to
UNIVERSITY OF PRETORIA, SOUTH AFRICA reconstruct modules of the SCW transcriptional
network and identify candidates affecting key
Biography
biological processes in xylem development.
Zander Myburg is a professor enetics at the
We also analysed sequence variation in the
University of Pretoria (UP) and holds the Chair
transcriptomes of the F2 progeny and identified
in Forest Genomics and Biotechnology at UP.
266,433 SNPs and 9,682 indels segregating
His research programme in the Forestry and
in the backcross progeny, half of which are
Agricultural Biotechnology Institute (FABI) and
predicted to affect protein sequence and
Genomics Research Institute (GRI) focuses
structure and thereby potentially contribute
on the genomics and molecular genetics
to phenotypic variation. A total of 1,183
of wood development in fast-growing forest
predicted loss-of-function (LoF) mutations are
trees and, in particular, the genetic regulation
polymorphic in the backcross and inactivate
of cellulose biosynthesis in trees. His research
alleles of 1,285 protein-coding genes including
group is pioneering the use of systems genetics
Eucalyptus orthologs of cell wall biosynthetic
approaches to unravel the genetic control
genes. Together these analyses provide the
of wood formation in Eucalyptus trees. He
foundation for gaining a systems genetics
has also been the lead investigator of the
understanding of the regulation of growth and
US Department of Energy (DOE) funded
development in Eucalyptus.
Eucalyptus Genome Project. He has supervised
30 postgraduate (MSc and PhD) students and
is author of 47 peer-reviewed papers and
book chapters in the field of plant molecular
genetics and genomics.
23
Handbook
2014 AGTA
Conference
24
Handbook
2014 AGTA
Conference
which define the RNA sequence recognised by The exotic fungal pathogen Myrtle rust
the protein. This recognition code has recently (Puccinia psidii) poses a serious biosecurity
been “cracked” (Barkan et al, 2012) allowing threat to ecosystems and forest-based
binding sites of PPRs to be predicted and then industries across the world, including Australia.
tested experimentally. In addition, reanalysis More than 140 Australian native species in
of data generated from small RNA sequencing the family Myrtaceae are susceptible to the
experiments supports the idea that short non- rust disease, including eucalypts (Eucalyptus
coding RNA fragments that accumulate in spp.) and tea trees (Melaleuca spp.). The rust
chloroplasts are likely to be “footprints” of RNA pathogen originated in South and Central
binding proteins (Ruwe & Schmitz-Linneweber, America, and has spread the length of the
2012). These two recent discoveries have the east coast since it was discovered in Australia
potential to rapidly accelerate our ability to in April 2010. The rust affects young leaves,
match PPRs with their RNA targets. flowers and fruits, causing lesions on the host
plant. Finding genetic markers for resistance
We have developed a method to analyse
against Myrtle rust in Myrtaceae may help
putative PPR “footprints” by sequencing small
identify resistant plants and provide a means
RNAs derived from wildtype and mutant plants.
of maintaining production in Myrtaceae-
As a proof-of-concept we have been able to
related rural industries (e.g. oils and culinary
show for 2 PPR proteins, SOT1 and EMB2564,
products such as lemon myrtle) as well as
that putative PPR footprints are lost or reduced
helping nurseries and contributing to restoring
in the mutant plants and that these observed
ecosystems.
changes are consistent with the PPR binding
sites predicted using computational methods. Understanding the molecular basis of
resistance and identifying potential genetic
markers of pathogen resistance in important
in plant breeding programs and we aim
1015 – 1030 to produce molecular markers that will
assist breeding of rust-resistant varieties. We
sequenced the transcriptomes of resistant and
DISCOVERING THE MOLECULAR BASIS OF
susceptible individuals of Melaleuca alternifolia
RESISTANCE IN AUSTRALIAN MYRTACEAE TO
THE EXOTIC FUNGAL PATHOGEN MYRTLE RUST and M. quinquenervia before and five days
(PUCCINIA PSIDII) post-inoculation with Myrtle rust spores. A large
proportion of the assembled transcripts showed
Ms Ji-Fan (Sarah) Hsieh differential expression. We investigated the
THE AUSTRALIAN NATIONAL UNIVERSITY, ACT GO categories of the differentially expressed
Biography transcripts and found many defence-related
Ji-Fan (Sarah) Hsieh is a PhD student in the categories up-regulated in the resistant
Research School of Biology at the Australian plant. One gene family that showed marked
National University. Her research focuses on differences between the resistant and
discovering of the molecular basis of resistance susceptible individuals was the nucleotide-
in Australian Myrtaceae to the exotic plant binding site leucine-rich repeat family (NBS-
pathogen Myrtle rust, and developing LRR), which is responsible for the recognition
molecular markers to help identifying resistant of pathogens and is the start of the signal
plants, providing a means of maintaining cascade that leads to the production of
production in Myrtaceae-related rural industries defence-related proteins and plant secondary
(e.g. oils and culinary products such as metabolites. We are currently developing
lemon myrtle) as well as helping nurseries and molecular markers for resistance based on this
contributing to restoring ecosystems. gene family.
25
Handbook
2014 AGTA
Conference
26
Handbook
2014 AGTA
Conference
Abstract
Cancer involves changes in gene expression
through a variety of mechanisms, including
SESSION 2 translocations and epigenetic alterations.
Recent genome-wide approaches have
revealed a plethora of distal regulatory
elements, some of which have been
EPIGENETICS associated with cancer initiation and
progression. Chromatin interactions, which
CHAIRS: EMMA WHITELAW are two or more distant genomic regions that
AND DAVID MARTINO come together in close spatial proximity in the
3-Dimensional environment of the cell nucleus,
may facilitate cancer through mechanisms
1115 – 1200 such as transcriptional control by connecting
distal regulatory elements with target genes.
3-DIMENSIONAL GENOME ORGANIZATION In this talk, I will present Chromatin Interaction
AND TRANSCRIPTIONAL CONTROL Analysis with Paired-End Tags (ChIA-PET)
sequencing, and present my lab’s recent
Professor Melissa Fullwood advances in refining the ChIA-PET method,
YALE-NUS COLLEGE, SINGAPORE analysis, and application of this method to
Biography better understand transcriptional control in
Dr Melissa J. Fullwood is a Junior Principal cancer cell lines, as well as future goals in terms
Investigator in the Cancer Science Institute with of translation of our results to the clinic.
a joint appointment as an Assistant Professor
at Yale-NUS and a Joint Principal Investigator
position at the Institute of Molecular and Cell 1200 – 1230
Biology, A*STAR, Singapore. Her research
focuses on chromatin interactions and
transcription in cancer. She completed her SHRNA SCREENS FOR NOVEL EPIGENETIC
undergraduate degree in Biological Sciences MODIFIERS AND CHARACTERISATION OF
in 2005 at Stanford University, and her PhD THEIR ROLES IN X INACTIVATION
at the Genome Institute under the NUS
Dr Marnie Blewitt
Graduate School for Integrative Sciences and
WALTER AND ELIZA HALL INSTITUTE, VIC
Engineering in 2009. She was a Lee Kuan Yew
Post-Doctoral Fellow in Duke-NUS. She has 17 Biography
publications, which have been cited over 1500 Marnie Blewitt performed her PhD with Emma
times, and 2 patents. She was one of three Whitelaw developing a mutagenesis screen for
recipients of the inaugural L’Oreal-UNESCO epigenetic modifiers in the mouse, graduating
for Women in Science National Fellowships in in 2005. Marnie worked with Douglas Hilton at
Singapore in 2009, and was the international The Walter and Eliza Hall Institute for her post-
winner of the GE and Science prize in 2010. doctoral studies, where she identified the role
She serves on the editorial board for Scientific ofthe novel protein Smchd1 in X inactivation,
Reports. In 2013, she received the National and studied the function of polycomb group
Research Foundation (NRF) fellowship, which proteins in hematopoietic stem cells. This
comes with a S$3.5 million grant over 5 years. earned her the AAS Gani medal and the
L’Oreal Australia Women in Science fellowship
2009. In 2010, Marnie established her own
group at The Walter and Eliza Hall Institute as
an ARC QEII fellow, to study the molecular
mechanisms of epigenetic control.
27
Handbook
2014 AGTA
Conference
Abstract
Epigenetic marks are frequently reported to 1230 – 1245
correlate with gene silencing, yet there are
few circumstances when the mechanistic role
of specific epigenetic marks in the process CHARACTERISATION OF THE NOVEL
EPIGENETIC MODIFIER RLF
of inactivation have been characterised. X
chromosome inactivation (XCI), the dosage Dr Harald Oey
compensation mechanism in female mammals, LA TROBE UNIVERSITY, VIC
provides a powerful system where the initiation,
Biography
establishment and maintenance of silencing
Harald completed his university studies at
can be studied for hundreds of genes in
Griffith University in 2007. he then worked
parallel. Each phase of silencing can be
for Professor John Mattick at the IMB in
monitored in vitro or in vivo, as female cells
Brisbane for two years before taking on a
transition from two active X chromosomes
postdoctoral research position with Professor
to one active and one inactive. While the
Emma Whitelaw, first at the QIMR in Brisbane
ontogeny of X chromosome inactivation (XCI)
and subsequently at La Trobe University in
is well described, much remains unknown
Melbourne.
about the molecular mechanisms governing
this process. Therefore we are using a bespoke Abstract
shRNA library to screen for epigenetic For many years our lab has been running
regulators of XCI, as in order to understand an ENU mutagenesis screen for modifiers
more about epigenetic silencing, we must first of transgene silencing in the mouse. The
identify all the key players. screen utilizes a GFP transgene expressed
in erythrocytes in a variegated manner
We are using our targeted library in multiple
(expressed in some erythrocytes, silent in
XCI screens: to identify epigenetic regulators
others), and genes that enhance or supress
involved in initiating, establishing or maintaining
variegation are called modifiers of murine
XCI. For high throughput screens, we need a
metastable epialleles (MommeD). Recently,
rapid system to survey XCI. We have produced
we reported the initial characterisation
2 knockin alleles of the X-linked house-keeping
of three novel enhancers of variegation,
gene Hprt: one which produces GFP and the
MommeD8, MommeD28 and MommeD34,
other mCherry under the control of the Hprt
all of which have mutations in a gene called
promoter. We have used these alleles and
Rearranged L-myc fusion (Rlf). Interestingly,
identified the H3K9 methyltransferase Setdb1
while several MommeD alleles cause reduced
as a new regulator of the maintenance of X
DNA methylation at the GFP transgene, and
inactivation. Using a combination of allele-
an increase in GFP expression, Rlf is currently
specific RNA-seq, ChIP-seq and RRBS we have
the only known MommeD causing increased
placed Setdb1 within the epigenetic hierarchy
methylation and a concomitant decrease in
of XCI, which has revealed the broader
GFP expression. To investigate the function
significance of H3K9 methylation in epigenetic
of Rlf further we have carried out a range of
silencing.
experiments to characterise the effect of Rlf
genome-wide. These include whole genome
bisulfite sequencing (~30x coverage) and
mRNA-seq in multiple mutant and wildtype
embryonic tissues, ChIP-seq to identify Rlf DNA
binding sites and co-immunoprecipitation,
followed by mass spectrometry, to identify
Rlf binding partners in-vivo. We find that Rlf
interacts with a histone demethylase in a
28
Handbook
2014 AGTA
Conference
29
Handbook
2014 AGTA
Conference
1600 – 1630
SESSION 3
The STATegra project: new statistical
tools for analysis and integration of
diverse omics data.
30
Handbook
2014 AGTA
Conference
Abstract
Next generation sequencing has speed
1630 – 1645
up genome analysis and brought omics
research closer to many organisms and CLASSIFICATION OF GENOMIC FUSIONS INTO
biological scenarios. Today an increasing STRUCTURAL VARIATION EVENTS
number of research projects propose the
combined use of different omics platforms Dr Jan Schroeder
to investigate diverse aspects of genome WALTER AND ELIZA HALL INSTITUTE, VIC
functioning. These proposals ideally seek to Abstract
provide complementary sources of molecular A precise understanding of the structural
information that eventually can be put variations (SVs) between samples of DNA is
together to obtain systems biology models important in the study of population diversity
of biological processes. Hence, it is not rare and disease causing mutations. There are
anymore to find experimental designs involving numerous methods that are designed to
the collection of genome, transcriptome, detect SVs in DNA sequencing data. Most of
epigenome and even metabolome data these methods do not provide a higher-level
on a particular system. However, standard interpretation of variants, instead identifying
methodologies for the integration of diverse a large set of coordinate pairs that are
omics data types are not yet ready and connected via candidate fusions – so-called
researchers frequently face post-experiment break points. While some SV callers include a
question on how to combine data of different description of underlying genomic events, they
nature, variability, and significance into an are currently limited to a set of simple variant
analysis routine that sheds more light than types (such as deletions or inversions).
the analysis of individual datasets separately.
Here we present an algorithm to process
The STATegra project has been conceived
the output of SV calling algorithms and infer
to address these problems and provide the
the rearrangement events that explain the
genomics community with user-friendly tools for
observed fusions. Our method creates a graph
the integration of different omics data types.
data structure from the fusion information
STATegra targets several sequencing based
and looks for patterns that are characteristic
functional genomics methods, proteomics
of more complex rearrangement types (such
and metabolomics. In this presentation I will
as balanced translocations). Our algorithm
describe current results of the project including
summarises fusion events that are seemingly
the STATegraEMS, an experiment management
independent in the original output into single,
system for storage and annotation of complex
more complex events (for example, an intra-
omics experiments, experimental design issues,
chromosomal translocation involves three
novel integrative visualization tools, statistical
fusions) thus reducing the output to a more
approaches to integrate RNA-seq data with
intelligible set. A better-categorised output
other omics technologies, and data mining
allows for better filtering of specific events
strategies to leverage public domain datasets
that are most relevant to the experimental
in the integrative effort. I will also present the
circumstances.
STATegRa, a new Bioconductor R package for
integrative omics data analysis. We demonstrate that our method can improve
the output of various SV algorithms (such as
CREST, BreakDancer, Delly, and Socrates),
identifying redundant events, and reducing
false positives. Augmenting SV calls in this
manner simplifies exploratory SV analysis of
samples in various experimental settings,
such as whole genome data, exome data, in
somatic or germline tissue.
31
Handbook
2014 AGTA
Conference
class of genomic lesion is a fusion gene - Liam graduated with a Bachelor of Science
whereby two genes are joined as a result (Biochem/Biomed) at the University of
of a translocation, producing a new and Queensland in 2007 and then completed his
potentially oncogenic gene. Many fusions Masters in Bioinformatics at the University of
are recurrent in patient tumour populations Queensland in 2014. Currently employed by SA
and their identification has both research pathology as a genetic informatician.
and clinical benefits. RNA sequencing (RNA- Abstract
Seq) is a powerful method for profiling cancer Next-generation sequencing (NGS) has
transcriptomes to identify fusions. However, reached a price point where it is viable for
many methods for identifying fusions have low clinical diagnostic application, however,
precision, require tuning of parameters or are the development of bioinformatics pipelines
not suitable for all RNA-Seq read lengths and remains an active area of research. We have
layouts. developed a read alignment and variant
Here, we present JAFFA (https://code. calling pipeline based on BWA-MEM and
google.com/p/jaffa-project/), a method that the GATK best practices and have adapted
overcomes these shortfalls. JAFFA is based on it to suit amplicon and enrichment based
the idea of comparing the sequenced cancer protocols for sequencing on the MiSeq. We
transcriptome to a reference transcriptome, have compared its performance against the
such as GENCODE. For short reads, the vendor pipeline, MiSeq Reporter, by assessing
cancer transcriptome is built through de novo recall rate and specificity using retrospective
32
Handbook
2014 AGTA
Conference
34
Handbook
2014 AGTA
Conference
35
Handbook
2014 AGTA
Conference
Abstract
Patients with neuroblastoma, due to the 1030 – 1045
amplification of a 130 kb genomic DNA region
containing the MYCN oncogene, have a
PATTERNS OF CLONAL EVOLUTION INVOLVED
survival rate of less than 10%. While MYCN
IN TREATMENT RESISTANCE IN DIFFUSE LARGE
has been extensively studied in the last three
B-CELL LYMPHOMA USING TUMOUR AND
decades, it is unknown whether other genetic
PLASMA SEQUENCING
elements within the 130 kb amplicon contribute
to tumorigenesis. The novel long noncoding Dr Ryan Morin
RNA up-stream of MYCN (lncUSMycN) in SIMON FRASER UNIVERSITY, CANADA
the 130 kb amplicon has been recently Biography
manually annotated by the Human and Dr Morin is an Assistant Professor, co-appointed
Vertebrate Analysis and Annotation (HAVANA) in the department of Molecular Biology and
bioinformatics team, but has not been Biochemistry at Simon Fraser University and
studied. Using 5’- and 3’- rapid amplification the BC Cancer Agency's Genome Sciences
of cDNA ends PCR, we experimentally Centre. He earned both M.Sc. and Ph.D.
identified 140bp extra nucleotides at the degrees from the University of British Columbia
5’- end of lncUSMycN. Analysis of 341 human (Canada) in the laboratory of Marco Marra.
neuroblastoma samples showed that the Dr Morin is best known for his early work using
lncUSMycN gene was co-amplified with the massively parallel sequencing to identify
MYCN oncogene. RNA-binding protein pull- driver mutations in leukemias and lymphomas
down assays and mass spectrometry analysis including histone modifiers such as EZH2 and
identified the RNA-binding protein NonO MLL2. Some of his recent recent work has been
as a lncUSMycN RNA binding protein, and published in world-class journals such as the
RNA-immunoprecipitation assays showed New England Journal of Medicine, Nature
NonO bound to both lncUSMycN and Genetics, Cancer Cell and Nature. He has
MYCN RNAs. While knocking-down NonO contributed to 42 peer-reviewed articles and
and lncUSMycN with small interfering RNAs 5 reviews/book chapters on various aspects
reduced MYCN RNA expression, forced over- of genomics and DNA sequence analysis and
expression of lncUSMycN resulted in MYCN his published work has been collectively cited
RNA up-regulation by binding to NonO. In 476 over 5500 times.
human neuroblastoma samples, high levels
Abstract
of lncUSMycN and NonO RNA expression
DLBCL is a common aggressive non-Hodgkin
correlated with high levels of MYCN RNA
lymphoma (NHL) and although sequencing
expression and independently predicted
diagnostic specimens has identified many
poor patient prognoses. Moreover, treatment
driver genes, it is unclear whether additional
with antisense oligonucleotides targeting
genes are involved in relapse or resistance
lncUSMycN in neuroblastoma-bearing mice
to standard treatments (R-CHOP). We
significantly reduced MYCN expression
performed a meta-analysis of all published
and blocked tumor progression. Our data
data interrogating the genomes of DLBCLs.
therefore demonstrate the important role of
To determine the extent to which clonal
the novel long noncoding RNA lncUSMycN in
evolution contributes to treatment resistance,
regulating MYCN oncogene expression and
we characterized DLBCLs following relapse
neuroblastoma oncogenesis, and provide
using exome and genome sequencing. We
the first evidence that amplification of long
have also implemented approaches to detect
noncoding RNA genes can contribute to
mutant cell-free DNA in the plasma known
tumorigenesis.
as circulating tumour DNA to monitor tumour
evolution. We analysed ctDNA using deep
amplicon sequencing and targeted capture
to non-invasively detect evolution of DLBCL
36
Handbook
2014 AGTA
Conference
37
Handbook
2014 AGTA
Conference
Abstract
Background: Characterization of genomic 1200 – 1230
changes driving an individual’s disease is
essential to inform rational use of targeted
INTEGRATING GENOMICS INTO CLINICAL
therapies and treatment planning for many
PRACTICE: A LOCAL AND INTERNATIONAL
patients with cancer. We have developed
PERSPECTIVE
a comprehensive, pan-solid tumor next
generation sequencing (NGS)-based Professor Kathryn North
diagnostic assay, optimized for routine FFPE MURDOCH CHILDRENS RESEARCH INSTITUTE, VIC
specimens including fine needle biopsies/
Biography
aspirates and applied it to 5,800 patients’
Kathryn is the Director of the Murdoch Childrens
tumors in a CLIA-certified, CAP-accredited lab.
Research Institute, Director of the Victoria Clinical
Methods: Hybridization capture of 3,769 Genetics Service and the David Danks Professor
exons from 236 cancer-related genes and 47 of Child Health Research at the University of
introns of 19 genes commonly rearranged was Melbourne.
applied to ≥20ng of DNA extracted from 5,800
Kathryn trained as a paediatrician, neurologist
consecutive FFPE specimens and sequenced to
and clinical geneticist in Sydney and Boston.
a unique median depth >1,000X.
Her laboratory research interests focus on gene
Results: 95% (5,508) of specimens were discovery and disease mechanism in inherited
successfully profiled, yielding 2,711 unique muscle disorders, as well as genes which
actionable alterations across 82.1% (4,520) influence normal skeletal muscle function and
of cases. 68.0% (3,744) of tumors harbored elite athletic performance. Her clinical research
≥1 actionable alteration not assayed by focuses on clinical trials for muscular dystrophy
available tumor-type specific tests or hotspot and neurofibromatosis type 1.
panels. Several novel, recurrent fusions
Kathryn is Chair of the NHMRC Research
were identified including NTRK1, FGFR2 and
Committee, and in 2012 was awarded the
BRAF fusions in pan-negative lung cancer,
Ramaciotti Medal for Excellence in Biomedical
cholangiocarcinoma, and melanoma,
Research and the Member of the Order of
respectively. Pan-negative melanoma was also
Australia (AM) for service to medicine in the
significantly enriched for NF1 loss. Constitutively
field of neuromuscular and neurogenetics
activating mutations in ESR1 were identified
research. Kathryn is one of the leaders of
as a novel resistance mechanism to hormonal
Melbourne Genomics Health Alliance, a
therapy in up to 20% of patients and ERBB2
collaborative approach to the integration of
activating mutations were enriched in 40% and
genomic information into everyday healthcare.
27.3% of micropapillary urothelial carcinoma
Melbourne Genomics is a partnership of the
and CDH1-mutated lobular breast cancer,
Royal Melbourne Hospital, Royal Children’s
respectively. Responses to targeted therapies
Hospital, Murdoch Childrens Research Institute,
against actionable alterations identified in
Walter and Eliza Hall Institute, CSIRO, Australian
patient tumors have been reported by treating
Genome Research Facility and the University
physicians.
of Melbourne. In 2013, Kathryn was invited to
Conclusions: Comprehensive genomic profiling participate in the Global Alliance for Genomics
identified actionable alterations in the majority and Health, an international consortium of more
of tumors profiled and identified potential than 130 institutions across 40 countries exploring
additional treatment options for 68% of patients the sharing of genomic and clinical data –
targeting alterations in genes not currently Kathryn is Vice Chair of the Steering Committee
assayed by available hotspot genotyping and Co-Chair of the Clinical Working Group. The
panels. Global Alliance will play a pivotal role in how
genomic data can be linked and shared on a
global level, and applied to the development of
targeted therapeutics.
38
Handbook
2014 AGTA
Conference
Abstract
disease and patient advocacy organizations.
Genomic technologies are proving
The Global Alliance aims to accelerate the
transformative, but ensuring full clinical
world-wide effort to responsibly aggregate,
and research benefit from their application
analyse and share large amounts of genomic
in clinical settings requires planning and
and clinical information to advance the
collaboration. The Melbourne Genomics
understanding, diagnosis, and treatment for
Health Alliance (MGHA) is a collaboration of
cancer, inherited diseases, infectious diseases,
10 research and healthcare organisations, with
and drug responses.
the goal of integrating genomic information
into everyday healthcare. MGHA are
conducting a clinically-led, prospective project
evaluating the feasibility of whole exome 1230 – 1300
sequencing (WES) as a single, first tier assay for
germline and somatic conditions. The Alliance
members have developed common ethics TRANSFERRING GENOMIC TESTING FROM
RESEARCH TO DIAGNOSTICS: SIMILARITIES
and consent, a common clinical bioinformatics
AND DIFFERENCES
pipeline and share genomic data via a
common clinical genomics data repository. Professor Hamish Scott
Decision-making is guided by structured CENTRE FOR CANCER BIOLOGY, SA
processes for clinician input and a Community
Biography
Advisory Group. Researchers are also an
Professor Hamish Scott did his PhD (1992) and
integral part of advisory and working groups.
first post-doc at the Women’s and Children’s
Mimicking usual clinical practice, patients Hospital and the University of Adelaide.
with one of five diverse germline or somatic During these 7 years he led the discovery of
conditions are being offered whole exome genes for 3 rare human diseases. After 11
sequencing in parallel to routine investigations. more years, with the persistence of Professor
Sequence data is generated by multiple John Hopwood and others in academia and
diagnostic laboratories, analysis is then industry, this resulted in either FDA approved
targeted to genes known to be related to the therapy (2003) or clinical trials of novel
clinical condition using a common analytic therapies for these diseases.
pipeline. Research results are returned to
In 1995, Hamish moved to the University of
clinicians and data is available to researchers.
Geneva Medical School in Switzerland. His
Exome data is linked to clinical data and a
focus was, and remains, the application
consolidated electronic view provided to
of genetic and genomic technologies to
clinicians and researchers. The pilot phase is
understand diseases processes to improve
being evaluated to determine the barriers,
diagnoses and treatment. He led international
feasibility, health economics and diagnostic
collaborations in identification of human genes
value of genomic sequencing. This systematic
causing Down syndrome and rare forms of
approach is designed to foster incremental
genetic deafness and autoimmunity (e.g.
change and future adoption, as well as
arthritis and multiple sclerosis). This continues to
ensure future implementation delivers a
have profound effects on our understanding
viable and sustainable system across multiple
of basic biology of Down syndrome, hearing
organisations.
and the immune system and lead to new
The Melbourne Genomics Health Alliance therapeutic strategies in these and related
is also a member of the Global Alliance diseases. This was also the start of his interest
for Genomics and Health (GA4GH), an in cancer and leukemia as children with
organisation of over 200 of the world’s leading Down syndrome have a low incidence of solid
biomedical research institutions, healthcare tumours and a high incidence of leukemia.
providers, information technology and life This is also when he started to work on familial
science companies, funders of research, and predisposition to leukemia.
39
Handbook
2014 AGTA
Conference
40
Handbook
2014 AGTA
Conference
pathogenic or potentially pathogenic variants into both lab bench work and informatics
by NGS technology. It will also briefly discuss pipelines. This is further confounded by the
the various applications and designs of NGS reality that many commercial components
and the best fit options for specific clinical care of such pipelines have not been designed
outcomes. strictly with diagnostics in mind nor are
they sold with that intention. This leaves the
laboratory implementing such technologies in
42
Handbook
2014 AGTA
Conference
43
Handbook
2014 AGTA
Conference
files) are often trimmed to discard poor quality high quality variants, which in some cases
bases at the ends of the reads, and then have 100% sensitivity, and as such have been
aligned to the human reference genome. clinically accredited without requiring Sanger
This produces SAM or BAM files which can be sequencing validation.
viewed in most modern genome browsers (eg In this presentation, I will run through the
IGV, UCSC genome browser). Short sequencing bioinformatics and computing aspects of a
reads are usually aligned individually, leading typical clinical genomics experiment, touching
to errors around regions of the genome that on both Mendelian and cancer genomics, and
are challenging to align to (eg known insertions some of the pitfalls.
and deletions), which are typically cleaned
up by realigning all reads around known
indels. Finally single nucleotide and small
insertion and deletion variants are identified 1640 – 1705
through a range of algorithms that identify
recurrent evidence for a variant across multiple
Curation and reporting of genetic and
overlapping sequencing reads. Some groups
genomic data in the modern era (2014!)
are performing copy number variant analysis
on NGS data, mostly based upon the idea of Professor Hamish Scott
read depth, however this is challenging due to CENTRE FOR CANCER BIOLOGY, SA
variability in sequencing depth and the patchy
Biography
sequencing coverage over only targeted NGS
Professor Hamish Scott did his PhD (1992) and
data. Variants are then extensively annotated
first post-doc at the Women’s and Children’s
and ultimately filtered to those that are rare in
Hospital and the University of Adelaide.
healthy individuals, known or predicted to be
During these 7 years he led the discovery of
damaging and associated with the patient’s
genes for 3 rare human diseases. After 11
phenotype. The final challenging task is the
more years, with the persistence of Professor
clinical interpretation, determining which of
John Hopwood and others in academia and
the one or two variants that match the patient
industry, this resulted in either FDA approved
phenotype from a list of often hundreds of
therapy (2003) or clinical trials of novel
variants. Due to the larger scope of these tests,
therapies for these diseases.
variants are usually orthogonally validated,
using Sanger sequencing, prior to return of In 1995, Hamish moved to the University of
results. Geneva Medical School in Switzerland. His focus
was, and remains, the application of genetic
In the context of cancer genomics, in the
and genomic technologies to understand
diagnostic setting, typically the DNA samples
diseases processes to improve diagnoses and
are far more challenging, usually with FFPE
treatment. He led international collaborations
material from either tumour blocks or biopsy,
in identification of human genes causing Down
resulting in low amounts of highly fragmented
syndrome and rare forms of genetic deafness
DNA. This has resulted in many groups adopting
and autoimmunity (e.g. arthritis and multiple
an amplicon based strategy around known
sclerosis). This continues to have profound
cancer hotspot mutations, micropanels of
effects on our understanding of basic biology
single genes (eg BRCA1 and BRCA2) or when
of Down syndrome, hearing and the immune
feasible, across panels of up to 300 cancer
system and lead to new therapeutic strategies
genes. Similar analysis strategies exist to
in these and related diseases. This was also
Mendelian genomics, though amplicon-based
the start of his interest in cancer and leukemia
sequencing data has nuances which must be
as children with Down syndrome have a low
accounted for in the sequence alignment and
incidence of solid tumours and a high incidence
variant calling stages. Typically the sequencing
of leukemia. This is also when he started to work
depth achieved is >500x, which results in very
on familial predisposition to leukemia.
44
Handbook
2014 AGTA
Conference
45
Handbook
2014 AGTA
Conference
46
Handbook
2014 AGTA
Conference
47
Handbook
2014 AGTA
Conference
48
Handbook
2014 AGTA
Conference
to transcript degradation via NMD. Thus, by if its predicted targets that lie in a common
coordinating inclusion of recursive exons with biological pathway have changed in the
the use of upstream promoters, RSS can control opposite direction. Utilising a database of
the stability of resulting transcripts. Collectively, predicted micro-RNA regulatory targets and
we demonstrate how exon definition and splice pathway annotations, pMimCor circumvents
site competition enable recursive splicing, much of the need to produce large datasets
a new mechanism for dictating promoter- in order to generate stable and accurate
dependent splicing and transcript stability in network estimation. Furthermore, by ranking
the brain. relevant micro-RNA and pathways concurrently
pMimCor also provides highly interpretable
biological information that can be used to
micro-RNA target and pathway databases Rima completed her PhD in Bioinformatics and
we propose an approach to stabilise the Computational Drug Design in the University
estimation and annotation of micro-RNA of Illinois at Chicago in 2010. Soon after, she
regulation that is suitable for datasets with moved to the Barcelona Supercomputing
small sample sizes. Specifically, we propose a Center and Institute for Research in
supervised framework, built upon concepts of Biomedicine as a postdoctoral candidate in
significance combination, for jointly ranking Spain to work on methodology development.
regulatory micro-RNA and their potential Next, she moved to Sydney, Australia and
functional impacts with respect to a condition joined the Garvan Institute for systems biology
of interest. and bioinformatic research in Diabetes and
Obesity.
Our proposed method, pMimCor, directly tests
if a micro-RNA is differentially expressed and
49
Handbook
2014 AGTA
Conference
Abstract
Insulin resistance (IR) is perhaps the best
predictor of future development of type 2 WEDNESDAY
diabetes (T2D). Identification of molecular 15 October 2014
signatures that can identify individuals with
higher risk of developing T2D could enable
early stage intervention. Gene expression data
could provide an ideal tool for identification
of such a molecular signature, but human SESSION 7
gene expression (GE) data is inherently noisy
and highly variable. To overcome the lack
of power, we integrated human data with
more stable data from a model organism and FUNCTIONAL GENOMICS
thereby developed a cross-species (human-
mouse) analysis platform. We applied this new CHAIRS: KAYLENE SIMPSON
approach to identify a gene signature that
AND AMEE GEORGE
could differentiate insulin resistant from insulin
sensitive individuals with improved prediction
accuracy (~20% increase) compared to 0900 – 0945
standard clinical measures (OCM). This
signature also identified beta-catenin and
JAK1 as novel ‘connection hubs’ to the Insulin EXPLORING IMMUNE SIGNALING SYSTEMS
Signaling Pathway whose GE pattern highly WITH HIGH THROUGHPUT, HIGH CONTENT
correlated to extensive metabolic phenotypic SCREENING
measures of insulin sensitivity in humans.
Dr Iain Fraser
Inhibiting these proteins impaired insulin-
NATIONAL INSTITUTES OF HEALTH, USA
stimulated glucose uptake in L6 myotubes,
confirming their role in insulin action. These Biography
data indicate the potential utility of using Iain Fraser is Chief of the Signaling Systems
systems biology approaches to segregate Unit in the Laboratory of Systems Biology at
individuals early based on differential diabetes the National Institute of Allergy and Infectious
risk. Disease, National Institutes of Health. He
received his B.S. in biochemistry from Heriot-
Watt University, Edinburgh, and his Ph.D.
in biochemistry from Imperial College,
University of London. He was a Wellcome Trust
International postdoctoral fellow at the Vollum
Institute and then co-director of the Alliance
for Cellular Signaling (AfCS) molecular biology
group at the California Institute of Technology.
He was appointed as an Investigator at the
NIH in 2008. His research has focused on the
mechanistic basis of cellular signaling, both in
G protein signaling systems and more recently
in Toll-like receptor signaling in innate immune
cells. He applies systems biology approaches
to decipher how mammalian cells integrate
stimuli in a complex environment to ensure
context-dependent cellular responses. This is
vital to understanding how a breakdown in
50
Handbook
2014 AGTA
Conference
51
Handbook
2014 AGTA
Conference
52
Handbook
2014 AGTA
Conference
1115 – 1145
Biography
Robert (Bobby) has extensive experience in
using PacBio for human and bacterial genome
sequencing. He has also implemented the
technology in a clinical setting
Abstract
Whether surveying the human genome for
pathologically relevant complex structural
variation or investigating various microbes
related to infectious disease, long read
length sequencing is rapidly enabling more
robust references that include previously
unachievable domains. Recent headway in
single molecule real time (SMRT) sequencing
has enabled faster, more comprehensive
diagnostics assays. Here, we focus specifically
on those methods related to uncovering
genetic structural variation towards a more
comprehensive sequencing pipeline. Despite
the progressive advances in next-generation
sequencing technologies over the past
decade, a variety of hard to sequence
repetitive human genetic loci remain
53
Handbook
2014 AGTA
Conference
54
Handbook
2014 AGTA
Conference
55
Handbook
2014 AGTA
Conference
56
Handbook
2014 AGTA
Conference
57
Handbook
2014 AGTA
Conference
POSTERS
58
Handbook
2014 AGTA
Conference
POSTER PRESENTATIONS
MONDAY 13 OCTOBER 2014
14 Controlling the error rate correctly when intersecting gene lists Aaron Lun
21 Why weight? Variance modelling for designed RNA-seq experiments Matthew Ritchie
High through-put Sequencing Guides the Design of DNA Minor Groove Anandhakumar
42 Binding Small Molecule with Improved Binding Specificity Chandran
POSTER PRESENTATIONS
TUESDAY 14 OCTOBER 2014
GENOME-WIDE MULTIMARKER MODEL TO PREDICT THE RISK OF BLADDER CANCER Evangelina Lopez De
57 RECURRENCE Maturana
69 A Handheld qPCR Device for Point-of-Care and In-Field Diagnostics Jo-Ann Stanton
74 Epilepsy: new genes and old pathways in a complex disorder Leanne Dibbens
76 Imperfect centered miRNA binding sites are common and functional Nicole Cloonan
The Clontech® SMARTer® UltraTM Low Input RNA Kit for Sequencing - v3
cDNA synthesis for single-cell transcriptome analysis that delivers:
Illumia is a regestered trademark of Illumia, Inc. For Research Use Only. Not for use in diagnostic or therapeutic procedures. Not for
resale. Your use of these products and technologoies is subject to compliance with any applicable licensing requirements described
on the product’s webpage at www.clontech.com. It is your responsibility to review, understand and adhere to any restrictions
imposed by such statements. Clontech, the Clontech logo, SMARTer and Ultra are trademarks of Clontech Laboratories, Inc. All
other marks are the property of their respective owners. Certain trademarks may not be registered in all jurisdictions.
Clontech is a Takara Bio Company © 2014.
66
Handbook
2014 AGTA
Conference
the kinases tested. When using the model to generation sequencing, proteomics and
complement sequence based kinase-specific ChIP-seq platforms. Currently the database
phosphorylation site prediction, we find contains mouse and human expression data,
that the additional information can greatly with support for additional model organisms
increase prediction performance at low false being planned. The user interface is designed
positive levels. Our results demonstrate that for rapid assessment of the expression of
the model can account for the short-falls in individual or multiple genes across and/or
sequence information and provide a robust within datasets. Users can store gene lists and
description of the regulation of kinase-protein share views with collaborators. Stemformatics
phosphorylation. is funded as part of the ARC special initiative
"Stem Cells Australia", has an international
community of users and is freely accessible to
Poster 3 all academic groups.
67
Handbook
2014 AGTA
Conference
68
Handbook
2014 AGTA
Conference
Poster 6 Poster 7
69
Handbook
2014 AGTA
Conference
Poster 8 Poster 9
70
Handbook
2014 AGTA
Conference
Poster 10 Poster 11
71
Handbook
2014 AGTA
Conference
Poster 12 Poster 13
72
Handbook
2014 AGTA
Conference
Poster 14 Poster 15
Mr Aaron Lun, Professor Gordon Smyth Paul Maclean, Vanessa Cave, Marlon Reis and
Christine Couldrey
Controlling the error rate correctly when
intersecting gene lists Integrating Next Generation Sequencing
Data with additional ‘omics datasets
Abstract
Intersection of lists of differentially expressed Abstract
(DE) genes is often performed to identify genes It is becoming more and more common to
that are changing across multiple comparisons. collect data from multiple ‘omics sources on
Knowing that a gene is potentially involved in an individual sample. Exploration of ‘omics
multiple biological processes can provide more data can shed light on the molecular factors
insights into its function. However, as with all that help determine the phenotype of an
statistical procedures, a degree of uncertainty individual, and a greater understanding of the
is involved in the intersection operation. This underlying biological system can be achieved
uncertainty is represented by the type I error by simultaneously looking at data from multiple
rate of the intersection, i.e., the probability of sources.
obtaining a gene in the intersected subset that This talk provides a brief overview of integrating
is actually non-DE in one or more comparisons. Next Generation Sequencing data with
A common strategy is to identify DE genes additional ‘omics datasets, and highlights some
in each comparison at a particular false methodology and considerations for a new
discovery rate (FDR), and then identify the biology.
shared subset of genes that are DE in each The talk is centred on research on livestock at
comparison. This approach will not control AgResearch Ltd.
the type I error rate, nor will it control the
corresponding FDR in the intersected subset.
Proper control of the type I error rate should
instead be maintained with an intersection-
union test (IUT). This is a very conservative test
in most cases, and will result in loss of detection
power.
Here, the inappropriateness of FDR-based
intersections and the need for the IUT are
demonstrated with some examples. A more
powerful version of the IUT is also proposed for
cases involving two DE comparisons that are
technically independent, i.e., using separate
datasets. The performance of this new method
is tested on both simulated and real data.
73
Handbook
2014 AGTA
Conference
Poster 16 Poster 17
Dr Yalchin Oytam, Dr Fariborz Sobhanmanesh, Dilmi Perera, Diego Chacon, Julie A. I. Thoms,
Dr Jason Ross, Dr Konsta Duesing Rebecca C. Poulos, Dominik Beck, Peter J.
Campbell, John E. Pimanda and Jason W. H.
Risk-conscious correction of batch
Wong,
effects: Maximising information extraction
from high-throughput genomic datasets OncoCis : An annotation tool for cis-
regulatory mutations in Cancer
Abstract
Identification and effective removal of batch Abstract
effects is a major challenge in genomic Whole genomes are being sequenced at
datasets. Batch effects are structured technical an accelerated pace but annotation of
noise which persist under best laboratory mutations and inference of their functional
practice, and typically account for more of the significance remains challenging. Whilst a
data variance than any of the experimental myriad of tools are available for the annotation
factors. Even with ideal experimental design, of protein coding mutations, few are suited
they artificially inflate within-group variances, for annotating non-coding mutations. To this
significantly reducing the power of statistical end, we have developed, OncoCis, to provide
tests, and resulting in actual effects going an easy to use web service for researchers to
undetected. Conversely, if “overcorrected” annotate and prioritize potential causal cis-
– i.e. genuine biological variance is removed regulatory mutations. In order to annotate
along with batch noise – the result is an artificial non-coding mutations, OncoCis integrates
deflation of within-group variances and as publicly available datasets from genome-
a consequence, false positives. Here, we wide DNase-seq and histone ChIP-seq data
present a novel technique which maximises the obtained from ENCODE and the Epigenome
removal of batch effects, with the constraint Atlas, to identify mutations that occur within
that the probability of overcorrection is kept potential cis-regulatory regions in a cell type
to a fraction which is set by the end-user. specific manner. These mutations are further
We benchmark the new technique against annotated with sequence conservation scores
the leading technique currently used, using and searched for possible removal or creation
data from published studies, with the two of transcription factor consensus binding
performance criteria being removal of (batch) motifs. Finally, FANTOM 5 data combined with
noise using guided-PCA as measure (Reese the GREAT tool is used to map mutations to
et. al, 2013), and preservation of (biological) the most likely gene that it may be regulating.
signal. The new method performs better If gene expression data is available, a fold
than the leading method on both fronts change is calculated for each of the mapped
simultaneously. For all datasets, while both genes. We have applied this method to
methods appear to remove batch noise whole genome sequencing data from 21
(guided-PCA p-values >> .05) our method breast cancer samples. Eighteen putative
does so with markedly lower probabilities for cis-regulatory mutations were identified
batch effect persistence, and manages to do in a number of samples that significantly
so while removing significantly less variance increased the expression of the mapped gene
from the datasets. Evidence suggests that when compared with other samples without
the novel method for batch correction, the mutation.This tool will be invaluable to
applicable to microarray and sequencing data researchers seeking to prioritise large lists
alike, is a notable advance on techniques of mutations for evaluating their functional
currently used. A software package has been relevance.transcription factor consensus
developed for public release. binding motifs. Finally, FANTOM 5 data
74
Handbook
2014 AGTA
Conference
combined with the GREAT tool is used to map We applied DiffVar to several publicly available
mutations to the most likely gene that it may be cancer datasets from The Cancer Genome
regulating. If gene expression data is available, Atlas, as well as a publicly available aging
a fold change is calculated for each of the dataset, all of which had samples profiled using
mapped genes. We have applied this method Illumina’s popular 450K human methylation
to whole genome sequencing data from 21 array. Unsurprisingly, the majority (>95%) of the
breast cancer samples. Eighteen putative significant differentially variable CpG sites in all
cis-regulatory mutations were identified cancer datasets tested showed more variability
in a number of samples that significantly in the cancer samples compared to the normal
increased the expression of the mapped gene samples. Interestingly, in the aging dataset,
when compared with other samples without 97% of the significant differentially variable
the mutation.This tool will be invaluable to CpGs were more variable in centenarians than
researchers seeking to prioritise large lists in newborns. A closer look at the genomic
of mutations for evaluating their functional composition of the significant CpGs revealed
relevance. that a greater proportion of the top ranked
differentially variable CpGs were found in
islands compared to differentially methylated
Poster 18 CpGs, a phenomenon that was consistently
observed across all datasets. DiffVar is available
in the Bioconductor R package missMethyl.
Dr Belinda Phipson, Dr Alicia Oshlack
75
Handbook
2014 AGTA
Conference
genotype pools with no prior of the individuals AGP is crucial for determining clinical strategies
constituting the pools. for intervention and prevention in at-risk
Two experiments were conducted each individuals. Analysis of a cytokine panel in a
using the same Sequenom iPLEX panel which longitudinal Finnish cohort (N = 2200) showed
generated results for 61 SNPs. The first larger an association between low-grade chronic
experiment genotyped 1041 individuals along inflammation with (future) elevated AGP. Next,
with 22 pools samples. The pool sizes varied we identified AGP-associated transcriptional
between 18 to 26 individuals all drawn from networks through weighted gene coexpression
the 1041 individuals. A second experiment network analysis on whole blood samples.
genotyped a sub set of 78 individuals from the Four transcriptional networks were found to
first experiment. We compared the performance be significantly associated with AGP, and
of the proposed calibration method with some of these sub-networks and their associations
the existing results. The RMSE error obtained using replicated in an independent cohort.
no calibration, commonly used piecewise linear Adjustment for lipid associated variation
transformation, and the proposed local–global identified one sub-network, enriched for host
learner fusion method are 0.135, 0.120, and response to bacterial and fungal infection.
0.080 respectively. Further evaluation revealed that its core genes
code for products released during neutrophil
degranulation, and Mendelian Randomisation
Poster 20 was used to infer causal or reactive
relationships between this sub-network and
AGP. Our study provides a general framework
Mr Scott Ritchie, Dr Peter Würtz , Ms Artika Nath, for integrative analysis of multi-omic data
Dr Gad Abraham, Dr Antti Kangas , Dr Aki sources, identifies molecular targets for follow-
Havulinna, Dr Pasi Soininen, Dr Kristiina Aalto, Dr up studies, and suggests prolonged innate
Johannes Kettunnen, Dr Michael Inouye immune response may be a driver of increased
mortality risk.
Integrative omics analysis elucidates
role for inflammation and neutrophil
degranulation in glycoprotein-associated
all-cause mortality risk
Abstract
Integration of molecular information captured
from high-throughput platforms has shown
immense promise for the elucidation of novel
biological processes underlying complex
disease phenotypes. Here, we integrate
genomic, transcriptomic, metabolomic,
and cytokine data in two Finnish population
cohorts to identify biological processes
underlying elevated circulating alpha-1-
acid glycoprotein (AGP) levels, a recently
discovered biomarker for 5-year all-cause
mortality risk. Circulating AGP is influenced by
a range of physiological effects; for example
AGP increases in response to inflammation,
infection, or tissue-injury, and interacts with a
variety of drugs. Identifying and understanding
the precise biological processes underlying
76
Handbook
2014 AGTA
Conference
Poster 21 Poster 22
Matthew E. Ritchie, Ruijie Liu, Gordon K. Smyth Alan F. Rubin, Lea M. Starita, Stanley Fields,
Terence P. Speed, and Douglas M. Fowler
Why weight? Variance modelling for
designed RNA-seq experiments Statistical methods and software for
functionally characterizing every single
Abstract mutation in Your Favorite Gene
Outlier samples are relatively common in
Abstract
RNA-seq experiments and the root cause
of such variation is generally unknown. In The use of high-throughput DNA sequencing
small experiments, the analyst is left with the has rapidly expanded catalogues of normal
difficult decision of what to do: removing the and disease-associated variation, but the
offending sample may reduce variation, but functional consequences of most mutations
at a cost of reducing power, which can limit are unknown. In deep mutational scanning,
our ability to detect biologically meaningful selection for protein function applied to
changes. A compromise is to use all of the a library of protein variants is combined
available data, but to down-weight the with high-throughput DNA sequencing to
observations from the outlier sample in the directly measure the activity of hundreds of
analysis. In this presentation we describe thousands of variants of the protein easily
a statistical approach that allows this by and cheaply. This approach helps to bridge
modelling heterogeneity at both the sample the gap between variant identification
and observational level in the differential and interpretation, allowing researchers to
expression analysis. Using both simulations and elucidate sequence-function relationships
real data, we tease apart scenarios where this at high resolution. The resulting data can
strategy leads to a more powerful analysis from be used in a variety of contexts, from aiding
those when it is better to discard the outlier the assessment of clinical variants to guiding
altogether. Our approach is implemented in protein engineering. Despite the growing
the open-source limma package available popularity of deep mutational scanning,
from Bioconductor (http://www.bioconductor. there are few tools and no formal statistical
org). methods available to help analyze these
complex datasets. Here we present a novel
method for assigning functional scores and
statistical significance to all variants in a
deep mutational scanning dataset based
on weighted regression. Our methods are
implemented as part of Enrich 2, a software
package that makes the initial data analysis
accessible to experimental biologists while
providing an extensible framework for
bioinformaticians manipulating these large
datasets. We share our results from applying
these methods to deep mutational scans of
BRCA1, the WW protein-binding domain of
YAP65, and other targets. Enrich 2 fully supports
barcode sequencing, complex experimental
designs involving controls and biological
replicates, and deep mutational scans of
noncoding sequences, such as structural RNAs
or regulatory elements.
77
Handbook
2014 AGTA
Conference
Poster 23 Poster 24
78
Handbook
2014 AGTA
Conference
79
Handbook
2014 AGTA
Conference
Using 12 extreme preterm birth cases and 12 A comparison of control samples for ChIP-
matched, term controls, we extracted DNA seq of histone modifications
from archived neonatal blood spots and blood
Abstract
collected in a similar way at 18 years of age.
The advent of high-throughput sequencing
DNA methylation was measured at 347,789
has allowed genome wide profiling of histone
autosomal locations throughout the genome
modifications by Chromatin ImmunoPrecipitation
using Infinium HM450 arrays. Representative
(ChIP) followed by sequencing (ChIP-seq). In
methylation differences were confirmed by
this assay the histone mark of interest is enriched
Sequenom MassArray EpiTyping.
through a chromatin pull-down assay using
At birth we found 1,555 sites with significant an antibody for the mark. Due to imperfect
differences in methylation between term and antibodies and other factors, many of the
preterm babies. At 18 years of age, these sequenced fragments do not originate from
differences had largely resolved, suggesting the histone mark of interest, and are referred
that DNA methylation differences at birth are to as background reads. Background reads
mainly driven by factors relating to gestational are not uniformly distributed and therefore
age, such as cell maturity and to some extent, control samples are usually used to estimate the
cell composition. Using matched longitudinal background distribution at any given genomic
samples, we found suggestive evidence for an position. The Encyclopedia of DNA Elements
epigenetic legacy of preterm birth, identifying (ENCODE) Consortium guidelines suggest
persistent methylation differences at several sequencing a whole cell extract (WCE, or “input”)
genomic loci. Longitudinal comparisons sample, or a mock ChIP reaction such as an IgG
uncovered a significant overlap between sites control, as a background sample. However, for
that were differentially-methylated at birth and a histone modification ChIP-seq investigation it is
those that changed with age. We found an also possible to use a Histone H3 (H3) pull-down to
unexpectedly high number of birth- and age- map the underlying distribution of histones.
80
Handbook
2014 AGTA
Conference
81
Handbook
2014 AGTA
Conference
82
Handbook
2014 AGTA
Conference
Poster 33 Poster 34
83
Handbook
2014 AGTA
Conference
The bovine gene expression atlas generated parasitic plants. In addition to advancing
by this study will advance our understanding our understanding stress in woody hosts it will
of the dynamics and complexity of bovine increase our understanding of the link between
transcription and provide a valuable resource belowground and aboveground events and
for future research on functional genomics in provide genetic information that can be used
the bovine species. for studying natural communities.
Poster 35 Poster 36
84
Handbook
2014 AGTA
Conference
transcript may interact differently with different Now, with the cost per marker continuing to
sense transcripts in the same locus. Overall, our drop, we investigate the feasibility of genome-
study improves the annotation of the soybean wide association and genomic selection for
genome by identifying numerous novel improving essential oil yield in E. polybractea.
transcripts and demonstrates the complex We compare traditional breeding and marker-
regulation of splicing variants with or without assisted selection for essential oil yield, and how
an association with antisense transcripts during the use of whole genome re-sequencing and
floral initiation. genomic selection may further improve selection
of elite individuals.
Poster 37
Poster 38
85
Handbook
2014 AGTA
Conference
Abstract
Mnasri Sameh, Saddoud Dabbebi Olfa, Ben
The rubber tree, Hevea brasiliensis, is an
Saleh Mohamed, Professeur Ferchichi Ali
important plant species that is commercially
Microsatellite marker-based identification grown to produce latex rubber in many
and genetic relationships of millennium countries. The rubber tree variety BPM 24
olive cultivars in Tunisia exhibits cytoplasmic male sterility, inherited
from the variety GT 1. We constructed the
Abstract
rubber tree mitochondrial genome of a
Microsatellite markers were used to characterize
cytoplasmic male sterile variety, BPM 24, using
the millennium olive cultivars localized in
454 sequencing, including 8 kb paired-end
nine different archeological sites in Tunisia.
libraries, plus Illumina paired-end sequencing.
Thirty genotypes were considered for genetic
We annotated this mitochondrial genome
fingerprinting using 10 pairs of microsatellite
with the aid of Illumina RNA-seq data and
primers. The number of alleles per locus ranged
performed comparative analysis. We then
from 3 to 5, with a mean of 3.7 alleles per
compared the sequence of BPM 24 to the
primer pair (a total of 37 alleles). The observed
contigs of the published rubber tree, variety
heterozygosity ranged from 0.4 to 1, while the
RRIM 600, and identified a rearrangement
expected heterozygosity varied between 0.37
that is unique to BPM 24 resulting in a novel
and 0.74. The polymorphism information content
transcript containing a portion of atp9. The
values ranged also from 0.37 to 0.74. The mean
novel transcript is consistent with changes
polymorphism information content value of 0.61
that cause cytoplasmic male sterility
for the SSR loci provided sufficient discriminating
through a slight reduction to ATP production
ability to evaluate the genetic diversity among
efficiency. The exhaustive nature of the search
the millennium cultivars. The UPGMA cluster
rules out alternative causes and supports
analyses using Jaccard’s index permitted a
previous findings of novel transcripts causing
segregation of the thirty millennium cultivars in
cytoplasmic male sterility.
three main groups and revealed that most of
the millennium cultivars grouped according
morphological parameters of the fruit and the
endocarp and no clear clustering trends were
observed according to their geographic origin.
As a sequel to the present work, new surveys
should be made in the archeological sites
localized in North and the Center of Tunisia to
sample more cultivars and to draw a clearer
picture of the diversity of the Tunisian millennium
olive germplasm.
86
Handbook
2014 AGTA
Conference
87
Handbook
2014 AGTA
Conference
88
Handbook
2014 AGTA
Conference
89
Handbook
2014 AGTA
Conference
90
Handbook
2014 AGTA
Conference
field cancerization theory. Oral epithelial biopsies tumour tissues, suggesting DPM accumulation is
can be separated based on their severity using a general feature of carcinogenesis and could
the number of mutations in their exome. be used as a novel cancer biomarker.
Poster 47 Poster 48
with accumulation of number of genetic Tumour cells undergo apoptosis release small
abnormalities. The vast majority of those amount of cell-free circulating tumour DNA
events, called passenger mutations, are non- (ctDNA) into blood throughout all stages
recurrent and dispersed across the cancer of cancer. These ctDNA carry the somatic
genome. A previous study has shown that mutations which have accumulated in the
deleterious passenger mutations (DPM) occur tumour cells, including copy number alteration
more frequently in tumour tissue compared (sCNA) of whole chromosomes (aneuploidy)
to surrounding non-tumour tissue. Therefore, and chromosome arms. Aneuploidy is also a
we hypothesise that there is a stepwise mark of high grade tumour development and
accumulation of DPM in HCC progression. To poor survival prognosis. We hypothesised that
address this, we aimed to characterise DPM in next generation sequencing of cell free DNA
whole exome sequencing (WES) datasets from in plasma may enable us to both detect the
patients in various stages liver disease. presence of ctDNA, and also characterise its
copy number aberration profile. Many existing
We analysed 3 publically-available WES
methods focus on using aberrant read-depth
datasets of HCC and paired non-tumour tissue
to infer the presence of copy number variation.
and generated our own datasets from donor
We have investigated the use of aberrant
liver tissue as a negative control. After filtering
allelic imbalance, as measured by the allele-
out probable germline mutations, we used two
specific read depths at SNP loci. We have
algorithms (PolyPhen and SIFT) to predict the
adapted methods which combine population
effect of the missense mutations on protein
level phasing with read-backed phasing in
function. We also determined whether these
order to maximise the information content in
SNPs were likely to alter genes expressed in the
allelic imbalance. We demonstrate that allelic
hepatocytes by filtering out those genes not
imbalance, modelled in this way, provides
detected in liver tissue based on our previous
sufficient power to detect tumour DNA even
data from microarray analysis.
when the mixture sample of tumour and normal
When more DPMs (but not benign SNPs) DNA at a minimum 5% tumour cellularity.
normalised to the total number of non- Moreover we are also able to pin down to
synonymous SNPs, significantly were found in the regions with copy number changes and
tumour tissue vs. non-tumour. This relationship predict possible aneuploidy.
was maintained when only liver-expressed
genes were analysed, suggesting that
damaging mutations affect cell phenotypes.
Accumulation was observed despite the
aetiology and known driver mutations in the
91
Handbook
2014 AGTA
Conference
92
Handbook
2014 AGTA
Conference
Poster 50 Poster 51
93
Handbook
2014 AGTA
Conference
94
Handbook
2014 AGTA
Conference
Our findings point to improvements in genetic and mRNA biomarkers. However, to date
screening for individuals at risk of early-onset there have been no such analyses applied
cancers and demonstrate the utility of germline to microRNA-based signatures of prognosis.
genetic variation in clinical decision-making. Our findings provide direct insights into the
translational potential of some signatures, while
highlighting critical limitations of others.
Poster 53
Poster 54
Ms Kaushala Jayawardana , Dr Sarah-Jane
Schramm , Dr Varsha Tembe, Dr Samuel Müller,
Professor Richard A. Scolyer, Professor Graham Anabel Kearney, Carrie Van Der Weyden &
J. Mann, Associate Professor Yee Hwa Yang Silvia Ling
95
Handbook
2014 AGTA
Conference
96
Handbook
2014 AGTA
Conference
97
Handbook
2014 AGTA
Conference
98
Handbook
2014 AGTA
Conference
99
Handbook
2014 AGTA
Conference
100
Handbook
2014 AGTA
Conference
101
Handbook
2014 AGTA
Conference
Poster 64 Poster 65
Dr Ruth MacKinnon, Dr Meaghan Wall, Gisela Mir, Tim Semple, Jason Li, Richard Tothill,
Associate Professor Lynda Campbell Andrew Fellowes
102
Handbook
2014 AGTA
Conference
103
Handbook
2014 AGTA
Conference
Poster 67 Poster 68
104
Handbook
2014 AGTA
Conference
Poster 69 Poster 70
Dr Jo-Ann Stanton, Dr Chris Rawle, Ms Christy Dr Siamak Tafavogh, Associate Professor Daniel
Rand, Mr Chris Mason, Dr Richard Hall, Dr Catchpoole, Associate Professor Paul Kennedy
Jeremie Langlet, Dr Joanne Hewitt, Dr Hugo
Generating Actionable Knowledge from
Strydom, Dr Angela Todd
Complex Genomic Data for Personalised
A Handheld qPCR Device for Point-of-Care Clinical Decisions
and In-Field Diagnostics
Abstract
Poster 71 Poster 72
Ms Kim Ton, Macmil, S, Miller, A.L,Chua, E.W, Dr Margaret Jordan, Dr Louise Laverick , Dr
Rand, C, Stanton, J.A, Kennedy, M.A Melissa Gresle, Dr Dragana Stanley, Ms Letitia
Smith, Dr Judith Field, Dr Laura Johnson, Dr Tim
Next generation sequencing for
Spelman, Professor Helmut Butzkueven, Professor
pharmacogenetic screening
Alan Baxter
Abstract
Causes of Multiple Sclerosis: a functional
Pharmacogenetics (PGx) is contributing to
genomics approach
the improvement of individualized medication
and the prediction of clinical treatment Abstract
outcome. With advances in this field, a Multiple Sclerosis (MS) is the most common
wide range of PGx assays have been being disabling neurological disease affecting young
developed and launched, ranging from adults in Western Society. Both environmental
low-throughput approaches that genotype and genetic associations have been confirmed.
single genetic mutations to high-throughput To date, 110 strongly associated single
approaches to evaluate hundreds of genes nucleotide polymorphisms (SNPs) have been
simultaneously. The rapid evolution of next discovered but these represent common
generation sequencing (NGS) makes it possible human genetic variations, and their individual
to sequence a large number of genes, up to association with MS is weak. The key challenge
the complete genome, in a single reaction. therefore is to identify the causal genes. One
Hi-plex is an approach for generating NGS clue lies in the fact that many SNPs are in
libraries via single multiplex PCR amplification close proximity to known immune genes that
(Nguyen-Dumont et al. 2014). This method could affect changes in immune cell function
uses automated primer design software and by changing their RNA expression and their
a simple, low-cost protocol for cost-effective encoded proteins. As it’s highly likely that eQTLs
and rapid performance. We are interested in and splice QTLs will be cell subset specific we
exploring the application of next generation wished to analyse expression levels in individual
DNA sequencing methods to understand the cell subtypes. While T cells likely play a critical
impact of pharmacogenetic variation on drug role in MS pathogenesis, successful drug trials
responses and risk of adverse drug reaction. In to date have targeted B-cells and monocytes
this study, we use this the Hi-plex approach to in addition to T cells and there is also strong
develop an throughput assay for genotyping clinical and in vitro data identifying NK cell
pharmacogenetic variants of CYP2D6, deficiencies in MS patients. We thus examined
CYP2C9, CYP2C19, DPYD, TMPT, UGT1A1, gene expression changes in 5 cell subtypes
SLCO1B1, VKORC1. These genes are involved (CD4 and CD8 T cells, B cells, monocyes and NK
in metabolism of a wide range of drugs, and cells) using Affymetrix Human Gene 1.0 ST arrays.
some are also implicated in adverse drug Expression values were standardized across
reactions. We have focused on those genes chips using RMA and quantile normalization and
with the best evidence of clinical utility (refer genes were ranked by expression difference
to PharmGKB VIP). Initially we are developing significance by MWU and ANOVA. Expression
this assay for the Ion Torrent PGM platform. differences of those genes adjacent to MS
Our current iteration of this assay contains 15 associated risk SNPs were considered. An
amplicons for these eight genes, with indexing investigation of candidate genes identified and
for 12 patients samples, and we aim to expand the pathways to which they contribute is being
this further for more routine application. undertaken to allow more targeted therapies to
Reference: Nguyen-Dumont, T., Pope, B. J., Hammet, be developed.
F., Southey, M. C., & Park, D. J. (2013). A high-plex PCR
approach for massively parallel sequencing. BioTechniques,
55(2), 69–74. doi:10.2144/000114052
106
Handbook
2014 AGTA
Conference
107
Handbook
2014 AGTA
Conference
108
Handbook
2014 AGTA
Conference
our protocol may have inflated the apparent on a collection of 86 M. tb clinical isolates from
prevalence of centered site-mediated South Africa, sequenced as part of another
interactions. study. We show that the lineages that the
Our results suggest that centered site-mediated spoligotypes associate to are consistent with
interactions are much more frequent than a phylogenetic tree derived from variation
previously thought. This may explain the across the whole genome. Furthermore,
evolutionary conservation of the central region spoligotype-defined lineage assignments
of miRNAs, and has significant implications for were consistent with those known for publicly
decoding miRNA-regulated genetic networks, available genome sequences included in the
and for predicting the functional effect of phylogenetic analysis. Our method is efficient
variants that do not alter protein sequence. and robust, meaning it is appropriate for
research deployment and public health use.
Poster 77
Poster 78
Abstract
Abstract
Spoligotyping is a genotyping assay Microsatellites are short (2-6bp) DNA
applicable for Mycobacterium tuberculosis sequences repeated in tandem, which make
(M. tb), the causative agent of tuberculosis. up approximately 3% of the human genome.
It is used for identifying the relationships These loci are prone to frequent mutations and
between isolates, tracking the spread and high polymorphism. Dozens of neurological and
evolution of M. tb. Traditional spoligotyping developmental disorders have been attributed
is a laboratory technique, whereby genomic to microsatellite expansions. Microsatellites
DNA is assayed for hybridisation to a set of have also been implicated in a range of
43 nucleotide probes, each 25 base pairs in functions such as DNA replication and repair,
length. Researchers are increasingly moving chromatin organisation and regulation of gene
away from traditional genotyping techniques expression.
in favour of high resolution next-generation Traditionally, microsatellite variation has been
sequencing (NGS), which is also generating measured using capillary gel electrophoresis.
much interest in the clinical and public health In addition to being time-consuming, and
microbiology space. However, derivation of expensive, this method fails to reveal the full
traditional genotyping information from these complexity at these loci because it cannot
data remains key for backwards compatibility detect SNP polymorphisms and compound
and contextualisation of results. microsatellites.
We present a method for deriving spoligotypes Next-generation sequencing has the potential
from NGS data. We use an alignment-free to address these problems. However,
method that avoids the cost and biases determining microsatellite lengths using next-
associated with mapping or de novo assembly, generation sequencing data is difficult. In
and results in an assay that closely mirrors the particular, polymerase slippage during PCR
original hybridization based assay. Isolates for amplification introduces stutter noise. A small
which experimental spoligotypes and NGS number of software tools claim to genotype
data are both present were not available, so simple microsatellites in next-generation
we demonstrate the efficacy of the method sequencing data, however they fail to address
109
Handbook
2014 AGTA
Conference
the issues of SNPs and compound repeats, genes associated with a disease often
and they tend to provide only approximate interact. Thus, we propose a novel kernel that
genotypes. incorporates the topology of pathways and
We have developed a microsatellite information on interactions. Using simulation
genotyping algorithm that addresses these studies, we demonstrate that the proposed
issues, providing high accuracy as well as method maintains the type I error correctly
more detailed analysis of microsatellite and can be more effective in the identification
loci. We have validated it using high depth of pathways associated with a disease than
amplicon sequencing data of microsatellites non-network-based methods. We apply our
near the AVPR1A gene. We found high approach to genome-wide association case-
concordance between our algorithm and control data on lung cancer and rheumatoid
repeat lengths obtained by electrophoresis, arthritis. We identify some promising new
manual inspection and Mendelian inheritance. pathways associated with these diseases,
By subsampling the reads, we found that our which may improve our current understanding
model is accurate to within one repeat unit of the genetic mechanisms.
down to coverages that we would expect in
standard exome sequencing.
Poster 80
Poster 79
Miriam Fanjul Fernandez, Peter Hickey, Vesna
Lukic, Natasha Brown, Greta Gillies, Sarah
Dr Saskia Freytag, Dr Juliane Manitz, Professor Wilson, Martin Delatycki, Ingrid Scheffer,
Dr Martin Schlather, Professor Dr Thomas Melanie Bahlo, Paul Lockhart
Kneib, Professor Dr Christopher Amo, Professor Gene discovery in large multiplex families
Dr Angela Risch Dr Jenny Chang-Claude, with autism spectrum disorders
Professor Dr Joachim Heinrich, Professor Dr
Heike Bickeböller Abstract
Autism spectrum disorders (ASD) are common
Identifying disease-associated pathways neuro-developmental disorders, defined
in genome-wide association studies via a by impairment in language and social
network-based kernel approach
interaction, that affect ~1:100 individuals. ASD
Abstract typically display complex inheritance with a
Biological pathways provide rich information multifactorial basis, but despite the significant
and biological context on the genetic causes heritability, the genetic basis underpinning the
of complex diseases. The logistic kernel disease remains largely unexplained. Moreover,
machine test integrates prior knowledge up to 20% of relatives of individuals with ASD
on pathways in order to analyze data display the “broader autism phenotype” (BAP),
from genome-wide association studies characterised by less severe deficits in one or
(GWAS). In this study, the kernel converts the more of the three core ASD domains.
genomic information of two individuals into Genetic studies have demonstrated that next-
a quantitative value reflecting their genetic generation sequencing (NGS) technology
similarity. With the selection of the kernel, one can be a powerful tool to identify the
implicitly chooses a genetic effect model. Like causative genes. However, the high genetic
many other pathway methods, none of the heterogeneity characterising ASD along with
standard kernels accounts for the topological studies focusing on small multiplex kindreds
structure of the pathway or gene-gene may explain previous unsuccessful molecular
interaction types. However, there is evidence studies. We have employed a novel approach
that connectivity and neighborhood of genes that analyses the BAP as marker of an ASD
are crucial in the context of GWAS, because genetic variant in large multiplex families that
110
Handbook
2014 AGTA
Conference
111
Handbook
2014 AGTA
Conference
Abstract
We set out to find the genetic cause of focal
epilepsy in a family with individuals that had
lesional and non-lesional focal epilepsy, as
112
Handbook
2014 AGTA
Conference
Poster 84 Poster 85
Lyndal Henden, Melanie Bahlo, Marie Shaw, Ms Tracey Van Stijn, Ms Shannon Clarke, Mr
Alison Gardner, Jozef Gecz Rudiger Brauning, Mr Ken Dodds, Mr John
McEwan
Inferring Identity by Descent on the X
Chromosome Genotyping-by-Sequencing (GBS):
Comparison to SNP chips and Whole
Abstract Genome Sequencing
Identity by descent (IBD) detection has played
a crucial role within the field of human disease Abstract
mapping. The idea behind IBD in disease Genotyping-by-Sequencing (GBS) has the
mapping is that if a disease is heritable, then potential to be a cost effective, reproducible
affected individuals in the same family are and high-throughput SNP genotyping method.
likely to have the same disease-causing Although we have been investigating GBS in a
variant. Affected individuals will then also number of livestock and aquaculture species,
have a shared region of IBD around the causal we have concentrated largely on sheep. To
variant that will not be present in unaffected establish the accuracy and reproducibility
individuals. Substantial progress has been of GBS, and to allow forward genotyping
made with regards to IBD methodology such capability, SNPs located in GBS targeted
as the inclusion of genotyping errors and genome regions were included on the sheep
accounting for linkage-disequilibrium, driven by High Density (HD) Illumina SNP chip. We have
the generation of denser marker data such as compared genotype calls between the two
that generated by next generation sequencing genotyping methods as well as whole genome
(NGS), however little to no work has been done sequencing (between 10x – 20x genome
on the X-chromosome. Our work focuses on the coverage).
“forgotten X” and as such we have developed
a Hidden Markov Model for inferring IBD sharing
on the X-chromosome. This method is designed
for SNP array genotype data, however it can
also be applied to NGS data once the data has
been formatted as SNP array data. Using exome
sequencing data from individuals with X-linked
intellectual disability, we show that our method
can narrow the search space for causal variants
on the X-chromosome as well as consolidate
evidence that a variant is causal.
113
Handbook
2014 AGTA
Conference
Poster 86 Poster 87
Ms Anita Goldinger , Dr Jian Yang, Dr Anjali Ms Beth Signal, Dr Brian Gloss, Mr Dominik
Henders, Dr Allan McRae, Professor Grant Kaczorowski, Mr Seth Cheetham, Ms Franziska
Montgomery, Professor Peter Visscher, Professor Gruhl, Dr Marcel Dinger
Greg Gibson, Dr Joseph Powell
The discovery of novel lincRNA with
Extensive genetic co-regulation drives potential roles in mESC development
highly correlated modular structures in
transcriptomic data Abstract
A large proportion of the genome can be
Abstract transcribed as long non-coding RNA (lncRNA),
Investigating the genetic co-regulation of RNA a class of non-coding RNA with demonstrated
transcripts can elucidate the complex genetic biological importance. Despite multiple
architecture governing the translation of gene available computational tools being available,
products into higher-order phenotypes. In a the categorisation of these transcripts as
previous study, we identified nine modules of non-coding from sequence information alone
correlated genes, which were representative remains difficult.
of distinct changes in immunological function
due to infection or geographical location. Here, It is also highly likely that novel lncRNAs
we report results from our examination into the remain to be found due to their characteristic
genetic co-regulation of these modules. The expression pattern of low levels combined with
phenotypic correlations between module probes high cell-type specificity and rapid turnover.
(ranging from 0.35 to 0.89) are significantly We hypothesised that employing higher
higher than estimates obtained from randomly temporal resolution RNA-Seq time courses
sampled probes. The probe heritability (h2), would allow us to capture novel, transiently
estimated using GREML analysis, was high with expressed or highly specific, lncRNA transcripts.
> 90% of the probes demonstrating a significant
h2 between 0.19 and 0.65. To calculate the We thus sought to identify potentially novel
amount of covariance attributed to common intergenic lncRNAs (lincRNAs) important
genomic control, we calculated the genetic in cellular differentiation using RNA-Seq
correlation (rG) using bi-variate GREML, for all data from a six-hourly time course over six
pair-wise combinations of module probes. The days of mouse embryonic stem cell (mESC)
average rG estimate of these probes was higher development. Genome guided transcriptome
than the phenotypic correlations, ranging from reconstruction was employed to generate a
0.56 to 0.96. The high rG suggests that common ‘novel transcriptome’ of mESC differentiation.
loci regulating module probes are driving the Transcripts from this were then subjected to
high phenotypic correlations. A significant a lincRNA detection pipeline established
number of trans-eQTLs were shared between by comparing the performance of multiple
probes in each module. These eQTLs were computational tools for classifying transcripts
located in genes involved in hematopoiesis as coding or non-coding on annotated
and were also enriched for previous GWAS hits mouse transcripts. This method facilitated
with various blood cell related disease. Our the discovery of hundreds of novel lincRNA
results demonstrate that these modules are with potential roles in mESC development.
under tight genetic regulation, which facilitates Insights into potential functions of these
their coordinated response to environmental novel transcripts were gained by comparing
pressures. This information has important expression profiles to those of known transcripts
application in the functional annotation of and detecting expression changes in other
genomic loci, building causal networks driving developmental transcriptome datasets.
disease and exploring the basis of disease
susceptibility between individuals.
114
Handbook
2014 AGTA
Conference
115
Handbook
2014 AGTA
Conference
genes (enriched for cell cycle, focal adhesion gene containing that exon under a certain
and transcription factors) responding to comparison. Similarly to the gene level analysis
expression changes of novel noncoding RNAs. of RNA-Seq data, we use negative binomial
Finally, this data is enabling unprecedented distribution to model exon counts. An empirical
capacity to unravel driver/carrier Bayes information sharing strategy is applied
transcriptional dynamics at complex gene loci to estimate the variation between biological
(bi-directional/antisense). replicates. Statistical tests are performed at
Fine temporal resolution has highlighted both the exon level and the gene level. The
important gene expression changes in early method has been implemented in the edgeR
differentiation and identified many noncoding package.
transcripts with a potential driver role in cellular
development. It is becoming possible to
elucidate the complexities of transcriptional Poster 91
control during cellular differentiation.
116
Handbook
2014 AGTA
Conference
117
Handbook
2014 AGTA
Conference
Poster 94
Poster 93
Dr Sonika Tyagi
Amanda Chamberlain, Christy Vander Jagt,
Michael Goddard, Leah Marrett, Ben Hayes Stranded versus non-stranded RNA-seq
118
Handbook
2014 AGTA
Conference
TO OUR GENEROUS
the opportunity to thank the following
organisations for their valuable contribution,
and ask that all delegates please consider
SPONSORS these organisations for future business
relationships.
Gold Sponsors
These revolutionary tools for DNA, RNA and SMRT technology enables real-time analysis of
protein analysis are enabling rapid advances in biomolecules with single molecule resolution,
disease research, drug development and the which has the potential to transform the
development of molecular tests in the clinic. understanding of biological systems by
providing a window into these systems that has
www.illumina.com
not previously been open for scientific study.
www.pacificbiosciences.com
119
Handbook
2014 AGTA
Conference
www.qfab.org
www.combine.org.au
120
Handbook
2014 AGTA
Conference
Illumina is advancing human health by unlocking Thermo Fisher Scientific is the world leader in
the power of the genome; from genome-wide serving science, enabling our customers to
discovery through to targeted validation and make the world healthier, cleaner and safer.
beyond. We are the leading provider of analytical
At Illumina, our goal is to apply innovative instruments, equipment, reagents and
technologies and revolutionary assays to the consumables, software and services for
analysis of genetic variation and function, research, analysis, discovery and diagnostics.
making studies possible that were not even In Australia and New Zealand, you can
imaginable just a few years ago. access the renowned Thermo Scientific
portfolio as well as premium, through Thermo
These revolutionary tools for DNA, RNA and
Fisher Scientific. Our comprehensive range
protein analysis are enabling rapid advances in
includes high-end analytical and process
disease research, drug development and the
instrumentation, laboratory equipment, a
development of molecular tests in the clinic.
complete range of consumables, chemicals
www.illumina.com and reagents as well as the service and
support to optimise your business.
www.thermofisher.com.au
Trade Booth 2
Trade Booth 4
122
Handbook
2014 AGTA
Conference
Trade Booth 12
To be your premier scientific partner, BGI Tech
is committed to providing high quality and
innovative products and solutions in biological,
agricultural and medical research.
Genesearch is the Australian distributor of
BGI Tech was founded in 2012 and is New England Biolabs products, including
headquartered in Shenzhen, China. Our service the NEBNext product range for NGS library
network covers more than 100 countries and preparation. NEBNext reagents provide
regions. Currently, BGI Tech collaborates with solutions for the major NGS platforms. Featuring
over 10,000 organizations and 30,000 partners, NEB’s NGS-optimized polymerase, Q5, for robust
providing highly reliable services for thousands performance with a wide range of amplicons,
of research projects, including several large- regardless of GC content. Please visit the NEB
scale international genome projects. To website to see the latest products including 96
date, BGI Tech has contributed to over 600 Dual Index Primers compatible with the Illumina
publications, with more than 100 articles platform.
appearing in top-tier journals such as Nature
www.genesearch.com.au
and Science.
www.bgitechsolutions.com
123
Handbook
2014 AGTA
Conference
Trade Booth 16
124
Handbook
2014 AGTA
Conference
On Q
software
OnQ Software is dedicated in providing
leading edge custom technology solutions with
the aim of enabling end to end integration
across your entire workplace using our suite
of state of the art software, tools and highly
trained experts who can tailor a solution for
your business needs.
www.onqsoft.com.au
125
Handbook
2014 AGTA
Conference
126
Handbook
2014 AGTA
Conference
127
Handbook
2014 AGTA
Conference
128
Handbook
2014 AGTA
Conference
EXHIBITION
FLOOR PLAN
EXHIBITOR BOOTH
Illumina 1
AGRF 2
Thermo Fisher 3
Scientific
Life Technologies 4
Tecan Australia 5
Roche 6
Agilent Technologies 7
Millennium Science 8&9
BGI Tech 10
Ramaciotti Centre 11
for Genomics
Genesearch 12
Bio-Strategy 13
AGTA Qiagen 14
REGISTRATION
Sarstedt 15
TrendBio 16
MIMR - PHI Institute 17
17 18 (Monash Health)
OnQ Software 18
16 19
Intergrated DNA 19
Scientifix LIFE 20
15 20 Eppendorf 21
14 21 SI Systems 22
PerkinElmer 23
DNA Genotek 24
13
Bio-Rad 25
12 22 Macrogen 26
GeneWorks 27
25 24 23
26
27 1
11
AGTA
2 SPEAKER
10
AGTA PREPARATION
PLENARY 9 3
8 4
AGTA
7 6 5 BREAKOUT
129
Handbook
2014 AGTA
Conference
AGTA 2014
DELEGATE LIST
130
Handbook
2014 AGTA
Conference
131
Handbook
2014 AGTA
Conference
132
Handbook
2014 AGTA
Conference
133
Handbook
2014 AGTA
Conference
134
Handbook
2014 AGTA
Conference
135
Handbook
2014 AGTA
Conference
136
Handbook
2014 AGTA
Conference
137
Handbook
The 2015 AGTA conference
will cover themes that include cancer genomics,
epigenomics, bioinformatics, transcriptomics, functional
genomics, clinical sequencing and plant genomics with
Save
a special focus on the genomics of good wine. We have
invited an outstanding list of international and national
speakers plus have engaged industry to showcase the best
in new genomic technologies. Please come and join us for
the Date!
some great networking and cutting edge science in the
heart of Australia’s top food and wine country.
We look forward to seeing you at AGTA 2015 in the Hunter
Valley to share the excitement of cutting-edge science!
11-14 October 2015
Co-Convenors
Marcel Dinger
Garvan Institute of Medical Research, Sydney
carsten Kulheim
Australian National University, Canberra
Any Questions?
Contact Brigitte at Leishman Associates
brigitte@leishman-associates.com.au
Phone: +61 3 6234 7844
AGTA 2014
SPONSORS AND EXHIBITORS
On Q
software