Molecular genetic diagnosis of kidney ciliopathies: Lessons from interpreting genomic sequencing data and the requirement for accurate phenotypic data.

INTRODUCTION: Massively parallel sequencing (MPS) techniques have made a major impact on the identification of the genetic basis of inherited kidney diseases such as the ciliopathy autosomal dominant polycystic kidney disease (ADPKD). Great care must be taken when analysing MPS data in isolation from accurate phenotypic information, as this can cause misdiagnosis. METHODS: Here, we describe a family trio, recruited to the Genomics England 100,000 Genomes Project, labelled as having cystic kidney disease, who were genetically unsolved following routine data analysis pipelines. We performed a bespoke reanalysis of Whole Genome Sequencing (WGS) data and coupled this with revised phenotypic data and targeted PCR and Sanger sequencing to provide a precise molecular genetic diagnosis. RESULTS: We detected a heterozygous PKD1 frameshift variant within the WGS data which segregated with the redefined ADPKD phenotypes. An additional heterozygous exon deletion in ALG8 was also found in affected and unaffected individuals, but its precise clinical significance remains unclear. CONCLUSION: This case illustrates that reanalysis of WGS data in unsolved cases of cystic kidney disease is valuable. Clinical phenotypes must be reassessed as these may have been incorrectly recorded and evolve over time. Undertaking additional studies including genotype-phenotype correlation in wider family members provides useful diagnostic information.

Use of whole genome sequencing to determine the genetic basis of visceral myopathies including Prune Belly syndrome.

Objectives/aims: The visceral myopathies (VM) are a group of disorders characterised by poorly contractile or acontractile smooth muscle. They manifest in both the GI and GU tracts, ranging from megacystis to Prune Belly syndrome. We aimed to apply a bespoke virtual genetic panel and describe novel variants associated with this condition using whole genome sequencing data within the Genomics England 100,000 Genomes Project. Methods: We screened the Genomics England 100,000 Genomes Project rare diseases database for patients with VM-related phenotypes. These patients were screened for sequence variants and copy number variants (CNV) in ACTG2, ACTA2, MYH11, MYLK, LMOD1, CHRM3, MYL9, FLNA and KNCMA1 by analysing whole genome sequencing data. The identified variants were analysed using variant effect predictor online tool, and any possible segregation in other family members and novel missense mutations was modelled using in silico tools. The VM cohort was also used to perform a genome-wide variant burden test in order to identify confirm gene associations in this cohort. Results: We identified 76 patients with phenotypes consistent with a diagnosis of VM. The range of presentations included megacystis/microcolon hypoperistalsis syndrome, Prune Belly syndrome and chronic intestinal pseudo-obstruction. Of the patients in whom we identified heterozygous ACTG2 variants, 7 had likely pathogenic variants including 1 novel likely pathogenic allele. There were 4 patients in whom we identified a heterozygous MYH11 variant of uncertain significance which leads to a frameshift and a predicted protein elongation. We identified one family in whom we found a heterozygous variant of uncertain significance in KCNMA1 which in silico models predicted to be disease causing and may explain the VM phenotype seen. We did not find any CNV changes in known genes leading to VM-related disease phenotypes. In this phenotype selected cohort, ACTG2 is the largest monogenic cause of VM-related disease accounting for 9% of the cohort, supported by a variant burden test approach, which identified ACTG2 variants as the largest contributor to VM-related phenotypes. Conclusions: VM are a group of disorders that are not easily classified and may be given different diagnostic labels depending on their phenotype. Molecular genetic analysis of these patients is valuable as it allows precise diagnosis and aids understanding of the underlying disease manifestations. We identified ACTG2 as the most frequent genetic cause of VM. We recommend a nomenclature change to 'autosomal dominant ACTG2 visceral myopathy' for patients with pathogenic variants in ACTG2 and associated VM phenotypes. Supplementary Information: The online version contains supplementary material available at 10.1007/s44162-023-00012-z.

Monoallelic intragenic POU3F2 variants lead to neurodevelopmental delay and hyperphagic obesity, confirming the gene's candidacy in 6q16.1 deletions.

While common obesity accounts for an increasing global health burden, its monogenic forms have taught us underlying mechanisms via more than 20 single-gene disorders. Among these, the most common mechanism is central nervous system dysregulation of food intake and satiety, often accompanied by neurodevelopmental delay (NDD) and autism spectrum disorder. In a family with syndromic obesity, we identified a monoallelic truncating variant in POU3F2 (alias BRN2) encoding a neural transcription factor, which has previously been suggested as a driver of obesity and NDD in individuals with the 6q16.1 deletion. In an international collaboration, we identified ultra-rare truncating and missense variants in another ten individuals sharing autism spectrum disorder, NDD, and adolescent-onset obesity. Affected individuals presented with low-to-normal birth weight and infantile feeding difficulties but developed insulin resistance and hyperphagia during childhood. Except for a variant leading to early truncation of the protein, identified variants showed adequate nuclear translocation but overall disturbed DNA-binding ability and promotor activation. In a cohort with common non-syndromic obesity, we independently observed a negative correlation of POU3F2 gene expression with BMI, suggesting a role beyond monogenic obesity. In summary, we propose deleterious intragenic variants of POU3F2 to cause transcriptional dysregulation associated with hyperphagic obesity of adolescent onset with variable NDD.

Biallelic variants in CEP164 cause a motile ciliopathy-like syndrome.

Ciliopathies may be classed as primary or motile depending on the underlying ciliary defect and are usually considered distinct clinical entities. Primary ciliopathies are associated with multisystem syndromes typically affecting the brain, kidney, and eye, as well as other organ systems such as the liver, skeleton, auditory system, and metabolism. Motile ciliopathies are a heterogenous group of disorders with defects in specialised motile ciliated tissues found within the lung, brain, and reproductive system, and are associated with primary ciliary dyskinesia, bronchiectasis, infertility and rarely hydrocephalus. Primary and motile cilia share defined core ultra-structures with an overlapping proteome, and human disease phenotypes can reflect both primary and motile ciliopathies. CEP164 encodes a centrosomal distal appendage protein vital for primary ciliogenesis. Human CEP164 mutations are typically described in patients with nephronophthisis-related primary ciliopathies but have also been implicated in motile ciliary dysfunction. Here we describe a patient with an atypical motile ciliopathy phenotype and biallelic CEP164 variants. This work provides further evidence that CEP164 mutations can contribute to both primary and motile ciliopathy syndromes, supporting their functional and clinical overlap, and informs the investigation and management of CEP164 ciliopathy patients.

Rare disease gene association discovery from burden analysis of the 100,000 Genomes Project data.

To discover rare disease-gene associations, we developed a gene burden analytical framework and applied it to rare, protein-coding variants from whole genome sequencing of 35,008 cases with rare diseases and their family members recruited to the 100,000 Genomes Project (100KGP). Following in silico triaging of the results, 88 novel associations were identified including 38 with existing experimental evidence. We have published the confirmation of one of these associations, hereditary ataxia with UCHL1 , and independent confirmatory evidence has recently been published for four more. We highlight a further seven compelling associations: hypertrophic cardiomyopathy with DYSF and SLC4A3 where both genes show high/specific heart expression and existing associations to skeletal dystrophies or short QT syndrome respectively; monogenic diabetes with UNC13A with a known role in the regulation of ß cells and a mouse model with impaired glucose tolerance; epilepsy with KCNQ1 where a mouse model shows seizures and the existing long QT syndrome association may be linked; early onset Parkinson's disease with RYR1 with existing links to tremor pathophysiology and a mouse model with neurological phenotypes; anterior segment ocular abnormalities associated with POMK showing expression in corneal cells and with a zebrafish model with developmental ocular abnormalities; and cystic kidney disease with COL4A3 showing high renal expression and prior evidence for a digenic or modifying role in renal disease. Confirmation of all 88 associations would lead to potential diagnoses in 456 molecularly undiagnosed cases within the 100KGP, as well as other rare disease patients worldwide, highlighting the clinical impact of a large-scale statistical approach to rare disease gene discovery.

The Idiopathic Pulmonary Fibrosis-Associated Single Nucleotide Polymorphism RS35705950 Is Transcribed in a MUC5B Promoter Associated Long Non-Coding RNA (AC061979.1).

LncRNAs are involved in regulatory processes in the human genome, including gene expression. The rs35705950 SNP, previously associated with IPF, overlaps with the recently annotated lncRNA AC061979.1, a 1712 nucleotide transcript located within the MUC5B promoter at chromosome 11p15.5. To document the expression pattern of the transcript, we processed 3.9 TBases of publicly available RNA-SEQ data across 27 independent studies involving lung airway epithelial cells. Epithelial lung cells showed expression of this putative pancRNA. The findings were independently validated in cell lines and primary cells. The rs35705950 is found within a conserved region (from fish to primates) within the expressed sequence indicating functional importance. These results implicate the rs35705950-containing AC061979.1 pancRNA as a novel component of the MUC5B expression control minicircuitry.

Monoallelic IFT140 pathogenic variants are an important cause of the autosomal dominant polycystic kidney-spectrum phenotype.

Autosomal dominant polycystic kidney disease (ADPKD), characterized by progressive cyst formation/expansion, results in enlarged kidneys and often end stage kidney disease. ADPKD is genetically heterogeneous; PKD1 and PKD2 are the common loci (∼78% and ∼15% of families) and GANAB, DNAJB11, and ALG9 are minor genes. PKD is a ciliary-associated disease, a ciliopathy, and many syndromic ciliopathies have a PKD phenotype. In a multi-cohort/-site collaboration, we screened ADPKD-diagnosed families that were naive to genetic testing (n = 834) or for whom no PKD1 and PKD2 pathogenic variants had been identified (n = 381) with a PKD targeted next-generation sequencing panel (tNGS; n = 1,186) or whole-exome sequencing (WES; n = 29). We identified monoallelic IFT140 loss-of-function (LoF) variants in 12 multiplex families and 26 singletons (1.9% of naive families). IFT140 is a core component of the intraflagellar transport-complex A, responsible for retrograde ciliary trafficking and ciliary entry of membrane proteins; bi-allelic IFT140 variants cause the syndromic ciliopathy, short-rib thoracic dysplasia (SRTD9). The distinctive monoallelic phenotype is mild PKD with large cysts, limited kidney insufficiency, and few liver cysts. Analyses of the cystic kidney disease probands of Genomics England 100K showed that 2.1% had IFT140 LoF variants. Analysis of the UK Biobank cystic kidney disease group showed probands with IFT140 LoF variants as the third most common group, after PKD1 and PKD2. The proximity of IFT140 to PKD1 (∼0.5 Mb) in 16p13.3 can cause diagnostic confusion, and PKD1 variants could modify the IFT140 phenotype. Importantly, our studies link a ciliary structural protein to the ADPKD spectrum.