ABSTRACT
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.
Subject(s)
Genome-Wide Association Study , Privacy , Humans , Genome-Wide Association Study/methods , Genotype , Software , GenomicsABSTRACT
To examine whether circulating interleukin-6 (IL-6) levels (CirIL6) have a causal effect on blood pressure using Mendelian randomization (MR) methods. We used data from genome-wide association studies (GWAS) of European ancestry to obtain genetic instruments for circulating IL-6 levels and blood pressure measurements. We applied several robust MR methods to estimate the causal effects and to test for heterogeneity and pleiotropy. We found that circulating IL-6 had a significant positive causal effect on systolic blood pressure (SBP) and pulmonary arterial hypertension (PAH), but not on diastolic blood pressure (DBP) or hypertension. We found that as CirIL6 genetically increased, SBP increased using Inverse Variance Weighted (IVW) method (for ukb-b-20175, ß = 0.082 with SE = 0.032, P = 0.011; for ukb-a-360, ß = 0.075 with SE = 0.031, P = 0.014) and weighted median (WM) method (for ukb-b-20175, ß = 0.061 with SE = 0.022, P = 0.006; for ukb-a-360, ß = 0.065 with SE = 0.027, P = 0.014). Moreover, CirIL6 may be associated with an increased risk of PAH using WM method (odds ratio (OR) = 15.503, 95% CI, 1.025-234.525, P = 0.048), but not with IVW method. Our study provides novel evidence that circulating IL-6 has a causal role in the development of SBP and PAH, but not DBP or hypertension. These findings suggest that IL-6 may be a potential therapeutic target for preventing or treating cardiovascular diseases and metabolic disorders. However, more studies are needed to confirm the causal effects of IL-6 on blood pressure and to elucidate the underlying mechanisms and pathways.
Subject(s)
Hypertension , Interleukin-6 , Humans , Blood Pressure/genetics , Interleukin-6/genetics , Genome-Wide Association Study , Mendelian Randomization Analysis , Hypertension/geneticsABSTRACT
BACKGROUND: Obesity and hypertension are major risk factors for cardiovascular diseases that affect millions of people worldwide. Both conditions are associated with chronic low-grade inflammation, which is mediated by adipokines such as adiponectin. Adiponectin is the most abundant adipokine that has a beneficial impact on metabolic and vascular biology, while high serum concentrations are associated with some syndromes. This "adiponectin paradox" still needs to be clarified in obesity-associated hypertension. The aim of this study was to investigate how adiponectin affects blood pressure, inflammation, and metabolic function in obesity hypertension using a Chinese adult case-control study. METHODS: A case-control study that had finished recruiting 153 subjects divided as four characteristic groups. Adiponectin serum levels were tested by ELISA in these subjects among these four characteristic Chinese adult physical examination groups. Waist circumference (WC), body mass index (BMI), systolic blood pressure (SB), diastolic blood pressure (DB), and other clinical laboratory data were collected. Analyzation of correlations between the research index and differences between groups was done by SPSS. RESULTS: Serum adiponectin levels in the| normal healthy group (NH group) were significantly higher than those in the newly diagnosed untreated just-obesity group (JO group), and negatively correlated with the visceral adiposity index. With multiple linear egression analysis, it was found that, for serum adiponectin, gender, serum albumin (ALB), alanine aminotransferase (ALT) and high-density lipoprotein cholesterol (HDLC) were the significant independent correlates, and for SB, age and HDLC were the significant independent correlates, and for DB, alkaline phosphatase (ALP) was the significant independent correlate. The other variables did not reach significance in the model. CONCLUSIONS: Our study reveals that adiponectin's role in obesity-hypertension is multifaceted and is influenced by the systemic metabolic homeostasis signaling axis. In obesity-related hypertension, compensatory effects, adiponectin resistance, and reduced adiponectin clearance from impaired kidneys and liver all contribute to the "adiponectin paradox".
Subject(s)
Adiponectin , Hypertension , Adult , Humans , Case-Control Studies , Hypertension/diagnosis , Obesity/complications , Obesity/diagnosis , Cholesterol, HDL , Inflammation , China/epidemiologyABSTRACT
Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome-wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined superenhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range 'mSE-miRNA' approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).
Subject(s)
Bone Density/genetics , Enhancer Elements, Genetic , Epigenomics , Gene Expression Regulation , Gene Regulatory Networks , Genomics , MicroRNAs/genetics , Computational Biology , Epigenomics/methods , Gene Expression Profiling , Genome-Wide Association Study , Genomics/methods , Humans , Osteoblasts/metabolism , Polymorphism, Single Nucleotide , Protein Interaction MapsABSTRACT
BACKGROUND: Obesity and hypertension are major risk factors for cardiovascular diseases that affect millions of people worldwide. Both conditions are associated with chronic low-grade inflammation, which is mediated by cytokines such as interleukin-6 (IL-6). IL-6 is a multifunctional cytokine that can have pro-inflammatory or anti-inflammatory effects depending on the context. The exact role of IL-6 in obesity-associated hypertension is unclear. OBJECTIVE: To investigate how IL-6 affects blood pressure, inflammation, and metabolic function in obesity-hypertension using a Chinese adult case-control study. METHODS: A total of 153 participants were sorted into four subgroups according to their body mass index (BMI) and blood pressure (BP): normal healthy group (NH), just obesity group (JO), just-hypertension group (JH), and obesity-hypertension group (OH). Serum IL-6 concentrations were measured by Enzyme-linked Immunosorbent Assay (ELISA) and their correlations with anthropometric and laboratory parameters and their differences across the subgroups were examined. Multiple linear regression analysis was performed to identify the predictors of serum IL-6 concentrations in each group. RESULTS: Serum IL-6 concentrations were higher in NH group than in JO group and correlated positively with diastolic blood pressure in NH and JO groups, but not in JH and OH groups. Serum IL-6 concentrations also correlated with albumin in NH group, alkaline phosphatase in JO group, serum creatinine and fasting blood glucose in JH group. The influencing factors of serum IL-6 concentrations varied among the four groups, with gender, diastolic blood pressure and albumin being significant predictors in NH group, alkaline phosphatase in JO group, age and serum creatinine in JH group, and none in OH group. CONCLUSIONS: These results suggest that IL-6 may play diverse effects in the pathogenesis of obesity- hypertension, depending on the presence or absence of obesity and hypertension. Further studies are needed to elucidate the underlying mechanisms of IL-6 signaling and function in these diseases.
Subject(s)
Hypertension , Interleukin-6 , Humans , Adult , Case-Control Studies , Alkaline Phosphatase , Creatinine , East Asian People , Obesity , Cytokines , Body Mass Index , Inflammation , AlbuminsABSTRACT
Coexistence and cooperation between dogs and humans over thousands of years have supported convergent evolutionary processes in the two species. Previous studies found that Eurasian dogs evolved into a distinct geographic cluster. In this study, we used the genomes of 242 European dogs, 38 Southeast Asian indigenous (SEAI) dogs, and 41 gray wolves to identify adaptation of European dogs . We report 86 unique positively selected genes in European dogs, among which is LCT (lactase). LCT encodes lactase, which is fundamental for the digestion of lactose. We found that an A-to-G mutation (chr19:38,609,592) is almost fixed in Middle Eastern and European dogs. The results of two-dimensional site frequency spectrum (2D SFS) support that the mutation is under soft sweep . We inferred that the onset of positive selection of the mutation is shorter than 6,535 years and behind the well-developed dairy economy in central Europe. It increases the expression of LCT by reducing its binding with ZEB1, which would enhance dog's ability to digest milk-based diets. Our study uncovers the genetic basis of convergent evolution between humans and dogs with respect to diet, emphasizing the import of the dog as a biomedical model for studying mechanisms of the digestive system.
Subject(s)
Lactase , Selection, Genetic , Animals , Dogs , Gene Frequency , Humans , Lactase/genetics , Lactase/metabolism , Lactose/metabolism , Polymorphism, Single Nucleotide , White PeopleABSTRACT
BACKGROUND: Birth weight is considered not only to undermine future growth, but also to induce lifelong diseases; the aim of this study is to explore the relationship between birth weight and adult bone mass. METHODS: We performed multivariable regression analyses to assess the association of birth weight with bone parameters measured by dual-energy X-ray absorptiometry (DXA) and by quantitative ultrasound (QUS), independently. We also implemented a systemic Mendelian randomization (MR) analysis to explore the causal association between them with both fetal-specific and maternal-specific instrumental variables. RESULTS: In the observational analyses, we found that higher birth weight could increase the adult bone area (lumbar spine, ß-coefficient= 0.17, P < 2.00 × 10-16; lateral spine, ß-coefficient = 0.02, P = 0.04), decrease bone mineral content-adjusted bone area (BMCadjArea) (lumbar spine, ß-coefficient= - 0.01, P = 2.27 × 10-14; lateral spine, ß-coefficient = - 0.05, P = 0.001), and decrease adult bone mineral density (BMD) (lumbar spine, ß-coefficient = - 0.04, P = 0.007; lateral spine; ß-coefficient = - 0.03, P = 0.02; heel, ß-coefficient = - 0.06, P < 2.00 × 10-16), and we observed that the effect of birth weight on bone size was larger than that on BMC. In MR analyses, the higher fetal-specific genetically determined birth weight was identified to be associated with higher bone area (lumbar spine; ß-coefficient = 0.15, P = 1.26 × 10-6, total hip, ß-coefficient = 0.15, P = 0.005; intertrochanteric area, ß-coefficient = 0.13, P = 0.0009; trochanter area, ß-coefficient = 0.11, P = 0.03) but lower BMD (lumbar spine, ß-coefficient = - 0.10, P = 0.01; lateral spine, ß-coefficient = - 0.12, P = 0.0003, and heel ß-coefficient = - 0.11, P = 3.33 × 10-13). In addition, we found that the higher maternal-specific genetically determined offspring birth weight was associated with lower offspring adult heel BMD (ß-coefficient = - 0.001, P = 0.04). CONCLUSIONS: The observational analyses suggested that higher birth weight was associated with the increased adult bone area but decreased BMD. By leveraging the genetic instrumental variables with maternal- and fetal-specific effects on birth weight, the observed relationship could be reflected by both the direct fetal and indirect maternal genetic effects.
Subject(s)
Bone Density , Lumbar Vertebrae , Absorptiometry, Photon , Adult , Birth Weight , Bone Density/genetics , Humans , Lumbar Vertebrae/diagnostic imaging , Mendelian Randomization AnalysisABSTRACT
BACKGROUND: Polydactyly is one of the most common congenital limb dysplasia in many animal species. Although preaxial polydactyly (PPD) has been comprehensively studied in humans as a common abnormality, the genetic variations in other animal species have not been fully understood. Herein, we focused on the pig, as an even-toed ungulate mammal model with its unique advantages in medical and genetic researches, two PPD families consisting of four affected and 20 normal individuals were sequenced. RESULTS: Our results showed that the PPD in the sampled pigs were not related to previously reported variants. A strong association was identified at ABCC4 and it encodes a transmembrane protein involved in ciliogenesis. We found that the affected and normal individuals were highly differentiated at ABCC4, and all the PPD individuals shared long haplotype stretches as compared with the unaffected individuals. A highly differentiated missense mutation (I85T) in ABCC4 was observed at a residue from a transmembrane domain highly conserved among a variety of organisms. CONCLUSIONS: This study reports ABCC4 as a new candidate gene and identifies a missense mutation for PPD in pigs. Our results illustrate a putative role of ciliogenesis process in PPD, coinciding with an earlier observation of ciliogenesis abnormality resulting in pseudo-thumb development in pandas. These results expand our knowledge on the genetic variations underlying PPD in animals.
Subject(s)
Multidrug Resistance-Associated Proteins/genetics , Polydactyly/genetics , Whole Genome Sequencing/methods , Animals , Female , Male , Mutation, Missense/genetics , Pedigree , SwineABSTRACT
The evolutionary and functional features of RNA editing are well studied in mammals, cephalopods, and insects, but not in birds. Here, we integrated transcriptomic and whole-genomic analyses to exhaustively characterize the expansive repertoire of adenosine-to-inosine (A-to-I) RNA editing sites (RESs) in the chicken. In addition, we investigated the evolutionary status of the chicken editome as a potential mechanism of domestication. We detected the lowest editing level in the liver of chickens, compared to muscles in humans, and found higher editing activity and specificity in the brain than in non-neural tissues, consistent with the brain's functional complexity. To a certain extent, specific editing activity may account for the specific functions of tissues. Our results also revealed that sequences critical to RES secondary structures remained conserved within avian evolution. Furthermore, the RNA editome was shaped by purifying selection during chicken domestication and most RESs may have served as a selection pool for a few functional RESs involved in chicken domestication, including evolution of nervous and immune systems. Regulation of RNA editing in chickens by adenosine deaminase acting on RNA (ADAR) enzymes may be affected by non-ADAR factors whose expression levels changed widely after ADAR knockdown. Collectively, we provide comprehensive lists of candidate RESs and non-ADAR-editing regulators in the chicken, thus contributing to our current understanding of the functions and evolution of RNA editing in animals.
Subject(s)
Adenosine Deaminase , Chickens , RNA Editing , Animals , Humans , Adenosine/genetics , Adenosine Deaminase/genetics , Adenosine Deaminase/metabolism , Chickens/genetics , Genomics , Inosine/genetics , RNA/genetics , TranscriptomeABSTRACT
We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.
Subject(s)
Asian People , Biological Specimen Banks , Asian People/genetics , China , Genomics , Humans , Pilot ProjectsABSTRACT
We combined conventional evidence from longitudinal data in UK Biobank and genetic evidence from Mendelian randomization (MR) approach to infer the causality between sleep behaviors and fracture risk. We found that participants with insomnia showed 6.4% higher risk of fracture (hazard ratio [HR] = 1.064, 95% CI = 1.038-1.090, P = 7.84 × 10-7), falls and bone mineral density (BMD) mediated 24.6% and 10.6% of the intermediary effect; the MR analyses provided the consistent evidence. A U-shape relationship was observed between sleep duration and fracture risk (P < 0.001) with the lowest risk at sleeping 7-8 h per day. The excessive daytime sleepiness and "evening" chronotype were associated with fracture risk in observational study, but the association between chronotype and fracture did not show in MR analyses. We further generated a sleep risk score (SRS) with potential risk factors (i.e., insomnia, sleep duration, chronotype, and daytime sleepiness). We found that the risk of fracture increased with an increasing SRS (HR = 1.087, 95% CI = 1.065-1.111, P = 1.27 × 10-14). Moreover, 17.4% of the fracture cases would be removed if all participants exhibited a healthy sleep pattern. In conclusion, insomnia had a causal effect on fracture, falls had a larger intermediary effect than BMD in this association. Individuals with fracture risk could benefit from the intervention on unhealthy sleep pattern.
Subject(s)
Fractures, Bone/epidemiology , Sleep Wake Disorders/epidemiology , Sleep , Adult , Aged , Female , Fractures, Bone/etiology , Fractures, Bone/genetics , Humans , Incidence , Male , Middle Aged , Sleep Wake Disorders/etiology , Sleep Wake Disorders/genetics , United Kingdom/epidemiologyABSTRACT
Piao chicken, a rare Chinese native poultry breed, lacks primary tail structures, such as pygostyle, caudal vertebra, uropygial gland, and tail feathers. So far, the molecular mechanisms underlying tail absence in this breed remain unclear. In this study, we comprehensively employed comparative transcriptomic and genomic analyses to unravel potential genetic underpinnings of rumplessness in Piao chicken. Our results reveal many biological factors involved in tail development and several genomic regions under strong positive selection in this breed. These regions contain candidate genes associated with rumplessness, including Irx4, Il18, Hspb2, and Cryab. Retrieval of quantitative trait loci (QTL) and gene functions implies that rumplessness might be consciously or unconsciously selected along with the high-yield traits in Piao chicken. We hypothesize that strong selection pressures on regulatory elements might lead to changes in gene activity in mesenchymal stem cells of the tail bud. The ectopic activity could eventually result in tail truncation by impeding differentiation and proliferation of the stem cells. Our study provides fundamental insights into early initiation and genetic basis of the rumpless phenotype in Piao chicken.
Subject(s)
Chickens , Transcriptome , Animals , Chickens/genetics , Genomics/methods , Phenotype , Polymorphism, Single Nucleotide , Quantitative Trait LociABSTRACT
The domestic dog, as a highly successful domestication model, is well known as a favored human companion. Exploring its domestication history should provide great insight into our understanding of the prehistoric development of human culture and productivity. Furthermore, investigation on the mechanisms underpinning the morphological and behavioral traits associated with canid domestication syndrome is of significance not only for scientific study but also for human medical research. Current development of a multidisciplinary canine genome database, which includes enormous omics data, has substantially improved our understanding of the genetic makeup of dogs. Here, we reviewed recent advances associated with the original history and genetic basis underlying environmental adaptations and phenotypic diversities in domestic dogs, which should provide perspectives on improving the communicative relationship between dogs and humans.
Subject(s)
Adaptation, Biological , Dogs , Domestication , Genetic Variation , Phenotype , AnimalsABSTRACT
Genomic imprinting often results in parent-of-origin specific differential expression of maternally and paternally inherited alleles and plays an essential role in mammalian development and growth. Mammalian genomic imprinting has primarily been studied in mice and humans, with only limited information available for pigs. To systematically characterize this phenomenon and evaluate imprinting status between different species, we investigated imprinted genes on a genome-wide scale in pig brain tissues. Specifically, we performed bioinformatics analysis of high-throughput sequencing results from parental genomes and offspring transcriptomes of hybrid crosses between Duroc and Diannan small-ear pigs. We identified 11 paternally and five maternally expressed imprinted genes in pigs with highly stringent selection criteria. Additionally, we found that the KCNQ1 and IGF2R genes, which are related to development, displayed a different imprinting status in pigs compared with that in mice and humans. This comprehensive research should help improve our knowledge on genomic imprinting in pigs and highlight the potential use of imprinted genes in the pig breeding field.
Subject(s)
Genomic Imprinting , Mammals/genetics , Swine/genetics , Alleles , Animals , Genome-Wide Association Study , Species SpecificityABSTRACT
The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East-Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.