Low codon bias and high rates of synonymous substitution in Drosophila hydei and D. melanogaster histone genes

We have evaluated codon usage bias in Drosophila histone genes and have obtained the nucleotide sequence of a 5,161-bp D. hydei histone gene repeat unit. This repeat contains genes for all five histone proteins (H1, H2a, H2b, H3, and H4) and differs from the previously reported one by a second EcoRI site. These D. hydei repeats have been aligned to each other and to the 5.0-kb (i.e., long) and 4.8-kb (i.e., short) histone repeat types from D. melanogaster. In each species, base composition at synonymous sites is similar to the average genomic composition and approaches that in the small intergenic spacers of the histone gene repeats. Accumulation of synonymous changes at synonymous sites after the species diverged is quite high. Both of these features are consistent with the relatively low codon usage bias observed in these genes when compared with other Drosophila genes. Thus, the generalization that abundantly expressed genes in Drosophila have high codon bias and low rates of silent substitution does not hold for the histone genes.      

Cloning, genomic structure, and expression profiles of TULIP1 (GARNL1), a brain-expressed candidate gene for 14q13-linked neurological phenotypes, and its murine homologue

Previously, we have described the clinical and molecular characterization of a de novo 14q13.1-q21.1 microdeletion, less than 3.5 Mb in size, in a patient with severe microcephaly, psychomotor retardation, and other clinical anomalies. Here we report the characterization of the genomic structure of the human tuberin-like protein gene 1 (TULIP1; approved gene symbol GARNL1), a CpGisland-associated, brain-expressed candidate gene for the neurological findings in our patient, and its murine homologue. The human TULIP1 gene was mapped to chromosome band 14q13.2 by fluorescence in situ hybridization of BAC clone RP11-355C3 (GenBank Accession No. AL160231), containing the 3' region of the gene. TULIP1 spans about 271 kb of human genomic DNA and is divided into 41 exons. An untranscribed, processed pseudogene of TULIP1 was found on human chromosome band 9q31.1. The active locus TULIP1, encoding a predicted protein of 2036 amino acids, is expressed ubiquitously in pre- and postnatal human tissues. The murine homologue Tulip1 spans about 220 kb of mouse genomic DNA and is also divided into 41 exons, encoding a predicted protein of 2035 amino acids. No pseudogene could be found in the available mouse sequence data. Several splicing variants were found. Considering the location, expression profile, and predicted function, TULIP1 is a strong candidate for several neurological features seen in 14q deletion patients. Additionally we searched for mutations in the coding region of TULIP1 in subjects from a family with idiopathic basal ganglia calcification (IBGC; Fahr disease), previously linked to chromosome 14q. We identified two novel SNPs in the intron-exon boundaries; however, they did not segregate only with affected subjects in the predicted model of an autosomal dominant disease such as IBGC.     

Neural progenitor genes. Germinal zone expression and analysis of genetic overlap in stem cell populations

The identification of the genes regulating neural progenitor cell (NPC) functions is of great importance to developmental neuroscience and neural repair. Previously, we combined genetic subtraction and microarray analysis to identify genes enriched in neural progenitor cultures. Here, we apply a strategy to further stratify the neural progenitor genes. In situ hybridization demonstrates expression in the central nervous system germinal zones of 54 clones so identified, making them highly relevant for study in brain and neural progenitor development. Using microarray analysis we find 73 genes enriched in three neural stem cell (NSC)-containing populations generated under different conditions. We use the custom microarray to identify 38 "stemness" genes, with enriched expression in the three NSC conditions and present in both embryonic stem cells and hematopoietic stem cells. However, comparison of expression profiles from these stem cell populations indicates that while there is shared gene expression, the amount of genetic overlap is no more than what would be expected by chance, indicating that different stem cells have largely different gene expression patterns. Taken together, these studies identify many genes not previously associated with neural progenitor cell biology and also provide a rational scheme for stratification of microarray data for functional analysis.     

A genetic analysis of neural progenitor differentiation

Genetic mechanisms regulating CNS progenitor function and differentiation are not well understood. We have used microarrays derived from a representational difference analysis (RDA) subtraction in a heterogeneous stem cell culture system to systematically study the gene expression patterns of CNS progenitors. This analysis identified both known and novel genes enriched in progenitor cultures. In situ hybridization in a subset of clones demonstrated that many of these genes were expressed preferentially in germinal zones, some showing distinct ventricular or subventricular zone labeling. Several genes were also enriched in hematopoietic stem cells, suggesting an overlap of gene expression in neural and hematopoietic progenitors. This combination of methods demonstrates the power of using custom microarrays derived from RDA-subtracted libraries for both gene discovery and gene expression analysis in the central nervous system.     

Thresholding rules for recovering a sparse signal from microarray experiments

We consider array experiments that compare expression levels of a high number of genes in two cell lines with few repetitions and with no subject effect. We develop a statistical model that illustrates under which assumptions thresholding is optimal in the analysis of such microarray data. The results of our model explain the success of the empirical rule of two-fold change. We illustrate a thresholding procedure that is adaptive to the noise level of the experiment, the amount of genes analyzed, and the amount of genes that truly change expression level. This procedure, in a world of perfect knowledge on noise distribution, would allow reconstruction of a sparse signal, minimizing the false discovery rate. Given the amount of information actually available, the thresholding rule described provides a reasonable estimator for the change in expression of any gene in two compared cell lines.     

Evidence for a Language Quantitative Trait Locus on Chromosome 7q in Multiplex Autism Families

Autism is a syndrome characterized by deficits in language and social skills and by repetitive behaviors. We hypothesized that potential quantitative trait loci (QTLs) related to component autism endophenotypes might underlie putative or significant regions of autism linkage. We performed nonparametric multipoint linkage analyses, in 152 families from the Autism Genetic Resource Exchange, focusing on three traits derived from the Autism Diagnostic Interview: "age at first word," "age at first phrase," and a composite measure of "repetitive and stereotyped behavior." Families were genotyped for 335 markers, and multipoint sib pair linkage analyses were conducted. Using nonparametric multipoint linkage analysis, we found the strongest QTL evidence for age at first word on chromosome 7q (nonparametric test statistic [Z] 2.98; P=.001), and subsequent linkage analyses of additional markers and association analyses in the same region supported the initial result (Z=2.85, P=.002; chi(2)=18.84, df 8, P=.016). Moreover, the peak fine-mapping result for repetitive behavior (Z=2.48; P=.007) localized to a region overlapping this language QTL. The putative autism-susceptibility locus on chromosome 7 may be the result of separate QTLs for the language and repetitive or stereotyped behavior deficits that are associated with the disorder.