Molecular Profiling of Human Induced Pluripotent Stem Cell-Derived Hypothalamic Neurones Provides Developmental Insights into Genetic Loci for Body Weight Regulation

Recent data suggest that common genetic risks for metabolic disorders such as obesity may be human-specific and exert effects via the central nervous system. To overcome the limitation of human tissue access for study, we have generated induced human pluripotent stem cell (hiPSC)-derived neuronal cultures that recapture many features of hypothalamic neurones within the arcuate nucleus. In the present study, we have comprehensively characterised this model across development, benchmarked these neurones to in vivo events, and demonstrate a link between obesity risk variants and hypothalamic development. The dynamic transcriptome across neuronal maturation was examined using microarray and RNA sequencing methods at nine time points. K-means clustering of the longitudinal data was conducted to identify co-regulation and microRNA control of biological processes. The transcriptomes were compared with those of 103 samples from 13 brain regions reported in the Genotype-Tissue Expression database (GTEx) using principal components analysis. Genes with proximity to body mass index (BMI)-associated genetic variants were mapped to the developmentally expressed genesets, and enrichment significance was assessed with Fisher\u27s exact test. The human neuronal cultures have a transcriptional and physiological profile of neuropeptide Y/agouti-related peptide arcuate nucleus neurones. The neuronal transcriptomes were highly correlated with adult hypothalamus compared to any other brain region from the GTEx. Also, approximately 25% of the transcripts showed substantial changes in expression across neuronal development and potential co-regulation of biological processes that mirror neuronal development in vivo. These developmentally expressed genes were significantly enriched for genes in proximity to BMI-associated variants. We confirmed the utility of this in vitro human model for studying the development of key hypothalamic neurones involved in energy balance and show that genes at loci associated with body weight regulation may share a pattern of developmental regulation. These data support the need to investigate early development to elucidate the human-specific central nervous system pathophysiology underlying obesity susceptibility

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

Benchmarking Relatedness Inference Methods with Genome-Wide Data from Thousands of Relatives

Inferring relatedness from genomic data is an essential component of genetic association studies, population genetics, forensics, and genealogy. While numerous methods exist for inferring relatedness, thorough evaluation of these approaches in real data has been lacking. Here, we report an assessment of 12 state-of-the-art pairwise relatedness inference methods using a data set with 2485 individuals contained in several large pedigrees that span up to six generations. We find that all methods have high accuracy (92–99%) when detecting first- and second-degree relationships, but their accuracy dwindles to \u3c43% for seventh-degree relationships. However, most identical by descent (IBD) segment-based methods inferred seventh-degree relatives correct to within one relatedness degree for \u3e76% of relative pairs. Overall, the most accurate methods are Estimation of Recent Shared Ancestry (ERSA) and approaches that compute total IBD sharing using the output from GERMLINE and Refined IBD to infer relatedness. Combining information from the most accurate methods provides little accuracy improvement, indicating that novel approaches, such as new methods that leverage relatedness signals from multiple samples, are needed to achieve a sizeable jump in performance

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

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

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 \u3e2000 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

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

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 \u3c 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

The Development of FVIII Inhibitor in Hispanic American Patients with Hemophilia A Critically Impacts Coagulation Potential

Background: Hemophilia A (HA) is caused by deficiencies in plasma-FVIII and heterogeneous factor-VIII-gene mutations that impair intrinsic coagulation amplification. In severe hemophilia A patients (HAPs), FVIII infusions are begun at toddlerhood to prevent hemarthrosis induced crippling. However, approximately 30% of these patients develop FVIII inhibitors. Gain-of-function mutations in the common pathway of coagulation increases coagulation potential and decreases bleeding and FVIII-utilization in HAPs which should decrease FVIII-inhibitor-risk. We identified loss-of-function mutations in this pathway which decrease coagulation-potential as they increase FVIII-inhibitor risk in HAPs. Methods: We screened Mexican-American-pedigrees of the South-Texas-Family-Study (STFS) for protein-altering-variants. Subjects were genotyped using Illumina-exome-24-chip. Protein-altering-variants were analyzed for associations with FII:C, PT, and aPTT. Linear-mixed-model-analyses was performed to estimate trait-heritability and examine single-nucleotide-variations (SNVs) for gene association. Significant associations’ p-values fell below Bonferroni-adjusted significance level. Results: Heritability-estimates for FII:C, aPTT, and PT were highly-significant with p-values of 0.49, 0.49, and 0.54 (for all, pT in the FII-gene (F2)—which encodes 543R\u3eL and has a large effect-size on each trait (for all, pT have lower FII:C levels but correspondingly prolonged aPTT and PT times. Conclusion: We hypothesize that FII-543R\u3eL (Prothrombin-RGV) likely contributes to the high-incidence of FVIII-inhibitor-development in HA-patients of Mexican-ancestry, resulting in higher risk of developing anti-tFVIII-antibodies than patients without the variant. Patients with the RGV variant are likely to bleed more which can require surgery, further increasing the development of FVIII inhibitor development

Prosaposin is a regulator of progranulin levels and oligomerization

Progranulin (GRN) loss-of-function mutations leading to progranulin protein (PGRN) haploinsufficiency are prevalent genetic causes of frontotemporal dementia. Reports also indicated PGRN-mediated neuroprotection in models of Alzheimer’s and Parkinson’s disease; thus, increasing PGRN levels is a promising therapeutic for multiple disorders. To uncover novel PGRN regulators, we linked whole-genome sequence data from 920 individuals with plasma PGRN levels and identified the prosaposin (PSAP) locus as a new locus significantly associated with plasma PGRN levels. Here we show that both PSAP reduction and overexpression lead to significantly elevated extracellular PGRN levels. Intriguingly, PSAP knockdown increases PGRN monomers, whereas PSAP overexpression increases PGRN oligomers, partly through a protein–protein interaction. PSAP-induced changes in PGRN levels and oligomerization replicate in human-derived fibroblasts obtained from a GRN mutation carrier, further supporting PSAP as a potential PGRN-related therapeutic target. Future studies should focus on addressing the relevance and cellular mechanism by which PGRN oligomeric species provide neuroprotection