Transcriptional profiling in microglia across physiological and pathological states identifies a transcriptional module associated with neurodegeneration.

Microglia are the resident immune cells of the central nervous system and are involved in brain development, homeostasis, and disease. New imaging and genomics technologies are revealing microglial complexity across developmental and functional states, brain regions, and diseases. We curated a set of publicly available gene expression datasets from human microglia spanning disease and health to identify sets of genes reflecting physiological and pathological microglial states. We also integrated multiple human microglial single-cell RNA-seq datasets in Alzheimer's disease (AD), multiple sclerosis (MS), and Parkinson's disease, and identified a distinct microglial transcriptional signature shared across diseases. Analysis of germ-line DNA identified genes with variants associated with AD and MS that are overrepresented in microglial gene sets, including the disease-associated transcriptional signature. This work points to genes that are dysregulated in disease states and provides a resource for the analysis of diseases in which microglia are implicated by genetic evidence.

NOTCH3 p.Arg1231Cys is markedly enriched in South Asians and associated with stroke.

The genetic factors of stroke in South Asians are largely unexplored. Exome-wide sequencing and association analysis (ExWAS) in 75 K Pakistanis identified NM_000435.3(NOTCH3):c.3691 C > T, encoding the missense amino acid substitution p.Arg1231Cys, enriched in South Asians (alternate allele frequency = 0.58% compared to 0.019% in Western Europeans), and associated with subcortical hemorrhagic stroke [odds ratio (OR) = 3.39, 95% confidence interval (CI) = [2.26, 5.10], p = 3.87 × 10-9), and all strokes (OR [CI] = 2.30 [1.77, 3.01], p = 7.79 × 10-10). NOTCH3 p.Arg231Cys was strongly associated with white matter hyperintensity on MRI in United Kingdom Biobank (UKB) participants (effect [95% CI] in SD units = 1.1 [0.61, 1.5], p = 3.0 × 10-6). The variant is attributable for approximately 2.0% of hemorrhagic strokes and 1.1% of all strokes in South Asians. These findings highlight the value of diversity in genetic studies and have major implications for genomic medicine and therapeutic development in South Asian populations.

Identification and functional validation of rare coding variants in genes linked to monogenic obesity.

OBJECTIVE: Rare cases of monogenic obesity, which may respond to specific therapeutics, can remain undetected in populations in which polygenic obesity is prevalent. This study examined rare DNA variation in established monogenic obesity genes within a community using whole-exome sequence data from 6803 longitudinally studied individuals. METHODS: Exome data across 15 monogenic obesity genes were analyzed for nonsynonymous variants observed in any child with a maximum BMI z score > 2 (N = 279) but not observed in a child with a maximum BMI z score ≤ 0 (n = 1542) or that occurred in adults in the top 5th percentile of BMI (n = 263) but not in adults below the median BMI (n = 2629). Variants were then functionally analyzed using luciferase assays. RESULTS: The comparisons between cases of obesity and controls identified eight missense variants in six genes: DYRK1B, KSR2, MC4R, NTRK2, PCSK1, and SIM1. Among these, MC4R p.A303P and p.R165G were previously shown to impair MC4R function. Functional analyses of the remaining six variants suggest that KSR2 p.I402F and p.T193I and NTRK2 p.S249Y alter protein function. CONCLUSIONS: In addition to MC4R, rare missense variants in KSR2 and NTRK2 may potentially explain the severe obesity observed for the carriers.

Genetic risk factors for COVID-19 and influenza are largely distinct.

Coronavirus disease 2019 (COVID-19) and influenza are respiratory illnesses caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza viruses, respectively. Both diseases share symptoms and clinical risk factors1, but the extent to which these conditions have a common genetic etiology is unknown. This is partly because host genetic risk factors are well characterized for COVID-19 but not for influenza, with the largest published genome-wide association studies for these conditions including >2 million individuals2 and about 1,000 individuals3-6, respectively. Shared genetic risk factors could point to targets to prevent or treat both infections. Through a genetic study of 18,334 cases with a positive test for influenza and 276,295 controls, we show that published COVID-19 risk variants are not associated with influenza. Furthermore, we discovered and replicated an association between influenza infection and noncoding variants in B3GALT5 and ST6GAL1, neither of which was associated with COVID-19. In vitro small interfering RNA knockdown of ST6GAL1-an enzyme that adds sialic acid to the cell surface, which is used for viral entry-reduced influenza infectivity by 57%. These results mirror the observation that variants that downregulate ACE2, the SARS-CoV-2 receptor, protect against COVID-19 (ref. 7). Collectively, these findings highlight downregulation of key cell surface receptors used for viral entry as treatment opportunities to prevent COVID-19 and influenza.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

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.

Acute pharmacodynamic responses to sitagliptin: Drug-induced increase in early insulin secretion in oral glucose tolerance test.

DPP4 inhibitors are widely prescribed as treatments for type 2 diabetes. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Sitagliptin (100 mg) was administered to 47 healthy volunteers. Several endpoints were measured to assess clinically relevant responses - including the effect of sitagliptin on glucose and insulin levels during an oral glucose tolerance test (OGTT). This pilot study confirmed that sitagliptin (100 mg) decreased the area under the curve for glucose during an OGTT (p = 0.0003). Furthermore, sitagliptin promoted insulin secretion during the early portion of the OGTT as reflected by an increase in the ratio of plasma insulin at 30 min divided by plasma insulin at 60 min (T30:T60) from mean ± SEM 0.87 ± 0.05 to 1.62 ± 0.36 mU/L (p = 0.04). The magnitude of sitagliptin's effect on insulin secretion (as judged by the increase in the T30:T60 ratio for insulin) was correlated with the magnitude of sitagliptin-induced increase in the area under the curve for intact plasma GLP1 levels during the first hour of the OGTT. This study confirmed previously reported sex differences in glucose and insulin levels during an OGTT. Specifically, females exhibited higher levels of glucose and insulin at the 90-180 min time points. However, we did not detect significant sex-associated differences in the magnitude of sitagliptin-induced changes in T30:T60 ratios for either glucose or insulin. In conclusion, T30:T60 ratios for insulin and glucose during an OGTT provide useful indices to assess pharmacodynamic responses to DPP4 inhibitors.

Massive Parallel DNA Sequencing of Patients with Inherited Cardiomyopathies in Cyprus and Suggestion of Digenic or Oligogenic Inheritance.

Inherited cardiomyopathies represent a highly heterogeneous group of cardiac diseases. DNA variants in genes expressed in cardiomyocytes cause a diverse spectrum of cardiomyopathies, ultimately leading to heart failure, arrythmias, and sudden cardiac death. We applied massive parallel DNA sequencing using a 72-gene panel for studying inherited cardiomyopathies. We report on variants in 25 families, where pathogenicity was predicted by different computational approaches, databases, and an in-house filtering analysis. All variants were validated using Sanger sequencing. Familial segregation was tested when possible. We identified 41 different variants in 26 genes. Analytically, we identified fifteen variants previously reported in the Human Gene Mutation Database: twelve mentioned as disease-causing mutations (DM) and three as probable disease-causing mutations (DM?). Additionally, we identified 26 novel variants. We classified the forty-one variants as follows: twenty-eight (68.3%) as variants of uncertain significance, eight (19.5%) as likely pathogenic, and five (12.2%) as pathogenic. We genetically characterized families with a cardiac phenotype. The genetic heterogeneity and the multiplicity of candidate variants are making a definite molecular diagnosis challenging, especially when there is a suspicion of incomplete penetrance or digenic-oligogenic inheritance. This is the first systematic study of inherited cardiac conditions in Cyprus, enabling us to develop a genetic baseline and precision cardiology.

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.

Acute pharmacodynamic responses to sitagliptin: Drug-induced increase in early insulin secretion in oral glucose tolerance test.

Aim: DPP4 inhibitors are widely prescribed as treatments for type 2 diabetes. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. Methods: Sitagliptin (100 mg) was administered to 47 healthy volunteers. Several endpoints were measured to assess clinically relevant responses - including the effect of sitagliptin on glucose and insulin levels during an oral glucose tolerance test (OGTT). Results: This pilot study confirmed that sitagliptin (100 mg) decreased the area under the curve for glucose during an OGTT (p=0.0003). Furthermore, sitagliptin promoted insulin secretion during the early portion of the OGTT as reflected by an increase in the ratio of plasma insulin at 30 min divided by plasma insulin at 60 min (T30:T60) from 0.87+/-0.05 to 1.62+/-0.36 mU/L (p=0.04). The magnitude of sitagliptin's effect on insulin secretion (as judged by the increase in the T30:T60 ratio for insulin) was correlated with the magnitude of sitagliptin-induced increase in the area under the curve for intact plasma GLP1 levels during the first hour of the OGTT. This study confirmed previously reported sex differences in glucose and insulin levels during an OGTT. Specifically, females exhibited higher levels of glucose and insulin at the 90-180 min time points. However, we did not detect significant sex-associated differences in the magnitude of sitagliptin-induced changes in T30:T60 ratios for either glucose or insulin. Conclusions: T30:T60 ratios for insulin and glucose during an OGTT provide useful indices to assess pharmacodynamic responses to DPP4 inhibitors.

Glucocorticoid Receptor Gene (NR3C1) Polymorphisms and Metabolic Syndrome: Insights from the Mennonite Population.

The regulation of the hypothalamic-pituitary-adrenal (HPA) axis is associated with polymorphisms and the methylation degree of the glucocorticoid receptor gene (NR3C1) and is potentially involved in the development of metabolic syndrome (MetS). In order to evaluate the association between MetS with the polymorphisms, methylation, and gene expression of the NR3C1 in the genetically isolated Brazilian Mennonite population, we genotyped 20 NR3C1 polymorphisms in 74 affected (MetS) and 138 unaffected individuals without affected first-degree relatives (Co), using exome sequencing, as well as five variants from non-exonic regions, in 70 MetS and 166 Co, using mass spectrometry. The methylation levels of 11 1F CpG sites were quantified using pyrosequencing (66 MetS and 141 Co), and the NR3C1 expression was evaluated via RT-qPCR (14 MetS and 25 Co). Age, physical activity, and family environment during childhood were associated with MetS. Susceptibility to MetS, independent of these factors, was associated with homozygosity for rs10482605*C (OR = 4.74, pcorr = 0.024) and the haplotype containing TTCGTTGATT (rs3806855*T_ rs3806854*T_rs10482605*C_rs10482614*G_rs6188*T_rs258813*T_rs33944801*G_rs34176759*A_rs17209258*T_rs6196*T, OR = 4.74, pcorr = 0.048), as well as for the CCT haplotype (rs41423247*C_ rs6877893*C_rs258763*T), OR = 6.02, pcorr = 0.030), but not to the differences in methylation or gene expression. Thus, NR3C1 polymorphisms seem to modulate the susceptibility to MetS in Mennonites, independently of lifestyle and early childhood events, and their role seems to be unrelated to DNA methylation and gene expression.

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.

Hemizygous variants in protein phosphatase 1 regulatory subunit 3F (PPP1R3F) are associated with a neurodevelopmental disorder characterized by developmental delay, intellectual disability and autistic features.

Protein phosphatase 1 regulatory subunit 3F (PPP1R3F) is a member of the glycogen targeting subunits (GTSs), which belong to the large group of regulatory subunits of protein phosphatase 1 (PP1), a major eukaryotic serine/threonine protein phosphatase that regulates diverse cellular processes. Here, we describe the identification of hemizygous variants in PPP1R3F associated with a novel X-linked recessive neurodevelopmental disorder in 13 unrelated individuals. This disorder is characterized by developmental delay, mild intellectual disability, neurobehavioral issues such as autism spectrum disorder, seizures and other neurological findings including tone, gait and cerebellar abnormalities. PPP1R3F variants segregated with disease in affected hemizygous males that inherited the variants from their heterozygous carrier mothers. We show that PPP1R3F is predominantly expressed in brain astrocytes and localizes to the endoplasmic reticulum in cells. Glycogen content in PPP1R3F knockout astrocytoma cells appears to be more sensitive to fluxes in extracellular glucose levels than in wild-type cells, suggesting that PPP1R3F functions in maintaining steady brain glycogen levels under changing glucose conditions. We performed functional studies on nine of the identified variants and observed defects in PP1 binding, protein stability, subcellular localization and regulation of glycogen metabolism in most of them. Collectively, the genetic and molecular data indicate that deleterious variants in PPP1R3F are associated with a new X-linked disorder of glycogen metabolism, highlighting the critical role of GTSs in neurological development. This research expands our understanding of neurodevelopmental disorders and the role of PP1 in brain development and proper function.

Rare coding variants in CHRNB2 reduce the likelihood of smoking.

Human genetic studies of smoking behavior have been thus far largely limited to common variants. Studying rare coding variants has the potential to identify drug targets. We performed an exome-wide association study of smoking phenotypes in up to 749,459 individuals and discovered a protective association in CHRNB2, encoding the ß2 subunit of the α4ß2 nicotine acetylcholine receptor. Rare predicted loss-of-function and likely deleterious missense variants in CHRNB2 in aggregate were associated with a 35% decreased odds for smoking heavily (odds ratio (OR) = 0.65, confidence interval (CI) = 0.56-0.76, P = 1.9 × 10-8). An independent common variant association in the protective direction ( rs2072659 ; OR = 0.96; CI = 0.94-0.98; P = 5.3 × 10-6) was also evident, suggesting an allelic series. Our findings in humans align with decades-old experimental observations in mice that ß2 loss abolishes nicotine-mediated neuronal responses and attenuates nicotine self-administration. Our genetic discovery will inspire future drug designs targeting CHRNB2 in the brain for the treatment of nicotine addiction.

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.

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.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

A deep catalog of protein-coding variation in 985,830 individuals.

Coding variants that have significant impact on function can provide insights into the biology of a gene but are typically rare in the population. Identifying and ascertaining the frequency of such rare variants requires very large sample sizes. Here, we present the largest catalog of human protein-coding variation to date, derived from exome sequencing of 985,830 individuals of diverse ancestry to serve as a rich resource for studying rare coding variants. Individuals of African, Admixed American, East Asian, Middle Eastern, and South Asian ancestry account for 20% of this Exome dataset. Our catalog of variants includes approximately 10.5 million missense (54% novel) and 1.1 million predicted loss-of-function (pLOF) variants (65% novel, 53% observed only once). We identified individuals with rare homozygous pLOF variants in 4,874 genes, and for 1,838 of these this work is the first to document at least one pLOF homozygote. Additional insights from the RGC-ME dataset include 1) improved estimates of selection against heterozygous loss-of-function and identification of 3,459 genes intolerant to loss-of-function, 83 of which were previously assessed as tolerant to loss-of-function and 1,241 that lack disease annotations; 2) identification of regions depleted of missense variation in 457 genes that are tolerant to loss-of-function; 3) functional interpretation for 10,708 variants of unknown or conflicting significance reported in ClinVar as cryptic splice sites using splicing score thresholds based on empirical variant deleteriousness scores derived from RGC-ME; and 4) an observation that approximately 3% of sequenced individuals carry a clinically actionable genetic variant in the ACMG SF 3.1 list of genes. We make this important resource of coding variation available to the public through a variant allele frequency browser. We anticipate that this report and the RGC-ME dataset will serve as a valuable reference for understanding rare coding variation and help advance precision medicine efforts.

Association of parity with body mass index and cardiometabolic risk in high-parous women.

OBJECTIVE: Postpregnancy weight retention contributes to obesity, but the long-term effect of parity on body mass index (BMI) and other cardiometabolic risk factors is unclear. We aimed to evaluate the relationship between parity and BMI among highly parous Amish women, both before and after menopause, and to evaluate the associations of parity with glucose, blood pressure, and lipids. METHODS: We conducted a cross-sectional study among 3,141 Amish women 18 years or older from Lancaster County, PA, who participated in our community-based Amish Research Program between 2003 and 2020. We evaluated the association between parity and BMI across different age groups, both before and after the menopausal transition. We further assessed associations between parity and cardiometabolic risk factors among the 1,128 postmenopausal women. Finally, we evaluated the association of change in parity with change in BMI in 561 women followed longitudinally. RESULTS: Approximately 62% of women in this sample (mean age, 45.2 y) reported having four or more children, and 36% reported having seven or more. A one-child increase in parity was associated with increased BMI in both premenopausal women (estimate [95% confidence interval], 0.4 kg/m 2 [0.2-0.5]) and to a lesser degree in postmenopausal women (0.2 kg/m 2 [0.02-0.3], Pint = 0.02), suggesting that the impact of parity on BMI decreases over time. Parity was not associated with glucose, blood pressure, total cholesterol, low-density lipoprotein, or triglycerides ( Padj > 0.05). CONCLUSIONS: Higher parity was associated with increased BMI in both premenopausal and postmenopausal women, but more so in younger/premenopausal women. Parity was not associated with other indices of cardiometabolic risk.