Seave: a comprehensive web platform for storing and interrogating human genomic variation

Filaggrin-Associated Atopic Skin, Eye, Airways, and Gut Disease, Modifying the Presentation of X-Linked Reticular Pigmentary Disorder (XLPDR)

BACKGROUND

X-linked reticular pigmentary disorder (XLPDR) is a rare condition characterized by skin hyperpigmentation, ectodermal features, multiorgan inflammation, and recurrent infections. All probands identified to date share the same intronic hemizygous POLA1 hypomorphic variant (NM_001330360.2(POLA1):c.1393-354A > G) on the X chromosome. Previous studies have supported excessive type 1 interferon (IFN) inflammation and natural killer (NK) cell dysfunction in disease pathogenesis. Common null polymorphisms in filaggrin (FLG) gene underlie ichthyosis vulgaris and atopic predisposition.

CASE

A 9-year-old boy born to non-consanguineous parents developed eczema with reticular skin hyperpigmentation in early infancy. He suffered recurrent chest infections with chronic cough, clubbing, and asthma, moderate allergic rhinoconjunctivitis with keratitis, multiple food allergies, and vomiting with growth failure. Imaging demonstrated bronchiectasis, while gastroscopy identified chronic eosinophilic gastroduodenitis. Interestingly, growth failure and bronchiectasis improved over time without specific treatment.

METHODS

Whole-genome sequencing (WGS) using Illumina short-read sequencing was followed by both manual and orthogonal automated bioinformatic analyses for single-nucleotide variants, small insertions/deletions (indels), and larger copy number variations. NK cell cytotoxic function was assessed using 51Cr release and degranulation assays. The presence of an interferon signature was investigated using a panel of six interferon-stimulated genes (ISGs) by QPCR.

RESULTS

WGS identified a de novo hemizygous intronic variant in POLA1 (NM_001330360.2(POLA1):c.1393-354A > G) giving a diagnosis of XLPDR, as well as a heterozygous nonsense FLG variant (NM_002016.2(FLG):c.441del, NP_0020.1:p.(Arg151Glyfs*43)). Compared to healthy controls, the IFN signature was elevated although the degree moderated over time with the improvement in his chest disease. NK cell functional studies showed normal cytotoxicity and degranulation.

CONCLUSION

This patient had multiple atopic manifestations affecting eye, skin, chest, and gut, complicating the presentation of XLPDR. This highlights that common FLG polymorphisms should always be considered when assessing genotype-phenotype correlations of other genetic variation in patients with atopic symptoms. Additionally, while the patient exhibited an enhanced IFN signature, he does not have an NK cell defect, suggesting this may not be a constant feature of XLPDR.

Atypical splicing variants in PKD1 explain most undiagnosed typical familial ADPKD

Autosomal dominant polycystic kidney disease (ADPKD) is the most common monogenic cause of kidney failure and is primarily associated with PKD1 or PKD2. Approximately 10% of patients remain undiagnosed after standard genetic testing. We aimed to utilise short and long-read genome sequencing and RNA studies to investigate undiagnosed families. Patients with typical ADPKD phenotype and undiagnosed after genetic diagnostics were recruited. Probands underwent short-read genome sequencing, PKD1 and PKD2 coding and non-coding analyses and then genome-wide analysis. Targeted RNA studies investigated variants suspected to impact splicing. Those undiagnosed then underwent Oxford Nanopore Technologies long-read genome sequencing. From over 172 probands, 9 met inclusion criteria and consented. A genetic diagnosis was made in 8 of 9 (89%) families undiagnosed on prior genetic testing. Six had variants impacting splicing, five in non-coding regions of PKD1. Short-read genome sequencing identified novel branchpoint, AG-exclusion zone and missense variants generating cryptic splice sites and a deletion causing critical intron shortening. Long-read sequencing confirmed the diagnosis in one family. Most undiagnosed families with typical ADPKD have splice-impacting variants in PKD1. We describe a pragmatic method for diagnostic laboratories to assess PKD1 and PKD2 non-coding regions and validate suspected splicing variants through targeted RNA studies.

TEFM variants impair mitochondrial transcription causing childhood-onset neurological disease

Mutations in the mitochondrial or nuclear genomes are associated with a diverse group of human disorders characterized by impaired mitochondrial respiration. Within this group, an increasing number of mutations have been identified in nuclear genes involved in mitochondrial RNA biology. The TEFM gene encodes the mitochondrial transcription elongation factor responsible for enhancing the processivity of mitochondrial RNA polymerase, POLRMT. We report for the first time that TEFM variants are associated with mitochondrial respiratory chain deficiency and a wide range of clinical presentations including mitochondrial myopathy with a treatable neuromuscular transmission defect. Mechanistically, we show muscle and primary fibroblasts from the affected individuals have reduced levels of promoter distal mitochondrial RNA transcripts. Finally, tefm knockdown in zebrafish embryos resulted in neuromuscular junction abnormalities and abnormal mitochondrial function, strengthening the genotype-phenotype correlation. Our study highlights that TEFM regulates mitochondrial transcription elongation and its defect results in variable, tissue-specific neurological and neuromuscular symptoms.

Biallelic pathogenic variants in COX11 are associated with an infantile-onset mitochondrial encephalopathy

Primary mitochondrial diseases are a group of genetically and clinically heterogeneous disorders resulting from oxidative phosphorylation (OXPHOS) defects. COX11 encodes a copper chaperone that participates in the assembly of complex IV and has not been previously linked to human disease. In a previous study, we identified that COX11 knockdown decreased cellular adenosine triphosphate (ATP) derived from respiration, and that ATP levels could be restored with coenzyme Q10 (CoQ10 ) supplementation. This finding is surprising since COX11 has no known role in CoQ10 biosynthesis. Here, we report a novel gene-disease association by identifying biallelic pathogenic variants in COX11 associated with infantile-onset mitochondrial encephalopathies in two unrelated families using trio genome and exome sequencing. Functional studies showed that mutant COX11 fibroblasts had decreased ATP levels which could be rescued by CoQ10 . These results not only suggest that COX11 variants cause defects in energy production but reveal a potential metabolic therapeutic strategy for patients with COX11 variants.

Whole exome and genome sequencing in mendelian disorders: a diagnostic and health economic analysis

Whole genome sequencing (WGS) improves Mendelian disorder diagnosis over whole exome sequencing (WES); however, additional diagnostic yields and costs remain undefined. We investigated differences between diagnostic and cost outcomes of WGS and WES in a cohort with suspected Mendelian disorders. WGS was performed in 38 WES-negative families derived from a 64 family Mendelian cohort that previously underwent WES. For new WGS diagnoses, contemporary WES reanalysis determined whether variants were diagnosable by original WES or unique to WGS. Diagnostic rates were estimated for WES and WGS to simulate outcomes if both had been applied to the 64 families. Diagnostic costs were calculated for various genomic testing scenarios. WGS diagnosed 34% (13/38) of WES-negative families. However, contemporary WES reanalysis on average 2 years later would have diagnosed 18% (7/38 families) resulting in a WGS-specific diagnostic yield of 19% (6/31 remaining families). In WES-negative families, the incremental cost per additional diagnosis using WGS following WES reanalysis was AU$36,710 (£19,407;US$23,727) and WGS alone was AU$41,916 (£22,159;US$27,093) compared to WES-reanalysis. When we simulated the use of WGS alone as an initial genomic test, the incremental cost for each additional diagnosis was AU$29,708 (£15,705;US$19,201) whereas contemporary WES followed by WGS was AU$36,710 (£19,407;US$23,727) compared to contemporary WES. Our findings confirm that WGS is the optimal genomic test choice for maximal diagnosis in Mendelian disorders. However, accepting a small reduction in diagnostic yield, WES with subsequent reanalysis confers the lowest costs. Whether WES or WGS is utilised will depend on clinical scenario and local resourcing and availability.

Use of Whole-Genome Sequencing for Mitochondrial Disease Diagnosis

BACKGROUND AND OBJECTIVES

Mitochondrial diseases (MDs) are the commonest group of heritable metabolic disorders. Phenotypic diversity can make molecular diagnosis challenging, and causative genetic variants may reside in either mitochondrial or nuclear DNA. A single comprehensive genetic diagnostic test would be highly useful and transform the field. We applied whole-genome sequencing (WGS) to evaluate the variant detection rate and diagnostic capacity of this technology with a view to simplifying and improving the MD diagnostic pathway.

METHODS

Adult patients presenting to a specialist MD clinic in Sydney, Australia, were recruited to the study if they satisfied clinical MD (Nijmegen) criteria. WGS was performed on blood DNA, followed by clinical genetic analysis for known pathogenic MD-associated variants and MD mimics.

RESULTS

Of the 242 consecutive patients recruited, 62 participants had "definite," 108 had "probable," and 72 had "possible" MD classification by the Nijmegen criteria. Disease-causing variants were identified for 130 participants, regardless of the location of the causative genetic variants, giving an overall diagnostic rate of 53.7% (130 of 242). Identification of causative genetic variants informed precise treatment, restored reproductive confidence, and optimized clinical management of MD.

DISCUSSION

Comprehensive bigenomic sequencing accurately detects causative genetic variants in affected MD patients, simplifying diagnosis, enabling early treatment, and informing the risk of genetic transmission.

Distinct diagnostic trajectories in NBAS-associated acute liver failure highlights the need for timely functional studies

Variants of uncertain significance (VUS) are commonly found following genomic sequencing, particularly in ethnically diverse populations that are underrepresented in large population databases. Functional characterization of VUS may assist in variant reclassification, however these studies are not readily available and often rely on research funding and good will. We present four individuals from three families at different stages of their diagnostic trajectory with recurrent acute liver failure (RALF) and biallelic NBAS variants, confirmed by either trio analysis or cDNA studies. Functional characterization was undertaken, measuring NBAS and p31 levels by Western blotting, demonstrating reduced NBAS levels in two of three families, and reduced p31 levels in all three families. These results provided functional characterization of the molecular impact of a missense VUS, allowing reclassification of the variant and molecular confirmation of NBAS-associated RALF. Importantly, p31 was decreased in all individuals, including an individual with two missense variants where NBAS protein levels were preserved. These results highlight the importance of access to timely functional studies after identification of putative variants, and the importance of considering a range of assays to validate variants whose pathogenicity is uncertain. We suggest that funding models for genomic sequencing should consider incorporating capabilities for adjunct RNA, protein, biochemical, and other specialized tests to increase the diagnostic yield which will lead to improved medical care, increased equity, and access to molecular diagnoses for all patients.

Decompensation of cardiorespiratory function and emergence of anemia during pregnancy in a case of mitochondrial myopathy, lactic acidosis, and sideroblastic anemia 2 with compound heterozygous YARS2 pathogenic variants

Myopathy, lactic acidosis, and sideroblastic anemia 2 (MLASA2) is an autosomal recessive mitochondrial disorder caused by pathogenic variants in YARS2. YARS2 variants confer heterogeneous phenotypes ranging from the full MLASA syndrome to a clinically unaffected state. Symptom onset is most common in the first decade of life but can occur in adulthood and has been reported following intercurrent illness. Early death can result from respiratory muscle weakness and cardiomyopathy. We report a case of MLASA2 with compound heterozygous YARS2 pathogenic variants; a known pathogenic nonsense variant [NM_001040436.3:c.98C>A (p.Ser33Ter)] and a likely pathogenic missense variant not previously associated with disease [NM_001040436.3:c.948G>T (p.Arg316Ser)]. The proband initially presented with a relatively mild phenotype of myopathy and lactic acidosis. During pregnancy, anemia emerged as an additional feature and in the postpartum period she experienced severe decompensation of cardiorespiratory function. This is the first reported case of pregnancy‐related complications in a patient with YARS2‐related mitochondrial disease. This case highlights the need for caution and careful counseling when considering pregnancy in mitochondrial disease, due to the risk of disease exacerbation and pregnancy complications.

Genomic diagnostics in polycystic kidney disease: an assessment of real-world use of whole-genome sequencing

Autosomal Dominant Polycystic Kidney Disease (ADPKD) is common, with a prevalence of 1/1000 and predominantly caused by disease-causing variants in PKD1 or PKD2. Clinical diagnosis is usually by age-dependent imaging criteria, which is challenging in patients with atypical clinical features, without family history, or younger age. However, there is increasing need for definitive diagnosis of ADPKD with new treatments available. Sequencing is complicated by six pseudogenes that share 97% homology to PKD1 and by recently identified phenocopy genes. Whole-genome sequencing can definitively diagnose ADPKD, but requires validation for clinical use. We initially performed a validation study, in which 42 ADPKD patients underwent sequencing of PKD1 and PKD2 by both whole-genome and Sanger sequencing, using a blinded, cross-over method. Whole-genome sequencing identified all PKD1 and PKD2 germline pathogenic variants in the validation study (sensitivity and specificity 100%). Two mosaic variants outside pipeline thresholds were not detected. We then examined the first 144 samples referred to a clinically-accredited diagnostic laboratory for clinical whole-genome sequencing, with targeted-analysis to a polycystic kidney disease gene-panel. In this unselected, diagnostic cohort (71 males :73 females), the diagnostic rate was 70%, including a diagnostic rate of 81% in patients with typical ADPKD (98% with PKD1/PKD2 variants) and 60% in those with atypical features (56% PKD1/PKD2; 44% PKHD1/HNF1B/GANAB/ DNAJB11/PRKCSH/TSC2). Most patients with atypical disease did not have clinical features that predicted likelihood of a genetic diagnosis. These results suggest clinicians should consider diagnostic genomics as part of their assessment in polycystic kidney disease, particularly in atypical disease.

Application of Genome Sequencing from Blood to Diagnose Mitochondrial Diseases

Mitochondrial diseases can be caused by pathogenic variants in nuclear or mitochondrial DNA-encoded genes that often lead to multisystemic symptoms and can have any mode of inheritance. Using a single test, Genome Sequencing (GS) can effectively identify variants in both genomes, but it has not yet been universally used as a first-line approach to diagnosing mitochondrial diseases due to related costs and challenges in data analysis. In this article, we report three patients with mitochondrial disease molecularly diagnosed through GS performed on DNA extracted from blood to demonstrate different diagnostic advantages of this technology, including the detection of a low-level heteroplasmic pathogenic variant, an intragenic nuclear DNA deletion, and a large mtDNA deletion. Current technical improvements and cost reductions are likely to lead to an expanded routine diagnostic usage of GS and of the complementary “Omic” technologies in mitochondrial diseases.

Diagnostic Yield of Whole Genome Sequencing After Nondiagnostic Exome Sequencing or Gene Panel in Developmental and Epileptic Encephalopathies

OBJECTIVE

To assess the benefits and limitations of whole genome sequencing (WGS) compared to exome sequencing (ES) or multigene panel (MGP) in the molecular diagnosis of developmental and epileptic encephalopathies (DEE).

METHODS

We performed WGS of 30 comprehensively phenotyped DEE patient trios that were undiagnosed after first-tier testing, including chromosomal microarray and either research ES (n = 15) or diagnostic MGP (n = 15).

RESULTS

Eight diagnoses were made in the 15 individuals who received prior ES (53%): 3 individuals had complex structural variants; 5 had ES-detectable variants, which now had additional evidence for pathogenicity. Eleven diagnoses were made in the 15 MGP-negative individuals (68%); the majority (n = 10) involved genes not included in the panel, particularly in individuals with postneonatal onset of seizures and those with more complex presentations including movement disorders, dysmorphic features, or multiorgan involvement. A total of 42% of diagnoses were autosomal recessive or X-chromosome linked.

CONCLUSION

WGS was able to improve diagnostic yield over ES primarily through the detection of complex structural variants (n = 3). The higher diagnostic yield was otherwise better attributed to the power of re-analysis rather than inherent advantages of the WGS platform. Additional research is required to assist in the assessment of pathogenicity of novel noncoding and complex structural variants and further improve diagnostic yield for patients with DEE and other neurogenetic disorders.

Recurrent SPECC1L-NTRK fusions in pediatric sarcoma and brain tumors

The identification of rearrangements driving expression of neurotrophic receptor tyrosine kinase (NTRK) family kinases in tumors has become critically important because of the availability of effective, specific inhibitor drugs. Whole-genome sequencing (WGS) combined with RNA sequencing (RNA-seq) can identify novel and recurrent expressed fusions. Here we describe three SPECC1L–NTRK fusions identified in two pediatric central nervous system cancers and an extracranial solid tumor using WGS and RNA-seq. These fusions arose either through a simple balanced rearrangement or in the context of a complex chromoplexy event. We cloned the SPECC1L–NTRK2 fusion directly from a patient sample and showed that enforced expression of this fusion is sufficient to promote cytokine-independent survival and proliferation. Cells transformed by SPECC1L–NTRK2 expression are sensitive to a TRK inhibitor drug. We report here that SPECC1L–NTRK fusions can arise in a range of pediatric cancers. Although WGS and RNA-seq are not required to detect NTRK fusions, these techniques may be of benefit when NTRK fusions are not suspected on clinical grounds or not identified by other methods.

Comprehensive Mutational and Phenotypic Characterization of New Metastatic Cutaneous Squamous Cell Carcinoma Cell Lines Reveal Novel Drug Susceptibilities

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer. Most patients who develop metastases (2–5%) present with advanced disease that requires a combination of radical surgery and adjuvant radiation therapy. There are few effective therapies for refractory disease. In this study, we describe novel patient-derived cell lines from cSCC metastases of the head and neck (designated UW-CSCC1 and UW-CSCC2). The cell lines genotypically and phenotypically resembled the original patient tumor and were tumorogenic in mice. Differences in cancer-related gene expression between the tumor and cell lines after various culturing conditions could be largely reversed by xenografting and reculturing. The novel drug susceptibilities of UW-CSCC1 and an irradiated subclone UW-CSCC1-R to drugs targeting cell cycle, PI3K/AKT/mTOR, and DNA damage pathways were observed using high-throughput anti-cancer and kinase-inhibitor compound libraries, which correlate with either copy number variations, targetable mutations and/or the upregulation of gene expression. A secondary screen of top hits in all three cell lines including PIK3CA-targeting drugs supports the utility of targeting the PI3K/AKT/mTOR pathway in this disease. UW-CSCC cell lines are thus useful preclinical models for determining targetable pathways and candidate therapeutics.

Whole genome, transcriptome and methylome profiling enhances actionable target discovery in high-risk pediatric cancer

The Zero Childhood Cancer Program is a precision medicine program to benefit children with poor-outcome, rare, relapsed or refractory cancer. Using tumor and germline whole genome sequencing (WGS) and RNA sequencing (RNAseq) across 252 tumors from high-risk pediatric patients with cancer, we identified 968 reportable molecular aberrations (39.9% in WGS and RNAseq, 35.1% in WGS only and 25.0% in RNAseq only). Of these patients, 93.7% had at least one germline or somatic aberration, 71.4% had therapeutic targets and 5.2% had a change in diagnosis. WGS identified pathogenic cancer-predisposing variants in 16.2% of patients. In 76 central nervous system tumors, methylome analysis confirmed diagnosis in 71.1% of patients and contributed to a change of diagnosis in two patients (2.6%). To date, 43 patients have received a recommended therapy, 38 of whom could be evaluated, with 31% showing objective evidence of clinical benefit. Comprehensive molecular profiling resolved the molecular basis of virtually all high-risk cancers, leading to clinical benefit in some patients.

NEMF mutations that impair ribosome-associated quality control are associated with neuromuscular disease

A hallmark of neurodegeneration is defective protein quality control. The E3 ligase Listerin (LTN1/Ltn1) acts in a specialized protein quality control pathway—Ribosome-associated Quality Control (RQC)—by mediating proteolytic targeting of incomplete polypeptides produced by ribosome stalling, and Ltn1 mutation leads to neurodegeneration in mice. Whether neurodegeneration results from defective RQC and whether defective RQC contributes to human disease have remained unknown. Here we show that three independently-generated mouse models with mutations in a different component of the RQC complex, NEMF/Rqc2, develop progressive motor neuron degeneration. Equivalent mutations in yeast Rqc2 selectively interfere with its ability to modify aberrant translation products with C-terminal tails which assist with RQC-mediated protein degradation, suggesting a pathomechanism. Finally, we identify NEMF mutations expected to interfere with function in patients from seven families presenting juvenile neuromuscular disease. These uncover NEMF’s role in translational homeostasis in the nervous system and implicate RQC dysfunction in causing neurodegeneration.

Defective protein quality control is a key feature of neurodegeneration. Here, the authors show that mutations in Nemf/NEMF, a component of the Ribosome-associated Quality Control complex, have a neurodegenerative effect in mice and may underlie neuromuscular disease in seven unrelated families.

Fatal perinatal mitochondrial cardiac failure caused by recurrent de novo duplications in the ATAD3 locus

Background

In about half of all patients with a suspected monogenic disease, genomic investigations fail to identify the diagnosis. A contributing factor is the difficulty with repetitive regions of the genome, such as those generated by segmental duplications. The ATAD3 locus is one such region, in which recessive deletions and dominant duplications have recently been reported to cause lethal perinatal mitochondrial diseases characterized by pontocerebellar hypoplasia or cardiomyopathy, respectively.

Methods

Whole exome, whole genome and long-read DNA sequencing techniques combined with studies of RNA and quantitative proteomics were used to investigate 17 subjects from 16 unrelated families with suspected mitochondrial disease.

Findings

We report six different de novo duplications in the ATAD3 gene locus causing a distinctive presentation including lethal perinatal cardiomyopathy, persistent hyperlactacidemia, and frequently corneal clouding or cataracts and encephalopathy. The recurrent 68 Kb ATAD3 duplications are identifiable from genome and exome sequencing but usually missed by microarrays. The ATAD3 duplications result in the formation of identical chimeric ATAD3A/ATAD3C proteins, altered ATAD3 complexes and a striking reduction in mitochondrial oxidative phosphorylation complex I and its activity in heart tissue.

Conclusions

ATAD3 duplications appear to act in a dominant-negative manner and the de novo inheritance infers a low recurrence risk for families, unlike most pediatric mitochondrial diseases. More than 350 genes underlie mitochondrial diseases. In our experience the ATAD3 locus is now one of the five most common causes of nuclear-encoded pediatric mitochondrial disease but the repetitive nature of the locus means ATAD3 diagnoses may be frequently missed by current genomic strategies.

Funding

Australian NHMRC, US Department of Defense, Japanese AMED and JSPS agencies, Australian Genomics Health Alliance and Australian Mito Foundation.

Development and validation of a targeted gene sequencing panel for application to disparate cancers

Next generation sequencing has revolutionised genomic studies of cancer, having facilitated the development of precision oncology treatments based on a tumour's molecular profile. We aimed to develop a targeted gene sequencing panel for application to disparate cancer types with particular focus on tumours of the head and neck, plus test for utility in liquid biopsy. The final panel designed through Roche/Nimblegen combined 451 cancer-associated genes (2.01 Mb target region). 136 patient DNA samples were collected for performance and application testing. Panel sensitivity and precision were measured using well-characterised DNA controls (n = 47), and specificity by Sanger sequencing of the Aryl Hydrocarbon Receptor Interacting Protein (AIP) gene in 89 patients. Assessment of liquid biopsy application employed a pool of synthetic circulating tumour DNA (ctDNA). Library preparation and sequencing were conducted on Illumina-based platforms prior to analysis with our accredited (ISO15189) bioinformatics pipeline. We achieved a mean coverage of 395x, with sensitivity and specificity of >99% and precision of >97%. Liquid biopsy revealed detection to 1.25% variant allele frequency. Application to head and neck tumours/cancers resulted in detection of mutations aligned to published databases. In conclusion, we have developed an analytically-validated panel for application to cancers of disparate types with utility in liquid biopsy.

Biparental inheritance of mitochondrial DNA in humans is not a common phenomenon

PURPOSE

A recent report has raised the possibility of biparental mitochondrial DNA (mtDNA) inheritance, which could lead to concerns by health-care professionals and patients regarding investigations and genetic counseling of families with pathogenic mitochondrial DNA variants. Our aim was to examine the frequency of this phenomenon by investigating a cohort of patients with suspected mitochondrial disease.

METHODS

We studied genome sequencing (GS) data of DNA extracted from blood samples of 41 pediatric patients with suspected mitochondrial disease and their parents.

RESULTS

All of the mtDNA variants in the probands segregated with their mother or were apparently de novo. There were no variants that segregated only with the father and none of these families showed evidence of biparental inheritance of their mtDNA.

CONCLUSION

Paternal mitochondrial transmission is unlikely to be a common occurrence and therefore at this point we would not recommend changes in clinical practice.

High Degree of Genetic Heterogeneity for Hereditary Cerebellar Ataxias in Australia

Genetic testing strategies such as next-generation sequencing (NGS) panels and whole genome sequencing (WGS) can be applied to the hereditary cerebellar ataxias (HCAs), but their exact role in the diagnostic pathway is unclear. We aim to determine the yield from genetic testing strategies and the genetic and phenotypic spectrum of HCA in Australia by analysing real-world data. We performed a retrospective review on 87 HCA cases referred to the Neurogenetics Clinic at the Royal North Shore Hospital, Sydney, Australia. Probands underwent triplet repeat expansion testing; those that tested negative had NGS-targeted panels and WGS testing when available. In our sample, 58.6% were male (51/87), with an average age at onset of 37.1 years. Individuals with sequencing variants had a prolonged duration of illness compared to those with a triplet repeat expansion. The detection rate in probands for routine repeat expansion panels was 13.8% (11/80). NGS-targeted panels yielded a further 11 individuals (11/32, 34.4%), with WGS yielding 1 more diagnosis (1/3, 33.3%). NGS panels and WGS improved the overall diagnostic rate to 28.8% (23/80) in 14 known HCA loci. The genetic findings included novel variants in ANO10, CACNA1A, PRKCG and SPG7. Our findings highlight the genetic heterogeneity of HCAs and support the use of NGS approaches for individuals who were negative on repeat expansion testing. In comparison to repeat disorders, individuals with sequencing variants may have a prolonged duration of illness, consistent with slower progression of disease.

Genomic stratification and liquid biopsy in a rare adrenocortical carcinoma (ACC) case, with dual lung metastases

Adrenocortical carcinoma is a rare malignancy with a poor prognosis and few treatment options. Molecular characterization of this cancer remains limited. We present a case of an adrenocortical carcinoma (ACC) in a 37-yr-old female, with dual lung metastases identified 1 yr following commencement of adjuvant mitotane therapy. As standard therapeutic regimens are often unsuccessful in ACC, we undertook a comprehensive genomic study into this case to identify treatment options and monitor disease progress. We performed targeted and whole-genome sequencing of germline, primary tumor, and both metastatic tumors from this patient and monitored recurrence over 2 years using liquid biopsy for ctDNA and steroid hormone measurements. Sequencing revealed the primary and metastatic tumors were hyperhaploid, with extensive loss of heterozygosity but few structural rearrangements. Loss-of-function mutations were identified in MSH2, TP53, RB1, and PTEN, resulting in tumors with mismatch repair signatures and microsatellite instability. At the cellular level, tumors were populated by mitochondria-rich oncocytes. Longitudinal ctDNA mutation and hormone profiles were unable to detect micrometastatic disease, consistent with clinical indicators of disease remission. The molecular signatures in our ACC case suggested immunotherapy in the event of disease progression; however, the patient remains free of cancer. The extensive molecular analysis presented here could be applied to other rare and/or poorly stratified cancers to identify novel or repurpose existing therapeutic options, thereby broadly improving diagnoses, treatments, and prognoses.

Deep multi-region whole-genome sequencing reveals heterogeneity and gene-by-environment interactions in treatment-naive, metastatic lung cancer

Our understanding of genomic heterogeneity in lung cancer is largely based on the analysis of early-stage surgical specimens. Here we used endoscopic sampling of paired primary and intrathoracic metastatic tumors from 11 lung cancer patients to map genomic heterogeneity inoperable lung cancer with deep whole-genome sequencing. Intra-patient heterogeneity in driver or targetable mutations was predominantly in the form of copy number gain. Private mutation signatures, including patterns consistent with defects in homologous recombination, were highly variable both within and between patients. Irrespective of histotype, we observed a smaller than expected number of private mutations, suggesting that ancestral clones accumulated large mutation burdens immediately prior to metastasis. Single-region whole-genome sequencing of from 20 patients showed that tumors in ever-smokers with the strongest tobacco signatures were associated with germline variants in genes implicated in the repair of cigarette-induced DNA damage. Our results suggest that lung cancer precursors in ever-smokers accumulate large numbers of mutations prior to the formation of frank malignancy followed by rapid metastatic spread. In advanced lung cancer, germline variants in DNA repair genes may interact with the airway environment to influence the pattern of founder mutations, whereas similar interactions with the tumor microenvironment may play a role in the acquisition of mutations following metastasis.

Reanalysis and optimisation of bioinformatic pipelines is critical for mutation detection

Rapid advances in genomic technologies have facilitated the identification pathogenic variants causing human disease. We report siblings with developmental and epileptic encephalopathy due to a novel, shared heterozygous pathogenic 13 bp duplication in SYNGAP1 (c.435_447dup, p.(L150Vfs*6)) that was identified by whole genome sequencing (WGS). The pathogenic variant had escaped earlier detection via two methodologies: whole exome sequencing and high-depth targeted sequencing. Both technologies had produced reads carrying the variant, however, they were either not aligned due to the size of the insertion or aligned to multiple major histocompatibility complex (MHC) regions in the hg19 reference genome, making the critical reads unavailable for variant calling. The WGS pipeline followed different protocols, including alignment of reads to the GRCh37 reference genome, which lacks the additional MHC contigs. Our findings highlight the benefit of using orthogonal clinical bioinformatic pipelines and all relevant inheritance patterns to re-analyze genomic data in undiagnosed patients.

Expanding the spectrum of PEX16 mutations and novel insights into disease mechanisms

Zellweger syndrome spectrum disorders are caused by mutations in any of at least 12 different PEX genes. This includes PEX16, an important regulator of peroxisome biogenesis. Using whole genome sequencing, we detected previously unreported, biallelic variants in PEX16 [NM_004813.2:c.658G>A, p.(Ala220Thr) and NM_004813.2:c.830G>A, p.(Arg277Gln)] in an individual with leukodystrophy, spastic paraplegia, cerebellar ataxia, and craniocervical dystonia with normal plasma very long chain fatty acids. Using olfactory-neurosphere derived cells, a population of neural stem cells, we showed patient cells had reduced peroxisome density and increased peroxisome size, replicating previously reported findings in PEX16 cell lines. Along with alterations in peroxisome morphology, patient cells also had impaired peroxisome function with reduced catalase activity. Furthermore, patient cells had reduced oxidative stress levels after exposure to hydrogen-peroxide (H2O2), which may be a result of compensation by H2O2 metabolising enzymes other than catalase to preserve peroxisome-related cell functions. Our findings of impaired catalase activity and altered oxidative stress response are novel. Our study expands the phenotype of PEX16 mutations by including dystonia and provides further insights into the pathological mechanisms underlying PEX16-associated disorders. Additional studies of the full spectrum of peroxisomal dysfunction could improve our understanding of the mechanism underlying PEX16-associated disorders.