Rapid, Phase-free Detection of Long Identity-by-Descent Segments Enables Effective Relationship Classification.

Identity-by-descent (IBD) segments are a useful tool for applications ranging from demographic inference to relationship classification, but most detection methods rely on phasing information and therefore require substantial computation time. As genetic datasets grow, methods for inferring IBD segments that scale well will be critical. We developed IBIS, an IBD detector that locates long regions of allele sharing between unphased individuals, and benchmarked it with Refined IBD, GERMLINE, and TRUFFLE on 3,000 simulated individuals. Phasing these with Beagle 5 takes 4.3 CPU days, followed by either Refined IBD or GERMLINE segment detection in 2.9 or 1.1 h, respectively. By comparison, IBIS finishes in 6.8 min or 7.8 min with IBD2 functionality enabled: speedups of 805-946× including phasing time. TRUFFLE takes 2.6 h, corresponding to IBIS speedups of 20.2-23.3×. IBIS is also accurate, inferring ≥7 cM IBD segments at quality comparable to Refined IBD and GERMLINE. With these segments, IBIS classifies first through third degree relatives in real Mexican American samples at rates meeting or exceeding other methods tested and identifies fourth through sixth degree pairs at rates within 0.0%-2.0% of the top method. While allele frequency-based approaches that do not detect segments can infer relationship degrees faster than IBIS, the fastest are biased in admixed samples, with KING inferring 30.8% fewer fifth degree Mexican American relatives correctly compared with IBIS. Finally, we ran IBIS on chromosome 2 of the UK Biobank dataset and estimate its runtime on the autosomes to be 3.3 days parallelized across 128 cores.

Genetic determinants of metabolic biomarkers and their associations with cardiometabolic traits in Hispanic/Latino adolescents.

BACKGROUND: Metabolic regulation plays a significant role in energy homeostasis, and adolescence is a crucial life stage for the development of cardiometabolic disease (CMD). This study aims to investigate the genetic determinants of metabolic biomarkers-adiponectin, leptin, ghrelin, and orexin-and their associations with CMD risk factors. METHODS: We characterized the genetic determinants of the biomarkers among Hispanic/Latino adolescents of the Santiago Longitudinal Study (SLS) and identified the cumulative effects of genetic variants on adiponectin and leptin using biomarker polygenic risk scores (PRS). We further investigated the direct and indirect effect of the biomarker PRS on downstream body fat percent (BF%) and glycemic traits using structural equation modeling. RESULTS: We identified putatively novel genetic variants associated with the metabolic biomarkers. A substantial amount of biomarker variance was explained by SLS-specific PRS, and the prediction was improved by including the putatively novel loci. Fasting blood insulin and insulin resistance were associated with PRS for adiponectin, leptin, and ghrelin, and BF% was associated with PRS for adiponectin and leptin. We found evidence of substantial mediation of these associations by the biomarker levels. CONCLUSIONS: The genetic underpinnings of metabolic biomarkers can affect the early development of CMD, partly mediated by the biomarkers. IMPACT: This study characterized the genetic underpinnings of four metabolic hormones and investigated their potential influence on adiposity and insulin biology among Hispanic/Latino adolescents. Fasting blood insulin and insulin resistance were associated with polygenic risk score (PRS) for adiponectin, leptin, and ghrelin, with evidence of some degree of mediation by the biomarker levels. Body fat percent (BF%) was also associated with PRS for adiponectin and leptin. This provides important insight on biological mechanisms underlying early metabolic dysfunction and reveals candidates for prevention efforts. Our findings also highlight the importance of ancestrally diverse populations to facilitate valid studies of the genetic architecture of metabolic biomarker levels.

Crossover interference and sex-specific genetic maps shape identical by descent sharing in close relatives.

Simulations of close relatives and identical by descent (IBD) segments are common in genetic studies, yet most past efforts have utilized sex averaged genetic maps and ignored crossover interference, thus omitting features known to affect the breakpoints of IBD segments. We developed Ped-sim, a method for simulating relatives that can utilize either sex-specific or sex averaged genetic maps and also either a model of crossover interference or the traditional Poisson model for inter-crossover distances. To characterize the impact of previously ignored mechanisms, we simulated data for all four combinations of these factors. We found that modeling crossover interference decreases the standard deviation of pairwise IBD proportions by 10.4% on average in full siblings through second cousins. By contrast, sex-specific maps increase this standard deviation by 4.2% on average, and also impact the number of segments relatives share. Most notably, using sex-specific maps, the number of segments half-siblings share is bimodal; and when combined with interference modeling, the probability that sixth cousins have non-zero IBD sharing ranges from 9.0 to 13.1%, depending on the sexes of the individuals through which they are related. We present new analytical results for the distributions of IBD segments under these models and show they match results from simulations. Finally, we compared IBD sharing rates between simulated and real relatives and find that the combination of sex-specific maps and interference modeling most accurately captures IBD rates in real data. Ped-sim is open source and available from https://github.com/williamslab/ped-sim.

Inferring Identical-by-Descent Sharing of Sample Ancestors Promotes High-Resolution Relative Detection.

As genetic datasets increase in size, the fraction of samples with one or more close relatives grows rapidly, resulting in sets of mutually related individuals. We present DRUID-deep relatedness utilizing identity by descent-a method that works by inferring the identical-by-descent (IBD) sharing profile of an ungenotyped ancestor of a set of close relatives. Using this IBD profile, DRUID infers relatedness between unobserved ancestors and more distant relatives, thereby combining information from multiple samples to remove one or more generations between the deep relationships to be identified. DRUID constructs sets of close relatives by detecting full siblings and also uses an approach to identify the aunts/uncles of two or more siblings, recovering 92.2% of real aunts/uncles with zero false positives. In real and simulated data, DRUID correctly infers up to 10.5% more relatives than PADRE when using data from two sets of distantly related siblings, and 10.7%-31.3% more relatives given two sets of siblings and their aunts/uncles. DRUID frequently infers relationships either correctly or within one degree of the truth, with PADRE classifying 43.3%-58.3% of tenth degree relatives in this way compared to 79.6%-96.7% using DRUID.

APOC3 genetic variation, serum triglycerides, and risk of coronary artery disease in Asian Indians, Europeans, and other ethnic groups.

BACKGROUND: Hypertriglyceridemia has emerged as a critical coronary artery disease (CAD) risk factor. Rare loss-of-function (LoF) variants in apolipoprotein C-III have been reported to reduce triglycerides (TG) and are cardioprotective in American Indians and Europeans. However, there is a lack of data in other Europeans and non-Europeans. Also, whether genetically increased plasma TG due to ApoC-III is causally associated with increased CAD risk is still unclear and inconsistent. The objectives of this study were to verify the cardioprotective role of earlier reported six LoF variants of APOC3 in South Asians and other multi-ethnic cohorts and to evaluate the causal association of TG raising common variants for increasing CAD risk. METHODS: We performed gene-centric and Mendelian randomization analyses and evaluated the role of genetic variation encompassing APOC3 for affecting circulating TG and the risk for developing CAD. RESULTS: One rare LoF variant (rs138326449) with a 37% reduction in TG was associated with lowered risk for CAD in Europeans (p = 0.007), but we could not confirm this association in Asian Indians (p = 0.641). Our data could not validate the cardioprotective role of other five LoF variants analysed. A common variant rs5128 in the APOC3 was strongly associated with elevated TG levels showing a p-value 2.8 × 10- 424. Measures of plasma ApoC-III in a small subset of Sikhs revealed a 37% increase in ApoC-III concentrations among homozygous mutant carriers than the wild-type carriers of rs5128. A genetically instrumented per 1SD increment of plasma TG level of 15 mg/dL would cause a mild increase (3%) in the risk for CAD (p = 0.042). CONCLUSIONS: Our results highlight the challenges of inclusion of rare variant information in clinical risk assessment and the generalizability of implementation of ApoC-III inhibition for treating atherosclerotic disease. More studies would be needed to confirm whether genetically raised TG and ApoC-III concentrations would increase CAD risk.

Role of miRNA-mRNA Interaction in Neural Stem Cell Differentiation of Induced Pluripotent Stem Cells.

miRNA regulates the expression of protein coding genes and plays a regulatory role in human development and disease. The human iPSCs and their differentiated progenies provide a unique opportunity to identify these miRNA-mediated regulatory mechanisms. To identify miRNA-mRNA regulatory interactions in human nervous system development, well characterized NSCs were differentiated from six validated iPSC lines and analyzed for differentially expressed (DE) miRNome and transcriptome by RNA sequencing. Following the criteria, moderated t statistics, FDR-corrected p-value ≤ 0.05 and fold change-absolute (FC-abs) ≥2.0, 51 miRNAs and 4033 mRNAs were found to be significantly DE between iPSCs and NSCs. The miRNA target prediction analysis identified 513 interactions between 30 miRNA families (mapped to 51 DE miRNAs) and 456 DE mRNAs that were paradoxically oppositely expressed. These 513 interactions were highly enriched in nervous system development functions (154 mRNAs; FDR-adjusted p-value range: 8.06 × 10-15-1.44 × 10-4). Furthermore, we have shown that the upregulated miR-10a-5p, miR-30c-5p, miR23-3p, miR130a-3p and miR-17-5p miRNA families were predicted to down-regulate several genes associated with the differentiation of neurons, neurite outgrowth and synapse formation, suggesting their role in promoting the self-renewal of undifferentiated NSCs. This study also provides a comprehensive characterization of iPSC-generated NSCs as dorsal neuroepithelium, important for their potential use in in vitro modeling of human brain development and disease.

Genetic and environmental (physical fitness and sedentary activity) interaction effects on cardiometabolic risk factors in Mexican American children and adolescents.

Knowledge on genetic and environmental (G × E) interaction effects on cardiometabolic risk factors (CMRFs) in children is limited.  The purpose of this study was to examine the impact of G × E interaction effects on CMRFs in Mexican American (MA) children (n = 617, ages 6-17 years). The environments examined were sedentary activity (SA), assessed by recalls from "yesterday" (SAy) and "usually" (SAu) and physical fitness (PF) assessed by Harvard PF scores (HPFS). CMRF data included body mass index (BMI), waist circumference (WC), fat mass (FM), fasting insulin (FI), homeostasis model of assessment-insulin resistance (HOMA-IR), high-density lipoprotein cholesterol (HDL-C), triglycerides (TG), systolic (SBP) and diastolic (DBP) blood pressure, and number of metabolic syndrome components (MSC). We examined potential G × E interaction in the phenotypic expression of CMRFs using variance component models and likelihood-based statistical inference. Significant G × SA interactions were identified for six CMRFs: BMI, WC, FI, HOMA-IR, MSC, and HDL, and significant G × HPFS interactions were observed for four CMRFs: BMI, WC, FM, and HOMA-IR. However, after correcting for multiple hypothesis testing, only WC × SAy, FM × SAy, and FI × SAu interactions became marginally significant. After correcting for multiple testing, most of CMRFs exhibited significant G × E interactions (Reduced G × E model vs. Constrained model). These findings provide evidence that genetic factors interact with SA and PF to influence variation in CMRFs, and underscore the need for better understanding of these relationships to develop strategies and interventions to effectively reduce or prevent cardiometabolic risk in children.

GWAS and transcriptional analysis prioritize ITPR1 and CNTN4 for a serum uric acid 3p26 QTL in Mexican Americans.

BACKGROUND: The variation in serum uric acid concentrations is under significant genetic influence. Elevated SUA concentrations have been linked to increased risk for gout, kidney stones, chronic kidney disease, and cardiovascular disease whereas reduced serum uric acid concentrations have been linked to multiple sclerosis, Parkinson's disease and Alzheimer's disease. Previously, we identified a novel locus on chromosome 3p26 affecting serum uric acid concentrations in Mexican Americans from San Antonio Family Heart Study. As a follow up, we examined genome-wide single nucleotide polymorphism data in an extended cohort of 1281 Mexican Americans from multigenerational families of the San Antonio Family Heart Study and the San Antonio Family Diabetes/Gallbladder Study. We used a linear regression-based joint linkage/association test under an additive model of allelic effect, while accounting for non-independence among family members via a kinship variance component. RESULTS: Univariate genetic analysis indicated serum uric acid concentrations to be significant heritable (h (2) = 0.50 ± 0.05, p < 4 × 10(-35)), and linkage analysis of serum uric acid concentrations confirmed our previous finding of a novel locus on 3p26 (LOD = 4.9, p < 1 × 10(-5)) in the extended sample. Additionally, we observed strong association of serum uric acid concentrations with variants in following candidate genes in the 3p26 region; inositol 1,4,5-trisphosphate receptor, type 1 (ITPR1), contactin 4 (CNTN4), decapping mRNA 1A (DCP1A); transglutaminase 4 (TGM4) and rho guanine nucleotide exchange factor (GEF) 26 (ARHGEF26) [p < 3 × 10(-7); minor allele frequencies ranged between 0.003 and 0.42] and evidence of cis-regulation for ITPR1 transcripts. CONCLUSION: Our results confirm the importance of the chromosome 3p26 locus and genetic variants in this region in the regulation of serum uric acid concentrations.

The Arg59Trp variant in ANGPTL8 (betatrophin) is associated with total and HDL-cholesterol in American Indians and Mexican Americans and differentially affects cleavage of ANGPTL3.

We previously identified a locus linked to total cholesterol (TC) concentration in Pima Indians on chromosome 19p. To characterize this locus, we genotyped >2000 SNPs in 1838 Pimas and assessed association with log(TC). We observed evidence for association with log(TC) with rs2278426 (3.5% decrease/copy of the T allele; P=5.045×10(-6)) in the ANGPTL8 (angiopoietin-like 8) gene. We replicated this association in 2413 participants of the San Antonio Mexican American Family Study (SAMAFS: 2.0% decrease per copy of the T allele; P=0.005842). In a meta-analysis of the combined data, we found the strongest estimated effect with rs2278426 (P=2.563×10(-7)). The variant T allele at rs2278426 predicts an Arg59Trp substitution and has previously been associated with LDL-C and HDL-C. In Pimas and SAMAFS participants, the T allele of rs2278426 was associated with reduced HDL-C levels (P=0.000741 and 0.00002, respectively), and the combined estimated effect for the two cohorts was -3.8% (P=8.526×10(-8)). ANGPTL8 transcript and protein levels increased in response to both glucose and insulin. The variant allele was associated with increased levels of cleaved ANGPTL3. We conclude that individuals with the variant allele may have lower TC and HDL-C levels due to increased activation of ANGPTL3 by ANGPTL8.

Linkage of type 2 diabetes on chromosome 9p24 in Mexican Americans: additional evidence from the Veterans Administration Genetic Epidemiology Study (VAGES).

OBJECTIVE: Type 2 diabetes (T2DM) is a complex metabolic disease and is more prevalent in certain ethnic groups such as the Mexican Americans. The goal of our study was to perform a genome-wide linkage (GWL) analysis to localize T2DM susceptibility loci in Mexican Americans. METHODS: We used the phenotypic and genotypic data from 1,122 Mexican-American individuals (307 families) who participated in the Veterans Administration Genetic Epidemiology Study (VAGES). GWL analysis was performed using the variance components approach. Data from 2 additional Mexican-American family studies, the San Antonio Family Heart Study (SAFHS) and the San Antonio Family Diabetes/Gallbladder Study (SAFDGS), were combined with the VAGES data to test for improved linkage evidence. RESULTS: After adjusting for covariate effects, T2DM was found to be under significant genetic influences (h2 = 0.62, p = 2.7 × 10(-6)). The strongest evidence for linkage of T2DM occurred between markers D9S1871 and D9S2169 on chromosome 9p24.2-p24.1 (LOD = 1.8). Given that we previously reported suggestive evidence for linkage of T2DM at this region also in SAFDGS, we found the significant and increased linkage evidence (LOD = 4.3, empirical p = 1.0 × 10(-5), genome-wide p = 1.6 × 10(-3)) for T2DM at the same chromosomal region, when we performed a GWL analysis of the VAGES data combined with the SAFHS and SAFDGS data. CONCLUSION: Significant T2DM linkage evidence was found on chromosome 9p24 in Mexican Americans. Importantly, the chromosomal region of interest in this study overlaps with several recent genome-wide association studies involving T2DM-related traits. Given its overlap with such findings and our own initial T2DM association findings in the 9p24 chromosomal region, high throughput sequencing of the linked chromosomal region could identify the potential causal T2DM genes.

Association of a low-frequency variant in HNF1A with type 2 diabetes in a Latino population.

IMPORTANCE: Latino populations have one of the highest prevalences of type 2 diabetes worldwide. OBJECTIVES: To investigate the association between rare protein-coding genetic variants and prevalence of type 2 diabetes in a large Latino population and to explore potential molecular and physiological mechanisms for the observed relationships. DESIGN, SETTING, AND PARTICIPANTS: Whole-exome sequencing was performed on DNA samples from 3756 Mexican and US Latino individuals (1794 with type 2 diabetes and 1962 without diabetes) recruited from 1993 to 2013. One variant was further tested for allele frequency and association with type 2 diabetes in large multiethnic data sets of 14,276 participants and characterized in experimental assays. MAIN OUTCOME AND MEASURES: Prevalence of type 2 diabetes. Secondary outcomes included age of onset, body mass index, and effect on protein function. RESULTS: A single rare missense variant (c.1522G>A [p.E508K]) was associated with type 2 diabetes prevalence (odds ratio [OR], 5.48; 95% CI, 2.83-10.61; P = 4.4 × 10(-7)) in hepatocyte nuclear factor 1-α (HNF1A), the gene responsible for maturity onset diabetes of the young type 3 (MODY3). This variant was observed in 0.36% of participants without type 2 diabetes and 2.1% of participants with it. In multiethnic replication data sets, the p.E508K variant was seen only in Latino patients (n = 1443 with type 2 diabetes and 1673 without it) and was associated with type 2 diabetes (OR, 4.16; 95% CI, 1.75-9.92; P = .0013). In experimental assays, HNF-1A protein encoding the p.E508K mutant demonstrated reduced transactivation activity of its target promoter compared with a wild-type protein. In our data, carriers and noncarriers of the p.E508K mutation with type 2 diabetes had no significant differences in compared clinical characteristics, including age at onset. The mean (SD) age for carriers was 45.3 years (11.2) vs 47.5 years (11.5) for noncarriers (P = .49) and the mean (SD) BMI for carriers was 28.2 (5.5) vs 29.3 (5.3) for noncarriers (P = .19). CONCLUSIONS AND RELEVANCE: Using whole-exome sequencing, we identified a single low-frequency variant in the MODY3-causing gene HNF1A that is associated with type 2 diabetes in Latino populations and may affect protein function. This finding may have implications for screening and therapeutic modification in this population, but additional studies are required.

Genetic epidemiology of cardiometabolic risk factors and their clustering patterns in Mexican American children and adolescents: the SAFARI Study.

Pediatric metabolic syndrome (MS) and its cardiometabolic components (MSCs) have become increasingly prevalent, yet little is known about the genetics underlying MS risk in children. We examined the prevalence and genetics of MS-related traits among 670 non-diabetic Mexican American (MA) children and adolescents, aged 6-17 years (49 % female), who were participants in the San Antonio Family Assessment of Metabolic Risk Indicators in Youth study. These children are offspring or biological relatives of adult participants from three well-established Mexican American family studies in San Antonio, TX, at increased risk of type 2 diabetes. MS was defined as ≥3 abnormalities among 6 MSC measures: waist circumference, systolic and/or diastolic blood pressure, fasting insulin, triglycerides, HDL-cholesterol, and fasting and/or 2-h OGTT glucose. Genetic analyses of MS, number of MSCs (MSC-N), MS factors, and bivariate MS traits were performed. Overweight/obesity (53 %), pre-diabetes (13 %), acanthosis nigricans (33 %), and MS (19 %) were strikingly prevalent, as were MS components, including abdominal adiposity (32 %) and low HDL-cholesterol (32 %). Factor analysis of MS traits yielded three constructs: adipo-insulin-lipid, blood pressure, and glucose factors, and their factor scores were highly heritable. MS itself exhibited 68 % heritability. MSC-N showed strong positive genetic correlations with obesity, insulin resistance, inflammation, and acanthosis nigricans, and negative genetic correlation with physical fitness. MS trait pairs exhibited strong genetic and/or environmental correlations. These findings highlight the complex genetic architecture of MS/MSCs in MA children, and underscore the need for early screening and intervention to prevent chronic sequelae in this vulnerable pediatric population.

Functional characterization of the disease-associated CCL2 rs1024611G-rs13900T haplotype: The role of the RNA-binding protein HuR.

CC-chemokine ligand 2 (CCL2) is involved in the pathogenesis of several diseases associated with monocyte/macrophage recruitment, such as HIV-associated neurocognitive disorder (HAND), tuberculosis, and atherosclerosis. The rs1024611 (alleles:A>G; G is the risk allele) polymorphism in the CCL2 cis-regulatory region is associated with increased CCL2 expression in vitro and ex vivo, leukocyte mobilization in vivo, and deleterious disease outcomes. However, the molecular basis for the rs1024611-associated differential CCL2 expression remains poorly characterized. It is conceivable that genetic variant(s) in linkage disequilibrium (LD) with rs1024611 could mediate such effects. Previously, we used rs13900 (alleles:_C>T) in the CCL2 3' untranslated region (3' UTR) that is in perfect LD with rs1024611 to demonstrate allelic expression imbalance (AEI) of CCL2 in heterozygous individuals. Here we tested the hypothesis that the rs13900 could modulate CCL2 expression by altering mRNA turnover and/or translatability. The rs13900 T allele conferred greater stability to the CCL2 transcript when compared to the rs13900 C allele. The rs13900 T allele also had increased binding to Human Antigen R (HuR), an RNA-binding protein, in vitro and ex vivo. The rs13900 alleles imparted differential activity to reporter vectors and influenced the translatability of the reporter transcript. We further demonstrated a role for HuR in mediating allele-specific effects on CCL2 expression in overexpression and silencing studies. The presence of the rs1024611G-rs13900T conferred a distinct transcriptomic signature related to inflammation and immunity. Our studies suggest that the differential interactions of HuR with rs13900 could modulate CCL2 expression and explain the interindividual differences in CCL2-mediated disease susceptibility.

Wide disparity in genetic admixture among Mexican Americans from San Antonio, TX.

We studied 706 participants of the San Antonio Family Diabetes Study (SAFDS) and 586 male samples from the San Antonio Center for Biomarkers of Risk of Prostate Cancer (SABOR) and used 64 ancestry informative markers to compare admixture proportions between both groups. Existence of population substructure was demonstrated by the excess association of unlinked markers. In the SAFDS sample, ancestral proportions were estimated at 50.2 ± 0.6% European, 46.4 ± 0.6% Native American, and 3.1 ± 0.2% West African. For the SABOR sample, the proportions were 58.9 ± 0.7%, 38.2 ± 0.7%, and 2.9 ± 0.2%, respectively. Additionally, in the SAFDS subjects a highly significant negative correlation was found between individual Native American ancestry and skin reflectance (R(2) = 0.07, P= 0.00006). The correlation was stronger in males than in females but clearly showed that ancestry only accounts for a small percentage of the variation in skin color and, conversely, that skin reflectance is not a robust surrogate for genetic admixture. Furthermore, a substantial difference in substructure is present in the two cohorts of Mexican American subjects from the San Antonio area in Texas, which emphasizes that genetic admixture estimates should be accounted for in association studies, even for geographically related subjects.

The Ser(326)Cys polymorphism of 8-oxoguanine glycosylase 1 (OGG1) is associated with type 2 diabetes in Mexican Americans.

OBJECTIVE: Human 8-oxoguanine glycosylase 1 (OGG1) excises oxidatively damaged promutagenic base 8-oxoguanine, a lesion previously observed in a rat model of type 2 diabetes (T2DM). The objective of the present study is to determine whether genetic variation in OGG1 is associated with type 2 diabetes (T2DM) in a Mexican American cohort. METHODS: Ten SNPs including two tagging SNPs (rs1052133, rs2072668) across the OGG1 gene region were selected from the Hapmap database and genotyped in the entire cohort (n = 670; 29% diabetes; 39 families) by TaqMan assay. Association analyses between the SNPs and T2DM were performed using the measured genotype approach as implemented in the program SOLAR. RESULTS: Of the ten SNPs genotyped, only five were polymorphic. The minor allele frequencies of these 5 SNPs ranged from 1-38%. Of the SNPs examined for association, the Ser(326)Cys (rs1052133) exhibited significant association with T2DM (p = 0.016) after accounting for age and sex effects. Another intronic variant (rs2072668), which was in strong linkage disequilibrium (r(2) = 0.96) with Ser(326)Cys also exhibited significant association with T2DM (p = 0.031). CONCLUSIONS: These results suggest for the first time that the variants in OGG1 could influence diabetes risk in these Mexican American families and support a role for alterations of OGG1 in the pathogenesis of T2DM.

A functional genomics pipeline identifies pleiotropy and cross-tissue effects within obesity-associated GWAS loci.

Genome-wide association studies (GWAS) have identified many disease-associated variants, yet mechanisms underlying these associations remain unclear. To understand obesity-associated variants, we generate gene regulatory annotations in adipocytes and hypothalamic neurons across cellular differentiation stages. We then test variants in 97 obesity-associated loci using a massively parallel reporter assay and identify putatively causal variants that display cell type specific or cross-tissue enhancer-modulating properties. Integrating these variants with gene regulatory information suggests genes that underlie obesity GWAS associations. We also investigate a complex genomic interval on 16p11.2 where two independent loci exhibit megabase-range, cross-locus chromatin interactions. We demonstrate that variants within these two loci regulate a shared gene set. Together, our data support a model where GWAS loci contain variants that alter enhancer activity across tissues, potentially with temporally restricted effects, to impact the expression of multiple genes. This complex model has broad implications for ongoing efforts to understand GWAS.

Association of HIV-1 Infection and Antiretroviral Therapy With Type 2 Diabetes in the Hispanic Population of the Rio Grande Valley, Texas, USA.

The Rio Grande Valley (RGV) in South Texas has one of the highest prevalence of obesity and type 2 diabetes (T2D) in the United States (US). We report for the first time the T2D prevalence in persons with HIV (PWH) in the RGV and the interrelationship between T2D, cardiometabolic risk factors, HIV-related indices, and antiretroviral therapies (ART). The PWH in this study received medical care at Valley AIDS Council (VAC) clinic sites located in Harlingen and McAllen, Texas. Henceforth, this cohort will be referred to as Valley AIDS Council Cohort (VACC). Cross-sectional analyses were conducted using retrospective data obtained from 1,827 registries. It included demographic and anthropometric variables, cardiometabolic traits, and HIV-related virological and immunological indices. For descriptive statistics, we used mean values of the quantitative variables from unbalanced visits across 20 months. Robust regression methods were used to determine the associations. For comparisons, we used cardiometabolic trait data obtained from HIV-uninfected San Antonio Mexican American Family Studies (SAMAFS; N = 2,498), and the Mexican American population in the National Health and Nutrition Examination Survey (HHANES; N = 5,989). The prevalence of T2D in VACC was 51% compared to 27% in SAMAFS and 19% in HHANES, respectively. The PWH with T2D in VACC were younger (4.7 years) and had lower BMI (BMI 2.43 units less) when compared to SAMAFS individuals. In contrast, VACC individuals had increased blood pressure and dyslipidemia. The increased T2D prevalence in VACC was independent of BMI. Within the VACC, ART was associated with viral load and CD4+ T cell counts but not with metabolic dysfunction. Notably, we found that individuals with any INSTI combination had higher T2D risk: OR 2.08 (95%CI 1.67, 2.6; p < 0.001). In summary, our results suggest that VACC individuals may develop T2D at younger ages independent of obesity. The high burden of T2D in these individuals necessitates rigorously designed longitudinal studies to draw potential causal inferences and develop better treatment regimens.

Serum carotenoids and Pediatric Metabolic Index predict insulin sensitivity in Mexican American children.

High concentrations of carotenoids are protective against cardiometabolic risk traits (CMTs) in adults and children. We recently showed in non-diabetic Mexican American (MA) children that serum α-carotene and ß-carotene are inversely correlated with obesity measures and triglycerides and positively with HDL cholesterol and that they were under strong genetic influences. Additionally, we previously described a Pediatric Metabolic Index (PMI) that helps in the identification of children who are at risk for cardiometabolic diseases. Here, we quantified serum lycopene and ß-cryptoxanthin concentrations in approximately 580 children from MA families using an ultraperformance liquid chromatography-photodiode array and determined their heritabilities and correlations with CMTs. Using response surface methodology (RSM), we determined two-way interactions of carotenoids and PMI on Matsuda insulin sensitivity index (ISI). The concentrations of lycopene and ß-cryptoxanthin were highly heritable [h2 = 0.98, P = 7 × 10-18 and h2 = 0.58, P = 1 × 10-7]. We found significant (P ≤ 0.05) negative phenotypic correlations between ß-cryptoxanthin and five CMTs: body mass index (- 0.22), waist circumference (- 0.25), triglycerides (- 0.18), fat mass (- 0.23), fasting glucose (- 0.09), and positive correlations with HDL cholesterol (0.29). In contrast, lycopene only showed a significant negative correlation with fasting glucose (- 0.08) and a positive correlation with HDL cholesterol (0.18). Importantly, we found that common genetic influences significantly contributed to the observed phenotypic correlations. RSM showed that increased serum concentrations of α- and ß-carotenoids rather than that of ß-cryptoxanthin or lycopene had maximal effects on ISI. In summary, our findings suggest that the serum carotenoids are under strong additive genetic influences and may have differential effects on susceptibility to CMTs in children.

Genotype phasing in pedigrees using whole-genome sequence data.

Phasing is the process of inferring haplotypes from genotype data. Efficient algorithms and associated software for accurate phasing in pedigrees are needed, especially for populations lacking reference panels of sequenced individuals. We present a novel method for phasing genotypes from whole-genome sequence data in pedigrees, called PULSAR (Phasing Using Lineage Specific Alleles/Rare variants). The method is based on the property that alleles specific to a single founding chromosome within a pedigree are highly informative for identifying haplotypes that are shared identical by descent. Simulation studies are used to assess the performance of PULSAR with various pedigree sizes and structures, and the effect of genotyping errors and the presence of nonsequenced individuals is investigated. In pedigrees with complete sequencing and realistic genotyping error rates, PULSAR correctly phases >99.9% of heterozygous genotypes, excluding sites at which all individuals are heterozygous, and does so with a switch error rate frequently below 10-4. PULSAR is highly accurate, capable of genotype error correction and imputation, and computationally competitive with alternative phasing software applicable to pedigrees. Our method has the significant advantage of not requiring reference panels that are essential for other population-based phasing algorithms. A software implementation of PULSAR is freely available.

Acanthosis nigricans as a composite marker of cardiometabolic risk and its complex association with obesity and insulin resistance in Mexican American children.

AIM: Acanthosis nigricans (AN) is a strong correlate of obesity and is considered a marker of insulin resistance (IR). AN is associated with various other cardiometabolic risk factors (CMRFs). However, the direct causal relationship of IR with AN in obesity has been debated. Therefore, we aimed to examine the complex causal relationships among the troika of AN, obesity, and IR in Mexican Americans (MAs). METHODS: We used data from 670 non-diabetic MA children, aged 6-17 years (49% girls). AN (prevalence 33%) severity scores (range 0-5) were used as a quasi-quantitative trait (AN-q) for analysis. We used the program SOLAR for determining phenotypic, genetic, and environmental correlations between AN-q and CMRFs (e.g., BMI, HOMA-IR, lipids, blood pressure, hs-C-reactive protein (CRP), and Harvard physical fitness score (PFS)). The genetic and environmental correlations were subsequently used in mediation analysis (AMOS program). Model comparisons were made using goodness-of-fit indexes. RESULTS: Heritability of AN-q was 0.75 (p<0.0001). It was positively/significantly (p<0.05) correlated with traits such as BMI, HOMA-IR, and CRP, and negatively with HDL-C and PFS. Of the models tested, indirect mediation analysis of BMI→HOMA-IR→AN-q yielded lower goodness-of-fit than a partial mediation model where BMI explained the relationship with both HOMA-IR and AN-q simultaneously. Using complex models, BMI was associated with AN-q and IR mediating most of the CMRFs; but no relationship between IR and AN-q. CONCLUSION: Our study suggests that obesity explains the association of IR with AN, but no causal relationship between IR and AN in Mexican American children.