Searching across-cohort relatives in 54,092 GWAS samples via encrypted genotype regression.

Explicitly sharing individual level data in genomics studies has many merits comparing to sharing summary statistics, including more strict QCs, common statistical analyses, relative identification and improved statistical power in GWAS, but it is hampered by privacy or ethical constraints. In this study, we developed encG-reg, a regression approach that can detect relatives of various degrees based on encrypted genomic data, which is immune of ethical constraints. The encryption properties of encG-reg are based on the random matrix theory by masking the original genotypic matrix without sacrificing precision of individual-level genotype data. We established a connection between the dimension of a random matrix, which masked genotype matrices, and the required precision of a study for encrypted genotype data. encG-reg has false positive and false negative rates equivalent to sharing original individual level data, and is computationally efficient when searching relatives. We split the UK Biobank into their respective centers, and then encrypted the genotype data. We observed that the relatives estimated using encG-reg was equivalently accurate with the estimation by KING, which is a widely used software but requires original genotype data. In a more complex application, we launched a finely devised multi-center collaboration across 5 research institutes in China, covering 9 cohorts of 54,092 GWAS samples. encG-reg again identified true relatives existing across the cohorts with even different ethnic backgrounds and genotypic qualities. Our study clearly demonstrates that encrypted genomic data can be used for data sharing without loss of information or data sharing barrier.

Causal impact of DNA methylation on refracture in elderly individuals with osteoporosis - a prospective cohort study.

BACKGROUND: Osteoporotic vertebral compression fractures (OVCF) in the elderly increase refracture risk post-surgery, leading to higher mortality rates. Genome-wide association studies (GWAS) have identified susceptibility genes for osteoporosis, but the phenotypic variance explained by these genes has been limited, indicating the need to explore additional causal factors. Epigenetic modifications, such as DNA methylation, may influence osteoporosis and refracture risk. However, prospective cohorts for assessing epigenetic alterations in Chinese elderly patients are lacking. Here, we propose to conduct a prospective cohort study to investigate the causal network of DNA polymorphisms, DNA methylation, and environmental factors on the development of osteoporosis and the risk of refracture. METHODS: We will collect vertebral and peripheral blood from 500 elderly OVCF patients undergoing surgery, extract DNA, and generate whole genome genotype data and DNA methylation data. Observation indicators will be collected and combined with one-year follow-up data. A healthy control group will be selected from a natural population cohort. Epigenome-wide association studies (EWAS) of osteoporosis and bone mineral density will be conducted. Differential methylation analysis will compare candidate gene methylation patterns in patients with and without refracture. Multi-omics prediction models using genetic variants and DNA methylation sites will be built to predict OVCF risk. DISCUSSION: This study will be the first large-scale population-based study of osteoporosis and bone mineral density phenotypes based on genome-wide data, multi-time point methylation data, and phenotype data. By analyzing methylation changes related to osteoporosis and bone mineral density in OVCF patients, the study will explore the feasibility of DNA methylation in evaluating postoperative osteoporosis intervention effects. The findings may identify new molecular markers for effective anti-osteoporosis treatment and inform individualized prevention and treatment strategies. TRIAL REGISTRATION: chictr.org.cn ChiCTR2200065316, 02/11/2022.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development.

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Investigating gene functions and single-cell expression profiles of de novo variants in orofacial clefts.

Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs' deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p < 0.0002), while only DNV genes were enriched in the endothelium (adjusted p = 0.032). Although we did not achieve significant findings using CPO gene sets, which was mainly due to the limited number of DNV genes, scRNA-seq analysis implicated various expression patterns among DNV and OMIM genes. Our results suggest that combinatory pathway and scRNA-seq data analyses are helpful for contextualizing genes in OFC development.

Analysis of blood methylation quantitative trait loci in East Asians reveals ancestry-specific impacts on complex traits.

Methylation quantitative trait loci (mQTLs) are essential for understanding the role of DNA methylation changes in genetic predisposition, yet they have not been fully characterized in East Asians (EAs). Here we identified mQTLs in whole blood from 3,523 Chinese individuals and replicated them in additional 1,858 Chinese individuals from two cohorts. Over 9% of mQTLs displayed specificity to EAs, facilitating the fine-mapping of EA-specific genetic associations, as shown for variants associated with height. Trans-mQTL hotspots revealed biological pathways contributing to EA-specific genetic associations, including an ERG-mediated 233 trans-mCpG network, implicated in hematopoietic cell differentiation, which likely reflects binding efficiency modulation of the ERG protein complex. More than 90% of mQTLs were shared between different blood cell lineages, with a smaller fraction of lineage-specific mQTLs displaying preferential hypomethylation in the respective lineages. Our study provides new insights into the mQTL landscape across genetic ancestries and their downstream effects on cellular processes and diseases/traits.

DNA methylation clocks for estimating biological age in Chinese cohorts.

Epigenetic clocks are accurate predictors of human chronological age based on the analysis of DNA methylation (DNAm) at specific CpG sites. However, a systematic comparison between DNA methylation data and other omics datasets has not yet been performed. Moreover, available DNAm age predictors are based on datasets with limited ethnic representation. To address these knowledge gaps, we generated and analyzed DNA methylation datasets from two independent Chinese cohorts, revealing age-related DNAm changes. Additionally, a DNA methylation aging clock (iCAS-DNAmAge) and a group of DNAm-based multi-modal clocks for Chinese individuals were developed, with most of them demonstrating strong predictive capabilities for chronological age. The clocks were further employed to predict factors influencing aging rates. The DNAm aging clock, derived from multi-modal aging features (compositeAge-DNAmAge), exhibited a close association with multi-omics changes, lifestyles, and disease status, underscoring its robust potential for precise biological age assessment. Our findings offer novel insights into the regulatory mechanism of age-related DNAm changes and extend the application of the DNAm clock for measuring biological age and aging pace, providing the basis for evaluating aging intervention strategies.