Association of the ROBO1 gene with reading disabilities in a family-based analysis.

Linkage studies have identified a locus on chromosome 3 as reading disabilities (RD) and speech and sound disorder (SSD) susceptibility region, with both RD and SSD sharing similar phonological processing and phonological memory difficulties. One gene in this region, roundabout homolog 1 (ROBO1), has been indicated as a RD candidate and has shown significant association with measures of phonological memory in a population-based sample. In this study, we conducted a family-based association analysis using two independent samples collected in Toronto and Calgary, Canada. Using the two samples, we tested for association between ROBO1 single nucleotide polymorphisms (SNPs) and RD, along with quantitative measures for reading, spelling and phonological memory. One SNP, rs331142, which was selected based on its correlation with ROBO1 expression in brain tissue, was found to be significantly associated with RD in the Toronto sample with over transmission of the minor C allele (P = 0.001), correlated with low expression. This SNP is located ~200 bp from a putative enhancer and results for a marker within the enhancer, rs12495133, showed evidence for association with the same allele in both the Toronto and Calgary samples (P = 0.005 and P = 0.007). These results support previous associations between ROBO1 and RD, as well as correlation with low gene expression, suggesting a possible mechanism of risk conferred by this gene.

A family-based association analysis and meta-analysis of the reading disabilities candidate gene DYX1C1.

Reading disabilities (RD) have a significant genetic basis and have shown linkage to multiple regions including chromosome 15q. Dyslexia susceptibility 1 candidate gene 1 (DYX1C1) on chromosome 15q21 was originally proposed as a candidate gene with two potentially functional polymorphisms at the -3G/A and 1249G/T positions showing association with RD. However, subsequent studies have yielded mixed results. We performed a literature review and meta-analysis of the -3G/A and 1249G/T polymorphisms, including new unpublished data from two family-based samples. Ten markers in DYX1C1 were genotyped in the two independently ascertained samples. Single marker and -3G/A:1249G/T haplotype analyses were performed for RD in both samples, and quantitative trait analyses using standardized reading-related measures was performed in one of the samples. For the meta-analysis, we used a random-effects model to summarize studies that tested for association between -3G/A or 1249G/T and RD. No significant association was found between the DYX1C1 SNPs and RD or any of the reading-related measures tested after correction for the number of tests performed. The previously reported risk haplotype (-3A:1249T) was not biased in transmission. A total of 9 and 10 study samples were included in the meta-analysis of the -3G/A and 1249G/T polymorphisms, respectively. Neither polymorphism reached statistical significance, but the heterogeneity for the 1249G/T polymorphism was high. The results of this study do not provide evidence for association between the putatively functional SNPs -3G/A and 1249G/T and RD.

Dense-map genome scan for dyslexia supports loci at 4q13, 16p12, 17q22; suggests novel locus at 7q36.

Analysis of genetic linkage to dyslexia was performed using 133,165 array-based SNPs genotyped in 718 persons from 101 dyslexia-affected families. Results showed five linkage peaks with lod scores >2.3 (4q13.1, 7q36.1-q36.2, 7q36.3, 16p12.1, and 17q22). Of these five regions, three have been previously implicated in dyslexia (4q13.1, 16p12.1, and 17q22), three have been implicated in attention-deficit hyperactivity disorder (ADHD, which highly co-occurs with dyslexia; 4q13.1, 7q36.3, 16p12.1) and four have been implicated in autism (a condition characterized by language deficits; 7q36.1-q36.2, 7q36.3, 16p12.1, and 17q22). These results highlight the reproducibility of dyslexia linkage signals, even without formally significant lod scores, and suggest dyslexia predisposing genes with relatively major effects and locus heterogeneity. The largest lod score (2.80) occurred at 17q22 within the MSI2 gene, involved in neuronal stem cell lineage proliferation. Interestingly, the 4q13.1 linkage peak (lod 2.34) occurred immediately upstream of the LPHN3 gene, recently reported both linked and associated with ADHD. Separate analyses of larger pedigrees revealed lods >2.3 at 1-3 regions per family; one family showed strong linkage (lod 2.9) to a known dyslexia locus (18p11) not detected in our overall data, demonstrating the value of analyzing single large pedigrees. Association analysis identified no SNPs with genome-wide significance, although a borderline significant SNP (P = 6 × 10(-7)) occurred at 5q35.1 near FGF18, involved in laminar positioning of cortical neurons during development. We conclude that dyslexia genes with relatively major effects exist, are detectable by linkage analysis despite genetic heterogeneity, and show substantial overlapping predisposition with ADHD and autism.

Association of AXIN2 with non-syndromic oral clefts in multiple populations.

We have previously shown the association of AXIN2 with oral clefts in a US population. Here, we expanded our study to explore the association of 11 AXIN2 markers in 682 cleft families from multiple populations. Alleles for each AXIN2 marker were tested for transmission distortion with clefts by means of the Family-based Association Test. We observed an association with SNP rs7224837 and all clefts in the combined populations (p = 0.001), and with SNP rs3923086 and cleft lip and palate in Asian populations (p = 0.004). We confirmed our association findings in an additional 528 cleft families from the United States (p < 0.009). We tested for gene-gene interaction between AXIN2 and additional cleft susceptibility loci. We assessed and detected Axin2 mRNA and protein expression during murine palatogenesis. In addition, we also observed co-localization of Axin2 with Irf6 proteins, particularly in the epithelium. Our results continue to support a role for AXIN2 in the etiology of human clefting. Additional studies should be performed to improve our understanding of the biological mechanisms linking AXIN2 to oral clefts.

Murine pancreatic beta TC3 cells show greater 2', 5'-oligoadenylate synthetase (2'5'AS) antiviral enzyme activity and apoptosis following IFN-alpha or poly(I:C) treatment than pancreatic alpha TC3 cells.

Type 1 diabetes is caused by autoimmune destruction of pancreatic beta cells, possibly virus initiated. Virus infection induces alpha-interferon (IFN-alpha), leading to upregulation of genes encoding double-stranded (ds) RNA-dependent antiviral enzymes 2', 5'-oligoadenylate synthetase (2'5'AS) and PKR (p68). To investigate whether beta cell specificity could be due to antiviral differences between beta and alpha cells, we treated beta and alpha TC3 cell lines with IFN-alpha and/or poly(I:C) (a synthetic dsRNA). Results showed that, following IFN-alpha stimulation, increases in 2'5'AS levels and activities were significantly higher in beta than alpha cells (P < .001), whereas increases in PKR level and activity were comparable in the two cell types. Poly(I:C) stimulated 2'5'AS activity in beta but not alpha cells, and co-transfection IFN-alpha plus poly(I:C) induced apoptosis in beta but not alpha cells. These findings suggest that the elevated 2'5'AS response of pancreatic beta cells could render them particularly vulnerable to damage and/or apoptosis during virus infection.

A specific requirement for PDGF-C in palate formation and PDGFR-alpha signaling.

PDGF-C is a member of the platelet-derived growth factor (PDGF) family, which signals through PDGF receptor (PDGFR) alphaalpha and alphabeta dimers. Here we show that Pdgfc(-/-) mice die in the perinatal period owing to feeding and respiratory difficulties associated with a complete cleft of the secondary palate. This phenotype was less severe than that of Pdgfra(-/-) embryos. Pdgfc(-/-) Pdgfa(-/-) embryos developed a cleft face, subepidermal blistering, deficiency of renal cortex mesenchyme, spina bifida and skeletal and vascular defects. Complete loss of function of both ligands, therefore, phenocopied the loss of PDGFR-alpha function, suggesting that both PDGF-A and PDGF-C signal through PDGFR-alpha to regulate the development of craniofacial structures, the neural tube and mesodermal organs. Our results also show that PDGF-C signaling is a new pathway in palatogenesis, different from, and independent of, those previously implicated.

TP63 mutation and clefting modifier genes in an EEC syndrome family.

Autosomal dominant EEC syndrome consists of ectrodactyly, ectodermal dysplasia, and cleft lip with or without cleft palate (CL/P). We investigated an EEC kindred with 10 affected persons in three generations in order to map the causative mutation in this family and to map modifier genes that contribute to the expression of facial clefting in the phenotype. DNA from 15 family members was genotyped for 388 genome screen markers. Analysis revealed maximal linkage between EEC and chromosome 3q27, which contains a known EEC gene - tumor protein 63 (TP63). Sequencing showed a CGT-->TGT missense mutation (R280C) in exon 7, previously reported to cause EEC in four families, and ectrodactyly alone (split hand-foot malformation) in one sporadic case and one large kindred. Analysis of the clefting phenotype in this EEC family demonstrated maximal linkage to two regions on chromosomes 4q and 14, which multiple studies have implicated in non-syndromic CL/P. In conclusion, this study demonstrates that the mutation of TP63 is the major (Mendelian) EEC gene in this kindred and suggests that additional minor modifying genes which predispose to non-syndromic CL/P could also contribute to the expression of the clefting component of the syndrome in this family.

Evidence for a locus (IDDM16) in the immunoglobulin heavy chain region on chromosome 14q32.3 producing susceptibility to type 1 diabetes.

Type 1 diabetes results from autoimmune destruction of pancreatic islet beta-cells, possibly initiated or exacerbated by viral infections. Recent studies have demonstrated that antibodies towards enterovirus and autoantibodies towards islet cell components develop in the long preclinical phase of type 1 diabetes. We therefore hypothesised that susceptibility to type 1 diabetes could be influenced by genetic factors controlling production of antiviral antibodies or autoantibodies or both. To search for evidence of linkage or association (linkage disequilibrium) between type 1 diabetes and the immunoglobulin heavy chain (IGH) region, 351 North American and British families with > or =2 diabetic children were genotyped for IGH region microsatellites. Using affected sibpair analysis, significant evidence for linkage was obtained for three markers close to the IGH gene cluster (P values 0.004, 0.002, 0.002). No evidence was found for association using family-based methods. To attempt to confirm these findings, a smaller dataset (241 families, 138 with > or =2 diabetic children) from Denmark, a more genetically-homogeneous population, was genotyped for one marker only. These families showed no linkage, but significant evidence for association (P = 0.019). This study suggests that a locus (assigned the symbol IDDM16) in the IGH region, possibly an IGH gene, influences susceptibility to type 1 diabetes.

Genetic linkage and association studies of Type I diabetes: challenges and rewards.

Family and twin studies have shown clearly that Type I (insulin-dependent) diabetes mellitus has a genetic basis. However, only within the past decade has it been possible to systematically attempt to identify the genes that increase susceptibility to this disorder using linkage and association analysis of genetic markers distributed across the genome. More than 20 putative diabetes-predisposing genes have been localised in addition to HLA region susceptibility genes detected more than 25 years ago. Unfortunately, the effects of most diabetes-predisposing genes are weak, with the exception of HLA region susceptibility genes (which contribute about half of the genetic risk). The overall effects of diabetes-predisposing genes could be weak because the susceptibility gene occurs in only a small proportion of diabetic patients or the susceptibility gene (although it might be common) produces only a modest increase in risk, probably by acting in concert with other such genes to cause disease. The weak effects of these genes have made them difficult to locate, difficult to confirm in independent studies and difficult to isolate by genetic procedures. This paper summarizes the major challenges that have faced geneticists in mapping Type I diabetes genes, and reviews the progress achieved to date. The rewards of the genetic studies will be twofold: increased understanding of the causes of Type I diabetes, facilitating creation of preventative therapies, and enabling clinicians in the future to use genetic information to predict which children are predisposed to diabetes in order to target them for preventative therapies.

Grade-of-membership sibpair linkage analysis maps IDDM11 to chromosome 14q24.3-q31.

We demonstrate the use of Grade-of-membership (GoM) (Manton et al. 1994) for sibpair linkage analysis: GoM was used to map the IDDM11 locus to the region of chromosome 14q24.3 identified by Field et al. (1996). Haplotype groups were constructed from sib pair information on the number of shared alleles. The sample consisted of 578 sibling pairs found in 246 multiplex IDDM families. Both siblings were diabetic in 53% of the pairs (AA). Pair members could share 0, 1 or 2 alleles IBS at each of eight linked marker loci spanning IDDM11. Three model-based groups best represented the data on allele sharing: the groups corresponded to 'No', 'One' and 'Two' shared haplotypes for the region. Group 'Two' was larger (37% vs. 25%, p < 0.0001) and more homogeneous (p < 0.0001) than expected by chance consistent with the IDDM11 locus being a determinant of diabetes in multiplex families. Genetic linkage of IDDM to the region was demonstrated by a 19% increase in the proportion of AA pairs over the haplotype groups: 'No', 42%; 'One', 49%; 'Two', 61%, p = 0.0005, representing a 43% relative increase.

Recombination in the pseudoautosomal region in a 47,XYY male.

Males with a 47,XYY karyotype generally have chromosomally normal children, despite the high theoretical risk of aneuploidy. Studies of sperm karyotypes or FISH analysis of sperm have demonstrated that the majority of sperm are chromosomally normal in 47,XYY men. There have been a number of meiotic studies of XYY males attempting to determine whether the additional Y chromosome is eliminated during spermatogenesis, with conflicting results regarding the pairing of the sex chromosomes and the presence of an additional Y. We analyzed recombination in the pseudoautosomal region of the XY bivalent to determine whether this is perturbed in a 47,XYY male. A recombination frequency similar to normal 46,XY men would indicate normal pairing within the XY bivalent, whereas a significantly altered frequency would suggest other types of pairing such as a YY bivalent or an XYY trivalent. Two DNA markers, STS/STS pseudogene and DXYS15, were typed in sperm from a heterozygous 47,XYY male. Individual sperm (23,X or Y) were isolated into PCR tubes using a FACStarPlus flow cytometer. Hemi-nested PCR analysis of the two DNA markers was performed to determine the frequency of recombination. A total of 108 sperm was typed with a 38% recombination frequency between the two DNA markers. This is very similar to the frequency of 38.3% that we have observed in 329 sperm from a normal 46,XY male. Thus our results suggest that XY pairing and recombination occur normally in this 47,XYY male. This could occur by the production of an XY bivalent and Y univalent (which is then lost in most cells) or by loss of the additional Y chromosome in some primitive germ cells or spermatogonia and a proliferative advantage of the normal XY cells.

Evidence for a susceptibility locus on chromosome 6q influencing phonological coding dyslexia.

A linkage study of 96 dyslexia families containing at least two affected siblings (totaling 877 individuals) has found evidence for a dyslexia susceptibility gene on chromosome 6q11.2-q12 (assigned the name DYX4). Using a qualitative phonological coding dyslexia (PCD) phenotype (affected, unaffected, or uncertain diagnoses), two-point parametric analyses found highly suggestive evidence for linkage between PCD and markers D6S254, D6S965, D6S280, and D6S251 (LOD(max) scores = 2.4 to 2.8) across an 11 cM region. Multipoint parametric analysis supported linkage of PCD to this region (peak HLOD = 1.6), as did multipoint nonparametric linkage analysis (P = 0.012). Quantitative trait linkage analyses of four reading measures (phonological awareness, phonological coding, spelling, and rapid automatized naming speed) also provided evidence for a dyslexia susceptibility locus on chromosome 6q. Using a variance-component approach, analysis of phonological coding and spelling measures resulted in peak LOD scores at D6S965 of 2.1 and 3.3, respectively, under 2 degrees of freedom. Furthermore, multipoint nonparametric quantitative trait sibpair analyses suggested linkage between the 6q region and phonological awareness, phonological coding, and spelling (P = 0.018, 0.017, 0.0005, respectively, for unweighted sibpairs < 18 years of age). Although conventional significance thresholds were not reached in the linkage analyses, the chromosome 6q11.2-q12 region clearly warrants investigation in other dyslexia family samples to attempt replication and confirmation of a dyslexia susceptibility gene in this region.

Single sperm typing demonstrates that reduced recombination is associated with the production of aneuploid 24,XY human sperm.

To account for the increased proportion of paternal nondisjunction in 47,XXY males as compared to other trisomies, it has been suggested that the XY bivalent, with its reduced region of homology, is particularly susceptible to nondisjunction. Molecular studies of liveborn Klinefelter syndrome (47,XXY) individuals have reported an association between the absence of recombination in the pseudoautosomal region and nondisjunction of the XY bivalent. In this study we examined single sperm from a normal 46,XY male to determine if there is any alteration in the recombination frequency of aneuploid disomic 24,XY sperm compared to unisomic sperm (23,X or Y). Two DNA markers STS/STS pseudogene and DXYS15 were typed in sperm from a heterozygous man to determine if recombination had occurred in the pseudoautosomal region. Individual unisomic sperm (23,X or Y) were isolated using a FACStar(Plus) flow cytometer into PCR tubes. To identify disomic 24,XY sperm, 3-colour FISH analysis was performed with probes for chromosomes X,Y and 1. The 24,XY cells were identified using fluorescence microscopy, each disomic sperm was scraped off the slide using a glass needle attached to a micromanipulator and then put into a PCR tube. Hemi-nested PCR analysis of the two markers was performed to determine the frequency of recombination. A total of 329 unisomic sperm and 150 disomic sperm have been typed. The frequency of recombination between the two DNA markers was 38.3% for the unisomic sperm, similar to frequencies previously reported. The 24,XY disomic sperm had an estimated recombination frequency of 25.3%, however, a highly significant decrease compared to the unisomic 23,X or 23,Y sperm (chi(2) = 10.7, P = 0.001). This direct analysis of human sperm indicates that lack of recombination in the pseudoautosomal region is a significant cause of XY nondisjunction and thus Klinefelter syndrome.

Absence of significant linkage between phonological coding dyslexia and chromosome 6p23-21.3, as determined by use of quantitative-trait methods: confirmation of qualitative analyses.

We recently reported the absence of significant linkage of phonological coding dyslexia (PCD) to chromosome 6p23-p21.3 in 79 families with at least two affected siblings, even though linkage of dyslexia to this region has been found in four other independent studies. Whereas, in our previous analyses, we used a qualitative (affected, unaffected, or uncertain) PCD phenotype, here we report a reanalysis of linkage to the chromosome 6p region, by use of four quantitative measures of reading disability: phonological awareness, phonological coding, spelling, and rapid-automatized-naming (RAN) speed. The phonological-coding and spelling measures were highly correlated with each other and with the qualitative PCD phenotype, whereas the phonological-awareness and RAN-speed measures were only moderately correlated with the other measures. Using two-point and multipoint quantitative-trait sib-pair linkage analyses and variance-components analyses, we were unable to detect significant evidence for a locus in the 6p23-p21.3 region influencing any of the quantitative reading measures, supporting our previous qualitative linkage results. The most likely explanation for our inability to detect linkage between dyslexia and this region is that families with subtypes of dyslexia linked to this region are underrepresented in our sample, because of either chance or varying ascertainment criteria.

Absence of linkage of phonological coding dyslexia to chromosome 6p23-p21.3 in a large family data set.

Previous studies have suggested that a locus predisposing to specific reading disability (dyslexia) resides on chromosome 6p23-p21.3. We investigated 79 families having at least two siblings affected with phonological coding dyslexia, the most common form of reading disability (617 people genotyped, 294 affected), and we tested for linkage with the genetic markers reported to be linked to dyslexia in those studies. No evidence for linkage was found by LOD score analysis or affected-sib-pair methods. However, using the affected-pedigree-member (APM) method, we detected significant evidence for linkage and/or association with some markers when we used published allele frequencies with weighting of rarer alleles. APM results were not significant when we used marker allele frequencies estimated from parents. Furthermore, results were not significant with the more robust SIMIBD method using either published or parental marker frequencies. Finally, family-based association analysis using the AFBAC program showed no evidence for association with any marker. We conclude that the APM method should be used only with extreme caution, because it appears to have generated false-positive results. In summary, using a large data set with high power to detect linkage, we were unable to find evidence for linkage or association between phonological coding dyslexia and chromosome 6p markers.

Specificity of eating disorders diagnoses in families of probands with anorexia nervosa and bulimia nervosa.

This study examines the prevalence of anorexia nervosa and bulimia nervosa in relatives of probands, and examines the probandwise specificity of any familial clustering. Data were collected from probands using the family history method. Probands were recruited in a sequential cohort fashion. Information collected from probands was rated semiblindly by two of the authors, and a diagnostic hierarchy applied to arrive at a diagnosis for each of the relatives assessed. Data are reported on 2,125 family members, collected from 93 probands. Diagnostic agreement between raters was high, with serious disagreement present in three of 167 possible cases of an eating disorder. Rates of anorexia nervosa, bulimia nervosa, major depression, and substance abuse declined from first- to third-degree relatives, which is consistent with genetic clustering, and there was evidence of a cohort effect operating for anorexia nervosa and bulimia nervosa. The rates of anorexia nervosa and bulimia nervosa in all family members were 5.1%, and 4.3% respectively. An analysis of maternal and paternal descent showed no evidence for X-linked dominant transmission in these families. Preliminary analysis of the clustering of diagnoses in relatives showed a tendency (chi 2 = 14.47, P = .006) for family members to be affected by the same diagnosis as was the proband. This trend was strongest for anorexia nervosa, but there was overlap when the proband had a lifetime diagnosis of bulimia nervosa, with or without anorexia nervosa. These results are compatible with the existence of genetic factors influencing predisposition to eating disorders, but do not prove such.

Localization of a gene for incomplete X-linked congenital stationary night blindness to the interval between DXS6849 and DXS8023 in Xp11.23.

Congenital stationary night blindness (CSNB) is a nonprogressive retinal disorder characterized by night blindness, nystagmus, myopia, a variable decrease in visual acuity, an abnormal electroretinographic response, and a disturbance in dark adaptation. Two forms of X-linked CSNB have been defined, complete CSNB in which rod function is extinguished, and incomplete CSNB in which rod function is reduced but not extinguished, as seen by electroretinography and dark adaptometry. In studying a large family of Mennonite ancestry, we have confirmed linkage between the locus (CSNB2) for incomplete CSNB and genetic markers in the Xp11 region. In particular, lod scores of 12.25 and 15.26 at zero recombination were observed between CSNB2 and the markers DXS573 and DXS255. Detailed analysis of critical recombinant chromosomes in this extended family have refined the minimal region for the CSNB2 locus to the interval between DXS6849 and DXS8023 in Xp11.23.

Association mapping of disease loci, by use of a pooled DNA genomic screen.

Genomic screening to map disease loci by association requires automation, pooling of DNA samples, and 3,000-6,000 highly polymorphic, evenly spaced microsatellite markers. Case-control samples can be used in an initial screen, followed by family-based data to confirm marker associations. Association mapping is relevant to genetic studies of complex diseases in which linkage analysis may be less effective and to cases in which multigenerational data are difficult to obtain, including rare or late-onset conditions and infectious diseases. The method can also be used effectively to follow up and confirm regions identified in linkage studies or to investigate candidate disease loci. Study designs can incorporate disease heterogeneity and interaction effects by appropriate subdivision of samples before screening. Here we report use of pooled DNA amplifications-the accurate determination of marker-disease associations for both case-control and nuclear family-based data-including application of correction methods for stutter artifact and preferential amplification. These issues, combined with a discussion of both statistical power and experimental design to define the necessary requirements for detecting of disease loci while virtually eliminating false positives, suggest the feasibility and efficiency of association mapping using pooled DNA screening.