Sequencing of candidate genes in Dominican families implicates both rare exonic and common non-exonic variants for carotid intima-media thickness at bifurcation

Rare variants in previously identified linkage regions associated with carotid plaque in Dominican Republic families

Carotid plaque is a subclinical measure of atherosclerosis. We have previously shown measures of carotid plaque to be heritable in a sample of 100 Dominican families and found evidence for linkage and association of common variants (CVs) on 7q36, 11p15, 14q32 and 15q23 with plaque presence. Our current study aimed to refine these regions further and identify rare variants (RVs) influencing plaque presence. Therefore, we performed targeted sequencing of the one LOD unit down region on 7q36, 11p15, 14q32 and 15q23 in 12 Dominican families with evidence for linkage to plaque presence. Gene-based RV analyses were performed using the Sequence Association Test for familial data (F-SKAT) under two filtering algorithms; 1. all exonic RVs and 2. non-synonymous RVs. Replication analyses were performed using a sample of 22 Dominican families and 556 unrelated Dominicans with Exome Array data. To identify additional non-synonymous RVs influencing plaque, we looked for co-segregation of RVs with plaque in each of the sequenced families. Our most strongly associated gene with evidence for replication was AMPD3 which showed suggestive association with plaque presence in the sequenced families (exonic RV p = 0.003, nonsynonymous RV p = 0.005) and replication families (exonic RV p = 0.04, nonsynonymous RV p = 0.02). Examination of the sequenced family pedigrees revealed two missense variants on chromosome 11 which co-segregated with plaque presence in one of our families; rs61751342 (located in DENND2B), and rs61760882 (located in RNF141). The rs61751342 missense variant is an eQTL for SCUBE2 in the atrial appendage. Notably, SCUBE2 encodes a protein which interacts with vascular endothelial growth factor (VEGF) receptor 2 to regulate VEGF-induced angiogenesis, thus providing biologic plausibility for this gene in atherosclerosis. In conclusion, using targeted sequencing of previously-identified linkage regions, we have identified suggestive evidence for the role of RVs in carotid plaque pathogenesis.

A host lipase prevents lipopolysaccharide-induced foam cell formation

Although microbe-associated molecular pattern (MAMP) molecules can promote cholesterol accumulation in macrophages, the existence of a host-derived MAMP inactivation mechanism that prevents foam cell formation has not been described. Here, we tested the ability of acyloxyacyl hydrolase (AOAH), the host lipase that inactivates gram-negative bacterial lipopolysaccharides (LPSs), to prevent foam cell formation in mice. Following exposure to small intraperitoneal dose(s) of LPSs, Aoah -/- macrophages produced more low-density lipoprotein receptor and less apolipoprotein E and accumulated more cholesterol than did Aoah +/+ macrophages. The Aoah -/- macrophages also maintained several pro-inflammatory features. Using a perivascular collar placement model, we found that Aoah -/- mice developed more carotid artery foam cells than did Aoah +/+ mice after they had been fed a high fat, high cholesterol diet, and received small doses of LPSs. This is the first demonstration that an enzyme that inactivates a stimulatory MAMP in vivo can reduce cholesterol accumulation and inflammation in arterial macrophages.

Targeted sequencing of linkage region in Dominican families implicates PRIMA1 and the SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus in carotid-intima media thickness and atherosclerotic events

Carotid intima-media thickness (cIMT) is a subclinical marker for atherosclerosis. Previously, we reported a quantitative trait locus (QTL) for total cIMT on chromosome 14q and identified PRiMA1, FOXN3 and CCDC88C as candidate genes using a common variants (CVs)-based approach. Herein, we further evaluated the genetic contribution of the QTL to cIMT by resequencing. We sequenced all exons within the QTL and genomic regions of PRiMA1, FOXN3 and CCDC88C in Dominican families with evidence for linkage to the QTL. Unrelated Dominicans from the Northern Manhattan Study (NOMAS) were used for validation. Single-variant-based and gene-based analyses were performed for CVs and rare variants (RVs). The strongest evidence for association with CVs was found in PRiMA1 (p = 8.2 × 10⁻⁵ in families, p = 0.01 in NOMAS at rs12587586), and in the five-gene cluster SPATA7-PTPN21-ZC3H14-EML5-TTC8 locus (p = 1.3 × 10⁻⁴ in families, p = 0.01 in NOMAS at rs2274736). No evidence for association with RVs was found in PRiMA1. The top marker from previous study in PRiMA1 (rs7152362) was associated with fewer atherosclerotic events (OR = 0.67; p = 0.02 in NOMAS) and smaller cIMT (β = −0.58, p = 2.8 × 10⁻⁴ in Family). Within the five-gene cluster, evidence for association was found for exonic RVs (p = 0.02 in families, p = 0.28 in NOMAS), which was enriched among RVs with higher functional potentials (p = 0.05 in NOMAS for RVs in the top functional tertile). In summary, targeted resequencing provided validation and novel insights into the genetic architecture of cIMT, suggesting stronger effects for RVs with higher functional potentials. Furthermore, our data support the clinical relevance of CVs associated with subclinical atherosclerosis.

Sequencing of Linkage Region on Chromosome 12p11 Identifies PKP2 as a Candidate Gene for Left Ventricular Mass in Dominican Families

Increased left ventricular mass (LVM) is an intermediate phenotype for cardiovascular disease (CVD) and a predictor of stroke. Using families from the Dominican Republic, we have previously shown LVM to be heritable and found evidence for linkage to chromosome 12p11. Our current study aimed to further characterize the QTL by sequencing the 1 LOD unit down region in 10 families from the Dominican Republic with evidence for linkage to LVM. Within this region, we tested 5477 common variants [CVs; minor allele frequency (MAF) ≥5%] using the Quantitative Transmission-Disequilibrium Test (QTDT). Gene-based analyses were performed to test rare variants (RVs; MAF < 5%) in 181 genes using the family-based sequence kernel association test. A sample of 618 unrelated Dominicans from the Northern Manhattan Study (NOMAS) and 12 Dominican families with Exome Array data were used for replication analyses. The most strongly associated CV with evidence for replication was rs1046116 (Discovery families P = 9.0 × 10-4; NOMAS P = 0.03; replication families P = 0.46), a missense variant in PKP2 In nonsynonymous RV analyses, PKP2 was one of the most strongly associated genes (P = 0.05) with suggestive evidence for replication in NOMAS (P = 0.05). PKP2 encodes the plakophilin 2 protein and is a desmosomal gene implicated in arrythmogenic right ventricular cardiomyopathy and recently in arrhythmogenic left ventricular cardiomyopathy, which makes PKP2 an excellent candidate gene for LVM. In conclusion, sequencing of our previously reported QTL identified common and rare variants within PKP2 to be associated with LVM. Future studies are necessary to elucidate the role these variants play in influencing LVM.

Rare Variants in NOD1 Associated with Carotid Bifurcation Intima-Media Thickness in Dominican Republic Families

Cardiovascular disorders including ischemic stroke (IS) and myocardial infarction (MI) are heritable; however, few replicated loci have been identified. One strategy to identify loci influencing these complex disorders is to study subclinical phenotypes, such as carotid bifurcation intima-media thickness (bIMT). We have previously shown bIMT to be heritable and found evidence for linkage and association with common variants on chromosome 7p for bIMT. In this study, we aimed to characterize contributions of rare variants (RVs) in 7p to bIMT. To achieve this aim, we sequenced the 1 LOD unit down region on 7p in nine extended families from the Dominican Republic (DR) with strong evidence for linkage to bIMT. We then performed the family-based sequence kernel association test (famSKAT) on genes within the 7p region. Analyses were restricted to single nucleotide variants (SNVs) with population based minor allele frequency (MAF) <5%. We first analyzed all exonic RVs and then the subset of only non-synonymous RVs. There were 68 genes in our analyses. Nucleotide-binding oligomerization domain (NOD1) was the most significantly associated gene when analyzing exonic RVs (famSKAT p = 9.2x10^-4; number of SNVs = 14). We achieved suggestive replication of NOD1 in an independent sample of twelve extended families from the DR (p = 0.055). Our study provides suggestive statistical evidence for a role of rare variants in NOD1 in bIMT. Studies in mice have shown Nod1 to play a role in heart function and atherosclerosis, providing biologic plausibility for a role in bIMT thus making NOD1 an excellent bIMT candidate.