Complex Genetic Disorders: Evaluating When Genetic Research Findings Are Applicable for Genetic Counseling Practice

Traditional genetic counseling processes and principles will be extended to a new realm-complex disorders. Although it may seem like a daunting task, understanding the methodologies used to study complex genetic disorders will enable genetic counselors to critically analyze research studies involving complex disorders. In this article, we explain newly evolving methodologies for genetic research, including case-control studies and transmission disequilibrium testing (TDT). Additionally, a framework is provided for evaluating original research findings and replication studies.

A functional alternative splicing mutation in human tryptophan hydroxylase-2

The brain serotonergic system has an essential role in the physiological functions of the central nervous system and dysregulation of serotonin (5-HT) homeostasis has been implicated in many neuropsychiatric disorders. The tryptophan hydroxylase-2 (TPH2) gene is the rate-limiting enzyme in brain 5-HT synthesis, and thus is an ideal candidate gene for understanding the role of dysregulation of brain serotonergic homeostasis. Here, we characterized a common, but functional single-nucleotide polymorphism (SNP rs1386493) in the TPH2 gene, which decreases efficiency of normal RNA splicing, resulting in a truncated TPH2 protein (TPH2-TR) by alternative splicing. TPH2-TR, which lacks TPH2 enzyme activity, dominant-negatively affects full-length TPH2 function, causing reduced 5-HT production. The predicted mRNA for TPH2-TR is present in postmortem brain of rs1386493 carriers. The rs13864923 variant does not appear to be overrepresented in either global or multiplex depression cohorts. However, in combination with other gene variants linked to 5-HT homeostasis, this variant may exhibit important epistatic influences.

Maternal fumonisin exposure as a risk factor for neural tube defects

Fumonisins are mycotoxins produced by the fungus F. verticillioides, a common contaminant of maize (corn) worldwide. Maternal consumption of fumonisin B(1)-contaminated maize during early pregnancy has recently been associated with increased risk for neural tube defects (NTDs) in human populations that rely heavily on maize as a dietary staple. Experimental administration of purified fumonisin to mice early in gestation also results in an increased incidence of NTDs in exposed offspring. Fumonisin inhibits the enzyme ceramide synthase in de novo sphingolipid biosynthesis, resulting in an elevation of free sphingoid bases and depletion of downstream glycosphingolipids. Increased sphingoid base metabolites (i.e., sphinganine-1-phosphate) may perturb signaling cascades involved in embryonic morphogenesis by functioning as ligands for sphingosine-1-P (S1P) receptors, a family of G-protein-coupled receptors that regulate key biological processes such as cell survival/proliferation, differentiation and migration. Fumonisin-induced depletion of glycosphingolipids impairs expression and function of the GPI-anchored folate receptor (Folr1), which may also contribute to adverse pregnancy outcomes. NTDs appear to be multifactorial in origin, involving complex gene-nutrient-environment interactions. Vitamin supplements containing folic acid have been shown to reduce the occurrence of NTDs, and may help protect the developing fetus from environmental teratogens. Fumonisins appear to be an environmental risk factor for birth defects, although other aspects of maternal nutrition and genetics play interactive roles in determining pregnancy outcome. Minimizing exposures to mycotoxins through enhanced agricultural practices, identifying biomarkers of exposure, characterizing mechanisms of toxicity, and improving maternal nutrition are all important strategies for reducing the NTD burden in susceptible human populations.

Novel congenital myopathy locus identified in Native American Indians at 12q13.13-14.1

BACKGROUND

Native American myopathy (NAM) is an autosomal recessive congenital myopathy first reported in the Lumbee Indian people. Features of NAM include congenital weakness, cleft palate, ptosis, short stature, and susceptibility to malignant hyperthermia provoked by anesthesia.

METHOD

We identified five individuals with NAM from the Lumbee population, and hypothesized that these affected individuals have disease alleles shared identical-by-descent inherited from common ancestry. To identify a NAM disease locus, homozygosity mapping methods were employed on a genomewide 10K single-nucleotide polymorphism (SNP) screen. To confirm regions of homozygosity identified in the SNP screen, microsatellite repeat markers were genotyped within those homozygous segments.

RESULTS

The SNP data demonstrated five regions of shared homozygosity in individuals with NAM. The additional genotyping data narrowed the region to one common segment of homozygosity spanning D12S398 to rs3842936 mapping to 12q13.13-14.1. Notably, loss of heterozygosity estimates from the SNP data also detected this same 12q region in the affected individuals.

CONCLUSION

This study reports the first gene mapping of Native American myopathy (NAM) using single-nucleotide polymorphism-based homozygosity mapping in only a few affected individuals from simplex families and identified a novel NAM locus. Identifying the genetic basis of NAM may suggest new genetic etiologies for other more common conditions such as congenital myopathy and malignant hyperthermia.

Genotype-phenotype study in an FSHD family with a proximal deletion encompassing p13E-11 and D4Z4

BACKGROUND

In the majority of facioscapulohumeral muscular dystrophy (FSHD) cases, the molecular basis of the disease is due to loss of subtelomeric D4Z4 repeat units at 4q35. Occasionally, an apparent absence of the contracted D4Z4 repeat is associated with FSHD. One explanation for this finding is a deletion in the region proximal to the D4Z4 repeat array that encompasses the p13E-11 (D4F104S1) probe-binding site used in the DNA diagnosis. The frequency of such proximally extended deletions is unknown, and to date, few patients have been described due to the difficulties in the molecular identification of such cases.

METHODS

We describe a family (DUK 2531) in which a contracted D4Z4 allele and a large proximal deletion of approximately 75 kb are segregating to 11 individuals. This is the largest deletion identified to date. Family DUK 2531 was initially thought to have normal D4Z4 fragment size and therefore unlinked to the 4q35 region (FSHD1B).

RESULTS

Further molecular analysis of DUK 2531 reveals the presence of 10 repeat units (33 kb). The extended deletion includes the probe p13E-11 and B31 binding sites, the inverted repeat D4S2463, and genes FRG2 and TUBB4Q.

CONCLUSION

Despite the length of the proximal deletion in this family, the range and severity of the clinical manifestations are typical for the disorder. Because such deletions can lead to misinterpretation in the diagnostic setting, this suggests the need for additional diagnostic tests in facioscapulohumeral muscular dystrophy.

Power calculations for likelihood ratio tests for offspring genotype risks, maternal effects, and parent-of-origin (POO) effects in the presence of missing parental genotypes when unaffected siblings are available

Genotype-based likelihood-ratio tests (LRT) of association that examine maternal and parent-of-origin effects have been previously developed in the framework of log-linear and conditional logistic regression models. In the situation where parental genotypes are missing, the expectation-maximization (EM) algorithm has been incorporated in the log-linear approach to allow incomplete triads to contribute to the LRT. We present an extension to this model which we call the Combined_LRT that incorporates additional information from the genotypes of unaffected siblings to improve assignment of incompletely typed families to mating type categories, thereby improving inference of missing parental data. Using simulations involving a realistic array of family structures, we demonstrate the validity of the Combined_LRT under the null hypothesis of no association and provide power comparisons under varying levels of missing data and using sibling genotype data. We demonstrate the improved power of the Combined_LRT compared with the family-based association test (FBAT), another widely used association test. Lastly, we apply the Combined_LRT to a candidate gene analysis in Autism families, some of which have missing parental genotypes. We conclude that the proposed log-linear model will be an important tool for future candidate gene studies, for many complex diseases where unaffected siblings can often be ascertained and where epigenetic factors such as imprinting may play a role in disease etiology.

Potential for expanded power in linkage studies using the ALLEGRO and MERLIN software programs

Multipoint linkage analysis in complex diseases requires the use of fast algorithms that can handle many markers and a large number of moderately sized pedigrees with unknown mode of inheritance. This need has led to the development of several competitive software programs. We recently completed a genomic screen of neural tube defects using GENEHUNTER-PLUS and the more recent ALLEGRO. The ALLEGRO software was found to offer expanded power for linkage studies, particularly for childhood onset diseases like neural tube defects, though the results must be treated with caution.

Whole genomewide linkage screen for neural tube defects reveals regions of interest on chromosomes 7 and 10

Neural tube defects (NTDs) are the second most common birth defects (1 in 1000 live births) in the world. Periconceptional maternal folate supplementation reduces NTD risk by 50-70%; however, studies of folate related and other developmental genes in humans have failed to definitively identify a major causal gene for NTD. The aetiology of NTDs remains unknown and both genetic and environmental factors are implicated. We present findings from a microsatellite based screen of 44 multiplex pedigrees ascertained through the NTD Collaborative Group. For the linkage analysis, we defined our phenotype narrowly by considering individuals with a lumbosacral level myelomeningocele as affected, then we expanded the phenotype to include all types of NTDs. Two point parametric analyses were performed using VITESSE and HOMOG. Multipoint parametric and nonparametric analyses were performed using ALLEGRO. Initial results identified chromosomes 7 and 10, both with maximum parametric multipoint lod scores (Mlod) >2.0. Chromosome 7 produced the highest score in the 24 cM interval between D7S3056 and D7S3051 (parametric Mlod 2.45; nonparametric Mlod 1.89). Further investigation demonstrated that results on chromosome 7 were being primarily driven by a single large pedigree (parametric Mlod 2.40). When this family was removed from analysis, chromosome 10 was the most interesting region, with a peak Mlod of 2.25 at D10S1731. Based on mouse human synteny, two candidate genes (Meox2, Twist1) were identified on chromosome 7. A review of public databases revealed three biologically plausible candidates (FGFR2, GFRA1, Pax2) on chromosome 10. The results from this screen provide valuable positional data for prioritisation of candidate gene assessment in future studies of NTDs.

The genes encoding for D4Z4 binding proteins HMGB2, YY1, NCL, and MYOD1 are excluded as candidate genes for FSHD1B

Facioscapulohumeral muscular dystrophy is a disease of skeletal muscle, with symptoms including both facial and shoulder girdle weakness and progression to involve the pelvic girdle and extremities in the majority of cases. For most cases of FSHD, the molecular basis of the disease can be identified as a partial deletion of the D4Z4 repeat array on the end of the long arm of chromosome 4. However, in up to 5% of FSHD families there is no linkage to 4q35. These cases are designated as FSHD1B. Proteins have been identified that bind to the D4Z4 repeats of chromosome 4q35. The genes encoding D4Z4 binding proteins YY1, HMGB2, NCL, and MYOD1 were investigated as candidate genes for FSHD1B. Coding sequences and promoter region were analyzed for HMBG2 and no sequence variations were detected. For YY1, all five exons were analyzed and a polymorphism was detected in both the unaffected and affected populations. In nucleolin (NCL), several SNPs were identified, including a SNP causing the non-synonymous change P515H; however, all polymorphisms either occurred in control samples or were previously reported. A novel polymorphism was also detected in MYOD1, but did not represent a disease-specific variation. These results suggest that HMBG2, YY1, NCL, and MYOD1 are unlikely to represent the genes responsible for FSHD in these families.

Cloning and characterization of an inversion breakpoint at 6q23.3 suggests a role for Map7 in sacral dysgenesis

Here we report on a male patient with sacral dysgenesis (SD) and constitutional pericentric inversion of chromosome 6 (p11.2;q23.3). SD is a heterogeneous group of congenital anomalies with complex genetic etiology. Previously, a patient with sacral abnormalities and an interstitial deletion of 6q23-->q25 region has been described. We speculated that a susceptibility gene for SD lies in 6q23.3 region (disrupted in both patients), and therefore, cloning of the breakpoint in our patient would lead to the identification of the disrupted gene. We performed FISH analysis followed by Southern blot analysis and inverse PCR to clone the breakpoint. The 6p11.2 breakpoint mapped very close to the centromere, and the 6q23.3 breakpoint localized in the ninth intron of the MAP7 gene. We then evaluated the involvement of MAP7 in SD by further screening of the gene in several patients with a similar phenotype. Two nucleotide changes causing Ile257Asn and Glu571Ala substitutions in the protein, both affecting amino acid residues conserved in the mouse homolog, were identified in two patients. Both changes are either very rare polymorphisms or true mutations, since they were not detected in 167 normal individuals nor found in the SNP database. Therefore, our study suggests MAP7 as a candidate gene for SD. However, we were unable to detect any sacral defects in the MAP7 knockout mice.

Updated investigations of the role of methylenetetrahydrofolate reductase in human neural tube defects

Folate supplementation appears to reduce the risk for neural tube defects (NTDs). Methylenetetrahydrofolate reductase (MTHFR) is a candidate gene in the folate metabolism pathway that has been extensively studied in different human populations. We examined the risk associated with having the thermolabile variant (TT) of MTHFR in a study of 175 American Caucasians with NTDs and their families. We found a significant association in patients compared with 195 unrelated controls [odds ratio (OR) = 2.13, 95% confidence interval (95% CI) = 1.11-4.09)], but not in mothers (OR = 1.29, 95% CI = 0.622-2.67) or in fathers (OR = 1.45, 95% CI = 0.681-3.09). We found no evidence for unequal transmission from parents to an affected child (p > 0.10). We failed to find a previously reported association for a combined haplotype for MTHFR and cystathionine beta-synthase, except in subjects with NTDs compared with 559 pooled controls (OR = 2.87, 95% CI = 1.03-8.03). We found no evidence for an association for a novel CA-repeat polymorphism identified in a gene closely linked to MTHFR (p > 0.10). Our studies continue to suggest that additional candidate genes other than MTHFR may be responsible for an increased risk to NTD in some American Caucasian families.

Life after the screen: making sense of many P-values

A multiple analytic approach may be useful for analyzing complex traits since different methods extract both similar and distinct, but complementary pieces of information from genome screen data on extended pedigrees. We examined the usefulness of combining p-values both across methods and across adjacent markers, taking into account the observed correlation structure among these p-values. To this end, we employed the recently proposed truncated product method [Zaykin et al., Genet Epidemiol, in press]. It appears that this approach is helpful for visualizing priority regions for follow-up analysis and reducing the number of false-positive linkage signals.

Additional glomangioma families link to chromosome 1p: no evidence for genetic heterogeneity

Venous malformations are a common abnormality of the vasculature that may occur sporadically or, more rarely, as an autosomal dominant trait. One familial form of venous malformations has previously been linked to chromosome 9p. Mutations in the gene encoding Tie2, an endothelial specific receptor tyrosine kinase, have been identified in four different families. Glomangiomas are a subtype of venous malformations with glomus cell involvement. These cutaneous lesions can be inherited as an autosomal dominant disease with reduced penetrance and variable expressivity. We present evidence of linkage to chromosome 1p21-1p22 using four new glomangioma families, with a combined maximum two-point lod score of 7.32 at marker D1S2804. Markers D1S2129 and D1S2881 define the 24-cM linkage interval determined by recombination within affected individuals. A recent report also showed linkage of the glomangioma locus to chromosome 1p. A total of 9 families now map to this region, suggesting a decreased likelihood of locus heterogenity in familial glomangiomas. Investigation of candidate genes within the interval should provide new insights into lesion formation in inherited venous malformations.

Myotilin is mutated in limb girdle muscular dystrophy 1A

We have identified a mutation in the myotilin gene in a large North American family of German descent expressing an autosomal dominant form of limb girdle muscular dystrophy (LGMD1A). We have previously mapped this gene to 5q31. Symptoms of this adult onset disease are progressive weakness of the hip and shoulder girdles, as well as a distinctive dysarthric pattern of speech. Muscle of affected individuals shows degeneration of myofibers, variations in fiber size, fiber splitting, centrally located myonuclei and a large number of autophagic vesicles. Affected muscle also exhibits disorganization and streaming of the Z-line similar to that seen in nemaline myopathy. We have identified a C450T missense mutation in the myotilin gene that is predicted to result in the conversion of residue 57 from threonine to isoleucine. This mutation has not been found in 396 control chromosomes. The mutant allele is transcribed and normal levels of correctly localized myotilin protein are seen in LGMD1A muscle. Myotilin is a sarcomeric protein that binds to alpha-actinin and is localized in the Z-line. The observed missense mutation does not disrupt binding to alpha-actinin.

GAW10: simulated family data for a common oligogenic disease with quantitative risk factors

GAW10 Problem 2 involves a simulated common disease defined by imposing a threshold, T, on a quantitative trait, Q1. Every individual with a value of Q1 > or = T (where T = 40) is defined as affected. Also thought to be associated with the disease as intervening variables are four other quantitative traits (Q2, Q3, Q4, and Q5) and an environmental factor (EF). Each individual has genotypes for 367 highly polymorphic markers on 10 chromosomes. The tasks for GAW10 were to characterize the genetic and environmental contribution to Q1 (and the disease), Q2, Q3, Q4, and Q5, and to localize any single genes that have detectable effects. Multiple replicates provided an opportunity to address questions of power and false positive rates.

Heterogeneity in Paget disease of the bone

Paget disease of the bone is a common skeletal disorder. Recently, a gene for Paget disease was localized to 18q with subsequent evidence for linkage heterogeneity. We report the identification and clinical characterization of a large pedigree of Paget disease and demonstrate that the Paget disease gene in this pedigree is not linked to the region on 18q, thus confirming linkage heterogeneity.

A Bethlem myopathy Gly to Glu mutation in the von Willebrand factor A domain N2 of the collagen alpha3(VI) chain interferes with protein folding

A single G1679E mutation in the amino-terminal globular domain N2 of the alpha3 chain of type VI collagen was found in a large family affected with Bethlem myopathy. Recombinant production of N2 ( approximately 200 residues) in transfected mammalian cells has now been used to examine the possibility that the mutation interfered with protein folding. The wild-type form and a G1679A mutant were produced at high levels and shown to fold into a stable globular structure. Only a small amount of secretion was observed for mutants G1679E and G1679Q, which apparently were efficiently degraded within the cells. Homology modeling onto the related von Willebrand factor A1 structure indicated that substitution of G1679 by the bulky E or Q cannot be accommodated without considerable changes in the folding pattern. This suggests protein misfolding as a molecular basis for this particular mutation in Bethlem myopathy, in agreement with radioimmunoassay data showing reduced levels of domain N2 in cultured fibroblasts from two patients.

Myelocystocele-cloacal exstrophy in a pedigree with a mitochondrial 12S rRNA mutation, aminoglycoside-induced deafness, pigmentary disturbances, and spinal anomalies

A large Filipino-American family with progressive matrilineal hearing loss, premature graying, depigmented patches, and digital anomalies was ascertained through a survey of a spina bifida clinic for neural crest disorders. Deafness followed a matrilineal pattern of inheritance and was associated with the A1555G mutation in the 12S rRNA gene (MTRNR1) in affected individuals as well as unaffected maternal relatives. Several other malformations were found in carriers of the mutation. The proband had a myelocystocele, Arnold-Chiari type I malformation, cloacal exstrophy, and severe early-onset hearing loss. Several family members had premature graying, white forelock, congenital leukoderma with or without telecanthus, somewhat suggestive of a Waardenburg syndrome variant. In addition to the patient with myelocystocele, two individuals had scoliosis and one had segmentation defects of spinal vertebrae. The syndromic characteristics reported here are novel for the mitochondrial A1555G substitution, and may result from dysfunction of mitochondrial genes during early development. However, the mitochondrial A1555G mutation is only rarely associated with neural tube defects as it was not found in a screen of 218 additional individuals with spina bifida, four of whom had congenital hearing loss.

A genetic hypothesis for Chiari I malformation with or without syringomyelia

In several reports the authors have suggested occasional familial aggregation of syringomyelia and/or Chiari 1 malformation (CM1). Familial aggregation is one characteristic of traits that have an underlying genetic basis. The authors provide evidence for familial aggregation of CM1 and syringomyelia (CM1/S) in a large series of families, establishing that there may be a genetic component to CM1/S in at least a subset of families. The authors observed no cases of isolated familial syringomyelia in their family studies, suggesting that familial syringomyelia is more accurately classified as familial CM1 with associated syringomyelia.

These data, together with the cosegregation of the trait with known genetic syndromes, support the authors' hypothesis of a genetic basis for some CM1/S cases.

Genetic and embryological approaches to studies of neural tube defects: a critical review. NTD Collaborative Group

Experimental embryological models have suggested that the morphology and quantity of neural tube defects may be governed by their position along the anteroposterior axis of the embryo. Inductive interactions and genetic regulation during axis development may play a role in the patterning of neural tube defects. A major challenge in the study of human neural tube defects is determining whether the spectrum of developmental neural tube anomalies found in individuals and their families mirror experimental models and are regulated by similar processes. We have found that the various neural tube defect phenotypes can be clustered according to their position along the anteroposterior axis. The findings correlate well to the pattern of early genes expression, inductive models of the embryonic axis, and mutant NTD animal models. We suggest that NTD should be studied by their location along the anteroposterior axis and that specific mutant genes may be identified by the observed pattern of NTD in an individual or a family.

Genetic studies in neural tube defects. NTD Collaborative Group

Neural tube defects (NTD) are one of the most common birth defects and are caused by both environmental and genetic factors. The approach to identifying the genes predisposing to NTD, through linkage analysis and candidate gene analysis, is reviewed along with characteristics of a large, nationally ascertained cohort of families. Results from specific assessments of p53, PAX3 and MTHFR failed to suggest that these genes play a major role in NTD development in these families. Advances in genetic laboratory and statistical techniques have made this a prime opportunity for investigation into the causes of complex disorders, such as NTD. However, traditional approaches may prove to be challenging due to the difficulty of ascertaining samplable multiplex families.

Complete genomic screen for disease susceptibility loci in nuclear families

We performed genome-wide model dependent and independent analyses on a simulated data set of 400 families segregating for a rare disorder. Regions on chromosomes 1, 3, and 5 were consistently indicated across the various analyses performed. Follow-up analyses included stratification for locus heterogeneity and clinical phenotype and studies of gene x gene and gene x environment interaction. The region around D1G024 was most notable, showing strong association and linkage with the trait. We also identified regions D3G043-46 and D5G037-39 by strong linkage and association findings and region D1G001-09 by linkage analysis. A complex statistical interaction was suggested between D1G024, D3G046 and environmental factor 1. This report suggests that traditional methods of analysis can be implemented to analyze and describe the mechanisms that may underlie the more complex genetic disorders.

Possible interaction of genotypes at cystathionine beta-synthase and methylenetetrahydrofolate reductase (MTHFR) in neural tube defects. NTD Collaborative Group

Neural tube defects are a common, complex disorder with genetic and environmental components to risk. We investigated the previously reported interaction between homozygosity for the thermolabile variant at the methylenetetrahydrofolate reductase and heterozygosity for the 844ins68 allele at the cystathionine beta-synthase loci in cases with lumbosacral myelomeningocele and their parents. Using control allele frequencies from our sample pooled with those published in the literature, we confirm a marginally significant interaction at these two loci. This finding suggests that additional, larger studies are warranted to investigate this possible interaction in more detail.

Allelic and locus heterogeneity in inherited venous malformations

Venous malformations are low-flow vascular lesions consisting of disorganized thin-walled vascular channels. These can occur sporadically but also as an autosomal dominant condition termed venous malformations, cutaneous and mucosal (VMCM; OMIM 600195). In two large unrelated kindreds mapping to chromosome 9, the identical R849W missense mutation was identified in the first kinase domain of Tie2, an endothelial cell-specific receptor tyrosine kinase. We report here the identification of four new kindreds with inherited venous malformations. Unlike the initial two families described, these four families demonstrate allelic and locus heterogeneity. In one of these families, the R849W mutation co-segregates with the disease phenotype. Three other families with venous malformations lack this mutation. One of these families is linked to markers near TIE2 on chromosome 9. In this family, we identified a novel mutation within the first kinase domain of Tie2 resulting in a Y897S change. Results from COS-1 cell transfections using expression constructs containing either the R849W or the Y897S mutation suggest that the receptors containing either mutation show ligand-independent hyperphosphorylation. These results suggest a gain-of-function mechanism for development of venous malformations in these families. Of the two remaining families, one excludes linkage to the TIE2 locus, establishing the existence of at least one additional locus for dominantly inherited venous malformations.

Chiari I malformation redefined: clinical and radiographic findings for 364 symptomatic patients

OBJECTIVE

Chiari malformations are regarded as a pathological continuum of hindbrain maldevelopments characterized by downward herniation of the cerebellar tonsils. The Chiari I malformation (CMI) is defined as tonsillar herniation of at least 3 to 5 mm below the foramen magnum. Increased detection of CMI has emphasized the need for more information regarding the clinical features of the disorder.

METHODS

We examined a prospective cohort of 364 symptomatic patients. All patients underwent magnetic resonance imaging of the head and spine, and some were evaluated using CINE-magnetic resonance imaging and other neurodiagnostic tests. For 50 patients and 50 age- and gender-matched control subjects, the volume of the posterior cranial fossa was calculated by the Cavalieri method. The families of 21 patients participated in a study of familial aggregation.

RESULTS

There were 275 female and 89 male patients. The age of onset was 24.9+/-15.8 years (mean +/- standard deviation), and 89 patients (24%) cited trauma as the precipitating event. Common associated problems included syringomyelia (65%), scoliosis (42%), and basilar invagination (12%). Forty-three patients (12%) reported positive family histories of CMI or syringomyelia. Pedigrees for 21 families showed patterns consistent with autosomal dominant or recessive inheritance. The clinical syndrome of CMI was found to consist of the following: 1) headaches, 2) pseudotumor-like episodes, 3) a Meniere's disease-like syndrome, 4) lower cranial nerve signs, and 5) spinal cord disturbances in the absence of syringomyelia. The most consistent magnetic resonance imaging findings were obliteration of the retrocerebellar cerebrospinal fluid spaces (364 patients), tonsillar herniation of at least 5 mm (332 patients), and varying degrees of cranial base dysplasia. Volumetric calculations for the posterior cranial fossa revealed a significant reduction of total volume (mean, 13.4 ml) and a 40% reduction of cerebrospinal fluid volume (mean, 10.8 ml), with normal brain volume.

CONCLUSION

These data support accumulating evidence that CMI is a disorder of the para-axial mesoderm that is characterized by underdevelopment of the posterior cranial fossa and overcrowding of the normally developed hindbrain. Tonsillar herniation of less than 5 mm does not exclude the diagnosis. Clinical manifestations of CMI seem to be related to cerebrospinal fluid disturbances (which are responsible for headaches, pseudotumor-like episodes, endolymphatic hydrops, syringomyelia, and hydrocephalus) and direct compression of nervous tissue. The demonstration of familial aggregation suggests a genetic component of transmission.

A radiation hybrid breakpoint map of the acute myeloid leukemia (AML) and limb-girdle muscular dystrophy 1A (LGMD1A) regions of chromosome 5q31 localizing 122 expressed sequences

We have constructed a high-resolution map of a 6-Mb interval of human chromosome 5, band q31, incorporating 175 sequence tagged sites, of which 33 are genetic polymorphisms and 122 are nonredundant expressed sequences. The map was assembled initially as a YAC contig, incorporating data from radiation hybrid maps. To improve resolution and to identify errors in the databases, a radiation hybrid breakpoint map was developed for the interval, which included hybrids from both Stanford G3 and GeneBridge 4 panels. This novel approach facilitated the integration of one RH panel with another and enabled the identification and localization of new, previously unmapped ESTs from the radiation hybrid databases. ESTs were assembled into overlapping transcription units and ordered with respect to polymorphic markers in the region, resulting in a comprehensive map that incorporates markers from multiple different types of maps. This map of 5q31 will facilitate gene discovery efforts for several disorders, including limb-girdle muscular dystrophy type 1A and the genes deleted in acute myeloid leukemias and myelodysplasia. The study demonstrates the utility of a radiation hybrid breakpoint panel for correction of map errors and for the efficient identification of new transcript units in a large genomic interval.

Identification of a new autosomal dominant limb-girdle muscular dystrophy locus on chromosome 7

We report the identification of a new locus for autosomal dominant limb-girdle muscular dystrophy (LGMD1) on 7q. Two of five families (1047 and 1701) demonstrate evidence in favor of linkage to this region. The maximum two-point LOD score for family 1047 was 3.76 for D7S427, and that for family 1701 was 2.63 for D7S3058. Flanking markers place the LGMD1 locus between D7S2423 and D7S427, with multipoint analysis slightly favoring the 9-cM interval spanned by D7S2546 and D7S2423. Three of five families appear to be unlinked to this new locus on chromosome 7, thus establishing further heterogeneity within the LGMD1 diagnostic classification.

Respiratory muscle involvement in Bethlem myopathy

We report a patient from a previously reported family with autosomal dominant Bethlem myopathy who demonstrated childhood onset, slowly progressive limb-girdle muscle weakness, contractures, and progressive respiratory compromise. Chest x-ray, pulmonary function tests, and electrophysiologic studies suggested respiratory muscle involvement, thus expanding the clinical repertoire in Bethlem myopathy.

Genetic mapping of a novel familial form of infantile hemangioma

Infantile hemangiomas are the most common tumor of infancy, occurring with an incidence of up to 10% of all births. They are benign but highly proliferative lesions involving aberrant localized growth of capillary endothelium. Although most hemangiomas occur sporadically and as single lesions, or in conjunction with pleiotropic genetic syndromes, we have previously identified six kindreds where hemangiomas appear to segregate as an autosomal dominant trait with high penetrance. Four such families contain affected individuals in three or more generations. In the current study, blood samples from five of these families were collected and used in a whole genome linkage search at 10-cM resolution. We established evidence for linkage to 5q in three families, and evidence for locus heterogeneity. The three 5q-linked families were further genotyped to generate haplotype information and narrow the candidate interval. Based on recombination breakpoint analysis, the interval exists between markers D5S2490 and D5S408, spanning 55 cM, and corresponding to 5q31-33. Using information from affected and unaffected individuals, the interval spans 38 cM between markers D5S1469 and D5S211. Three candidate genes involved with blood vessel growth map to this region: fibroblast growth factor receptor-4 (FGFR4), platelet-derived growth factor receptor-beta (PDG-FRB), and fms-related tyrosine kinase-4 (FLT4). The genes and gene products associated with familial hemangiomas may be involved somatically in the more common sporadic cases.

Use of a CEPH meiotic breakpoint panel to refine the locus of limb-girdle muscular dystrophy type 1A (LGMD1A) to a 2-Mb interval on 5q31

Limb-girdle muscular dystrophy type 1A (LGMD1A) is an autosomal dominant disease characterized by progressive weakness of the hip and shoulder girdle. The gene for LGMD1A had been localized to a 7-cM interval at 5q31 in a single large family (Family 39). To refine the localization of LGMD1A further and to aid in its identification, a high-resolution physical map of the locus was used to identify and provisionally localize 25 polymorphic markers. A subset of these markers was then ordered genetically, using a CEPH meiotic breakpoint panel, resulting in an integrated physical-genetic map of the locus. Relevant markers were genotyped on the members of Family 39 who contained informative recombination events, resulting in a further narrowing of LGMD1A to an interval bounded by D5S479 and D5S594, estimated to be 2 Mb in size. Integration of the genetic and physical map permits the identification of several transcription units from within the narrowed LGMD1A interval, including one that is muscle specific, representing candidate genes for this familial dystrophy.

A PHEX gene mutation is responsible for adult-onset vitamin D-resistant hypophosphatemic osteomalacia: evidence that the disorder is not a distinct entity from X-linked hypophosphatemic rickets

Previous investigators described a kindred with an X-linked dominant form of phosphate wasting in which affected children did not have radiographic evidence of rickets, whereas older individuals were progressively disabled by severe bowing. They proposed that this kindred suffered from a distinct disorder that they referred to as adult-onset vitamin D-resistant hypophosphatemic osteomalacia (AVDRR). We recently identified a gene, PHEX, that is responsible for the disorder X-linked hypophosphatemic rickets. To determine whether AVDRR is a distinct form of phosphate wasting, we searched for PHEX mutations in affected members of the original AVDRR kindred. We found that affected individuals have a missense mutation in PHEX exon 16 that results in an amino acid change from leucine to proline in residue 555. Clinical evaluation of individuals from this family indicates that some of these individuals display classic features of X-linked hypophosphatemic rickets, and we were unable to verify progressive bowing in adults. In light of the variability in the clinical spectrum of X-linked hypophosphatemic rickets and the presence of a PHEX mutation in affected members of this kindred, we conclude that there is only one form of X-linked dominant phosphate wasting.

Exclusion of identified LGMD1 loci from four dominant limb-girdle muscular dystrophy families

The limb-girdle muscular dystrophies are a clinically and genetically heterogeneous group of disorders. Recent linkage analyses and positional cloning studies have identified numerous loci responsible for the recessive and dominant forms, underscoring the inherent heterogeneity. In this report, we investigate four large autosomal dominant limb-girdle pedigrees and exclude these pedigrees from linkage to these loci. In addition, there is no evidence for linkage to any of the seven recessive LGMD loci.

CA repeat polymorphism of the COL6A3 gene on chromosome 2q37

A CA dinucleotide repeat has been identified in an intron of the human alpha3(VI) collagen gene (COL6A3) located on chromosome 2q37. Analysis of 100 chromosomes in unrelated Caucasians has demonstrated the existence of eight alleles, and the allelic frequencies have been determined.

Missense mutation in a von Willebrand factor type A domain of the alpha 3(VI) collagen gene (COL6A3) in a family with Bethlem myopathy

The Bethlem myopathy is a rare autosomal dominant proximal myopathy characterized by early childhood onset and joint contractures. Evidence for linkage and genetic heterogeneity has been established, with the majority of families linked to 21q22.3 and one large family linked to 2q37, implicating the three type VI collagen subunit genes, COL6A1 (chromosome 21), COL6A2 (chromosome 21) and COL6A3 (chromosome 2) as candidate genes. Mutations of the invariant glycine residues in the triple-helical domain-coding region of COL6A1 and COL6A2 have been reported previously in the chromosome 21-linked families. We report here the identification of a G-->A mutation in the N-terminal globular domain-coding region of COL6A3 in a large American pedigree (19 affected, 12 unaffected), leading to the substitution of glycine by glutamic acid in the N2 motif, which is homologous to the type A domains of the von Willebrand factor. This mutation segregated to all affected family members, to no unaffected family members, and was not identified in 338 unrelated Caucasian control chromosomes. Thus mutations in either the triple-helical domain or the globular domain of type VI collagen appear to cause Bethlem myopathy.

Promoter sequence, expression, and fine chromosomal mapping of the human gene (MLP) encoding the MARCKS-like protein: identification of neighboring and linked polymorphic loci for MLP and MACS and use in the evaluation of human neural tube defects

The MARCKS-like protein (MLP), also known as F52, MacMARCKS, or MARCKS-related protein, is a widely distributed substrate for protein kinase C (PKC). Recent studies using gene disruption in vivo have demonstrated the importance of both MARCKS and MLP to the development of the central nervous system; specifically, mice lacking either protein exhibit a high frequency of neural tube defects. We isolated a genomic clone for human MLP and discovered a directly linked polymorphism (MLP1) useful for genetic linkage analysis. The MLP promoter was 71% identical over 433 bp to that of the corresponding mouse gene, Mlp, with conservation of many putative transcription factor-binding sites; it was only 36% identical over 433 bp to the promoter of the human gene, MACS, which encodes the MLP homologue MARCKS. This 433-bp fragment drove expression of an MLP-beta-galactosidase transgene in a tissue-specific and developmental expression pattern that was similar to that observed for the endogenous gene, as shown by in situ hybridization histochemistry. In contrast to MACS, the MLP and Mlp promoters contain a TATA box approximately 40 bp 5' of the presumed transcription initiation site. MLP was localized to chromosome 1p34-->1pter by analysis of human-mouse somatic cell hybrid DNA and to 1p34 by fluorescence in situ hybridization. Radiation hybrid mapping of MLP placed it between genetic markers D1S511 (LOD > 3.0) and WI9232. MACS was localized to 6q21 between D6S266 (LOD > 3.0) and AFM268uh5 by the same technique. We tested the novel MLP1 polymorphism and the MACS flanking markers in a series of 43 Caucasian simplex families in which the affected child had a lumbosacral myelomeningocele. We found no evidence of linkage disequilibrium, suggesting that these loci were not major genes for spina bifida in these families. Nonetheless, the identification of linked and neighboring polymorphisms for MACS and MLP should permit similar genetic studies in other groups of patients with neural tube defects and other neurodevelopmental abnormalities.

Evidence for anticipation in autosomal dominant limb-girdle muscular dystrophy

Anticipation, an increase in severity or decrease in age of onset (AO) inherent in the transmission of the disease gene from affected parent to affected child, has been increasingly described in human disease. To assess anticipation in a large kindred in which autosomal dominant limb-girdle muscular dystrophy (LGMD1A) is segregating, age of disease onset was collected from patient interviews of affected family members. A total of 25 parent-offspring pairs, in which the parents are three (3R), four (4R), or five (5R) generations removed from a common founding ancestor, were available for analysis. Life table analyses showed significant decreases in age at first reported symptoms in the offspring of the 3R (chi2=5.55, p=0.02) and 4R (chi2=7.81, p=0.005) parents. Pairwise analyses confirmed this decrease with a median decrease of 13 years in transmission to offspring from 3R parents and 18 years in transmission to offspring from 4R parents. The finding of anticipation in this pedigree suggests that the mutation in LGMD1A may be the result of the expansion of an unstable trinucleotide repeat.

Using neural networks as an aid in the determination of disease status: comparison of clinical diagnosis to neural-network predictions in a pedigree with autosomal dominant limb-girdle muscular dystrophy

Studies of the genetics of certain inherited diseases require expertise in the determination of disease status even for single-locus traits. For example, in the diagnosis of autosomal dominant limb-girdle muscular dystrophy (LGMD1A), it is not always possible to make a clear-cut determination of disease, because of variability in the diagnostic criteria, age at onset, and differential presentation of disease. Mapping such diseases is greatly simplified if the data present a homogeneous genetic trait and if disease status can be reliably determined. Here, we present an approach to determination of disease status, using methods of artificial neural-network analysis. The method entails "training" an artificial neural network, with input facts (based on diagnostic criteria) and related results (based on disease diagnosis). The network contains weight factors connecting input "neurons" to output "neurons," and these connections are adjusted until the network can reliably produce the appropriate outputs for the given input facts. The trained network can be "tested" with a second set of facts, in which the outcomes are known but not provided to the network, to see how well the training has worked. The method was applied to members of a pedigree with LGMD1A, now mapped to chromosome 5q. We used diagnostic criteria and disease status to train a neural network to classify individuals as "affected" or "not affected." The trained network reproduced the disease diagnosis of all individuals of known phenotype, with 98% reliability. This approach defined an appropriate choice of clinical factors for determination of disease status. Additionally, it provided insight into disease classification of those considered to have an "unknown" phenotype on the basis of standard clinical diagnostic methods.

Quantitative DNA pooling to increase the efficiency of linkage analysis in autosomal dominant disease

DNA pooling is an efficient method to rapidly perform genome-wide linkage scans in autosomal recessive diseases in inbred populations where affected individuals are likely to be homozygous for alleles near the disease gene locus. We wanted to examine whether this approach would detect linkage in autosomal dominant (AD) disorders where affected individuals may share one allele identical by descent at loci tightly linked to the disease. Two large outbred pedigrees in which the AD diseases familial venous malformation (FVM) and hereditary hemorrhagic telangiectasia (HHT1), linked to 9p and 9q, respectively, were investigated. Separate pools of DNA from affected (n = 21 for FVM and 17 for HHT1) and unaffected family members (n = 9 FVM and HHT1), and 25 unrelated population controls were established. Polymorphic markers spanning chromosome 9 at approximately 13.5-cM intervals were amplified using standard PCR. Allele quantitation was performed with a fluorimager. Visual inspection of allele intensities and frequency distributions suggested a shift in frequency of the most common allele in the affecteds lane when compared to control lanes for markers within 30 cM of the FVM and HHT1 loci. These subjective assessments were confirmed statistically by testing for the difference between two proportions (one-sided; P < or = 0.05). When using population controls, the true-positive rates for FVM and HHT1 were 5/5 and 2/5 markers, respectively. False-positive rates for FVM and HHT1 were 3/9 and 2/9, respectively. In both AD diseases investigated, quantitative DNA pooling detected shifts in allele frequency, thus identifying areas of known linkage in most cases. The utility of this technique depends on the size of the pedigree, frequency of the disease-associated allele in the population, and the choice of appropriate controls. Although the false-positive rate appears to be high, this approach still serves to reduce the amount of overall genotyping by about 60%. DNA pooling merits further investigation as a potential strategy in increasing the efficiency of genomic linkage scans.

Autosomal dominant hypophosphatemic rickets is linked to chromosome 12p13

Autosomal dominant hypophosphatemic rickets (ADHR) is an inherited disorder of isolated renal phosphate wasting, the pathogenesis of which is unknown. We performed a genome-wide linkage study in a large kindred to determine the chromosome location of the ADHR gene. Two-point LOD scores indicate that the gene is linked to the markers D12S314 [Z(theta) = 3.15 at theta = 0.0], vWf [Z(theta) = 5.32 at theta = 0.0], and CD4 [Z(theta) = 3.53 at theta = 0.0]. Moreover, multilocus analysis indicates that the ADHR gene locus is located on chromosome 12p13 in the 18-cM interval between the flanking markers D12S100 and D12S397. These data are the first to establish a chromosomal location for the ADHR locus and to provide a framework map to further localize the gene. Such studies will permit ultimate identification of the ADHR gene and provide further insight into phosphate homeostasis.

Linkage of familial dilated cardiomyopathy with conduction defect and muscular dystrophy to chromosome 6q23

Inherited cardiomyopathies may arise from mutations in genes that are normally expressed in both heart and skeletal muscle and therefore may be accompanied by skeletal muscle weakness. Phenotypically, patients with familial dilated cardiomyopathy (FDC) show enlargement of all four chambers of the heart and develop symptoms of congestive heart failure. Inherited cardiomyopathies may also be accompanied by cardiac conduction-system defects that affect the atrioventricular node, resulting in bradycardia. Several different chromosomal regions have been linked with the development of autosomal dominant FDC, but the gene defects in these disorders remain unknown. We now characterize an autosomal dominant disorder involving dilated cardiomyopathy, cardiac conduction-system disease, and adult-onset limb-girdle muscular dystrophy (FDC, conduction disease, and myopathy [FDC-CDM]). Genetic linkage was used to exclude regions of the genome known to be linked to dilated cardiomyopathy and muscular dystrophy phenotypes and to confirm genetic heterogeneity of these disorders. A genomewide scan identified a region on the long arm of chromosome 6 that is significantly associated with the presence of myopathy (D6S262; maximum LOD score [Z(max)] 4.99 at maximum recombination fraction [theta(max)] .00), identifying FDC-CDM as a genetically distinct disease. Haplotype analysis refined the interval containing the genetic defect, to a 3-cM interval between D6S1705 and D6S1656. This haplotype analysis excludes a number of striated muscle-expressed genes present in this region, including laminin alpha2, laminin alpha4, triadin, and phospholamban.

The thermolabile variant of methylenetetrahydrofolate reductase (MTHFR) is not a major risk factor for neural tube defect in American Caucasians. The NTD Collaborative Group

Mutations in the gene for methylenetetrahydrofolate reductase (MTHFR) have been implicated as a risk factor in the formation of neural tube defects. We investigated this gene in a series of 65 sporadic American Caucasian patients with lumbosacral NTD and their unaffected parents, using both case-control design and assessment of linkage disequilibrium. We found no evidence to support mutations in MTHFR as a significant risk factor for NTD in this population.

Further exclusion of FSHD1B from the telomeric region of 10q

The localization of the gene for the majority of cases of facioscapulohumeral muscular dystrophy is established as 4q35-qter (FSHD1A), although locus heterogeneity has been demonstrated with a minority of families unlinked to 4q. In FSHD1A, the disease is associated with a deletion of 3.3 kb repeats from a tandem repeat located near the as-yet-unidentified gene. This repeat cross-hybridizes with a telomeric region on 10q, making this cross-hybridizing region a feasible candidate gene for FSHD1B. We have tested the most telomeric marker on 10q (sAVA4) and excluded approximately 17 cM on either side of this marker as harboring the FSHD1B gene.

Confirmation of a second locus for CMT2 and evidence for additional genetic heterogeneity

The Charcot-Marie-Tooth (CMT) neuropathies are a group of disorders exhibiting neurophysical, pathological and genetic heterogeneity. CMT2 is a diagnostic subtype of this group of disorders characterized by variable expression and age-of-onset and normal or slightly diminished nerve conduction velocities. Previously, linkage and heterogeneity had been reported in CMT2 with linked families localizing to chromosome 1p (CMT2A). Recently a second CMT2 locus has been described on chromosome 7 in a single large CMT2 family (CMT2D). We have performed pedigree linkage analysis on 15 CMT2 families (N = 371 individuals, 106 affected family members) and have confirmed linkage to chromosome 7. Furthermore, using both admixture and multipoint linkage analysis we show conclusive evidence for additional heterogeneity within this clinical subtype with evidence of families that exclude linkage to both the CMT2D and CMT2A regions. In addition, unlike the previous report we found no obvious consistent clinical differences between the linked family types.

False positive rates in a genomic screen for complex quantitative traits

We conducted a genomic screen for genes associated with Q1, Q2, and Q3 in 239 nuclear pedigrees from replicate 115, Problem Set 2A. We compared false positive (FP) and true positive (TP) rates for three significance levels and two map densities. Using the 2 cM genetic map and alpha = 0.05 produced the most FP but detected the greatest number of major genes. Following up only 31 plateaus (two or more adjacent markers with significant results) from the 2 cM screen eliminated some FP, but failed to detect MG3 for Q3. Multipoint analysis reduced the number of priority regions from 31 to seven; only two of these regions were TP. Replication of the two-point analysis of plateau markers in replicate 80 detected all of the genes associated with Q1 and Q2, but not Q3. Multipoint analysis in replicate 80 failed to replicate any genes associated with Q1, Q2, or Q3, but "replicated" two FP regions. While FP may be reduced by decreasing map density, considering only plateaus for follow up and decreasing significance levels, such adjustments may also fail to detect weak TP. Multipoint analysis and replication in independent data sets may not be reliable methods of distinguishing FP from TP.

Type VI collagen mutations in Bethlem myopathy, an autosomal dominant myopathy with contractures

Among the diverse family of collagens, the widely expressed microfibrillar type VI collagen is believed to play a role in bridging cells with the extracellular matrix. Several observations imply substrate properties for cell attachment as well as association with major collagen fibers. Previously, we have established genetic linkage between the genes encoding the three constituent alpha-chains of type VI collagen and Bethlem myopathy. A distinctive feature of this autosomal dominant disorder consists of contractures of multiple joints in addition to generalized muscular weakness and wasting. Nine kindreds show genetic linkage to the COL6A1-COL6A2 cluster on chromosome 21q22.3 (refs 3,4; manuscript submitted) whereas one family shows linkage to markers on chromosome 2q37 close to COL6A3 (ref. 5). Sequence analysis in four families reveals a mutation in COL6A1 in one and a COL6A2 mutation in two other kindreds. Both mutations disrupt the Gly-X-Y motif of the triple helical domain by substitution of Gly for either Val or Ser. Analogous to the putative perturbation of the anchoring function of the dystrophin-associated complex in congenital muscular dystrophy with mutations in the alpha 2-subunit of laminin, our observations suggest a similar mechanism in Bethlem myopathy.

Evidence for locus heterogeneity in the Bethlem myopathy and linkage to 2q37

The Bethlem myopathy, a childhood onset autosomal dominant myopathy with joint contractures, has recently been localized to 21q in a series of Dutch families and the alpha 1 and alpha 2 subunits of type VI collagen (COL6A1 and COL6A2) have been postulated as candidate genes. We investigate a large family of French Canadian descent (family 1489) in which the Bethlem myopathy is segregating. Family 1489 is unlinked to the region of interest on 21q, thus demonstrating locus heterogeneity within the Bethlem myopathy classification. In view of the localization of the genes coding the alpha 1 and alpha 2 subunits of type VI collagen on chromosome 21q, we carried out linkage analysis on chromosome 2q where the alpha 3 subunit of type VI collagen has been localized. We demonstrate linkage to markers in this region, define the region of disease gene localization, and confirm by FISH analysis that COL6A3 is located within the interval of interest making COL6A3 a feasible candidate gene for the Bethlem myopathy.