Genes related to sex steroids, neural growth, and social-emotional behavior are associated with autistic traits, empathy, and Asperger syndrome.

Genetic studies of autism spectrum conditions (ASC) have mostly focused on the "low functioning" severe clinical subgroup, treating it as a rare disorder. However, ASC is now thought to be relatively common ( approximately 1%), and representing one end of a quasi-normal distribution of autistic traits in the general population. Here we report a study of common genetic variation in candidate genes associated with autistic traits and Asperger syndrome (AS). We tested single nucleotide polymorphisms in 68 candidate genes in three functional groups (sex steroid synthesis/transport, neural connectivity, and social-emotional responsivity) in two experiments. These were (a) an association study of relevant behavioral traits (the Empathy Quotient (EQ), the Autism Spectrum Quotient (AQ)) in a population sample (n=349); and (b) a case-control association study on a sample of people with AS, a "high-functioning" subgroup of ASC (n=174). 27 genes showed a nominally significant association with autistic traits and/or ASC diagnosis. Of these, 19 genes showed nominally significant association with AQ/EQ. In the sex steroid group, this included ESR2 and CYP11B1. In the neural connectivity group, this included HOXA1, NTRK1, and NLGN4X. In the socio-responsivity behavior group, this included MAOB, AVPR1B, and WFS1. Fourteen genes showed nominally significant association with AS. In the sex steroid group, this included CYP17A1 and CYP19A1. In the socio-emotional behavior group, this included OXT. Six genes were nominally associated in both experiments, providing a partial replication. Eleven genes survived family wise error rate (FWER) correction using permutations across both experiments, which is greater than would be expected by chance. CYP11B1 and NTRK1 emerged as significantly associated genes in both experiments, after FWER correction (P<0.05). This is the first candidate-gene association study of AS and of autistic traits. The most promising candidate genes require independent replication and fine mapping.

Loss-of-function mutations in growth differentiation factor-1 (GDF1) are associated with congenital heart defects in humans.

Congenital heart defects (CHDs) are among the most common birth defects in humans (incidence 8-10 per 1,000 live births). Although their etiology is often poorly understood, most are considered to arise from multifactorial influences, including environmental and genetic components, as well as from less common syndromic forms. We hypothesized that disturbances in left-right patterning could contribute to the pathogenesis of selected cardiac defects by interfering with the extrinsic cues leading to the proper looping and vessel remodeling of the normally asymmetrically developed heart and vessels. Here, we show that heterozygous loss-of-function mutations in the human GDF1 gene contribute to cardiac defects ranging from tetralogy of Fallot to transposition of the great arteries and that decreased TGF- beta signaling provides a framework for understanding their pathogenesis. These findings implicate perturbations of the TGF- beta signaling pathway in the causation of a major subclass of human CHDs.

Essential structural and functional determinants within the forkhead domain of FOXC1.

The forkhead domain (FHD)-containing developmental transcription factor FOXC1 is mutated in patients presenting with Axenfeld-Rieger malformations. In this paper, we report the introduction of positive, negative or neutral charged amino acids into critical positions within the forkhead domain of FOXC1 in an effort to better understand the essential structural and functional determinants within the FHD. We found that FOXC1 is intolerant of mutations at I87. Additionally, alterations of amino acids within alpha-helix 1 of the FOXC1 FHD affected both nuclear localization and transactivation. Amino acids within alpha-helix 3 were also found to be necessary for transactivation and can have roles in correct localization. Interestingly, changing amino acids within alpha-helix 3, particularly R127, resulted in altered DNA-binding specificity and granted FOXC1 the ability to bind to a novel DNA sequence. Given the limited topological variation of FHDs, due to the high conservation of residues, we anticipate that models of forkhead domain function derived from these data will be relevant to other members of the FOX family of transcription factors.

Molecular evolution of the homeodomain family of transcription factors.

The homeodomain family of transcription factors plays a fundamental role in a diverse set of functions that include body plan specification, pattern formation and cell fate determination during metazoan development. Members of this family are characterized by a helix-turn-helix DNA-binding motif known as the homeodomain. Homeodomain proteins regulate various cellular processes by specifically binding to the transcriptional control region of a target gene. These proteins have been conserved across a diverse range of species, from yeast to human. A number of inherited human disorders are caused by mutations in homeodomain-containing proteins. In this study, we present an evolutionary classification of 129 human homeodomain proteins. Phylogenetic analysis of these proteins, whose sequences were aligned based on the three-dimensional structure of the homeodomain, was performed using a distance matrix approach. The homeodomain proteins segregate into six distinct classes, and this classification is consistent with the known functional and structural characteristics of these proteins. An ancestral sequence signature that accurately describes the unique sequence characteristics of each of these classes has been derived. The phylogenetic analysis, coupled with the chromosomal localization of these genes, provides powerful clues as to how each of these classes arose from the ancestral homeodomain.

Analyses of the effects that disease-causing missense mutations have on the structure and function of the winged-helix protein FOXC1.

Five missense mutations of the winged-helix FOXC1 transcription factor, found in patients with Axenfeld-Rieger (AR) malformations, were investigated for their effects on FOXC1 structure and function. Molecular modeling of the FOXC1 forkhead domain predicted that the missense mutations did not alter FOXC1 structure. Biochemical analyses indicated that, whereas all mutant proteins correctly localize to the cell nucleus, the I87M mutation reduced FOXC1-protein levels. DNA-binding experiments revealed that, although the S82T and S131L mutations decreased DNA binding, the F112S and I126M mutations did not. However, the F112S and I126M mutations decrease the transactivation ability of FOXC1. All the FOXC1 mutations had the net effect of reducing FOXC1 transactivation ability. These results indicate that the FOXC1 forkhead domain contains separable DNA-binding and transactivation functions. In addition, these findings demonstrate that reduced stability, DNA binding, or transactivation, all causing a decrease in the ability of FOXC1 to transactivate genes, can underlie AR malformations.

The Homeodomain Resource: sequences, structures, DNA binding sites and genomic information.

The Homeodomain Resource is an annotated collection of non-redundant protein sequences, three-dimensional structures and genomic information for the homeodomain protein family. Release 3.0 contains 795 full-length homeodomain-containing sequences, 32 experimentally-derived structures and 143 homeo-box loci implicated in human genetic disorders. Entries are fully hyperlinked to facilitate easy retrieval of the original records from source databases. A simple search engine with a graphical user interface is provided to query the component databases and assemble customized data sets. A new feature for this release is the addition of DNA recognition sites for all human homeodomain proteins described in the literature. The Homeodomain Resource is freely available through the World Wide Web at http://genome.nhgri.nih.gov/homeodomain.

Mutation of a gene encoding a putative chaperonin causes McKusick-Kaufman syndrome.

McKusick-Kaufman syndrome (MKKS, MIM 236700) is a human developmental anomaly syndrome comprising hydrometrocolpos (HMC), postaxial polydactyly (PAP) and congenital heart disease (CHD). MKKS has been mapped in the Old Order Amish population to 20p12, between D20S162 and D20S894 (ref. 3). Here we describe the identification of a gene mutated in MKKS. We analysed the approximately 450-kb candidate region by sample sequencing, which revealed the presence of several known genes and EST clusters. We evaluated candidate transcripts by northern-blot analysis of adult and fetal tissues. We selected one transcript with widespread expression, MKKS, for analysis in a patient from the Amish pedigree and a sporadic, non-Amish case. The Old Order Amish patient was found to be homozygous for an allele that had two missense substitutions and the non-Amish patient was a compound heterozygote for a frameshift mutation predicting premature protein truncation and a distinct missense mutation. The MKKS predicted protein shows amino acid similarity to the chaperonin family of proteins, suggesting a role for protein processing in limb, cardiac and reproductive system development. We believe that this is the first description of a human disorder caused by mutations affecting a putative chaperonin molecule.

MLH3: a DNA mismatch repair gene associated with mammalian microsatellite instability.

DNA mismatch repair is important because of its role in maintaining genomic integrity and its association with hereditary non-polyposis colon cancer (HNPCC). To identify new human mismatch repair proteins, we probed nuclear extracts with the conserved carboxy-terminal MLH1 interaction domain. Here we describe the cloning and complete genomic sequence of MLH3, which encodes a new DNA mismatch repair protein that interacts with MLH1. MLH3 is more similar to mismatch repair proteins from yeast, plants, worms and bacteria than to any known mammalian protein, suggesting that its conserved sequence may confer unique functions in mice and humans. Cells in culture stably expressing a dominant-negative MLH3 protein exhibit microsatellite instability. Mlh3 is highly expressed in gastrointestinal epithelium and physically maps to the mouse complex trait locus colon cancer susceptibility I (Ccs1). Although we were unable to identify a mutation in the protein-coding region of Mlh3 in the susceptible mouse strain, colon tumours from congenic Ccs1 mice exhibit microsatellite instability. Functional redundancy among Mlh3, Pms1 and Pms2 may explain why neither Pms1 nor Pms2 mutant mice develop colon cancer, and why PMS1 and PMS2 mutations are only rarely found in HNPCC families.

The homeodomain resource: a prototype database for a large protein family.

The Homeodomain Resource is an annotated collection of non-redundant protein sequences, three-dimensional structures and genomic information for the homeodomain protein family. Release 2.0 contains 765 full-length homeodomain-containing sequences, 29 experimentally derived structures and 116 homeobox loci implicated in human genetic disorders. Entries are fully hyperlinked to facilitate easy retrieval of the original records from source databases. A simple search engine with a graphical user interface is provided to query the component databases and assemble customized data sets. A new feature for this release is the addition of more automated methods for database searching, maintenance and implementation of efficient data management. The Homeodomain Resource is freely available through the WWW at http://genome.nhgri.nih.gov/homeodomain

Threading analysis of the Pitx2 homeodomain: predicted structural effects of mutations causing Rieger syndrome and iridogoniodysgenesis.

Mutations in the homeobox gene PITX2 are responsible for a range of clinical phenotypes involving ocular and craniofacial development. Several mutations within the Pitx2 homeodomain region are specifically responsible for the development of the related autosomal-dominant disorders Rieger syndrome and iridogoniodysgenesis. To address the question of the structural effect of disease-causing mutations on the Pitx2 homeodomain, we used threading techniques to examine the tertiary structure of the Pitx2 wild-type and mutant homeodomain, using the crystal structure of Drosophila engrailed homeodomain bound with DNA as a template [Kissinger et al., 1990]. The threading analysis reveals that the wild-type Pitx2 homeodomain is indeed capable of forming the typical three-helical bundle-fold characteristic of homeodomain proteins. Energy calculations indicate that the homeodomain structure is stabilized primarily by hydrophobic interactions between residues at the helical interface. Point mutations responsible for the development of these genetic disorders were also examined; the results suggest that these mutations lead to the inability of Pitx2 to adopt its proper structure and bind to the regulatory sequences of its target gene(s), which in turn affects its metabolic role in the cell. Published 1999 Wiley-Liss, Inc.

Threading analysis of prospero-type homeodomains.

The homeodomain is a common structural motif found in many transcription factors involved in cell fate determination during development. We have used threading analysis techniques to predict whether the atypical homeodomain of prospero (pros) family members could form the three-helical homeodomain structural motif, even though these proteins are not statistically similar to canonical homeodomains as assessed by BLAST searches. Amino acid sequences of these divergent homeodomain proteins were threaded through the X-ray coordinates of the Drosophila engrailed homeodomain protein [23]. The analysis confirms that the prospero class of homeodomain proteins is indeed capable of forming the homeodomain structure despite its low degree of sequence identity to the canonical homeodomain. Energy calculations indicate that the homeodomain structure is stabilized primarily by hydrophobic interactions between residues at the helical interfaces. Although the atypical prospero-type homeodomain shows very little sequence similarity when compared to other homeodomain proteins, the critical amino acids responsible for maintaining the three-dimensional structure are highly conserved. A number of other homeodomain proteins, such as PHO2p from Saccharomyces and Pax6 from human, were also included in the threading analysis and were shown to be able to form the engrailed structure, indicating that there are no rigid overall sequence requirements for the formation of the homeodomain structural motif. Based on the threading experiments and the subsequent structural alignment, a new amino acid signature that unambiguously identifies the prospero-type proteins was deduced.

The Homeodomain Resource: sequences, structures and genomic information.

The Homeodomain Resource is a comprehensive collection of sequence, structure and genomic information on the homeodomain protein family. Available through the Resource are both full-length and domain-only sequence data, as well as X-ray and NMR structural data for proteins and protein-DNA complexes. Also available is information on human genetic diseases and disorders in which proteins from the homeodomain family play an important role; genomic information includes relevant gene symbols, cytogenetic map locations, and specific mutation data. Search engines are provided to allow users to easily query the component databases and assemble specialized data sets. The Homeodomain Resource is available through the World Wide Web at http://genome.nhgri.nih.gov/homeodomain

Chicken homeobox gene Prox 1 related to Drosophila prospero is expressed in the developing lens and retina.

Prox 1 is the vertebrate homolog of Drosophila prospero, a gene known to be expressed in the lens-secreting cone cells of fly ommatidia. Chicken Prox 1 cDNAs were isolated from 14 day embryonic chicken lenses, and a complete open reading frame encoding an 83 kDa protein was elucidated. The homeodomains of chicken and mouse Prox 1 are identical at the amino acid level and are 65-67% similar to the homeodomains of Drosophila and C. elegans prospero. The homology between these proteins extends beyond the homeodomain. There is 56% identity between chicken Prox 1 and Drosophila prospero in the C-terminal region downstream of the homeodomain, whereas there is little similarity upstream of the homeodomain. Prox 1 is expressed most actively in the developing lens and midgut and at lower levels in the developing brain, heart, muscle, and retina. cDNA sequencing has established that there are alternatively spliced forms of the single Prox 1 gene, which probably account for the two abundant RNAs of about 2 and 8 kb and two less abundant RNAs close to 3.5 kb in length in the lens. In the lens fibers, only the shortest mRNA was present, whereas, in the epithelial cells, both short and long mRNAs were detected. By using in situ hybridization, expression of the Prox 1 gene was first detected at stage 14 in the early lens placode and slightly preceded the expression of delta 1-crystallin, the first crystallin gene expressed in the developing chicken lens. At later stages of development, Prox 1 mRNA was observed throughout the lens, but it appeared more abundant around the bow region of the equator than in the anterior epithelium or the fibers. In the retina, expression of the Prox 1 gene was detected mainly in the inner nuclear layer during later stages of histogenesis. The conserved pattern of Prox 1/prospero gene expression in vertebrates and Drosophila suggests that Prox 1, like Pax-6, may be essential for eye development in different systematic groups.

Common core sequences are found in skeletal muscle slow- and fast-fiber-type-specific regulatory elements.

The molecular mechanisms generating muscle diversity during development are unknown. The phenotypic properties of slow- and fast-twitch myofibers are determined by the selective transcription of genes coding for contractile proteins and metabolic enzymes in these muscles, properties that fail to develop in cultured muscle. Using transgenic mice, we have identified regulatory elements in the evolutionarily related troponin slow (TnIs) and fast (TnIf) genes that confer specific transcription in either slow or fast muscles. Analysis of serial deletions of the rat TnIs upstream region revealed that sequences between kb -0.95 and -0.5 are necessary to confer slow-fiber-specific transcription; the -0.5-kb fragment containing the basal promoter was inactive in five transgenic mouse lines tested. We identified a 128-bp regulatory element residing at kb -0.8 that, when linked to the -0.5-kb TnIs promoter, specifically confers transcription to slow-twitch muscles. To identify sequences directing fast-fiber-specific transcription, we generated transgenic mice harboring a construct containing the TnIs kb -0.5 promoter fused to a 144-bp enhancer derived from the quail TnIf gene. Mice harboring the TnIf/TnIs chimera construct expressed the transgene in fast but not in slow muscles, indicating that these regulatory elements are sufficient to confer fiber-type-specific transcription. Alignment of rat TnIs and quail TnIf regulatory sequences indicates that there is a conserved spatial organization of core elements, namely, an E box, a CCAC box, a MEF-2-like sequence, and a previously uncharacterized motif. The core elements were shown to bind their cognate factors by electrophoretic mobility shift assays, and their mutation demonstrated that the TnIs CCAC and E boxes are necessary for transgene expression. Our results suggest that the interaction of closely related transcriptional protein-DNA complexes is utilized to specify fiber type diversity.

Sequence and expression of chicken beta A2- and beta B3-crystallins.

Crystallins are a diverse group of proteins that contribute to the transparency and refractive properties of the eye lens. Previously, the chicken orthologs of four out of the six known bovine beta-crystallin genes have been cloned and sequenced. In the present study, cDNAs corresponding to the chicken orthologs of beta A2- and beta B3-crystallin, the two previously unidentified chicken beta-crystallins, have been isolated. In addition, sequence analysis of three independent chicken beta B2-crystallin cDNAs yielded a deduced connecting peptide sequence which is considerably shorter than that reported previously. Thus, direct homologs of all of the known bovine beta-crystallins are expressed in the chicken lens. This demonstrates that the duplications giving rise to the known vertebrate beta-crystallins occurred over 300 million years ago. beta B2- and beta B3/A1-crystallin are the most highly conserved of the beta-crystallins suggesting that these genes may be important for other functions besides their refractive role in the lens. By Northern blot hybridization analysis, both beta A2- and beta B3-crystallin were shown to be lens-specific in the chicken embryo. The relative levels of beta A2-crystallin remained stable from five days of embryogenesis until adulthood, while the relative amounts of beta B3-crystallin increased until hatching and were appreciably lower in the adult lens. Approximately equal relative amounts of beta A2-crystallin mRNA were found in the lens epithelia and fibers of 5 day embryonic chicken embryos; by contrast, beta B3-crystallin mRNA was detected preferentially in the lens fibers. These data in combination with previous studies suggest that beta-crystallin genes are regulated independently from each other in the developing chicken lens. The elucidation of the primary structures for all seven chicken beta-crystallin polypeptides will facilitate future studies on the structure/function relationships responsible for lens transparency and on the molecular basis for beta-crystallin gene expression during development.

Isolation and structure of the rat gene encoding troponin I(slow).

Troponin I is a myofibrillar protein involved in the Ca(2+)-mediated regulation of actomyosin ATPase activity. We report here the isolation and characterization of the gene coding for the slow-muscle-specific isoform of the rat troponin I polypeptide (TpnI). Using restriction mapping, PCR mapping and partial DNA sequencing, we have determined the exon/intron arrangement of this gene. The transcription unit is 10.5-kb long and contains nine exons ranging in size from 4 bp to 330 bp. The rat TpnI(slow) gene is interrupted by large intervening sequences; a 3.3-kb intron separates the 5' untranslated exons from the protein-coding exons. Comparison of the structure of rat TpnI(slow) with that of quail TpnI(fast) reveals that they have a similar intron/exon organization. The 5' untranslated region of the rat gene contains an additional exon, otherwise, the positions of introns and coding exons map to essentially identical regions in both genes.

cis-acting sequences of the rat troponin I slow gene confer tissue- and development-specific transcription in cultured muscle cells as well as fiber type specificity in transgenic mice.

Transcription of the genes coding for troponin I slow (TnIslow) and other contractile proteins is activated during skeletal muscle differentiation, and their expression is later restricted to specific fiber types during maturation. We have isolated and characterized the rat TnIslow gene in order to begin elucidating its regulation during myogenesis. Transcriptional regulatory regions were delineated by using constructs, containing TnIslow gene sequences driving the expression of the chloramphenicol acetyltransferase (CAT) reporter gene, that were transiently transfected into undifferentiated and differentiated C2C12 cells. TnIslow 5'-flanking sequences directed transcription specifically in differentiated cells. However, transcription rates were approximately 10-fold higher in myotubes transfected with constructs containing the 5'-flanking sequences plus the intragenic region residing upstream of the translation initiation site (introns 1 and 2), indicative of interactions between elements residing upstream and in the introns of the gene. Deletion analysis of the 5' region of the TnIslow gene showed that the 200 bp upstream of the transcription initiation site is sufficient to confer differentiation-specific transcription in C2C12 myocytes. MyoD consensus binding sites were found both in the upstream 200-bp region and in a region residing in the second intron that is highly homologous to the quail TnIfast enhancer. Transactivation experiments using transfected NIH 3T3 fibroblasts with TnI-CAT constructs containing intragenic and/or upstream sequences and with the myogenic factors MyoD, myogenin, and MRF4 showed different potentials of these factors to induce transcription. Transgenic mice harboring the rat TnI-CAT fusion gene expressed the reporter specifically in the skeletal muscle. Furthermore, CAT levels were approximately 50-fold higher in the soleus than in the extensor digitorum longus, gastrocnemius, or tibialis muscle, indicating that the regulatory elements that restrict TnI transcription to slow-twitch myofibers reside in the sequences we have analyzed.