Exome Sequencing Identifies Genetic Variants Associated with Circulating Lipid Levels in Mexican Americans: The Insulin Resistance Atherosclerosis Family Study (IRASFS)

Association of Apolipoprotein A5 Gene Variants with Hyperlipidemic Acute Pancreatitis in Southeastern China

Background: Apolipoprotein A5 (APOA5) is involved in serum triglyceride (TG) regulation. Several studies have reported that the rs651821 locus in the APOA5 gene is associated with serum TG levels in the Chinese population. However, no research has been performed regarding the association between the variants of rs651821 and the risk of hyperlipidemic acute pancreatitis (HLAP). Methods: A case-control study was conducted and is reported following the STROBE guidelines. We enrolled a total of 88 participants in this study (60 HLAP patients and 28 controls). APOA5 was genotyped using PCR and Sanger sequencing. Logistic regression models were conducted to calculate odds ratios and a 95% confidence interval. Results: The genotype distribution of the rs651821 alleles in both groups follow the Hardy-Weinberg distribution. The frequency of the "C" allele in rs651821 was increased in HLAP patients compared to controls. In the recessive model, subjects with the "CC" genotype had an 8.217-fold higher risk for HLAP (OR = 8.217, 95% CI: 1.023-66.01, p = 0.046) than subjects with the "TC+TT" genotypes. After adjusting for sex, the association remained significant (OR = 9.898, 95% CI: 1.176-83.344, p = 0.035). Additionally, the "CC" genotype was related to an increased TG/apolipoprotein B (APOB) ratio and fasting plasma glucose (FPG) levels. Conclusions: Our findings suggest that the C allele of rs651821 in APOA5 increases the risk of HLAP in persons from Southeastern China.

Exome Sequencing Data Analysis and a Case-Control Study in Mexican Population Reveals Lipid Trait Associations of New and Known Genetic Variants in Dyslipidemia-Associated Loci

Background: Plasma lipid levels are a major risk factor for cardiovascular diseases. Although international efforts have identified a group of loci associated with the risk of dyslipidemia, Latin American populations have been underrepresented in these studies.

Objective: To know the genetic variation occurring in lipid-related loci in the Mexican population and its association with dyslipidemia.

Methods: We searched for single-nucleotide variants in 177 lipid candidate genes using previously published exome sequencing data from 2838 Mexican individuals belonging to three different cohorts. With the extracted variants, we performed a case-control study. Logistic regression and quantitative trait analyses were implemented in PLINK software. We used an LD pruning using a 50-kb sliding window size, a 5-kb window step size and a r2 threshold of 0.1.

Results: Among the 34251 biallelic variants identified in our sample population, 33% showed low frequency. For case-control study, we selected 2521 variants based on a minor allele frequency ≥1% in all datasets. We found 19 variants in 9 genes significantly associated with at least one lipid trait, with the most significant associations found in the APOA1/C3/A4/A5-ZPR1-BUD13 gene cluster on chromosome 11. Notably, all 11 variants associated with hypertriglyceridemia were within this cluster; whereas variants associated with hypercholesterolemia were located at chromosome 2 and 19, and for low high density lipoprotein cholesterol were in chromosomes 9, 11, and 19. No significant associated variants were found for low density lipoprotein. We found several novel variants associated with different lipemic traits: rs3825041 in BUD13 with hypertriglyceridemia, rs7252453 in CILP2 with decreased risk to hypercholesterolemia and rs11076176 in CETP with increased risk to low high density lipoprotein cholesterol.

Conclusions: We identified novel variants in lipid-regulation candidate genes in the Mexican population, an underrepresented population in genomic studies, demonstrating the necessity of more genomic studies on multi-ethnic populations to gain a deeper understanding of the genetic structure of the lipemic traits.

Metabolomic architecture of obesity implicates metabolonic lactone sulfate in cardiometabolic disease

OBJECTIVE

Identify and characterize circulating metabolite profiles associated with adiposity to inform precision medicine.

METHODS

Untargeted plasma metabolomic profiles in the Insulin Resistance Atherosclerosis Family Study (IRASFS) Mexican American cohort (n = 1108) were analyzed for association with anthropometric (body mass index, BMI; waist circumference, WC; waist-to-hip ratio, WHR) and computed tomography measures (visceral adipose tissue, VAT; subcutaneous adipose tissue, SAT; visceral-to-subcutaneous ratio, VSR) of adiposity. Genetic data, inclusive of genome-wide array-based genotyping, whole exome sequencing (WES) and whole genome sequencing (WGS), were evaluated to identify the genetic contributors. Phenotypic and genetic association signals were replicated across ancestries. Transcriptomic data were analyzed to explore the relationship between genetic and metabolomic data.

RESULTS

A partially characterized metabolite, tentatively named metabolonic lactone sulfate (X-12063), was consistently associated with BMI, WC, WHR, VAT, and SAT in IRASFS Mexican Americans (P_MA <2.02 × 10^-27). Trait associations were replicated in IRASFS African Americans (P_AA < 1.12 × 10^-07). Expanded analyses revealed associations with multiple phenotypic measures of cardiometabolic health, e.g. insulin sensitivity (S_I), triglycerides (TG), diastolic blood pressure (DBP) and plasminogen activator inhibitor-1 (PAI-1) in both ancestries. Metabolonic lactone sulfate levels were heritable (h² > 0.47), and a significant genetic signal at the ZSCAN25/CYP3A5 locus (P_MA = 9.00 × 10^-41, P_AA = 2.31 × 10^-10) was observed, highlighting a putative functional variant (rs776746, CYP3A5∗3). Transcriptomic analysis in the African American Genetics of Metabolism and Expression (AAGMEx) cohort supported the association of CYP3A5 with metabolonic lactone sulfate levels (P_FDR = 6.64 × 10^-07).

CONCLUSIONS

Variant rs776746 is associated with a decrease in the transcript levels of CYP3A5, which in turn is associated with increased metabolonic lactone sulfate levels and poor cardiometabolic health.

Genome-Wide Association Study Identifies a Functional SIDT2 Variant Associated With HDL-C (High-Density Lipoprotein Cholesterol) Levels and Premature Coronary Artery Disease

Objective

Low HDL-C (high-density lipoprotein cholesterol) is the most frequent dyslipidemia in Mexicans, but few studies have examined the underlying genetic basis. Our purpose was to identify genetic variants associated with HDL-C levels and cardiovascular risk in the Mexican population.

Approach and Results

A genome-wide association studies for HDL-C levels in 2335 Mexicans, identified four loci associated with genome-wide significance: CETP, ABCA1, LIPC, and SIDT2. The SIDT2 missense Val636Ile variant was associated with HDL-C levels and was replicated in 3 independent cohorts (P=5.9×10−18 in the conjoint analysis). The SIDT2/Val636Ile variant is more frequent in Native American and derived populations than in other ethnic groups. This variant was also associated with increased ApoA1 and glycerophospholipid serum levels, decreased LDL-C (low-density lipoprotein cholesterol) and ApoB levels, and a lower risk of premature CAD. Because SIDT2 was previously identified as a protein involved in sterol transport, we tested whether the SIDT2/Ile636 protein affected this function using an in vitro site-directed mutagenesis approach. The SIDT2/Ile636 protein showed increased uptake of the cholesterol analog dehydroergosterol, suggesting this variant affects function. Finally, liver transcriptome data from humans and the Hybrid Mouse Diversity Panel are consistent with the involvement of SIDT2 in lipid and lipoprotein metabolism.

Conclusions

This is the first genome-wide association study for HDL-C levels seeking associations with coronary artery disease in the Mexican population. Our findings provide new insight into the genetic architecture of HDL-C and highlight SIDT2 as a new player in cholesterol and lipoprotein metabolism in humans.

A multi-ethnic epigenome-wide association study of leukocyte DNA methylation and blood lipids

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Genetic Epidemiology in Latin America: Identifying Strong Genetic Proxies for Complex Disease Risk Factors

Epidemiology seeks to determine the causal effects of exposures on outcomes related to the health and wellbeing of populations. Observational studies, one of the most commonly used designs in epidemiology, can be biased due to confounding and reverse causation, which makes it difficult to establish causal relationships. In recent times, genetically informed methods, like Mendelian randomization (MR), have been developed in an attempt to overcome these disadvantages. MR relies on the association of genetic variants with outcomes of interest, where the genetic variants are proxies or instruments for modifiable exposures. Because genotypes are sorted independently and at random at the time of conception, they are less prone to confounding and reverse causation. Implementation of MR depends on, among other things, a strong association of the genetic variants with the exposure, which has usually been defined via genome-wide association studies (GWAS). Because GWAS have been most often carried out in European populations, the limited identification of strong instruments in other populations poses a major problem for the application of MR in Latin America. We suggest potential solutions that can be realized with the resources at hand and others that will have to wait for increased funding and access to technology.

Genome-wide association study of blood lipids in Indians confirms universality of established variants

Lipids foster energy production and their altered levels have been coupled with metabolic ailments. Indians feature high prevalence of metabolic diseases, yet uncharacterized for genes regulating lipid homeostasis. We performed first GWAS for quantitative lipids (total cholesterol, LDL, HDL, and triglycerides) exclusively in 5271 Indians. Further to corroborate our genetic findings, we investigated DNA methylation marks in peripheral blood in Indians at the identified loci (N = 233) and retrieved gene regulatory features from public domains. Recurrent GWAS loci-CELSR2, CETP, LPL, ZNF259, and BUD13 cropped up as lead signals in Indians, reflecting their universal applicability. Besides established variants, we found certain unreported variants at sub-genome-wide level-QKI, REEP3, TMCC2, FAM129C, FAM241B, and LOC100506207. These variants though failed to attain GWAS significance in Indians, but largely turned out to be active CpG sites in human subcutaneous adipose tissue and showed robust association to two or more lipid traits. Of which, QKI variants showed significant association to all four lipid traits and their designated region was observed to be a key gene regulatory segment denoting active transcription particularly in human subcutaneous adipose tissue. Both established and novel loci were observed to be significantly associated with altered DNA methylation in Indians for specific CpGs that resided in key regulatory elements. Further, gene-based association analysis pinpointed novel GWAS loci-LINC01340 and IQCJ-SCHIP1 for TC; IFT27, IFT88, and LINC02141 for HDL; and TEX26 for TG. Present study ascertains universality of selected known genes and also identifies certain novel loci for lipids in Indians by integrating data from various levels of gene regulation.

An exome wide association study of pulmonary tuberculosis patients and their asymptomatic household contacts

Tuberculosis is a leading cause of death in India. To identify genetic variants associated with susceptibility or resistance to Mycobacterium tuberculosis infection, we have performed an exome-wide association study with 0.2 million exonic variants among 119 pairs of tuberculosis patients and their clinically asymptomatic household contacts. The strongest association was identified for rs61104666[A], a synonymous variant (p.E292E) of exon 5 of the gene SIGLEC15 (OR = 2.4, p = 1.49 × 10^-5). We also found association of non-coding variants in the 3'UTR region of a gene encoding the class II human leukocyte antigens (HLAs), HLA-DRA. rs13209234[A] (minor allele frequency (MAF) = 13.8%) (OR = 0.35, P = 2.5 × 10^-4) and rs3177928[A] (minor allele frequency (MAF) = 13.7%) (OR = 0.35, P = 3.3 × 10^-4) were associated with protection from tuberculosis. These two SNPs, rs13209234 and rs3177928, are in complete linkage disequilibrium. These associations remained valid when additional data on freshly recruited individuals were jointly analyzed on 250 patient-control pairs. The identified gene, HLA-DRA, suggest involvement of immune regulation, indicating pathways associated with antigen presentation in tuberculosis infection.