Incomplete Distal Renal Tubular Acidosis from a Heterozygous Mutation of the V-ATPase B1 Subunit

Congenital distal renal tubular acidosis (dRTA) from mutations of the B1 subunit of the V-ATPase is considered an autosomal recessive disease. We analyzed a dRTA kindred with a truncation-mutation of B1 (p.Phe468fsX487) previously shown to have failure of assembly into the V1 domain of the V-ATPase. All heterozygous carriers in this kindred have normal plasma bicarbonate concentrations, thus evaded the diagnosis of RTA. However, inappropriately high urine pH, hypocitraturia, and hypercalciuria are present either individually or in combination in the heterozygotes at baseline. Two of the heterozygotes studied also have inappropriate urinary acidification with acute ammonium chloride loading and impaired urine-blood pCO2 gradient during bicarbonaturia indicating presence of H+ gradient and flux defects. In normal human renal papillae, wild type B1 is located primarily on the plasma membrane but papilla from one of the heterozygote who had kidney stones had renal tissue secured from surgery showed B1 in both plasma membrane as well as a diffuse intracellular staining. Titrating increasing amounts of the mutant B1 subunit did not exhibit negative dominance over the expression, cellular distribution, or H+-pump activity of the wild type B1 in mammalian HEK293 cells and in V-ATPase-deficient S. cerevisiae. This is the first demonstration of renal acidification defects and nephrolithiasis in heterozygous carriers of mutant B1 subunit; which cannot be attributable to negative dominance. We propose that heterozygosity may lead to mild real acidification defects due to haploinsufficiency. B1 heterozygosity should be considered in patients with calcium nephrolithiasis and urinary abnormalities such as alkalinuria or hypocitraturia.

An ARHGEF10 deletion is highly associated with a juvenile-onset inherited polyneuropathy in Leonberger and Saint Bernard dogs

An inherited polyneuropathy (PN) observed in Leonberger dogs has clinical similarities to a genetically heterogeneous group of peripheral neuropathies termed Charcot-Marie-Tooth (CMT) disease in humans. The Leonberger disorder is a severe, juvenile-onset, chronic, progressive, and mixed PN, characterized by exercise intolerance, gait abnormalities and muscle atrophy of the pelvic limbs, as well as inspiratory stridor and dyspnea. We mapped a PN locus in Leonbergers to a 250 kb region on canine chromosome 16 (Praw = 1.16×10-10, Pgenome, corrected = 0.006) utilizing a high-density SNP array. Within this interval is the ARHGEF10 gene, a member of the rho family of GTPases known to be involved in neuronal growth and axonal migration, and implicated in human hypomyelination. ARHGEF10 sequencing identified a 10 bp deletion in affected dogs that removes four nucleotides from the 3'-end of exon 17 and six nucleotides from the 5'-end of intron 17 (c.1955_1958+6delCACGGTGAGC). This eliminates the 3'-splice junction of exon 17, creates an alternate splice site immediately downstream in which the processed mRNA contains a frame shift, and generates a premature stop codon predicted to truncate approximately 50% of the protein. Homozygosity for the deletion was highly associated with the severe juvenile-onset PN phenotype in both Leonberger and Saint Bernard dogs. The overall clinical picture of PN in these breeds, and the effects of sex and heterozygosity of the ARHGEF10 deletion, are less clear due to the likely presence of other forms of PN with variable ages of onset and severity of clinical signs. This is the first documented severe polyneuropathy associated with a mutation in ARHGEF10 in any species.

Echinococcus P29 antigen: molecular characterization and implication on post-surgery follow-up of CE patients infected with different species of the Echinococcus granulosus complex.

The protein P29 is a potential serological marker for post-treatment monitoring of cystic echinococcosis (CE) especially in young patients. We now have demonstrated that P29 is encoded in the Echinococcus genus by a single gene consisting of 7 exons spanning 1.2 kb of DNA. Variability of the p29 gene at inter- and intra-species level was assessed with 50 cDNA and 280 genomic DNA clones isolated from different E. granulosus s.l. isolates (E. granulosus sensu stricto (G1), E. equinus (G4), E. ortleppi (G5), E. canadensis (G6), E. canadensis (G7) and E. canadensis (G10)) as well as four E. multilocularis isolates. Scarce interspecies polymorphism at the p29 locus was observed and affected predominantly E. granulosus s.s. (G1), where we identified two alleles (A1 and A2) coding for identical P29 proteins and yielding in three genotypes (A1/A1, A2/A2 and A1/A2). Genotypic frequencies expected under Hardy-Weinberg equilibrium revealed a high rate of heterozygosity (47%) that strongly supports the hypothesis that E. granulosus s.s. (G1) is predominantly outbreeding. Comparative sequence analyses of the complete p29 gene showed that phylogenetic relationships within the genus Echinococcus were in agreement with those of previous nuclear gene studies. At the protein level, the deduced P29 amino acid (AA) sequences exhibited a high level of conservation, ranging from 97.9% AA sequence identity among the whole E. granulosus s.l. group to 99.58% identity among E. multilocularis isolates. We showed that P29 proteins of these two species differ by three AA substitutions without implication for antigenicity. In Western-blot analyses, serum antibodies from a human CE patient infected with E. canadensis (G6) strongly reacted with recombinant P29 from E. granulosus s.s. (G1) (recEg(G1)P29). In the same line, human anti-Eg(G1)P29 antibodies bound to recEcnd(G6)P29. Thus, minor AA sequence variations appear not to impair the prognostic serological use of P29.

Polymorphous oligodendroglioma of Zülch revisited: a genetically heterogeneous group of anaplastic gliomas including tumors of bona fide oligodendroglial differentiation.

A polymorphous variant of oligodendroglioma was described by K.J. Zülch half a century ago, and is only very sporadically referred to in the subsequent literature. In particular, no comprehensive analysis with respect to clinical or genetic features of these tumors is available. From a current perspective, the term polymorphous oligodendroglioma (pO) may appear as contradictory in terms, as nuclear monotony is a histomorphological hallmark of oligodendrogliomas. For the purpose of this study, we defined pO as diffusely infiltrating gliomas felt to be of oligodendroglial rather than astrocytic differentiation and characterized by the presence of multinucleate tumor giant cells and/or nuclear pleomorphism. In a total of nine patients, we identified tumors consistent with this working definition. All tumors were high-grade. We characterized these with respect to clinical, histomorphological and genetic features. Despite clinical and genetic heterogeneity, we identified a subset of tumors of bona fide oligodendroglial differentiation as characterized by combined loss of heterozygosity of chromosome arms 1p and 19q (LOH 1p19q). Those tumors that lacked LOH 1p19q showed a high frequency of IDH1 mutations and loss of alpha thalassemia/mental retardation syndrome X-linked gene (ATRX) immunoreactivity, indicating a possible phenotypic convergence of true oligodendrogliomas and gliomas of the alternative lengthening of telomeres (ALT) pathway. p53 alterations were common irrespective of the 1p19q status. Histomorphologically, the tumors featured interspersed bizarre multinucleate giant tumor cells, while the background population varied from monotonous to significantly pleomorphic. Our findings indicate, that a rare polymorphous - or "giant cell" - variant of oligodendroglioma does indeed exist.

Combined ATRX/IDH1 immunohistochemistry predicts genotype of oligoastrocytomas

AIMS To assess whether in oligoastrocytomas ATRX deficiency - as a surrogate of the alternative lengthening of telomeres (ALT) pathway - has a role in predicting the presence or absence of loss of heterozygosity of 1p and 19q (LOH), the genetic signature of oligodendroglial differentiation and a favourable prognostic marker. METHODS AND RESULTS A series of 54 oligoastrocytomas were investigated by immunohistochemistry as well as microsatellite analysis for LOH 1p19q. Genetic findings were correlated with morphological assessment. CONCLUSIONS ATRX deficiency was mutually exclusive with LOH. Conversely, ATRX-proficient tumours immunoreactive for R132H-mutant isocitrate dehydrogenase 1 (IDH1) showed a high rate (85%) of LOH. A more oligodendroglioma-like morphology was associated with a higher rate of LOH even in the morphologically ambiguous group of oligoastrocytomas. Our findings support the concept that oligoastrocytomas represent a morphological grey zone rather than a group of truly "mixed" or "intermediate" tumours. More precise classification of diffuse gliomas may also improve grading of borderline cases. We propose an immunohistochemical algorithm for classification of morphologically ambiguous diffuse gliomas. This article is protected by copyright. All rights reserved.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

Thrombophilia and risk of VTE recurrence according to the age at the time of first VTE manifestation

BACKGROUND Whether screening for thrombophilia is useful for patients after a first episode of venous thromboembolism (VTE) is a controversial issue. However, the impact of thrombophilia on the risk of recurrence may vary depending on the patient's age at the time of the first VTE. PATIENTS AND METHODS Of 1221 VTE patients (42 % males) registered in the MAISTHRO (MAin-ISar-THROmbosis) registry, 261 experienced VTE recurrence during a 5-year follow-up after the discontinuation of anticoagulant therapy. RESULTS Thrombophilia was more common among patients with VTE recurrence than those without (58.6 % vs. 50.3 %; p = 0.017). Stratifying patients by the age at the time of their initial VTE, Cox proportional hazards analyses adjusted for age, sex and the presence or absence of established risk factors revealed a heterozygous prothrombin (PT) G20210A mutation (hazard ratio (HR) 2.65; 95 %-confidence interval (CI) 1.71 - 4.12; p < 0.001), homozygosity/double heterozygosity for the factor V Leiden and/or PT mutation (HR 2.35; 95 %-CI 1.09 - 5.07, p = 0.030), and an antithrombin deficiency (HR 2.12; 95 %-CI 1.12 - 4.10; p = 0.021) to predict recurrent VTE in patients aged 40 years or older, whereas lupus anticoagulants (HR 3.05; 95%-CI 1.40 - 6.66; p = 0.005) increased the risk of recurrence in younger patients. Subgroup analyses revealed an increased risk of recurrence for a heterozygous factor V Leiden mutation only in young females without hormonal treatment whereas the predictive value of a heterozygous PT mutation was restricted to males over the age of 40 years. CONCLUSIONS Our data do not support a preference of younger patients for thrombophilia testing after a first venous thromboembolic event.

The Vacuolar H+-ATPase B1 Subunit Polymorphism p.E161K Associates with Impaired Urinary Acidification in Recurrent Stone Formers

Mutations in the vacuolar–type H+-ATPase B1 subunit gene ATP6V1B1 cause autosomal–recessive distal renal tubular acidosis (dRTA). We previously identified a single-nucleotide polymorphism (SNP) in the human B1 subunit (c.481G.A; p.E161K) that causes greatly diminished pump function in vitro. To investigate the effect of this SNP on urinary acidification, we conducted a genotype-phenotype analysis of recurrent stone formers in theDallas and Bern kidney stone registries. Of 555 patients examined, 32 (5.8%) were heterozygous for the p.E161K SNP, and the remaining 523 (94.2%) carried two wild–type alleles. After adjustment for sex, age, body mass index, and dietary acid and alkali intake, p.E161K SNP carriers had a nonsignificant tendency to higher urinary pH on a random diet (6.31 versus 6.09; P=0.09). Under an instructed low–Ca and low–Na diet, urinary pH was higher in p.E161K SNP carriers (6.56 versus 6.01; P,0.01). Kidney stones of p.E161K carriers were more likely to contain calcium phosphate than stones of wild-type patients. In acute NH4Cl loading, p.E161K carriers displayed a higher trough urinary pH (5.34 versus 4.89; P=0.01) than wild-type patients. Overall, 14.6% of wild-type patients and 52.4% of p.E161K carriers were unable to acidify their urine below pH 5.3 and thus, can be considered to have incomplete dRTA. In summary, our data indicate that recurrent stone formers with the vacuolar H+-ATPase B1 subunit p.E161K SNP exhibit a urinary acidification deficit with an increased prevalence of calcium phosphate– containing kidney stones. The burden of E161K heterozygosity may be a forme fruste of dRTA.

A significant effect of the KIR ligand HLA-C on fibrosis progression in chronic C hepatitis with or without liver transplantation.

BACKGROUND & AIMS The interaction of KIR with their HLA ligands drives the activation and inhibition of natural killer (NK) cells. NK cells could be implicated in the development of liver fibrosis in chronic hepatitis C. METHODS We analysed 206 non-transplanted and 53 liver transplanted patients, selected according to their Metavir fibrosis stage. Several variables such as the number of activator KIR or the HLA ligands were considered in multinomial and logistic regression models. Possible confounding variables were also investigated. RESULTS The KIRs were not significant predictors of the fibrosis stage. Conversely, a significant reduction of the HLA-C1C2 genotype was observed in the most advanced fibrosis stage group (F4) in both cohorts. Furthermore, the progression rate of fibrosis was almost 10 times faster in the subgroup of patients after liver transplantation and HLA-C1C2 was significantly reduced in this cohort compared to non-transplanted patients. CONCLUSION This study suggests a possible role of KIR and their ligands in the development of liver damage. The absence of C1 and C2 ligands heterozygosity could lead to less inhibition of NK cells and a quicker progression to a high level of fibrosis in patients infected by HCV, especially following liver transplantation. This article is protected by copyright. All rights reserved.

IGHD II: A Novel GH-1 Gene Mutation (GH-L76P) Severely Affects GH Folding, Stability, and Secretion

CONTEXT The autosomal dominant form of GH deficiency (IGHD II) is characterized by markedly reduced GH secretion combined with low concentrations of IGF-1 leading to short stature. OBJECTIVE Structure-function analysis of a missense mutation in the GH-1 gene converting codon 76 from leucine (L) to proline (P) yielding a mutant GH-L76P peptide. DESIGN, SETTINGS, AND PATIENTS Heterozygosity for GH-L76P/wt-GH was identified in a nonconsanguineous Spanish family. The index patients, two siblings, a boy and a girl, were referred for assessment of their short stature (-3.2 and -3.8 SD). Their grandmother, father, and aunt were also carrying the same mutation and showed severe short stature; therefore, IGHD II was diagnosed. INTERVENTIONS AND RESULTS AtT-20 cells coexpressing both wt-GH and GH-L76P showed a reduced GH secretion (P < .001) after forskolin stimulation when compared with the cells expressing only wt-GH. In silico mutagenesis and molecular dynamics simulations presented alterations of correct folding and mutant stability compared with wt-GH. Therefore, further structural analysis of the GH-L76P mutant was performed using expressed and purified proteins in Escherichia coli by thermofluor assay and fast degradation proteolysis assay. Both assays revealed that the GH-L76P mutant is unstable and misfolded compared to wt-GH confirming the bioinformatic model prediction. CONCLUSIONS This is the first report of a family suffering from short stature caused by IGHD II, which severely affects intracellular GH folding and stability as well as secretion, highlighting the necessity of functional analysis of any GH variant for defining new mechanisms as a cause for IGHD II.

Expansion load: recessive mutations and the role of standing genetic variation

Expanding populations incur a mutation burden – the so-called expansion load. Previous studies of expansion load have focused on codominant mutations. An important consequence of this assumption is that expansion load stems exclusively from the accumulation of new mutations occurring in individuals living at the wave front. Using individual-based simulations, we study here the dynamics of standing genetic variation at the front of expansions, and its consequences on mean fitness if mutations are recessive. We find that deleterious genetic diversity is quickly lost at the front of the expansion, but the loss of deleterious mutations at some loci is compensated by an increase of their frequencies at other loci. The frequency of deleterious homozygotes therefore increases along the expansion axis, whereas the average number of deleterious mutations per individual remains nearly constant across the species range. This reveals two important differences to codominant models: (i) mean fitness at the front of the expansion drops much faster if mutations are recessive, and (ii) mutation load can increase during the expansion even if the total number of deleterious mutations per individual remains constant. We use our model to make predictions about the shape of the site frequency spectrum at the front of range expansion, and about correlations between heterozygosity and fitness in different parts of the species range. Importantly, these predictions provide opportunities to empirically validate our theoretical results. We discuss our findings in the light of recent results on the distribution of deleterious genetic variation across human populations and link them to empirical results on the correlation of heterozygosity and fitness found in many natural range expansions.

Implications of Longterm Diameter-Limit Harvesting: Effects on Radial Growth of Red Spruce (Picea rubens) and Genetic Diversity of White Pine (Pinus strobus)

For over 3 centuries, diameter-limit harvesting has been a predominant logging method in the northeastern United States. Silvicultural theory asserts that such intensively selective harvesting can lead to genetic degradation. A decrease in softwood productivity has recently been reported in Maine - has a long history of dysgenic selection degraded the genetic resources of Maine softwoods, contributing to a decrease in growth and productivity? This study examines two aspects of potential implications of diameter-limit harvesting: effects on residual phenotypes of red spruce and impacts on genetic diversity of white pine. Radial growth of residual red spruce trees in stands experiencing 50 years of fixed diameter-limit harvesting was measured using annual increment rings and compared with residual red spruce trees in positive selection stands. Trees remaiaing after several rounds of diameter-limit harvesting exhibited sigdicantl y smaller radial sizes throughout their lives, and displayed significantly slower growth rates for the first 80 years of measured growth. These results strongly suggest that the largest and fastest-growing genotypes and their respective gene complexes determining good radial growth have been removed from the diameter-limit stand. Dysgenic selection can be observed in fixed diarneter-limit stands, resulting in a diminished genetic resource and decreased residual stand value. To examine more direct genetic implications of long-term diameter-limit harvesting, microsatellite DNA markers were implemented to study genetic diversity of eastern white pine in Maine. Three age groups of trees were studied: mature trees older than 200 years, juvenile trees 5-30 years old, and embryos. Trees were genotyped at 10 microsatellite loci. Overall genetic diversity levels of eastern white pine in Maine were extremely high, with an average observed heterozygosity of 0.762. Genetic differentiation was minimal among and between all three age groups, although an excess of heterozygotes was shown in the mature and juvenile groups that was not reflected in the embryo group, which actually had a slight heterozygote deficiency. Allele frequencies did not differ significantly between age groups, but did reveal more rare and low frequency alleles in the embryo groups than in the mature group. Overall, low frequency alleles comprise the largest portion of alleles in the sample population, with no common alleles evident overall. These results suggest that significant genetic degradation has either not occurred for white pine, or that the results of dysgenic selection have not yet emerged. It is clear, however, that selective harvesting could result in a loss of low frequency alleles, which are a primary reserve of evolutionary potential in a species. Implications of these studies affect industrial forestry, regional economics, and ecological concerns for the northeast. Long-term diameter-limit harvesting can lead to a degradation of residual phenotypes, and an overall decrease in stand quality. Potentially, a loss of low frequency, locally adapted alleles could result in a decrease of allelic richness and degradation of the regidnal genetic resource. Decreased genetic variation can lead to seriously limited evolutionary potential of species and ecosystems, particularly in rapidly changing environments. Based on these findings, I recommend a reassessment of any harvesting prescription that includes fixed diameter-limit removals, particularly for species that have low natural genetic diversity levels or a limited natural range, such as red spruce. Maintenance of a healthy genetic reserve can avoid effects of dysgenic harvesting.

The localization and identification of candidate tumor suppressor genes involved in prostate cancer

Frequent loss of heterozygosity (LOH) at specific chromosomal regions are highly associated with the inactivation of tumor suppressor genes (TSGs) (Weinberg, 1991; Bishop, 1989). Chromosome 8p is the most frequently reported site of LOH (∼60%) in prostate cancer (PC), suggesting that there may be inactivated TSG(s) involved in PC on chromosome 8p. (Bergerheim et. al., 1991; Kagan et. al., 1995). In order to identify the smallest common regions of frequent LOH (SCLs) on chromosome 8, we screened 52 PC patient/tumor samples with 39 polymorphic markers in successive screenings. In the course of refining the SCLs, we identified 3 tumors with >6 Mb homozygous deletions (HZDs) at 8p22 and 8p21, suggesting the presence of candidate TSGs at both loci. These HZDs spanned the two SCLs at 8p22 (46%) and 8p21 (45%). The SCLs were narrowed to 3.2 cM at 8p22 and less than 3 cM at 8p21. ^ In order to identify candidate TSGs within the SCLs on 8p, two approaches were used. In the candidate gene approach, thirty genes that mapped to the SCLs were evaluated for expression in normal prostate and in PC cell lines. One of the candidate genes, Clusterin, showed decreased expression in 4/7 (57%) prostate cancer cell lines by Northern blot analysis. Clusterin will be further examined as a candidate TSG. ^ The second approach involved utilizing subtractive hybridization and hybrid affinity capture to generate pools of expressed sequence tags (ESTs) enriched for genes that are downregulated or deleted in PC and that map to specific regions of interest. We took advantage of a prostate cancer cell line (PC3) with a known HZD of a candidate TSG, CTNNA1 on 5q31, to develop and validate a model system. We then developed subtracted libraries enriched for 8p22 and 8p21 ESTs by this method, using two cell lines, MDAPCa-2b and PC3. The ESTs were cloned, and 40 were sequenced and evaluated for expression in normal prostate and PC cell lines. Three ESTs from the subtracted libraries, C2, C17 and F12, showed decreased expression in 29–57% of the prostate tumor cell lines studied, and will be further examined as candidate TSGs. ^

Evaluation of candidate genes for susceptibility to spina bifida

Neural tube defects (NTDs) are the most common severely disabling birth defects in the United States, with a frequency of approximately 1–2 of every 1,000 births. This text includes the identification and evaluation of candidate susceptibility genes that confer risk for the development of neural tube defects (NTDs). The project focused on isolated meningomyelocele, also termed spina bifida (SB). ^ Spina bifida is a complex disease with multifactorial inheritance, therefore the subject population (consisting of North American Caucasians and Hispanics of Mexicali-American descent) was composed of 459 simplex SB families who were tested for genetic associations utilizing the transmission disequilibrium test (TDT), a nonparametric linkage technique. Three categories of candidate genes were studied, including (1) human equivalents of genes determined in mouse models to cause NTDs, (2) HOX and PAX genes, and (3) the MTHFR gene involved in the metabolic pathway of folate. ^ The C677T variant of the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene was the first mutation in this gene to be implicated as a risk factor for NTDs. Our evaluation of the MTHFR gene provides evidence that maternal C677T homozygosity is a risk factor for upper level spina bifida defects in Hispanics [OR = 2.3, P = 0.02]. This observed risk factor is of great importance due to the high prevalence of this homozygous genotype in the Hispanic population. Additionally, maternal C677T/A1298C compound heterozygosity is a risk factor for upper level spina bifida defects in non-Hispanic whites [OR = 3.6, P = 0.03]. ^ For TDT analysis, our total population of 1128 subjects were genotyped for 54 markers from within and/or flanking the 20 candidate genes/gene regions of interest. Significant TDT findings were obtained for 3 of the 54 analyzed markers: d20s101 flanking the PAX1 gene (P = 0.019), d1s228 within the PAX7 gene (P = 0.011), and d2s110 within the PAX8 gene (P = 0.013). These results were followed-up by testing the genes directly for mutations utilizing single-strand conformational analysis (SSCA) and direct sequencing. Multiple variations were detected in each of these PAX genes; however, these variations were not passed from parent to child in phase with the positively transmitted alleles. Therefore, these variations do not contribute to the susceptibility of spina bifida, but rather are previously unreported single nucleotide polymorphisms. ^