A deep catalogue of protein-coding variation in 983,578 individuals.

Rare coding variants that substantially affect function provide insights into the biology of a gene1-3. However, ascertaining the frequency of such variants requires large sample sizes4-8. Here we present a catalogue of human protein-coding variation, derived from exome sequencing of 983,578 individuals across diverse populations. In total, 23% of the Regeneron Genetics Center Million Exome (RGC-ME) data come from individuals of African, East Asian, Indigenous American, Middle Eastern and South Asian ancestry. The catalogue includes more than 10.4 million missense and 1.1 million predicted loss-of-function (pLOF) variants. We identify individuals with rare biallelic pLOF variants in 4,848 genes, 1,751 of which have not been previously reported. From precise quantitative estimates of selection against heterozygous loss of function (LOF), we identify 3,988 LOF-intolerant genes, including 86 that were previously assessed as tolerant and 1,153 that lack established disease annotation. We also define regions of missense depletion at high resolution. Notably, 1,482 genes have regions that are depleted of missense variants despite being tolerant of pLOF variants. Finally, we estimate that 3% of individuals have a clinically actionable genetic variant, and that 11,773 variants reported in ClinVar with unknown significance are likely to be deleterious cryptic splice sites. To facilitate variant interpretation and genetics-informed precision medicine, we make this resource of coding variation from the RGC-ME dataset publicly accessible through a variant allele frequency browser.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Exome sequencing and characterization of 49,960 individuals in the UK Biobank.

The UK Biobank is a prospective study of 502,543 individuals, combining extensive phenotypic and genotypic data with streamlined access for researchers around the world1. Here we describe the release of exome-sequence data for the first 49,960 study participants, revealing approximately 4 million coding variants (of which around 98.6% have a frequency of less than 1%). The data include 198,269 autosomal predicted loss-of-function (LOF) variants, a more than 14-fold increase compared to the imputed sequence. Nearly all genes (more than 97%) had at least one carrier with a LOF variant, and most genes (more than 69%) had at least ten carriers with a LOF variant. We illustrate the power of characterizing LOF variants in this population through association analyses across 1,730 phenotypes. In addition to replicating established associations, we found novel LOF variants with large effects on disease traits, including PIEZO1 on varicose veins, COL6A1 on corneal resistance, MEPE on bone density, and IQGAP2 and GMPR on blood cell traits. We further demonstrate the value of exome sequencing by surveying the prevalence of pathogenic variants of clinical importance, and show that 2% of this population has a medically actionable variant. Furthermore, we characterize the penetrance of cancer in carriers of pathogenic BRCA1 and BRCA2 variants. Exome sequences from the first 49,960 participants highlight the promise of genome sequencing in large population-based studies and are now accessible to the scientific community.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

In-Depth Mass Spectrometry Analysis Reveals the Plasma Proteomic and N-Glycoproteomic Impact of an Amish-Enriched Cardioprotective Variant in B4GALT1.

B4GALT1 encodes ß-1,4-galactosyltransferase 1, an enzyme that plays a major role in glycan synthesis in the Golgi apparatus by catalyzing the addition of terminal galactose. Studies increasingly suggest that B4GALT1 may be involved in the regulation of lipid metabolism pathways. Recently, we discovered a single-site missense variant Asn352Ser (N352S) in the functional domain of B4GALT1 in an Amish population, which decreases the level of LDL-cholesterol (LDL-c) as well as the protein levels of ApoB, fibrinogen, and IgG in the blood. To systematically evaluate the effects of this missense variant on protein glycosylation, expression, and secretion, we developed a nano-LC-MS/MS-based platform combined with TMT-labeling for in-depth quantitative proteomic and glycoproteomic analyses in the plasma of individuals homozygous for the B4GALT1 missense variant N352S versus non-carriers (n = 5 per genotype). A total of 488 secreted proteins in the plasma were identified and quantified, 34 of which showed significant fold changes in protein levels between N352S homozygotes and non-carriers. We determined N-glycosylation profiles from 370 glycosylation sites in 151 glycoproteins and identified ten proteins most significantly associated with decreased galactosylation and sialyation in B4GALT1 N352S homozygotes. These results further support that B4GALT1 N352S alters the glycosylation profiles of a variety of critical target proteins, thus governing the functions of these proteins in multiple pathways, such as those involved in lipid metabolism, coagulation, and the immune response.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Missense variants in SORT1 are associated with LDL-C in an Amish population.

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

Associations of genome-wide and regional autozygosity with 96 complex traits in old order Amish.

BACKGROUND: Autozygosity, the proportion of the genome that is homozygous by descent, has been associated with variation in multiple health-related traits impacting evolutionary fitness. Autozygosity (FROH) is typically measured from runs of homozygosity (ROHs) that arise when identical-by-descent (IBD) haplotypes are inherited from each parent. Population isolates with a small set of common founders have elevated autozygosity relative to outbred populations. METHODS: In this study, we examined whether degree of autozygosity was associated with variation in 96 cardiometabolic traits among 7221 Old Order Amish individuals residing in Lancaster County, PA. We estimated the average length of an ROH segment to be 6350 KB, with each individual having on average 17.2 segments 1.5 KB or larger. Measurements of genome-wide and regional FROH were used as the primary predictors of trait variation in association analysis. RESULTS: In genome-wide FROH analysis, we did not identify any associations that withstood Bonferroni-correction (p = 0.0005). However, on regional FROH analysis, we identified associations exceeding genome-wide thresholds for two traits: serum bilirubin levels, which were significantly associated with a region on chromosome 2 localized to a region surrounding UGT1A10 (p = 1 × 10- 43), and HbA1c levels, which were significantly associated with a region on chromosome 8 localized near CHRNB3 (p = 8 × 10- 10). CONCLUSIONS: These analyses highlight the potential value of autozygosity mapping in founder populations.

Characterization of Exome Variants and Their Metabolic Impact in 6,716 American Indians from the Southwest US.

Applying exome sequencing to populations with unique genetic architecture has the potential to reveal novel genes and variants associated with traits and diseases. We sequenced and analyzed the exomes of 6,716 individuals from a Southwestern American Indian (SWAI) population with well-characterized metabolic traits. We found that the SWAI population has distinct allelic architecture compared to populations of European and East Asian ancestry, and there were many predicted loss-of-function (pLOF) and nonsynonymous variants that were highly enriched or private in the SWAI population. We used pLOF and nonsynonymous variants in the SWAI population to evaluate gene-burden associations of candidate genes from European genome-wide association studies (GWASs) for type 2 diabetes, body mass index, and four major plasma lipids. We found 19 significant gene-burden associations for 11 genes, providing additional evidence for prioritizing candidate effector genes of GWAS signals. Interestingly, these associations were mainly driven by pLOF and nonsynonymous variants that are unique or highly enriched in the SWAI population. Particularly, we found four pLOF or nonsynonymous variants in APOB, APOE, PCSK9, and TM6SF2 that are private or enriched in the SWAI population and associated with low-density lipoprotein (LDL) cholesterol levels. Their large estimated effects on LDL cholesterol levels suggest strong impacts on protein function and potential clinical implications of these variants in cardiovascular health. In summary, our study illustrates the utility and potential of exome sequencing in genetically unique populations, such as the SWAI population, to prioritize candidate effector genes within GWAS loci and to find additional variants in known disease genes with potential clinical impact.

Pharmacogenetics of sodium-glucose co-transporter-2 inhibitors: Validation of a sex-agnostic pharmacodynamic biomarker.

AIM: To validate pharmacodynamic responses to sodium-glucose co-transporter-2 (SGLT2) inhibitors and test for association with genetic variants in SLC5A4, SLC5A9, and SLC2A9. METHODS: Canagliflozin (300 mg), a SGLT2 inhibitor, was administered to 30 healthy volunteers. Several endpoints were measured to assess clinically relevant responses, including drug-induced increases in urinary excretion of glucose, sodium and uric acid. RESULTS: This pilot study confirmed that canagliflozin (300 mg) triggered acute changes in mean levels of several biomarkers: fasting plasma glucose (-4.1 mg/dL; P = 6 × 10-5 ), serum creatinine (+0.05 mg/dL; P = 8 × 10-4 ) and serum uric acid (-0.90 mg/dL; P = 5 × 10-10 ). The effects of sex on glucosuria depended upon how data were normalized. Whereas males' responses were ~60% greater when data were normalized to body surface area, males and females exhibited similar responses when glucosuria was expressed as grams of urinary glucose per gram-creatinine. The magnitude of glucosuria was not significantly correlated with fasting plasma glucose, estimated glomerular filtration rate or age in those healthy individuals without diabetes with an estimated glomerular filtration rate of more than 60 mL/min/1.73m2 . CONCLUSIONS: Normalizing data relative to creatinine excretion will facilitate including data from males and females in a single analysis. Furthermore, because our ongoing pharmacogenomic study (NCT02891954) is conducted in healthy individuals, this will facilitate detection of genetic associations with limited confounding by other factors such as HbA1c and renal function.

Acute pharmacodynamic responses to exenatide: Drug-induced increases in insulin secretion and glucose effectiveness.

AIM: Glucagon-like peptide-1 receptor agonists provide multiple benefits to patients with type 2 diabetes, including improved glycaemic control, weight loss and decreased risk of major adverse cardiovascular events. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. METHODS: Exenatide (5 µg, subcutaneously) or saline (0.2 ml, subcutaneously) was administered to 62 healthy volunteers. Frequently sampled intravenous glucose tolerance tests were conducted to assess the impact of exenatide on insulin secretion and insulin action. This pilot study was a crossover design in which participants received exenatide and saline in random order. RESULTS: Exenatide increased first phase insulin secretion 1.9-fold (p = 1.9 × 10-9 ) and accelerated the rate of glucose disappearance 2.4-fold (p = 2 × 10-10 ). Minimal model analysis showed that exenatide increased glucose effectiveness (Sg ) by 32% (p = .0008) but did not significantly affect insulin sensitivity (Si ). The exenatide-induced increase in insulin secretion made the largest contribution to interindividual variation in exenatide-induced acceleration of glucose disappearance while interindividual variation in the drug effect on Sg contributed to a lesser extent (ß = 0.58 or 0.27, respectively). CONCLUSIONS: This pilot study provides validation for the value of a frequently sampled intravenous glucose tolerance test (including minimal model analysis) to provide primary data for our ongoing pharmacogenomic study of pharmacodynamic effects of semaglutide (NCT05071898). Three endpoints provide quantitative assessments of the effects of glucagon-like peptide-1 receptor agonists on glucose metabolism: first phase insulin secretion, glucose disappearance rates and glucose effectiveness.

Genome-wide association analysis of serum alanine and aspartate aminotransferase, and the modifying effects of BMI in 388k European individuals.

Serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) are biomarkers for liver health. Here we report the largest genome-wide association analysis to date of serum ALT and AST levels in over 388k people of European ancestry from UK biobank and DiscovEHR. Eleven million imputed markers with a minor allele frequency (MAF) ≥ 0.5% were analyzed. Overall, 300 ALT and 336 AST independent genome-wide significant associations were identified. Among them, 81 ALT and 61 AST associations are reported for the first time. Genome-wide interaction study identified 9 ALT and 12 AST independent associations significantly modified by body mass index (BMI), including several previously reported potential liver disease therapeutic targets, for example, PNPLA3, HSD17B13, and MARC1. While further work is necessary to understand the effect of ALT and AST-associated variants on liver disease, the weighted burden of significant BMI-modified signals is significantly associated with liver disease outcomes. In summary, this study identifies genetic associations which offer an important step forward in understanding the genetic architecture of serum ALT and AST levels. Significant interactions between BMI and genetic loci not only highlight the important role of adiposity in liver damage but also shed light on the genetic etiology of liver disease in obese individuals.

Cerebral blood flow and cardiovascular risk effects on resting brain regional homogeneity.

Regional homogeneity (ReHo) is a measure of local functional brain connectivity that has been reported to be altered in a wide range of neuropsychiatric disorders. Computed from brain resting-state functional MRI time series, ReHo is also sensitive to fluctuations in cerebral blood flow (CBF) that in turn may be influenced by cerebrovascular health. We accessed cerebrovascular health with Framingham cardiovascular risk score (FCVRS). We hypothesize that ReHo signal may be influenced by regional CBF; and that these associations can be summarized as FCVRS→CBF→ReHo. We used three independent samples to test this hypothesis. A test-retest sample of N = 30 healthy volunteers was used for test-retest evaluation of CBF effects on ReHo. Amish Connectome Project (ACP) sample (N = 204, healthy individuals) was used to evaluate association between FCVRS and ReHo and testing if the association diminishes given CBF. The UKBB sample (N = 6,285, healthy participants) was used to replicate the effects of FCVRS on ReHo. We observed strong CBF→ReHo links (p<2.5 × 10-3) using a three-point longitudinal sample. In ACP sample, marginal and partial correlations analyses demonstrated that both CBF and FCVRS were significantly correlated with the whole-brain average (p<10-6) and regional ReHo values, with the strongest correlations observed in frontal, parietal, and temporal areas. Yet, the association between ReHo and FCVRS became insignificant once the effect of CBF was accounted for. In contrast, CBF→ReHo remained significantly linked after adjusting for FCVRS and demographic covariates (p<10-6). Analysis in N = 6,285 replicated the FCVRS→ReHo effect (p = 2.7 × 10-27). In summary, ReHo alterations in health and neuropsychiatric illnesses may be partially driven by region-specific variability in CBF, which is, in turn, influenced by cardiovascular factors.

UGT1A1 genetic variants are associated with increases in bilirubin levels in rheumatoid arthritis patients treated with sarilumab.

Sarilumab is a human monoclonal antibody against interleukin (IL)-6Rα that has been approved for the treatment of adult patients with moderately to severely active rheumatoid arthritis (RA) and an inadequate response or intolerance to one or more disease-modifying antirheumatic drugs (DMARDs). Mild liver function test abnormalities have been observed in patients treated with sarilumab. We describe a genome-wide association study of bilirubin elevations in RA patients treated with sarilumab. Array genotyping and exome sequencing were performed on DNA samples from 1075 patients. Variants in the UGT1A1 gene were strongly associated with maximum bilirubin elevations in sarilumab-treated patients (rs4148325; p = 2.88 × 10-41) but were not associated with aminotransferase elevations. No other independent loci showed evidence of association with bilirubin elevations after sarilumab treatment. These findings suggest that most bilirubin increases during sarilumab treatment are related to genetic variation in UGT1A1 rather than underlying liver injury.

A missense variant Arg611Cys in LIPE which encodes hormone sensitive lipase decreases lipolysis and increases risk of type 2 diabetes in American Indians.

AIMS: Hormone sensitive lipase (HSL), encoded by the LIPE gene, is involved in lipolysis. Based on prior animal and human studies, LIPE was analysed as a candidate gene for the development of type 2 diabetes (T2D) in a community-based sample of American Indians. MATERIALS AND METHODS: Whole-exome sequence data from 6782 participants with longitudinal clinical measures were used to identify variation in LIPE. RESULTS: Amongst the 16 missense variants identified, an Arg611Cys variant (rs34052647; Cys-allele frequency = 0.087) significantly associated with T2D (OR [95% CI] = 1.38 [1.17-1.64], p = 0.0002, adjusted for age, sex, birth year, and the first five genetic principal components) and an earlier onset age of T2D (HR = 1.22 [1.09-1.36], p = 0.0005). This variant was further analysed for quantitative traits related to T2D. Amongst non-diabetic American Indians, those with the T2D risk Cys-allele had increased insulin levels during an oral glucose tolerance test (0.07 SD per Cys-allele, p = 0.04) and a mixed meal test (0.08 log10 µU/ml per Cys-allele, p = 0.003), and had increased lipid oxidation rates post-absorptively and during insulin infusion (0.07 mg [kg estimated metabolic body size {EMBS}]-1  min-1 per Cys-allele for both, p = 0.01 and 0.009, respectively), compared to individuals with the non-risk Arg-allele. In vitro functional studies showed that cells expressing the Cys-allele had a 17.2% decrease in lipolysis under isoproterenol stimulation (p = 0.03) and a 21.3% decrease in lipase enzyme activity measured by using p-nitrophenyl butyrate as a substrate (p = 0.04) compared to the Arg-allele. CONCLUSION: The Arg611Cys variant causes a modest impairment in lipolysis, thereby affecting glucose homoeostasis and risk of T2D.

Impact of parental relatedness on reproductive outcomes among the Old Order Amish of Lancaster County.

Genetically isolated populations that arise due to recent bottleneck events have reduced genetic variation reflecting the common set of founders. Increased genetic relatedness among members of isolated populations puts them at increased risk for some recessive disorders that are rare in outbred populations. To assess the burden on reproductive health, we compared frequencies of adverse reproductive outcomes between Amish couples who were both heterozygous carriers of a highly penetrant pathogenic or likely pathogenic variant and noncarrier couples from the same Amish community. In addition, we evaluated whether overall genetic relatedness of parents was associated with reproductive outcomes. Of the 1824 couples included in our study, 11.1% were at risk of producing a child with an autosomal recessive disorder. Carrier couples reported a lower number of miscarriages compared to noncarrier couples (p = 0.02), although the number of stillbirths (p = 0.3), live births (p = 0.9), and number of pregnancies (p = 0.9) did not differ significantly between groups. In contrast, higher overall relatedness between spouses was positively correlated with number of live births (p < 0.0001), pregnancies (p < 0.0001), and stillbirths (p = 0.03), although not with the number of miscarriages (p = 0.4). These results highlight a complex association between relatedness of parents and reproductive health outcomes in this community.

Biallelic truncating variants in the muscular A-type lamin-interacting protein (MLIP) gene cause myopathy with hyperCKemia.

BACKGROUND AND PURPOSE: Muscular A-type lamin-interacting protein (MLIP) is most abundantly expressed in cardiac and skeletal muscle. In vitro and animal studies have shown its regulatory role in myoblast differentiation and in organization of myonuclear positioning in skeletal muscle, as well as in cardiomyocyte adaptation and cardiomyopathy. We report the association of biallelic truncating variation in the MLIP gene with human disease in five individuals from two unrelated pedigrees. METHODS: Clinical evaluation and exome sequencing were performed in two unrelated families with elevated creatine kinase level. RESULTS: Family 1. A 6-year-old girl born to consanguineous parents of Arab-Muslim origin presented with myalgia, early fatigue after mild-to-moderate physical exertion, and elevated creatine kinase levels up to 16,000 U/L. Exome sequencing revealed a novel homozygous nonsense variant, c.2530C>T; p.Arg844Ter, in the MLIP gene. Family 2. Three individuals from two distantly related families of Old Order Amish ancestry presented with elevated creatine kinase levels, one of whom also presented with abnormal electrocardiography results. On exome sequencing, all showed homozygosity for a novel nonsense MLIP variant c.1825A>T; p.Lys609Ter. Another individual from this pedigree, who had sinus arrhythmia and for whom creatine kinase level was not available, was also homozygous for this variant. CONCLUSIONS: Our findings suggest that biallelic truncating variants in MLIP result in myopathy characterized by hyperCKemia. Moreover, these cases of MLIP-related disease may indicate that at least in some instances this condition is associated with muscle decompensation and fatigability during low-to-moderate intensity muscle exertion as well as possible cardiac involvement.