Multi-gene panel sequencing in highly consanguineous families and patients with congenital forms of skeletal dysplasias.

Skeletal dysplasias (SKDs) are a heterogeneous group of more than 750 genetic disorders characterized by abnormal development, growth, and maintenance of bones or cartilage in the human skeleton. SKDs are often caused by variants in early patterning genes and in many cases part of multiple malformation syndromes and occur in combination with non-skeletal phenotypes. The aim of this study was to investigate the underlying genetic cause of congenital SKDs in highly consanguineous Pakistani families, as well as in sporadic and familial SKD cases from India using multigene panel sequencing analysis. Therefore, we performed panel sequencing of 386 bone-related genes in 7 highly consanguineous families from Pakistan and 27 cases from India affected with SKDs. In the highly consanguineous families, we were able to identify the underlying genetic cause in five out of seven families, resulting in a diagnostic yield of 71%. Whereas, in the sporadic and familial SKD cases, we identified 12 causative variants, corresponding to a diagnostic yield of 44%. The genetic heterogeneity in our cohorts was very high and we were able to detect various types of variants, including missense, nonsense, and frameshift variants, across multiple genes known to cause different types of SKDs. In conclusion, panel sequencing proved to be a highly effective way to decipher the genetic basis of SKDs in highly consanguineous families as well as sporadic and or familial cases from South Asia. Furthermore, our findings expand the allelic spectrum of skeletal dysplasias.

Quality of Life After Deep Brain Stimulation of Pediatric Patients with Dyskinetic Cerebral Palsy: A Prospective, Single-Arm, Multicenter Study with a Subsequent Randomized Double-Blind Crossover (STIM-CP).

BACKGROUND: Patients with dyskinetic cerebral palsy are often severely impaired with limited treatment options. The effects of deep brain stimulation (DBS) are less pronounced than those in inherited dystonia but can be associated with favorable quality of life outcomes even in patients without changes in dystonia severity. OBJECTIVE: The aim is to assess DBS effects in pediatric patients with pharmacorefractory dyskinetic cerebral palsy with focus on quality of life. METHODS: The method used is a prospective, single-arm, multicenter study. The primary endpoint is improvement in quality of life (CPCHILD [Caregiver Priorities & Child Health Index of Life with Disabilities]) from baseline to 12 months under therapeutic stimulation. The main key secondary outcomes are changes in Burke-Fahn-Marsden Dystonia Rating Scale, Dyskinesia Impairment Scale, Gross Motor Function Measure-66, Canadian Occupational Performance Measure (COPM), and Short-Form (SF)-36. After 12 months, patients were randomly assigned to a blinded crossover to receive active or sham stimulation for 24 hours each. Severity of dystonia and chorea were blindly rated. Safety was assessed throughout. The trial was registered at ClinicalTrials.gov, number NCT02097693. RESULTS: Sixteen patients (age: 13.4 ± 2.9 years) were recruited by seven clinical sites. Primary outcome at 12-month follow-up is as follows: mean CPCHILD increased by 4.2 ± 10.4 points (95% CI [confidence interval] -1.3 to 9.7; P = 0.125); among secondary outcomes: improvement in COPM performance measure of 1.1 ± 1.5 points (95% CI 0.2 to 1.9; P = 0.02) and in the SF-36 physical health component by 5.1 ± 6.2 points (95% CI 0.7 to 9.6; P = 0.028). Otherwise, there are no significant changes. CONCLUSION: Evidence to recommend DBS as routine treatment to improve quality of life in pediatric patients with dyskinetic cerebral palsy is not yet sufficient. Extended follow-up in larger cohorts will determine the impact of DBS further to guide treatment decisions in these often severely disabled patients. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features.

PURPOSE: The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10. METHODS: Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry. RESULTS: A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein. CONCLUSION: We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

The murine ortholog of Kaufman oculocerebrofacial syndrome protein Ube3b regulates synapse number by ubiquitinating Ppp3cc.

Kaufman oculocerebrofacial syndrome (KOS) is a severe autosomal recessive disorder characterized by intellectual disability, developmental delays, microcephaly, and characteristic dysmorphisms. Biallelic mutations of UBE3B, encoding for a ubiquitin ligase E3B are causative for KOS. In this report, we characterize neuronal functions of its murine ortholog Ube3b and show that Ube3b regulates dendritic branching in a cell-autonomous manner. Moreover, Ube3b knockout (KO) neurons exhibit increased density and aberrant morphology of dendritic spines, altered synaptic physiology, and changes in hippocampal circuit activity. Dorsal forebrain-specific Ube3b KO animals show impaired spatial learning, altered social interactions, and repetitive behaviors. We further demonstrate that Ube3b ubiquitinates the catalytic γ-subunit of calcineurin, Ppp3cc, the overexpression of which phenocopies Ube3b loss with regard to dendritic spine density. This work provides insights into the molecular pathologies underlying intellectual disability-like phenotypes in a genetically engineered mouse model.

SQSTM1/p62 variants in 486 patients with familial ALS from Germany and Sweden.

Several studies reported amyotrophic lateral sclerosis (ALS)-linked mutations in TBK1, OPTN, VCP, UBQLN2, and SQSTM1 genes encoding proteins involved in autophagy. SQSTM1 was originally identified by a candidate gene approach because it encodes p62, a multifunctional protein involved in protein degradation both through proteasomal regulation and autophagy. Both p62 and optineurin (encoded by OPTN) are direct interaction partners and substrates of TBK1, and these 3 proteins form the core of a genetic and functional network that may connect autophagy with ALS. Considering the molecular and conceptual relevance of the TBK1/OPTN/SQSTM1 "triangle," we here performed a targeted screen for SQSTM1 variants in 486 patients with familial ALS from Germany and Sweden by analyzing whole-exome sequencing data. We report 9 novel and 5 previously reported rare variants in SQSTM1 and discuss the current evidence for SQSTM1 as a primary disease gene for ALS. We conclude that the evidence for causality remains vague for SQSTM1 and is weaker than for the other autophagy genes, for example, TBK1 and OPTN.

Variants in KIAA0825 underlie autosomal recessive postaxial polydactyly.

Postaxial polydactyly (PAP) is a common limb malformation that often leads to cosmetic and functional complications. Molecular evaluation of polydactyly can serve as a tool to elucidate genetic and signaling pathways that regulate limb development, specifically, the anterior-posterior specification of the limb. To date, only five genes have been identified for nonsyndromic PAP: FAM92A, GLI1, GLI3, IQCE and ZNF141. In this study, two Pakistani multiplex consanguineous families with autosomal recessive nonsyndromic PAP were clinically and molecularly evaluated. From both pedigrees, a DNA sample from an affected member underwent exome sequencing. In each family, we identified a segregating frameshift (c.591dupA [p.(Q198Tfs*21)]) and nonsense variant (c.2173A > T [p.(K725*)]) in KIAA0825 (also known as C5orf36). Although KIAA0825 encodes a protein of unknown function, it has been demonstrated that its murine ortholog is expressed during limb development. Our data contribute to the establishment of a catalog of genes important in limb patterning, which can aid in diagnosis and obtaining a better understanding of the biology of polydactyly.

Next-generation sequencing of 32 genes associated with hereditary aortopathies and related disorders of connective tissue in a cohort of 199 patients.

PURPOSE: Heritable factors play an important etiologic role in connective tissue disorders (CTD) with vascular involvement, and a genetic diagnosis is getting increasingly important for gene-tailored, personalized patient management. METHODS: We analyzed 32 disease-associated genes by using targeted next-generation sequencing and exome sequencing in a clinically relevant cohort of 199 individuals. We classified and refined sequence variants according to their likelihood for pathogenicity. RESULTS: We identified 1 pathogenic variant (PV; in FBN1 or SMAD3) in 15 patients (7.5%) and ≥1 likely pathogenic variant (LPV; in COL3A1, FBN1, FBN2, LOX, MYH11, SMAD3, TGFBR1, or TGFBR2) in 19 individuals (9.6%), together resulting in 17.1% diagnostic yield. Thirteen PV/LPV were novel. Of PV/LPV-negative patients 47 (23.6%) showed ≥1 variant of uncertain significance (VUS). Twenty-five patients had concomitant variants. In-depth evaluation of reported/calculated variant classes resulted in reclassification of 19.8% of variants. CONCLUSION: Variant classification and refinement are essential for shaping mutational spectra of disease genes, thereby improving clinical sensitivity. Obligate stringent multigene analysis is a powerful tool for identifying genetic causes of clinically related CTDs. Nonetheless, the relatively high rate of PV/LPV/VUS-negative patients underscores the existence of yet unknown disease loci and/or oligogenic/polygenic inheritance.

Unknown mutations and genotype/phenotype correlations of autosomal recessive congenital ichthyosis in patients from Saudi Arabia and Pakistan.

BACKGROUND: Autosomal recessive congenital ichthyosis (ARCI) is a genetically and phenotypically heterogeneous skin disease, associated with defects in the skin permeability barrier. Several but not all genes with underlying mutations have been identified, but a clear correlation between genetic causes and clinical picture has not been described to date. METHODS: Our study included 19 families from Saudi Arabia, Yemen, and Pakistan. All patients were born to consanguineous parents and diagnosed with ARCI. Mutations were analyzed by homozygosity mapping and direct sequencing. RESULTS: We have detected mutations in all families in five different genes: TGM1, ABCA12, CYP4F22, NIPAL4, and ALOXE3. Five likely pathogenic variants were unknown so far, a splice site and a missense variant in TGM1, a splice site variant in NIPAL4, and missense variants in ABCA12 and CYP4F22. We attributed TGM1 and ABCA12 mutations to the most severe forms of lamellar and erythematous ichthyoses, respectively, regardless of treatment. Other mutations highlighted the presence of a phenotypic spectrum in ARCI. CONCLUSION: Our results contribute to expanding the mutational spectrum of ARCI and revealed new insights into genotype/phenotype correlations. The findings are instrumental for a faster and more precise diagnosis, a better understanding of the pathophysiology, and the definition of targets for more specific therapies for ARCI.

Dominant Noonan syndrome-causing LZTR1 mutations specifically affect the Kelch domain substrate-recognition surface and enhance RAS-MAPK signaling.

Noonan syndrome (NS), the most common RASopathy, is caused by mutations affecting signaling through RAS and the MAPK cascade. Recently, genome scanning has discovered novel genes implicated in NS, whose function in RAS-MAPK signaling remains obscure, suggesting the existence of unrecognized circuits contributing to signal modulation in this pathway. Among these genes, leucine zipper-like transcriptional regulator 1 (LZTR1) encodes a functionally poorly characterized member of the BTB/POZ protein superfamily. Two classes of germline LZTR1 mutations underlie dominant and recessive forms of NS, while constitutional monoallelic, mostly inactivating, mutations in the same gene cause schwannomatosis, a cancer-prone disorder clinically distinct from NS. Here we show that dominant NS-causing LZTR1 mutations do not affect significantly protein stability and subcellular localization. We provide the first evidence that these mutations, but not the missense changes occurring as biallelic mutations in recessive NS, enhance stimulus-dependent RAS-MAPK signaling, which is triggered, at least in part, by an increased RAS protein pool. Moreover, we document that dominant NS-causing mutations do not perturb binding of LZTR1 to CUL3, a scaffold coordinating the assembly of a multimeric complex catalyzing protein ubiquitination but are predicted to affect the surface of the Kelch domain mediating substrate binding to the complex. Collectively, our data suggest a model in which LZTR1 contributes to the ubiquitinationof protein(s) functioning as positive modulator(s) of the RAS-MAPK signaling pathway. In this model, LZTR1 mutations are predicted to variably impair binding of these substrates to the multi-component ligase complex and their efficient ubiquitination and degradation, resulting in MAPK signaling upregulation.

Functional and Phenotypic Characteristics of Human Leptin Receptor Mutations.

Several case series of extreme early-onset obesity due to mutations in the human leptin receptor (LEPR) gene have been reported. In this review we summarize published functional and phenotypic data on mutations in the human LEPR gene causing severe early-onset obesity. Additionally, we included data on six new cases from our obesity center. Literature research was performed using PubMed and OMIM. Functional relevance of mutations was estimated based on reported functional analysis, mutation size, and location, as well as phenotypic characteristics of affected patients. We identified 57 cases with 38 distinct LEPR mutations. We found severe early-onset obesity, hyperphagia, and hypogonadotropic hypogonadism as cardinal features of a complete loss of LEPR function. Other features, for example, metabolic disorders and recurring infections, were variable in manifestation. Obesity degree or other manifestations did not aggregate by genotype. Few patients underwent bariatric surgery with variable success. Most mutations occurred in the fibronectin III and cytokine receptor homology II domains, whereas none was found in cytoplasmic domain. In silico data were available for 25 mutations and in vitro data were available for four mutations, revealing residual activity in one case. By assessing provided information on the clinical phenotype, functional analysis, and character of the 38 mutations, we assume residual LEPR activity for five additional mutations. Functional in vitro analysis is necessary to confirm this assumption.

Elucidating the genetic architecture of Adams-Oliver syndrome in a large European cohort.

Adams-Oliver syndrome (AOS) is a rare developmental disorder, characterized by scalp aplasia cutis congenita (ACC) and transverse terminal limb defects (TTLD). Autosomal dominant forms of AOS are linked to mutations in ARHGAP31, DLL4, NOTCH1 or RBPJ, while DOCK6 and EOGT underlie autosomal recessive inheritance. Data on the frequency and distribution of mutations in large cohorts are currently limited. The purpose of this study was therefore to comprehensively examine the genetic architecture of AOS in an extensive cohort. Molecular diagnostic screening of 194 AOS/ACC/TTLD probands/families was conducted using next-generation and/or capillary sequencing analyses. In total, we identified 63 (likely) pathogenic mutations, comprising 56 distinct and 22 novel mutations, providing a molecular diagnosis in 30% of patients. Taken together with previous reports, these findings bring the total number of reported disease variants to 63, with a diagnostic yield of 36% in familial cases. NOTCH1 is the major contributor, underlying 10% of AOS/ACC/TTLD cases, with DLL4 (6%), DOCK6 (6%), ARHGAP31 (3%), EOGT (3%), and RBPJ (2%) representing additional causality in this cohort. We confirm the relevance of genetic screening across the AOS/ACC/TTLD spectrum, highlighting preliminary but important genotype-phenotype correlations. This cohort offers potential for further gene identification to address missing heritability.

Comprehensive analysis of the mutation spectrum in 301 German ALS families.

OBJECTIVES: Recent advances in amyotrophic lateral sclerosis (ALS) genetics have revealed that mutations in any of more than 25 genes can cause ALS, mostly as an autosomal-dominant Mendelian trait. Detailed knowledge about the genetic architecture of ALS in a specific population will be important for genetic counselling but also for genotype-specific therapeutic interventions. METHODS: Here we combined fragment length analysis, repeat-primed PCR, Southern blotting, Sanger sequencing and whole exome sequencing to obtain a comprehensive profile of genetic variants in ALS disease genes in 301 German pedigrees with familial ALS. We report C9orf72 mutations as well as variants in consensus splice sites and non-synonymous variants in protein-coding regions of ALS genes. We furthermore estimate their pathogenicity by taking into account type and frequency of the respective variant as well as segregation within the families. RESULTS: 49% of our German ALS families carried a likely pathogenic variant in at least one of the earlier identified ALS genes. In 45% of the ALS families, likely pathogenic variants were detected in C9orf72, SOD1, FUS, TARDBP or TBK1, whereas the relative contribution of the other ALS genes in this familial ALS cohort was 4%. We identified several previously unreported rare variants and demonstrated the absence of likely pathogenic variants in some of the recently described ALS disease genes. CONCLUSIONS: We here present a comprehensive genetic characterisation of German familial ALS. The present findings are of importance for genetic counselling in clinical practice, for molecular research and for the design of diagnostic gene panels or genotype-specific therapeutic interventions in Europe.

Early childhood BMI trajectories in monogenic obesity due to leptin, leptin receptor, and melanocortin 4 receptor deficiency.

OBJECTIVE: To evaluate whether early childhood body mass index (BMI) is an appropriate indicator for monogenic obesity. METHODS: A cohort of n = 21 children living in Germany or Austria with monogenic obesity due to congenital leptin deficiency (group LEP, n = 6), leptin receptor deficiency (group LEPR, n = 6) and primarily heterozygous MC4 receptor deficiency (group MC4R, n = 9) was analyzed. A control group (CTRL) was defined that consisted of n = 22 obese adolescents with no mutation in the above mentioned genes. Early childhood (0-5 years) BMI trajectories were compared between the groups at selected time points. RESULTS: The LEP and LEPR group showed a tremendous increase in BMI during the first 2 years of life with all patients displaying a BMI >27 kg/m2 (27.2-38.4 kg/m2) and %BMIP95 (percentage of the 95th percentile BMI for age and sex) >140% (144.8-198.6%) at the age of 2 years and a BMI > 33 kg/m2 (33.3-45.9 kg/m2) and %BMIP95 > 184% (184.1-212.6%) at the age of 5 years. The MC4R and CTRL groups had a later onset of obesity with significantly lower BMI values at both time points (p < 0.01). CONCLUSION: As result of the investigation of early childhood BMI trajectories in this pediatric cohort with monogenic obesity we suggest that BMI values >27.0 kg/m2 or %BMIP95 > 140% at the age of 2 years and BMI values >33.0 kg/m2 or %BMIP95 > 184% at the age of 5 years may be useful cut points to identify children who should undergo genetic screening for monogenic obesity due to functionally relevant mutations in the leptin gene or leptin receptor gene.

Autosomal recessive Noonan syndrome associated with biallelic LZTR1 variants.

PURPOSE: To characterize the molecular genetics of autosomal recessive Noonan syndrome. METHODS: Families underwent phenotyping for features of Noonan syndrome in children and their parents. Two multiplex families underwent linkage analysis. Exome, genome, or multigene panel sequencing was used to identify variants. The molecular consequences of observed splice variants were evaluated by reverse-transcription polymerase chain reaction. RESULTS: Twelve families with a total of 23 affected children with features of Noonan syndrome were evaluated. The phenotypic range included mildly affected patients, but it was lethal in some, with cardiac disease and leukemia. All of the parents were unaffected. Linkage analysis using a recessive model supported a candidate region in chromosome 22q11, which includes LZTR1, previously shown to harbor mutations in patients with Noonan syndrome inherited in a dominant pattern. Sequencing analyses of 21 live-born patients and a stillbirth identified biallelic pathogenic variants in LZTR1, including putative loss-of-function, missense, and canonical and noncanonical splicing variants in the affected children, with heterozygous, clinically unaffected parents and heterozygous or normal genotypes in unaffected siblings. CONCLUSION: These clinical and genetic data confirm the existence of a form of Noonan syndrome that is inherited in an autosomal recessive pattern and identify biallelic mutations in LZTR1.

A novel homozygous missense variant in NECTIN4 (PVRL4) causing ectodermal dysplasia cutaneous syndactyly syndrome.

Ectodermal dysplasia syndactyly syndrome 1 (EDSS1) is a rare form of ectodermal dysplasia including anomalies of hair, nails, and teeth along with bilateral cutaneous syndactyly of hands and feet. In the present report, we performed a clinical and genetic characterization of a consanguineous Pakistani family with four individuals affected by EDSS1. We performed exome sequencing using DNA of one affected individual. Exome data analysis identified a novel homozygous missense variant (c.242T>C; p.(Leu81Pro)) in NECTIN4 (PVRL4). Sanger sequencing validated this variant and confirmed its cosegregation with the disease phenotype in the family members. Thus, our report adds a novel variant to the NECTIN4 mutation spectrum and contributes to the NECTIN4-related clinical characterization.

Clinical and genetic spectrum of AMPD2-related pontocerebellar hypoplasia type 9.

Pontocerebellar hypoplasia (PCH) represents a group of autosomal-recessive progressive neurodegenerative disorders of prenatal onset. Eleven PCH subtypes are classified according to clinical, neuroimaging and genetic findings. Individuals with PCH type 9 (PCH9) have a unique combination of postnatal microcephaly, hypoplastic cerebellum and pons, and hypoplastic or absent corpus callosum. PCH9 is caused by biallelic variants in AMPD2 encoding adenosine monophosphate deaminase 2; however, a homozygous AMPD2 frameshift variant has recently been reported in two family members with spastic paraplegia type 63 (SPG63). We identified homozygous or compound heterozygous AMPD2 variants in eight PCH-affected individuals from six families. The eight variants likely affect function and comprise one frameshift, one nonsense and six missense variants; seven of which were novel. The main clinical manifestations in the eight new patients and 17 previously reported individuals with biallelic AMPD2 variants were postnatal microcephaly, severe global developmental delay, spasticity, and central visual impairment. Brain imaging data identified hypomyelination, hypoplasia of the cerebellum and pons, atrophy of the cerebral cortex, complete or partial agenesis of the corpus callosum and the "figure 8" shape of the hypoplastic midbrain as consistent features. We broaden the AMPD2-related clinical spectrum by describing one individual without microcephaly and absence of the characteristic "figure 8" shape of the midbrain. The existence of various AMPD2 isoforms with different functions possibly explains the variability in phenotypes associated with AMPD2 variants: variants leaving some of the isoforms intact may cause SPG63, while those affecting all isoforms may result in the severe and early-onset PCH9.

Hot-spot KIF5A mutations cause familial ALS.

Heterozygous missense mutations in the N-terminal motor or coiled-coil domains of the kinesin family member 5A (KIF5A) gene cause monogenic spastic paraplegia (HSP10) and Charcot-Marie-Tooth disease type 2 (CMT2). Moreover, heterozygous de novo frame-shift mutations in the C-terminal domain of KIF5A are associated with neonatal intractable myoclonus, a neurodevelopmental syndrome. These findings, together with the observation that many of the disease genes associated with amyotrophic lateral sclerosis disrupt cytoskeletal function and intracellular transport, led us to hypothesize that mutations in KIF5A are also a cause of amyotrophic lateral sclerosis. Using whole exome sequencing followed by rare variant analysis of 426 patients with familial amyotrophic lateral sclerosis and 6137 control subjects, we detected an enrichment of KIF5A splice-site mutations in amyotrophic lateral sclerosis (2/426 compared to 0/6137 in controls; P = 4.2 × 10-3), both located in a hot-spot in the C-terminus of the protein and predicted to affect splicing exon 27. We additionally show co-segregation with amyotrophic lateral sclerosis of two canonical splice-site mutations in two families. Investigation of lymphoblast cell lines from patients with KIF5A splice-site mutations revealed the loss of mutant RNA expression and suggested haploinsufficiency as the most probable underlying molecular mechanism. Furthermore, mRNA sequencing of a rare non-synonymous missense mutation (predicting p.Arg1007Gly) located in the C-terminus of the protein shortly upstream of the splice donor of exon 27 revealed defective KIF5A pre-mRNA splicing in respective patient-derived cell lines owing to abrogation of the donor site. Finally, the non-synonymous single nucleotide variant rs113247976 (minor allele frequency = 1.00% in controls, n = 6137), also located in the C-terminal region [p.(Pro986Leu) in exon 26], was significantly enriched in familial amyotrophic lateral sclerosis patients (minor allele frequency = 3.40%; P = 1.28 × 10-7). Our study demonstrates that mutations located specifically in a C-terminal hotspot of KIF5A can cause a classical amyotrophic lateral sclerosis phenotype, and underline the involvement of intracellular transport processes in amyotrophic lateral sclerosis pathogenesis.

Association Between Germline Mutations in BRF1, a Subunit of the RNA Polymerase III Transcription Complex, and Hereditary Colorectal Cancer.

BACKGROUND & AIMS: Although there is a genetic predisposition to colorectal cancer (CRC), few of the genes that affect risk have been identified. We performed whole-exome sequence analysis of individuals in a high-risk family without mutations in genes previously associated with CRC risk to identify variants associated with inherited CRC. METHODS: We collected blood samples from 3 relatives with CRC in Spain (65, 62, and 40 years old at diagnosis) and performed whole-exome sequence analyses. Rare missense, truncating or splice-site variants shared by the 3 relatives were selected. We used targeted pooled DNA amplification followed by next generation sequencing to screen for mutations in candidate genes in 547 additional hereditary and/or early-onset CRC cases (502 additional families). We carried out protein-dependent yeast growth assays and transfection studies in the HT29 human CRC cell line to test the effects of the identified variants. RESULTS: A total of 42 unique or rare (population minor allele frequency below 1%) nonsynonymous genetic variants in 38 genes were shared by all 3 relatives. We selected the BRF1 gene, which encodes an RNA polymerase III transcription initiation factor subunit for further analysis, based on the predicted effect of the identified variant and previous association of BRF1 with cancer. Previously unreported or rare germline variants in BRF1 were identified in 11 of 503 CRC families, a significantly greater proportion than in the control population (34 of 4300). Seven of the identified variants (1 detected in 2 families) affected BRF1 mRNA splicing, protein stability, or expression and/or function. CONCLUSIONS: In an analysis of families with a history of CRC, we associated germline mutations in BRF1 with predisposition to CRC. We associated deleterious BRF1 variants with 1.4% of familial CRC cases, in individuals without mutations in high-penetrance genes previously associated with CRC. Our findings add additional evidence to the link between defects in genes that regulate ribosome synthesis and risk of CRC.