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.

LVM in the discovery cohorts of each study, none achieved replication in their respective validation cohorts. Additional variants have been implicated in LVM through other genome-wide association studies, although none achieved statistical significance (Vasan et al. 2007;Arnett et al. 2009;Harper et al. 2013).
In a previous genome-wide linkage study, we identified a novel region on chromosome 12p11 linked to LVM using extended families from the Dominican Republic (MLOD = 3.11), particularly in families with higher waist circumference (MLOD = 4.45) (Wang et al. 2009). We further characterized this region by investigating common variants in a population-based cohort of Caribbean Hispanics and identified several candidates, including SOX5 (Della- Morte et al. 2011). For this study we aimed to identify additional common and rare variants contributing to the 12p11 linkage signal by sequencing the 1 LOD down region on 12p11 in 10 extended Dominican families with strong evidence for linkage to this region (LOD $ 0.1). We validated the observed associations using both genome-wide and exome-wide genotype data from unrelated Dominicans in the Northern Manhattan Study (NOMAS) as well as exome-wide genotype data from 12 additional extended Dominican families with suggestive evidence for linkage to the 12p11 region.

Study samples
Individuals from both NOMAS, a population-based cohort, and the Family Study of Stroke Risk and Carotid Atherosclerosis, a family study consisting of select probands from NOMAS and their family members (Sacco et al. 2004), were used for the current study. Details of NOMAS and the Family Study have been published previously (Sacco et al. 2004(Sacco et al. , 2007. Briefly, a total of 3298 stroke-free population-based participants were enrolled in NOMAS from 1993 to 2001. The family study enrolled a subset of Caribbean Hispanic probands from NOMAS with a high risk of cardiovascular disease (Sacco et al. 2007) who could provide a family history, obtain family members' permission for research staff to contact them, and have at least three first-degree relatives able to participate. All probands were identified in northern Manhattan and family members were enrolled in New York at Columbia University and in the Dominican Republic (DR) at the Clinicas Corazones Unidos in Santo Domingo. All participants provided written informed consent and the study was approved by the Institutional Review Boards of Columbia University, University of Miami, the National Bioethics Committee, and the Independent Ethics Committee of Instituto Oncologico Regional del Cibao in the DR (20070478 and 20072012).
Twenty-two families with a family-specific LOD score .0.1 at the chromosome 12p11 QTL for LVM were selected for our current study; 10 families were sequenced as part of our discovery analyses (n = 180 individuals) ( Table 1) and 12 families were genotyped on the Exome Array for our replication analyses (n = 143 individuals) (Supplemental Material, Table S1). Details of the family selection have been published previously (Dueker et al. 2016). Additional replication analyses were performed using a sample of 618 unrelated Dominicans from NOMAS (Sacco et al. 2004) who were not enrolled in the family study.
Vascular risk factors, including body mass index (BMI), systolic blood pressure, and smoking status, were collected during a standardized interview (Elkind et al. 2006). Smoking status was defined as never vs. ever, systolic blood pressure was defined as the average of two separate systolic blood pressure measurements taken after rest, and presence of diabetes was defined as fasting blood glucose level $126 mg/dl or self-reported history of diabetes. All data collection procedures were standardized and identical across NOMAS and the Family Study.

Discovery sample sequencing and quality control
In the family sample, genomic DNA was isolated from whole blood. Targeted sequencing of the exons in 181 genes within the 1 LOD unit down region on 12p11 (chr12:24-51 Mb), as well as sequencing beyond the exons for candidate genes identified previously (ARID2, BICD1, BIN2, c12orf68, FAR2, RACGAP1, SLC38A1, SOX5) (Della-Morte et al. 2011), was performed using a customized Agilent SureSelect Enrichment kit. A detailed description of sequencing methods has been previously published (Wang et al. 2015). Briefly, DNA libraries were sequenced on an Illumina HiSeq2000 and the raw sequencing reads were aligned to the human reference sequence hg19 with the Burrows-Wheeler Aligner (Li and Durbin 2010). Variant calling was performed using the Genome Analysis ToolKit and potential functions of variants were annotated using ANNOVAR v. 2016Feb01 (Wang et al. 2010) and SeattleSequation 138.
Quality control was conducted at both variant and sample levels, as described previously (Wang et al. 2015;Dueker et al. 2016). Within each individual sample, variants with a depth ,4 or Phred-Like score ,100 were set as missing. Variants with VQSLOD ,24 and variants with call rate ,75% were removed from further analysis. Individuals with low concordance (,95%) between the sequencing data and available genotype data were removed (n = 3). Additionally, individuals missing LVM measures (n = 9) and/or covariate values (n = 1) were removed. For the remaining family study samples, pedigree structure was confirmed using the Graphical Relationship Representation software v. 1.2.1.41. Mendelian error checking was performed, and Mendelian errors were set to missing for all the variants called using PLATO v. 0.84 (Grady et al. 2010).

Replication sample genotyping and quality control
Exome Array: NOMAS participants and our 12 replication families were genotyped using the Illumina HumanExome-24v1_B Beadchip, at the Hussman Institute for Human Genomics in the Center for Genome Technology (Miami, FL). Our Exome Array included custom exonic variants selected on the basis of sequencing data obtained in the discovery family data set. Details of the variant selection have been described previously (Wang et al. 2015). A total of 4128 single nucleotide variants (SNVs) within our 12p11 region were available for quality control analyses.
Of the 659 NOMAS participants and 150 replication family members genotyped on the Exome Array, 99.8% had a genotype call rate .98%. A subset of NOMAS participants genotyped on the Exome Array also had Affymetrix 6.0 whole-genome genotype data available and all participants had high concordance with this additional data set ($96%). We removed seven NOMAS individuals and five replication family members due to unexpected duplication or relatedness, gender discrepancy, and low call rate (,98%). For our current study, individuals missing LVM measures and/or covariate values were removed (n = 38 NOMAS participants; n = 2 replication family members), leaving us with a final sample of 618 Dominican NOMAS participants and 143 replication family members. At the variant level, we removed SNVs with call rate ,95% (n = 3 in NOMAS; n = 7 in the replication families) and monomorphic SNVs, leaving us with 1419 exonic rare single nucleotide variants (RVs) in our region for Exome Array analysis in NOMAS and 842 exonic RVs in the replication families. Mendelian error checking was performed in the replication families and Mendelian errors were set to missing for all the variants called using PLATO v. 0.84 (Grady et al. 2010).
Affymetrix 6.0 genotyping chip: In addition to Exome Array data, many of our NOMAS participants also had Affymetrix 6.0 wholegenome genotyping data available. These data were used for our common single nucleotide variant (CV) analyses since only 10.5% of CVs identified in the discovery families were available on the Exome Array (n = 574). Details of our genotyping and QC have been reported previously (Della-Morte et al. 2011). These data were imputed using the 1000 Genomes phase 1, version 3 reference panel with IMPUTE2 v.2.2.2 (Howie et al. 2009). Variants with INFO #0.4 were removed from analyses.

Statistical methods
Family-based discovery analyses: As in previous analyses, LVM was natural log transformed and multiplied by 10 to ensure it was normally distributed and properly scaled for analyses in SOLAR 6.6.2 (Wang et al. 2009;Della-Morte et al. 2011). Common and rare SNVs were defined based on frequencies from Dominican NOMAS participants, as described previously (Wang et al. 2015). SNVs were classified as common if they had MAF $5% and rare if they had MAF ,5% or could not be imputed efficiently (INFO #0.4) in NOMAS Dominicans. Analyses in the families were performed using the sequencing data. Single-variant analyses were performed for CVs using the Quantitative Transmission-Disequilibrium Test (QTDT), implemented in SOLAR. Adjustment was made for sex, BMI, systolic blood pressure, and smoking status. Covariates were identified using a polygenic screen implemented in SOLAR, with covariates having P , 0.1 included in analyses. CVs with P , 9.13 · 10 26 were considered significant based on a Bonferroni correction of 5477 tests.
To evaluate the contribution of our most strongly associated CV to our linkage results, linkage analyses were performed in our combined sample of sequenced and replication families in SOLAR following our previously detailed protocol (Wang et al. 2009). Briefly, linkage analyses were run with and without CV genotype as a covariate and a likelihood ratio test was performed to determine if the LOD score significantly decreased after conditioning on the CV.
Gene-based analyses were performed for RVs using the Family SNPset (Sequence) Kernel Association Test (Fam-SKAT) v. 1.8 (Chen et al. 2013), adjusting for the same covariates included in our single-variant analyses. We employed two different gene-based analyses based on annotation from ANNOVAR v. 2016Feb01 (Wang et al. 2010) and Seat-tleSequation 138: exonic RVs (UTR3, UTR5, synonymous, missense, nonsense, or splice-site variants) and a subset of nonsynonymous RVs only (missense, nonsense, or splice-site variants). These analyses were restricted to genes with $2 polymorphic variants and a P , 1.83 · 10 24 was considered significant based on a Bonferroni correction of 273 tests (181 exonic RV genes and 92 nonsynonymous RV genes).
n Using SAS v. 9.3, we computed the residual LVM value after adjusting for the associated risk factors to better visualize the distribution of RVs in relation to LVM.
Replication analyses: LVM was natural log transformed and multiplied by 10 to ensure a normal distribution. We additionally removed outliers falling 3 SD above or below the mean LVM value. Single-variant CV analyses were performed in the NOMAS sample using linear regression, implemented in PLINK v. 1.7. Variants were coded additively. Replication was defined as CVs with P , 0.05 and having the same direction of effect. Affy 6.0 genotyping data were used for these analyses.
Gene-based RV association analyses using Exome Array data were performed using SKAT-O v. 1.1.2 in the NOMAS sample and Fam-SKAT v. 1.8 in the replication family sample. These analyses were performed using the same two filtering algorithms employed in the discovery analyses; all exonic RVs and then a subset of nonsynonymous RVs only. Analyses were restricted to genes with $2 polymorphic variants. Replication was defined as genes with P , 0.05. All association analyses were adjusted for the same covariates included in our discovery analyses. NOMAS analyses were additionally adjusted for factors that were associated with LVM (P , 0.1) in this sample: diabetes, time between LVM measurement and baseline patient assessment, waist hip ratio, and the first principal component obtained via principal components analyses implemented in Eigenstrat to account for population substructure. Details of our Eigenstrat analysis and resulting principal components can be found in Figure S1 and File S1.

Data availability
The data that support the findings of this study are currently being uploaded to dbGaP.

Participant and SNV characteristics
A total of 180 individuals in 10 Dominican families were included in our discovery analyses. Characteristics of these 10 families are summarized in Table 1. Family size ranged from 10 to 35 individuals and familyspecific LOD scores for the 12p11 region ranged from 0.14 to 0.94. LVM residual values ranged from 25.62 to 6.89. Mean age was between 40 and 54 yr and mean BMI was similar across the 10 families. Variability was seen with respect to percent of participants with diabetes and hypertension in each family (Table 1).
Within the 10 families, sequencing identified 5473 CVs and 10,167 RVs. Among RVs, 24.6% (n = 2503) were novel and 22.8% (n = 2318) were classified as exonic. A total of 20.2% of exonic RVs (n = 468 variants) were either missense, nonsense, or splice-site variants ( Table  2). In the NOMAS sample, 96.5% of the CVs and 44.4% of the exonic RVs identified through sequencing were available for replication analyses. An additional 572 exonic RVs on the Exome Array were polymorphic in the NOMAS sample but were monomorphic in the discovery family sample and were included in gene-based replication analyses. In the replication families, 29.6% of the exonic RVs identified through sequencing were available for analysis in addition to 153 exonic RVs that were on the Exome Array and polymorphic in the replication families but monomorphic in the discovery family sample.

CV association results
Results from the CV analyses are shown in Figure 1 and our most strongly associated variants with suggestive evidence for replication are listed in Table 3. While no CV reached peak-wide significance (P , 9.2 · 10 26 ), rs1046116 achieved a P-value of 9.0 · 10 24 and showed evidence for replication in NOMAS (P = 0.03). This variant was associated with decreased LVM and is a missense variant located within PKP2. To investigate this association further, we tested rs1046116 for association with LVM in our replication families and found no evidence for replication (P = 0.46). To evaluate the contribution of rs1046116 to our linkage signal, we performed linkage analyses in our combined discovery and replication family sample and observed a change in LOD score from 9.09 to 8.30 (P , 0.0001) with rs1046116 adjustment, indicating that this CV accounted for some of the linkage signal at the 12p11 QTL. Eight additional variants showed suggestive evidence for association in both the families and NOMAS (P , 0.05). The minor allele of these variants, except for rs11168985 (located in CPNE8), was associated with increased LVM. n

RV association results
Gene-based RV analyses were performed using two filtering algorithms; first, analyzing all exonic RVs (n = 181 genes) and second, restricting analyses to nonsynonymous RVs (n = 92 genes). A total of 54 genes showed suggestive evidence for association in the discovery families under at least one filtering algorithm (P , 0.05). Results for these genes are shown in Table S2. Among these top genes, we observed at least suggestive evidence for replication of five genes in NOMAS and three genes in the replication families (P , 0.10). Summary results for these genes are shown in Table 4. In exonic RV analyses, NELL2 was our most strongly associated gene (P = 2.2 · 10 24 ), although this finding did not replicate in NOMAS (P = 0.53) or the replication families (P = 0.18) (Table S2). While no gene achieved replication in both replication samples, our most strongly associated gene with evidence for replication in NOMAS was ALG10B (P = 0.02 in the discovery families; P = 0.04 in NOMAS; P = 0.86 in the replication families). When analyses were restricted to nonsynonymous RVs, this association became slightly attenuated in the discovery families (P = 0.06), and was even stronger in NOMAS (P = 0.03). GXYLT1 was our most strongly associated gene with evidence for replication in the replication families (P = 0.04 in the discovery families; P = 0.004 in the replication families; P = 0.36 in NOMAS). In nonsynonymous RV analyses, IFLTD1 showed the strongest evidence for association in the discovery families (P = 0.003). However, since no variants met our inclusion criteria in NOMAS or the replication families, replication was not possible. NELL2 was our second most strongly associated gene (P = 0.006), and replicated in the replication families (P = 0.02) but not NOMAS (P = 1.00) (Table S2). Two genes, PKP2 and SLC2A13, were moderately associated in the discovery families and showed at least suggestive evidence for replication in NOMAS (P , 0.10), although they did not replicate in the replication families. Details of the nonsynonymous RVs in PKP2, SLC2A13, and NELL2 are shown in Table 5. Within the discovery families, six nonsynonymous variants in PKP2 were observed, all of which were missense variants. These variants were carried by individuals within five families: F2783, F3561, F3719, F5103, and F5275. In families F2783 and F3719, carriers of a rare PKP2 allele had higher average LVM residual compared to noncarriers. In the remaining families, carriers of a rare PKP2 allele had lower average LVM residual compared to noncarriers (Figure 2). When examining each missense variant individually, we observed that carriers of the rare allele of rs143004808, rs146882581, rs151264959, and rs62001016 had higher average LVM residual compared to noncarriers. In contrast, carriers of the rare allele of rs200069860 and rs75909145 had lower average LVM residual compared to noncarriers (Figure 3).
Since our most strongly associated CV was a missense variant in PKP2, we performed an additional gene-based analysis on all nonsynonymous variants in PKP2 (common and rare) which resulted in a slightly stronger association in the discovery families (P = 0.04, seven variants) and no association in NOMAS (P = 0.41, 12 variants) or the replication families (P = 0.98, seven variants).

Candidate variant and gene association results
We additionally tested whether candidate CVs identified in our previous study of LVM (Della-Morte et al. 2011) (n = 33 CVs in 14 genes), a n n Table 4 Gene-based association results for rare variant analysis in the chr12p11 region, for genes with P < 0.05 in the Dominican families and P < 0.10 in the replication families or NOMAS  study which primarily investigated CV associations with LVM in all Caribbean Hispanics, were associated in our current study which focused specifically on Dominicans. Association results for these CVs are summarized in Table S3. Eleven of the previously reported CVs, located in SLC38A1, BICD1, RACGAP1, ARID2, FAR2, and c12orf68, were found in our discovery families. These CVs included three intronic variants in SLC38A1 which showed suggestive evidence for association in the discovery families and NOMAS, although the direction of effect differed (rs11183394, rs6582621, rs7133522). One other CV in SLC38A1, rs7956629, was moderately associated in the discovery families (P = 0.02), but did not replicate in NOMAS (P = 0.29). No other CVs were associated with LVM in the families. We also performed gene-based RV analyses on 14 genes containing the previously reported candidate CVs. Results are summarized in Table S4. Exonic RV analyses identified PDZRN4 (P = 0.03) and ARID2 (P = 0.04) to be moderately associated in the discovery families, although neither replicated in NOMAS or the replication families (P . 0.05). One additional gene, c12orf68, was moderately associated in the discovery families (P = 0.03); however, replication was not possible as no variants met our inclusion criteria in either replication sample. In nonsynonymous RV analyses, SLC2A13 was the most strongly associated gene and was a top gene in our peak-wide nonsynonymous RV analyses.

DISCUSSION
Building upon our previous studies which identified a region on 12p11 to be linked with LVM in Dominicans, we performed targeted rese-quencing in 10 extended families with evidence for linkage to refine the region and identify LVM susceptibility genes. Through these resequencing efforts, we found suggestive evidence for both common and rare variants within 12p11 to influence LVM in Dominican families. Common variant analyses revealed rs1046116, a missense variant in PKP2, to be the most strongly associated variant with evidence for replication in our population-based sample of Dominicans. When investigating this association further by testing rs1046116 for association in our 12 additional families with suggestive evidence for linkage to the region, we observed no association; however, these families generally showed weaker evidence for linkage to this region than the 10 families sequenced for our discovery analyses. Further, although this association could not be replicated in our replication families, linkage analyses in our combined family sample revealed that rs1046116 significantly contributed to the original linkage signal as evidenced by a decrease in LOD score from 9.09 to 8.30 in these families. This variant was associated with decreased LVM, suggesting that the G allele of rs1046116 may be protective against the development of left ventricular hypertrophy.
Interestingly, we also found suggestive evidence that RVs within PKP2 influence LVM. In nonsynonymous RV analyses, PKP2 was one of the most strongly associated genes with suggestive evidence for replication, thus suggesting that both common and rare nonsynonymous variants, particularly missense variants, may influence LVM. Indeed, when a gene-based analysis was performed on all nonsynonymous variants, common and rare, the association between PKP2 and LVM became slightly stronger in the families, although the association diminished in NOMAS.
To explore the associations underlying PKP2 further, we performed an in-depth characterization of the individual rare nonsynonymous variants within this gene. Three missense variants were predicted to be probably-damaging by PolyPhen and showed strong evidence of being conserved (rs200069860, rs151264959, rs143004808). Of particular interest was rs143004808, which was observed in F3719 and had a Combined Annotation Dependent Depletion (CADD) score of 33, indicating that this variant is in the top 0.1% of variants with respect to deleteriousness. Two nonsense variants (located at base pair positions 32,975,421 and 33,021,968) also had similar CADD scores and were observed in the NOMAS sample.
Missense variants in PKP2, including several of the ones identified in our study, have been previously implicated in cardiac phenotypes (Cerrone et al. 2014;Peters 2014;van der Zwaag et al. 2009;Saguner et al. 2015), making PKP2 an excellent candidate gene for LVM. PKP2 encodes plakophilin 2, a member of one of three major protein families found in desmosomes (Getsios et al. 2004). Desmosomes are protein structures in cell membranes that maintain adhesion between neighboring cells and serve as anchoring sites for intermediate filaments. They are found in tissues that experience mechanical stress, including the myocardium (Getsios et al. 2004;Desai et al. 2009). Mutations in PKP2 are known to play a role in arrythmogenic cardiomyopathies (AC), most notably arrythmogenic right ventricular dysplasia/cardiomyopathy (ARD/C) (Gerull et al. 2004;Fernandez-Rosado et al. 2015;Groeneweg et al. 2015;Trenkwalder et al. 2015). Several studies have shown that PKP2 haploinsufficiency contributes to pathogenesis in AC (Joshi-Mukherjee et al. 2008;Hall et al. 2009;Kirchner et al. 2012;Rasmussen et al. 2014).
In characterizing the clinical course of ARD/C, one study found that 9% of ARD/C patients with a PKP2 mutation had left ventricular dysfunction and patients carrying more than one mutation were three times more likely to have left ventricular dysfunction compared to patients with one mutation . Studies have also provided further evidence suggesting that PKP2 mutations may impact left ventricular cardiomyopathies (Horimoto et al. 2000;Hamid et al. 2002;Sen-Chowdhry et al. 2007;Saguner et al. 2015), including a recent case report that identified a pathogenic PKP2 deletion in two siblings with left ventricular noncompaction cardiomyopathy (Ramond et al. 2017). Together, these studies support the biologic plausibility for a role of PKP2 in LVM.
In addition to PKP2, our analyses revealed nonsynonymous variants within SLC2A13 to be suggestively associated in both our discovery families and the NOMAS sample. Included in our nonsynonymous analyses were two missense variants observed in both the families and NOMAS (rs139518863 and rs186341127), as well as one missense variant observed in NOMAS only (rs146020551). Rs139518863 has strong evidence of being evolutionarily conserved (Genomic Evolutionary Rate Profiling (GERP) score = 4), is predicted to be possibly-damaging according to PolyPhen, and has a CADD score of 19.66 indicating that this variant is in the top 10-1% of all variants with respect to deleteriousness. Similar characteristics were observed for rs146020551. Interestingly, our previous association study found CVs in SLC2A13 to be associated with LVM in Caribbean Hispanics (Della-Morte et al. 2011). These findings, in combination with results from our current study, suggest that variants in SLC2A13 may be involved in LVM in Dominicans and, more broadly, Caribbean Hispanics.
Additional candidate genes identified through our resequencing study were GXYLT1 and NELL2. GXYLT1 was identified in our exonic RV analyses and encodes an enzyme that adds xylose to O-glucose residues bound to epidermal growth factor repeats of Notch proteins (Sethi et al. 2010). Notch signaling plays a role in cardiac development and mutations within Notch signaling genes have been associated with cardiac structural abnormalities including left ventricular outflow tract abnormalities, making GXYLT1 an excellent candidate gene for LVM (Nemir and Pedrazzini 2008;Padang et al. 2012;Penton et al. 2012). NELL2 was one of our most strongly associated genes in RV analyses and encodes the neural epidermal growth factor-like 2 protein which is largely expressed in brain but is also expressed in hematopoietic cells (Luce and Burrows 1999). Its role in LVM is unknown.
When evaluating the results from our study, there are several limitations which should be noted. First, only a subset of families was sequenced; however, these were the families with the strongest evidence for linkage to the 12p11 region. Second, sequencing was performed primarily on exons, therefore missing noncoding regions affecting LVM. However, for candidate genes previously implicated in LVM, sequencing beyond the exons was performed to allow for identification of noncoding variants. Third, the RV replication analyses were performed using Exome Array data, thereby limiting our analyses to only those variants included on the Exome Array. However, variants on the Exome Array were selected to be functional (primarily missense variants) and had to be observed at least three times in at least two people (http://genome.sph. umich.edu/wiki/Exome_Chip_Design), and we added custom content within the 12p11 region to our Exome Array. Fourth, our study was limited to individuals from the Dominican Republic and may not be generalizable to other Hispanic or non-Hispanic populations. Fifth, due to our sample sizes, we had limited power to detect RVs in our study.
In conclusion, our current targeted resequencing study, in combination with our previous studies, shows evidence of a role for both common and rare variants within the 12p11 region in LVM pathogenesis, particularly missense variants within PKP2. Functional studies are needed to elucidate the mechanism underlying the association of the implicated genes with LVM.