Tracing genetic diversity captures the molecular basis of misfolding disease

Genetic variation in human populations can result in the misfolding and aggregation of proteins, giving rise to systemic and neurodegenerative diseases that require management by proteostasis. Here, we define the role of GRP94, the endoplasmic reticulum Hsp90 chaperone paralog, in managing alpha-1-antitrypsin deficiency on a residue-by-residue basis using Gaussian process regression-based machine learning to profile the spatial covariance relationships that dictate protein folding arising from sequence variants in the population. Covariance analysis suggests a role for the ATPase activity of GRP94 in controlling the N- to C-terminal cooperative folding of alpha-1-antitrypsin responsible for the correction of liver aggregation and lung-disease phenotypes of alpha-1-antitrypsin deficiency. Gaussian process-based spatial covariance profiling provides a standard model built on covariant principles to evaluate the role of proteostasis components in guiding information flow from genome to proteome in response to genetic variation, potentially allowing us to intervene in the onset and progression of complex multi-system human diseases.

The prevalence and clinical features of MYO7A-related hearing loss including DFNA11, DFNB2 and USH1B

The MYO7A gene is known to be responsible for both syndromic hearing loss (Usher syndrome type1B:USH1B) and non-syndromic hearing loss including autosomal dominant and autosomal recessive inheritance (DFNA11, DFNB2). However, the prevalence and detailed clinical features of MYO7A-associated hearing loss across a large population remain unclear. In this study, we conducted next-generation sequencing analysis for a large cohort of 10,042 Japanese hearing loss patients. As a result, 137 patients were identified with MYO7A-associated hearing loss so that the prevalence among Japanese hearing loss patients was 1.36%. We identified 70 disease-causing candidate variants in this study, with 36 of them being novel variants. All variants identified in autosomal dominant cases were missense or in-frame deletion variants. Among the autosomal recessive cases, all patients had at least one missense variant. On the other hand, in patients with Usher syndrome, almost half of the patients carried biallelic null variants (nonsense, splicing, and frameshift variants). Most of the autosomal dominant cases showed late-onset progressive hearing loss. On the other hand, cases with autosomal recessive inheritance or Usher syndrome showed congenital or early-onset hearing loss. The visual symptoms in the Usher syndrome cases developed between age 5-15, and the condition was diagnosed at about 6-15 years of age.

Expanding the clinical spectrum of biglycan-related Meester-Loeys syndrome

Pathogenic loss-of-function variants in BGN, an X-linked gene encoding biglycan, are associated with Meester-Loeys syndrome (MRLS), a thoracic aortic aneurysm/dissection syndrome. Since the initial publication of five probands in 2017, we have considerably expanded our MRLS cohort to a total of 18 probands (16 males and 2 females). Segregation analyses identified 36 additional BGN variant-harboring family members (9 males and 27 females). The identified BGN variants were shown to lead to loss-of-function by cDNA and Western Blot analyses of skin fibroblasts or were strongly predicted to lead to loss-of-function based on the nature of the variant. No (likely) pathogenic missense variants without additional (predicted) splice effects were identified. Interestingly, a male proband with a deletion spanning the coding sequence of BGN and the 5' untranslated region of the downstream gene (ATP2B3) presented with a more severe skeletal phenotype. This may possibly be explained by expressional activation of the downstream ATPase ATP2B3 (normally repressed in skin fibroblasts) driven by the remnant BGN promotor. This study highlights that aneurysms and dissections in MRLS extend beyond the thoracic aorta, affecting the entire arterial tree, and cardiovascular symptoms may coincide with non-specific connective tissue features. Furthermore, the clinical presentation is more severe and penetrant in males compared to females. Extensive analysis at RNA, cDNA, and/or protein level is recommended to prove a loss-of-function effect before determining the pathogenicity of identified BGN missense and non-canonical splice variants. In conclusion, distinct mechanisms may underlie the wide phenotypic spectrum of MRLS patients carrying loss-of-function variants in BGN.

The aldehyde dehydrogenase 2 rs671 variant enhances amyloid β pathology

In the ALDH2 rs671 variant, a guanine changes to an adenine, resulting in a dramatic decrease in the catalytic activity of the enzyme. Population-based data are contradictory about whether this variant increases the risk of Alzheimer's disease. In East Asian populations, the prevalence of the ALDH2 rs671 variant is 30-50%, making the National Human Brain Bank for Development and Function (the largest brain bank in East Asia) an important resource to explore the link between the ALDH2 rs671 polymorphism and Alzheimer's disease pathology. Here, using 469 postmortem brains, we find that while the ALDH2 rs671 variant is associated with increased plaque deposits and a higher Aβ40/42 ratio, it is not an independent risk factor for Alzheimer's disease. Mechanistically, we show that lower ALDH2 activity leads to 4-HNE accumulation in the brain. The (R)-4-HNE enantiomer adducts to residue Lys53 of C99, favoring Aβ40 generation in the Golgi apparatus. Decreased ALDH2 activity also lowers inflammatory factor secretion, as well as amyloid β phagocytosis and spread in brains of patients with Alzheimer's disease. We thus define the relationship between the ALDH2 rs671 polymorphism and amyloid β pathology, and find that ALDH2 rs671 is a key regulator of Aβ40 or Aβ42 generation.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

DAG1 haploinsufficiency is associated with sporadic and familial isolated or pauci-symptomatic hyperCKemia

DAG1 encodes for dystroglycan, a key component of the dystrophin-glycoprotein complex (DGC) with a pivotal role in skeletal muscle function and maintenance. Biallelic loss-of-function DAG1 variants cause severe muscular dystrophy and muscle-eye-brain disease. A possible contribution of DAG1 deficiency to milder muscular phenotypes has been suggested. We investigated the genetic background of twelve subjects with persistent mild-to-severe hyperCKemia to dissect the role of DAG1 in this condition. Genetic testing was performed through exome sequencing (ES) or custom NGS panels including various genes involved in a spectrum of muscular disorders. Histopathological and Western blot analyses were performed on muscle biopsy samples obtained from three patients. We identified seven novel heterozygous truncating variants in DAG1 segregating with isolated or pauci-symptomatic hyperCKemia in all families. The variants were rare and predicted to lead to nonsense-mediated mRNA decay or the formation of a truncated transcript. In four cases, DAG1 variants were inherited from similarly affected parents. Histopathological analysis revealed a decreased expression of dystroglycan subunits and Western blot confirmed a significantly reduced expression of beta-dystroglycan in muscle samples. This study supports the pathogenic role of DAG1 haploinsufficiency in isolated or pauci-symptomatic hyperCKemia, with implications for clinical management and genetic counseling.

Identification of a DLG3 stop mutation in the MRX20 family

Here, we identified the causal mutation in the MRX20 family, one of the larger X-linked pedigrees that have been described in which no gene had been identified up till now. In 1995, the putative disease gene had been mapped to the pericentromeric region on the X chromosome, but no follow-up studies were performed. Here, whole exome sequencing (WES) on two affected and one unaffected family member revealed the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) that segregated with the affected individuals in the family. DLG3 mutations have been consequently associated with intellectual disability and are a plausible explanation for the clinical abnormalities observed in this family. In addition, we identified two other variants co-segregating with the phenotype: a stop gain mutation in SSX1 (c.358G>T/p.(Glu120Ter)) (NM_001278691.2) and a nonsynonymous SNV in USP27X (c.56 A>G/p.(Gln19Arg)) (NM_001145073.3). RNA sequencing revealed 14 differentially expressed genes (p value < 0.1) in 7 affected males compared to 4 unaffected males of the family, including four genes known to be associated with neurological disorders. Thus, in this paper we identified the c.195del/p.(Thr66ProfsTer55) mutation in the DLG3 gene (NM_021120.4) as likely responsible for the phenotype observed in the MRX20 family.

A deletion variant in LMX1B causing nail-patella syndrome in Japanese twins

Nail-patella syndrome (NPS) is a hereditary disease caused by pathogenic variants in LMX1B and characterized by nail, limb, and renal symptoms. This study revealed a likely pathogenic LMX1B variant, NM_002316.4: c.723_726delinsC (p.Ser242del), in Japanese twins with clubfoot. The patients' mother, who shared this variant, developed proteinuria after delivery. p.Ser242del is located in the homeodomain of the protein, in which variants that cause renal disease tend to cluster. Our findings highlight p.Ser242del as a likely pathogenic variant, expanding our knowledge of NPS.

The frequency of pathogenic variation in the All of Us cohort reveals ancestry-driven disparities

Disparities in data underlying clinical genomic interpretation is an acknowledged problem, but there is a paucity of data demonstrating it. The All of Us Research Program is collecting data including whole-genome sequences, health records, and surveys for at least a million participants with diverse ancestry and access to healthcare, representing one of the largest biomedical research repositories of its kind. Here, we examine pathogenic and likely pathogenic variants that were identified in the All of Us cohort. The European ancestry subgroup showed the highest overall rate of pathogenic variation, with 2.26% of participants having a pathogenic variant. Other ancestry groups had lower rates of pathogenic variation, including 1.62% for the African ancestry group and 1.32% in the Latino/Admixed American ancestry group. Pathogenic variants were most frequently observed in genes related to Breast/Ovarian Cancer or Hypercholesterolemia. Variant frequencies in many genes were consistent with the data from the public gnomAD database, with some notable exceptions resolved using gnomAD subsets. Differences in pathogenic variant frequency observed between ancestral groups generally indicate biases of ascertainment of knowledge about those variants, but some deviations may be indicative of differences in disease prevalence. This work will allow targeted precision medicine efforts at revealed disparities.

Unraveling the genetic architecture of congenital vertebral malformation with reference to the developing spine

Congenital vertebral malformation, affecting 0.13-0.50 per 1000 live births, has an immense locus heterogeneity and complex genetic architecture. In this study, we analyze exome/genome sequencing data from 873 probands with congenital vertebral malformation and 3794 control individuals. Clinical interpretation identifies Mendelian etiologies in 12.0% of the probands and reveals a muscle-related disease mechanism. Gene-based burden test of ultra-rare variants identifies risk genes with large effect sizes (ITPR2, TBX6, TPO, H6PD, and SEC24B). To further investigate the biological relevance of the genetic association signals, we perform single-nucleus RNAseq on human embryonic spines. The burden test signals are enriched in the notochord at early developmental stages and myoblast/myocytes at late stages, highlighting their critical roles in the developing spine. Our work provides insights into the developmental biology of the human spine and the pathogenesis of spine malformation.

Clinically accessible amplitude-based multiplex ddPCR assay for tryptase genotyping

Hereditary α tryptasemia (HαT) is an autosomal dominant trait characterized by increased TPSAB1 copy number (CN) encoding α-tryptase. The determination of HαT is being discussed as an important biomarker to be included in risk assessment models and future diagnostic algorithms for patients with mastocytosis and anaphylaxis. Due to the complex genetic structure at the human tryptase locus, genetic testing for tryptase gene composition is presently notably limited and infrequently pursued. This study aimed to develop, optimise and validate a multiplex droplet digital PCR (ddPCR) assay that can reliably quantify α- and β-tryptase encoding sequences in a single reaction. To optimise the ddPCR conditions and establish an amplitude-based multiplex ddPCR assay, additional primers and probes, a thermal gradient with varying annealing temperatures, different primers/probe concentrations, and various initial DNA quantities were tested. Results obtained from all 114 samples analysed using multiplex ddPCR were identical to those obtained through the use of original duplex assays. Utilizing this multiplex ddPCR assay, in contrast to conducting distinct duplex ddPCRs, presents noteworthy benefits for tryptase genotyping. These advantages encompass a substantial threefold decrease in material costs and considerable time savings. Consequently, this approach exhibits high suitability and particularly captures interest for routine clinical implementation.

Episodic ataxia type 2 with a novel missense variant (Leu602Arg) in CACNA1A

Autosomal dominant episodic ataxia type 2 (EA2) is caused by variants in CACNA1A. We examined a 20-year-old male with EA symptoms from a Japanese family with hereditary EA. Cerebellar atrophy was not evident, but single photon emission computed tomography showed cerebellar hypoperfusion. We identified a novel nonsynonymous variant in CACNA1A, NM_001127222.2:c.1805T>G (p.Leu602Arg), which is predicted to be functionally deleterious; therefore, this variant is likely responsible for EA2 in this pedigree.

Large scale plasma proteomics identifies novel proteins and protein networks associated with heart failure development

Heart failure (HF) causes substantial morbidity and mortality but its pathobiology is incompletely understood. The proteome is a promising intermediate phenotype for discovery of novel mechanisms. We measured 4877 plasma proteins in 13,900 HF-free individuals across three analysis sets with diverse age, geography, and HF ascertainment to identify circulating proteins and protein networks associated with HF development. Parallel analyses in Atherosclerosis Risk in Communities study participants in mid-life and late-life and in Trøndelag Health Study participants identified 37 proteins consistently associated with incident HF independent of traditional risk factors. Mendelian randomization supported causal effects of 10 on HF, HF risk factors, or left ventricular size and function, including matricellular (e.g. SPON1, MFAP4), senescence-associated (FSTL3, IGFBP7), and inflammatory (SVEP1, CCL15, ITIH3) proteins. Protein co-regulation network analyses identified 5 modules associated with HF risk, two of which were influenced by genetic variants that implicated trans hotspots within the VTN and CFH genes.

Possibilities and limitations of antisense oligonucleotide therapies for the treatment of monogenic disorders

Antisense oligonucleotides (ASOs) are incredibly versatile molecules that can be designed to specifically target and modify RNA transcripts to slow down or halt rare genetic disease progression. They offer the potential to target groups of patients or can be tailored for individual cases. Nonetheless, not all genetic variants and disorders are amenable to ASO-based treatments, and hence, it is important to consider several factors before embarking on the drug development journey. Here, we discuss which genetic disorders have the potential to benefit from a specific type of ASO approach, based on the pathophysiology of the disease and pathogenic variant type, as well as those disorders that might not be suitable for ASO therapies. We further explore additional aspects, such as the target tissues, intervention time points, and potential clinical benefits, which need to be considered before developing a compound. Overall, we provide an overview of the current potentials and limitations of ASO-based therapeutics for the treatment of monogenic disorders.

Fatty acid synthesis suppresses dietary polyunsaturated fatty acid use

Dietary polyunsaturated fatty acids (PUFA) are increasingly recognized for their health benefits, whereas a high production of endogenous fatty acids - a process called de novo lipogenesis (DNL) - is closely linked to metabolic diseases. Determinants of PUFA incorporation into complex lipids are insufficiently understood and may influence the onset and progression of metabolic diseases. Here we show that fatty acid synthase (FASN), the key enzyme of DNL, critically determines the use of dietary PUFA in mice and humans. Moreover, the combination of FASN inhibition and PUFA-supplementation decreases liver triacylglycerols (TAG) in mice fed with high-fat diet. Mechanistically, FASN inhibition causes higher PUFA uptake via the lysophosphatidylcholine transporter MFSD2A, and a diacylglycerol O-acyltransferase 2 (DGAT2)-dependent incorporation of PUFA into TAG. Overall, the outcome of PUFA supplementation may depend on the degree of endogenous DNL and combining PUFA supplementation and FASN inhibition might be a promising approach to target metabolic disease.

Bi-allelic truncating variants in CASP2 underlie a neurodevelopmental disorder with lissencephaly

Lissencephaly (LIS) is a malformation of cortical development due to deficient neuronal migration and abnormal formation of cerebral convolutions or gyri. Thirty-one LIS-associated genes have been previously described. Recently, biallelic pathogenic variants in CRADD and PIDD1, have associated with LIS impacting the previously established role of the PIDDosome in activating caspase-2. In this report, we describe biallelic truncating variants in CASP2, another subunit of PIDDosome complex. Seven patients from five independent families presenting with a neurodevelopmental phenotype were identified through GeneMatcher-facilitated international collaborations. Exome sequencing analysis was carried out and revealed two distinct novel homozygous (NM_032982.4:c.1156delT (p.Tyr386ThrfsTer25), and c.1174 C > T (p.Gln392Ter)) and compound heterozygous variants (c.[130 C > T];[876 + 1 G > T] p.[Arg44Ter];[?]) in CASP2 segregating within the families in a manner compatible with an autosomal recessive pattern. RNA studies of the c.876 + 1 G > T variant indicated usage of two cryptic splice donor sites, each introducing a premature stop codon. All patients from whom brain MRIs were available had a typical fronto-temporal LIS and pachygyria, remarkably resembling the CRADD and PIDD1-related neuroimaging findings. Other findings included developmental delay, attention deficit hyperactivity disorder, hypotonia, seizure, poor social skills, and autistic traits. In summary, we present patients with CASP2-related ID, anterior-predominant LIS, and pachygyria similar to previously reported patients with CRADD and PIDD1-related disorders, expanding the genetic spectrum of LIS and lending support that each component of the PIDDosome complex is critical for normal development of the human cerebral cortex and brain function.

A population-based survey of FBN1 variants in Iceland reveals underdiagnosis of Marfan syndrome

Marfan syndrome (MFS) is an autosomal dominant condition characterized by aortic aneurysm, skeletal abnormalities, and lens dislocation, and is caused by variants in the FBN1 gene. To explore causes of MFS and the prevalence of the disease in Iceland we collected information from all living individuals with a clinical diagnosis of MFS in Iceland (n = 32) and performed whole-genome sequencing of those who did not have a confirmed genetic diagnosis (27/32). Moreover, to assess a potential underdiagnosis of MFS in Iceland we attempted a genotype-based approach to identify individuals with MFS. We interrogated deCODE genetics' database of 35,712 whole-genome sequenced individuals to search for rare sequence variants in FBN1. Overall, we identified 15 pathogenic or likely pathogenic variants in FBN1 in 44 individuals, only 22 of whom were previously diagnosed with MFS. The most common of these variants, NM_000138.4:c.8038 C > T p.(Arg2680Cys), is present in a multi-generational pedigree, and was found to stem from a single forefather born around 1840. The p.(Arg2680Cys) variant associates with a form of MFS that seems to have an enrichment of abdominal aortic aneurysm, suggesting that this may be a particularly common feature of p.(Arg2680Cys)-associated MFS. Based on these combined genetic and clinical data, we show that MFS prevalence in Iceland could be as high as 1/6,600 in Iceland, compared to 1/10,000 based on clinical diagnosis alone, which indicates underdiagnosis of this actionable genetic disorder.

Biallelic ATP2B1 variants as a likely cause of a novel neurodevelopmental malformation syndrome with primary hypoparathyroidism

ATP2B1 encodes plasma membrane calcium-transporting-ATPase1 and plays an essential role in maintaining intracellular calcium homeostasis that regulates diverse signaling pathways. Heterozygous de novo missense and truncating ATP2B1 variants are associated with a neurodevelopmental phenotype of variable expressivity. We describe a proband with distinctive craniofacial gestalt, Pierre-Robin sequence, neurodevelopmental and growth deficit, periventricular heterotopia, brachymesophalangy, cutaneous syndactyly, and persistent hypocalcemia from primary hypoparathyroidism. Proband-parent trio exome sequencing identified compound heterozygous ATP2B1 variants: a maternally inherited splice-site (c.3060+2 T > G) and paternally inherited missense c.2938 G > T; p.(Val980Leu). Reverse-transcription-PCR on the proband's fibroblast-derived mRNA showed aberrantly spliced ATP2B1 transcripts targeted for nonsense-mediated decay. All correctly-spliced ATP2B1 mRNA encoding p.(Val980Leu) functionally causes decreased cellular Ca2+ extrusion. Immunoblotting showed reduced fibroblast ATP2B1. We conclude that biallelic ATP2B1 variants are the likely cause of the proband's phenotype, strengthening the association of ATP2B1 as a neurodevelopmental gene and expanding the phenotypic characterization of a biallelic loss-of-function genotype.

vcfdist: accurately benchmarking phased small variant calls in human genomes

Accurately benchmarking small variant calling accuracy is critical for the continued improvement of human whole genome sequencing. In this work, we show that current variant calling evaluations are biased towards certain variant representations and may misrepresent the relative performance of different variant calling pipelines. We propose solutions, first exploring the affine gap parameter design space for complex variant representation and suggesting a standard. Next, we present our tool vcfdist and demonstrate the importance of enforcing local phasing for evaluation accuracy. We then introduce the notion of partial credit for mostly-correct calls and present an algorithm for clustering dependent variants. Lastly, we motivate using alignment distance metrics to supplement precision-recall curves for understanding variant calling performance. We evaluate the performance of 64 phased Truth Challenge V2 submissions and show that vcfdist improves measured insertion and deletion performance consistency across variant representations from R2 = 0.97243 for baseline vcfeval to 0.99996 for vcfdist.

Deep phenotypic characterization of the retinal dystrophy in patients with RNU4ATAC-associated Roifman syndrome

PURPOSE

To characterize the retinal phenotype in RNU4ATAC-associated Roifman syndrome.

METHODS

Ten patients (including 8 males) with molecularly confirmed Roifman syndrome underwent detailed ophthalmologic evaluation including fundus imaging, fundus autofluorescence (FAF) imaging, spectral-domain optical coherence tomography (SD-OCT), and electroretinography (ERG). Six patients had follow-up eye exams. All patients also underwent comprehensive examination for features of extra-retinal Roifman syndrome.

RESULTS

All patients had biallelic RNU4ATAC variants. Nyctalopia was common (7/10). Visual acuity at presentation ranged from 20/20 to 20/200 (Age Range: 5-41 years). Retinal exam revealed features of generalized retinopathy with mid-peripheral pigment epithelial changes. A para or peri-foveal ring of hyper-autofluorescence was the commonest FAF abnormality noted (6/8). The SD-OCT demonstrated relative preservation of the foveal ellipsoid zone in six cases; associated features included cystoid changes (5/10) and posterior staphyloma (3/10). The ERG was abnormal in all patients; nine showed generalized rod-cone dystrophy, whilst one patient with sectoral retinal involvement only had isolated rod dystrophy (20 years old). On follow-up examination (Mean duration: 8.16 years), progressive loss of visual acuity (2/6), mid-peripheral retinal atrophy (3/6) or shortening of ellipsoid zone width (1/6) were observed.

CONCLUSION

This study has characterized the retinal phenotype in RNU4ATAC-associated Roifman syndrome. Retinal involvement is universal, early-onset, and overall, the retinal and FAF features are consistent with rod-cone degeneration that is slowly progressive over time. The sub-foveal retinal ultrastructure is relatively preserved in majority of patients. Phenotypic variability independent of age exists, and more study of allelic- and sex-based determinants of disease severity are necessary.

Börjeson-Forssman-Lehmann syndrome: delineating the clinical and allelic spectrum in 14 new families

Börjeson-Forssman-Lehmann syndrome (BFLS) is an X-linked intellectual disability syndrome caused by variants in the PHF6 gene. We ascertained 19 individuals from 15 families with likely pathogenic or pathogenic PHF6 variants (11 males and 8 females). One family had previously been reported. Six variants were novel. We analysed the clinical and genetic findings in our series and compared them with reported BFLS patients. Affected males had classic features of BFLS including intellectual disability, distinctive facies, large ears, gynaecomastia, hypogonadism and truncal obesity. Carrier female relatives of affected males were unaffected or had only mild symptoms. The phenotype of affected females with de novo variants overlapped with the males but included linear skin hyperpigmentation and a higher frequency of dental, retinal and cortical brain anomalies. Complications observed in our series included keloid scarring, digital fibromas, absent vaginal orifice, neuropathy, umbilical hernias, and talipes. Our analysis highlighted sex-specific differences in PHF6 variant types and locations. Affected males often have missense variants or small in-frame deletions while affected females tend to have truncating variants or large deletions/duplications. Missense variants were found in a minority of affected females and clustered in the highly constrained PHD2 domain of PHF6. We propose recommendations for the evaluation and management of BFLS patients. These results further delineate and extend the genetic and phenotypic spectrum of BFLS.

A biallelic variant in COX18 cause isolated Complex IV deficiency associated with neonatal encephalo-cardio-myopathy and axonal sensory neuropathy

Pathogenic variants impacting upon assembly of mitochondrial respiratory chain Complex IV (Cytochrome c Oxidase or COX) predominantly result in early onset mitochondrial disorders often leading to CNS, skeletal and cardiac muscle manifestations. The aim of this study is to describe a molecular defect in the COX assembly factor gene COX18 as the likely cause of a neonatal form of mitochondrial encephalo-cardio-myopathy and axonal sensory neuropathy. The proband is a 19-months old female displaying hypertrophic cardiomyopathy at birth and myopathy with axonal sensory neuropathy and failure to thrive developing in the first months of life. Serum lactate was consistently increased. Whole exome sequencing allowed the prioritization of the unreported homozygous substitution NM_001297732.2:c.667 G > C p.(Asp223His) in COX18. Patient's muscle biopsy revealed severe and diffuse COX deficiency and striking mitochondrial abnormalities. Biochemical and enzymatic studies in patient's myoblasts and in HEK293 cells after COX18 silencing showed a severe impairment of both COX activity and assembly. The biochemical defect was partially rescued by delivery of wild-type COX18 cDNA into patient's myoblasts. Our study identifies a novel defect of COX assembly and expands the number of nuclear genes involved in a mitochondrial disorder due to isolated COX deficiency.

SLC35D3 promotes white adipose tissue browning to ameliorate obesity by NOTCH signaling

White adipose tissue browning can promote lipid burning to increase energy expenditure and improve adiposity. Here, we show that Slc35d3 expression is significantly lower in adipose tissues of obese mice. While adipocyte-specific Slc35d3 knockin is protected against diet-induced obesity, adipocyte-specific Slc35d3 knockout inhibits white adipose tissue browning and causes decreased energy expenditure and impaired insulin sensitivity in mice. Mechanistically, we confirm that SLC35D3 interacts with the NOTCH1 extracellular domain, which leads to the accumulation of NOTCH1 in the endoplasmic reticulum and thus inhibits the NOTCH1 signaling pathway. In addition, knockdown of Notch1 in mouse inguinal white adipose tissue mediated by orthotopic injection of AAV8-adiponectin-shNotch1 shows considerable improvement in obesity and glucolipid metabolism, which is more pronounced in adipocyte-specific Slc35d3 knockout mice than in knockin mice. Overall, in this study, we reveal that SLC35D3 is involved in obesity via NOTCH1 signaling, and low adipose SLC35D3 expression in obesity might be a therapeutic target for obesity and associated metabolic disorders.

Clinical, genetic and structural delineation of RPL13-related spondyloepimetaphyseal dysplasia suggest extra-ribosomal functions of eL13

Spondyloepimetaphyseal dysplasia with severe short stature, RPL13-related (SEMD-RPL13), MIM#618728), is a rare autosomal dominant disorder characterized by short stature and skeletal changes such as mild spondylar and epimetaphyseal dysplasia affecting primarily the lower limbs. The genetic cause was first reported in 2019 by Le Caignec et al., and six disease-causing variants in the gene coding for a ribosomal protein, RPL13 (NM_000977.3) have been identified to date. This study presents clinical and radiographic data from 12 affected individuals aged 2-64 years from seven unrelated families, showing highly variable manifestations. The affected individuals showed a range from mild to severe short stature, retaining the same radiographic pattern of spondylar- and epi-metaphyseal dysplasia, but with varying severity of the hip and knee deformities. Two new missense variants, c.548 G>A, p.(Arg183His) and c.569 G>T, p.(Arg190Leu), and a previously known splice variant c.477+1G>A were identified, confirming mutational clustering in a highly specific RNA binding motif. Structural analysis and interpretation of the variants' impact on the protein suggests that disruption of extra-ribosomal functions of the protein through binding of mRNA may play a role in the skeletal phenotype of SEMD-RPL13. In addition, we present gonadal and somatic mosaicism for the condition.

Recessive dystrophic epidermolysis bullosa caused by a novel COL7A1 variant with isodisomy

Recessive dystrophic epidermolysis bullosa is a genetic collagen disorder characterized by skin fragility that leads to generalized severe blistering, wounds, and scarring. In this report, we present a patient with a novel COL7A1 homozygous nonsense variant, c.793C>T p.(Gln265*). Although the parents were not consanguineous, both were heterozygous carriers of the variant. Single nucleotide polymorphism (SNP) array analysis revealed an isodisomy area on 3p22.1p21.1, encompassing COL7A1, suggesting that the variant originated from a common ancestor.

Allelic hierarchy for USH2A influences auditory and visual phenotypes in South Korean patients

When medical genetic syndromes are influenced by allelic hierarchies, mutant alleles have distinct effects on clinical phenotypes. Genotype-phenotype correlations for Usher syndrome type 2 (USH2) suggest that the USH2A gene exhibits an allelic hierarchy. Here, we analyzed the phenotypes and genotypes of 16 South Korean patients with USH2A biallelic variants to investigate an allelic hierarchy from audiological and ophthalmological perspectives. Using whole exome and genome sequencing, 18 mutant alleles, including 4 novel alleles, were identified and implicated in USH2A-related disorders. Truncated alleles were linked to earlier onset of subjective hearing loss and more severe thresholds; biallelic truncated alleles had more severe effects. Truncated alleles were also associated with retinal structure degeneration and severe functional deterioration. However, younger patients (aged < 16 years) did not exhibit overt retinitis pigmentosa even when they had biallelic truncated alleles, suggesting that USH2A-related USH2 can mimic nonsyndromic hearing loss. For truncated alleles, there was a clear correlation between mean hearing threshold and 30-Hz flicker electroretinography implicit time. This study provides the first evidence of an USH2A-related allelic hierarchy among South Korean patients; our data yield valuable insights concerning the natural courses of clinical phenotypes and how genotype-based therapies may be used.

KiT-GENIE, the French genetic biobank of kidney transplantation

KiT-GENIE is a monocentric DNA biobank set up to consolidate the very rich and homogeneous DIVAT French cohort of kidney donors and recipients (D/R) in order to explore the molecular factors involved in kidney transplantation outcomes. We collected DNA samples for kidney transplantations performed in Nantes, and we leveraged GWAS genotyping data for securing high-quality genetic data with deep SNP and HLA annotations through imputations and for inferring D/R genetic ancestry. Overall, the biobank included 4217 individuals (n = 1945 D + 2,272 R, including 1969 D/R pairs), 7.4 M SNPs and over 200 clinical variables. KiT-GENIE represents an accurate snapshot of kidney transplantation clinical practice in Nantes between 2002 and 2018, with an enrichment in living kidney donors (17%) and recipients with focal segmental glomerulosclerosis (4%). Recipients were predominantly male (63%), of European ancestry (93%), with a mean age of 51yo and 86% experienced their first graft over the study period. D/R pairs were 93% from European ancestry, and 95% pairs exhibited at least one HLA allelic mismatch. The mean follow-up time was 6.7 years with a hindsight up to 25 years. Recipients experienced biopsy-proven rejection and graft loss for 16.6% and 21.3%, respectively. KiT-GENIE constitutes one of the largest kidney transplantation genetic cohorts worldwide to date. It includes homogeneous high-quality clinical and genetic data for donors and recipients, hence offering a unique opportunity to investigate immunogenetic and genetic factors, as well as donor-recipient interactions and mismatches involved in rejection, graft survival, primary disease recurrence and other comorbidities.

Epidemiology of the disorders of the Pik3ca-related overgrowth spectrum (Pros)

PIK3CA pathogenic variants are responsible for a group of overgrowth syndromes, collectively known as PIK3CA-Related Overgrowth Spectrum (PROS). These gain-of-function variants arise postzygotically, and, according to time of onset, kind of embryonal tissue affected and regional body extension, give rise to heterogeneous phenotypes. PROS rarity and heterogeneity hamper the correct estimation of its epidemiology. Our work represents the first attempt to define the prevalence of PROS according to the established diagnostic criteria and molecular analysis and based on solid demographic data. We assessed the prevalence in Piedmont Region (Italy), including in the study all participants diagnosed with PROS born there from 1998 to 2021. The search identified 37 cases of PROS born across the 25-year period, providing a prevalence of 1:22,313 live births. Molecular analysis was positive in 81.0% of participants. Taking into account the cases with a detected variant in PIK3CA (n = 30), prevalence of molecularly positive PROS was 1:27,519.

Skewed X-chromosome inactivation in unsolved neurodevelopmental disease cases can guide re-evaluation For X-linked genes

Despite major advances in genome technology and analysis, >50% of patients with a neurodevelopmental disorder (NDD) remain undiagnosed after extensive evaluation. A point in case is our clinically heterogeneous cohort of NDD patients that remained undiagnosed after FRAXA testing, chromosomal microarray analysis and trio exome sequencing (ES). In this study, we explored the frequency of non-random X chromosome inactivation (XCI) in the mothers of male patients and affected females, the rationale being that skewed XCI might be masking previously discarded genetic variants found on the X chromosome. A multiplex fluorescent PCR-based assay was used to analyse the pattern of XCI after digestion with HhaI methylation-sensitive restriction enzyme. In families with skewed XCI, we re-evaluated trio-based ES and identified pathogenic variants and a deletion on the X chromosome. Linkage analysis and RT-PCR were used to further study the inactive X chromosome allele, and Xdrop long-DNA technology was used to define chromosome deletion boundaries. We found skewed XCI (>90%) in 16/186 (8.6%) mothers of NDD males and in 12/90 (13.3%) NDD females, far beyond the expected rate of XCI in the normal population (3.6%, OR = 4.10; OR = 2.51). By re-analyzing ES and clinical data, we solved 7/28 cases (25%) with skewed XCI, identifying variants in KDM5C, PDZD4, PHF6, TAF1, OTUD5 and ZMYM3, and a deletion in ATRX. We conclude that XCI profiling is a simple assay that targets a subgroup of patients that can benefit from re-evaluation of X-linked variants, thus improving the diagnostic yield in NDD patients and identifying new X-linked disorders.

Systematic analysis of paralogous regions in 41,755 exomes uncovers clinically relevant variation

The short lengths of short-read sequencing reads challenge the analysis of paralogous genomic regions in exome and genome sequencing data. Most genetic variants within these homologous regions therefore remain unidentified in standard analyses. Here, we present a method (Chameleolyser) that accurately identifies single nucleotide variants and small insertions/deletions (SNVs/Indels), copy number variants and ectopic gene conversion events in duplicated genomic regions using whole-exome sequencing data. Application to a cohort of 41,755 exome samples yields 20,432 rare homozygous deletions and 2,529,791 rare SNVs/Indels, of which we show that 338,084 are due to gene conversion events. None of the SNVs/Indels are detectable using regular analysis techniques. Validation by high-fidelity long-read sequencing in 20 samples confirms >88% of called variants. Focusing on variation in known disease genes leads to a direct molecular diagnosis in 25 previously undiagnosed patients. Our method can readily be applied to existing exome data.

Clinical and genetic analysis further delineates the phenotypic spectrum of ALDH1A3-related anophthalmia and microphthalmia

Biallelic pathogenic variants in ALDH1A3 are responsible for approximately 11% of recessively inherited cases of severe developmental eye anomalies. Some individuals can display variable neurodevelopmental features, but the relationship to the ALDH1A3 variants remains unclear. Here, we describe seven unrelated families with biallelic pathogenic ALDH1A3 variants: four compound heterozygous and three homozygous. All affected individuals had bilateral anophthalmia/microphthalmia (A/M), three with additional intellectual or developmental delay, one with autism and seizures and three with facial dysmorphic features. This study confirms that individuals with biallelic pathogenic ALDH1A3 variants consistently manifest A/M, but additionally display neurodevelopmental features with significant intra- and interfamilial variability. Furthermore, we describe the first case with cataract and highlight the importance of screening ALDH1A3 variants in nonconsanguineous families with A/M.

FBXL4 mutations cause excessive mitophagy via BNIP3/BNIP3L accumulation leading to mitochondrial DNA depletion syndrome

Mitochondria are essential organelles found in eukaryotic cells that play a crucial role in ATP production through oxidative phosphorylation (OXPHOS). Mitochondrial DNA depletion syndrome (MTDPS) is a group of genetic disorders characterized by the reduction of mtDNA copy number, leading to deficiencies in OXPHOS and mitochondrial functions. Mutations in FBXL4, a substrate-binding adaptor of Cullin 1-RING ubiquitin ligase complex (CRL1), are associated with MTDPS, type 13 (MTDPS13). Here, we demonstrate that, FBXL4 directly interacts with the mitophagy cargo receptors BNIP3 and BNIP3L, promoting their degradation through the ubiquitin-proteasome pathway via the assembly of an active CRL1FBXL4 complex. However, MTDPS13-associated FBXL4 mutations impair the assembly of an active CRL1FBXL4 complex. This results in a notable accumulation of BNIP3/3L proteins and robust mitophagy even at basal levels. Excessive mitophagy was observed in Knockin (KI) mice carrying a patient-derived FBXL4 mutation and cortical neurons (CNs)-induced from MTDPS13 patient human induced pluripotent stem cells (hiPSCs). In summary, our findings suggest that abnormal activation of BNIP3/BNIP3L-dependent mitophagy impairs mitochondrial homeostasis and underlies FBXL4-mutated MTDPS13.

Double heterozygous pathogenic mutations in KIF3C and ZNF513 cause hereditary gingival fibromatosis

Hereditary gingival fibromatosis (HGF) is a rare inherited condition with fibromatoid hyperplasia of the gingival tissue that exhibits great genetic heterogeneity. Five distinct loci related to non-syndromic HGF have been identified; however, only two disease-causing genes, SOS1 and REST, inducing HGF have been identified at two loci, GINGF1 and GINGF5, respectively. Here, based on a family pedigree with 26 members, including nine patients with HGF, we identified double heterozygous pathogenic mutations in the ZNF513 (c.C748T, p.R250W) and KIF3C (c.G1229A, p.R410H) genes within the GINGF3 locus related to HGF. Functional studies demonstrated that the ZNF513 p.R250W and KIF3C p.R410H variants significantly increased the expression of ZNF513 and KIF3C in vitro and in vivo. ZNF513, a transcription factor, binds to KIF3C exon 1 and participates in the positive regulation of KIF3C expression in gingival fibroblasts. Furthermore, a knock-in mouse model confirmed that heterozygous or homozygous mutations within Zfp513 (p.R250W) or Kif3c (p.R412H) alone do not led to clear phenotypes with gingival fibromatosis, whereas the double mutations led to gingival hyperplasia phenotypes. In addition, we found that ZNF513 binds to the SOS1 promoter and plays an important positive role in regulating the expression of SOS1. Moreover, the KIF3C p.R410H mutation could activate the PI3K and KCNQ1 potassium channels. ZNF513 combined with KIF3C regulates gingival fibroblast proliferation, migration, and fibrosis response via the PI3K/AKT/mTOR and Ras/Raf/MEK/ERK pathways. In summary, these results demonstrate ZNF513 + KIF3C as an important genetic combination in HGF manifestation and suggest that ZNF513 mutation may be a major risk factor for HGF.

Mechanosensitive ion channel gene survey suggests potential roles in primary open angle glaucoma

Although glaucoma is a disease modulated by eye pressure, the mechanisms of pressure sensing in the eye are not well understood. Here, we investigated associations between mechanosensitive ion channel gene variants and primary open-angle glaucoma (POAG). Common (minor allele frequency > 5%) single nucleotide polymorphisms located within the genomic regions of 20 mechanosensitive ion channel genes in the K2P, TMEM63, PIEZO and TRP channel families were assessed using genotype data from the NEIGHBORHOOD consortium of 3853 cases and 33,480 controls. Rare (minor allele frequency < 1%) coding variants were assessed using exome array genotyping data for 2606 cases and 2606 controls. Association with POAG was analyzed using logistic regression adjusting for age and sex. Two rare PIEZO1 coding variants with protective effects were identified in the NEIGHBOR dataset: R1527H, (OR 0.17, P = 0.0018) and a variant that alters a canonical splice donor site, g.16-88737727-C-G Hg38 (OR 0.38, P = 0.02). Both variants showed similar effects in the UK Biobank and the R1527H also in the FinnGen database. Several common variants also reached study-specific thresholds for association in the NEIGHBORHOOD dataset. These results identify novel variants in several mechanosensitive channel genes that show associations with POAG, suggesting that these channels may be potential therapeutic targets.

Curated incidence of lysosomal storage diseases from the Taiwan Biobank

Lysosomal storage diseases (LSDs) are a group of metabolic disorders resulting from a deficiency in one of the lysosomal hydrolases. Most LSDs are inherited in an autosomal or X-linked recessive manner. As LSDs are rare, their true incidence in Taiwan remains unknown. In this study, we used high-coverage whole-genome sequencing data from 1,495 Taiwanese individuals obtained from the Taiwan Biobank. We found 3826 variants in 71 genes responsible for autosomal recessive LSDs. We first excluded benign variants by allele frequency and other criteria. As a result, 270 variants were considered disease-causing. We curated these variants using published guidelines from the American College of Medical Genetics and Genomics (ACMG). Our results revealed a combined incidence rate of 13 per 100,000 (conservative estimation by pathologic and likely pathogenic variants; 95% CI 6.92-22.23) to 94 per 100,000 (extended estimation by the inclusion of variants of unknown significance; 95% CI 75.96-115.03) among 71 autosomal recessive disease-associated genes. The conservative estimations were similar to those in published clinical data. No disease-causing mutations were found for 18 other diseases; thus, these diseases are likely extremely rare in Taiwan. The study results are important for designing screening and treatment methods for LSDs in Taiwan and demonstrate the importance of mutation curation to avoid overestimating disease incidences from genomic data.

Whole genome sequencing resolves 10 years diagnostic odyssey in familiar myxoma

Carney complex (CNC) is an ultrarare disorder causing cutaneous and cardiac myxomas, primary pigmented nodular adrenocortical disease, hypophyseal adenoma, and gonadal tumours. Genetic alterations are often missed under routine genetic testing. Pathogenic variants in PRKAR1A are identified in most cases, while large exonic or chromosomal deletions have only been reported in a few cases. Our aim was to identify the causal genetic alteration in our kindred with a clinical diagnosis of CNC and prove its pathogenic role by functional investigation. Targeted testing of PRKAR1A gene, whole exome and whole genome sequencing (WGS) were performed in the proband, one clinically affected and one unaffected relative. WGS identified a novel, large, 10,662 bp (10.6 kbp; LRG_514t1:c.-10403_-7 + 265del; hg19, chr17:g.66498293_66508954del) deletion in the promoter of PRKAR1A in heterozygous form in the affected family members. The exact breakpoints and the increased enzyme activity in deletion carriers compared to wild type carrier were proved. Segregation analysis and functional evaluation of PKA activity confirmed the pathogenic role of this alteration. A novel deletion upstream of the PRKAR1A gene was proved to be the cause of CNC. Our study underlines the need for WGS in molecular genetic testing of patients with monogenic disorders where conventional genetic analysis fails.

Somatic and germinal mosaicism in a Han Chinese family with laminopathies

"Laminopathies" refers to a wide spectrum of myopathies caused by mutations in the LMNA gene. These myopathies include limb girdle muscular dystrophy type 1B (LGMD1B) and dilated cardiomyopathy 1 A (DCM1A), which are both autosomal dominant neurogenetic diseases. There have been few studies on mosaicism in laminopathies. Herein, a Han Chinese family with laminopathies was enrolled in our study. Genetic analysis revealed that the proband carried a novel splice site mutation, c. 1158-3 C > T, in the LMNA gene due to her mother having de novo somatic and gonadal mosaicism. Reverse-transcription polymerase chain reaction (RT-PCR) analysis revealed reduced levels of LMNA mRNA in the proband, which were probably due to nonsense-mediated mRNA decay (NMD). Western blotting revealed reduced lamin A/C protein levels in the skeletal muscle tissue of the proband. In this family, the clinical phenotypes of the proband's mother were normal, and the c. 1158-3 C > T splicing mutation was identified in the blood sample of the proband's mother. Thus, the mutation could be easily considered to be nonpathogenic. Our study emphasizes the importance of mosaicism in the identification of pathogenic variants and genetic counseling.

BRAT1-related disorders: phenotypic spectrum and phenotype-genotype correlations from 97 patients

BRAT1 biallelic variants are associated with rigidity and multifocal seizure syndrome, lethal neonatal (RMFSL), and neurodevelopmental disorder associating cerebellar atrophy with or without seizures syndrome (NEDCAS). To date, forty individuals have been reported in the literature. We collected clinical and molecular data from 57 additional cases allowing us to study a large cohort of 97 individuals and draw phenotype-genotype correlations. Fifty-nine individuals presented with BRAT1-related RMFSL phenotype. Most of them had no psychomotor acquisition (100%), epilepsy (100%), microcephaly (91%), limb rigidity (93%), and died prematurely (93%). Thirty-eight individuals presented a non-lethal phenotype of BRAT1-related NEDCAS phenotype. Seventy-six percent of the patients in this group were able to walk and 68% were able to say at least a few words. Most of them had cerebellar ataxia (82%), axial hypotonia (79%) and cerebellar atrophy (100%). Genotype-phenotype correlations in our cohort revealed that biallelic nonsense, frameshift or inframe deletion/insertion variants result in the severe BRAT1-related RMFSL phenotype (46/46; 100%). In contrast, genotypes with at least one missense were more likely associated with NEDCAS (28/34; 82%). The phenotype of patients carrying splice variants was variable: 41% presented with RMFSL (7/17) and 59% with NEDCAS (10/17).

Myeloid/lymphoid neoplasms with eosinophilia and tyrosine kinase gene fusions: reevaluation of the defining characteristics in a registry-based cohort

In a registry-based analysis of 135 patients with "myeloid/lymphoid neoplasms with eosinophilia and tyrosine kinase gene fusions" (MLN-TK; FIP1L1::PDGFRA, n = 78; PDGFRB, diverse fusions, n = 26; FGFR1, diverse, n = 9; JAK2, diverse, n = 11; ETV6::ABL1, n = 11), we sought to evaluate the disease-defining characteristics. In 81/135 (60%) evaluable patients, hypereosinophilia (>1.5 × 109/l) was observed in 40/44 (91%) FIP1L1::PDGFRA and 7/7 (100%) ETV6::ABL1 positive patients but only in 13/30 (43%) patients with PDGFRB, FGFR1, and JAK2 fusion genes while 9/30 (30%) patients had no eosinophilia. Monocytosis >1 × 109/l was identified in 27/81 (33%) patients, most frequently in association with hypereosinophilia (23/27, 85%). Overall, a blast phase (BP) was diagnosed in 38/135 (28%) patients (myeloid, 61%; lymphoid, 39%), which was at extramedullary sites in 18 (47%) patients. The comparison between patients with PDGFRA/PDGFRB vs. FGFR1, JAK2, and ETV6::ABL1 fusion genes revealed a similar occurrence of primary BP (17/104, 16% vs. 8/31 26%, p = 0.32), a lower frequency (5/87, 6% vs. 8/23, 35%, p = 0.003) of and a later progression (median 87 vs. 19 months, p = 0.053) into secondary BP, and a better overall survival from diagnosis of BP (17.1 vs. 1.7 years, p < 0.0008). We conclude that hypereosinophilia with or without monocytosis and various phenotypes of BP occur at variable frequencies in MLN-TK.

A common allele of HLA is associated with asymptomatic SARS-CoV-2 infection

Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic1-4. Although most global efforts have focused on severe illness in COVID-19, examining asymptomatic infection provides a unique opportunity to consider early immunological features that promote rapid viral clearance. Here, postulating that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection, we enrolled 29,947 individuals, for whom high-resolution HLA genotyping data were available, in a smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n = 1,428) comprised unvaccinated individuals who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci with disease course and identified a strong association between HLA-B*15:01 and asymptomatic infection, observed in two independent cohorts. Suggesting that this genetic association is due to pre-existing T cell immunity, we show that T cells from pre-pandemic samples from individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF. The majority of the reactive T cells displayed a memory phenotype, were highly polyfunctional and were cross-reactive to a peptide derived from seasonal coronaviruses. The crystal structure of HLA-B*15:01-peptide complexes demonstrates that the peptides NQKLIANQF and NQKLIANAF (from OC43-CoV and HKU1-CoV) share a similar ability to be stabilized and presented by HLA-B*15:01. Finally, we show that the structural similarity of the peptides underpins T cell cross-reactivity of high-affinity public T cell receptors, providing the molecular basis for HLA-B*15:01-mediated pre-existing immunity.

High throughput human genotyping for variants associated with malarial disease outcomes using custom targeted amplicon sequencing

Malaria has exhibited the strongest known selective pressure on the human genome in recent history and is the evolutionary driving force behind genetic conditions, such as sickle-cell disease, glucose-6-phosphatase deficiency, and some other erythrocyte defects. Genomic studies (e.g., The 1000 Genomes project) have provided an invaluable baseline for human genetics, but with an estimated two thousand ethno-linguistic groups thought to exist across the African continent, our understanding of the genetic differences between indigenous populations and their implications on disease is still limited. Low-cost sequencing-based approaches make it possible to target specific molecular markers and genes of interest, leading to potential insights into genetic diversity. Here we demonstrate the versatility of custom dual-indexing technology and Illumina next generation sequencing to generate a genetic profile of human polymorphisms associated with malaria pathology. For 100 individuals diagnosed with severe malaria in Northeast Tanzania, variants were successfully characterised on the haemoglobin subunit beta (HBB), glucose-6-phosphate dehydrogenase (G6PD), atypical chemokine receptor 1 (ACKR1) genes, and the intergenic Dantu genetic blood variant, then validated using pre-existing genotyping data. High sequencing coverage was observed across all amplicon targets in HBB, G6PD, ACKR1, and the Dantu blood group, with variants identified at frequencies previously observed within this region of Tanzania. Sequencing data exhibited high concordance rates to pre-existing genotyping data (> 99.5%). Our work demonstrates the potential utility of amplicon sequencing for applications in human genetics, including to personalise medicine and understand the genetic diversity of loci linked to important host phenotypes, such as malaria susceptibility.

Clinical profiling of MRD48 and functional characterization of two novel pathogenic RAC1 variants

RAC1 is a member of the Rac/Rho GTPase subfamily within the RAS superfamily of small GTP-binding proteins, comprising 3 paralogs playing a critical role in actin cytoskeleton remodeling, cell migration, proliferation and differentiation. De novo missense variants in RAC1 are associated with a rare neurodevelopmental disorder (MRD48) characterized by DD/ID and brain abnormalities coupled with a wide range of additional features. Structural and functional studies have documented either a dominant negative or constitutively active behavior for a subset of mutations. Here, we describe two individuals with previously unreported de novo missense RAC1 variants. We functionally demonstrate their pathogenicity proving a gain-of-function (GoF) effect for both. By reviewing the clinical features of these two individuals and the previously published MRD48 subjects, we further delineate the clinical profile of the disorder, confirming its phenotypic variability. Moreover, we compare the main features of MRD48 with the neurodevelopmental disease caused by GoF variants in the paralog RAC3, highlighting similarities and differences. Finally, we review all previously reported variants in RAC proteins and in the closely related CDC42, providing an updated overview of the spectrum and hotspots of pathogenic variants affecting these functionally related GTPases.

CDK13-related disorder: a deep characterization of speech and language abilities and addition of 33 novel cases

Speech and language impairments are central features of CDK13-related disorder. While pathogenic CDK13 variants have been associated with childhood apraxia of speech (CAS), a systematic characterisation of communication has not been conducted. Here we examined speech, language, non-verbal communication skills, social behaviour and health and development in 41 individuals with CDK13-related disorder from 10 countries (male = 22, median-age 7 years 1 month, range 1-25 years; 33 novel). Most participants used augmentative and alternative communication (AAC) in early childhood (24/41). CAS was common (14/22). Performance varied widely across intellectual ability, social behaviour and expressive language skills, with participants ranging from within average through to the severely impaired range. Receptive language was significantly stronger than expressive language ability. Social motivation was a relative strength. In terms of a broader health phenotype, a quarter had one or more of: renal, urogenital, musculoskeletal, and cardiac malformations, vision impairment, ear infections and/or sleep disturbance. All had gross and fine motor impairments (41/41). Other conditions included mild-moderate intellectual disability (16/22) and autism (7/41). No genotype-phenotype correlations were found. Recognition of CAS, a rare speech disorder, is required to ensure appropriately targeted therapy. The high prevalence of speech and language impairment underscores the importance of tailored speech therapy, particularly early access to AAC supports.