Innovative approach to identify multigenomic and environmental interactions associated with birth defects in family-based hybrid designs.

Genes, including those with transgenerational effects, work in concert with behavioral, environmental, and social factors via complex biological networks to determine human health. Understanding complex relationships between causal factors underlying human health is an essential step towards deciphering biological mechanisms. We propose a new analytical framework to investigate the interactions between maternal and offspring genetic variants or their surrogate single nucleotide polymorphisms (SNPs) and environmental factors using family-based hybrid study design. The proposed approach can analyze diverse genetic and environmental factors and accommodate samples from a variety of family units, including case/control-parental triads, and case/control-parental dyads, while minimizing potential bias introduced by population admixture. Comprehensive simulations demonstrated that our innovative approach outperformed the log-linear approach, the best available method for case-control family data. The proposed approach had greater statistical power and was capable to unbiasedly estimate the maternal and child genetic effects and the effects of environmental factors, while controlling the Type I error rate against population stratification. Using our newly developed approach, we analyzed the associations between maternal and fetal SNPs and obstructive and conotruncal heart defects, with adjustment for demographic and lifestyle factors and dietary supplements. Fourteen and 11 fetal SNPs were associated with obstructive and conotruncal heart defects, respectively. Twenty-seven and 17 maternal SNPs were associated with obstructive and conotruncal heart defects, respectively. In addition, maternal body mass index was a significant risk factor for obstructive defects. The proposed approach is a powerful tool for interrogating the etiological mechanism underlying complex traits.

A parent-of-origin analysis of paternal genetic variants and increased risk of conotruncal heart defects.

The association between conotruncal heart defects (CTHDs) and maternal genetic and environmental exposures is well studied. However, little is known about paternal genetic or environmental exposures and risk of CTHDs. We assessed the effect of paternal genetic variants in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs in offspring. We utilized National Birth Defects Prevention Study data to conduct a family-based case only study using 616 live-born infants with CTHDs, born October 1997-August 2008. Maternal, paternal and infant DNA was genotyped using an Illumina® Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR) and 95% confidence intervals (CI) from log-linear models determined parent of origin effects for 921 SNPs in 60 candidate genes involved in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs. The risk of CTHD among children who inherited a paternally derived copy of the A allele on GLRX (rs17085159) or the T allele of GLRX (rs12109442) was 0.23 (95%CI: 0.12, 0.42; p = 1.09 × 10-6 ) and 0.27 (95%CI: 0.14, 0.50; p = 2.06 × 10-5 ) times the risk among children who inherited a maternal copy of the same allele. The paternally inherited copy of the GSR (rs7818511) A allele had a 0.31 (95%CI: 0.18, 0.53; p = 9.94 × 10-6 ] risk of CTHD compared to children with the maternal copy of the same allele. The risk of CTHD is less influenced by variants in paternal genes involved in the folate, homocysteine, or transsulfuration pathways than variants in maternal genes in those pathways.

Testing Allele Transmission of an SNP Set Using a Family-Based Generalized Genetic Random Field Method.

Family-based association studies are commonly used in genetic research because they can be robust to population stratification (PS). Recent advances in high-throughput genotyping technologies have produced a massive amount of genomic data in family-based studies. However, current family-based association tests are mainly focused on evaluating individual variants one at a time. In this article, we introduce a family-based generalized genetic random field (FB-GGRF) method to test the joint association between a set of autosomal SNPs (i.e., single-nucleotide polymorphisms) and disease phenotypes. The proposed method is a natural extension of a recently developed GGRF method for population-based case-control studies. It models offspring genotypes conditional on parental genotypes, and, thus, is robust to PS. Through simulations, we presented that under various disease scenarios the FB-GGRF has improved power over a commonly used family-based sequence kernel association test (FB-SKAT). Further, similar to GGRF, the proposed FB-GGRF method is asymptotically well-behaved, and does not require empirical adjustment of the type I error rates. We illustrate the proposed method using a study of congenital heart defects with family trios from the National Birth Defects Prevention Study (NBDPS).

A Three-Way Interaction among Maternal and Fetal Variants Contributing to Congenital Heart Defects.

Congenital heart defects (CHDs) develop through a complex interplay between genetic variants, epigenetic modifications, and maternal environmental exposures. Genetic studies of CHDs have commonly tested single genetic variants for association with CHDs. Less attention has been given to complex gene-by-gene and gene-by-environment interactions. In this study, we applied a recently developed likelihood-ratio Mann-Whitney (LRMW) method to detect joint actions among maternal variants, fetal variants, and maternal environmental exposures, allowing for high-order statistical interactions. All subjects are participants from the National Birth Defect Prevention Study, including 623 mother-offspring pairs with CHD-affected pregnancies and 875 mother-offspring pairs with unaffected pregnancies. Each individual has 872 single nucleotide polymorphisms encoding for critical enzymes in the homocysteine, folate, and trans-sulfuration pathways. By using the LRMW method, three variants (fetal rs625879, maternal rs2169650, and maternal rs8177441) were identified with a joint association to CHD risk (nominal P-value = 1.13e-07). These three variants are located within genes BHMT2, GSTP1, and GPX3, respectively. Further examination indicated that maternal SNP rs2169650 may interact with both fetal SNP rs625879 and maternal SNP rs8177441. Our findings suggest that the risk of CHD may be influenced by both the intragenerational interaction within the maternal genome and the intergenerational interaction between maternal and fetal genomes.

Krüppel-like factor 9 (KLF9) prevents colorectal cancer through inhibition of interferon-related signaling.

UNLABELLED: Expression of the transcription factor Krüppel-like factor 9 (KLF9) is frequently reduced in colorectal cancers, although a tumor suppressive role has not been established. To determine if KLF9 suppresses intestinal adenoma formation, we generated mice of distinct Klf9 genotypes in the background of the Apc (Min/+) mouse and compared their adenoma burdens at 16 weeks of age. While small intestine adenoma burden remained unchanged among Klf9 genotypes, male and female Apc(Min/+)/Klf9(-/-) and Apc(Min/+)/Klf9(+/-) mice exhibited significantly more colon adenomas than their Apc(Min/+)/Klf9(+/+) counterparts. Microarray analysis showed significant increases in the expression of interferon-induced genes in the colon mucosa of female Apc (Min/+)/Klf9(+/-) and Apc(Min/+)/Klf9(-/-) compared to Apc(Min/+)/Klf9(+/+) mice, prior to overt adenoma occurrence. Gene upregulation was confirmed by qPCR of colon mucosa and by siRNA knockdown of KLF9 in human HT29 colorectal cancer cells. Increases in expression of these genes were further augmented by supplementation with Interferon ß1. Circulating levels of the cytokine, interferon-stimulated gene 15 (ISG15) were increased in Apc(Min/+)/Klf9(+/-) and Apc(Min/+)/Klf9(-/-) mice relative to Apc(Min/+)/Klf9(+/+). Additionally, colon mucosal levels of ISG15 were increased in Apc(Min/+)/Klf9(+/-) mice. Chromatin immunoprecipitation demonstrated KLF9 recruitment to the ISG15 promoter. Lastly, treatment with ISG15 suppressed apoptosis in HT29 cells, in the presence and absence of 5-fluorouracil (5FU). Results show KLF9 to be a haploinsufficient suppressor of colon tumorigenesis in Apc(Min/+) mice in part, by repression of ISG15 and the latter's antiapoptotic function. SUMMARY: Krüppel-like factor 9 (KLF9) is a haploinsufficient tumor suppressor in the ApcMin/+ mouse colon by suppressing expression of ISG15, an apoptosis-inhibiting cytokine.

Detecting maternal-fetal genotype interactions associated with conotruncal heart defects: a haplotype-based analysis with penalized logistic regression.

Nonsyndromic congenital heart defects (CHDs) develop during embryogenesis as a result of a complex interplay between environmental exposures, genetics, and epigenetic causes. Genetic factors associated with CHDs may be attributed to either independent effects of maternal or fetal genes, or the intergenerational interactions between maternal and fetal genes. Detecting gene-by-gene interactions underlying complex diseases is a major challenge in genetic research. Detecting maternal-fetal genotype (MFG) interactions and differentiating them from the maternal/fetal main effects has presented additional statistical challenges due to correlations between maternal and fetal genomes. Traditionally, genetic variants are tested separately for maternal/fetal main effects and MFG interactions on a single-locus basis. We conducted a haplotype-based analysis with a penalized logistic regression framework to dissect the genetic effect associated with the development of nonsyndromic conotruncal heart defects (CTD). Our method allows simultaneous model selection and effect estimation, providing a unified framework to differentiate maternal/fetal main effect from the MFG interaction effect. In addition, the method is able to test multiple highly linked SNPs simultaneously with a configuration of haplotypes, which reduces the data dimensionality and the burden of multiple testing. By analyzing a dataset from the National Birth Defects Prevention Study (NBDPS), we identified seven genes (GSTA1, SOD2, MTRR, AHCYL2, GCLC, GSTM3, and RFC1) associated with the development of CTDs. Our findings suggest that MFG interactions between haplotypes in three of seven genes, GCLC, GSTM3, and RFC1, are associated with nonsyndromic conotruncal heart defects.

Obstructive heart defects associated with candidate genes, maternal obesity, and folic acid supplementation.

Right-sided and left-sided obstructive heart defects (OHDs) are subtypes of congenital heart defects, in which the heart valves, arteries, or veins are abnormally narrow or blocked. Previous studies have suggested that the development of OHDs involved a complex interplay between genetic variants and maternal factors. Using the data from 569 OHD case families and 1,644 control families enrolled in the National Birth Defects Prevention Study (NBDPS) between 1997 and 2008, we conducted an analysis to investigate the genetic effects of 877 single nucleotide polymorphisms (SNPs) in 60 candidate genes for association with the risk of OHDs, and their interactions with maternal use of folic acid supplements, and pre-pregnancy obesity. Applying log-linear models based on the hybrid design, we identified a SNP in methylenetetrahydrofolate reductase (MTHFR) gene (C677T polymorphism) with a main genetic effect on the occurrence of OHDs. In addition, multiple SNPs in betaine-homocysteine methyltransferase (BHMT and BHMT2) were also identified to be associated with the occurrence of OHDs through significant main infant genetic effects and interaction effects with maternal use of folic acid supplements. We also identified multiple SNPs in glutamate-cysteine ligase, catalytic subunit (GCLC) and DNA (cytosine-5-)-methyltransferase 3 beta (DNMT3B) that were associated with elevated risk of OHDs among obese women. Our findings suggested that the risk of OHDs was closely related to a combined effect of variations in genes in the folate, homocysteine, or glutathione/transsulfuration pathways, maternal use of folic acid supplements and pre-pregnancy obesity.

Genetic variation affects congenital heart defect susceptibility in offspring exposed to maternal tobacco use.

BACKGROUND: Congenital heart defects (CHDs) are among the most prevalent and serious birth defects, occurring in 8 to 10 of every 1000 live births in the United States. Epidemiologic studies have reported an association between CHDs and maternal smoking, but it remains unknown how genes impact the susceptibility of offspring to CHDs in the presence of maternal tobacco use. METHODS: Using data from 403 case- and 219 control-parental triads enrolled in the National Birth Defects Prevention Study between 1998 and 2008, we investigated the association between CHDs and maternal and infant genetic variants involved in the tobacco metabolism and DNA repair pathways among mothers who smoked prenatally. RESULTS: The maternal genotypes of single nucleotide polymorphisms in the excision repair cross-complementation group 1 (ERCC1), poly (ADP-ribose) polymerase 2 (PARP2), and ERCC5 genes were identified to be significantly associated with the occurrence of CHDs in the presence of maternal tobacco use. Our analysis also revealed a moderate association between the infant genotypes of polymorphisms in the O-sialoglycoprotein endopeptidase (OSGEP) gene and increased risk of CHDs among mothers who smoked. CONCLUSION: Our study provides evidence that maternal and infant polymorphisms within the ERCC1, PARP2, ERCC5, and OSGEP genes are associated with CHD risk in the presence of maternal tobacco use. These results may provide insight into the susceptibility of having a pregnancy affected by CHDs among women who smoke.

Developments in our understanding of the genetic basis of birth defects.

Birth defects are a major cause of morbidity and mortality worldwide. There has been much progress in understanding the genetic basis of familial and syndromic forms of birth defects. However, the etiology of nonsydromic birth defects is not well-understood. Although there is still much work to be done, we have many of the tools needed to accomplish the task. Advances in next-generation sequencing have introduced a sea of possibilities, from disease-gene discovery to clinical screening and diagnosis. These advances have been fruitful in identifying a host of candidate disease genes, spanning the spectrum of birth defects. With the advent of CRISPR-Cas9 gene editing, researchers now have a precise tool for characterizing this genetic variation in model systems. Work in model organisms has also illustrated the importance of epigenetics in human development and birth defects etiology. Here we review past and current knowledge in birth defects genetics. We describe genotyping and sequencing methods for the detection and analysis of rare and common variants. We remark on the utility of model organisms and explore epigenetics in the context of structural malformation. We conclude by highlighting approaches that may provide insight into the complex genetics of birth defects.

A genetic association study detects haplotypes associated with obstructive heart defects.

The development of congenital heart defects (CHDs) involves a complex interplay between genetic variants, epigenetic variants, and environmental exposures. Previous studies have suggested that susceptibility to CHDs is associated with maternal genotypes, fetal genotypes, and maternal-fetal genotype (MFG) interactions. We conducted a haplotype-based genetic association study of obstructive heart defects (OHDs), aiming to detect the genetic effects of 877 SNPs involved in the homocysteine, folate, and transsulfuration pathways. Genotypes were available for 285 mother-offspring pairs with OHD-affected pregnancies and 868 mother-offspring pairs with unaffected pregnancies. A penalized logistic regression model was applied with an adaptive least absolute shrinkage and selection operator (lasso), which dissects the maternal effect, fetal effect, and MFG interaction effects associated with OHDs. By examining the association between 140 haplotype blocks, we identified 9 blocks that are potentially associated with OHD occurrence. Four haplotype blocks, located in genes MGMT, MTHFS, CBS, and DNMT3L, were statistically significant using a Bayesian false-discovery probability threshold of 0.8. Two blocks in MGMT and MTHFS appear to have significant fetal effects, while the CBS and DNMT3L genes may have significant MFG interaction effects.

Conotruncal heart defects and common variants in maternal and fetal genes in folate, homocysteine, and transsulfuration pathways.

BACKGROUND: We investigated the association between conotruncal heart defects (CTDs) and maternal and fetal single nucleotide polymorphisms (SNPs) in 60 genes in the folate, homocysteine, and transsulfuration pathways. We also investigated whether periconceptional maternal folic acid supplementation modified associations between CTDs and SNPs METHODS: Participants were enrolled in the National Birth Defects Prevention Study between 1997 and 2008. DNA samples from 616 case-parental triads affected by CTDs and 1645 control-parental triads were genotyped using an Illumina® Golden Gate custom SNP panel. A hybrid design analysis, optimizing data from case and control trios, was used to identify maternal and fetal SNPs associated with CTDs RESULTS: Among 921 SNPs, 17 maternal and 17 fetal SNPs had a Bayesian false-discovery probability of <0.8. Ten of the 17 maternal SNPs and 2 of the 17 fetal SNPs were found within the glutamate-cysteine ligase, catalytic subunit (GCLC) gene. Fetal SNPs with the lowest Bayesian false-discovery probability (rs2612101, rs2847607, rs2847326, rs2847324) were found within the thymidylate synthetase (TYMS) gene. Additional analyses indicated that the risk of CTDs associated with candidate SNPs was modified by periconceptional folic acid supplementation. Nineteen maternal and nine fetal SNPs had a Bayesian false-discovery probability <0.8 for gene-by-environment (G × E) interactions with maternal folic acid supplementation. CONCLUSION: These results support previous studies suggesting that maternal and fetal SNPs within folate, homocysteine, and transsulfuration pathways are associated with CTD risk. Maternal use of supplements containing folic acid may modify the impact of SNPs on the developing heart.

Gene-Folic Acid Interactions and Risk of Conotruncal Heart Defects: Results from the National Birth Defects Prevention Study.

Conotruncal heart defects (CTDs) are heart malformations that affect the cardiac outflow tract and typically cause significant morbidity and mortality. Evidence from epidemiological studies suggests that maternal folate intake is associated with a reduced risk of heart defects, including CTD. However, it is unclear if folate-related gene variants and maternal folate intake have an interactive effect on the risk of CTDs. In this study, we performed targeted sequencing of folate-related genes on DNA from 436 case families with CTDs who are enrolled in the National Birth Defects Prevention Study and then tested for common and rare variants associated with CTD. We identified risk alleles in maternal MTHFS (ORmeta = 1.34; 95% CI 1.07 to 1.67), maternal NOS2 (ORmeta = 1.34; 95% CI 1.05 to 1.72), fetal MTHFS (ORmeta = 1.35; 95% CI 1.09 to 1.66), and fetal TCN2 (ORmeta = 1.38; 95% CI 1.12 to 1.70) that are associated with an increased risk of CTD among cases without folic acid supplementation. We detected putative de novo mutations in genes from the folate, homocysteine, and transsulfuration pathways and identified a significant association between rare variants in MGST1 and CTD risk. Results suggest that periconceptional folic acid supplementation is associated with decreased risk of CTD among individuals with susceptible genotypes.

Associations between maternal genotypes and metabolites implicated in congenital heart defects.

BACKGROUND: The development of non-syndromic congenital heart defects (CHDs) involves a complex interplay of genetics, metabolism, and lifestyle. Previous studies have implicated maternal single nucleotide polymorphisms (SNPs) and altered metabolism in folate-related pathways as CHD risk factors. OBJECTIVE: We sought to discover associations between maternal SNPs and metabolites involved in the homocysteine, folate, and transsulfuration pathways, and determine if these associations differ between CHD cases and controls. DESIGN: Genetic, metabolic, demographic, and lifestyle information was available for 335 mothers with CHD-affected pregnancies and 263 mothers with unaffected pregnancies. Analysis was conducted on 1160 SNPs, 13 plasma metabolites, and 2 metabolite ratios. A two-stage multiple linear regression was fitted to each combination of SNP and metabolite/ratio. RESULTS: We identified 4 SNPs in the methionine adenosyltransferase II alpha (MAT2A) gene that were associated with methionine levels. Three SNPs in tRNA aspartic acid methyltransferase 1 (TRDMT1) gene were associated with total plasma folate levels. Glutamylcysteine (GluCys) levels were associated with multiple SNPs within the glutathione peroxidase 6 (GPX6) and O-6-methylguanine-DNA methyltransferase (MGMT) genes. The regression model revealed interactions between genotype and case-control status in the association of total plasma folate, total glutathione (GSH), and free GSH, to SNPs within the MGMT, 5,10-methenyltetrahydrofolate synthetase (MTHFS), and catalase (CAT) genes, respectively. CONCLUSIONS: Our study provides further evidence that genetic variation within folate-related pathways accounts for inter-individual variability in key metabolites. We identified specific SNP-metabolite relationships that differed in mothers with CHD-affected pregnancies, compared to controls. Our results underscore the importance of multifactorial studies to define maternal CHD risk.

Cheek swabs, SNP chips, and CNVs: assessing the quality of copy number variant calls generated with subject-collected mail-in buccal brush DNA samples on a high-density genotyping microarray.

BACKGROUND: Multiple investigators have established the feasibility of using buccal brush samples to genotype single nucleotide polymorphisms (SNPs) with high-density genome-wide microarrays, but there is currently no consensus on the accuracy of copy number variants (CNVs) inferred from these data. Regardless of the source of DNA, it is more difficult to detect CNVs than to genotype SNPs using these microarrays, and it therefore remains an open question whether buccal brush samples provide enough high-quality DNA for this purpose. METHODS: To demonstrate the quality of CNV calls generated from DNA extracted from buccal samples, compared to calls generated from blood samples, we evaluated the concordance of calls from individuals who provided both sample types. The Illumina Human660W-Quad BeadChip was used to determine SNPs and CNVs of 39 Arkansas participants in the National Birth Defects Prevention Study (NBDPS), including 16 mother-infant dyads, who provided both whole blood and buccal brush DNA samples. RESULTS: We observed a 99.9% concordance rate of SNP calls in the 39 blood-buccal pairs. From the same dataset, we performed a similar analysis of CNVs. Each of the 78 samples was independently segmented into regions of like copy number using the Optimal Segmentation algorithm of Golden Helix SNP & Variation Suite 7.Across 640,663 loci on 22 autosomal chromosomes, segment-mean log R ratios had an average correlation of 0.899 between blood-buccal pairs of samples from the same individual, while the average correlation between all possible blood-buccal pairs of samples from unrelated individuals was 0.318. An independent analysis using the QuantiSNP algorithm produced average correlations of 0.943 between blood-buccal pairs from the same individual versus 0.332 between samples from unrelated individuals.Segment-mean log R ratios had an average correlation of 0.539 between mother-offspring dyads of buccal samples, which was not statistically significantly different than the average correlation of 0.526 between mother-offspring dyads of blood samples (p=0.302). CONCLUSIONS: We observed performance from the subject-collected mail-in buccal brush samples comparable to that of blood. These results show that such DNA samples can be used for genome-wide scans of both SNPs and CNVs, and that high rates of CNV concordance were achieved whether using a change-point-based algorithm or one based on a hidden Markov model (HMM).

Congenital heart defects and maternal genetic, metabolic, and lifestyle factors.

BACKGROUND: The purpose of this study was to identify metabolic, genetic, and lifestyle factors that discriminate between women who have pregnancies affected by congenital heart defects (CHDs) from those who have unaffected pregnancies. METHODS: We analyzed the concentrations of 13 serum biomarkers, 3 functional genetic variants, and 4 lifestyle factors among 417 women with CHD-affected pregnancies and 250 controls. To identify risk factors that discriminated between cases and controls, we used logistic regression followed by recursive partitioning to identify non-linear interactions. A receiver operating characteristic (ROC) curve was constructed to evaluate the discriminatory accuracy of the final model. RESULTS: A combination of risk factors discriminated women who had pregnancies affected by CHDs from those who had unaffected pregnancies. Among 21 possible determinants, serum concentrations of homocysteine and methionine, and reduced:oxidized glutathione ratios (GSH:GSSG) had the greatest discriminatory power. Recursive partition modeling resulted in five terminal nodes each illustrating the interplay of these three biomarkers. Women with elevated homocysteine and low GSH:GSSG had the highest risk of having CHD-affected pregnancy, whereas women with low homocysteine, high methionine, and high GSH:GSSG had the lowest risk. The corresponding area under the ROC curve was 81.6% (95% confidence interval [CI], 78.1-85.2%), indicating high ability to discriminate between cases and controls. CONCLUSION: High homocysteine, low methionine, and a reduced GSH:GSSG ratio were the strongest discriminating factors between cases and controls. Measurement of total homocysteine, methionine, and total and reduced glutathione in reproductive aged women may play a role in primary prevention strategies targeted at CHDs.

Maternal DNA hypomethylation and congenital heart defects.

BACKGROUND: Congenital heart defects (CHDs) are among the most prevalent and serious of birth defects. Multiple maternal factors are thought to contribute to CHD development, including folate intake. Maternal DNA methylation, which is dependent on folate metabolism, may impact the risk of CHDs. Our study was designed to determine whether maternal long interspersed nucleotide elements-1 (LINE-1) DNA hypomethylation is associated with increased occurrence of non-syndromic CHDs and whether maternal folate-dependent metabolites are correlated with DNA methylation status. METHODS: Using a case-control study design, we measured global DNA methylation status among mothers whose pregnancies were affected by non-syndromic CHDs (n = 180) and mothers of unaffected pregnancies (n = 187). Methylation of LINE-1 was used as a surrogate marker of global DNA methylation status. The association between DNA methylation and CHD risk was determined while adjusting for selected lifestyle factors. RESULTS: LINE-1 DNA methylation was significantly lower in cases compared to controls (p = 0.049). After covariate adjustments, a significant difference between cases and controls remained (p = 0.010). Among women with LINE-1 methylation in the lowest decile of DNA methylation, the estimated risk of having a CHD-affected pregnancy was almost twice that of women in all other deciles (odds ratio [OR], 1.91; 95% confidence interval [CI], 1.03-3.58). CONCLUSIONS: Our findings indicate that maternal LINE-1 DNA hypomethylation is associated with an increased risk of CHDs. Future studies investigating the association between maternal DNA methylation patterns and CHDs should be pursued.

Association between selected folate pathway polymorphisms and nonsyndromic limb reduction defects: a case-parental analysis.

Inadequate folate status resulting from either genetic variation or nutritional deficiencies has been associated with an increased risk of congenital malformations including orofacial clefting, limb, cardiac and neural tube defects. Few epidemiological studies have examined the association between limb reduction defects (LRDs) and folate-related genetic polymorphisms other than MTHFR 677C→T. We conducted a case-parental analysis of 148 families who participated in the National Birth Defects Prevention Study to examine the association between nonsyndromic transverse and longitudinal LRDs with five single nucleotide polymorphisms (SNPs) in genes encoding enzymes in folate and methionine pathways. Log-linear Poisson regression, adapted for analysis of case-parental data assuming an additive genetic model, was used to estimate genetic relative risks and 95% confidence intervals for the association between LRDs and each SNP. Among women who did not take multivitamin supplements, the MTHFR 677T variant acts via the offspring's genome to increase the risk of LRDs. No association between LRDs and any fetal SNP was found among women who used multivitamin supplements. These results suggest the possibility that initiating folic acid supplementation prior to pregnancy may reduce the risk of having a LRD-affected pregnancy, especially in women whose offspring inherit one or two copies of the MTHFR 677T variant.

Mapping methylation quantitative trait loci in cardiac tissues nominates risk loci and biological pathways in congenital heart disease.

BACKGROUND: Most congenital heart defects (CHDs) result from complex interactions among genetic susceptibilities, epigenetic modifications, and maternal environmental exposures. Characterizing the complex relationship between genetic, epigenetic, and transcriptomic variation will enhance our understanding of pathogenesis in this important type of congenital disorder. We investigated cis-acting effects of genetic single nucleotide polymorphisms (SNPs) on local DNA methylation patterns within 83 cardiac tissue samples and prioritized their contributions to CHD risk by leveraging results of CHD genome-wide association studies (GWAS) and their effects on cardiac gene expression. RESULTS: We identified 13,901 potential methylation quantitative trait loci (mQTLs) with a false discovery threshold of 5%. Further co-localization analyses and Mendelian randomization indicated that genetic variants near the HLA-DRB6 gene on chromosome 6 may contribute to CHD risk by regulating the methylation status of nearby CpG sites. Additional SNPs in genomic regions on chromosome 10 (TNKS2-AS1 gene) and chromosome 14 (LINC01629 gene) may simultaneously influence epigenetic and transcriptomic variations within cardiac tissues. CONCLUSIONS: Our results support the hypothesis that genetic variants may influence the risk of CHDs through regulating the changes of DNA methylation and gene expression. Our results can serve as an important source of information that can be integrated with other genetic studies of heart diseases, especially CHDs.

Gene-by-gene interactions associated with the risk of conotruncal heart defects.

BACKGROUND: The development of conotruncal heart defects (CTDs) involves a complex relationship among genetic variants and maternal lifestyle factors. In this article, we focused on the interactions between 13 candidate genes within folate, homocysteine, and transsulfuration pathways for potential association with CTD risk. METHODS: Targeted sequencing was used for 328 case-parental triads enrolled in the National Birth Defects Prevention Study (NBDPS). To evaluate the interaction of two genes, we applied a conditional logistic regression model for all possible SNP pairs within two respective genes by contrasting the affected infants with their pseudo-controls. The findings were replicated in an independent sample of 86 NBDPS case-parental triads genotyped by DNA microarrays. The results of two studies were further integrated by a fixed-effect meta-analysis. RESULTS: One SNP pair (i.e., rs4764267 and rs6556883) located in gene MGST1 and GLRX, respectively, was found to be associated with CTD risk after multiple testing adjustment using simpleM, a modified Bonferroni correction approach (nominal p-value of 4.62e-06; adjusted p-value of .04). Another SNP pair (i.e., rs11892646 and rs56219526) located in gene DNMT3A and MTRR, respectively, achieved marginal significance after multiple testing adjustment (adjusted p-value of .06). CONCLUSION: Further studies with larger sample sizes are needed to confirm and elucidate these potential interactions.

Neonatal Diet Impacts Circulatory miRNA Profile in a Porcine Model.

microRNAs (miRNAs) are conserved non-coding small nucleotide molecules found in nearly all species and breastmilk. miRNAs present in breastmilk are very stable to freeze-thaw, RNase treatment, and low pH as they are protected inside exosomes. They are involved in regulating several physiologic and pathologic processes, including immunologic pathways, and we have demonstrated better immune response to vaccines in piglets fed with human milk (HM) in comparison to dairy-based formula (MF). To understand if neonatal diet impacts circulatory miRNA expression, serum miRNA expression was evaluated in piglets fed HM or MF while on their neonatal diet at postnatal day (PND) 21 and post-weaning to solid diet at PND 35 and 51. MF fed piglets showed increased expression of 14 miRNAs and decreased expression of 10 miRNAs, relative to HM fed piglets at PND 21. At PND 35, 9 miRNAs were downregulated in the MF compared to the HM group. At PND 51, 10 miRNAs were decreased and 17 were increased in the MF relative to HM suggesting the persistent effect of neonatal diet. miR-148 and miR-181 were decreased in MF compared to HM at PND 21. Let-7 was decreased at PND 35 while miR-199a and miR-199b were increased at PND 51 in MF compared to HM. Pathway analysis suggested that many of the miRNAs are involved in immune function. In conclusion, we observed differential expression of blood miRNAs at both PND 21 and PND 51. miRNA found in breastmilk were decreased in the serum of the MF group, suggesting that diet impacts circulating miRNA profiles at PND 21. The miRNAs continue to be altered at PND 51 suggesting a persistent effect of the neonatal diet. The sources of miRNAs in circulation need to be evaluated, as the piglets were fed the same solid diet leading up to PND 51 collections. In conclusion, the HM diet appears to have an immediate and persistent effect on the miRNA profile and likely regulates the pathways that impact the immune system and pose benefits to breastfed infants.