Genetic modifiers of body mass index in individuals with cystic fibrosis.

To identify modifier loci underlying variation in body mass index (BMI) in persons with cystic fibrosis (pwCF), we performed a genome-wide association study (GWAS). Utilizing longitudinal height and weight data, along with demographic information and covariates from 4,393 pwCF, we calculated AvgBMIz representing the average of per-quarter BMI Z scores. The GWAS incorporated 9.8M single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) > 0.005 extracted from whole-genome sequencing (WGS) of each study subject. We observed genome-wide significant association with a variant in FTO (FaT mass and Obesity-associated gene; rs28567725; p value = 1.21e-08; MAF = 0.41, ß = 0.106; n = 4,393 individuals) and a variant within ADAMTS5 (A Disintegrin And Metalloproteinase with ThromboSpondin motifs 5; rs162500; p value = 2.11e-10; MAF = 0.005, ß = -0.768; n = 4,085 pancreatic-insufficient individuals). Notably, BMI-associated variants in ADAMTS5 occur on a haplotype that is much more common in African (AFR, MAF = 0.183) than European (EUR, MAF = 0.006) populations (1000 Genomes project). A polygenic risk score (PRS) calculated using 924 SNPs (excluding 17 in FTO) showed significant association with AvgBMIz (p value = 2.2e-16; r2 = 0.03). Association between variants in FTO and the PRS correlation reveals similarities in the genetic architecture of BMI in CF and the general population. Inclusion of Black individuals in whom the single-gene disorder CF is much less common but genomic diversity is greater facilitated detection of association with variants that are in LD with functional SNPs in ADAMTS5. Our results illustrate the importance of population diversity, particularly when attempting to identify variants that manifest only under certain physiologic conditions.

Whole-exome sequencing uncovers the genetic complexity of bicuspid aortic valve in families with early-onset complications.

Bicuspid aortic valve (BAV) is the most common congenital heart lesion with an estimated population prevalence of 1%. We hypothesize that specific gene variants predispose to early-onset complications of BAV (EBAV). We analyzed whole-exome sequences (WESs) to identify rare coding variants that contribute to BAV disease in 215 EBAV-affected families. Predicted damaging variants in candidate genes with moderate or strong supportive evidence to cause developmental cardiac phenotypes were present in 107 EBAV-affected families (50% of total), including genes that cause BAV (9%) or heritable thoracic aortic disease (HTAD, 19%). After appropriate filtration, we also identified 129 variants in 54 candidate genes that are associated with autosomal-dominant congenital heart phenotypes, including recurrent deleterious variation of FBN2, MYH6, channelopathy genes, and type 1 and 5 collagen genes. These findings confirm our hypothesis that unique rare genetic variants drive early-onset presentations of BAV disease.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Expanding the genetic and phenotypic landscape of replication factor C complex-related disorders: RFC4 deficiency is linked to a multisystemic disorder.

The precise regulation of DNA replication is vital for cellular division and genomic integrity. Central to this process is the replication factor C (RFC) complex, encompassing five subunits, which loads proliferating cell nuclear antigen onto DNA to facilitate the recruitment of replication and repair proteins and enhance DNA polymerase processivity. While RFC1's role in cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS) is known, the contributions of RFC2-5 subunits on human Mendelian disorders is largely unexplored. Our research links bi-allelic variants in RFC4, encoding a core RFC complex subunit, to an undiagnosed disorder characterized by incoordination and muscle weakness, hearing impairment, and decreased body weight. We discovered across nine affected individuals rare, conserved, predicted pathogenic variants in RFC4, all likely to disrupt the C-terminal domain indispensable for RFC complex formation. Analysis of a previously determined cryo-EM structure of RFC bound to proliferating cell nuclear antigen suggested that the variants disrupt interactions within RFC4 and/or destabilize the RFC complex. Cellular studies using RFC4-deficient HeLa cells and primary fibroblasts demonstrated decreased RFC4 protein, compromised stability of the other RFC complex subunits, and perturbed RFC complex formation. Additionally, functional studies of the RFC4 variants affirmed diminished RFC complex formation, and cell cycle studies suggested perturbation of DNA replication and cell cycle progression. Our integrated approach of combining in silico, structural, cellular, and functional analyses establishes compelling evidence that bi-allelic loss-of-function RFC4 variants contribute to the pathogenesis of this multisystemic disorder. These insights broaden our understanding of the RFC complex and its role in human health and disease.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions.

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis.

Individuals with cystic fibrosis (CF) develop complications of the gastrointestinal tract influenced by genetic variants outside of CFTR. Cystic fibrosis-related diabetes (CFRD) is a distinct form of diabetes with a variable age of onset that occurs frequently in individuals with CF, while meconium ileus (MI) is a severe neonatal intestinal obstruction affecting ∼20% of newborns with CF. CFRD and MI are slightly correlated traits with previous evidence of overlap in their genetic architectures. To better understand the genetic commonality between CFRD and MI, we used whole-genome-sequencing data from the CF Genome Project to perform genome-wide association. These analyses revealed variants at 11 loci (6 not previously identified) that associated with MI and at 12 loci (5 not previously identified) that associated with CFRD. Of these, variants at SLC26A9, CEBPB, and PRSS1 associated with both traits; variants at SLC26A9 and CEBPB increased risk for both traits, while variants at PRSS1, the higher-risk alleles for CFRD, conferred lower risk for MI. Furthermore, common and rare variants within the SLC26A9 locus associated with MI only or CFRD only. As expected, different loci modify risk of CFRD and MI; however, a subset exhibit pleiotropic effects indicating etiologic and mechanistic overlap between these two otherwise distinct complications of CF.

MagicalRsq: Machine-learning-based genotype imputation quality calibration.

Whole-genome sequencing (WGS) is the gold standard for fully characterizing genetic variation but is still prohibitively expensive for large samples. To reduce costs, many studies sequence only a subset of individuals or genomic regions, and genotype imputation is used to infer genotypes for the remaining individuals or regions without sequencing data. However, not all variants can be well imputed, and the current state-of-the-art imputation quality metric, denoted as standard Rsq, is poorly calibrated for lower-frequency variants. Here, we propose MagicalRsq, a machine-learning-based method that integrates variant-level imputation and population genetics statistics, to provide a better calibrated imputation quality metric. Leveraging WGS data from the Cystic Fibrosis Genome Project (CFGP), and whole-exome sequence data from UK BioBank (UKB), we performed comprehensive experiments to evaluate the performance of MagicalRsq compared to standard Rsq for partially sequenced studies. We found that MagicalRsq aligns better with true R2 than standard Rsq in almost every situation evaluated, for both European and African ancestry samples. For example, when applying models trained from 1,992 CFGP sequenced samples to an independent 3,103 samples with no sequencing but TOPMed imputation from array genotypes, MagicalRsq, compared to standard Rsq, achieved net gains of 1.4 million rare, 117k low-frequency, and 18k common variants, where net gains were gained numbers of correctly distinguished variants by MagicalRsq over standard Rsq. MagicalRsq can serve as an improved post-imputation quality metric and will benefit downstream analysis by better distinguishing well-imputed variants from those poorly imputed. MagicalRsq is freely available on GitHub.

Genetic variation in severe cystic fibrosis liver disease is associated with novel mechanisms for disease pathogenesis.

BACKGROUND AND AIMS: It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. APPROACH AND RESULTS: Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. CONCLUSION: These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD.

A mutation in SLC37A4 causes a dominantly inherited congenital disorder of glycosylation characterized by liver dysfunction.

SLC37A4 encodes an endoplasmic reticulum (ER)-localized multitransmembrane protein required for transporting glucose-6-phosphate (Glc-6P) into the ER. Once transported into the ER, Glc-6P is subsequently hydrolyzed by tissue-specific phosphatases to glucose and inorganic phosphate during times of glucose depletion. Pathogenic variants in SLC37A4 cause an established recessive disorder known as glycogen storage disorder 1b characterized by liver and kidney dysfunction with neutropenia. We report seven individuals who presented with liver dysfunction multifactorial coagulation deficiency and cardiac issues and were heterozygous for the same variant, c.1267C>T (p.Arg423∗), in SLC37A4; the affected individuals were from four unrelated families. Serum samples from affected individuals showed profound accumulation of both high mannose and hybrid type N-glycans, while N-glycans in fibroblasts and undifferentiated iPSC were normal. Due to the liver-specific nature of this disorder, we generated a CRISPR base-edited hepatoma cell line harboring the c.1267C>T (p.Arg423∗) variant. These cells replicated the secreted abnormalities seen in serum N-glycosylation, and a portion of the mutant protein appears to relocate to a distinct, non-Golgi compartment, possibly ER exit sites. These cells also show a gene dosage-dependent alteration in the Golgi morphology and reduced intraluminal pH that may account for the altered glycosylation. In summary, we identify a recurrent mutation in SLC37A4 that causes a dominantly inherited congenital disorder of glycosylation characterized by coagulopathy and liver dysfunction with abnormal serum N-glycans.

A dyadic approach to the delineation of diagnostic entities in clinical genomics.

The delineation of disease entities is complex, yet recent advances in the molecular characterization of diseases provide opportunities to designate diseases in a biologically valid manner. Here, we have formalized an approach to the delineation of Mendelian genetic disorders that encompasses two distinct but inter-related concepts: (1) the gene that is mutated and (2) the phenotypic descriptor, preferably a recognizably distinct phenotype. We assert that only by a combinatorial or dyadic approach taking both of these attributes into account can a unitary, distinct genetic disorder be designated. We propose that all Mendelian disorders should be designated as "GENE-related phenotype descriptor" (e.g., "CFTR-related cystic fibrosis"). This approach to delineating and naming disorders reconciles the complexity of gene-to-phenotype relationships in a simple and clear manner yet communicates the complexity and nuance of these relationships.

Targeted long-read sequencing identifies missing disease-causing variation.

Despite widespread clinical genetic testing, many individuals with suspected genetic conditions lack a precise diagnosis, limiting their opportunity to take advantage of state-of-the-art treatments. In some cases, testing reveals difficult-to-evaluate structural differences, candidate variants that do not fully explain the phenotype, single pathogenic variants in recessive disorders, or no variants in genes of interest. Thus, there is a need for better tools to identify a precise genetic diagnosis in individuals when conventional testing approaches have been exhausted. We performed targeted long-read sequencing (T-LRS) using adaptive sampling on the Oxford Nanopore platform on 40 individuals, 10 of whom lacked a complete molecular diagnosis. We computationally targeted up to 151 Mbp of sequence per individual and searched for pathogenic substitutions, structural variants, and methylation differences using a single data source. We detected all genomic aberrations-including single-nucleotide variants, copy number changes, repeat expansions, and methylation differences-identified by prior clinical testing. In 8/8 individuals with complex structural rearrangements, T-LRS enabled more precise resolution of the mutation, leading to changes in clinical management in one case. In ten individuals with suspected Mendelian conditions lacking a precise genetic diagnosis, T-LRS identified pathogenic or likely pathogenic variants in six and variants of uncertain significance in two others. T-LRS accurately identifies pathogenic structural variants, resolves complex rearrangements, and identifies Mendelian variants not detected by other technologies. T-LRS represents an efficient and cost-effective strategy to evaluate high-priority genes and regions or complex clinical testing results.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing.

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing, the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare diseases. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery, which should, in turn, increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks such as Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, and researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Bi-allelic variants in MYH3 cause recessively-inherited arthrogryposis.

Arthrogryposis is a clinical feature defined by congenital joint contractures in two or more different body areas which occurs in between 1/3000 and 1/5000 live births. Variants in multiple genes have been associated with distal arthrogryposis syndromes. Heterozygous variants in MYH3 have been identified to cause the dominantly-inherited distal arthrogryposis conditions, Freeman-Sheldon syndrome, Sheldon-Hall syndrome, and multiple pterygium syndrome. In contrast, MYH3 variants underlie both dominantly and recessively inherited Contractures, Pterygia, and Spondylocarpotarsal Fusion syndromes (CPSFS) which are characterized by extensive bony abnormalities in addition to congenital contractures. Here we report two affected sibs with distal arthrogryposis born to unaffected, distantly related parents. Sequencing revealed that both sibs were homozygous for two ultra-rare MYH3 variants, c.3445G>A (p.Glu1149Lys) and c.4760T>C (p.Leu1587Pro). Sequencing and deletion/duplication analysis of 169 other arthrogryposis genes yielded no other compelling candidate variants. This is the first report of biallelic variants in MYH3 being implicated in a distal arthrogryposis phenotype without the additional features of CPSFS. Thus, akin to CPSFS, both dominant and recessively inherited distal arthrogryposis can be caused by variants in MYH3.

The phenotypic and genotypic spectrum of individuals with mono- or biallelic ANK3 variants.

ANK3 encodes ankyrin-G, a protein involved in neuronal development and signaling. Alternative splicing gives rise to three ankyrin-G isoforms comprising different domains with distinct expression patterns. Mono- or biallelic ANK3 variants are associated with non-specific syndromic intellectual disability in 14 individuals (seven with monoallelic and seven with biallelic variants). In this study, we describe the clinical features of 13 additional individuals and review the data on a total of 27 individuals (16 individuals with monoallelic and 11 with biallelic ANK3 variants) and demonstrate that the phenotype for biallelic variants is more severe. The phenotypic features include language delay (92%), autism spectrum disorder (76%), intellectual disability (78%), hypotonia (65%), motor delay (68%), attention deficit disorder (ADD) or attention deficit hyperactivity disorder (ADHD) (57%), sleep disturbances (50%), aggressivity/self-injury (37.5%), and epilepsy (35%). A notable phenotypic difference was presence of ataxia in three individuals with biallelic variants, but in none of the individuals with monoallelic variants. While the majority of the monoallelic variants are predicted to result in a truncated protein, biallelic variants are almost exclusively missense. Moreover, mono- and biallelic variants appear to be localized differently across the three different ankyrin-G isoforms, suggesting isoform-specific pathological mechanisms.

SAMHD1 compound heterozygous rare variants associated with moyamoya and mitral valve disease in the absence of other features of Aicardi-Goutières syndrome.

Aicardi-Goutières syndrome (AGS) is an autosomal recessive inflammatory syndrome that manifests as an early-onset encephalopathy with both neurologic and extraneurologic clinical findings. AGS has been associated with pathogenic variants in nine genes: TREX1, RNASEH2B, RNASEH2C, RNASEH2A, SAMHD1, ADAR, IFIH1, LSM11, and RNU7-1. Diagnosis is established by clinical findings (encephalopathy and acquired microcephaly, intellectual and physical impairments, dystonia, hepatosplenomegaly, sterile pyrexia, and/or chilblains), characteristic abnormalities on cranial CT (calcification of the basal ganglia and white matter) and MRI (leukodystrophic changes), or the identification of pathogenic/likely pathogenic variants in the known genes. One of the genes associated with AGS, SAMHD1, has also been associated with a spectrum of cerebrovascular diseases, including moyamoya disease (MMD). In this report, we describe a 31-year-old male referred to genetics for MMD since childhood who lacked the hallmark features of AGS patients but was found to have compound heterozygous SAMHD1 variants. He later developed mitral valve insufficiency due to recurrent chordal rupture and ultimately underwent a heart transplant at 37 years of age. Thus, these data suggest that SAMHD1 pathogenic variants can cause MMD without typical AGS symptoms and support that SAMHD1 should be assessed in MMD patients even in the absence of AGS features.

Rare variants in ANO1, encoding a calcium-activated chloride channel, predispose to moyamoya disease.

Moyamoya disease, a cerebrovascular disease leading to strokes in children and young adults, is characterized by progressive occlusion of the distal internal carotid arteries and the formation of collateral vessels. Altered genes play a prominent role in the aetiology of moyamoya disease, but a causative gene is not identified in the majority of cases. Exome sequencing data from 151 individuals from 84 unsolved families were analysed to identify further genes for moyamoya disease, then candidate genes assessed in additional cases (150 probands). Two families had the same rare variant in ANO1, which encodes a calcium-activated chloride channel, anoctamin-1. Haplotype analyses found the families were related, and ANO1 p.Met658Val segregated with moyamoya disease in the family with an LOD score of 3.3. Six additional ANO1 rare variants were identified in moyamoya disease families. The ANO1 rare variants were assessed using patch-clamp recordings, and the majority of variants, including ANO1 p.Met658Val, displayed increased sensitivity to intracellular Ca2+. Patients harbouring these gain-of-function ANO1 variants had classic features of moyamoya disease, but also had aneurysm, stenosis and/or occlusion in the posterior circulation. Our studies support that ANO1 gain-of-function pathogenic variants predispose to moyamoya disease and are associated with unique involvement of the posterior circulation.

The Genetic Landscape of Familial Pulmonary Fibrosis.

Rationale and Objectives: Up to 20% of idiopathic interstitial lung disease is familial, referred to as familial pulmonary fibrosis (FPF). An integrated analysis of FPF genetic risk was performed by comprehensively evaluating for genetic rare variants (RVs) in a large cohort of FPF kindreds. Methods: Whole-exome sequencing and/or candidate gene sequencing from affected individuals in 569 FPF kindreds was performed, followed by cosegregation analysis in large kindreds, gene burden analysis, gene-based risk scoring, cell-type enrichment analysis, and coexpression network construction. Measurements and Main Results: It was found that 14.9-23.4% of genetic risk in kindreds could be explained by RVs in genes previously linked to FPF, predominantly telomere-related genes. New candidate genes were identified in a small number of families-including SYDE1, SERPINB8, GPR87, and NETO1-and tools were developed for evaluation and prioritization of RV-containing genes across kindreds. Several pathways were enriched for RV-containing genes in FPF, including focal adhesion and mitochondrial complex I assembly. By combining single-cell transcriptomics with prioritized candidate genes, expression of RV-containing genes was discovered to be enriched in smooth muscle cells, type II alveolar epithelial cells, and endothelial cells. Conclusions: In the most comprehensive FPF genetic study to date, the prevalence of RVs in known FPF-related genes was defined, and new candidate genes and pathways relevant to FPF were identified. However, new RV-containing genes shared across multiple kindreds were not identified, thereby suggesting that heterogeneous genetic variants involving a variety of genes and pathways mediate genetic risk in most FPF kindreds.

Genetic Modifiers of Cystic Fibrosis Lung Disease Severity: Whole-Genome Analysis of 7,840 Patients.

Rationale: Lung disease is the major cause of morbidity and mortality in persons with cystic fibrosis (pwCF). Variability in CF lung disease has substantial non-CFTR (CF transmembrane conductance regulator) genetic influence. Identification of genetic modifiers has prognostic and therapeutic importance. Objectives: Identify genetic modifier loci and genes/pathways associated with pulmonary disease severity. Methods: Whole-genome sequencing data on 4,248 unique pwCF with pancreatic insufficiency and lung function measures were combined with imputed genotypes from an additional 3,592 patients with pancreatic insufficiency from the United States, Canada, and France. This report describes association of approximately 15.9 million SNPs using the quantitative Kulich normal residual mortality-adjusted (KNoRMA) lung disease phenotype in 7,840 pwCF using premodulator lung function data. Measurements and Main Results: Testing included common and rare SNPs, transcriptome-wide association, gene-level, and pathway analyses. Pathway analyses identified novel associations with genes that have key roles in organ development, and we hypothesize that these genes may relate to dysanapsis and/or variability in lung repair. Results confirmed and extended previous genome-wide association study findings. These whole-genome sequencing data provide finely mapped genetic information to support mechanistic studies. No novel primary associations with common single variants or rare variants were found. Multilocus effects at chr5p13 (SLC9A3/CEP72) and chr11p13 (EHF/APIP) were identified. Variant effect size estimates at associated loci were consistently ordered across the cohorts, indicating possible age or birth cohort effects. Conclusions: This premodulator genomic, transcriptomic, and pathway association study of 7,840 pwCF will facilitate mechanistic and postmodulator genetic studies and the development of novel therapeutics for CF lung disease.

Multiplexed Functional Assessment of Genetic Variants in CARD11.

Genetic testing has increased the number of variants identified in disease genes, but the diagnostic utility is limited by lack of understanding variant function. CARD11 encodes an adaptor protein that expresses dominant-negative and gain-of-function variants associated with distinct immunodeficiencies. Here, we used a "cloning-free" saturation genome editing approach in a diploid cell line to simultaneously score 2,542 variants for decreased or increased function in the region of CARD11 associated with immunodeficiency. We also described an exon-skipping mechanism for CARD11 dominant-negative activity. The classification of reported clinical variants was sensitive (94.6%) and specific (88.9%), which rendered the data immediately useful for interpretation of seven coding and splicing variants implicated in immunodeficiency found in our clinic. This approach is generalizable for variant interpretation in many other clinically actionable genes, in any relevant cell type.

Mutations in MYLPF Cause a Novel Segmental Amyoplasia that Manifests as Distal Arthrogryposis.

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.