Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Association of seafood consumption with cardiovascular disease among adults in Qingdao, China.

BACKGROUND AND AIMS: The relationship between seafood consumption and cardiovascular disease (CVD) is controversial, and studies have not considered competing risk events. Our study examined the association between a full range of seafood consumption and CVD incidence and mortality based on the Qingdao Diabetes Prevention Program. METHODS AND RESULTS: We followed up 5285 participants without CVD at baseline until December 31, 2021. CVD cases and deaths were identified through record linkage with the Qingdao CVD Surveillance System and the Qingdao Death Surveillance System, respectively. Information on seafood consumption was obtained using a food frequency questionnaire. We used the Cox proportional hazard model and the competing risk model to evaluate the association between all types of seafood consumption and CVD incidence and mortality. During a median follow-up of 11.4 years, 122 CVD cases and 75 deaths occurred. After adjustment for potential confounders, compared with nonconsumers, seafood consumption of 300-500 and > 500 g/week was associated with a lower risk of CVD incidence [hazards ratio and 95 % confidence interval (CI): 0.54 (0.29-0.99) and 0.49 (0.26-0.91), respectively]. However, seafood consumption of >500 g/week had a significantly lower risk of CVD mortality [subdistribution hazard ratio and 95 % CI: 0.40 (0.17-0.95)], but it was insignificant in other groups. CONCLUSION: Seafood consumption of 300-500 g/week and >500 g/week was associated with a lower CVD incidence and mortality. Our findings provide evidence of the recommendations of the 2022 Dietary Guidelines for Chinese residents and may guide the promotion of strategies for CVD prevention.

Accuracy of subtraction fractional flow reserve with computed tomography in identifying early revascularization in patients with coronary artery disease.

OBJECTIVES: The diagnostic performance of fractional flow reserve with computed tomography (FFR-CT) is affected by the presence of calcified plaque. Subtraction can remove the influence of calcification in coronary computed tomography angiography (CCTA) to increase confidence in the diagnosis of coronary artery stenosis. Our purpose is to investigate the accuracy of post-subtraction FFR-CT in predicting early revascularization. DESIGN: Based on CCTA data of 237 vessels from 79 patients with coronary artery disease, subtraction CCTA images were obtained at a local post-processing workstation, and the conventional and post-subtraction FFR-CT measurements and the difference in proximal and distal FFR-CT values of the narrowest segment of the vessel (ΔFFR-CT) were analyzed for their accuracy in predicting early coronary artery hemodynamic reconstruction. RESULTS: With FFR-CT ≤ 0.8 as the criterion, the accuracy of conventional and post-subtraction FFR-CT measurements in predicting early revascularization was 73.4% and 77.2% at the patient level, and 64.6% and 72.2% at the vessel level, respectively. The specificity of post-subtraction FFR-CT measurements was significantly higher than that of conventional FFR-CT at both the patient and vessel levels (P of 0.013 and 0.015, respectively). At the vessel level, the area under the curve of receiver operating characteristic was 0.712 and 0.797 for conventional and post-subtraction ΔFFR-CT, respectively, showing a difference (P = 0.047), with optimal cutoff values of 0.07 and 0.11, respectively. CONCLUSION: The post-subtraction FFR-CT measurements enhance the specificity in predicting early revascularization. The post-subtraction ΔFFR-CT value of the stenosis segment > 0.11 may be an important indicator for early revascularization.

Association between famine exposure during infancy and childhood and the risk of chronic kidney disease in adulthood.

BACKGROUND: Famine exposure in childhood is proven to be associated with multiple chornic disease in adult but has not been studied with chronic kidney disease (CKD). AIMS: This study was conducted to identify the relationship between famine exposure during infancy and childhood - specifically, the Chinese famine of 1959-1961 - and the risk of adult-onset chronic kidney disease (CKD) among Chinese individuals. METHODS: This study included 2937 individuals from the Qingdao Diabetes Prevention Program. They were stratified by birth year into infancy-exposed (1956-1958), childhood-exposed (1950-1955) and unexposed (1963-1971) groups. The estimated glomerular filtration rate (eGFR) was calculated using the Chronic Kidney Disease Epidemiology Collaboration equation. CKD was defined as an eGFR of <90 mL/min/1.73 m2. RESULTS: The mean eGFR values for the infancy-exposed and childhood-exposed groups were 107.23 ± 12.53 and 103.23 ± 12.44 mL/min/1.73 m2, respectively, both of which were lower than that of the unexposed group (114.82 ± 13.39 mL/min/1.73 m2; P < 0.05). In the crude model, the odds ratio (OR) for CKD was 2.00 (95% confidence interval (CI): 1.39-2.88) in the infancy-exposed group and 2.92 (95% CI: 2.17-3.93) in the childhood-exposed group. Further adjustments for urban/rural residence, body mass index, age, current smoking, type 2 diabetes, systolic blood pressure, diastolic blood pressure and total cholesterol did not significantly alter the association between famine exposure and CKD. The corresponding ORs were 1.71 (95% CI: 1.17-2.50) and 2.48 (95% CI: 1.81-3.40) for the infancy-exposed and childhood-exposed groups respectively. CONCLUSIONS: Famine exposure during infancy and childhood is associated with a long-term decline in eGFR and an increased adult-onset CKD risk. Early intervention for high-risk individuals may mitigate the risk of adult-onset CKD.

Taraxasterol attenuates zearalenone-induced kidney damage in mice by modulating oxidative stress and endoplasmic reticulum stress.

Taraxasterol is one of the bioactive ingredients from traditional Chinese herb Taraxacum, which exhibits multiple pharmacological activities and protective effects. However, the underlying influence and mechanism of its use against kidney damage caused from zearalenone (ZEA) remain unexplored. The ZEA-induced kidney damage model of mice was established by feeding diets containing ZEA (2 mg/kg), and taraxasterol (5 and 10 mg/kg) was administered by gavage for 28 days. Results demonstrated taraxasterol increased average daily gain (ADG) and average daily feed intake (ADFI), reduced feed-to-gain ratio (F/G) and kidney index of mice induced by ZEA. Taraxasterol alleviated histopathological changes of kidney, reduced ZEA residue and the levels of blood urea nitrogen (BUN), uric acid (UA), and creatinine (CRE). Concurrently, taraxasterol reduced the contents of oxidative stress indicator reactive oxygen species (ROS) and malondialdehyde (MDA), and increased the activities of antioxidant enzymes catalase (CAT), total superoxide dismutase (T-SOD), and glutathione peroxidase (GSH-Px). Further, taraxasterol up-regulated the mRNA and protein expression of nuclear factor erythroid-2-related factor 2 (Nrf2), GSH-Px, NAD(P)H quinone oxidoreductase 1 (NQO1), and heme oxygenase-1 (HO-1), and down-regulated the mRNA and protein expression of KELCH like ECH associated protein (Keap1) in Nrf2/Keap1 pathway. Taraxasterol down-regulated the mRNA and protein expression of immunoglobulin binding protein (Bip), C/EBP homologous protein (CHOP), Bcl-2 associated X (Bax), cysteine protease (Caspase)-12, and Caspase-3, and up-regulated B-cell lymphoma 2 (Bcl-2) expression in endoplasmic reticulum stress pathway. This study suggests that taraxasterol attenuates ZEA-induced mouse kidney damage through the modulation of Nrf2/Keapl pathway to play antioxidant role and endoplasmic reticulum stress pathway to enhance anti-apoptotic ability. It will provide a basis for taraxasterol as a potential drug to prevent and treat ZEA-induced kidney damage.

Does obesity mediate the relationship between diabetes and osteoporosis in Chinese elderly population?

OBJECTIVES: Diabetes mellitus (DM), osteoporosis (OP), and obesity (OB) are three complex diseases. OB is associated with both DM and OP, but it is unclear whether OB mediates association between DM and OP. The study aimed to investigate the potential mediation effects of OB on association between DM and bone mineral density (BMD) by the causal inference tests (CIT). METHODS: A total of 5682 Chinese aged over 65 years were enrolled in an ongoing cohort: Osteoporosis Preventive Project (OPP). Obesity-related indexes, including body mass index (BMI), waist circumference, and waist circumference-hip circumference-ratio (WHR), and BMD at total hip (TH) and femur neck (FN) were measured. RESULTS: Subjects with DM had significant greater values of age, weight, BMI, waist circumference, WHR, and BMD than non-DM subjects. BMD at TH and FN was significantly associated with DM (p < 0.05) with adjustment of age both in males and females. Further CIT showed that OB-related indexes (BMI, waist circumference, and WHR) are significantly mediators in the associations between DM and BMD in females, but not in males. Furthermore, the mediation effects of waist circumference were detected on DM and TH BMD in the females of normal-weight group. CONCLUSIONS: Obesity-related indexes, especially waist circumference, serve as significant mediator(s) between DM and OP in Chinese female elderly. Diabetes increases BMD by increasing obesity-related indexes. The findings established the intermediate role of OB underlying the association between DM and OP in human population.

The obesity indices mediate the relationships of blood lipids and bone mineral density in Chinese elders.

OBJECTIVE: It is unclear what role does obesity (OB) index play between blood lipid and bone mineral density (BMD). SUBJECTS: This study recruited a total of 4,558 Chinese elders >65 years. OB indices: waist circumference (WC), body mass index (BMI), waist-hip-ratio (WHR); blood lipid parameters: low density lipoprotein (LDL); total cholesterol (TC), triglyceride (TG), and BMDs at femur neck (FN), total hip (TH), and lumbar spine (LS) were measured. The t-test and multiple linear regression analysis were used to detect the differences of variables. Casual inference test (CIT) were performed to test potential mediators underlying the associations between blood lipid and BMD. RESULTS: The blood lipids were positively associated with BMD (p < 0.05) after adjustment of age and sex (Model 1) both in total subjects and in sex-stratified subjects. The CIT showed that OB indices had significant mediation effects on the associations between blood lipid (TG and LDL) and BMD in total subjects and males. Comparably, the correlations of TG and BMD are most likely mediated by BMI and WC. CONCLUSIONS: This study represented the first effort to report that OB indices, especially BMI and WC, served as significant mediators between blood lipid (TG and LDL) and BMD in Chinese elderly.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Newly developed hip geometry parameters are associated with hip fracture.

BACKGROUND: The conventional and newly developed geometry parameters at the femoral neck have formed a large geometry profile and their relationship with hip fracture was largely unknown. We aimed to evaluate the relationship between the geometry profile and hip fracture in Chinese. METHODS: The hip geometry profile contains seven conventional geometry parameters at femoral neck (FN) and thirty newly developed parameters at three sub-regions (Narrow Neck, NN; Intertrochanter, IT; Femoral shaft, FS) of the total hip. Based on 6294 recruited Chinese (≥65 years), 97 subjects with osteoporotic fracture (OF) history and 388 matched controls were selected. The t test, Chi-square statistic, conditional logistic regression model were used. RESULTS: Three geometric parameters (endocortical diameter, ED; cortical thickness, CT; buckling ratio, BR) have consistent differences at all the sites between the cases and controls (p < 0.01). Conventional geometry parameters (e.g., cross-sectional area, CSA; BR) and the newly developed parameters (e.g., NN_ED, NN_Outer Diameter, IT_ED) were identified as the risk factors of hip fracture independent of BMD. The additional predictive ability of the hip geometric parameters, over BMD alone, (receiver operating characteristic curve (ROC) analysis) was improved. Especially at NN, the area of ROC used single NN_BMD is only 0.596, but increased rapidly at 0.705 when combined with the hip geometric parameters. CONCLUSIONS: This study found that three newly developed hip geometry parameters are associated with hip fracture. The results will increase our understanding of the determinants of fracture and provide potential clue for future prevention of fracture in Chinese Population.

Association of Plasma Irisin with Bone Mineral Density in a Large Chinese Population Using an Extreme Sampling Design.

Irisin, a myokine produced by skeletal muscle in response to physical exercise, promotes trans-differentiation of white adipose tissue into brown adipose tissue. Recent evidences suggested that irisin also plays an important role in the control of bone metabolism. This study aimed to ascertain the relationship between plasma irisin and bone mineral density (BMD) in Chinese population by adoption of an extreme sampling method. Based on a large and screened Chinese elderly population (N = 6308), two subgroups with extremely high and low hip BMD were selected for discovery (N = 80, high vs. low BMD = 44:36) and validation (N = 60, high vs. low BMD = 30:30), respectively. Plasma irisin, P1NP, and ß-CTx were measured using commercially available ELISA kits. Other metabolic parameters (e.g., blood glucose, total cholesterol and triglycerides) were collected. Student's t test and Spearman correlation analyses were conducted in SPSS. Significant difference was discovered for plasma irisin between females and age-matched males (N = 80, male vs. female = 42:38, P = 0.002). The plasma irisin levels were significantly higher in high BMD subjects than in low BMD subjects, which was observed in both discovery (P = 0.012) and validation samples (P = 0.022). However, such observation was limited to males only. Further correlation analyses in males showed that plasma irisin was correlated with BMD (r = 0.362, P = 0.025) and triglyceride (r = - 0.354, P = 0.032). Plasma irisin levels were associated with hip BMD in Chinese elderly men. This study represented the first effort of investigating the relationship of plasma irisin and BMD in elderly population. The positive correlation between plasma irisin and BMD hints intrinsic communication between muscle and bone.

Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: ß = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Genetic risk factors for decreased bone mineral accretion in children with asthma receiving multiple oral corticosteroid bursts.

BACKGROUND: Long-term intermittent oral corticosteroid (OCS) use in children with asthma leads to significant decreases in bone mineral accretion (BMA). OBJECTIVE: We aimed to identify genetic factors influencing OCS dose effects on BMA in children with asthma. METHODS: We first performed a gene-by-OCS interaction genome-wide association study (GWAS) of BMA in 489 white participants in the Childhood Asthma Management Program trial who took short-term oral prednisone bursts when they experienced acute asthma exacerbations. We selected the top-ranked 2000 single nucleotide polymorphisms (SNPs) in the GWAS and determined whether these SNPs also had cis-regulatory effects on dexamethasone-induced gene expression in osteoblasts. RESULTS: We identified 2 SNPs (rs9896933 and rs2074439) associated with decreased BMA and related to the tubulin γ pathway. The rs9896933 variant met the criteria for genome-wide significance (P = 3.15 × 10(-8) in the GWAS) and is located on the intron of tubulin folding cofactor D (TBCD) gene. The rs2074439 variant (P = 2.74 × 10(-4) in the GWAS) showed strong cis-regulatory effects on dexamethasone-induced tubulin γ gene expression in osteoblasts (P = 8.64 × 10(-4)). Interestingly, we found that BMA worsened with increasing prednisone dose as the number of mutant alleles of the 2 SNPs increased. CONCLUSIONS: We have identified 2 novel tubulin γ pathway SNPs, rs9896933 and rs2074439, showing independent interactive effects with cumulative corticosteroid dose on BMA in children with asthma receiving multiple OCS bursts.

Allelic expression mapping across cellular lineages to establish impact of non-coding SNPs.

Most complex disease-associated genetic variants are located in non-coding regions and are therefore thought to be regulatory in nature. Association mapping of differential allelic expression (AE) is a powerful method to identify SNPs with direct cis-regulatory impact (cis-rSNPs). We used AE mapping to identify cis-rSNPs regulating gene expression in 55 and 63 HapMap lymphoblastoid cell lines from a Caucasian and an African population, respectively, 70 fibroblast cell lines, and 188 purified monocyte samples and found 40-60% of these cis-rSNPs to be shared across cell types. We uncover a new class of cis-rSNPs, which disrupt footprint-derived de novo motifs that are predominantly bound by repressive factors and are implicated in disease susceptibility through overlaps with GWAS SNPs. Finally, we provide the proof-of-principle for a new approach for genome-wide functional validation of transcription factor-SNP interactions. By perturbing NFκB action in lymphoblasts, we identified 489 cis-regulated transcripts with altered AE after NFκB perturbation. Altogether, we perform a comprehensive analysis of cis-variation in four cell populations and provide new tools for the identification of functional variants associated to complex diseases.

Catecholic amides as potential selective phosphodiesterase 4D inhibitors: Design, synthesis, pharmacological evaluation and structure-activity relationships.

In this study, a series of catechol-based amides (8a-n) with different amide linkers linking the catecholic moiety to the terminal phenyl ring was designed and synthesized as potent phosphodiesterase (PDE) 4D inhibitors. The inhibitory activities of these compounds were evaluated against the core catalytic domains of human PDE4 (PDE4CAT), full-length PDE4B1 and PDE4D7 enzymes, and other PDE family members. The results indicated the majority of compounds 8a-n displayed moderate to good inhibitory activities against PDE4CAT. Among these compounds, compound 8 j with a short amide linker (-CONHCH2-) displayed comparable PDE4CAT inhibitory activity (IC50=410 nM) with rolipram. More interestingly, compound 8 g, a potent and selective PDE4D inhibitor (IC50=94 nM), exhibited a 10-fold selectivity over the PDE4B subtypes and an over 1000-fold selectivity against other PDE family members. Docking simulations suggested that 8 g forms three extra H-bonds with the N-H of residue Asn487 and two water molecules.

Global analysis of the impact of environmental perturbation on cis-regulation of gene expression.

Genetic variants altering cis-regulation of normal gene expression (cis-eQTLs) have been extensively mapped in human cells and tissues, but the extent by which controlled, environmental perturbation influences cis-eQTLs is unclear. We carried out large-scale induction experiments using primary human bone cells derived from unrelated donors of Swedish origin treated with 18 different stimuli (7 treatments and 2 controls, each assessed at 2 time points). The treatments with the largest impact on the transcriptome, verified on two independent expression arrays, included BMP-2 (t = 2h), dexamethasone (DEX) (t = 24 h), and PGE2 (t = 24 h). Using these treatments and control, we performed expression profiling for 18,144 RefSeq transcripts on biological replicates of the complete study cohort of 113 individuals (n(total) = 782) and combined it with genome-wide SNP-genotyping data in order to map treatment-specific cis-eQTLs (defined as SNPs located within the gene ± 250 kb). We found that 93% of cis-eQTLs at 1% FDR were observed in at least one additional treatment, and in fact, on average, only 1.4% of the cis-eQTLs were considered as treatment-specific at high confidence. The relative invariability of cis-regulation following perturbation was reiterated independently by genome-wide allelic expression tests where only a small proportion of variance could be attributed to treatment. Treatment-specific cis-regulatory effects were, however, 2- to 6-fold more abundant among differently expressed genes upon treatment. We further followed-up and validated the DEX-specific cis-regulation of the MYO6 and TNC loci and found top cis-regulatory variants located 180 kb and 250 kb upstream of the transcription start sites, respectively. Our results suggest that, as opposed to tissue-specificity of cis-eQTLs, the interactions between cellular environment and cis-variants are relatively rare (∼1.5%), but that detection of such specific interactions can be achieved by a combination of functional genomic approaches as described here.

Identification, replication, and fine-mapping of Loci associated with adult height in individuals of african ancestry.

Adult height is a classic polygenic trait of high heritability (h(2) approximately 0.8). More than 180 single nucleotide polymorphisms (SNPs), identified mostly in populations of European descent, are associated with height. These variants convey modest effects and explain approximately10% of the variance in height. Discovery efforts in other populations, while limited, have revealed loci for height not previously implicated in individuals of European ancestry. Here, we performed a meta-analysis of genome-wide association (GWA) results for adult height in 20,427 individuals of African ancestry with replication in up to 16,436 African Americans. We found two novel height loci (Xp22-rs12393627, P = 3.4×10(-12) and 2p14-rs4315565, P = 1.2×10(-8)). As a group, height associations discovered in European-ancestry samples replicate in individuals of African ancestry (P = 1.7×10(-4) for overall replication). Fine-mapping of the European height loci in African-ancestry individuals showed an enrichment of SNPs that are associated with expression of nearby genes when compared to the index European height SNPs (P<0.01). Our results highlight the utility of genetic studies in non-European populations to understand the etiology of complex human diseases and traits.

Short-term lifestyle education on obesity reduction in adolescents.

Backgrounds: Obesity is increasing in adolescents in China. However, the awareness of obesity and prevention on related risk factors were not well known. We aim to assess the effectiveness of short-term health education intervention on obesity in Chinese adolescents. Methods: In this study, 42 primary and secondary schools from Qingdao were randomly divided into the education and control groups. A total of 11,739 adolescents was included in the current study. The logistic regression was employed to assess odds ratio (OR) of education intervention on overweight and obesity prevalence adjusting for covariates. Results: The baseline prevalence of overweight and obesity was significantly higher in urban than in rural areas and in boys than in girls. After 1 year lifestyle intervention, the proportion of students with awareness of obesity was higher, meanwhile age-adjusted mean values of weight, body mass index, duration of watching TV and doing homework were lower in education group than control group. The corresponding figures were 43.6 [95% CI (confidence intervals); 43.3-43.9] kg versus 44.3 (95% CI; 44.0-44.6) kg, 18.6 (95% CI; 18.5-18.7) kg/m2 versus 18.9 (95% CI; 18.8-19.1) kg/m2, 1.3 (95% CI; 1.2-1.3) hours/d versus 1.4 (95% CI; 1.3-1.4) hours/d, and 1.5 (95% CI; 1.4-1.5) hours/d versus 1.8 (95% CI, 1.7-1.8) hours/d. The multivariable adjusted OR for combined prevalence of overweight and obesity was 0.85 (95% CI, 0.76-0.96) in education group as compared with control group. Conclusion: Short-term health education intervention results in significantly higher reductions in obesity parameters and improvement in awareness in Chinese adolescents.

Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores.

Recent studies have revealed the pervasive landscape of rare structural variants (rSVs) present in human genomes. rSVs can have extreme effects on the expression of proximal genes and, in a rare disease context, have been implicated in patient cases where no diagnostic single nucleotide variant (SNV) was found. Approaches for integrating rSVs to date have focused on targeted approaches in known Mendelian rare disease genes. This approach is intractable for rare diseases with many causal loci or patients with complex, multi-phenotype syndromes. We hypothesized that integrating trait-relevant polygenic scores (PGS) would provide a substantial reduction in the number of candidate disease genes in which to assess rSV effects. We further implemented a method for ranking PGS genes to define a set of core/key genes where a rSV has the potential to exert relatively larger effects on disease risk. Among a subset of patients enrolled in the Genomic Answers for Kids (GA4K) rare disease program (N=497), we used PacBio HiFi long-read whole genome sequencing (lrWGS) to identify rSVs intersecting genes in trait-relevant PGSs. Illustrating our approach in Autism (N=54 cases), we identified 22, 019 deletions, 2,041 duplications, 87,826 insertions, and 214 inversions overlapping putative core/key PGS genes. Additionally, by integrating genomic constraint annotations from gnomAD, we observed that rare duplications overlapping putative core/key PGS genes were frequently in higher constraint regions compared to controls (P = 1×10-03). This difference was not observed in the lowest-ranked gene set (P = 0.15). Overall, our study provides a framework for the annotation of long-read rSVs from lrWGS data and prioritization of disease-linked genomic regions for downstream functional validation of rSV impacts. To enable reuse by other researchers, we have made SV allele frequencies and gene associations freely available.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - in the onset and phenotypic presentation of rare diseases. Here, we comprehensively map individual polygenic liability for 1102 open-source PGS in a cohort of 3059 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. Using this resource, we demonstrate increased polygenic liability in probands with an inherited candidate disease variant (VUS) compared to unaffected carrier parents. Further, we show an enrichment for large-effect rare variants in putative core PGS genes for associated complex traits. Overall, our study supports and expands on previous findings of complex trait associations in rare diseases, implicates polygenic liability as a potential mechanism underlying variable penetrance of candidate causal variants, and provides a framework for identifying novel candidate rare disease genes.