Rapid Paediatric Sequencing (RaPS): comprehensive real-life workflow for rapid diagnosis of critically ill children.

BACKGROUND: Rare genetic conditions are frequent risk factors for, or direct causes of, paediatric intensive care unit (PICU) admission. Such conditions are frequently suspected but unidentified at PICU admission. Compassionate and effective care is greatly assisted by definitive diagnostic information. There is therefore a need to provide a rapid genetic diagnosis to inform clinical management.To date, whole genome sequencing (WGS) approaches have proved successful in diagnosing a proportion of children with rare diseases, but results may take months to report. Our aim was to develop an end-to-end workflow for the use of rapid WGS for diagnosis in critically ill children in a UK National Health Service (NHS) diagnostic setting. METHODS: We sought to establish a multidisciplinary Rapid Paediatric Sequencing team for case selection, trio WGS, rapid bioinformatics sequence analysis and a phased analysis and reporting system to prioritise genes with a high likelihood of being causal. RESULTS: Trio WGS in 24 critically ill children led to a molecular diagnosis in 10 (42%) through the identification of causative genetic variants. In 3 of these 10 individuals (30%), the diagnostic result had an immediate impact on the individual's clinical management. For the last 14 trios, the shortest time taken to reach a provisional diagnosis was 4 days (median 8.5 days). CONCLUSION: Rapid WGS can be used to diagnose and inform management of critically ill children within the constraints of an NHS clinical diagnostic setting. We provide a robust workflow that will inform and facilitate the rollout of rapid genome sequencing in the NHS and other healthcare systems globally.

Distinctive characters of Nostoc genomes in cyanolichens.

BACKGROUND: Cyanobacteria of the genus Nostoc are capable of forming symbioses with a wide range of organism, including a diverse assemblage of cyanolichens. Only certain lineages of Nostoc appear to be able to form a close, stable symbiosis, raising the question whether symbiotic competence is determined by specific sets of genes and functionalities. RESULTS: We present the complete genome sequencing, annotation and analysis of two lichen Nostoc strains. Comparison with other Nostoc genomes allowed identification of genes potentially involved in symbioses with a broad range of partners including lichen mycobionts. The presence of additional genes necessary for symbiotic competence is likely reflected in larger genome sizes of symbiotic Nostoc strains. Some of the identified genes are presumably involved in the initial recognition and establishment of the symbiotic association, while others may confer advantage to cyanobionts during cohabitation with a mycobiont in the lichen symbiosis. CONCLUSIONS: Our study presents the first genome sequencing and genome-scale analysis of lichen-associated Nostoc strains. These data provide insight into the molecular nature of the cyanolichen symbiosis and pinpoint candidate genes for further studies aimed at deciphering the genetic mechanisms behind the symbiotic competence of Nostoc. Since many phylogenetic studies have shown that Nostoc is a polyphyletic group that includes several lineages, this work also provides an improved molecular basis for demarcation of a Nostoc clade with symbiotic competence.

Metagenomic natural product discovery in lichen provides evidence for a family of biosynthetic pathways in diverse symbioses.

Bacteria are a major source of natural products that provide rich opportunities for both chemical and biological investigation. Although the vast majority of known bacterial metabolites derive from free-living organisms, increasing evidence supports the widespread existence of chemically prolific bacteria living in symbioses. A strategy based on bioinformatic prediction, symbiont cultivation, isotopic enrichment, and advanced analytics was used to characterize a unique polyketide, nosperin, from a lichen-associated Nostoc sp. cyanobacterium. The biosynthetic gene cluster and the structure of nosperin, determined from 30 µg of compound, are related to those of the pederin group previously known only from nonphotosynthetic bacteria associated with beetles and marine sponges. The presence of this natural product family in such highly dissimilar associations suggests that some bacterial metabolites may be specific to symbioses with eukaryotes and encourages exploration of other symbioses for drug discovery and better understanding of ecological interactions mediated by complex bacterial metabolites.

Unravelling morphoea aetiopathogenesis by next-generation sequencing of paired skin biopsies.

BACKGROUND: Morphoea can have a significant disease burden. Aetiopathogenesis remains poorly understood, with very limited existing genetic studies. Linear morphoea (LM) may follow Blascho's lines of epidermal development, providing potential pathogenic clues. OBJECTIVE: The first objective of this study was to identify the presence of primary somatic epidermal mosaicism in LM. The second objective was tTo explore differential gene expression in morphoea epidermis and dermis to identify potential pathogenic molecular pathways and tissue layer cross-talk. METHODOLOGY: Skin biopsies from paired affected and contralateral unaffected skin were taken from 16 patients with LM. Epidermis and dermis were isolated using a 2-step chemical-physical separation protocol. Whole Genome Sequencing (WGS; n = 4 epidermal) and RNA-seq (n = 5-epidermal, n = 5-dermal) with gene expression analysis via GSEA-MSigDBv6.3 and PANTHER-v14.1 pathway analyses, were performed. RTqPCR and immunohistochemistry were used to replicate key results. RESULTS: Sixteen participants (93.8% female, mean age 27.7 yrs disease-onset) were included. Epidermal WGS identified no single affected gene or SNV. However, many potential disease-relevant pathogenic variants were present, including ADAMTSL1 and ADAMTS16. A highly proliferative, inflammatory and profibrotic epidermis was seen, with significantly-overexpressed TNFα-via-NFkB, TGFß, IL6/JAKSTAT and IFN-signaling, apoptosis, p53 and KRAS-responses. Upregulated IFI27 and downregulated LAMA4 potentially represent initiating epidermal 'damage' signals and enhanced epidermal-dermal communication. Morphoea dermis exhibited significant profibrotic, B-cell and IFN-signatures, and upregulated morphogenic patterning pathways such as Wnt. CONCLUSION: This study supports the absence of somatic epidermal mosaicism in LM, and identifies potential disease-driving epidermal mechanisms, epidermal-dermal interactions and disease-specific dermal differential-gene-expression in morphoea. We propose a potential molecular narrative for morphoea aetiopathogenesis which could help guide future targeted studies and therapies.

Neuroinflammation, autoinflammation, splenomegaly and anemia caused by bi-allelic mutations in IRAK4.

We describe a novel, severe autoinflammatory syndrome characterized by neuroinflammation, systemic autoinflammation, splenomegaly, and anemia (NASA) caused by bi-allelic mutations in IRAK4. IRAK-4 is a serine/threonine kinase with a pivotal role in innate immune signaling from toll-like receptors and production of pro-inflammatory cytokines. In humans, bi-allelic mutations in IRAK4 result in IRAK-4 deficiency and increased susceptibility to pyogenic bacterial infections, but autoinflammation has never been described. We describe 5 affected patients from 2 unrelated families with compound heterozygous mutations in IRAK4 (c.C877T (p.Q293*)/c.G958T (p.D320Y); and c.A86C (p.Q29P)/c.161 + 1G>A) resulting in severe systemic autoinflammation, massive splenomegaly and severe transfusion dependent anemia and, in 3/5 cases, severe neuroinflammation and seizures. IRAK-4 protein expression was reduced in peripheral blood mononuclear cells (PBMC) in affected patients. Immunological analysis demonstrated elevated serum tumor necrosis factor (TNF), interleukin (IL) 1 beta (IL-1ß), IL-6, IL-8, interferon α2a (IFN-α2a), and interferon ß (IFN-ß); and elevated cerebrospinal fluid (CSF) IL-6 without elevation of CSF IFN-α despite perturbed interferon gene signature. Mutations were located within the death domain (DD; p.Q29P and splice site mutation c.161 + 1G>A) and kinase domain (p.Q293*/p.D320Y) of IRAK-4. Structure-based modeling of the DD mutation p.Q29P showed alteration in the alignment of a loop within the DD with loss of contact distance and hydrogen bond interactions with IRAK-1/2 within the myddosome complex. The kinase domain mutation p.D320Y was predicted to stabilize interactions within the kinase active site. While precise mechanisms of autoinflammation in NASA remain uncertain, we speculate that loss of negative regulation of IRAK-4 and IRAK-1; dysregulation of myddosome assembly and disassembly; or kinase active site instability may drive dysregulated IL-6 and TNF production. Blockade of IL-6 resulted in immediate and complete amelioration of systemic autoinflammation and anemia in all 5 patients treated; however, neuroinflammation has, so far proven recalcitrant to IL-6 blockade and the janus kinase (JAK) inhibitor baricitinib, likely due to lack of central nervous system penetration of both drugs. We therefore highlight that bi-allelic mutation in IRAK4 may be associated with a severe and complex autoinflammatory and neuroinflammatory phenotype that we have called NASA (neuroinflammation, autoinflammation, splenomegaly and anemia), in addition to immunodeficiency in humans.

ZSWIM7 Is Associated With Human Female Meiosis and Familial Primary Ovarian Insufficiency.

BACKGROUND: Primary ovarian insufficiency (POI) affects 1% of women and is associated with significant medical consequences. A genetic cause for POI can be found in up to 30% of women, elucidating key roles for these genes in human ovary development. OBJECTIVE: We aimed to identify the genetic mechanism underlying early-onset POI in 2 sisters from a consanguineous pedigree. METHODS: Genome sequencing and variant filtering using an autosomal recessive model was performed in the 2 affected sisters and their unaffected family members. Quantitative reverse transcriptase PCR (qRT-PCR) and RNA sequencing were used to study the expression of key genes at critical stages of human fetal gonad development (Carnegie Stage 22/23, 9 weeks post conception (wpc), 11 wpc, 15/16 wpc, 19/20 wpc) and in adult tissue. RESULTS: Only 1 homozygous variant cosegregating with the POI phenotype was found: a single nucleotide substitution in zinc finger SWIM-type containing 7 (ZSWIM7), NM_001042697.2: c.173C > G; resulting in predicted loss-of-function p.(Ser58*). qRT-PCR demonstrated higher expression of ZSWIM7 in the 15/16 wpc ovary compared with testis, corresponding to peak meiosis in the fetal ovary. RNA sequencing of fetal gonad samples showed that ZSWIM7 has a similar temporal expression profile in the developing ovary to other homologous recombination genes. MAIN CONCLUSIONS: Disruption of ZSWIM7 is associated with POI in humans. ZSWIM7 is likely to be important for human homologous recombination; these findings expand the range of genes associated with POI in women.

Pathogenic variants in the human m6A reader YTHDC2 are associated with primary ovarian insufficiency.

Primary ovarian insufficiency (POI) affects 1% of women and carries significant medical and psychosocial sequelae. Approximately 10% of POI has a defined genetic cause, with most implicated genes relating to biological processes involved in early fetal ovary development and function. Recently, Ythdc2, an RNA helicase and N6-methyladenosine reader, has emerged as a regulator of meiosis in mice. Here, we describe homozygous pathogenic variants in YTHDC2 in 3 women with early-onset POI from 2 families: c. 2567C>G, p.P856R in the helicase-associated (HA2) domain and c.1129G>T, p.E377*. We demonstrated that YTHDC2 is expressed in the developing human fetal ovary and is upregulated in meiotic germ cells, together with related meiosis-associated factors. The p.P856R variant resulted in a less flexible protein that likely disrupted downstream conformational kinetics of the HA2 domain, whereas the p.E377* variant truncated the helicase core. Taken together, our results reveal that YTHDC2 is a key regulator of meiosis in humans and pathogenic variants within this gene are associated with POI.

Activating mutations in BRAF disrupt the hypothalamo-pituitary axis leading to hypopituitarism in mice and humans.

Germline mutations in BRAF and other components of the MAPK pathway are associated with the congenital syndromes collectively known as RASopathies. Here, we report the association of Septo-Optic Dysplasia (SOD) including hypopituitarism and Cardio-Facio-Cutaneous (CFC) syndrome in patients harbouring mutations in BRAF. Phosphoproteomic analyses demonstrate that these genetic variants are gain-of-function mutations leading to activation of the MAPK pathway. Activation of the MAPK pathway by conditional expression of the BrafV600E/+ allele, or the knock-in BrafQ241R/+ allele (corresponding to the most frequent human CFC-causing mutation, BRAF p.Q257R), leads to abnormal cell lineage determination and terminal differentiation of hormone-producing cells, causing hypopituitarism. Expression of the BrafV600E/+ allele in embryonic pituitary progenitors leads to an increased expression of cell cycle inhibitors, cell growth arrest and apoptosis, but not tumour formation. Our findings show a critical role of BRAF in hypothalamo-pituitary-axis development both in mouse and human and implicate mutations found in RASopathies as a cause of endocrine deficiencies in humans.

Novel missense variants in the RNF213 gene from a European family with Moyamoya disease.

In this report, we present a European family with six individuals affected with Moyamoya disease (MMD). We detected two novel missense variants in the Moyamoya susceptibility gene RNF213, c.12553A>G (p.(Lys4185Glu)) and c.12562G>A (p.(Ala4188Thr)). Cosegregation of the variants with MMD, as well as a previous report of a variant affecting the same amino acid residue in unrelated MMD patients, supports the role of RNF213 in the pathogenesis of MMD.

Sotos Syndrome Presenting with Neonatal Hyperinsulinaemic Hypoglycaemia, Extensive Thrombosis, and Multisystem Involvement.

Initially described as an uncommon presenting feature of Sotos syndrome (SoS), over the last decades, congenital hyperinsulinaemic hypoglycaemia (CHI) has been increasingly reported in association with this condition. The mechanism responsible for CHI in SoS is unclear. We report the case of a neonate presenting with CHI and extensive venous and arterial thrombosis associated with kidney, heart, liver, skeleton, and brain abnormalities and finally diagnosed with SoS on whole genome sequencing. Our case describes an extended phenotype associated with SoS presenting with CHI (including thrombosis and liver dysfunction) and reinforces the association of transient CHI with SoS. The case also shows that an early neonatal diagnosis of rare genetic conditions is challenging, especially in acutely unwell patients, and that in complex cases with incomplete, atypical, or overlapping phenotypes, broad genomic testing by whole exome or whole genome sequencing may be a useful diagnostic strategy.

Cloning and heterologous transcription of a polyketide synthase gene from the lichen Solorina crocea.

Lichens and most ascomycete fungi produce polyketide secondary metabolites often with valuable biological activities. Their biosynthesis is primarily governed by large iterative multifunctional type I polyketide synthases. Although there has been good progress studying filamentous non-lichenized fungi, there is limited information on polyketide biosynthesis in lichens and their mycobionts, due to their slow growth, difficulties in establishing pure cultures, and the absence of methods for direct genetic manipulation. However, heterologous expression in a surrogate host offers an alternative approach for exploring lichen polyketide biosynthesis. Here, we report cloning of a type I polyketide synthase gene from the foliose lichen Solorina crocea and its heterologous transcription in the filamentous fungus Aspergillus oryzae, including processing of the transcript. No new polyketide product was detected. The lichen polyketide synthase showed greatest homology with uncharacterized genes from filamentous fungi and lower homology with proteins catalysing biosynthesis of the decaketide alternapyrone and the tetraketide side-chain of squalestatin. The technology platform utilized here presents a useful tool for functional characterization of fungal biosynthetic genes and provides a means for novel production of valuable compounds.