Genome-Wide DNA Methylation Analysis and Epigenetic Variations Associated with Congenital Aortic Valve Stenosis (AVS)

Abstract

Congenital heart defect (CHD) is the most common cause of death from congenital anomaly. Among several candidate epigenetic mechanisms, DNA methylation may play an important role in the etiology of CHDs. We conducted a genome-wide DNA methylation analysis using an Illumina Infinium 450k human methylation assay in a cohort of 24 newborns who had aortic valve stenosis (AVS), with gestational-age matched controls. The study identified significantly-altered CpG methylation at 59 sites in 52 genes in AVS subjects as compared to controls (either hypermethylated or demethylated). Gene Ontology analysis identified biological processes and functions for these genes including positive regulation of receptor-mediated endocytosis. Consistent with prior clinical data, the molecular function categories as determined using DAVID identified low-density lipoprotein receptor binding, lipoprotein receptor binding and identical protein binding to be over-represented in the AVS group. A significant epigenetic change in the APOA5 and PCSK9 genes known to be involved in AVS was also observed. A large number CpG methylation sites individually demonstrated good to excellent diagnostic accuracy for the prediction of AVS status, thus raising possibility of molecular screening markers for this disorder. Using epigenetic analysis we were able to identify genes significantly involved in the pathogenesis of AVS.

About the Authors:

Uppala Radhakrishna

* E-mail: [email protected]

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Samet Albayrak

Affiliation: Department of Obstetrics and Gynecology, Wayne State University School of Medicine, Detroit, Michigan, United States of America

Zeynep Alpay-Savasan

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Amna Zeb

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Onur Turkoglu

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Paul Sobolewski

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Ray O. Bahado-Singh

Affiliation: Department of Obstetrics and Gynecology, Oakland University William Beaumont School of Medicine, Royal Oak, Michigan, United States of America

Introduction

Congenital heart defect (CHD) is the most common type of birth defect, and affects close to 1% of all live births [1, 2]. Over 50% of children born with CHD will have at least one invasive surgery in their lifetime, while many will have multiple surgeries. In the US, approximately 20–30% of all infants with CHD have valve defects [3, 4].

The etiology of CHD is multifactorial and complex. Approximately 30% of children with a chromosomal abnormality have associated CHD [5]; however, a majority of CHD cases have no identifiable genetic abnormality. Families with autosomal dominant [6], recessive [7] and X-linked [8] CHD have been reported; however, in the majority of cases genetic mechanisms have not heretofore been identified. It is postulated that sporadic CHD, which accounts for around 80% of all cases, arises from an interaction of genetic and environmental factors [9]. Gestational diabetes, dietary deficiency, medication use, maternal age, cigarette smoking, obesity, febrile illnesses in pregnancy, alcohol intake and viral infections have been reported as maternal-environmental risk factors for the development of CHD [10]. The exact mechanisms by which these factors cause CHD remain unknown.

At least 18 distinct types of congenital heart defects are recognized; among these aortic valve stenosis (AVS) is one of the most common and most serious valve disease that has strong genetic basis [11, 12]. Congenital AVS, defined as incomplete obstruction of the valve orifice, is an important category of structural heart defect, and occurs in 3–6% of such cases [13]. There is variability in both the site of obstruction and severity of the obstruction. Sites of obstruction are sub-classified as valvular, subvalvular and supravalvular [14]. About half of infants with severe AVS require surgery [15]. Mild aortic stenosis is difficult to detect prenatally; however, critical aortic stenosis can lead to left ventricular myocardial dysfunction with endocardial fibroelastosis, left atrial dilation and narrowing of the aortic root [16]. These changes can be a prelude to the development of hypoplastic left heart syndrome. The present standard for prenatal screening is mostly based the use of ultrasonography and there is no reliable biologic screening marker for any type of CHD [17].

Recently, some studies have identified an association between dietary folate supplementation and decreased CHD risk [18]. Deficiency in folate is known to result in hyperhomocysteinemia, another risk factor for CHD development [19]. Folate is a B-vitamin that supplies the methyl group for DNA methylation and for other methylation reactions in the cell. Such findings are consistent with a role for epigenetic DNA methylation in the development of CHD. In a recently conducted genome-wide DNA-methylation analysis, we reported evidence of significantly altered CpG methylation levels genome-wide in multiple categories of CHD compared to controls. The study included cases of hypoplastic left heart syndrome, Teratology of Fallot, ventricular septal defect, atrial septal defect and coarctation of the aorta [20].

The prenatal and newborn detection of CHD remains a significant challenge [21, 22]. We undertook a study to examine genome-wide DNA methylation patterns in newborns with AVS to identify genomic regions containing disease-related genes and epigenetic changes that may contribute to CHD pathophysiology. An important objective of the study was to identify DNA methylation biomarkers, serum molecules that could potentially be used in the future for risk estimation and detection of AVS.

Materials and Methods

Genomic DNA was obtained from neonatal dried blood spots using commercial DNA extraction kits (Qiagen QIAamp^®) according to manufacturer’s protocol. Many studies have reported genome-wide DNA methylation profiles from archived dried blood spots using the Infinium HumanMethylation450 BeadChip as a suitable template [23, 24]. Blood spot specimens were collected previously for the mandated newborn screening and treatment program run by the Michigan Department of Community Health in the State of Michigan (MDCH). All specimens were collected between 24 and 79 hours after birth. Parents/legal guardians were aware at the time of blood collection that residual blood spots after clinical testing may be utilized for research pending review of such study requests by the MDCH. This study was approved by both the institutional review boards from William Beaumont Hospital and the MDCH. Limited demographic information was available for each subject including date of sample collection, maternal age and race, gestational age at delivery and newborn sex along with the type of CHD anomaly. Suspected or diagnosis-unknown AVS cases were excluded. Unaffected controls had no reported medical disorder and were matched for birth weight, gestational age at delivery, ethnicity, year of birth, and interval from specimen collection to testing. Our cohort included 24 AVS subjects and 24 controls. All specimens were de-identified by removal of further protected health information and researchers were masked to subject identity. Details of the case control cohort are available in S1 Table.

Genome-wide methylation analysis using the HumanMethylation450

Genome-wide methylation analysis was performed for 48 individuals (24 AVS subjects and 24 controls) using the HumanMethylation450, Illumina’s newest Infinium® HD BeadChip assay for methylation (Illumina, Inc., California, USA), which contains 485,577 methylation sites and requires only 500 ng of genomic DNA. These sites are equally distributed in the genome and represent 96% of RefSeq genes, 95% of CpG islands and an average of 17 CpG sites per covered gene region including the promoter, 5’UTR, coding, and 3’UTR regions. DNA methylation profiling using Illumina Infinium technology with peripheral blood lymphocytes has been used to identify CpG sites associated with disease states [25, 26].

The DNA samples were bisulfite converted using the EZ DNA Methylation-Direct Kit (Zymo Research, Orange, CA) according to the manufacturer´s protocol. The fluorescently stained BeadChips were imaged by the Illumina iScan. Prior to detailed bioinformatic and statistical analysis, data preprocessing and quality control was performed including examination of the background signal intensity of both affected and negative controls, the methylated and unmethylated signals, and ratio of the methylated and unmethylated signal intensities. The processing is done fully according to manufacturer's protocol and 99% of the CpG loci are determined unequivocally.

Statistical and Bioinformatic analysis

Genome-wide, gene-specific DNA methylation was measured using the Genome Studio methylation analysis package (Illumina). Following the pre-processing described above, a DNA methylation ß-value was assigned to each CpG site. Differential methylation was assessed by comparing the ß-values per individual nucleotide at each CpG site between AVS subjects and controls. In order to avoid potential confounding factors, probes associated with sex chromosomes and/or containing SNPs in the probe sequence (listing dbSNP entries near or within the probe sequence, i.e., within 10 bp of the CpG site) were excluded from further analysis [27–29]. Probes targeting CpG loci associated with SNPs near or within the probe sequence may influence corresponding methylated probes [30]. The remaining CpG sites were taken forward for analysis.

The most discriminating 59 CpG sites were selected based on the pre-set cutoff criteria of >1.5-fold increase and/or >1.5-fold decrease with p < 0.05. These sites include the highest 46 hypermethylated and top 13 demethylated CpG sites for further analysis. The identified genes are listed in Table 1. The 59 CpG sites, corresponding to 52 genes, were differentially methylated either in the coding and/or promoter regions and were subsequently used to generate a heatmap using the ComplexHeatmap (v1.6.0) R package (v3.2.2). We have used ward distance for the hierarchical clustering of samples [31] (Fig 1). Subsequently these 59 CpG sites were used to calculate receiver operating characteristic (ROC) curves and the area under the ROC curves (ROC AUC).

[Figure omitted. See PDF.]

Fig 1. Heatmap of unsupervised hierarchical clustering of AVS data on the basis of 59 differentially methylated CpG sites.

Unsupervised hierarchical clustering analysis is very popular in identifying methylation-defined patient sub-groups. Fig 1 displays CpG sites that are at least either 1.5 fold demethylated or 1.5 fold hypermethylated in the disease (AVS) condition (red colored squares) compared to normal subjects (blue colored squares). Differentially methylated CpG sites have been displayed in three clusters (row-wise). The figure also displays direction, fold change in disease, probe relationship and probe annotation of differentially methylated CpG sites. Red color in the heatmap indicates hyper DNA-methylation, and blue hypo DNA-methylation values.

https://doi.org/10.1371/journal.pone.0154010.g001

[Figure omitted. See PDF.]

Table 1. Highly Differentially Methylated CpG sites in Aortic Stenosis subjects.

Differentially methylated genes with Target ID, Gene ID, chromosome location, % methylation change and FDR p-value for each gene methylated. CpG sites with significant FDR p-value indicating methylation status and ROC AUC ≥0.75 appear to have a strong potential as diagnostic biomarkers for AVS.

https://doi.org/10.1371/journal.pone.0154010.t001

Validation of methylation status by bisulfite sequencing

To confirm the results obtained from the Illumina HumanMethylation450 arrays, we performed pyrosequencing analysis of bisulphite-converted DNA with appropriate oligos using the PyroMark Q24 System and advanced CpG Reagents (Qiagen ^®).

The p-value for methylation differences between case and normal groups at each locus was calculated as previously described [32]. Filtering criteria for p-values were set at <0.05 and also <0.01 in order to identify the most differentiating cytosines. P-values were calculated with and without Benjamini and Hochberg False Discovery Rate (FDR) correction for multiple testing [33]. The Benjamini and Hochberg correction tolerates more false positive genes than the Bonferroni correction.

Further analysis of the differentially methylated genes was conducted for potential biological significance. ROC curves and ROC AUC were calculated to determine the diagnostic accuracy of specific cytosine loci to differentiate AVS from control groups. Data were normalized using the Controls Normalization Method.

Gene ontology analysis and functional enrichment

The genes found to be differentially methylated (at FDR p-value < 0.01) were uploaded to the web-based functional annotation tool DAVID V67 (DAVID/EASE, WebGestalt) for Gene Ontology analysis [34, 35] including gene ID conversion, bio-pathways analysis, and the molecular functions of methylated and unmethylated regions. Literature data mining for co-occurrence of gene names and keywords of interest was performed using Chilibot. Only genes for which Entrez identifiers were available were further analyzed. Pathway analysis was carried out using Ingenuity pathway analysis (Ingenuity Systems). Over-represented canonical pathways, biological processes and molecular processes were identified.

Web Resources

The URLs for data presented herein are as follows:

Illumina: http://www.illumina.com/

Genome Studio: http://www.solexa.com/gsp/genomestudio_software.ilmn

Ensemble: http://www.ensembl.org/

UCSC: http://genome.ucsc.edu/

NCBI: http://www.ncbi.nih.gov/

DAVID/EASE: http://david.abcc.ncifcrf.gov/

Web Gestalt: http://genereg.ornl.gov/webgestalt/

Chilibot: www.chilibot.net

Ingenuity Systems: www.ingenuity.com

Online Mendelian Inheritance in Man (OMIM): http://www.ncbi.nlm.nih.gov/Omim/

Results

There were no differences in gestational age at birth in weeks: mean (SD) 38.75 (1.42) in AVS subjects vs. 38.88 (1.19) in controls (p = 0.743), nor in the timing of specimen collection after birth (in hours): mean (SD) 31.042 (11.86) in AVS subjects vs. 32.46 (8.62) in controls (p = 0.638). There were no variations in maternal age: mean (SD) 29.87 (4.56) in AVS subjects vs. 29.87 (4.56) years in controls (p-value 1.00). Finally, maternal race and newborn gender were matched for analysis. Of the 52 genes identified during genome-wide methylation analysis, hierarchical clustering analysis demonstrated ~10 as novel principal candidate genes that are commonly methylated and whose methylation was associated with altered gene expression in AVS individuals. Table 1 includes the most significantly differentially methylated CpG sites based on FDR-corrected p-values. The methylation status is represented as percentage methylation for a given probe in the sample. A positive ‘% m Change’ value indicates an average increase in methylation status in AVS subjects compared to control samples. Similarly, a negative ‘% m Change’ value indicates a decrease in methylation status in AVS subjects compared to controls. The p-value indicates the significance of the differential methylation levels. The University of California Santa Cruz (UCSC) gene name and genomic location of the C in the CG dinucleotide and the chromosome on which it is located as provided by Illumina are also shown in Table 1. The results obtained from the DAVID Pathway and Gene Ontology over-representation analysis for canonical pathways and for biological processes are presented in Table 2. Gene Set Enrichment analysis using multiple computational tools showed no significant functional enrichment due to the relatively small size of the gene list. Therefore Gene Ontology information for all genes given in the list was obtained and classified. DAVID pathway analysis software was used to identify molecular pathways associated with genes having differentially methylated CpG sites between AVS subjects and controls. Analysis was done on genes with at least one differentially methylated CpG site based on the uncorrected p-values.

[Figure omitted. See PDF.]

Table 2. Over-represented Gene Ontology Molecular Function and Biological Process categories as determined using DAVID Categories: AVS.

Biological Process and Metabolic Function categories for over-represented pathways determined using DAVID Pathway and Gene Ontology analysis

https://doi.org/10.1371/journal.pone.0154010.t002

In combination with the FDR p-value indicating methylation status, the area under the ROC curves can be used to distinguish AVS subjects from controls. All data cleaning and analysis was performed using R (version 3.2.3) and RStudio (version 0.99.489). The CpG sites corresponding to the 52 differentially methylated genes have a ROC AUC ≥0.50 including six CpG sites with ROC AUC ≥0.75: cg16464924 (AUC 075; 95% CI, 0.62 to 0.89), cg05890887 (AUC 0.80; 95% CI, 0.68 to 0.93), cg26196700 (AUC 0.75; 95% CI, 061 to 081), cg22469274 (AUC 0.76; 95% CI, 0.63 to 0.90), cg00994804 (AUC 0.76; 0.62 to 0.89), cg04836786 (AUC 0.76; 0.63 to 0.90). At each locus, the FDR p-value for methylation difference between AVS subjects and controls was highly significantly different. ROC curve analysis (Fig 2) narrowed down the number of markers commonly methylated in AVS for the development and further consideration of their role and validation and also possible use as biomarkers.

[Figure omitted. See PDF.]

Fig 2. Receiver operating characteristic (ROC) curve analysis of methylation profiles for four specific markers associated with aortic valve stenosis.

AUC: Area Under Curve; 95% CI: 95% Confidence Interval. Lower and upper confidence interval was given in parentheses. We have identified six CpG sites (cg16464924, cg05890887, cg26196700, cg22469274 cg00994804 cg04836786) among 52 differentially methylated genes that have ROC AUC ≥0.75. At each locus, the FDR p-value for methylation difference between AVS subjects and controls was highly significantly different. Due to figure resolution, we have included only for four markers.

https://doi.org/10.1371/journal.pone.0154010.g002

Biological processes and molecular function determination for these genes are shown in Table 2. Genes were further grouped according to their Gene Ontology-characterized function. Two genes were identified whose biological processes (APOA5 PCSK9) have known molecular function (AIMP1, PCSK9, TXNRD2, CLDN11, RUNX1, MCM6, APOA5 and PCSK9).

Discussion

We have presented genome-wide differential methylation analysis of AVS subjects. To our knowledge, this is the first study that has performed an array-based genome-wide high resolution DNA methylation in a cohort of AVS newborns. We found that epigenetic alteration of CpG sites exhibits a relationship with AVS phenotype, not only of those located within the CpG islands in the promoter region of genes, but also of those distributed throughout the gene including 5’ UTR, coding, and 3’UTR regions. Recent work by our group [20] found a strong association between CpG methylation changes and the presence of multiple different types of CHD. Another recently published study reported DNA methylation changes in patients with Tetralogy of Fallot (TOF) using MassARRAY-based quantitative methylation analysis [36]. That study was, however, limited to targeted promoter regions. Taken together, these studies provide evidence that DNA methylation may be critical in the genetic and environmental interactions underlying cardiac morphogenesis. We have identified 59 CpG sites based on either at least 1.5 fold demethylation or 1.5 fold hypermethylation in disease samples. We hypothesize from the present study that altered expression of one or a combination of the 52 genes corresponding to the 59 CpG sites identified by differential methylation analysis is likely to be responsible for AVS.

There is a well-established and known association between defects in receptor-mediated endocytosis and many diseases such as hypercholesterolemia, a condition characterized by very high levels of cholesterol in the blood. Higher levels of total cholesterol increase the risk of cardiovascular disease [37, 38] including aortic valve stenosis [39] and a high risk for heart disease even in childhood. The report by Chui et al [40] also showed significantly higher serum total cholesterol concentrations in patients with aortic stenosis than in controls. However, molecular mechanisms responsible for the association of hypercholesterolemia with CHD in offspring remain unknown. In the present study, the identification of multiple genes modified in AVS indicates that the phenotype is a complex trait.

DUSP27 (Dual Specificity Phosphatase 27 (Putative)) is a protein coding gene that belongs to the dual-specificity phosphatase (DUSPs) family of proteins. These proteins play important roles in a number of cell types and events, but are preferentially expressed in the heart and skeletal muscles. Li [41] conducted a Dusp27 knockdown experiment using siRNA to study its effect on myogenic differentiation. These experiments identified this gene’s important role during muscle and heart development. Together, these findings indicate that Dusp27 may have a more important function in these two tissues. Arrington et al. [42] conducted exome sequencing of DNA samples from multiple affected family members with diverse CHDs and identified gene variations in DUSP27 among the 18 variations they noted.

APOA5 (OMIM 606368) or apolipoprotein A5 is associated with plasma triglyceride levels, a major risk factor for coronary heart disease. Earlier association studies were conducted to examine the links between high-density lipoprotein genetics of various genes, including the APOA5 gene, and aortic valve stenosis risk [43]. However due to limited sample size, their study did not yield strong associations. A significant correlation between the APOA5 gene polymorphism and the levels of plasma high-density lipoprotein-cholosteal and increased risk for cardiovascular disease was previously identified [44, 45].

Another important gene identified by Gene Ontology analysis is proprotein convertase, subtilisin/Kexin-Type 9 (PCSK9, OMIM 607786). PCSK9 is a serine protease that reduces both hepatic and extrahepatic low-density lipoprotein (LDL) receptor levels and increases plasma LDL cholesterol. Mutation in this gene leads to an increase in the LDL cholesterol level. Cohen et al [46] studied a large population and identified PCSK9 gene variations associated with differing plasma levels of LDL cholesterol and the risk of coronary heart disease, including myocardial infarction and mortality from coronary disease.

Two other genes displaying altered methylation are Runt-related transcription factor 1 (RUNX1; OMIM 151385) on chromosome 21q22.12, and Thioredoxin reductase 2 (TXNRD2; OMIM 606448) on chromosome 22q11.21. Deletion of chromosome 21q21.1–22.12 including the RUNX1 gene has been associated with multiple congenital anomalies and congenital heart defects including aortic stenosis [47, 48]. Thioredoxin reductase is directly involved in reducing proteins such as insulin. Experiments in mice demonstrated TrxR2 is required for proper heart development [49]. Chromosome micro-deletions (del22q11) involving TXNRD2 have also been associated with various congenital heart defects including aortic valve stenosis [50]. Helicase-Like transcription factor (HLTF; OMIM 603257), a gene that maps to chromosome 3, also plays an important role in the mouse heart development and function [51].

A number of other unrecognized genes identified in the genome-wide analysis are also potentially implicated in the pathogenesis of AVS, however the association between these genes and heart defects and putative mechanisms of action remain to be further explored.

In summary, we have identified alterations in DNA methylation in a number of genes involved in biological processes important to cardiac development in newborns with AVS, supporting the hypothesis that epigenetic modification may play a crucial role in the development of congenital heart defects such as AVS.

In the present study we have identified novel DNA-methylated genes and critical biological pathways that correlate with human aortic valve stenosis. We have demonstrated profound methylation differences in multiple CpG sites in these genes in AVS subjects. The methylation levels of individual CpG sites were used to calculate the area under the ROC curves as a measure of the accuracy of a putative diagnostic test with 59 CpG sites. Many of the identified CpG sites have a higher ROC AUC than 0.5; among them, the six CpGs with the highest AUC of more than 0.75 display strong diagnostic potential. Our overall results raise the possibility of using a large number of different marker combinations for effective detection of AVS. The results suggest that DNA methylation analysis is a molecular technique that might have the potential to be developed into a test for screening and diagnosis of AVS.

Finally the present study also determined that among 52 genes identified, many have not previously been reported in association with AVS. These results also suggest that the selected genes may be involved in the development of AVS by DNA methylation. The functional role of these genes or their potential as novel DNA methylation markers remains to be examined. An important limitation of the present study is that the findings are based on a relatively small number of subjects; however, these results can form the basis for specific hypotheses regarding AVS pathophysiology. Future larger studies based on this preliminary data may contain further clues for new approaches to non-invasive prenatal diagnosis and even prevention of this serious congenital heart malformation.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Details of the AVS subject cohort and controls used in the present analysis (online only).

https://doi.org/10.1371/journal.pone.0154010.s001

(DOCX)

Acknowledgments

We thank the Applied Genomics Technology Center, Wayne State University School of Medicine, Detroit, Ml for bioinformatics support. We thank the Michigan Department of Community Health for specimens. We also thank Rose Callahan, Beaumont Health, for critical manuscript review and editing.

Author Contributions

Conceived and designed the experiments: UR. Performed the experiments: UR. Analyzed the data: UR AZ. Contributed reagents/materials/analysis tools: SA. Wrote the paper: UR SA ZA-S PS OT. Contributed to design of the study, supervised the clinical, experimental data interpretation and critically revised the manuscript: ROB-S.

Citation: Radhakrishna U, Albayrak S, Alpay-Savasan Z, Zeb A, Turkoglu O, Sobolewski P, et al. (2016) Genome-Wide DNA Methylation Analysis and Epigenetic Variations Associated with Congenital Aortic Valve Stenosis (AVS). PLoS ONE 11(5): e0154010. https://doi.org/10.1371/journal.pone.0154010

References

1. Khoshnood B, Lelong N, Houyel L, Thieulin AC, Jouannic JM, Magnier S, et al. Prevalence, timing of diagnosis and mortality of newborns with congenital heart defects: a population-based study. Heart. 2012;98(22):1667–73. pmid:22888161.

2. Reller MD, Strickland MJ, Riehle-Colarusso T, Mahle WT, Correa A. Prevalence of congenital heart defects in metropolitan Atlanta, 1998–2005. The Journal of pediatrics. 2008;153(6):807–13. pmid:18657826; PubMed Central PMCID: PMC2613036.

3. Loffredo CA. Epidemiology of cardiovascular malformations: prevalence and risk factors. American journal of medical genetics. 2000;97(4):319–25. pmid:11376444.

4. Roberts WC, Ko JM. Frequency by decades of unicuspid, bicuspid, and tricuspid aortic valves in adults having isolated aortic valve replacement for aortic stenosis, with or without associated aortic regurgitation. Circulation. 2005;111(7):920–5. pmid:15710758.

5. Pierpont ME, Basson CT, Benson DW Jr., Gelb BD, Giglia TM, Goldmuntz E, et al. Genetic basis for congenital heart defects: current knowledge: a scientific statement from the American Heart Association Congenital Cardiac Defects Committee, Council on Cardiovascular Disease in the Young: endorsed by the American Academy of Pediatrics. Circulation. 2007;115(23):3015–38. pmid:17519398.

6. Sheffield VC, Pierpont ME, Nishimura D, Beck JS, Burns TL, Berg MA, et al. Identification of a complex congenital heart defect susceptibility locus by using DNA pooling and shared segment analysis. Human molecular genetics. 1997;6(1):117–21. pmid:9002679.

7. McGregor TL, Misri A, Bartlett J, Orabona G, Friedman RD, Sexton D, et al. Consanguinity mapping of congenital heart disease in a South Indian population. PloS one. 2010;5(4):e10286. pmid:20422016; PubMed Central PMCID: PMC2858208.

8. Monteleone PL, Fagan LF. Possible X-linked congenital heart disease. Circulation. 1969;39(5):611–4. pmid:5787314.

9. Hoffman JI, Kaplan S. The incidence of congenital heart disease. Journal of the American College of Cardiology. 2002;39(12):1890–900. pmid:12084585.

10. Patel SS, Burns TL. Nongenetic risk factors and congenital heart defects. Pediatric cardiology. 2013;34(7):1535–55. pmid:23963188.

11. Probst V, Le Scouarnec S, Legendre A, Jousseaume V, Jaafar P, Nguyen JM, et al. Familial aggregation of calcific aortic valve stenosis in the western part of France. Circulation. 2006;113(6):856–60. pmid:16461814.

12. Le Gal G, Bertault V, Bezon E, Cornily JC, Barra JA, Blanc JJ. Heterogeneous geographic distribution of patients with aortic valve stenosis: arguments for new aetiological hypothesis. Heart. 2005;91(2):247–9. pmid:15657257; PubMed Central PMCID: PMC1768671.

13. Friedman WF. Aortic stenosis. In Adama FH, Emmnoulides GC (Eds) Heart diseases in infants, children and adolescents including the fetus and young adults 5th edition Baltimore MD Williams & Wilkins 1995.

14. McBride KL, Marengo L, Canfield M, Langlois P, Fixler D, Belmont JW. Epidemiology of noncomplex left ventricular outflow tract obstruction malformations (aortic valve stenosis, coarctation of the aorta, hypoplastic left heart syndrome) in Texas, 1999–2001. Birth defects research Part A, Clinical and molecular teratology. 2005;73(8):555–61. pmid:16007587; PubMed Central PMCID: PMC1361303.

15. Csaky-Szunyogh M, Vereczkey A, Kosa Z, Gerencser B, Czeizel AE. Risk factors in the origin of congenital left-ventricular outflow-tract obstruction defects of the heart: a population-based case-control study. Pediatric cardiology. 2014;35(1):108–20. pmid:23843102.

16. Abuhamad A, Chaoui RC. A practical guide to fetal echocardiography. Normal and abnormal hearts. Second edition 2010: LippincottWilliams and Wilkins Philadelphia, PA. 2010.

17. Kehler L, Biro O, Lazar L, Rigo J Jr., Nagy B. Elevated hsa-miR-99a levels in maternal plasma may indicate congenital heart defects. Biomed Rep. 2015;3(6):869–73. pmid:26623032; PubMed Central PMCID: PMC4660590.

18. Czeizel AE, Dobo M, Vargha P. Hungarian cohort-controlled trial of periconceptional multivitamin supplementation shows a reduction in certain congenital abnormalities. Birth defects research Part A, Clinical and molecular teratology. 2004;70(11):853–61. pmid:15523663.

19. Marti-Carvajal AJ, Sola I, Lathyris D, Salanti G. Homocysteine lowering interventions for preventing cardiovascular events. The Cochrane database of systematic reviews. 2009;(4):CD006612. pmid:19821378.

20. Bahado-Singh RO, Zaffra R, Albayarak S, Chelliah A, Bolinjkar R, Turkoglu O, et al. Epigenetic markers for newborn congenital heart defect (CHD). J Matern Fetal Neonatal Med. 2015:1–7.

21. Ailes EC, Gilboa SM, Riehle-Colarusso T, Johnson CY, Hobbs CA, Correa A, et al. Prenatal diagnosis of nonsyndromic congenital heart defects. Prenatal diagnosis. 2014;34(3):214–22. pmid:24222433.

22. Knowles R, Griebsch I, Dezateux C, Brown J, Bull C, Wren C. Newborn screening for congenital heart defects: a systematic review and cost-effectiveness analysis. Health technology assessment. 2005;9(44):1–152, iii-iv. pmid:16297355.

23. Joo JE, Wong EM, Baglietto L, Jung CH, Tsimiklis H, Park DJ, et al. The use of DNA from archival dried blood spots with the Infinium HumanMethylation450 array. BMC Biotechnol. 2013;13:23. pmid:23497093; PubMed Central PMCID: PMC3610215.

24. Wong N, Morley R, Saffery R, Craig J. Archived Guthrie blood spots as a novel source for quantitative DNA methylation analysis. Biotechniques. 2008;45(4):423–4, 6, 8 passim. pmid:18855769.

25. Docherty LE, Rezwan FI, Poole RL, Jagoe H, Lake H, Lockett GA, et al. Genome-wide DNA methylation analysis of patients with imprinting disorders identifies differentially methylated regions associated with novel candidate imprinted genes. Journal of medical genetics. 2014;51(4):229–38. pmid:24501229; PubMed Central PMCID: PMC3963529.

26. Pan H, Chen L, Dogra S, Teh AL, Tan JH, Lim YI, et al. Measuring the methylome in clinical samples: improved processing of the Infinium Human Methylation450 BeadChip Array. Epigenetics: official journal of the DNA Methylation Society. 2012;7(10):1173–87. pmid:22964528; PubMed Central PMCID: PMC3469459.

27. Liu Y, Aryee MJ, Padyukov L, Fallin MD, Hesselberg E, Runarsson A, et al. Epigenome-wide association data implicate DNA methylation as an intermediary of genetic risk in rheumatoid arthritis. Nat Biotechnol. 2013;31(2):142–7. pmid:23334450; PubMed Central PMCID: PMCPMC3598632.

28. Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9. pmid:23314698; PubMed Central PMCID: PMCPMC3592906.

29. Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, et al. Review of processing and analysis methods for DNA methylation array data. Br J Cancer. 2013;109(6):1394–402. pmid:23982603; PubMed Central PMCID: PMC3777004.

30. Daca-Roszak P, Pfeifer A, Zebracka-Gala J, Rusinek D, Szybinska A, Jarzab B, et al. Impact of SNPs on methylation readouts by Illumina Infinium HumanMethylation450 BeadChip Array: implications for comparative population studies. BMC Genomics. 2015;16(1):1003. pmid:26607064; PubMed Central PMCID: PMC4659175.

31. Gu. Z. ComplexHeatmap: Making Complex Heatmaps. R package version 1.6.0. https://githubcom/jokergoo/ComplexHeatmap. 2015.

32. Altorok N, Tsou PS, Coit P, Khanna D, Sawalha AH. Genome-wide DNA methylation analysis in dermal fibroblasts from patients with diffuse and limited systemic sclerosis reveals common and subset-specific DNA methylation aberrancies. Annals of the rheumatic diseases. 2014. pmid:24812288.

33. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Statist Soc B. 1995;57:289–300.

34. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nature protocols. 2009;4(1):44–57. pmid:19131956.

35. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic acids research. 2009;37(1):1–13. pmid:19033363; PubMed Central PMCID: PMC2615629.

36. Sheng W, Qian Y, Zhang P, Wu Y, Wang H, Ma X, et al. Association of promoter methylation statuses of congenital heart defect candidate genes with Tetralogy of Fallot. Journal of translational medicine. 2014;12:31. pmid:24479926; PubMed Central PMCID: PMC3915623.

37. Hutter CM, Austin MA, Humphries SE. Familial hypercholesterolemia, peripheral arterial disease, and stroke: a HuGE minireview. Am J Epidemiol. 2004;160(5):430–5. pmid:15321839.

38. Austin MA, Hutter CM, Zimmern RL, Humphries SE. Familial hypercholesterolemia and coronary heart disease: a HuGE association review. Am J Epidemiol. 2004;160(5):421–9. pmid:15321838.

39. Rabus MB, Kayalar N, Sareyyupoglu B, Erkin A, Kirali K, Yakut C. Hypercholesterolemia association with aortic stenosis of various etiologies. J Card Surg. 2009;24(2):146–50. pmid:19267822.

40. Chui MC, Newby DE, Panarelli M, Bloomfield P, Boon NA. Association between calcific aortic stenosis and hypercholesterolemia: is there a need for a randomized controlled trial of cholesterol-lowering therapy? Clin Cardiol. 2001;24(1):52–5. pmid:11195607.

41. Li Y. Functional analysis of Dusp27, a novel target gene of the JAK1/STAT1 pathway, in myogenesis. Thesis (PhD)—Hong Kong University of Science and Technology,. 2011.

42. Arrington CB, Bleyl SB, Matsunami N, Bonnell GD, Otterud BE, Nielsen DC, et al. Exome analysis of a family with pleiotropic congenital heart disease. Circ Cardiovasc Genet. 2012;5(2):175–82. pmid:22337856; PubMed Central PMCID: PMC3329568.

43. Arsenault BJ, Dube MP, Brodeur MR, de Oliveira Moraes AB, Lavoie V, Kernaleguen AE, et al. Evaluation of links between high-density lipoprotein genetics, functionality, and aortic valve stenosis risk in humans. Arterioscler Thromb Vasc Biol. 2014;34(2):457–62. pmid:24334872.

44. Liu HK, Wang CT, Zhang SZ, Xiao CY, Li XF, Zhang KL, et al. [Association of APOA5 gene single nucleotide polymorphism with levels of lipids and coronary heart disease in Chinese]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2004;21(4):335–8. pmid:15300628.

45. Wang J, Ban MR, Kennedy BA, Anand S, Yusuf S, Huff MW, et al. APOA5 genetic variants are markers for classic hyperlipoproteinemia phenotypes and hypertriglyceridemia. Nat Clin Pract Cardiovasc Med. 2008;5(11):730–7. pmid:18779834.

46. Cohen JC, Boerwinkle E, Mosley TH Jr., Hobbs HH. Sequence variations in PCSK9, low LDL, and protection against coronary heart disease. N Engl J Med. 2006;354(12):1264–72. pmid:16554528.

47. Click ES, Cox B, Olson SB, Grompe M, Akkari Y, Moreau LA, et al. Fanconi anemia-like presentation in an infant with constitutional deletion of 21q including the RUNX1 gene. Am J Med Genet A. 2011;155A(7):1673–9. pmid:21626672.

48. Lindstrand A, Malmgren H, Sahlen S, Schoumans J, Nordgren A, Ergander U, et al. Detailed molecular and clinical characterization of three patients with 21q deletions. Clin Genet. 2010;77(2):145–54. pmid:19863549.

49. Conrad M, Jakupoglu C, Moreno SG, Lippl S, Banjac A, Schneider M, et al. Essential role for mitochondrial thioredoxin reductase in hematopoiesis, heart development, and heart function. Mol Cell Biol. 2004;24(21):9414–23. pmid:15485910; PubMed Central PMCID: PMC522221.

50. Lindsay EA. Chromosomal microdeletions: dissecting del22q11 syndrome. Nat Rev Genet. 2001;2(11):858–68. pmid:11715041.

51. Helmer RA, Martinez-Zaguilan R, Dertien JS, Fulford C, Foreman O, Peiris V, et al. Helicase-like transcription factor (Hltf) regulates G2/M transition, Wt1/Gata4/Hif-1a cardiac transcription networks, and collagen biogenesis. PLoS One. 2013;8(11):e80461. pmid:24278285; PubMed Central PMCID: PMC3835564.