Critically ill infants and children with rare diseases need equitable access to rapid and accurate diagnosis to direct clinical management. Over 2 years, the Acute Care Genomics program provided whole-genome sequencing to 290 families whose critically ill infants and children were admitted to hospitals throughout Australia with suspected genetic conditions. The average time to result was 2.9 d and diagnostic yield was 47%. We performed additional bioinformatic analyses and transcriptome sequencing in all patients who remained undiagnosed. Long-read sequencing and functional assays, ranging from clinically accredited enzyme analysis to bespoke quantitative proteomics, were deployed in selected cases. This resulted in an additional 19 diagnoses and an overall diagnostic yield of 54%. Diagnostic variants ranged from structural chromosomal abnormalities through to an intronic retrotransposon, disrupting splicing. Critical care management changed in 120 diagnosed patients (77%). This included major impacts, such as informing precision treatments, surgical and transplant decisions and palliation, in 94 patients (60%). Our results provide preliminary evidence of the clinical utility of integrating multi-omic approaches into mainstream diagnostic practice to fully realize the potential of rare disease genomic testing in a timely manner.

Chromobox protein homolog 7 (CBX7), a member of the Polycomb repressor complex, is a potent epigenetic regulator and gene silencer. Our group has previously reported that CBX7 functions as a tumor suppressor in ovarian cancer cells and its loss accelerated formation of carcinomatosis and drove tumor progression in an ovarian cancer mouse model. The goal of this study is to identify specific signaling pathways in the ovarian tumor microenvironment that down-regulate CBX7. Given that adipocytes are an integral component of the peritoneal cavity and the ovarian tumor microenvironment, we hypothesize that the adipose microenvironment is an important regulator of CBX7 expression.

Lipopolysaccharide (LPS) phase variation is a critical aspect of virulence in many Gram-negative bacteria. It is of particular importance to Coxiella burnetii, the biothreat pathogen that causes Q fever, as in vitro propagation of this organism leads to LPS truncation, which is associated with an attenuated and exempted from select agent status (Nine Mile II, NMII). Here, we demonstrate that NMII was recovered from the spleens of infected guinea pigs. Moreover, these strains exhibit a previously unrecognized form of elongated LPS and display increased virulence in comparison with the initial NMII strain. The reversion of a 3-bp mutation in the gene cbu0533 directly leads to LPS elongation. To address potential safety concerns, we introduce a modified NMII strain unable to produce elongated LPS.

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.

Disorganization of the valve extracellular matrix (ECM) is a hallmark of calcific aortic valve disease (CAVD). However, while microarchitectural features of the ECM can strongly influence the biological and mechanical behavior of tissues, little is known about the ECM microarchitecture in CAVD. In this work, we apply advanced imaging techniques to quantify spatially heterogeneous changes in collagen microarchitecture in CAVD. Human aortic valves were obtained from individuals between 50 and 75 years old with no evidence of valvular disease (healthy) and individuals who underwent valve replacement surgery due to severe stenosis (diseased). Second Harmonic Generation microscopy and subsequent image quantification revealed layer-specific changes in fiber characteristics in healthy and diseased valves. Specifically, the majority of collagen fiber changes in CAVD were found to occur in the spongiosa, where collagen fiber number increased by over 2-fold, and fiber width and density also significantly increased. Relatively few fibrillar changes occurred in the fibrosa in CAVD, where fibers became significantly shorter, but did not otherwise change in terms of number, width, density, or alignment. Immunohistochemical staining for lysyl oxidase showed localized increased expression in the diseased fibrosa. These findings reveal a more complex picture of valvular collagen enrichment and arrangement in CAVD than has previously been described using traditional analysis methods. Changes in fiber architecture may play a role in regulating the pathobiological events and mechanical properties of valves during CAVD. Additionally, characterization of the ECM microarchitecture can inform the design of fibrous scaffolds for heart valve tissue engineering.

Heparan sulfate (HS) is a polysaccharide fundamentally important for biologically activities. T/Tn antigens are universal carbohydrate cancer markers. Here, we report the specific imaging of these carbohydrates using a mesenchymal stem cell line and human umbilical vein endothelial cells (HUVEC). The staining specificities were demonstrated by comparing imaging of different glycans and validated by either removal of target glycans, which results in loss of signal, or installation of target glycans, which results in gain of signal. As controls, representative key glycans including O-GlcNAc, lactosaminyl glycans and hyaluronan were also imaged. HS staining revealed novel architectural features of the extracellular matrix (ECM) of HUVEC cells. Results from T/Tn antigen staining suggest that O-GalNAcylation is a rate-limiting step for O-glycan synthesis. Overall, these highly specific approaches for HS and T/Tn antigen imaging should greatly facilitate the detection and functional characterization of these biologically important glycans.

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the β2 subunit of the α4β2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10-8). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10-6) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that β2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction.

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.

Oesophageal adenocarcinoma (EAC) incidence is rapidly increasing in Western countries. A better understanding of EAC underpins efforts to improve early detection and treatment outcomes. While large EAC exome sequencing efforts to date have found recurrent loss-of-function mutations, oncogenic driving events have been underrepresented. Here we use a combination of whole-genome sequencing (WGS) and single-nucleotide polymorphism-array profiling to show that genomic catastrophes are frequent in EAC, with almost a third (32%, n=40/123) undergoing chromothriptic events. WGS of 22 EAC cases show that catastrophes may lead to oncogene amplification through chromothripsis-derived double-minute chromosome formation (MYC and MDM2) or breakage-fusion-bridge (KRAS, MDM2 and RFC3). Telomere shortening is more prominent in EACs bearing localized complex rearrangements. Mutational signature analysis also confirms that extreme genomic instability in EAC can be driven by somatic BRCA2 mutations. These findings suggest that genomic catastrophes have a significant role in the malignant transformation of EAC.

Tendons are ineffective at repairing sub-rupture fatigue injuries. Accordingly, we evaluated whether an exercise protocol that we have previously found to decrease structural damage kinks in fatigue damaged tendons, leads to improvement in mechanical properties. We hypothesized that exercise that promotes repair of fatigue damage will decrease apoptosis and increase the population of myofibroblasts. Rat patellar tendons underwent in vivo fatigue loading for 500 or 7200 cycles. Animals resumed cage activity for 2-weeks, then either remained cage active or began treadmill running until sacrifice at 4- or 10-weeks post-fatigue loading. Exercise following fatigue damage increased the stiffness back towards naïve levels, decreased apoptosis and increased the population of myofibroblasts. Next, proteins associated with inhibition of apoptosis (Collagen VI) or activation of myofibroblast (pSmad 2/3, fibrillin, integrin subunits αV and α5) were evaluated. Data suggests that collagen VI may not be integral to inhibition of apoptosis in this context. Exercise increased pSmad 2/3 and fibrillin in the insertion region for the 7200-cycles group. In addition, exercise decreased integrin αV and increased integrin α5 in fatigue damaged tendons. Data suggests that a decrease in apoptosis and an increase in population of myofibroblasts may be integral to remodeling of fatigue damaged tendons.

The cerebellum is a brain region responsible for motor coordination and for refining motor programs. While a great deal is known about the structure and connectivity of the mammalian cerebellum, fundamental questions regarding its function in behavior remain unanswered. Recently, the zebrafish has emerged as a useful model organism for cerebellar studies, owing in part to the similarity in cerebellar circuits between zebrafish and mammals. While the cell types composing their cerebellar cortical circuits are generally conserved with mammals, zebrafish lack deep cerebellar nuclei, and instead a majority of cerebellar output comes from a single type of neuron: the eurydendroid cell. To describe spatial patterns of cerebellar output in zebrafish, we have used genetic techniques to label and trace eurydendroid cells individually and en masse. We have found that cerebellar output targets the thalamus and optic tectum, and have confirmed the presence of pre-synaptic terminals from eurydendroid cells in these structures using a synaptically targeted GFP. By observing individual eurydendroid cells, we have shown that different medial-lateral regions of the cerebellum have eurydendroid cells projecting to different targets. Finally, we found topographic organization in the connectivity between the cerebellum and the optic tectum, where more medial eurydendroid cells project to the rostral tectum while lateral cells project to the caudal tectum. These findings indicate that there is spatial logic underpinning cerebellar output in zebrafish with likely implications for cerebellar function.

Chromobox protein homolog 7 (CBX7), a member of the Polycomb repressor complex, is a potent epigenetic regulator and gene silencer. Our group has previously reported that CBX7 functions as a tumor suppressor in ovarian cancer cells and its loss accelerated formation of carcinomatosis and drove tumor progression in an ovarian cancer mouse model. The goal of this study is to identify specific signaling pathways in the ovarian tumor microenvironment that down-regulate CBX7. Given that adipocytes are an integral component of the peritoneal cavity and the ovarian tumor microenvironment, we hypothesize that the adipose microenvironment is an important regulator of CBX7 expression.

Tumour cellularity, the relative proportion of tumour and normal cells in a sample, affects the sensitivity of mutation detection, copy number analysis, cancer gene expression and methylation profiling. Tumour cellularity is traditionally estimated by pathological review of sectioned specimens; however this method is both subjective and prone to error due to heterogeneity within lesions and cellularity differences between the sample viewed during pathological review and tissue used for research purposes. In this paper we describe a statistical model to estimate tumour cellularity from SNP array profiles of paired tumour and normal samples using shifts in SNP allele frequency at regions of loss of heterozygosity (LOH) in the tumour. We also provide qpure, a software implementation of the method. Our experiments showed that there is a medium correlation 0.42 ([Formula: see text]-value=0.0001) between tumor cellularity estimated by qpure and pathology review. Interestingly there is a high correlation 0.87 ([Formula: see text]-value [Formula: see text] 2.2e-16) between cellularity estimates by qpure and deep Ion Torrent sequencing of known somatic KRAS mutations; and a weaker correlation 0.32 ([Formula: see text]-value=0.004) between IonTorrent sequencing and pathology review. This suggests that qpure may be a more accurate predictor of tumour cellularity than pathology review. qpure can be downloaded from https://sourceforge.net/projects/qpure/.

Somatic mutation calling from next-generation sequencing data remains a challenge due to the difficulties of distinguishing true somatic events from artifacts arising from PCR, sequencing errors or mis-mapping. Tumor cellularity or purity, sub-clonality and copy number changes also confound the identification of true somatic events against a background of germline variants. We have developed a heuristic strategy and software (http://www.qcmg.org/bioinformatics/qsnp/) for somatic mutation calling in samples with low tumor content and we show the superior sensitivity and precision of our approach using a previously sequenced cell line, a series of tumor/normal admixtures, and 3,253 putative somatic SNVs verified on an orthogonal platform.

Variants of microRNAs (miRNAs), called isomiRs, are commonly reported in deep-sequencing studies; however, the functional significance of these variants remains controversial. Observational studies show that isomiR patterns are non-random, hinting that these molecules could be regulated and therefore functional, although no conclusive biological role has been demonstrated for these molecules.

Proteins of the mammalian PYHIN (IFI200/HIN-200) family are involved in defence against infection through recognition of foreign DNA. The family member absent in melanoma 2 (AIM2) binds cytosolic DNA via its HIN domain and initiates inflammasome formation via its pyrin domain. AIM2 lies within a cluster of related genes, many of which are uncharacterised in mouse. To better understand the evolution, orthology and function of these genes, we have documented the range of PYHIN genes present in representative mammalian species, and undertaken phylogenetic and expression analyses.

The return of research results (RoR) remains a complex and well-debated issue. Despite the debate, actual data related to the experience of giving individual results back, and the impact these results may have on clinical care and health outcomes, is sorely lacking. Through the work of the Australian Pancreatic Cancer Genome Initiative (APGI) we: (1) delineate the pathway back to the patient where actionable research data were identified; and (2) report the clinical utilisation of individual results returned. Using this experience, we discuss barriers and opportunities associated with a comprehensive process of RoR in large-scale genomic research that may be useful for others developing their own policies.

Our understanding of the importance of transcriptional regulation for biological function is continuously improving. We still know, however, comparatively little about how environmentally induced stress affects gene expression in vertebrates, and the consistency of transcriptional stress responses to different types of environmental stress. In this study, we used a multi-stressor approach to identify components of a common stress response as well as components unique to different types of environmental stress. We exposed individuals of the coral reef fish Pomacentrus moluccensis to hypoxic, hyposmotic, cold and heat shock and measured the responses of approximately 16,000 genes in liver. We also compared winter and summer responses to heat shock to examine the capacity for such responses to vary with acclimation to different ambient temperatures.

We describe a sibling pair displaying an early infantile-onset, progressive neurodegenerative phenotype, with symptoms of developmental delay and epileptic encephalopathy developing from 12 to 14 months of age. Using whole exome sequencing, compound heterozygous variants were identified in SLC5A6, which encodes the sodium-dependent multivitamin transporter (SMVT) protein. SMVT is an important transporter of the B-group vitamins biotin, pantothenate, and lipoate. The protein is ubiquitously expressed and has major roles in vitamin uptake in the digestive system, as well as transport of these vitamins across the blood-brain barrier. Pathogenicity of the identified variants was demonstrated by impaired biotin uptake of mutant SMVT. Identification of this vitamin transporter as the genetic basis of this disorder guided targeted therapeutic intervention, resulting clinically in improvement of the patient's neurocognitive and neuromotor function. This is the second report of biallelic mutations in SLC5A6 leading to a neurodegenerative disorder due to impaired biotin, pantothenate and lipoate uptake. The genetic and phenotypic overlap of these cases confirms mutations in SLC5A6 as the genetic cause of this disease phenotype. Recognition of the genetic disorder caused by SLC5A6 mutations is essential for early diagnosis and to facilitate timely intervention by triple vitamin (biotin, pantothenate, and lipoate) replacement therapy.

Non-caloric artificial sweeteners (NCAS) are widely used as a substitute for dietary sugars to control body weight or glycemia. Paradoxically, some interventional studies in humans and rodents have shown unfavorable changes in glucose homeostasis in response to NCAS consumption. The causative mechanisms are largely unknown, but adverse changes in gut microbiota have been proposed to mediate these effects. These findings have raised concerns about NCAS safety and called into question their broad use, but further physiological and dietary considerations must be first addressed before these results are generalized. We also reasoned that, since NCAS are bona fide ligands for sweet taste receptors (STRs) expressed in the intestine, some metabolic effects associated with NCAS use could be attributed to a common mechanism involving the host.