Genes linked to energy metabolism and immunoregulatory mechanisms are associated with subcutaneous adipose tissue distribution in HIV-infected men

Admixture/fine-mapping in Brazilians reveals a West African associated potential regulatory variant (rs114066381) with a strong female-specific effect on body mass and fat mass indexes

BACKGROUND/OBJECTIVES

Admixed populations are a resource to study the global genetic architecture of complex phenotypes, which is critical, considering that non-European populations are severely underrepresented in genomic studies. Here, we study the genetic architecture of BMI in children, young adults, and elderly individuals from the admixed population of Brazil.

SUBJECTS/METHODS

Leveraging admixture in Brazilians, whose chromosomes are mosaics of fragments of Native American, European, and African origins, we used genome-wide data to perform admixture mapping/fine-mapping of body mass index (BMI) in three Brazilian population-based cohorts from Northeast (Salvador), Southeast (Bambuí), and South (Pelotas).

RESULTS

We found significant associations with African-associated alleles in children from Salvador (PALD1 and ZMIZ1 genes), and in young adults from Pelotas (NOD2 and MTUS2 genes). More importantly, in Pelotas, rs114066381, mapped in a potential regulatory region, is significantly associated only in females (p = 2.76e-06). This variant is rare in Europeans but with frequencies of ~3% in West Africa and has a strong female-specific effect (95% CI: 2.32-5.65 kg/m² per each A allele). We confirmed this sex-specific association and replicated its strong effect for an adjusted fat mass index in the same Pelotas cohort, and for BMI in another Brazilian cohort from São Paulo (Southeast Brazil). A meta-analysis confirmed the significant association. Remarkably, we observed that while the frequency of rs114066381-A allele ranges from 0.8 to 2.1% in the studied populations, it attains ~9% among women with morbid obesity from Pelotas, São Paulo, and Bambuí. The effect size of rs114066381 is at least five times higher than the FTO SNPs rs9939609 and rs1558902, already emblematic for their high effects.

CONCLUSIONS

We identified six candidate SNPs associated with BMI. rs114066381 stands out for its high effect that was replicated and its high frequency in women with morbid obesity. We demonstrate how admixed populations are a source of new relevant phenotype-associated genetic variants.

Comparison of Growth Performance and Meat Quality Traits of Commercial Cross-Bred Pigs versus the Large Black Pig Breed

The meat quality of different pig breeds is associated with their different muscle tissue physiological processes, which involves a large variety of genes related with muscle fat and energy metabolism. Understanding the differences of biological processes of muscle after slaughter is helpful to reveal the meat quality development of different breeds. Therefore, eight native Large Black pigs (BP), with high fat content in meat, and seven cross-bred commercial pigs (CP), which had a high feed efficiency with high lean meat, were used to investigate the differences in their meat quality and RNA transcriptomes. The average daily gain (ADG) and hot carcass weight (HCW) of CP were higher than BP, but the back-fat thickness of BP was higher than CP (p < 0.05). The CP had higher a* (redness) but lower h (hue angle) than BP (p < 0.05). The metmyoglobin (MMb) percentage of CP was higher (p < 0.05) than BP. The fat content and oxygen consumption of longissimus dorsi (LD) muscles in BP were higher (p < 0.05) than CP. BP had higher monounsaturated fatty acids (MUFA) content, but CP had higher polyunsaturated fatty acids (PUFA) content (p < 0.05). The RNA-seq data highlighted 201 genes differentially expressed between the two groups (corrected false discovery rate (FDR) p < 0.05), with 75 up-regulated and 126 down-regulated genes in BP compared with CP using the fold change (FC). The real-time PCR was used to validate the results of RNA-seq for eight genes, and the genes related to lipid and energy metabolism were highly expressed in BP (p < 0.05). Based on the results, BP had superior intramuscular fat content to CP, while the growth performance of CP was better, and the transcriptomic differences between these two groups of pigs may cause the meat quality and growth performance variance.

Genetic architecture of cardiometabolic risks in people living with HIV

BACKGROUND

Advances in antiretroviral therapies have greatly improved the survival of people living with human immunodeficiency virus (HIV) infection (PLWH); yet, PLWH have a higher risk of cardiovascular disease than those without HIV. While numerous genetic loci have been linked to cardiometabolic risk in the general population, genetic predictors of the excessive risk in PLWH are largely unknown.

METHODS

We screened for common and HIV-specific genetic variants associated with variation in lipid levels in 6284 PLWH (3095 European Americans [EA] and 3189 African Americans [AA]), from the Centers for AIDS Research Network of Integrated Clinical Systems cohort. Genetic hits found exclusively in the PLWH cohort were tested for association with other traits. We then assessed the predictive value of a series of polygenic risk scores (PRS) recapitulating the genetic burden for lipid levels, type 2 diabetes (T2D), and myocardial infarction (MI) in EA and AA PLWH.

RESULTS

We confirmed the impact of previously reported lipid-related susceptibility loci in PLWH. Furthermore, we identified PLWH-specific variants in genes involved in immune cell regulation and previously linked to HIV control, body composition, smoking, and alcohol consumption. Moreover, PLWH at the top of European-based PRS for T2D distribution demonstrated a > 2-fold increased risk of T2D compared to the remaining 95% in EA PLWH but to a much lesser degree in AA. Importantly, while PRS for MI was not predictive of MI risk in AA PLWH, multiethnic PRS significantly improved risk stratification for T2D and MI.

CONCLUSIONS

Our findings suggest that genetic loci involved in the regulation of the immune system and predisposition to risky behaviors contribute to dyslipidemia in the presence of HIV infection. Moreover, we demonstrate the utility of the European-based and multiethnic PRS for stratification of PLWH at a high risk of cardiometabolic diseases who may benefit from preventive therapies.

Minding the gap in HIV host genetics: opportunities and challenges

Genome-wide association studies (GWAS) have been successful in identifying and confirming novel genetic variants that are associated with diverse HIV phenotypes. However, these studies have predominantly focused on European cohorts. HLA molecules have been consistently associated with HIV outcomes, some of which have been found to be population specific, underscoring the need for diversity in GWAS. Recently, there has been a concerted effort to address this gap that leads to health care (disease prevention, diagnosis, treatment) disparities with marginal improvement. As precision medicine becomes more utilized, non-European individuals will be more and more disadvantaged, as the genetic variants identified in genomic research based on European populations may not accurately reflect that of non-European individuals. Leveraging pre-existing, large, multiethnic cohorts, such as the UK Biobank, 23andMe, and the National Institute of Health's All of Us Research Program, can contribute in raising genomic research in non-European populations and ultimately lead to better health outcomes.

BMPR2 promotes fatty acid oxidation and protects white adipocytes from cell death in mice

Adipocyte cell death is pathologically involved in both obesity and lipodystrophy. Inflammation and pro-inflammatory cytokines are generally regarded as inducers for adipocyte apoptosis, but whether some innate defects affect their susceptibility to cell death has not been extensively studied. Here, we found bone morphogenetic protein receptor type 2 (BMPR2) knockout adipocytes were prone to cell death, which involved both apoptosis and pyroptosis. BMPR2 deficiency in adipocytes inhibited phosphorylation of perilipin, a lipid-droplet-coating protein, and impaired lipolysis when stimulated by tumor necrosis factor (TNFα), which lead to failure of fatty acid oxidation and oxidative phosphorylation. In addition, impaired lipolysis was associated with mitochondria-mediated apoptosis and pyroptosis as well as elevated inflammation. These results suggest that BMPR2 is important for maintaining the functional integrity of adipocytes and their ability to survive when interacting with inflammatory factors, which may explain why adipocytes among individuals show discrepancy for death responses in inflammatory settings.

Qian, Pan, Su et al. show that bone morphogenetic protein receptor type 2 (BMPR2) deficiency makes white adipocytes prone to death. This study suggests that BMPR2 is required for fat cells to effectively generate energy from the stored fat for their survival during energetically demanding situations such as inflammation.

Exploratory analysis of genetic variants influencing molecular traits in cerebral cortex of suicide completers

Genetic factors have been implicated in suicidal behavior. It has been suggested that one of the roles of genetic factors in suicide could be represented by the effect of genetic variants on gene expression regulation. Alteration in the expression of genes participating in multiple biological systems in the suicidal brain has been demonstrated, so it is imperative to identify genetic variants that could influence gene expression or its regulatory mechanisms. In this study, we integrated DNA methylation, gene expression, and genotype data from the prefrontal cortex of suicides to identify genetic variants that could be factors in the regulation of gene expression, generally called quantitative trait locus (xQTLs). We identify 6,224 methylation quantitative trait loci and 2,239 expression quantitative trait loci (eQTLs) in the prefrontal cortex of suicide completers. The xQTLs identified influence the expression of genes involved in neurodevelopment and cell organization. Two of the eQTLs identified (rs8065311 and rs1019238) were previously associated with cannabis dependence, highlighting a candidate genetic variant for the increased suicide risk in subjects with substance use disorders. Our findings suggest that genetic variants may regulate gene expression in the prefrontal cortex of suicides through the modulation of promoter and enhancer activity, and to a lesser extent, binding transcription factors.

Important candidate genes for abdominal fat content identified by linkage disequilibrium and fixation index information

Selection for rapid growth in chickens has always been accompanied by increased fat deposition and excessive fat deposition, especially abdominal fat, cannot only decrease feed efficiency but also cause many diseases. Finding the candidate genes associated with abdominal fat deposition is essential for breeding. To identify these candidate genes, we applied linkage disequilibrium and selection signature analysis using chicken 60 k single nucleotide polymorphism (SNP) chips in two broiler lines divergently selected for abdominal fat content for 11 generations. After quality control, 46,033 SNPs were left for analysis. Using these SNPs, we found that r2 was 0.06 to 0.14 in the lean line and 0.07 to 0.13 in the fat line for all 28 chromosomes (except GGA16). Pairwise SNP distances <25 kb showed a mean r2 = 0.33 in the lean line and r2 = 0.32 in the fat line. The fixation index (FST) analysis was carried out and 46 SNPs with the top 0.1% of the FST value was detected as the loci with selection signatures. Besides FST, hapFLK was also used to detect selection signatures for abdominal fat content. A total of 11 genes, including transient receptor potential cation channel subfamily C member 4, estrogen related receptor gamma, fibroblast growth factor 13, G-protein-signaling modulator 2, RAR related orphan receptor A, phospholipase A2 group X, mitochondrial ribosomal protein L28, metadherin, calcitonin receptor like receptor, serine/threonine kinase 39, and nuclear factor I A, were detected as the important candidate genes for abdominal fat deposition based on their basic functions. The results of the present study may benefit the understanding of genetic mechanism of abdominal fat deposition in chicken.

Cord blood DNA methylation and adiposity measures in early and mid-childhood

Background

Excess adiposity in childhood is associated with numerous adverse health outcomes. As this condition is difficult to treat once present, identification of risk early in life can help inform and implement strategies to prevent the onset of the condition. We performed an epigenome-wide association study to prospectively investigate the relationship between cord blood DNA methylation and adiposity measurements in childhood.

Methods

We measured genome-wide DNA methylation from 478 children in cord blood and measured overall and central adiposity via skinfold caliper measurements in early (range 3.1–3.3 years) and mid-childhood (age range 7.3–8.3 years) and via dual X-ray absorptiometry (DXA) in mid-childhood. Final models were adjusted for maternal age at enrollment, pre-pregnancy body mass index, education, folate intake during pregnancy, smoking during pregnancy, and gestational weight gain, and child sex, race/ethnicity, current age, and cord blood cell composition.

Results

We identified four promoter proximal CpG sites that were associated with adiposity as measured by subscapular (SS) and triceps (TR) ratio (SS:TR) in early childhood, in the genes KPRP, SCL9A10, MYLK2, and PRLHR. We additionally identified one gene body CpG site associated with early childhood SS + TR on PPAPDC1A; this site was nominally associated with SS + TR in mid-childhood. Higher methylation at one promoter proximal CpG site in MMP25 was also associated with SS:TR in mid-childhood. In regional analyses, methylation at an exonal region of GFPT2 was positively associated with SS:TR in early childhood. Finally, we identified regions of two long, non-coding RNAs which were associated with SS:TR (LOC100049716) and fat-free mass index (LOC102723493) in mid-childhood.

Conclusion

This analysis identified novel CpG loci associated with adiposity outcomes. However, our results suggest little consistency across the various adiposity outcomes tested, particularly among the more accurate DXA measurements of body composition. We recommend using caution when interpreting these associations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0384-9) contains supplementary material, which is available to authorized users.

Host genetic variation and HIV disease: from mapping to mechanism

This review aims to provide a summary of current knowledge of host genetic effects on human immunodeficiency virus (HIV) disease. Mapping of simple single nucleotide polymorphisms (SNP) has been largely successful in HIV, but more complex genetic associations involving haplotypic or epigenetic variation, for example, remain elusive. Mechanistic insights explaining SNP associations are incomplete, but continue to be forthcoming. The number of robust immunogenetic correlates of HIV is modest and their discovery mostly predates the genome-wide era. Nevertheless, genome-wide evaluations have nicely validated the impact of HLA and CCR5 variants on HIV disease, and importantly, made clear the many false positive associations that were previously suggested by studies using the candidate gene approach. We describe how multiple HIV outcome measures such as acquisition, viral control, and immune decline have been studied in adults and in children, but that collectively these identify only the two replicable loci responsible for modifying HIV disease, CCR5, and HLA. Recent heritability estimates in this disease corroborate the modest impact of genetic determinants and their oligogenic nature. While the mechanism of protection afforded by genetic variants that diminish CCR5 expression is clear, new aspects of HLA class I-mediated protection continue to be uncovered. We describe how these genetic findings have enhanced insights into immunobiology, been clinically translated into CCR5 antagonists, allowed prioritization of antigens for vaccination efforts, and identified targets for genome-editing interventions. Finally, we describe how studies of genetically complex parts of the genome using new tools may begin revealing additional correlates.

Comparative Analyses of QTLs Influencing Obesity and Metabolic Phenotypes in Pigs and Humans

The pig is a well-known animal model used to investigate genetic and mechanistic aspects of human disease biology. They are particularly useful in the context of obesity and metabolic diseases because other widely used models (e.g. mice) do not completely recapitulate key pathophysiological features associated with these diseases in humans. Therefore, we established a F2 pig resource population (n = 564) designed to elucidate the genetics underlying obesity and metabolic phenotypes. Segregation of obesity traits was ensured by using breeds highly divergent with respect to obesity traits in the parental generation. Several obesity and metabolic phenotypes were recorded (n = 35) from birth to slaughter (242 ± 48 days), including body composition determined at about two months of age (63 ± 10 days) via dual-energy x-ray absorptiometry (DXA) scanning. All pigs were genotyped using Illumina Porcine 60k SNP Beadchip and a combined linkage disequilibrium-linkage analysis was used to identify genome-wide significant associations for collected phenotypes. We identified 229 QTLs which associated with adiposity- and metabolic phenotypes at genome-wide significant levels. Subsequently comparative analyses were performed to identify the extent of overlap between previously identified QTLs in both humans and pigs. The combined analysis of a large number of obesity phenotypes has provided insight in the genetic architecture of the molecular mechanisms underlying these traits indicating that QTLs underlying similar phenotypes are clustered in the genome. Our analyses have further confirmed that genetic heterogeneity is an inherent characteristic of obesity traits most likely caused by segregation or fixation of different variants of the individual components belonging to cellular pathways in different populations. Several important genes previously associated to obesity in human studies, along with novel genes were identified. Altogether, this study provides novel insight that may further the current understanding of the molecular mechanisms underlying human obesity.

CNV Analysis Associates AKNAD1 with Type-2 Diabetes in Jordan Subpopulations

Previous studies have identified a number of single nucleotide polymorphisms (SNPs) associated with type-2 diabetes (T2D), but copy number variation (CNV) association has rarely been addressed, especially in populations from Jordan. To investigate CNV associations for T2D in populations in Jordan, we conducted a CNV analysis based on intensity data from genome-wide SNP array, including 34 T2D cases and 110 healthy controls of Chechen ethnicity, as well as 34 T2D cases and 106 healthy controls of Circassian ethnicity. We found a CNV region in protein tyrosine phosphatase receptor type D (PTPRD) with significant association with T2D. PTPRD has been reported to be associated with T2D in genome-wide association studies (GWAS). We additionally identified 16 CNV regions associated with T2D which overlapped with gene exons. Of particular interest, a CNV region in the gene AKNA Domain Containing 1 (AKNAD1) surpassed the experiment-wide significance threshold. Endoplasmic reticulum (ER)-related pathways were significantly enriched among genes which are predicted to be functionally associated with human or mouse homologues of AKNAD1. This is the first CNV analysis of a complex disease in populations of Jordan. We identified and experimentally validated a significant CNVR in gene AKNAD1 associated with T2D.

Differential subcutaneous adipose tissue gene expression patterns in a randomized clinical trial of efavirenz or lopinavir-ritonavir in antiretroviral-naive patients

Gene expression studies of subcutaneous adipose tissue may help to better understand the mechanisms behind body fat changes in HIV-infected patients who initiate antiretroviral therapy (ART). Here, we evaluated early changes in adipose tissue gene expression and their relationship to fat changes in ART-naive HIV-infected patients randomly assigned to initiate therapy with emtricitabine/tenofovir plus efavirenz (EFV) or ritonavir-boosted lopinavir (LPV/r). Patients had abdominal subcutaneous adipose tissue biopsies at baseline and week 16 and dual-energy-X-ray absorptiometry at baseline and weeks 16 and 48. mRNA changes of 11 genes involved in adipogenesis, lipid and glucose metabolism, mitochondrial energy, and inflammation were assessed through reverse transcription-quantitative PCR (RT-qPCR). Additionally, correlations between gene expression changes and fat changes were evaluated. Fat increased preferentially in the trunk with EFV and in the limbs with LPV/r (P < 0.05). After 16 weeks of exposure to the drug regimen, transcripts of CEBP/A, ADIPOQ, GLUT4, LPL, and COXIV were significantly down-regulated in the EFV arm compared to the LPV/r arm (P < 0.05). Significant correlations were observed between LPL expression change and trunk fat change at week 16 in both arms and between CEBP/A or COXIV change and trunk fat change at the same time point only in the EFV arm and not in the LPV/r arm. When combined with emtricitabine/tenofovir as standard backbone therapy, EFV and LPV/r induced differential early expression of genes involved in adipogenesis and energy metabolism. Moreover, these mRNA expression changes correlated with trunk fat change in the EFV arm. (This was a substudy of a randomized clinical trial [LIPOTAR study] registered at ClinicalTrials.gov under identifier NCT00759070.).

A comparison of multivariate genome-wide association methods

Joint association analysis of multiple traits in a genome-wide association study (GWAS), i.e. a multivariate GWAS, offers several advantages over analyzing each trait in a separate GWAS. In this study we directly compared a number of multivariate GWAS methods using simulated data. We focused on six methods that are implemented in the software packages PLINK, SNPTEST, MultiPhen, BIMBAM, PCHAT and TATES, and also compared them to standard univariate GWAS, analysis of the first principal component of the traits, and meta-analysis of univariate results. We simulated data (N = 1000) for three quantitative traits and one bi-allelic quantitative trait locus (QTL), and varied the number of traits associated with the QTL (explained variance 0.1%), minor allele frequency of the QTL, residual correlation between the traits, and the sign of the correlation induced by the QTL relative to the residual correlation. We compared the power of the methods using empirically fixed significance thresholds (α = 0.05). Our results showed that the multivariate methods implemented in PLINK, SNPTEST, MultiPhen and BIMBAM performed best for the majority of the tested scenarios, with a notable increase in power for scenarios with an opposite sign of genetic and residual correlation. All multivariate analyses resulted in a higher power than univariate analyses, even when only one of the traits was associated with the QTL. Hence, use of multivariate GWAS methods can be recommended, even when genetic correlations between traits are weak.

A computational study identifies HIV progression-related genes using mRMR and shortest path tracing

Since statistical relationships between HIV load and CD4+ T cell loss have been demonstrated to be weak, searching for host factors contributing to the pathogenesis of HIV infection becomes a key point for both understanding the disease pathology and developing treatments. We applied Maximum Relevance Minimum Redundancy (mRMR) algorithm to a set of microarray data generated from the CD4+ T cells of viremic non-progressors (VNPs) and rapid progressors (RPs) to identify host factors associated with the different responses to HIV infection. Using mRMR algorithm, 147 gene had been identified. Furthermore, we constructed a weighted molecular interaction network with the existing protein-protein interaction data from STRING database and identified 1331 genes on the shortest-paths among the genes identified with mRMR. Functional analysis shows that the functions relating to apoptosis play important roles during the pathogenesis of HIV infection. These results bring new insights of understanding HIV progression.

Genomic copy number variants: evidence for association with antibody response to anthrax vaccine adsorbed

Background

Anthrax and its etiologic agent remain a biological threat. Anthrax vaccine is highly effective, but vaccine-induced IgG antibody responses vary widely following required doses of vaccinations. Such variation can be related to genetic factors, especially genomic copy number variants (CNVs) that are known to be enriched among genes with immunologic function. We have tested this hypothesis in two study populations from a clinical trial of anthrax vaccination.

Methods

We performed CNV-based genome-wide association analyses separately on 794 European Americans and 200 African-Americans. Antibodies to protective antigen were measured at week 8 (early response) and week 30 (peak response) using an enzyme-linked immunosorbent assay. We used DNA microarray data (Affymetrix 6.0) and two CNV detection algorithms, hidden markov model (PennCNV) and circular binary segmentation (GeneSpring) to determine CNVs in all individuals. Multivariable regression analyses were used to identify CNV-specific associations after adjusting for relevant non-genetic covariates.

Results

Within the 22 autosomal chromosomes, 2,943 non-overlapping CNV regions were detected by both algorithms. Genomic insertions containing HLA-DRB5, DRB1 and DQA1/DRA genes in the major histocompatibility complex (MHC) region (chromosome 6p21.3) were moderately associated with elevated early antibody response (β = 0.14, p = 1.78×10⁻³) among European Americans, and the strongest association was observed between peak antibody response and a segmental insertion on chromosome 1, containing NBPF4, NBPF5, STXMP3, CLCC1, and GPSM2 genes (β = 1.66, p = 6.06×10⁻⁵). For African-Americans, segmental deletions spanning PRR20, PCDH17 and PCH68 genes on chromosome 13 were associated with elevated early antibody production (β = 0.18, p = 4.47×10⁻⁵). Population-specific findings aside, one genomic insertion on chromosome 17 (containing NSF, ARL17 and LRRC37A genes) was associated with elevated peak antibody response in both populations.

Conclusion

Multiple CNV regions, including the one consisting of MHC genes that is consistent with earlier research, can be important to humoral immune responses to anthrax vaccine adsorbed.