Human myelomeningocele risk and ultra-rare deleterious variants in genes associated with cilium, WNT-signaling, ECM, cytoskeleton and cell migration.

Myelomeningocele (MMC) affects one in 1000 newborns annually worldwide and each surviving child faces tremendous lifetime medical and caregiving burdens. Both genetic and environmental factors contribute to disease risk but the mechanism is unclear. This study examined 506 MMC subjects for ultra-rare deleterious variants (URDVs, absent in gnomAD v2.1.1 controls that have Combined Annotation Dependent Depletion score ≥ 20) in candidate genes either known to cause abnormal neural tube closure in animals or previously associated with human MMC in the current study cohort. Approximately 70% of the study subjects carried one to nine URDVs among 302 candidate genes. Half of the study subjects carried heterozygous URDVs in multiple genes involved in the structure and/or function of cilium, cytoskeleton, extracellular matrix, WNT signaling, and/or cell migration. Another 20% of the study subjects carried heterozygous URDVs in candidate genes associated with gene transcription regulation, folate metabolism, or glucose metabolism. Presence of URDVs in the candidate genes involving these biological function groups may elevate the risk of developing myelomeningocele in the study cohort.

Common variants in the Na(+)-coupled bicarbonate transporter genes and salt sensitivity of blood pressure: the GenSalt study.

The current study comprehensively examined the association between common variants in the Na(+)-coupled bicarbonate transporter (NCBT) genes and blood pressure (BP) responses to dietary sodium intervention. A 7-day low-sodium followed by a 7-day high-sodium dietary intervention was conducted among 1906 Han participants from rural areas of northern China. Nine BP measurements were obtained at baseline and each intervention using a random-zero sphygmomanometer. A mixed-effect model was used to assess the additive associations of 76 common variants in five NCBT genes, including SLC4A4, SLC4A5, SLC4A7, SLC4A8 and SLC4A10, with salt sensitivity phenotypes. The Bonferroni method was used to adjust for multiple testing. SLC4A4 marker rs4254735 was significantly associated with diastolic BP (DBP) response to low-sodium intervention (P=5.05 × 10(-4)), with mean (95% confidence interval (CI)) response of -2.91 (-3.21, -2.61) and -0.40 (-1.84, 1.05) mmHg for genotype AA and AG, respectively. In addition, BP responses to high-sodium intervention significantly increased with the number of minor C alleles of SLC4A4 marker rs10022637. Mean systolic BP responses among those with genotypes TT, CT and CC were 4.62 (4.29, 4.99), 5.94 (5.31, 6.58) and 6.00 (3.57, 8.43) mmHg (P=1.14 × 10(-4)); mean DBP responses were 1.72 (1.41, 2.03), 3.22 (2.52, 3.92) and 3.94 (1.88, 5.99) mmHg (P=2.26 × 10(-5)) and mean arterial pressure responses were 2.69 (2.40, 2.97), 4.13 (3.57, 4.70) and 4.61 (2.51, 6.71) mmHg (P=2.07 × 10(-6)), respectively. In brief, the present study indicated that common variants in the SLC4A4 gene might contribute to the variation of BP responses to dietary sodium intake in Han Chinese population.

Associations of epithelial sodium channel genes with blood pressure: the GenSalt study.

In order to investigate the associations of SCNN1A, SCNN1G and SCNN1B genes with blood pressure (BP) in the Han Chinese population, we included 2880 participants did not use antihypertensive medication in the month prior to the baseline survey in the current analysis. Forty-four tag-single-nucleotide polymorphisms (SNPs) in epithelial sodium channel (ENaC) genes were selected and genotyped, and nine BP measurements were obtained during the 3-day examination. In the single-marker analyses, we identified significant associations of SCNN1A marker rs13306613 with diastolic BP (DBP) and SCNN1B marker rs12447134 with systolic BP (SBP) under codominant model after Bonferroni's correction (P=2.82 × 10(-5) and 4.63 × 10(-4), respectively). In addition, five SNPs in SCNN1G and four SNPs in SCNN1B achieved nominal significance for SBP, DBP or mean arterial pressure (MAP) under the additive model. For example, the minor C allele of rs5735 in SCNN1G gene was associated with decreased SBP, DBP and MAP (P=0.016, 5.41 × 10(-3) and 4.36 × 10(-3), respectively). Gene-based results showed significant associations of SCNN1G and SCNN1B with BP levels. This study suggested that ENaC genes have important roles in BP regulation in the Han Chinese population. Future studies are warranted to replicate these findings, and functional studies are needed to identify true causal variants in ENaC genes.

The PPAR alpha gene is associated with triglyceride, low-density cholesterol and inflammation marker response to fenofibrate intervention: the GOLDN study.

As a peroxisome proliferator-activated receptor alpha (PPARα) agonist, fenofibrate favorably modulates dyslipidemia and inflammation markers, which are associated with cardiovascular risk. To determine whether variation in the PPARα receptor gene was associated with lipid and inflammatory marker response, we conducted a 3-week trial of fenofibrate in 861 men and women. Mixed linear models that controlled for age and sex, as well as family pedigree and study center, were constructed using single-nucleotide polymorphisms (SNPs) in the PPARα gene as predictors and changes in fasting triglycerides (TGs), cholesterol and inflammatory markers as outcomes. Significant associations with low-density cholesterol and interleukin-2 (P<0.001) responses to fenofibrate were found. Although there were suggestive associations with tumor necrosis factor-alpha and TG responses (P<0.05), these did not survive the correction for multiple testing. We conclude that variants in the PPARα gene may contribute to future pharmacogenomic paradigms seeking to predict fenofibrate responders from both an anti-dyslipidemic and anti-inflammatory perspective.

Genetic correlation of blood pressure responses to dietary sodium and potassium intervention and cold pressor test in Chinese population.

We examined the genetic association between blood pressure (BP) responses to dietary sodium and potassium intervention and to cold pressor test (CPT) in a large family-based dietary feeding study. The dietary intervention and CPT were conducted among 1906 participants in rural China. The dietary intervention included three 7-day periods of low-sodium feeding (51.3 mmol per day), high-sodium feeding (307.8 mmol per day) and high-sodium feeding plus potassium supplementation (60 mmol per day). BP responses to high-sodium intervention had strong genetic correlations (ρ(G)) with both BP responses to low sodium (ρ(G)=-0.43 to -0.54, P-values=0.0005 to 0.03) and to potassium supplementation (ρ(G)=-0.41 to -0.49, P-values=0.001 to 0.005) interventions. Most environmental correlations between BP responses to various dietary interventions were significant. The ρ(G) between BP responses to CPT and to high-sodium intervention and potassium supplementation were statistically significant. For example, the ρ(G) between maximum BP responses to CPT and BP responses to high-sodium intervention was 0.37 (P=0.006) for systolic BP (SBP) and 0.41 (P=0.002) for diastolic BP (DBP). The ρ(G) between maximum BP responses to CPT and BP responses to potassium intervention was -0.42 (P=0.001) for SBP and -0.46 (P=0.001) for SBP. Our study suggests that there are common genetic determinants that influence BP responses to dietary sodium and potassium interventions and to CPT.

Circulating MCP-1 levels shows linkage to chemokine receptor gene cluster on chromosome 3: the NHLBI family heart study follow-up examination.

Atherogenesis is a chronic inflammatory process. Critical in the inflammation process is monocyte chemoattractant protein-1 (MCP-1). To locate genomic regions that affect circulating MCP-1 levels, a genome-wide linkage scan was conducted in a sample of whites and blacks. Phenotype and genetic marker data were available for 2501 white and 513 black participants in the National Heart Lung Blood Institute Family Heart Study follow-up examination. Heritability for MCP-1 was 0.37 in whites and 0.47 in blacks after adjusting for the effects of sex, age, age-sex interaction, smoking status, lifetime smoking exposure (pack-years) and field center. Significant linkage was observed for MCP-1 in a combined black and white sample on chromosome 3 (logarithm of the odds ratio (LOD)=3.5 at 78 cM, P=0.0001) and suggestive linkage was observed in whites on chromosome 5 (LOD=1.8 at 128 cM, P=0.002). Located under the linkage peak on chromosome 3 is the chemokine receptor gene cluster, including CCR2, the receptor for MCP-1. This study provides preliminary evidence linking genetic variation in a receptor to circulating levels of its ligand, as previously demonstrated for the low-density lipoprotein receptor. Further characterization of these chromosomal regions is needed to identify the functional mutations associated with circulating levels of MCP-1.

QTL-specific genotype-by-smoking interaction and burden of calcified coronary atherosclerosis: the NHLBI Family Heart Study.

BACKGROUND: Calcified coronary plaque (CCP) is a complex trait influenced by both genes and environment, and plausibly an interaction between the two. Because the familial aggregation of CCP has been demonstrated and smoking is a significant, independent predictor of CCP, we assessed the evidence for genotype-by-smoking interaction and conducted linkage analysis of quantitative Agatston CCP scores in participants of the NHLBI Family Heart Study (FHS). METHODS: During standardized clinical exams smoking habits were ascertained and CCP was quantified with cardiac computed tomography (CT). Among 4387 relationship pairs from 2128 Caucasian examinees variance component analysis was implemented in SOLAR to examine: (1) additive genotype-by-smoking status interaction using a variance component approach; (2) linkage analysis in the full sample and among smoking subsets defined by individual smoking exposure; (3) QTL-specific genotype-by-smoking interaction in the regions that appeared to differentiate between smoking strata. RESULTS: The prevalence of CCP (and median Agatston score) was 75% (184.6) in men and 48% (51.0) in women. We detected four genome-wide significant logarithm of odds (LOD) scores in samples stratified by individual smoking exposure: chromosome 4 at 122cM (nearest marker D4S2297; robust adjusted LOD=3.1; q=0.053), chromosome 6 at 99cM (nearest marker D6S1056; robust adjusted LOD=3.3; q=0.053), chromosome 11 at 19cM (nearest marker D11S199; robust adjusted LOD=4.0; q=0.02) and chromosome 13 at 77cM (nearest marker D13S892; robust adjusted LOD=3.1; q=0.053). Additive and QTL-specific genotype-by-smoking interaction was detected on chromosomes 4, 6, 11 and 13; all P<0.05. Three of the four QTLs identified in this report have been previously linked to atherosclerosis and harbor interesting candidate genes. CONCLUSIONS: These findings demonstrate the importance of considering complex interactions in the search for genes that influence the pathogenesis of CCP.

The contribution of individual and pairwise combinations of SNPs in the APOA1 and APOC3 genes to interindividual HDL-C variability.

Apolipoproteins (apo) A-I and C-III are components of high-density lipoprotein-cholesterol (HDL-C), a quantitative trait negatively correlated with risk of cardiovascular disease (CVD). We analyzed the contribution of individual and pairwise combinations of single nucleotide polymorphisms (SNPs) in the APOA1/APOC3 genes to HDL-C variability to evaluate (1) consistency of published single-SNP studies with our single-SNP analyses; (2) consistency of single-SNP and two-SNP phenotype-genotype relationships across race-, gender-, and geographical location-dependent contexts; and (3) the contribution of single SNPs and pairs of SNPs to variability beyond that explained by plasma apo A-I concentration. We analyzed 45 SNPs in 3,831 young African-American (N=1,858) and European-American (N=1,973) females and males ascertained by the Coronary Artery Risk Development in Young Adults (CARDIA) study. We found three SNPs that significantly impact HDL-C variability in both the literature and the CARDIA sample. Single-SNP analyses identified only one of five significant HDL-C SNP genotype relationships in the CARDIA study that was consistent across all race-, gender-, and geographical location-dependent contexts. The other four were consistent across geographical locations for a particular race-gender context. The portion of total phenotypic variance explained by single-SNP genotypes and genotypes defined by pairs of SNPs was less than 3%, an amount that is miniscule compared to the contribution explained by variability in plasma apo A-I concentration. Our findings illustrate the impact of context-dependence on SNP selection for prediction of CVD risk factor variability.

Linkage analysis of plasma ApoE in three ethnic groups: multiple genes with context-dependent effects.

We performed variance component-based linkage analysis in four samples (two of non-Hispanic European-Americans from Rochester, MN; African-Americans from Jackson, MS; and Mexican-Americans from Starr County, TX) to identify chromosomal regions containing genes influencing plasma apolipoprotein E (apoE) levels. The APOE gene region on chromosome (chr) 19 was identified with a LOD > or = 2.00 in both samples from Rochester and the sample from Jackson. Adjustment of apoE levels for differences among means of genotypes defined by the APOE epsilon2/3/4 alleles reduced evidence of linkage, indicating that the APOE gene was responsible for the majority of the linkage signal. In stratified linkage analyses, there was a LOD of 1.70 in the Starr County sibships with average total cholesterol (TC) above the median level for all sibships in that population. Adjustment for APOE genotype did not remove this LOD score, suggesting a second gene in this region may influence apoE variation. Evidence of linkage (LOD= 3.32) on chr 17 was observed in the Starr County sibships with average TC below the median. Inter-individual variation in plasma apoE level may be influenced by variations in the structural gene, and at least one other gene whose effects differ among populations and are dependent on the influence of unmeasured genetic and environmental factors indexed by correlated measures of lipid metabolism.

A quantitative trait locus on chromosome 16q influences variation in plasma HDL-C levels in Mexican Americans.

OBJECTIVE: We conducted a whole-genome, multipoint linkage screen to localize a previously reported major locus accounting for 56% to 67% of the additive genetic effects on covariate-adjusted plasma HDL cholesterol (HDL-C) levels in Mexican Americans from the San Antonio Family Heart Study (SAFHS). METHODS AND RESULTS: After using complex segregation analysis to recover the major locus in 472 SAFHS participants from 10 genotyped families, we incorporated covariates required to detect that major locus, including plasma levels of triglycerides and apolipoprotein A-I, in a maximum-likelihood-based variance-components linkage screen. Only chromosome 16 exhibited convincing evidence for a quantitative trait locus (QTL), with a peak multipoint log of the odds (LOD)=3.73 (P=0.000034). Subsequent penetrance model-based linkage analysis, incorporating genotypes at the marker locus nearest the multipoint peak (D16S518) into the segregation model, detected linkage with the previously detected major locus (LOD=2.73, P=0.000642). Initial estimates place this QTL within a 15-cM region of chromosome 16q near the structural loci for lecithin:cholesterol acyltransferase (LCAT) and cholesteryl ester transfer protein (CETP). CONCLUSIONS: A QTL influencing plasma levels of HDL-C in Mexican Americans from San Antonio maps to a region of human chromosome 16q near LCAT and CETP.

Interactions between variants in the beta3-adrenergic receptor and peroxisome proliferator-activated receptor-gamma2 genes and obesity.

OBJECTIVE: Previous studies have reported modest associations between measures of obesity and the Trp64-Arg variant of the beta3-adrenergic receptor (ADRbeta3) and the Pro12Ala variant of the peronisome proliferator-activated receptor (PPAR)-gamma2. We hypothesized that these single gene variants may mark mutations that act through convergent pathways to produce synergistic effects on obesity. RESEARCH DESIGN AND METHODS: The sample included 453 subjects from 10 large Mexican-American families participating in the population-based San Antonio Family Heart Study. The effects of each gene variant singly and jointly were estimated as fixed effects using the measured genotype approach framework. Analyses were conditioned on the pedigree structures to account for the correlations among family members. Statistical significance was evaluated by the likelihood ratio test with adjustment for age, sex and diabetes status. RESULTS: The allele frequencies for the ADRbeta3 Trp64Arg and PPARgamma2 Pro12Ala variants were 18 and 12%, respectively. The ADRbeta3 variant was not significantly associated with any of the obesity-related traits, but subjects with the PPAR-gamma2 variant (n = 98) had significantly higher levels of lasting insulin (P = 0.03), leptin (P = 0.009), and waist circumference (P = 0.03) than those without. Subjects with the gene variants (n = 32) had significantly higher BMI, insulin, and leprtin levels than those with only the PPARgamma2 variant (n = 66) (P for interaction: 0.04, 0.02, and 0.01 for BMI, fasting insulin, and leptin, respectively). CONCLUSIONS: Our results suggest that epistatic models with genes that have modest individual effects may be useful in understanding the genetic underpinnings of typical obesity in humans.

Genome-wide linkage analysis of blood pressure in Mexican Americans.

The genetic mechanisms that control variation in blood pressure level are largely unknown. One of the first steps in understanding those mechanisms is the localization of the genes that have a significant effect on blood pressure. We performed genome scans of systolic (SBP) and diastolic blood pressure (DBP) on a population-based sample of families in the San Antonio Family Heart Study. A likelihood-based Mendelian model incorporating genotype-specific effects of sex, age, age(2), BMI, and blood pressure (SBP or DBP, as appropriate) as covariates was used to perform two-point lodscore (Z) linkage on 399 polymorphic markers. Results showed that the genotype-specific covariate effects were highly significant for both SBP and DBP. Linkage results showed that a quantitative trait locus (QTL) influencing DBP was significantly linked to D2S1790 (Z = 3.92, theta = 0.00) and showed suggestive linkage to D8S373 (Z = 1.92, theta = 0.00). A QTL influencing SBP showed suggestive linkage to D21S1440 (Z = 2.82, theta = 0.00) and D18S844 (Z = 2.09, theta = 0.11). Without the genotype-specific effects in the model, the linkage to D2S1790 was not even suggestive (Z = 1.33, theta = 0.09); thus genotype-specific modeling was crucial in detecting this linkage. A comparison with linkage studies based in other populations showed that the significant linkage to D2S1790 has been replicated at the same marker in the Quebec Family Study. The replicated significant linkage at D2S1790 may begin to establish the locations of the genes that significantly affect blood pressure across several human ethnic groups.

Genome-wide linkage analysis of pulse pressure in Mexican Americans.

Pulse pressure, a measure of aortic stiffness, is a strong predictor of cardiovascular mortality. To locate genes that affect pulse pressure, we performed genetic analysis on randomly ascertained families in the San Antonio Family Heart Study. Pulse pressure was defined as the difference between systolic and diastolic blood pressures. Likelihood methods were used to construct a model that had both single-locus and polygenic components for 46 families (1308 individuals). The single-locus component included sex-specific and genotype-specific effects of both age and body mass index. Using this model, we then performed 2-point linkage analysis in 10 families (440 individuals) that were among the largest of the 46 families and that had been genotyped for 399 polymorphic markers. The model that contained only the polygenic component and simple effects of the covariates showed pulse pressure heritability of 0.21. When the single-locus component was added, the sex-specific and genotype-specific effects of age and body mass index were highly significant (P<0.002). The full model accounted for 73% of the total variation of pulse pressure. Linkage analysis using this model with each marker revealed 4 markers with lod scores >1.9, which is the Lander-Kruglyak suggestive linkage standard. D21S1440 had a lod score of 2.78 with a recombination fraction (theta) of 0.02. D7S1799 had a lod score of 2.04 (theta=0.01), D8S1100 had a lod score of 1.98 (theta=0.08), and D18S844 had a lod score of 1.95 (theta=0.11). These results are highly correlated with results involving systolic blood pressure, indicating that pulse pressure may not be genetically distinct from systolic blood pressure.

A quantitative trait locus influencing estrogen levels maps to a region homologous to human chromosome 20.

Estrogen, a steroid hormone, regulates reproduction and has been implicated in several diseases. We performed a genome-wide scan using multipoint linkage analysis implemented in a general pedigree-based variance component approach to identify genes with measurable effects on variation in estrogen levels in baboons. A microsatellite polymorphism, D20S171, located on human chromosome 20q13.11, showed strong evidence of linkage with a LOD score of 3.06 (P = 0.00009). This region contains several potential candidate genes including melanocortin 3 receptor (MC3R), cytochrome P-450 subfamily XXIV (CYP24), and breast carcinoma amplified sequence (BCAS1). This is the first evidence of a quantitative trait locus with a significant effect on estrogen.

A quantitative trait locus influencing activin-to-estrogen ratio in pedigreed baboons maps to a region homologous to human chromosome 19.

Activin is a multifunctional hormone playing a major role in the regulation of reproduction and growth and development. We performed a genomewide scan using multipoint linkage analysis implemented in a general pedigree-based variance component approach to identify genes with measurable effects on variation in the activin-to-estrogen ratio in baboons. A microsatellite polymorphism, D19S714, which maps to human chromosome 19p13.2, showed marginal evidence of linkage with a lod (log10 of the odds in favor of genetic linkage) score of 1.95 (0.014). This region contains several potential candidate genes including PKA (protein kinase, cAMP-dependent, catalytic alpha) and the gene pair JUN-B and JUN-D. This is the first evidence of a quantitative trait locus with a significant effect on the activin-to-estrogen ratio.

Linkage exclusion analysis of the chromosome 11 region containing UCP2 and UCP3 with obesity-related phenotypes in Mexican Americans.

OBJECTIVE: To investigate whether the region of chromosome 11 (11q13) containing the genes UCP2 and UCP3 could be excluded for linkage with a variety of obesity-related phenotypes in humans. DESIGN: Exclusion mapping using a variance component approach in extended pedigrees. SUBJECTS: Four-hundred and fifty eight individuals (195 females, 263 males) distributed in 10 Mexican American families of probands randomly ascertained with respect to any disease state and who are participating in the San Antonio Family Heart Study. Ages range from 18 to 87 (mean age 35 y). MEASUREMENTS: Serum leptin levels, fat mass (FM), body mass index (BMI), and waist circumference. RESULTS: We were able to exclude the chromosomal region containing UCP2/UCP3 as having an effect on this set of obesity-related phenotypes at relative effect sizes of 10% or greater (P-values<0.05). CONCLUSIONS: These results suggest that variation in these genes is unlikely to have a substantial effect on the expression of obesity-related phenotypes in the Mexican American population.

Genomic searches for genes that influence atherosclerosis and its risk factors.

We are performing genomic searches in randomly ascertained families to identify new quantitative trait loci (QTLs) that influence atherosclerosis and its risk factors. The genetic markers used for genomic searches are random microsatellite markers distributed throughout the human chromosomes. These markers are used for linkage analysis with variance component methods to identify QTLs for measured phenotypes related to lipid metabolism and atherosclerosis. We conducted such a genomic search in 477 participants of the San Antonio Family Heart Study. This genomic search identified QTLs on chromosomes 3 and 4 that influence LDL size class, an important risk factor of atherosclerosis. In addition to lipid risk factors, we measured a variety of gene products involved in atherogenesis in the arterial wall (such as adhesion molecules and components of hemostasis). We found QTLs for serum levels of soluble P-selectin on chromosome 15 (LOD = 3.8) and chromosome 12 (LOD = 2.6).

Linkage of serum insulin concentrations to chromosome 3p in Mexican Americans.

Hyperinsulinemia predicts the development of type 2 diabetes, and family studies suggest that insulin levels are regulated in part by genes. We conducted a genome-wide scan to detect genes influencing variation in fasting serum insulin concentrations in 391 nondiabetic individuals from 10 large multigenerational families. Approximately 380 microsatellite markers with an average spacing of 10 cM were genotyped in all study subjects. Insulin concentrations measured by radioimmunoassay were transformed by their natural logarithms before analysis. In multipoint analysis, peak evidence for linkage occurred on chromosome 3p approximately 109 cM from pter in the region of 3p14.2-p14.1. The multipoint logarithm of odds (LOD) score was 3.07, occurring in the region flanked by markers D3S1600 and D3S1285 (P value by simulation <0.0001). In a two-point analysis, LOD scores ranged from 0.75 to 2.52 for the nine markers typed in the region spanning 88-143 cM from pter. The fasting insulin resistance index was highly correlated with fasting insulin concentrations in this sample and also provided strong evidence for linkage to this region (LOD = 2.99). There was no evidence in our genome-wide scan for linkage of insulin levels to any other chromosome. These results provide evidence that a gene-influencing variation in insulin concentrations exists on chromosome 3p. Possible candidate genes in this region include GBE1 and ACOX2, which encode enzymes involved in glycogen and fatty acid metabolism, respectively.