A universal molecular control for DNA, mRNA and protein expression

The expression of genes encompasses their transcription into mRNA followed by translation into protein. In recent years, next-generation sequencing and mass spectrometry methods have profiled DNA, RNA and protein abundance in cells. However, there are currently no reference standards that are compatible across these genomic, transcriptomic and proteomic methods, and provide an integrated measure of gene expression. Here, we use synthetic biology principles to engineer a multi-omics control, termed pREF, that can act as a universal molecular standard for next-generation sequencing and mass spectrometry methods. The pREF sequence encodes 21 synthetic genes that can be in vitro transcribed into spike-in mRNA controls, and in vitro translated to generate matched protein controls. The synthetic genes provide qualitative controls that can measure sensitivity and quantitative accuracy of DNA, RNA and peptide detection. We demonstrate the use of pREF in metagenome DNA sequencing and RNA sequencing experiments and evaluate the quantification of proteins using mass spectrometry. Unlike previous spike-in controls, pREF can be independently propagated and the synthetic mRNA and protein controls can be sustainably prepared by recipient laboratories using common molecular biology techniques. Together, this provides a universal synthetic standard able to integrate genomic, transcriptomic and proteomic methods.

Quantitative ctDNA Detection in Hepatoblastoma: Implications for Precision Medicine

Hepatoblastoma is characterized by driver mutations in CTNNB1, making it an attractive biomarker for a liquid biopsy approach utilizing circulating tumor DNA (ctDNA). This prospective observational study sought to ascertain the feasibility of ctDNA detection in patients with hepatoblastoma and explore its associations with established clinical indicators and biomarkers, including serum Alpha-fetoprotein (AFP). We obtained 38 plasma samples and 17 tumor samples from 20 patients with hepatoblastoma. These samples were collected at various stages: 10 at initial diagnosis, 17 during neoadjuvant chemotherapy, 6 post-operatively, and 5 at disease recurrence. Utilizing a bespoke sequencing assay we developed called QUENCH, we identified single nucleotide variants and deletions in CTNNB1 ctDNA. Our study demonstrated the capability to quantitate ctDNA down to a variant allele frequency of 0.3%, achieving a sensitivity of 90% for patients at initial diagnosis, and a specificity of 100% at the patient level. Notably, ctDNA positivity correlated with tumor burden, and ctDNA levels exhibited associations with macroscopic residual disease and treatment response. Our findings provide evidence for the utility of quantitative ctDNA detection in hepatoblastoma management. Given the distinct detection targets, ctDNA and AFP-based stratification and monitoring approaches could synergize to enhance clinical decision-making. Further research is needed to elucidate the interplay between ctDNA and AFP and determine the optimal clinical applications for both methods in risk stratification and residual disease detection.

Circulating Tumour DNA Biomarkers Associated with Outcomes in Metastatic Prostate Cancer Treated with Lutetium-177-PSMA-617

Background

Lutetium-177-prostate-specific membrane antigen- 617 (Lu-PSMA) is an effective therapy for metastatic castration-resistant prostate cancer (mCRPC). However, treatment responses are heterogeneous despite stringent positron emission tomography (PET)-based imaging selection criteria. Molecularly based biomarkers have potential to refine patient selection and optimise outcomes.

Objective

To identify circulating tumour DNA (ctDNA) features associated with treatment outcomes for men treated with Lu-PSMA.

Design setting and participants

ctDNA from men treated with Lu-PSMA in combination with idronoxil for progressive mCRPC were analysed using an 85-gene customised sequencing assay. ctDNA fractions, molecular profiles, and the presence of alterations in aggressive-variant prostate cancer (AVPC) genes were analysed at baseline, cycle 3 and at disease progression.

Intervention

Men received Lu-PSMA with idronoxil every 6 wk for up to six cycles.

Outcome measurements and statistical analysis

Baseline and exit PSMA and fluorodeoxyglucose PET/computed tomography (CT) imaging was conducted at baseline and study exit. Single-photon emission CT (SPECT) scans were performed 24 h after Lu-PSMA. Blood samples were collected at baseline,cycle 3 and at disease progression. Cox proportional-hazards models were used to assess associations and derive hazard ratios (HRs) and confidence intervals (CIs) for associations between molecular factors, imaging features, and clinical outcomes.

Results and limitations

Sixty samples from 32 men were sequenced (32 at baseline, 24 at cycle 3, four from patients with disease progression); two samples (baseline, on-treatment) from one individual were excluded from analysis owing to poor quality of the baseline sequencing data. Alterations in AVPC genes were associated with shorter prostate-specific antigen (PSA) progression-free survival (PFS) and overall survival (OS) in univariate (HR 3.4, 95% CI 1.5-7.7; p = 0.0036; and HR 3.3, 95% CI 1.4-7.7; p = 0.0063, respectively) and multivariate analyses (HR 4.8, 95% CI 1.8-13; p = 0.0014; and HR 4.1, 95% CI 1.6-11; p = 0.004).

Conclusions

ctDNA alterations in AVPC genes were associated with shorter PSA PFS and OS among men treated with Lu-PSMA and intermittent idronoxil. These candidate molecular biomarkers warrant further study to determine whether they have predictive value and potential to guide synergistic combination strategies to enhance outcomes for men treated with Lu-PSMA for mCRPC.

Patient summary

Certain DNA/gene changes detected in the blood of men with advanced prostate cancer were associated with shorter benefit from lutetium PSMA, a targeted radioactive therapy. This information may be useful in determining which men may benefit most from this treatment, but additional research is needed.

Community Genetics screening in a pandemic: solutions for pre-test education, informed consent, and specimen collection

A Community Genetics carrier screening program for the Jewish community has operated on-site in high schools in Sydney (Australia) for 25 years. During 2020, in response to the COVID-19 pandemic, government-mandated social-distancing, 'lock-down' public health orders, and laboratory supply-chain shortages prevented the usual operation and delivery of the annual testing program. We describe development of three responses to overcome these challenges: (1) pivoting to online education sufficient to ensure informed consent for both genetic and genomic testing; (2) development of contactless telehealth with remote training and supervision for collecting genetic samples using buccal swabs; and (3) a novel patient and specimen identification 'GeneTrustee' protocol enabling fully identified clinical-grade specimens to be collected and DNA extracted by a research laboratory while maintaining full participant confidentiality and privacy. These telehealth strategies for education, consent, specimen collection and sample processing enabled uninterrupted delivery and operation of complex genetic testing and screening programs even amid pandemic restrictions. These tools remain available for future operation and can be adapted to other programs.

Low disease risk and penetrance in Leber hereditary optic neuropathy

The risk of Leber hereditary optic neuropathy (LHON) has largely been extrapolated from disease cohorts, which underestimate the population prevalence of pathogenic primary LHON variants as a result of incomplete disease penetrance. Understanding the true population prevalence of primary LHON variants, alongside the rate of clinical disease, provides a better understanding of disease risk and variant penetrance. We identified pathogenic primary LHON variants in whole-genome sequencing data of a well-characterized population-based control cohort and found that the prevalence is far greater than previously estimated, as it occurs in approximately 1 in 800 individuals. Accordingly, we were able to more accurately estimate population risk and disease penetrance in LHON variant carriers, validating our findings by using other large control datasets. These findings will inform accurate counseling in relation to the risk of vision loss in LHON variant carriers and disease manifestation in their family. This Matters Arising paper is in response to Lopez Sanchez et al. (2021), published in The American Journal of Human Genetics. See also the response by Mackey et al. (2022), published in this issue.

Pharmacogenomic testing: perception of clinical utility, enablers and barriers to adoption in Australian hospitals

Background

Despite healthcare professionals (HCP) endorsing the clinical utility of pharmacogenomics testing, use in clinical practice is limited.

Aims

To assess HCP' perceptions of pharmacogenomic testing and identify barriers to implementation.

Methods

HCP involved in prescribing decisions at three hospitals in Sydney, Australia, were invited to participate. The online survey assessed perceptions of pharmacogenomic testing, including: (i) demographic and practice variables; (ii) use, knowledge and confidence; (iii) perceived benefits; (iv) barriers to implementation; and (v) operational and/or system changes and personnel required to implement on site.

Results

HCP were predominantly medical practitioners (75/107) and pharmacists (25/107). HCP perceived pharmacogenomic testing was beneficial to identify reasons for drug intolerance (85/95) and risk of side‐effects (86/95). Although testing was considered relevant to their practice (79/100), few HCP (23/100) reported past or intended future use (26/100). Few HCP reported confidence in their ability to identify indications for pharmacogenomic testing (14/107), order tests (19/106) and communicate results with patients (16/107). Lack of clinical practice guidelines (62/79) and knowledge (54/77) were identified as major barriers to implementation of pharmacogenomics. Comprehensive reimbursement for testing and clinical practice guidelines, alongside models‐of‐care involving multidisciplinary teams and local clinical champions were suggested as strategies to facilitate implementation of pharmacogenomic testing into practice.

Conclusions

Pharmacogenomic testing was considered important to guide drug selection and dosing decisions. However, limited knowledge, low confidence and an absence of guidelines impede the use of pharmacogenomic testing. Establishment of local resources including multidisciplinary models‐of‐care was suggested to facilitate implementation of pharmacogenomics.

Genetic Pathways Associated With Sleep Problems in Children With Autism Spectrum Disorder

AIMS

Children on the autism spectrum are more likely to have sleep problems than non-autistic children. Sleep disturbance may exacerbate emotional and behavioral problems of children on the autism spectrum. A better understanding of the biological mechanisms underlying sleep disturbance provide clues to better management for this co-morbid condition in autism. The goal of the current study is to identify genetic variants associated with sleep disturbance and melatonin levels in autistic children.

METHODS

A total of 969 children on the autism spectrum were genotyped using the Global Screening Array v1 or Global Screening Array v2. Sleep problems were assessed using the Children's Sleep Habits Questionnaire (CSHQ). Melatonin levels were measured using the urine samples of 219 probands. The relationship between the melatonin level and CSHQ score was examined using the general linear model. The genetic variants associated with the CSHQ score and melatonin level as two separate quantitative traits were determined using genomewide association studies.

RESULTS

The data indicates that urine melatonin levels were positively associated with CSHQ scores, suggesting that autistic children with a poorer sleep qualiy could has higher melatonin level. Furthermore, genetic assocication studies suggest that genetic pathways involved in pro-inflammatory responses might be involved in sleep disturbance, while genetic pathways involved in catecholamine-secreting PC12 cells and Schwann cells could be associated with melatonin levels.

CONCLUSIONS

Taken together, our findings indicate that sleep disturbance and melatonin metabolism could be attributable to distinct biological mechanisms in autistic children since they might not share genetic contributors.

Abstract 582: Utility of CTNNB1 ctDNA as a biomarker for hepatoblastoma

Driver mutations in CTNNB1 are a hallmark of hepatoblastoma (HB) and offer a common biomarker for a liquid biopsy approach that is based on the presence of CTNNB1 circulating tumor DNA (ctDNA). Initial results 1 showed an association of CTNNB1 ctDNA variant allele frequency (VAF) identified through digital droplet PCR (ddPCR) and the levels of serum Alpha-fetoprotein (AFP, current HB biomarker) throughout the course of treatment of three patients with HB. However, applying custom probe-designed ddPCR assays may pose challenges for real-time data at diagnosis. Our primary objective is to investigate the utility of a universal next-generation sequencing (NGS) assay to detect CTNNB1 ctDNA in patients with HB. Our secondary objective is to compare the levels of: (1) NGS ctDNA, (2) ddPCR ctDNA, (3) AFP, and verify if they correlate with tumor burden and treatment response. We developed and tested a custom NGS assay covering exons 2 to 4 of CTNNB1 for detection of somatic mutations in plasma. We used the QIASeq platform which: (1) is based on a single primer extension and, thus, captures both single nucleotide variants (SNV) and structural variants (SV), and (2) integrates unique molecular indexing (UMI), thus allowing determination of ctDNA levels from VAF. Our cohort included 18 patients, 13 with Sanger sequencing confirmed CTNNB1 somatic mutations in the matched tumors: 6 missense SNVs and 7 deletions ranging from 96 to 348 base pairs. Our NGS assay was able to accurately detect mutations in 9/13 (69%, 7 SV and 2 SNV) of the Sanger confirmed cases, and SV and missense SNV in 2/5 (40%) cases pending orthogonal Sanger confirmation. Available ddPCR data for 3 patients (8 samples) 1 and corresponding NGS data showed similar ctDNA levels (deviation ranging between 1.2-8.5%). NGS CTNNB1 variants were detected in samples taken at initial diagnosis in 8/9 (89%), following neoadjuvant chemotherapy in 4/9 (44%), post-operatively for fully resected localized disease in 0/4 (0%), and metastatic recurrence in 1/5 (20%), suggesting a positive correlation between CTNNB1 variants and tumor burden. AFP levels were available for 24/27 samples, with abnormal levels detected in 23 samples. The ctDNA levels determined from NGS did not correlate with AFP levels; however, longitudinal plasma samples (n=9) showed similar dynamics of both ctDNA and AFP, reflecting the response to treatment. To conclude, we show that ctDNA detectable with our universal HB NGS assay is a good surrogate marker of tumor burden and shows correlation with treatment response. SV, prevalent in approximately 50% of cases of HB, may be inadequately captured and thereby misrepresented in probe-capture NGS platforms but were highly detectable in our assay. Further work is needed to understand the interplay of both AFP and ctDNA in HB monitoring. Reference: 1. Kahana-Edwin S, McCowage G, Cain L, et al. Exploration of CTNNB1 ctDNA as a putative biomarker for hepatoblastoma. Pediatr Blood Cancer. 2020. doi:10.1002/pbc.28594

Citation Format: Smadar Kahana-Edwin, Andre E. Minoche, Lucy E. Cain, Geoffrey McCowage, Sarah E. Woodfield, Sanjeev A. Vasudevan, Sarah Kummerfeld, Jonathan Karpelowsky. Utility of CTNNB1 ctDNA as a biomarker for hepatoblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 582.

ClinSV: clinical grade structural and copy number variant detection from whole genome sequencing data

Whole genome sequencing (WGS) has the potential to outperform clinical microarrays for the detection of structural variants (SV) including copy number variants (CNVs), but has been challenged by high false positive rates. Here we present ClinSV, a WGS based SV integration, annotation, prioritization, and visualization framework, which identified 99.8% of simulated pathogenic ClinVar CNVs > 10 kb and 11/11 pathogenic variants from matched microarrays. The false positive rate was low (1.5-4.5%) and reproducibility high (95-99%). In clinical practice, ClinSV identified reportable variants in 22 of 485 patients (4.7%) of which 35-63% were not detectable by current clinical microarray designs. ClinSV is available at https://github.com/KCCG/ClinSV .

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.

Whole genome sequencing for the genetic diagnosis of heterogenous dystonia phenotypes

INTRODUCTION

Dystonia is a clinically and genetically heterogeneous disorder and a genetic cause is often difficult to elucidate. This is the first study to use whole genome sequencing (WGS) to investigate dystonia in a large sample of affected individuals.

METHODS

WGS was performed on 111 probands with heterogenous dystonia phenotypes. We performed analysis for coding and non-coding variants, copy number variants (CNVs), and structural variants (SVs). We assessed for an association between dystonia and 10 known dystonia risk variants.

RESULTS

A genetic diagnosis was obtained for 11.7% (13/111) of individuals. We found that a genetic diagnosis was more likely in those with an earlier age at onset, younger age at testing, and a combined dystonia phenotype. We identified pathogenic/likely-pathogenic variants in ADCY5 (n = 1), ATM (n = 1), GNAL (n = 2), GLB1 (n = 1), KMT2B (n = 2), PRKN (n = 2), PRRT2 (n = 1), SGCE (n = 2), and THAP1 (n = 1). CNVs were detected in 3 individuals. We found an association between the known risk variant ARSG rs11655081 and dystonia (p = 0.003).

CONCLUSION

A genetic diagnosis was found in 11.7% of individuals with dystonia. The diagnostic yield was higher in those with an earlier age of onset, younger age at testing, and a combined dystonia phenotype. WGS may be particularly relevant for dystonia given that it allows for the detection of CNVs, which accounted for 23% of the genetically diagnosed cases.

IRF2 transcriptionally induces GSDMD expression for pyroptosis

Gasdermin-D (GSDMD) is cleaved by caspase-1, caspase-4, and caspase-11 in response to canonical and noncanonical inflammasome activation. Upon cleavage, GSDMD oligomerizes and forms plasma membrane pores, resulting in interleukin-1β (IL-1β) secretion, pyroptotic cell death, and inflammatory pathologies, including periodic fever syndromes and septic shock—a plague on modern medicine. Here, we showed that IRF2, a member of the interferon regulatory factor (IRF) family of transcription factors, was essential for the transcriptional activation of GSDMD. A forward genetic screen with N-ethyl-N-nitrosourea (ENU)–mutagenized mice linked IRF2 to inflammasome signaling. GSDMD expression was substantially attenuated in IRF2-deficient macrophages, endothelial cells, and multiple tissues, which corresponded with reduced IL-1β secretion and inhibited pyroptosis. Mechanistically, IRF2 bound to a previously uncharacterized but unique site within the GSDMD promoter to directly drive GSDMD transcription for the execution of pyroptosis. Disruption of this single IRF2-binding site abolished signaling by both the canonical and noncanonical inflammasomes. Together, our data illuminate a key transcriptional mechanism for expression of the gene encoding GSDMD, a critical mediator of inflammatory pathologies.

Intrinsic apoptosis shapes the tumor spectrum linked to inactivation of the deubiquitinase BAP1

Malignancies arising from mutation of tumor suppressors have unexplained tissue proclivity. For example, BAP1 encodes a widely expressed deubiquitinase for histone H2A, but germline mutations are predominantly associated with uveal melanomas and mesotheliomas. We show that BAP1 inactivation causes apoptosis in mouse embryonic stem cells, fibroblasts, liver, and pancreatic tissue but not in melanocytes and mesothelial cells. Ubiquitin ligase RNF2, which silences genes by monoubiquitinating H2A, promoted apoptosis in BAP1-deficient cells by suppressing expression of the prosurvival genes Bcl2 and Mcl1. In contrast, BAP1 loss in melanocytes had little impact on expression of prosurvival genes, instead inducing Mitf Thus, BAP1 appears to modulate gene expression by countering H2A ubiquitination, but its loss only promotes tumorigenesis in cells that do not engage an RNF2-dependent apoptotic program.

Caspase-11 cleaves gasdermin D for non-canonical inflammasome signalling

Intracellular lipopolysaccharide from Gram-negative bacteria including Escherichia coli, Salmonella typhimurium, Shigella flexneri, and Burkholderia thailandensis activates mouse caspase-11, causing pyroptotic cell death, interleukin-1β processing, and lethal septic shock. How caspase-11 executes these downstream signalling events is largely unknown. Here we show that gasdermin D is essential for caspase-11-dependent pyroptosis and interleukin-1β maturation. A forward genetic screen with ethyl-N-nitrosourea-mutagenized mice links Gsdmd to the intracellular lipopolysaccharide response. Macrophages from Gsdmd(-/-) mice generated by gene targeting also exhibit defective pyroptosis and interleukin-1β secretion induced by cytoplasmic lipopolysaccharide or Gram-negative bacteria. In addition, Gsdmd(-/-) mice are protected from a lethal dose of lipopolysaccharide. Mechanistically, caspase-11 cleaves gasdermin D, and the resulting amino-terminal fragment promotes both pyroptosis and NLRP3-dependent activation of caspase-1 in a cell-intrinsic manner. Our data identify gasdermin D as a critical target of caspase-11 and a key mediator of the host response against Gram-negative bacteria.

Gene expression patterns of Th2 inflammation and intercellular communication in asthmatic airways

Asthma is canonically thought of as a disorder of excessive Th2-driven inflammation in the airway, although recent studies have described heterogeneity with respect to asthma pathophysiology. We have previously described distinct phenotypes of asthma based on the presence or absence of a three-gene "Th2 signature" in bronchial epithelium, which differ in terms of eosinophilic inflammation, mucin composition, subepithelial fibrosis, and corticosteroid responsiveness. In the present analysis, we sought to describe Th2 inflammation in human asthmatic airways quantitatively with respect to known mediators of inflammation and intercellular communication. Using whole-genome microarray and quantitative real-time PCR analysis of endobronchial biopsies from 27 mild-to-moderate asthmatics and 13 healthy controls with associated clinical and demographic data, we found that asthmatic Th2 inflammation is expressed over a variable continuum, correlating significantly with local and systemic measures of allergy and eosinophilia. We evaluated a composite metric describing 79 coexpressed genes associated with Th2 inflammation against the biological space comprising cytokines, chemokines, and growth factors, identifying distinctive patterns of inflammatory mediators as well as Wnt, TGF-β, and platelet-derived growth factor family members. This integrated description of the factors regulating inflammation, cell migration, and tissue remodeling in asthmatic airways has important consequences for the pathophysiological and clinical impacts of emerging asthma therapeutics targeting Th2 inflammation.

The genome of the social amoeba Dictyostelium discoideum

The social amoebae are exceptional in their ability to alternate between unicellular and multicellular forms. Here we describe the genome of the best-studied member of this group, Dictyostelium discoideum. The gene-dense chromosomes of this organism encode approximately 12,500 predicted proteins, a high proportion of which have long, repetitive amino acid tracts. There are many genes for polyketide synthases and ABC transporters, suggesting an extensive secondary metabolism for producing and exporting small molecules. The genome is rich in complex repeats, one class of which is clustered and may serve as centromeres. Partial copies of the extrachromosomal ribosomal DNA (rDNA) element are found at the ends of each chromosome, suggesting a novel telomere structure and the use of a common mechanism to maintain both the rDNA and chromosomal termini. A proteome-based phylogeny shows that the amoebozoa diverged from the animal-fungal lineage after the plant-animal split, but Dictyostelium seems to have retained more of the diversity of the ancestral genome than have plants, animals or fungi.