Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

An Epigenome-Wide Association Study of DNA Methylation and Proliferative Retinopathy over 28 Years in Type 1 Diabetes.

Purpose: To perform a prospective epigenome-wide association study of DNA methylation (DNAm) and 28-year proliferative diabetic retinopathy (PDR) incidence in type 1 diabetes (T1D). Design: Prospective observational cohort study. Participants: The Pittsburgh Epidemiology of Diabetes Complications (EDC) study of childhood-onset (< 17 years) T1D. Methods: Stereoscopic fundus photographs were taken in fields 1, 2, and 4 at baseline, 2, 4, 6, 8, 16, 23, and 28 years after DNAm measurements. The photos were graded using the modified Airlie House System. In those free of PDR at baseline (n = 265; mean T1D duration of 18 years at baseline), whole blood DNAm (EPIC array) at 683 597 CpGs was analyzed in Cox models for time to event. Associations between significant CpGs and clinical risk factors were assessed; genetic variants associated with DNAm were identified (methylation quantitative trait loci [meQTLs]). Mendelian randomization was used to examine evidence of causal associations between DNAm and PDR. Post hoc regional and functional analyses were performed. Main Outcome Measures: Proliferative diabetic retinopathy was defined as the first instance of a grade of ≥ 60 in at least 1 eye or pan-retinal photocoagulation for PDR. Follow-up time was calculated from the study visit at which DNAm data were available (baseline) until PDR incidence or censoring (December 31, 2018 or last follow-up). Results: PDR incidence was 53% over 28-years' follow-up. Greater DNAm of cg27512687 (KIF16B) was associated with reduced PDR incidence (P = 6.3 × 10-9; false discovery rate [FDR]: < 0.01); 113 cis-meQTLs (P < 5 × 10-8) were identified. Mendelian randomization analysis using the sentinel meQTL as the instrumental variable supported a potentially causal association between cg27512687 and PDR. Cg27512687 was also associated with lower pulse rate and albumin excretion rate and higher estimated glomerular filtration rate, but its association with PDR remained independently significant after adjustment for those factors. In regional analyses, DNAm of FUT4, FKBP1A, and RIN2 was also associated with PDR incidence. Conclusions: DNA methylation of KIF16B, FUT4, FKBP1A, and RIN2 was associated with PDR incidence, supporting roles for epigenetic regulation of iron clearance, developmental pathways, and autophagy in PDR pathogenesis. Further study of those loci may provide insight into novel targets for interventions to prevent or delay PDR in T1D. Financial Disclosures: Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Mosaic chromosomal alterations in blood across ancestries using whole-genome sequencing.

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis.

Myocardial Fibrosis and Cardiomyopathy Risk: A Genetic Link in the MESA.

BACKGROUND: Common genetic variants are associated with risk for hypertrophic cardiomyopathy and dilated cardiomyopathy and with left ventricular (LV) traits. Whether these variants are associated with myocardial fibrosis, an important pathophysiological mediator of cardiomyopathy, is unknown. METHODS: Multi-Ethnic Study of Atherosclerosis participants with T1-mapping cardiac magnetic resonance imaging in-whom extracellular volume was assessed, and genotyping information was available were included (N=1255). Log extracellular volume (%) was regressed on 50 candidate single nucleotide polymorphisms (previously identified to be associated with hypertrophic cardiomyopathy, dilated cardiomyopathy, and LV traits) adjusting for age, sex, diabetes, blood pressure, and principal components of ancestry. Ancestry-specific results were pooled by fixed-effect meta-analyses. Gene knockdown experiments were performed in human cardiac fibroblasts. RESULTS: The SMARCB1 rs2186370 intronic variant (minor allele frequency: 0.18 in White and 0.50 in Black participants), previously identified as a risk variant for dilated cardiomyopathy and hypertrophic cardiomyopathy, was significantly associated with increased extracellular volume (P=0.0002) after adjusting for confounding clinical variables. The SMARCB1 rs2070458 locus previously associated with increased LV wall thickness and mass was similarly significantly associated with increased extracellular volume (P=0.0002). The direction of effect was similar in all 4 ancestry groups, but the effect was strongest in Black participants. The variants are strong expression quantitative loci in human LV tissue and associated with genotype-dependent decreased expression of SMARCB1 (P=7.3×10-22). SMARCB1 knockdown in human cardiac fibroblasts resulted in increased TGF (transforming growth factor)-ß1-mediated α-smooth muscle actin and collagen expression. CONCLUSIONS: Common genetic variation in SMARCB1 previously associated with risk for cardiomyopathies and increased LV wall thickness is associated with increased cardiac magnetic resonance imaging-based myocardial fibrosis and increased TGF-ß1 mediated myocardial fibrosis in vitro. Whether these findings suggest a pathophysiologic link between myocardial fibrosis and cardiomyopathy risk remains to be proven.

DNA methylation and 28-year cardiovascular disease risk in type 1 diabetes: the Epidemiology of Diabetes Complications (EDC) cohort study.

BACKGROUND: The potential for DNA methylation (DNAm) as an early marker for cardiovascular disease (CVD) and how such an association might differ by glycemic exposure has not been examined in type 1 diabetes, a population at increased CVD risk. We thus performed a prospective epigenome-wide association study of blood leukocyte DNAm (EPIC array) and time to CVD incidence over 28 years in a childhood-onset (< 17 years) type 1 diabetes cohort, the Pittsburgh Epidemiology of Diabetes Complications (EDC) study (n = 368 with DNA and no CVD at baseline), both overall and separately by glycemic exposure, as measured by HbA1c at baseline (split at the median: < 8.9% and ≥ 8.9%). We also assessed whether DNAm-CVD associations were independent of established cardiometabolic risk factors, including body mass index, estimated glucose disposal rate, cholesterol, triglycerides, blood pressure, pulse rate, albumin excretion rate, and estimated glomerular filtration rate. RESULTS: CVD (first instance of CVD death, myocardial infarction, coronary revascularization, ischemic ECG, angina, or stroke) developed in 172 participants (46.7%) over 28 years. Overall, in Cox regression models for time to CVD, none of the 683,597 CpGs examined reached significance at a false discovery rate (FDR) ≤ 0.05. In participants with HbA1c < 8.9% (n = 180), again none reached FDR ≤ 0.05, but three were associated at the a priori nominal significance level FDR ≤ 0.10: cg07147033 in MIB2, cg12324048 (intergenic, chromosome 3), and cg15883830 (intergenic, chromosome 1). In participants with HbA1c ≥ 8.9% (n = 188), two CpGs in loci involved in calcium channel activity were significantly associated with CVD (FDR ≤ 0.05): cg21823999 in GPM6A and cg23621817 in CHRNA9; four additional CpGs were nominally associated (FDR ≤ 0.10). In participants with HbA1c ≥ 8.9%, DNAm-CVD associations were only modestly attenuated after cardiometabolic risk factor adjustment, while attenuation was greater in those with HbA1c < 8.9%. No pathways were enriched in those with HbA1c < 8.9%, while pathways for calcium channel activity and integral component of synaptic membrane were significantly enriched in those with HbA1c ≥ 8.9%. CONCLUSIONS: These results provide novel evidence that DNAm at loci involved in calcium channel activity and development may contribute to long-term CVD risk beyond known risk factors in type 1 diabetes, particularly in individuals with greater glycemic exposure, warranting further study.

Leveraging type 1 diabetes human genetic and genomic data in the T1D knowledge portal.

To address the challenge of translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight, we have developed the T1D Knowledge Portal (T1DKP), an open-access resource for hypothesis development and target discovery in T1D.

Associations between Ambient Air Pollutants and Clonal Hematopoiesis of Indeterminate Potential.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP) is an age-related somatic mutation associated with incident hematologic cancer. Environmental stressors which, like air pollution, generate oxidative stress at the cellular level, may induce somatic mutations and some mutations may provide a selection advantage for persistence and expansion of specific clones. METHODS: We used data from the Multi-Ethnic Study of Atherosclerosis (MESA) N = 4,379 and the Women's Health Initiative (WHI) N = 7,701 to estimate cross-sectional associations between annual average air pollution concentrations at participant address the year before blood draw using validated spatiotemporal models. We used covariate-adjusted logistic regression to estimate risk of CHIP per interquartile range increases in particulate matter (PM2.5; 4 µg/m3) and nitrogen dioxide (NO2; 10 ppb) as ORs (95% confidence intervals). RESULTS: Prevalence of CHIP at blood draw (variant allele fraction > 2%) was 4.4% and 8.7% in MESA and WHI, respectively. The most common CHIP driver mutation was in DNMT3A. Neither pollutant was associated with CHIP: ORMESA PM2.5 = 1.00 (0.68-1.45), ORMESA NO2 = 1.05 (0.69-1.61), ORWHI PM2.5 = 0.97 (0.86-1.09), ORWHI NO2 = 0.98 (0.88-1.10); or with DNMT3A-driven CHIP. CONCLUSIONS: We did not find evidence that air pollution contributes to CHIP prevalence in two large observational cohorts. IMPACT: This is the first study to estimate associations between air pollution and CHIP.

Leveraging type 1 diabetes human genetic and genomic data in the T1D Knowledge Portal.

Translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight can reveal novel biology and therapeutic targets but remains a major challenge. We developed the T1D Knowledge Portal (T1DKP), a disease-specific resource of genetic and functional annotation data that enables users to develop hypotheses for T1D-based research and target discovery. The T1DKP can be used to query genes and genomic regions for genetic associations, identify epigenomic features, access results of bioinformatic analyses, and obtain expert-curated resources. The T1DKP is available at http://t1d.hugeamp.org .

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

TXNIP DNA methylation is associated with glycemic control over 28 years in type 1 diabetes: findings from the Pittsburgh Epidemiology of Diabetes Complications (EDC) study.

INTRODUCTION: DNA methylation (DNAme) has been cross-sectionally associated with type 2 diabetes and hemoglobin A1c (HbA1c) in the general population. However, longitudinal data and data in type 1 diabetes are currently very limited. Thus, we performed an epigenome-wide association study (EWAS) in an observational type 1 diabetes cohort to identify loci with DNAme associated with concurrent and future HbA1cs, as well as other clinical risk factors, over 28 years. RESEARCH DESIGN AND METHODS: Whole blood DNAme in 683 597 CpGs was analyzed in the Pittsburgh Epidemiology of Diabetes Complications study of childhood onset (<17 years) type 1 diabetes (n=411). An EWAS of DNAme beta values and concurrent HbA1c was performed using linear models adjusted for diabetes duration, sex, pack years of smoking, estimated cell type composition variables, and technical/batch covariates. A longitudinal EWAS of subsequent repeated HbA1c measures was performed using mixed models. We further identified methylation quantitative trait loci (meQTLs) for significant CpGs and conducted a Mendelian randomization. RESULTS: DNAme at cg19693031 (Chr 1, Thioredoxin-Interacting Protein (TXNIP)) and cg21534330 (Chr 17, Casein Kinase 1 Isoform Delta) was significantly inversely associated with concurrent HbA1c. In longitudinal analyses, hypomethylation of cg19693031 was associated with consistently higher HbA1c over 28 years, and with higher triglycerides, pulse rate, and albumin:creatinine ratio (ACR) independently of HbA1c. We further identified 34 meQTLs in SLC2A1/SLC2A1-AS1 significantly associated with cg19693031 DNAme. CONCLUSIONS: Our results extend prior findings that TXNIP hypomethylation relates to worse glycemic control in type 1 diabetes by demonstrating the association persists over the long term. Additionally, the associations with triglycerides, pulse rate, and ACR suggest TXNIP DNAme could play a role in vascular damage independent of HbA1c. These findings strengthen potential for interventions targeting TXNIP to improve glycemic control in type 1 diabetes through its role in SLC2A1/glucose transporter 1-mediated glucose regulation.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Circulating Soluble CD163, Associations With Cardiovascular Outcomes and Mortality, and Identification of Genetic Variants in Older Individuals: The Cardiovascular Health Study.

Background Monocytes/macrophages participate in cardiovascular disease. CD163 (cluster of differentiation 163) is a monocyte/macrophage receptor, and the shed sCD163 (soluble CD163) reflects monocyte/macrophage activation. We examined the association of sCD163 with incident cardiovascular disease events and performed a genome-wide association study to identify sCD163-associated variants. Methods and Results We measured plasma sCD163 in 5214 adults (aged ≥65 years, 58.7% women, 16.2% Black) of the CHS (Cardiovascular Health Study). We used Cox regression models (associations of sCD163 with incident events and mortality); median follow-up was 26 years. Genome-wide association study analyses were stratified on race. Adjusted for age, sex, and race and ethnicity, sCD163 levels were associated with all-cause mortality (hazard ratio [HR], 1.08 [95% CI, 1.04-1.12] per SD increase), cardiovascular disease mortality (HR, 1.15 [95% CI, 1.09-1.21]), incident coronary heart disease (HR, 1.10 [95% CI, 1.04-1.16]), and incident heart failure (HR, 1.18 [95% CI, 1.12-1.25]). When further adjusted (eg, cardiovascular disease risk factors), only incident coronary heart disease lost significance. In European American individuals, genome-wide association studies identified 38 variants on chromosome 2 near MGAT5 (top result rs62165726, P=3.3×10-18),19 variants near chromosome 17 gene ASGR1 (rs55714927, P=1.5×10-14), and 18 variants near chromosome 11 gene ST3GAL4. These regions replicated in the European ancestry ADDITION-PRO cohort, a longitudinal cohort study nested in the Danish arm of the Anglo-Danish-Dutch study of Intensive Treatment Intensive Treatment In peOple with screeNdetcted Diabetes in Primary Care. In Black individuals, we identified 9 variants on chromosome 6 (rs3129781 P=7.1×10-9) in the HLA region, and 3 variants (rs115391969 P=4.3×10-8) near the chromosome 16 gene MYLK3. Conclusions Monocyte function, as measured by sCD163, may be predictive of overall and cardiovascular-specific mortality and incident heart failure.

Hematopoietic loss of Y chromosome leads to cardiac fibrosis and heart failure mortality.

Hematopoietic mosaic loss of Y chromosome (mLOY) is associated with increased risk of mortality and age-related diseases in men, but the causal and mechanistic relationships have yet to be established. Here, we show that male mice reconstituted with bone marrow cells lacking the Y chromosome display increased mortality and age-related profibrotic pathologies including reduced cardiac function. Cardiac macrophages lacking the Y chromosome exhibited polarization toward a more fibrotic phenotype, and treatment with a transforming growth factor ß1-neutralizing antibody ameliorated cardiac dysfunction in mLOY mice. A prospective study revealed that mLOY in blood is associated with an increased risk for cardiovascular disease and heart failure-associated mortality. Together, these results indicate that hematopoietic mLOY causally contributes to fibrosis, cardiac dysfunction, and mortality in men.

Oral rotavirus vaccine shedding as a marker of mucosal immunity.

Group A rotaviruses (RVA) remain a leading cause of pediatric diarrhea worldwide, in part due to underperformance of currently approved live-attenuated, oral vaccines in low-and-middle income countries. Improved immune correlates of protection (CoP) for existing oral vaccines and novel strategies to evaluate the performance of next-generation vaccines are needed. Use of oral vaccines as challenge agents in controlled human infection models is a potential approach to CoP discovery that remains underexplored. In a live-attenuated, oral rotavirus vaccine (Rotarix, GlaxoSmithKline) efficacy trial conducted among infants in Dhaka, Bangladesh, we explored the potential for the second dose of the two-dose series to be considered a challenge agent through which RVA immunity could be explored, using fecal virus shedding post-dose 2 as a marker of mucosal immunity. Among 180 vaccinated infants who completed the parent study per protocol, the absence of fecal vaccine shedding following the second dose of Rotarix suggested intestinal mucosal immunity generated by the first dose and a decreased risk of RVA diarrhea through 2 years of life (RR 0.616, 95% CI 0.392-0.968). Further development of controlled human infection models for group A rotaviruses, especially in prospective studies with larger sample sizes, may be a promising tool to assess rotavirus vaccine efficacy and CoPs.

Population sequencing data reveal a compendium of mutational processes in the human germ line.

Biological mechanisms underlying human germline mutations remain largely unknown. We statistically decompose variation in the rate and spectra of mutations along the genome using volume-regularized nonnegative matrix factorization. The analysis of a sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. We provide a biological interpretation for seven of these processes. We associate one process with bulky DNA lesions that are resolved asymmetrically with respect to transcription and replication. Two processes track direction of replication fork and replication timing, respectively. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions and a mutagenic effect of long interspersed nuclear elements. We localize a mutagenic process specific to oocytes from population sequencing data. This process appears transcriptionally asymmetric.

Association of Coding Variants in Hydroxysteroid 17-beta Dehydrogenase 14 (HSD17B14) with Reduced Progression to End Stage Kidney Disease in Type 1 Diabetes.

BACKGROUND: Rare variants in gene coding regions likely have a greater impact on disease-related phenotypes than common variants through disruption of their encoded protein. We searched for rare variants associated with onset of ESKD in individuals with type 1 diabetes at advanced kidney disease stage. METHODS: Gene-based exome array analyses of 15,449 genes in five large incidence cohorts of individuals with type 1 diabetes and proteinuria were analyzed for survival time to ESKD, testing the top gene in a sixth cohort (n=2372/1115 events all cohorts) and replicating in two retrospective case-control studies (n=1072 cases, 752 controls). Deep resequencing of the top associated gene in five cohorts confirmed the findings. We performed immunohistochemistry and gene expression experiments in human control and diseased cells, and in mouse ischemia reperfusion and aristolochic acid nephropathy models. RESULTS: Protein coding variants in the hydroxysteroid 17-ß dehydrogenase 14 gene (HSD17B14), predicted to affect protein structure, had a net protective effect against development of ESKD at exome-wide significance (n=4196; P value=3.3 × 10-7). The HSD17B14 gene and encoded enzyme were robustly expressed in healthy human kidney, maximally in proximal tubular cells. Paradoxically, gene and protein expression were attenuated in human diabetic proximal tubules and in mouse kidney injury models. Expressed HSD17B14 gene and protein levels remained low without recovery after 21 days in a murine ischemic reperfusion injury model. Decreased gene expression was found in other CKD-associated renal pathologies. CONCLUSIONS: HSD17B14 gene is mechanistically involved in diabetic kidney disease. The encoded sex steroid enzyme is a druggable target, potentially opening a new avenue for therapeutic development.

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci.

BACKGROUND: DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach. METHODS: The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses. RESULTS: We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development. CONCLUSIONS: We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Genetic landscape of Gullah African Americans.

OBJECTIVES: Gullah African Americans are descendants of formerly enslaved Africans living in the Sea Islands along the coast of the southeastern U.S., from North Carolina to Florida. Their relatively high numbers and geographic isolation were conducive to the development and preservation of a unique culture that retains deep African features. Although historical evidence supports a West-Central African ancestry for the Gullah, linguistic and cultural evidence of a connection to Sierra Leone has led to the suggestion of this country/region as their ancestral home. This study sought to elucidate the genetic structure and ancestry of the Gullah. MATERIALS AND METHODS: We leveraged whole-genome genotype data from Gullah, African Americans from Jackson, Mississippi, African populations from Sierra Leone, and population reference panels from Africa and Europe to infer population structure, ancestry proportions, and global estimates of admixture. RESULTS: Relative to non-Gullah African Americans from the Southeast US, the Gullah exhibited higher mean African ancestry, lower European admixture, a similarly small Native American contribution, and increased male-biased European admixture. A slightly tighter bottleneck in the Gullah 13 generations ago suggests a largely shared demographic history with non-Gullah African Americans. Despite a slightly higher relatedness to populations from Sierra Leone, our data demonstrate that the Gullah are genetically related to many West African populations. DISCUSSION: This study confirms that subtle differences in African American population structure exist at finer regional levels. Such observations can help to inform medical genetics research in African Americans, and guide the interpretation of genetic data used by African Americans seeking to explore ancestral identities.

Whole genome sequence analyses of eGFR in 23,732 people representing multiple ancestries in the NHLBI trans-omics for precision medicine (TOPMed) consortium.

BACKGROUND: Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants. METHODS: We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity. FINDINGS: When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10-11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10-9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10-9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10-9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10-9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants. INTERPRETATION: This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

Maternal Secretor Status Affects Oral Rotavirus Vaccine Response in Breastfed Infants in Bangladesh.

Secretor status controls mucosal histo-blood group antigen expression and is associated with susceptibility to rotavirus (RV) diarrhea, with nonsecretors less susceptible to symptomatic infection. The role of breast milk secretor status on oral live-attenuated RV vaccine response in breastfed infants has not been explored. In a monovalent G1P[8] RV vaccine (Rotarix) trial in Bangladesh, RV-specific plasma immunoglobulin A antibody seroconversion rates were higher among infants of maternal nonsecretors (39%) than infants of maternal secretors (23%; P = .001). Maternal status remained a significant predictor when correcting for infant status (P = .002). Maternal secretor status should be considered when interpreting oral RV vaccine responses in low- and middle-income settings. Clinical Trials Registration. NCT01375647.