Genomic analysis of childhood hearing loss in the Yoruba population of Nigeria

Although variant alleles of hundreds of genes are associated with sensorineural deafness in children, the genes and alleles involved remain largely unknown in the Sub-Saharan regions of Africa. We ascertained 56 small families mainly of Yoruba ethno-lingual ancestry in or near Ibadan, Nigeria, that had at least one individual with nonsyndromic, severe-to-profound, prelingual-onset, bilateral hearing loss not attributed to nongenetic factors. We performed a combination of exome and Sanger sequencing analyses to evaluate both nuclear and mitochondrial genomes. No biallelic pathogenic variants were identified in GJB2, a common cause of deafness in many populations. Potential causative variants were identified in genes associated with nonsyndromic hearing loss (CIB2, COL11A1, ILDR1, MYO15A, TMPRSS3, and WFS1), nonsyndromic hearing loss or Usher syndrome (CDH23, MYO7A, PCDH15, and USH2A), and other syndromic forms of hearing loss (CHD7, OPA1, and SPTLC1). Several rare mitochondrial variants, including m.1555A>G, were detected in the gene MT-RNR1 but not in control Yoruba samples. Overall, 20 (33%) of 60 independent cases of hearing loss in this cohort of families were associated with likely causal variants in genes reported to underlie deafness in other populations. None of these likely causal variants were present in more than one family, most were detected as compound heterozygotes, and 77% had not been previously associated with hearing loss. These results indicate an unusually high level of genetic heterogeneity of hearing loss in Ibadan, Nigeria and point to challenges for molecular genetic screening, counseling, and early intervention in this population.Subject terms: Genetics research, Medical genomics     

PRICKLE2 revisited—further evidence implicating PRICKLE2 in neurodevelopmental disorders

PRICKLE2 encodes a member of a highly conserved family of proteins that are involved in the non-canonical Wnt and planar cell polarity signaling pathway. Prickle2 localizes to the post-synaptic density, and interacts with post-synaptic density protein 95 and the NMDA receptor. Loss-of-function variants in prickle2 orthologs cause seizures in flies and mice but evidence for the role of PRICKLE2 in human disease is conflicting. Our goal is to provide further evidence for the role of this gene in humans and define the phenotypic spectrum of PRICKLE2-related disorders. We report a cohort of six subjects from four unrelated families with heterozygous rare PRICKLE2 variants ({"type":"entrez-nucleotide","attrs":{"text":"NM_198859.4","term_id":"1519316443","term_text":"NM_198859.4"}}NM_198859.4). Subjects were identified through an international collaboration. Detailed phenotypic and genetic assessment of the subjects were carried out and in addition, we assessed the variant pathogenicity using bioinformatic approaches. We identified two missense variants (c.122 C > T; p.(Pro41Leu), c.680 C > G; p.(Thr227Arg)), one nonsense variant (c.214 C > T; p.(Arg72*) and one frameshift variant (c.1286_1287delGT; p.(Ser429Thrfs*56)). While the p.(Ser429Thrfs*56) variant segregated with disease in a family with three affected females, the three remaining variants occurred de novo. Subjects shared a mild phenotype characterized by global developmental delay, behavioral difficulties ± epilepsy, autistic features, and attention deficit hyperactive disorder. Computational analysis of the missense variants suggest that the altered amino acid residues are likely to be located in protein regions important for function. This paper demonstrates that PRICKLE2 is involved in human neuronal development and that pathogenic variants in PRICKLE2 cause neurodevelopmental delay, behavioral difficulties and epilepsy in humans.Subject terms: Autism spectrum disorders, Genetics of the nervous system, Paediatric neurological disorders     

A detailed clinical and molecular survey of subjects with nonsyndromic USH2A retinopathy reveals an allelic hierarchy of disease-causing variants

Defects in USH2A cause both isolated retinal disease and Usher syndrome (ie, retinal disease and deafness). To gain insights into isolated/nonsyndromic USH2A retinopathy, we screened USH2A in 186 probands with recessive retinal disease and no hearing complaint in childhood (discovery cohort) and in 84 probands with recessive retinal disease (replication cohort). Detailed phenotyping, including retinal imaging and audiological assessment, was performed in individuals with two likely disease-causing USH2A variants. Further genetic testing, including screening for a deep-intronic disease-causing variant and large deletions/duplications, was performed in those with one likely disease-causing change. Overall, 23 of 186 probands (discovery cohort) were found to harbour two likely disease-causing variants in USH2A. Some of these variants were predominantly associated with nonsyndromic retinal degeneration (‘retinal disease-specific''); these included the common c.2276 G>T, p.(Cys759Phe) mutation and five additional variants: c.2802 T>G, p.(Cys934Trp); c.10073 G>A, p.(Cys3358Tyr); c.11156 G>A, p.(Arg3719His); c.12295-3 T>A; and c.12575 G>A, p.(Arg4192His). An allelic hierarchy was observed in the discovery cohort and confirmed in the replication cohort. In nonsyndromic USH2A disease, retinopathy was consistent with retinitis pigmentosa and the audiological phenotype was variable. USH2A retinopathy is a common cause of nonsyndromic recessive retinal degeneration and has a different mutational spectrum to that observed in Usher syndrome. The following model is proposed: the presence of at least one ‘retinal disease-specific'' USH2A allele in a patient with USH2A-related disease results in the preservation of normal hearing. Careful genotype–phenotype studies such as this will become increasingly important, especially now that high-throughput sequencing is widely used in the clinical setting.     

Novel biallelic variants expand the SLC5A6-related phenotypic spectrum

The sodium (Na⁺):multivitamin transporter (SMVT), encoded by SLC5A6, belongs to the sodium:solute symporter family and is required for the Na⁺-dependent uptake of biotin (vitamin B7), pantothenic acid (vitamin B5), the vitamin-like substance α-lipoic acid, and iodide. Compound heterozygous SLC5A6 variants have been reported in individuals with variable multisystemic disorder, including failure to thrive, developmental delay, seizures, cerebral palsy, brain atrophy, gastrointestinal problems, immunodeficiency, and/or osteopenia. We expand the phenotypic spectrum associated with biallelic SLC5A6 variants affecting function by reporting five individuals from three families with motor neuropathies. We identified the homozygous variant c.1285 A > G [p.(Ser429Gly)] in three affected siblings and a simplex patient and the maternally inherited c.280 C > T [p.(Arg94*)] variant and the paternally inherited c.485 A > G [p.(Tyr162Cys)] variant in the simplex patient of the third family. Both missense variants were predicted to affect function by in silico tools. 3D homology modeling of the human SMVT revealed 13 transmembrane helices (TMs) and Tyr162 and Ser429 to be located at the cytoplasmic facing region of TM4 and within TM11, respectively. The SLC5A6 missense variants p.(Tyr162Cys) and p.(Ser429Gly) did not affect plasma membrane localization of the ectopically expressed multivitamin transporter suggesting reduced but not abolished function, such as lower catalytic activity. Targeted therapeutic intervention yielded clinical improvement in four of the five patients. Early molecular diagnosis by exome sequencing is essential for timely replacement therapy in affected individuals.Subject terms: Peripheral neuropathies, Genetics research     

Progressive hearing loss and vestibular dysfunction caused by a homozygous nonsense mutation in CLIC5

In a consanguineous Turkish family diagnosed with autosomal recessive nonsyndromic hearing impairment (arNSHI), a homozygous region of 47.4 Mb was shared by the two affected siblings on chromosome 6p21.1-q15. This region contains 247 genes including the known deafness gene MYO6. No pathogenic variants were found in MYO6, neither with sequence analysis of the coding region and splice sites nor with mRNA analysis. Subsequent candidate gene evaluation revealed CLIC5 as an excellent candidate gene. The orthologous mouse gene is mutated in the jitterbug mutant that exhibits progressive hearing impairment and vestibular dysfunction. Mutation analysis of CLIC5 revealed a homozygous nonsense mutation c.96T>A (p.(Cys32Ter)) that segregated with the hearing loss. Further analysis of CLIC5 in 213 arNSHI patients from mostly Dutch and Spanish origin did not reveal any additional pathogenic variants. CLIC5 mutations are thus not a common cause of arNSHI in these populations. The hearing loss in the present family had an onset in early childhood and progressed from mild to severe or even profound before the second decade. Impaired hearing is accompanied by vestibular areflexia and in one of the patients with mild renal dysfunction. Although we demonstrate that CLIC5 is expressed in many other human tissues, no additional symptoms were observed in these patients. In conclusion, our results show that CLIC5 is a novel arNSHI gene involved in progressive hearing impairment, vestibular and possibly mild renal dysfunction in a family of Turkish origin.     

A novel HRAS c.466C>T p.(Phe156Leu) variant in two patients with attenuated features of Costello syndrome

Costello syndrome (CS) is caused by heterozygous HRAS germline mutations. Most patients share the HRAS variant p.Gly12Ser that is associated with a typical, homogeneous phenotype. Rarer pathogenic HRAS variants (e.g., p.Thr56Ile) were identified in individuals with attenuated CS phenotypes. The obvious phenotypical variability reflects different dysfunctional consequences of distinct HRAS variants. We report on two boys with the novel de novo HRAS variant c.466 C > T p.(Phe156Leu). Both had severe feeding difficulties, airway obstruction and developmental delay, which are typical findings in CS. They showed subtle facial and dermatologic features consistent with attenuated CS. They significantly differed in their musculoskeletal, cardiovascular and endocrinologic manifestations underscoring the clinical variability of individuals with identical, in particular rarer pathogenic HRAS variants. Functional studies revealed enhanced effector-binding, increased downstream signaling activation and impaired growth factor-induced signaling dynamics in cells expressing HRAS^Phe156Leu. Our data further illustrate the molecular and phenotypic variability of CS.Subject terms: Genetics research, Disease genetics     

Targeted Resequencing of Deafness Genes Reveals a Founder MYO15A Variant in Northeastern Brazil

Identifying the genetic etiology in a person with hearing loss (HL) is challenging due to the extreme genetic heterogeneity in HL and the population‐specific variability. In this study, after excluding GJB2 variants, targeted resequencing of 180 deafness‐related genes revealed the causative variants in 11 of 19 (58%) Brazilian probands with autosomal recessive HL. Identified pathogenic variants were in MYO15A (10 families) and CLDN14 (one family). Remarkably, the MYO15A p.(Val1400Met) variant was identified in eight families from the city of Monte Santo in the northeast region of Brazil. Haplotype analysis of this variant was consistent with a single founder. No other cases with this variant were detected among 105 simplex cases from other cities of northeastern Brazil, suggesting that this variant is confined to a geographical region. This study suggests that it is feasible to develop population‐specific screening for deafness variants once causative variants are identified in different geographical groups.     

Sequestosome-1 (SQSTM1) sequence variants in ALS cases in the UK: prevalence and coexistence of SQSTM1 mutations in ALS kindred with PDB

Mutations in the SQSTM1 gene have been reported to be associated with amyotrophic lateral sclerosis (ALS). We sought to determine the frequency of these mutations in a UK familial ALS (FALS) cohort. Sequences of all eight exons of the SQSTM1 gene were analysed in index cases from 61 different FALS kindred lacking known FALS mutations. Six exonic variants c.463G>A, p.(Glu155Lys), c.822G>C, p.(Glu274Asp), c.888G>T, p.(=), c.954C>T, p.(=), c.1038G>A, p.(=) and c.1175C>T, p.(Pro392Leu) were identified in five FALS index cases, three of which were non-synonymous and three were synonymous. One index case harboured three variants (c.822G>C, c.888G>T and c.954C>T), and a second index case harboured two variants (c.822G>C and c.954C>T). Only the p.(Pro392Leu) and p.(Glu155Lys) mutations were predicted to be pathogenic. In one p.(Pro392Leu) kindred, the carrier developed both ALS and Paget''s disease of bone (PDB), and, in the p.(Glu155Lys) kindred, the father of the proband developed PDB. All p.(Pro392Leu) carriers were heterozygous for a previously reported founder haplotype for PDB, where this mutation has an established causal effect. The frequency of the p.(Pro392Leu) mutation in this UK FALS cohort was 2.3% and 0.97% overall including three previously screened FALS cohorts. Our results confirm the presence of the p.(Pro392Leu) SQSTM1 mutation in FALS. This mutation is the most common SQSTM1 mutation found in ALS to date, and a likely pathogenicity is supported by having an established causal role in PDB. The occurrence of the same mutation in ALS and PDB is indicative of a common pathogenic pathway that converges on protein homeostasis.     

A systematic review assessing the existence of pneumothorax-only variants of FLCN. Implications for lifelong surveillance of renal tumours

Individuals with Birt–Hogg–Dubé syndrome (BHDS) may develop fibrofolliculomas, pneumothorax and/or renal cell carcinoma (RCC). Currently, all patients with pathogenic FLCN variants are recommended to have renal surveillance. It has however been suggested that some FLCN variants only cause pneumothorax, which would make surveillance unnecessary in certain cases. This review assesses this possibility. We provide an up-to-date analysis of clinical and genetic features of BHDS. The PUBMED database was systematically searched to find all articles describing patients with pathogenic FLCN variants. The relevant clinical and genetic features of these patients were recorded and analysed. The prevalence of pneumothorax, pulmonary cysts, RCC and characteristic skin lesions in BHDS were 50.9% (n = 1038), 91.9% (n = 720), 22.5% (n = 929) and 47.9% (n = 989), respectively. There was a higher prevalence of pneumothoraces (p < 0.0001) but lower prevalence of dermatological findings (p < 0.0001) in patients from East Asia compared to North America or Europe. Of the 194 pathogenic FLCN variants, 76 could be defined as ‘pneumothorax-only’. Pneumothorax only pathogenic variants (POPVs) were distributed throughout the gene, and there were no statistical differences in variant type. The majority of POPVs (65/76) affected no more than three individuals. Individuals with ‘POPVs’ also tended to be younger (45 vs. 47 years, p < 0.05). Many apparent POPVs in the literature could result from variable expressivity, age-related penetrance and other confounding factors. We therefore recommend that all individuals found to carry a pathogenic FLCN variant be enroled in lifelong surveillance for RCC.Subject terms: Clinical genetics, Respiratory tract diseases     

Functional characterization of BRCA1 gene variants by mini-gene splicing assay

Mutational screening of the breast cancer susceptibility gene BRCA1 leads to the identification of numerous pathogenic variants such as frameshift and nonsense variants, as well as large genomic rearrangements. The screening moreover identifies a large number of variants, for example, missense, silent, and intron variants, which are classified as variants of unknown clinical significance owing to the lack of causal evidence. Variants of unknown clinical significance can potentially have an impact on splicing and therefore functional examinations are warranted to classify whether these variants are pathogenic or benign. Here we validate a mini-gene splicing assay by comparing the results of 24 variants with previously published data from RT-PCR analysis on RNA from blood samples/lymphoblastoid cell lines. The analysis showed an overall concordance of 100%. In addition, we investigated 13 BRCA1 variants of unknown clinical significance or putative variants affecting splicing by in silico analysis and mini-gene splicing assay. Both the in silico analysis and mini-gene splicing assay classified six BRCA1 variants as pathogenic (c.80+1G>A, c.132C>T (p.=), c.213−1G>A, c.670+1delG, c.4185+1G>A, and c.5075−1G>C), whereas six BRCA1 variants were classified as neutral (c.-19-22_-19-21dupAT, c.302−15C>G, c.547+14delG, c.4676−20A>G, c.4987−21G>T, and c.5278−14C>G) and one BRCA1 variant remained unclassified (c.670+16G>A). In conclusion, our study emphasizes that in silico analysis and mini-gene splicing assays are important for the classification of variants, especially if no RNA is available from the patient. This knowledge is crucial for proper genetic counseling of patients and their family members.     

DNA Diagnostics of Hereditary Hearing Loss: A Targeted Resequencing Approach Combined with a Mutation Classification System

Although there are nearly 100 different causative genes identified for nonsyndromic hearing loss (NSHL), Sanger sequencing‐based DNA diagnostics usually only analyses three, namely, GJB2, SLC26A4, and OTOF. As this is seen as inadequate, there is a need for high‐throughput diagnostic methods to detect disease‐causing variations, including single‐nucleotide variations (SNVs), insertions/deletions (Indels), and copy‐number variations (CNVs). In this study, a targeted resequencing panel for hearing loss was developed including 79 genes for NSHL and selected forms of syndromic hearing loss. One‐hundred thirty one presumed autosomal‐recessive NSHL (arNSHL) patients of Western‐European ethnicity were analyzed for SNVs, Indels, and CNVs. In addition, we established a straightforward variant classification system to deal with the large number of variants encountered. We estimate that combining prescreening of GJB2 with our panel leads to a diagnosis in 25%–30% of patients. Our data show that after GJB2, the most commonly mutated genes in a Western‐European population are TMC1, MYO15A, and MYO7A (3.1%). CNV analysis resulted in the identification of causative variants in two patients in OTOA and STRC. One of the major challenges for diagnostic gene panels is assigning pathogenicity for variants. A collaborative database collecting all identified variants from multiple centers could be a valuable resource for hearing loss diagnostics.     

Targeted and Genomewide NGS Data Disqualify Mutations in MYO1A,the “DFNA48 Gene”, as a Cause of Deafness

MYO1A is considered the gene underlying autosomal dominant nonsyndromic hearing loss DFNA48, based on six missense variants, one small in‐frame insertion, and one nonsense mutation. Results from NGS targeting 66 deafness genes in 109 patients identified three families challenging this assumption: two novel nonsense (p.Tyr740* and p.Arg262*) and a known missense variant were identified heterozygously not only in index patients, but also in unaffected relatives. Deafness in these families clearly resulted from mutations in other genes (MYO7A, EYA1, and CIB2). Most of the altogether 10 MYO1A mutations are annotated in dbSNP, and population frequencies (dbSNP, 1000 Genomes, Exome Sequencing Project) above 0.1% contradict pathogenicity under a dominant model. One healthy individual was even homozygous for p.Arg262*, compatible with homozygous Myo1a knockout mice lacking any overt pathology. MYO1A seems dispensable for hearing and overall nonessential. MYO1A adds to the list of “erroneous disease genes”, which will expand with increasing availability of large‐scale sequencing data.     

PNPT1, MYO15A,PTPRQ, and SLC12A2-associated genetic and phenotypic heterogeneity among hearing impaired assortative mating families in Southern India

The study was conducted between 2018 and 2020. From a cohort of 113 hearing impaired (HI), five non-DFNB12 probands identified with heterozygous CDH23 variants were subjected to exome analysis. This resolved the etiology of hearing loss (HL) in four South Indian assortative mating families. Six variants, including three novel ones, were identified in four genes: PNPT1 p.(Ala46Gly) and p.(Asn540Ser), MYO15A p.(Leu1485Pro) and p.(Tyr1891Ter), PTPRQ p.(Gln1336Ter), and SLC12A2 p.(Pro988Ser). Compound heterozygous PNPT1 variants were associated with DFNB70 causing prelingual profound sensorineural hearing loss (SNHL), vestibular dysfunction, and unilateral progressive vision loss in one family. In the second family, MYO15A variants in the myosin motor domain, including a novel variant, causing DFNB3, were found to be associated with prelingual profound SNHL. A novel PTPRQ variant was associated with postlingual progressive sensorineural/mixed HL and vestibular dysfunction in the third family with DFNB84A. In the fourth family, the SLC12A2 novel variant was found to segregate with severe-to-profound HL causing DFNA78, across three generations. Our results suggest a high level of allelic, genotypic, and phenotypic heterogeneity of HL in these families. This study is the first to report the association of PNPT1, PTPRQ, and SLC12A2 variants with HL in the Indian population.     

Consolidating biallelic SDHD variants as a cause of mitochondrial complex II deficiency

Isolated mitochondrial complex II deficiency is a rare cause of mitochondrial respiratory chain disease. To date biallelic variants in three genes encoding mitochondrial complex II molecular components have been unequivocally associated with mitochondrial disease (SDHA/SDHB/SDHAF1). Additionally, variants in one further complex II component (SDHD) have been identified as a candidate cause of isolated mitochondrial complex II deficiency in just two unrelated affected individuals with clinical features consistent with mitochondrial disease, including progressive encephalomyopathy and lethal infantile cardiomyopathy. We present clinical and genomic investigations in four individuals from an extended Palestinian family with clinical features consistent with an autosomal recessive mitochondrial complex II deficiency, in which our genomic studies identified a homozygous NM_003002.3:c.[205 G > A];[205 G > A];p.[(Glu69Lys)];[(Glu69Lys)] SDHD variant as the likely cause. Reviewing previously published cases, these findings consolidate disruption of SDHD function as a cause of mitochondrial complex II deficiency and further define the phenotypic spectrum associated with SDHD gene variants.Subject terms: Disease genetics, Genetic counselling, Metabolic disorders, Genetics research     

Variant pathogenic prediction by locus variability: the importance of the current picture of evolution

Accurate detection of pathogenic single nucleotide variants (SNVs) is a key challenge in whole exome and whole genome sequencing studies. To date, several in silico tools have been developed to predict deleterious variants from this type of data. However, these tools have limited power to detect new pathogenic variants, especially in non-coding regions. In this study, we evaluate the use of a new metric, the Shannon Entropy of Locus Variability (SELV), calculated as the Shannon entropy of the variant frequencies reported in genome-wide population studies at a given locus, as a new predictor of potentially pathogenic variants in non-coding nuclear and mitochondrial DNA and also in coding regions with a selective pressure other than that imposed by the genetic code, e.g splice-sites. For benchmarking, SELV was compared to predictors of pathogenicity in different genomic contexts. In nuclear non-coding DNA, SELV outperformed CDTS (AUC_SELV = 0.97 in ROC curve and PR-AUC_SELV = 0.96 in Precision-recall curve). For non-coding mitochondrial variants (AUC_SELV = 0.98 in ROC curve and PR-AUC_SELV = 1.00 in Precision-recall curve) SELV outperformed HmtVar. Moreover, SELV was compared against two state-of-the-art ensemble predictors of pathogenicity in splice-sites, ada-score, and rf-score, matching their overall performance both in ROC (AUC_SELV = 0.95) and Precision-recall curves (PR-AUC = 0.97), with the advantage that SELV can be easily calculated for every position in the genome, as opposite to ada-score and rf-score. Therefore, we suggest that the information about the observed genetic variability in a locus reported from large scale population studies could improve the prioritization of SNVs in splice-sites and in non-coding regions.Subject terms: Computational biology and bioinformatics, High-throughput screening     

Autosomal recessive cataract (CTRCT18) in the Yakut population isolate of Eastern Siberia: a novel founder variant in the FYCO1 gene

Congenital autosomal recessive cataract with unknown genetic etiology is one of the most common Mendelian diseases among the Turkic-speaking Yakut population (Eastern Siberia, Russia). To identify the genetic cause of congenital cataract spread in this population, we performed whole-exome sequencing (Illumina NextSeq 500) in one Yakut family with three affected siblings whose parents had preserved vision. We have revealed the novel homozygous c.1621C>T transition leading to premature stop codon p.(Gln541*) in exon 8 of the FYCO1 gene ({"type":"entrez-nucleotide","attrs":{"text":"NM_024513.4","term_id":"1519312749"}}NM_024513.4). Subsequent screening of c.1621C>T p.(Gln541*) revealed this variant in a homozygous state in 25 out of 29 Yakut families with congenital cataract (86%). Among 424 healthy individuals from seven populations of Eastern Siberia (Russians, Yakuts, Evenks, Evens, Dolgans, Chukchi, and Yukaghirs), the highest carrier frequency of c.1621C>T p.(Gln541*) was found in the Yakut population (7.9%). DNA samples of 25 homozygous for c.1621C>T p.(Gln541*) patients with congenital cataract and 114 unaffected unrelated individuals without this variant were used for a haplotype analysis based on the genotyping of six STR markers (D3S3512, D3S3685, D3S3582, D3S3561, D3S1289, and D3S3698). The structure of the identified haplotypes indicates a common origin for all of the studied mutant chromosomes bearing c.1621C>T p.(Gln541*). The age of the с.1621C>T p.(Gln541*) founder haplotype was estimated to be approximately 260 ± 65 years (10 generations). These findings characterize Eastern Siberia as the region of the world with the most extensive accumulation of the unique variant c.1621C>T p.(Gln541*) in the FYCO1 gene as a result of the founder effect.Subject terms: Visual system, Genetics research, Clinical genetics     

Clinical and molecular characterization of the potential CF disease modifier syntaxin 1A

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator gene (CFTR). Disease severity in CF varies greatly, and sibling studies strongly indicate that genes other than CFTR modify disease outcome. Syntaxin 1A (STX1A) has been reported as a negative regulator of CFTR and other ion channels. We hypothesized that STX1A variants act as a CF modifier by influencing the remaining function of mutated CFTR. We identified STX1A variants by genomic resequencing patients from the Bernese CF Patient Data Registry and applied linear mixed model analysis to establish genotype–phenotype correlations, revealing STX1A rs4363087 (c.467−38A>G) to significantly influence lung function. The same STX1A risk allele was recognized in the European CF Twin and Sibling Study (P=0.0027), demonstrating that the genotype–phenotype association of STX1A to CF disease severity is robust enough to allow replication in two independent CF populations. rs4363087 is in linkage disequilibrium to the exonic variant rs2228607 (c.204C>T). Considering that neither rs4363087 nor rs2228607 changes the amino-acid sequence of STX1A, we investigated their effects on mRNA level. We show that rs2228607 reinforces aberrant splicing of STX1A mRNA, leading to nonsense-mediated mRNA decay. In conclusion, we demonstrate the clinical relevance of STX1A variants in CF, and evidence the functional relevance of STX1A variant rs2228607 at molecular level. Our findings show that genes interacting with CFTR can modify CF disease progression.