Germline ERG haploinsufficiency defines a new syndrome with cytopenia and hematological malignancy predisposition.

The genomics era has facilitated discovery of new genes predisposing to bone marrow failure (BMF) and hematological malignancy (HM). We report the discovery of ERG as a novel autosomal dominant BMF/HM predisposition gene. ERG is a highly constrained transcription factor critical for definitive hematopoiesis, stem cell function and platelet maintenance. ERG colocalizes with other transcription factors including RUNX1 and GATA2 on promoters/enhancers of genes orchestrating hematopoiesis. We identified a rare heterozygous ERG missense variant in 3 thrombocytopenic individuals from one family and 14 additional ERG variants in unrelated individuals with BMF/HM including 2 de novo cases and 3 truncating variants. Phenotypes associated with pathogenic germline ERG variants included cytopenias (thrombocytopenia, neutropenia, pancytopenia) and HMs (acute myeloid leukemia, myelodysplastic syndrome, acute lymphoblastic leukemia) with onset before 40 years. Twenty ERG variants (19 missense, 1 truncating) including 3 missense population variants were functionally characterized. Thirteen potentially pathogenic ETS domain missense variants displayed loss-of-function characteristics disrupting transcriptional transactivation, DNA-binding and/or nuclear localization. Selected variants overexpressed in mouse fetal liver cells failed to drive myeloid differentiation and cytokine-independent growth in culture, and to promote acute erythroleukemia when transplanted into mice, concordant with these variants being loss-of-function. Four individuals displayed somatic genetic rescue by copy neutral loss of heterozygosity. Identification of predisposing germline ERG variants has clinical implications for patient/family diagnosis, counselling, surveillance, and treatment strategies including selection of bone marrow donors or cell/gene therapy.

Shariant platform: Enabling evidence sharing across Australian clinical genetic-testing laboratories to support variant interpretation.

Sharing genomic variant interpretations across laboratories promotes consistency in variant assertions. A landscape analysis of Australian clinical genetic-testing laboratories in 2017 identified that, despite the national-accreditation-body recommendations encouraging laboratories to submit genotypic data to clinical databases, fewer than 300 variants had been shared to the ClinVar public database. Consultations with Australian laboratories identified resource constraints limiting routine application of manual processes, consent issues, and differences in interpretation systems as barriers to sharing. This information was used to define key needs and solutions required to enable national sharing of variant interpretations. The Shariant platform, using both the GRCh37 and GRCh38 genome builds, was developed to enable ongoing sharing of variant interpretations and associated evidence between Australian clinical genetic-testing laboratories. Where possible, two-way automated sharing was implemented so that disruption to laboratory workflows would be minimized. Terms of use were developed through consultation and currently restrict access to Australian clinical genetic-testing laboratories. Shariant was designed to store and compare structured evidence, to promote and record resolution of inter-laboratory classification discrepancies, and to streamline the submission of variant assertions to ClinVar. As of December 2021, more than 14,000 largely prospectively curated variant records from 11 participating laboratories have been shared. Discrepant classifications have been identified for 11% (28/260) of variants submitted by more than one laboratory. We have demonstrated that co-design with clinical laboratories is vital to developing and implementing a national variant-interpretation sharing effort. This approach has improved inter-laboratory concordance and enabled opportunities to standardize interpretation practices.

Extending the new era of genomic testing into pregnancy management: A proposed model for Australian prenatal services.

BACKGROUND: Trio exome sequencing can be used to investigate congenital abnormalities identified on pregnancy ultrasound, but its use in an Australian context has not been assessed. AIMS: Assess clinical outcomes and changes in management after expedited genomic testing in the prenatal period to guide the development of a model for widespread implementation. MATERIALS AND METHODS: Forty-three prospective referrals for whole exome sequencing, including 40 trios (parents and pregnancy), two singletons and one duo were assessed in a tertiary hospital setting with access to a state-wide pathology laboratory. Diagnostic yield, turn-around time (TAT), gestational age at reporting, pregnancy outcome, change in management and future pregnancy status were assessed for each family. RESULTS: A clinically significant genomic diagnosis was made in 15/43 pregnancies (35%), with an average TAT of 12 days. Gestational age at time of report ranged from 16 + 5 to 31 + 6 weeks (median 21 + 3 weeks). Molecular diagnoses included neuromuscular and skeletal disorders, RASopathies and a range of other rare Mendelian disorders. The majority of families actively used the results in pregnancy decision making as well as in management of future pregnancies. CONCLUSIONS: Rapid second trimester prenatal genomic testing can be successfully delivered to investigate structural abnormalities in pregnancy, providing crucial guidance for current and future pregnancy management. The time-sensitive nature of this testing requires close laboratory and clinical collaboration to ensure appropriate referral and result communication. We found the establishment of a prenatal coordinator role and dedicated reporting team to be important facilitators. We propose this as a model for genomic testing in other prenatal services.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Truncating Variants in NAA15 Are Associated with Variable Levels of Intellectual Disability, Autism Spectrum Disorder, and Congenital Anomalies.

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.

Whole-genome characterization of chemoresistant ovarian cancer.

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

Whole genomes redefine the mutational landscape of pancreatic cancer.

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

Hypermutation In Pancreatic Cancer.

Pancreatic cancer is molecularly diverse, with few effective therapies. Increased mutation burden and defective DNA repair are associated with response to immune checkpoint inhibitors in several other cancer types. We interrogated 385 pancreatic cancer genomes to define hypermutation and its causes. Mutational signatures inferring defects in DNA repair were enriched in those with the highest mutation burdens. Mismatch repair deficiency was identified in 1% of tumors harboring different mechanisms of somatic inactivation of MLH1 and MSH2. Defining mutation load in individual pancreatic cancers and the optimal assay for patient selection may inform clinical trial design for immunotherapy in pancreatic cancer.

Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Computational approaches to identify functional genetic variants in cancer genomes.

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

Integrated genomic and transcriptomic analysis of human brain metastases identifies alterations of potential clinical significance.

Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2) = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.

Case report of whole genome sequencing in the XY female: identification of a novel SRY mutation and revision of a misdiagnosis of androgen insensitivity syndrome.

BACKGROUND: The 46,XY female is characterised by a male karyotype and female phenotype arising due to any interruption in the sexual development pathways in utero. The cause is usually genetic and various genes are implicated. CASE PRESENTATION: Herein we describe a 46,XY woman who was first diagnosed with androgen insensitivity syndrome (testicular feminisation) at 18 years; however, this was later questioned due to the presence of intact Müllerian structures. The clinical phenotype suggested several susceptibility genes including SRY, DHH, NR5A1, NR0B1, AR, AMH, and AMHR2. To study candidate genes simultaneously, we performed whole genome sequencing. This revealed a novel and likely pathogenic missense variant (p.Arg130Pro, c.389G>C) in SRY, one of the major genes implicated in complete gonadal dysgenesis, hence securing this condition over androgen insensitivity syndrome as the cause of the patient's disorder of sexual development. CONCLUSION: This case highlights the emerging clinical utility of whole genome sequencing as a tool in differentiating disorders of sexual development.

Presentation of m.3243A>G (MT-TL1; tRNALeu) variant with focal neurology in infancy.

The Mitochondrial tRNALeu (MT-TL1) mutation, m.3243A>G constitutes the commonest identified mitochondrial genome mutation. Characteristically, giving rise to MELAS (mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes), a phenotypic spectrum associated with this genetic variant is now apparent. We report on the first patient with infantile hemiparesis, without comorbid encephalopathy, attributed to this variant. This further expands the recognized disease spectrum and highlights the need to consider mitochondrial genomic mutations in cases of cryptogenic focal neurological deficit in infancy. The potential for genetic disease modifiers is additionally discussed.

Sleeping Beauty mutagenesis reveals cooperating mutations and pathways in pancreatic adenocarcinoma.

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

Integrating massively parallel sequencing into diagnostic workflows and managing the annotation and clinical interpretation challenge.

Massively parallel sequencing has become a powerful tool for the clinical management of patients with applications in diagnosis, guidance of treatment, prediction of drug response, and carrier screening. A considerable challenge for the clinical implementation of these technologies is the management of the vast amount of sequence data generated, in particular the annotation and clinical interpretation of genomic variants. Here, we describe annotation steps that can be automated and common strategies employed for variant prioritization. The definition of best practice standards for variant annotation and prioritization is still ongoing; at present, there is limited consensus regarding an optimal clinical sequencing pipeline. We provide considerations to help define these. For the first time, clinical genetics and genomics is not limited by our ability to sequence, but our ability to clinically interpret and use genomic information in health management. We argue that the development of standardized variant annotation and interpretation approaches and software tools implementing these warrants further support. As we gain a better understanding of the significance of genomic variation through research, patients will be able to benefit from the full scope that these technologies offer.

Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling.

The importance of epigenetic modifications such as DNA methylation in tumorigenesis is increasingly being appreciated. To define the genome-wide pattern of DNA methylation in pancreatic ductal adenocarcinomas (PDAC), we captured the methylation profiles of 167 untreated resected PDACs and compared them to a panel of 29 adjacent nontransformed pancreata using high-density arrays. A total of 11,634 CpG sites associated with 3,522 genes were significantly differentially methylated (DM) in PDAC and were capable of segregating PDAC from non-malignant pancreas, regardless of tumor cellularity. As expected, PDAC hypermethylation was most prevalent in the 5' region of genes (including the proximal promoter, 5'UTR and CpG islands). Approximately 33% DM genes showed significant inverse correlation with mRNA expression levels. Pathway analysis revealed an enrichment of aberrantly methylated genes involved in key molecular mechanisms important to PDAC: TGF-ß, WNT, integrin signaling, cell adhesion, stellate cell activation and axon guidance. Given the recent discovery that SLIT-ROBO mutations play a clinically important role in PDAC, the role of epigenetic perturbation of axon guidance was pursued in more detail. Bisulfite amplicon deep sequencing and qRT-PCR expression analyses confirmed recurrent perturbation of axon guidance pathway genes SLIT2, SLIT3, ROBO1, ROBO3, ITGA2 and MET and suggests epigenetic suppression of SLIT-ROBO signaling and up-regulation of MET and ITGA2 expression. Hypomethylation of MET and ITGA2 correlated with high gene expression, which was associated with poor survival. These data suggest that aberrant methylation plays an important role in pancreatic carcinogenesis affecting core signaling pathways with potential implications for the disease pathophysiology and therapy.

Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells.

Recent RNA-sequencing studies have shown remarkable complexity in the mammalian transcriptome. The ultimate impact of this complexity on the predicted proteomic output is less well defined. We have undertaken strand-specific RNA sequencing of multiple cellular RNA fractions (>20 Gb) to uncover the transcriptional complexity of human embryonic stem cells (hESCs). We have shown that human embryonic stem (ES) cells display a high degree of transcriptional diversity, with more than half of active genes generating RNAs that differ from conventional gene models. We found evidence that more than 1000 genes express long 5' and/or extended 3'UTRs, which was confirmed by "virtual Northern" analysis. Exhaustive sequencing of the membrane-polysome and cytosolic/untranslated fractions of hESCs was used to identify RNAs encoding peptides destined for secretion and the extracellular space and to demonstrate preferential selection of transcription complexity for translation in vitro. The impact of this newly defined complexity on known gene-centric network models such as the Plurinet and the cell surface signaling machinery in human ES cells revealed a significant expansion of known transcript isoforms at play, many predicting possible alternative functions based on sequence alterations within key functional domains.

PINA v2.0: mining interactome modules.

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

Reanalysis of genomic data in rare disease: current practice and attitudes among Australian clinical and laboratory genetics services.

Reanalyzing stored genomic data over time is highly effective in increasing diagnostic yield in rare disease. Automation holds the promise of delivering the benefits of reanalysis at scale. Our study aimed to understand current reanalysis practices among Australian clinical and laboratory genetics services and explore attitudes towards large-scale automated re-analysis. We collected audit data regarding testing and reanalysis volumes, policies and procedures from all Australian diagnostic laboratories providing rare disease genomic testing. A genetic health professionals' survey explored current practices, barriers to reanalysis, preferences and attitudes towards automation. Between 2018 and 2021, Australian diagnostic laboratories performed over 25,000 new genomic tests and 950 reanalyses, predominantly in response to clinician requests. Laboratory and clinical genetic health professionals (N = 134) identified workforce capacity as the principal barrier to reanalysis. No specific laboratory or clinical guidelines for genomic data reanalysis or policies were identified nationally. Perceptions of acceptability and feasibility of automating reanalysis were positive, with professionals emphasizing clinical and workflow benefits. In conclusion, there is a large and rapidly growing unmet need for reanalysis of existing genomic data. Beyond developing scalable automated reanalysis pipelines, leadership and policy are needed to successfully transform service delivery models and maximize clinical benefit.