Search

Background and ObjectivesMitochondrial 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.MethodsAdult 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.ResultsOf 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.DiscussionComprehensive bigenomic sequencing accurately detects causative genetic variants in affected MD patients, simplifying diagnosis, enabling early treatment, and informing the risk of genetic transmission.

Cerebral palsy (CP) is the most common cause of childhood physical disability, with incidence between 1/500 and 1/700 births in the developed world. Despite increasing evidence for a major contribution of genetics to CP aetiology, genetic testing is currently not performed systematically. We assessed the diagnostic rate of genome sequencing (GS) in a clinically unselected cohort of 150 singleton CP patients, with CP confirmed at >4 years of age. Clinical grade GS was performed on the proband and variants were filtered, and classified according to American College of Medical Genetics and Genomics–Association for Molecular Pathology (ACMG-AMP) guidelines. Variants classified as pathogenic or likely pathogenic (P/LP) were further assessed for their contribution to CP. In total, 24.7% of individuals carried a P/LP variant(s) causing or increasing risk of CP, with 4.7% resolved by copy number variant analysis and 20% carrying single nucleotide or indel variants. A further 34.7% carried one or more rare, high impact variants of uncertain significance (VUS) in variation intolerant genes. Variants were identified in a heterogeneous group of genes, including genes associated with hereditary spastic paraplegia, clotting and thrombophilic disorders, small vessel disease, and other neurodevelopmental disorders. Approximately 1/2 of individuals were classified as likely to benefit from changed clinical management as a result of genetic findings. In addition, no significant association between genetic findings and clinical factors was detectable in this cohort, suggesting that systematic sequencing of CP will be required to avoid missed diagnoses.

Driver mutations in CTNNB1 are a hallmark of hepatoblastoma and offer a common biomarker for a liquid biopsy approach based on the presence of CTNNB1 circulating tumor DNA (ctDNA). We developed and investigated the utility of a quantitative universal next-generation sequencing (NGS) ctDNA assay for hepatoblastoma (QUENCH) to detect CTNNB1 ctDNA and assessed the links between ctDNA and current clinical indicators/biomarkers in hepatoblastoma. Applied to patients with hepatoblastoma, we demonstrate quantitation of various variants including single base substitutions and deletions down to 0.3% variant allele frequency, with 65% sensitivity and 100% specificity at the patient level, to allow biopsy-free tumor genotyping and sensitive ctDNA quantitation. CtDNA positivity correlates with tumor burden and ctDNA levels correlate with macroscopic residual disease and treatment response, thus providing promising evidence for the utility of quantitative ctDNA detection in hepatoblastoma.

Chard (Beta vulgaris ssp. vulgaris var. cicla) is a member of one of four different cultigroups of beets. While the genome of sugar beet, the most prominent beet crop, has been studied extensively, molecular data on other beet cultivars is scant. Here, we present a genome assembly of chard, a vegetable crop grown for its fleshy leaves. We report a de novo genome assembly of 604 Mbp, slightly larger than sugar beet assemblies presented so far. About 57 % of the assembly was annotated as repetitive sequence, of which LTR retrotransposons were the most abundant. Based on the presence of conserved genes, the chard assembly was estimated to be at least 96 % complete regarding its gene space. We predicted 34,521 genes of which 27,582 genes were supported by evidence from transcriptomic sequencing reads, and 5503 of the evidence-supported genes had multiple isoforms. We compared the chard gene set with gene sets from sugar beet and two wild beets (i.e. Beta vulgaris ssp. maritima and Beta patula) to find orthology relationships and identified genome-wide syntenic regions between chard and sugar beet. Lastly, we determined genomic variants that distinguish sugar beet and chard. Assessing the variation distribution along the chard chromosomes, we found extensive haplotype sharing between the two cultivars. In summary, our work provides a foundation for the molecular analysis of Beta vulgaris cultigroups as a basis for chard genomics and to unravel the domestication history of beet crops.

ObjectiveTo 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).MethodsWe 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).ResultsEight 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.ConclusionWGS 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.

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.

Dogs are uniquely associated with human dispersal and bring transformational insight into the domestication process. Dingoes represent an intriguing case within canine evolution being geographically isolated for thousands of years. Here, we present a high-quality de novo assembly of a pure dingo (CanFam_DDS). We identified large chromosomal differences relative to the current dog reference (CanFam3.1) and confirmed no expanded pancreatic amylase gene as found in breed dogs. Phylogenetic analyses using variant pairwise matrices show that the dingo is distinct from five breed dogs with 100% bootstrap support when using Greenland wolf as the outgroup. Functionally, we observe differences in methylation patterns between the dingo and German shepherd dog genomes and differences in serum biochemistry and microbiome makeup. Our results suggest that distinct demographic and environmental conditions have shaped the dingo genome. In contrast, artificial human selection has likely shaped the genomes of domestic breed dogs after divergence from the dingo.

Interpretation of the significance of maternally inherited X chromosome variants in males with neurocognitive phenotypes continues to present a challenge to clinical geneticists and diagnostic laboratories. Here we report 14 males from 9 families with duplications at the Xq13.2-q13.3 locus with a common facial phenotype, intellectual disability (ID), distinctive behavioral features, and a seizure disorder in two cases. All tested carrier mothers had normal intelligence. The duplication arose de novo in three mothers where grandparental testing was possible. In one family the duplication segregated with ID across three generations. RLIM is the only gene common to our duplications. However, flanking genes duplicated in some but not all the affected individuals included the brain-expressed genes NEXMIF, SLC16A2, and the long non-coding RNA gene FTX. The contribution of the RLIM-flanking genes to the phenotypes of individuals with different size duplications has not been fully resolved. Missense variants in RLIM have recently been identified to cause X-linked ID in males, with heterozygous females typically having normal intelligence and highly skewed X chromosome inactivation. We detected consistent and significant increase of RLIM mRNA and protein levels in cells derived from seven affected males from five families with the duplication. Subsequent analysis of MDM2, one of the targets of the RLIM E3 ligase activity, showed consistent downregulation in cells from the affected males. All the carrier mothers displayed normal RLIM mRNA levels and had highly skewed X chromosome inactivation. We propose that duplications at Xq13.2-13.3 including RLIM cause a recognizable but mild neurocognitive phenotype in hemizygous males.

Purpose: The utility of genome sequencing (GS) in the diagnosis of suspected pediatric mitochondrial disease (MD) was investigated. Methods: An Australian cohort of 40 pediatric patients with clinical features suggestive of MD were classified using the modified Nijmegen mitochondrial disease severity scoring into definite (17), probable (17), and possible (6) MD groups. Trio GS was performed using DNA extracted from patient and parent blood. Data were analyzed for single-nucleotide variants, indels, mitochondrial DNA variants, and structural variants. Results: A definitive MD gene molecular diagnosis was made in 15 cases and a likely MD molecular diagnosis in a further five cases. Causative mitochondrial DNA (mtDNA) variants were identified in four of these cases. Three potential novel MD genes were identified. In seven cases, causative variants were identified in known disease genes with no previous evidence of causing a primary MD. Diagnostic rates were higher in patients classified as having definite MD. Conclusion: GS efficiently identifies variants in MD genes of both nuclear and mitochondrial origin. A likely molecular diagnosis was identified in 67% of cases and a definitive molecular diagnosis achieved in 55% of cases. This study highlights the value of GS for a phenotypically and genetically heterogeneous disorder like MD.

The inherited retinal dystrophies (IRDs) are a clinically and genetically complex group of disorders primarily affecting the rod and cone photoreceptors or other retinal neuronal layers, with emerging therapies heralding the need for accurate molecular diagnosis. Targeted capture and panel-based strategies examining the partial or full exome deliver molecular diagnoses in many IRD families tested. However, approximately one in three families remain unsolved and unable to obtain personalised recurrence risk or access to new clinical trials or therapy. In this study, we investigated whole genome sequencing (WGS), focused assays and functional studies to assist with unsolved IRD cases and facilitate integration of these approaches to a broad molecular diagnostic clinical service. The WGS approach identified variants not covered or underinvestigated by targeted capture panel-based clinical testing strategies in six families. This included structural variants, with notable benefit of the WGS approach in repetitive regions demonstrated by a family with a hybrid gene and hemizygous missense variant involving the opsin genes, OPN1LW and OPN1MW. There was also benefit in investigation of the repetitive GC-rich ORF15 region of RPGR. Further molecular investigations were facilitated by focused assays in these regions. Deep intronic variants were identified in IQCB1 and ABCA4, with functional RNA based studies of the IQCB1 variant revealing activation of a cryptic splice acceptor site. While targeted capture panel-based methods are successful in achieving an efficient molecular diagnosis in a proportion of cases, this study highlights the additional benefit and clinical value that may be derived from WGS, focused assays and functional genomics in the highly heterogeneous IRDs.

Monogenic diseases can often manifest diverse clinical phenotypes and cause diagnostic dilemmas. While monoallelic loss-of-function variants in TNFAIP3 (Haploinsufficiency of A20; HA20) cause a highly penetrant autoinflammatory disease, the variable expressivity suggest a role for additional genetic and environmental disease modifiers. Here, we identify critically ill children who inherited a family-specific TNFAIP3 deletion from one of their otherwise healthy parents. Each of the probands also inherited in trans a subtle loss-of-function I207L TNFAIP3 variant that is common in Oceania, originally introgressed from Denisovans. Modelling this compound heterozgous state in mice under specific pathogen free conditions demonstrated a reduced threshold to break immune tolerance. Exaggerated immune responses were precipitated by inheriting the two genetic hits on the TNFAIP3 checkpoint coupled with increasing the microbial challenge to immune tolerance, either by co-housing with pet store mice carrying a wild microbial burden or by transient dietary exposure to a chemical that diminishes the intestinal mucin barrier separating gut microbes from immune sensing systems. These data illuminate second-hit genetic and environmental modifiers contributing to complex inflammatory and autoimmune disease. Increased mechanistic understanding of the presence and contribution of disease modifiers will aid diagnostic and prognostic patient stratification and potentially reveal novel therapeutic opportunities.

Additional file 1: Fig. S1. number of citations from PubMed for breast cancer cell lines. Number of citations for each of the cell lines obtained from PubMed, data retrieved from Pubmed in May 2022. Figure S2. Heatmap of correlation of genotyping calls of breast cancer cell lines from three different studies. Correlation analysis of genotyping calls from this study to compare with SNP calls from Heiser et al and CCLE. Figure S3. barplots of number of different types of variants identified in the cell lines. Figure S4. IGV plot of a non-coding gene MALAT1 in MCF7 and MDAMB231 respectively, showing a mutation in MCF7 but not in MDAMB231. Figure S5. Circos plot of structural variations in breast cancer cell lines, MDA-MB-231, T47D, MDA-MB-468 and SKBR3. Arcs connecting two loci of difference chromosomes indicate inter-chromosomal structural variations. Figure S6. GREAT analysis of genes affected by SV variants in the breast cancer cell lines. Figure S7. Representative IGV plot showing copy number gains in FOXA1, GATA3 and ESR1 in MCF7. Tracks from top to bottom: depth of coverage in an NA12878 control (control), all reads in the sample (all reads), or reads with mapping quality >=20 (MQ>20), the average mapping quality of aligned reads from the sample (MQ, if no reads align MQ=0), coverage standard deviation from 500 controls (Coverage SD, indicating common CNV), overlapping segmental duplications published by Bailey JA et al. 2002 (SEG-DUP, used as control for germline CNVs), discordant pairs (DP), split reads (SR), variants from the Database of Genomic Variants (DGV), and RefSeq genes (Genes). Figure S8. barplots of number of different types of variants identified in the six PDX models. Figure S9. Kaplan-Meier survival analysis of METABRIC samples stratified by CSMD1 expression status by four different breast cancer PAM50 subtypes. Top 25% samples with high CSMD1 is in red, showing a better survival outcome compared to those with low CSMD1 expression. (PDF 8335 kb).

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.