Increased variability of stimulus-driven cortical responses is associated with genetic variability in children with and without dyslexia

Individuals with dyslexia exhibit increased brainstem variability in response to sound. It is unknown as to whether increased variability extends to neocortical regions associated with audition and reading, extends to visual stimuli, and whether increased variability characterizes all children with dyslexia or, instead, a specific subset of children. We evaluated the consistency of stimulus-evoked neural responses in children with (N = 20) or without dyslexia (N = 12) as measured by magnetoencephalography (MEG). Approximately half of the children with dyslexia had significantly higher levels of variability in cortical responses to both auditory and visual stimuli in multiple nodes of the reading network. There was a significant and positive relationship between the number of risk alleles at rs6935076 in the dyslexia-susceptibility gene KIAA0319 and the degree of neural variability in primary auditory cortex across all participants. This gene has been linked with neural variability in rodents and in typical readers. These findings indicate that unstable representations of auditory and visual stimuli in auditory and other reading-related neocortical regions are present in a subset of children with dyslexia and support the link between the gene KIAA0319 and the auditory neural variability across children with or without dyslexia.

Multipoint genome-wide linkage scan for nonword repetition in a multigenerational family further supports chromosome 13q as a locus for verbal trait disorders

Verbal trait disorders encompass a wide range of conditions and are marked by deficits in five domains that impair a person's ability to communicate: speech, language, reading, spelling, and writing. Nonword repetition is a robust endophenotype for verbal trait disorders that is sensitive to cognitive processes critical to verbal development, including auditory processing, phonological working memory, and motor planning and programming. In the present study, we present a six-generation extended pedigree with a history of verbal trait disorders. Using genome-wide multipoint variance component linkage analysis of nonword repetition, we identified a region spanning chromosome 13q14-q21 with LOD = 4.45 between 52 and 55 cM, spanning approximately 5.5 Mb on chromosome 13. This region overlaps with SLI3, a locus implicated in reading disability in families with a history of specific language impairment. Our study of a large multigenerational family with verbal trait disorders further implicates the SLI3 region in verbal trait disorders. Future studies will further refine the specific causal genetic factors in this locus on chromosome 13q that contribute to language traits.

The DYX2 locus and neurochemical signaling genes contribute to speech sound disorder and related neurocognitive domains

A major milestone of child development is the acquisition and use of speech and language. Communication disorders, including speech sound disorder (SSD), can impair a child's academic, social and behavioral development. Speech sound disorder is a complex, polygenic trait with a substantial genetic component. However, specific genes that contribute to SSD remain largely unknown. To identify associated genes, we assessed the association of the DYX2 dyslexia risk locus and markers in neurochemical signaling genes (e.g., nicotinic and dopaminergic) with SSD and related endophenotypes. We first performed separate primary associations in two independent samples - Cleveland SSD (210 affected and 257 unaffected individuals in 127 families) and Denver SSD (113 affected individuals and 106 unaffected individuals in 85 families) - and then combined results by meta-analysis. DYX2 markers, specifically those in the 3' untranslated region of DCDC2 (P = 1.43 × 10(-4) ), showed the strongest associations with phonological awareness. We also observed suggestive associations of dopaminergic-related genes ANKK1 (P = 1.02 × 10(-2) ) and DRD2 (P = 9.22 × 10(-3) ) and nicotinic-related genes CHRNA3 (P = 2.51 × 10(-3) ) and BDNF (P = 8.14 × 10(-3) ) with case-control status and articulation. Our results further implicate variation in putative regulatory regions in the DYX2 locus, particularly in DCDC2, influencing language and cognitive traits. The results also support previous studies implicating variation in dopaminergic and nicotinic neural signaling influencing human communication and cognitive development. Our findings expand the literature showing genetic factors (e.g., DYX2) contributing to multiple related, yet distinct neurocognitive domains (e.g., dyslexia, language impairment, and SSD). How these factors interactively yield different neurocognitive and language-related outcomes remains to be elucidated.

Genome-wide association study of shared components of reading disability and language impairment

Written and verbal languages are neurobehavioral traits vital to the development of communication skills. Unfortunately, disorders involving these traits-specifically reading disability (RD) and language impairment (LI)-are common and prevent affected individuals from developing adequate communication skills, leaving them at risk for adverse academic, socioeconomic and psychiatric outcomes. Both RD and LI are complex traits that frequently co-occur, leading us to hypothesize that these disorders share genetic etiologies. To test this, we performed a genome-wide association study on individuals affected with both RD and LI in the Avon Longitudinal Study of Parents and Children. The strongest associations were seen with markers in ZNF385D (OR = 1.81, P = 5.45 × 10(-7) ) and COL4A2 (OR = 1.71, P = 7.59 × 10(-7) ). Markers within NDST4 showed the strongest associations with LI individually (OR = 1.827, P = 1.40 × 10(-7) ). We replicated association of ZNF385D using receptive vocabulary measures in the Pediatric Imaging Neurocognitive Genetics study (P = 0.00245). We then used diffusion tensor imaging fiber tract volume data on 16 fiber tracts to examine the implications of replicated markers. ZNF385D was a predictor of overall fiber tract volumes in both hemispheres, as well as global brain volume. Here, we present evidence for ZNF385D as a candidate gene for RD and LI. The implication of transcription factor ZNF385D in RD and LI underscores the importance of transcriptional regulation in the development of higher order neurocognitive traits. Further study is necessary to discern target genes of ZNF385D and how it functions within neural development of fluent language.

Polymorphisms in the gene encoding thymus-specific serine protease in the extended HLA complex: a potential candidate gene for autoimmune and HLA-associated diseases

Positive selection plays a role, together with negative selection, in the prevention of autoimmunity. Thymus-specific serine protease is highly expressed in the thymus and is believed to be involved in positive selection of T cells. The gene encoding thymus-specific serine protease (PRSS16) maps to the extended HLA complex, which harbours several genes predisposing for autoimmune diseases. Here we report the results of scanning the genetic region containing PRSS16 for polymorphisms. Twenty-two polymorphisms were identified, including one missense polymorphism, one deletion leading to elimination of five amino acids, as well as several SNPs in the promoter region.

Evidence for linkage and association with reading disability on 6p21.3-22

Reading disability (RD), or dyslexia, is a common heterogeneous syndrome with a large genetic component. Several studies have consistently found evidence for a quantitative-trait locus (QTL) within the 17 Mb (14.9 cM) that span D6S109 and D6S291 on chromosome 6p21.3-22. To characterize further linkage to the QTL, to define more accurately the location and the effect size, and to identify a peak of association, we performed Haseman-Elston and DeFries-Fulker linkage analyses, as well as transmission/disequilibrium, total-association, and variance-components analyses, on 11 quantitative reading and language phenotypes. One hundred four families with RD were genotyped with a new panel of 29 markers that spans 9 Mb of this region. Linkage results varied widely in degree of statistical significance for the different linkage tests, but multipoint analysis suggested a peak near D6S461. The average 6p QTL heritability for the 11 reading and language phenotypes was 0.27, with a maximum of 0.66 for orthographic choice. Consistent with the region of linkage described by these studies and others, there was a peak of transmission disequilibrium with a QTL centered at JA04 (chi2=9.48; empirical P=.0033; orthographic choice), and there was strong evidence for total association at this same marker (chi2=11.49; P=.0007; orthographic choice). Although the boundaries of the peak could not be precisely defined, the most likely location of the QTL is within a 4-Mb region surrounding JA04.

Peaks of linkage are localized by a BAC/PAC contig of the 6p reading disability locus

A gene for reading disability has been localized by nonparametric linkage to 6p21.3-p22 in several published reports. However, the lack of an uninterrupted genomic clone contig has made it difficult to determine accurate intermarker distances, precise marker order, and genetic boundaries and hinders direct comparisons of linkage. The search and discovery of the hemochromatosis gene (HFE) led to the creation of a bacterial artificial chromosome (BAC) and P-1 derived artificial chromosome (PAC) contig that extended physical maps 4 Mb from the MHC toward pter and localized new markers in that region [10-12]. Using this contig, we localized 124 sequence tagged sites, expressed sequence tags, and short tandem repeats including most of the markers in linkage with reading disability phenotypes, succinic semialdehyde dehydrogenase, GPLD1, prolactin, and 18 uncharacterized genes. This new contig joins and extends previously published physical maps to span the entire chromosome 6 reading disability genetic locus. Physical mapping data from the complete contig show overlap of the published linkage peaks for reading disability, provide accurate intermarker distances and order, and offer resources for generating additional markers and candidate genes for high resolution genetic studies in this region.

Human MHC class III and IV genes and disease associations

The major histocompatibility complex (MHC) was initially defined as the genetic locus encoding the Class I and Class II highly polymorphic cell surface antigens that are now known to present antigen to matched sets of T cell receptors. Genes for several diverse complement components, specifically Bf, C2, and C4 were found between the Class I and II genes, in a region later dubbed Class III. More recently, several genes have been described that are encoded in the telomeric end of the Class III region and that appear to be involved in both global and specific inflammatory responses. Due to this commonality of function this gene-rich region was dubbed Class IV, and includes the TNF family, AIF1, and HSP70. The genes of the Class III and Class IV regions are sufficiently divergent in sequence and structure so that clustering is not explicable in terms of gene duplication or divergence. We present some of the newer pertinent information and puzzling features of the genes embraced in the Class IV region and discuss possible roles in specific autoimmune diseases linked to this region.

Human GABA(B) receptor 1 gene: eight novel sequence variants

GABA (gamma-aminobutyric acid) is the principal inhibitory neurotransmitter in the brain. The human GABA(B) receptor (GABBR1) maps to the human leukocyte antigen (HLA) region of chromosome 6. Its function and location in a susceptibility region for schizophrenia, epilepsy, and dyslexia make GABBR1 a candidate gene for neurobehavioral disorders. We report the characterization of GABBR1 gene mutations in 100 chromosomes from a mixed American population. Eleven distinct mutations were found, including two previously reported missense mutations (A20V and G489S) and a previously reported silent 1977 T>C transition. Here, we report four novel silent substitutions (39C>T, 1473T>C, 1476T>C, 1545T>C) and four novel intron variants. These DNA variants may be useful in association and linkage studies of neurobehavioral disorders, and in pharmacogenetic studies of drugs targeting GABBR1.

Cloning of a novel MHC-encoded serine peptidase highly expressed by cortical epithelial cells of the thymus

Antigen presentation by cortical thymic epithelial cells (cTEC) during the positive selection of T cells has been shown to differ from that of other antigen-presenting cells. In the case of MHC class II presentation, cathepsin L as opposed to cathepsin S is responsible at least in part for the degradation of invariant chain. Other proteases, however, must be involved. We have identified a putative serine protease that is specifically expressed in the thymus. Encoded within the class I major histocompatibility complex (MHC) region, this gene has sequence homology with lysosomal prolylcarboxypeptidase, suggesting that it is a serine protease. We have, therefore, designated this gene thymus-specific serine protease (TSSP). In situ hybridization and immunofluorescence staining reveal that TSSP is expressed exclusively by cortical thymic epithelial cells, with the strongest staining noted around vessels and the thymic capsule. The identification of TSSP further supports the theory that MHC class II antigen processing and presentation in the thymic cortex involves a proteolytic milieu that differs from that of other antigen-presenting cells.

Human gamma-aminobutyric acid B receptor gene: complementary DNA cloning, expression, chromosomal location, and genomic organization

BACKGROUND

The 6p21.3 region of human chromosome 6 is a genetic locus for schizophrenia, juvenile myoclonic epilepsy, and dyslexia.

METHODS

Due to our interest in these disorders we performed complementary DNA (cDNA) hybridization selection on genomic DNA clones spanning this region to identify potential positional-candidate genes.

RESULTS

We identified a full-length cDNA with an open reading frame of 2883 bp corresponding to a predicted protein of 961 amino acids that shares greater than 95% homology with the rat gamma-aminobutyric acid B (GABAB) receptor. Northern blot hybridization identified a 4.4-kb transcript in human brain. The human gene mapped to two sites on 6p21.3 separated by 2 Mb. Sequence analysis of both sites showed that the centromeric gene is transcribed, whereas the telomeric site is likely a pseudogene. The transcribed gene is distributed over 22 exons spanning 18 kb of genomic DNA.

CONCLUSIONS

The genomic location, tissue expression, and function of the human GABAB receptor gene suggest that it is an important positional-candidate for the neurobehavioral disorders with a genetic locus on 6p21.3. In addition, delineation of the genomic organization will now permit it to be integrated as part of pharmacogenetic studies in trials of anxiolytic, narcotic, antiepileptic, and fluoxetine therapies.

Haplotype analysis of hemochromatosis: evaluation of different linkage-disequilibrium approaches and evolution of disease chromosomes

We applied several types of linkage-disequilibrium calculations to analyze the hereditary hemochromatosis (hh) locus. Twenty-four polymorphic markers in the major histocompatibility complex (MHC) class I region were used to haplotype hh and normal chromosomes. A total of 169 hh and 161 normal chromosomes were analyzed. Disequilibrium values were found to be high over an unusually large region beginning 150 kb centromeric of HLA-A and extending nearly 5 Mb telomeric of it. Recombination in this region was approximately 28% of the expected value. This low level of recombination contributes to the unusually broad region of linkage disequilibrium found with hh. The strongest disequilibrium was found at locus HLA-H (delta = .84) and at locus D6S2239 (delta = .85), a marker approximately 10 kb telomeric to HLA-H. All disequilibrium methods employed in this study found peak disequilibrium at HLA-H or D6S2239. The cys282tyr mutation in HLA-H, a candidate gene for hh, was found in 85% of disease chromosomes. A haplotype phylogeny for hh chromosomes was constructed and suggests that the mutation associated with the most common haplotype occurred relatively recently. The age of the hh mutation was estimated to be approximately 60-70 generations. Disequilibrium was maintained over a greater distance for hh-carrying chromosomes, consistent with a recent mutation for hh. Our data provide a reasonable explanation for previous difficulties in localizing the hh locus and provide an evolutionary history for disease chromosomes.

Recombinations defining centromeric and telomeric borders for the hereditary haemochromatosis locus

Hereditary haemochromatosis (HFE) is a common inherited disorder, affecting approximately five per thousand white people of northern European descent. Genetic linkage and linkage disequilibrium studies indicate that the disease locus is tightly linked to HLA-A and D6S105. Recombination between HFE and HLA class I loci is known to be rare. We report here two pedigrees in which recombinations telomeric of HLA-A occurred. These recombinant events define new centromeric and telomeric borders for the HFE locus.

FISH-Mapped CEPH YACs spanning 0 to 46 cM on human chromosome 6

Seventy-six CEPH YACs were mapped by fluorescence in situ hybridization (FISH) to human metaphase chromosomes. These clones have been ordered from pter to 46 cM by combining the results of FISH with sequence-tagged site content mapping using data from the public databases. This created a minimal tiling path containing at least 37 Mb of human genomic DNA from 0 to 46 cM on chromosome 6 that contains up to four gaps not greater than 200 kb. These data provide an integration of the FLpter physical map values with cytogenetic band localization and markers on the genetic and radiation hybrid maps. We also assessed YAC chimerism and placed three additional Whitehead contigs (WC952, WC799, WC436) within the integrated map.

A transcription map of the major histocompatibility complex (MHC) class I region

We have applied cDNA hybridization selection to nine YACs spanning 3 Mb of genomic DNA from a region centromeric to HLA-A to the histone cluster that lies telomeric to the human major histocompatibility complex (MHC). In addition to Class I genes and pseudogenes, we describe over 63 genes and 23 additional expressed sequence tags distributed throughout the region. Many of the full-length genes belong to gene families. Prominent among these are a group of genes encoding proteins showing homology to the carboxyl-terminal sequences of butyrophilin and an additional group of zinc finger genes. We also detected several previously undefined genes that are specifically expressed in cells of the immune system, indicating a more complex role of the MHC in the immune response than has been appreciated.

Cloning of a new "finger" protein gene (ZNF173) within the class I region of the human MHC

The human major histocompatability complex contains genes of both immune and nonimmune importance. Recently, several genes encoding novel, non-HLA products have been described in this area. We have performed positional cloning of short fragment cDNA sequences from the class I region of the human MHC using a hybridization selection approach. This report describes isolation of full-length cDNA clones and partial genomic clones that encode a protein that contains two domains rich in cysteine and histidine similar to those characteristic of metal-dependent DNA binding proteins (C3HC4). The predicted protein also contains a domain thought to form a coiled-coil that may promote dimerization. A third feature is a polyglutamic acid region near the carboxyl terminus of the conceptual protein. Because of these properties, we have named this gene product acid finger protein (AFP). Although the biological role of AFP is unknown at present, one potential function is binding of nucleic acids. The gene (ZNF173) is expressed in multiple tissues and is conserved among mammals. In particular, the mouse and human coding regions are highly conserved. In addition to AFP, other related sequences have been localized to the MHC, suggesting that multiple AFP-like genes exist in this area.

Isolation of novel non-HLA gene fragments from the hemochromatosis region (6p21.3) by cDNA hybridization selection

It has previously been shown that cDNA hybridization selection can identify and recover novel genes from large cloned genomic DNA such as cosmids or YACs. In an effort to identify candidate genes for hemochromatosis, this technique was applied to a 320-kb YAC containing the HLA-A gene. A short fragment cDNA library derived from human duodenum was selected with the YAC DNA. Ten novel gene fragments were isolated, characterized, and localized on the physical map of the YAC.

Physical and genetic mapping of the telomeric major histocompatibility complex region in man and relevance to the primary hemochromatosis gene (HFE)

We performed pulsed-field gel electrophoresis (PFGE) on genomic DNA from a radiation hybrid (RH) cell line and constructed a high-resolution physical map of the major histocompatibility complex class I region in 6p21.3, where the gene for primary hemochromatosis (HFE) is believed to be located. Due to the intact microsegment of hemizygous human genomic DNA preserved in the RH cell line, simplified and distinct restriction fragment banding patterns were generated. Using the RH cell line, we were able to extend the physical map of the HLA class I region to about 3000 kb, order the known HLA class I genes from centromere to telomere: HLA-B, -C, -E, (-A, -H, -G), and -F, and orient the HLA-F gene along the chromosome. The proximity of HLA-F to HLA-A was confirmed by linkage and linkage disequilibrium analysis. This study shows that RH cell lines can be useful for constructing long-range physical maps in specific regions of the human genome with PFGE. Physical and genetic mapping studies of this region are consistent with a localization of the HFE gene proximal or distal to HLA-A.

Differential expression of novel Gs alpha signal transduction protein cDNA species

The Gs alpha guanine nucleotide-binding signal transduction protein is part of a heterotrimeric complex that is involved in the stimulation of adenylate cyclase upon activation of membrane receptors. We report the characterization of 16 Gs alpha cDNA clones isolated from the human adult retina and fetal eye libraries. Molecular heterogeneity in the 5'-region defines four novel Gs alpha cDNA species which are generated either by alternate splicing or by using alternative promoter. The novel exons upstream of exon 2 interrupt the highly conserved 'region A' in the Gs alpha polypeptide. Non-AUG codons in the novel 5'-exon can initiate translation of these Gs alpha species in vitro. Reverse transcription of total RNA coupled with polymerase chain reaction (RTPCR) using specific primers and in situ hybridization to mRNA in baboon tissue sections with a specific oligonucleotide probe show a high level of expression of these species in retina and brain but not in liver. Differential expression of alternatively spliced Gs alpha species suggests novel signal transducing pathways.