Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA Microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design microarray probes with high specificities and sensitivities. Highly specific probes correspond to probes having unique DNA sequences; whereas highly sensitive probes correspond to those with melting temperature within a desired range and having no secondary structure. The selection of these probes from a set of functional DNA sequences (exons) constitutes a computationally expensive discrete non-linear search problem. We delegate the search task to a simple yet effective Evolution Strategy algorithm. The computational efficiency is also greatly improved by making use of an available bioinformatics tool.

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design probes with high specificities, i.e. uniqueness, and also sensitivities, i.e., suitable melting temperature and no secondary structure. We make use of available biology tools to gain necessary sequence information of human chromosome 12, and combined with evolutionary strategy (ES) to find unique subsequences representing all predicted exons. The results are presented and discussed.

Ultraconserved elements in the human genome

There are 481 segments longer than 200 base pairs (bp) that are absolutely conserved (100% identity with no insertions or deletions) between orthologous regions of the human, rat, and mouse genomes. Nearly all of these segments are also conserved in the chicken and dog genomes, with an average of 95 and 99% identity, respectively. Many are also significantly conserved in fish. These ultraconserved elements of the human genome are most often located either overlapping exons in genes involved in RNA processing or in introns or nearby genes involved in the regulation of transcription and development. Along with more than 5000 sequences of over 100 bp that are absolutely conserved among the three sequenced mammals, these represent a class of genetic elements whose functions and evolutionary origins are yet to be determined, but which are more highly conserved between these species than are proteins and appear to be essential for the ontogeny of mammals and other vertebrates.

Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function

The mammalian transcriptome contains many nonprotein-coding RNAs (ncRNAs), but most of these are of unclear significance and lack strong sequence conservation, prompting suggestions that they might be non-functional. However, certain long functional ncRNAs such as Air and Xist are also poorly conserved. In this article, we systematically analyzed the conservation of several groups of functional ncRNAs, including miRNAs, snoRNAs and longer ncRNAs whose function has been either documented or confidently predicted. As expected, miRNAs and snoRNAs were highly conserved. By contrast, the longer functional non-micro, non-sno ncRNAs were much less conserved with many displaying rapid sequence evolution. Our findings suggest that longer ncRNAs are under the influence of different evolutionary constraints and that the lack of conservation displayed by the thousands of candidate ncRNAs does not necessarily signify an absence of function.

Segmenting eukaryotic genomes with the generalized Gibbs sampler

Eukaryotic genomes display segmental patterns of variation in various properties, including GC content and degree of evolutionary conservation. DNA segmentation algorithms are aimed at identifying statistically significant boundaries between such segments. Such algorithms may provide a means of discovering new classes of functional elements in eukaryotic genomes. This paper presents a model and an algorithm for Bayesian DNA segmentation and considers the feasibility of using it to segment whole eukaryotic genomes. The algorithm is tested on a range of simulated and real DNA sequences, and the following conclusions are drawn. Firstly, the algorithm correctly identifies non-segmented sequence, and can thus be used to reject the null hypothesis of uniformity in the property of interest. Secondly, estimates of the number and locations of change-points produced by the algorithm are robust to variations in algorithm parameters and initial starting conditions and correspond to real features in the data. Thirdly, the algorithm is successfully used to segment human chromosome 1 according to GC content, thus demonstrating the feasibility of Bayesian segmentation of eukaryotic genomes. The software described in this paper is available from the author's website (www.uq.edu.au/similar to uqjkeith/) or upon request to the author.

Computational analysis of base composition pattern and promoter elements in the putative promoter regions in relation to expression profiles of 682 human genes on chromosome 22

The base composition pattern (BCP) in the putative promoter region (PPRs) up to 5 Kb lengths of 682 human genes on Chromosome 22 (Chr22) was examined. Two-dimensional (2D) and three-dimensional (3D) functions were designed to delineate the DNA base composition, with four major patterns identified. It is found that 17.6% genes include TATA box, 28.0% GC box, 18.9% CAAT box and 38.4% CpG islands, and approximately 10% genes have one of four putative initiator (Inr) motifs. The occurrence of the promoter elements is tightly associated with the base composition features in the promoter regions, and the associations of the base composition features with occurrence of the promoter elements in the promoter regions mediate tissue-wide expression of the genes in human. The occurrence of two or more promoter elements in the promoter regions is required for the medium- and wide-range expression profiles of the human genes on Chr22. Thus, the reported data shed light on the characteristics of the PPRs of the human genes on Chr22, which may improve our understanding of regulatory roles of the PPRs with occurrence of the promoter elements in gene expression.

Mapping of autosomal recessive chronic distal spinal muscular atrophy to chromosome 11q13

Distal spinal muscular atrophy is a heterogeneous group of neuromuscular disorders caused by progressive anterior born cell degeneration and characterized by progressive motor weakness and muscular atrophy, predominantly in the distal parts of the limbs. Here we report on chronic autosomal recessive distal spinal muscular atrophy in a large, inbred family with onset at various ages. Because this condition had some of the same clinical features as spinal muscular atrophy with respiratory distress, we tested the disease gene for linkage to chromosome 11q and mapped the disease locus to chromosome 11q13 in the genetic interval that included the spinal muscular atrophy with respiratory distress gene (D11S1889-D11S1321, Z(max) = 4.59 at theta = 0 at locus D11S4136). The sequencing of IGHMBP2, the human homologue of the mouse neuromuscular degeneration gene (nmd) that accounts for spinal muscular atrophy with respiratory distress, failed to detect any mutation in our chronic distal spinal muscular atrophy patients, suggesting that spinal muscular atrophy with respiratory distress and chronic distal spinal muscular atrophy are caused by distinct genes located in the so-me chromosomal region. In addition, the high intrafamilial variability in age at onset raises the question of whether nonallelic modifying genes could be involved in chronic distal spinal muscular atrophy.

Expression of SHOX in human fetal and childhood growth plate

Abnormalities in the growth plate may lead to short stature and skeletal deformity including Leri Weil syndrome, which has been shown to result from deletions or mutations in the SHOX gene, a homeobox gene located at the pseudoautosomal region of the X and Y chromosome. We studied the expression of SHOX protein, by immunohistochemistry, in human fetal and childhood growth plates and mRNA by in situ hybridization in childhood normal and Leri Weil growth plate. SHOX protein was found in reserve, proliferative, and hypertrophic zones of fetal growth plate from 12 wk to term and childhood control and Leri Weil growth plates. The pattern of immunostaining in the proliferative zone of childhood growth plate was patchy, with more intense uniform immunostaining in the hypertrophic zone. In situ hybridization studies of childhood growth plate demonstrated SHOX mRNA expression throughout the growth plate. No difference in the pattern of SHOX protein or mRNA expression was seen between the control and Leri Weil growth plate. These findings suggest that SHOX plays a role in chondrocyte function in the growth plate.

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia - and, more recently, for bipolar disorder - on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P = .0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.

Multicenter linkage study of schizophrenia loci on chromosome 22q

The hypothesis of the existence of one or more schizophrenia susceptibility loci on chromosome 22q is supported by reports of genetic linkage and association, meta-analyses of linkage, and the observation of elevated risk for psychosis in people with velocardiofacial syndrome, caused by 22q11 microdeletions. We tested this hypothesis by evaluating 10 microsatellite markers spanning 22q in a multicenter sample of 779 pedigrees. We also incorporated age at onset and sex into the analysis as covariates. No significant evidence for linkage to schizophrenia or for linkage associated with earlier age at onset, gender, or heterogeneity across sites was observed. We interpret these findings to mean that the population-wide effects of putative 22q schizophrenia susceptibility loci are too weak to detect with linkage analysis even in large samples.

Human SULT1A genes: Cloning and activity assays of the SULT1A promoters

The three human SULT1A sulfotransferase enzymes are closely related in amino acid sequence (>90%), yet differ in their substrate preference and tissue distribution. SULT1A1 has a broad tissue distribution and metabolizes a range of xenobiotics as well as endogenous substrates such as estrogens and iodothyronines. While the localization of SULT1A2 is poorly understood, it has been shown to metabolize a number of aromatic amines. SULT1A3 is the major catecholamine sulfonating form, which is consistent with it being expressed principally in the gastrointestinal tract. SULT1A proteins are encoded by three separate genes, located in close proximity to each other on chromosome 16. The presence of differential 5′-untranslated regions identified upon cloning of the SULT1A cDNAs suggested the utilization of differential transcriptional start sites and/or differential splicing. This chapter describes the methods utilized by our laboratory to clone and assay the activity of the promoters flanking these different untranslated regions found on SULT1A genes. These techniques will assist investigators in further elucidating the differential mechanisms that control regulation of the human SULT1A genes. They will also help reveal how different cellular environments and polymorphisms affect the activity of SULT1A gene promoters.

Familial aggregation of human infection with Schistosoma japonicum in the Poyang Lake region, China

Despite the success of extensive control measures that have been implemented in China for over 50 years, the number of individuals infected with Schistosoma japonicum remains high in the existing endemic areas. A variance components analysis was undertaken to estimate the heritable and environmental components that contribute to S. japonicum infection in the Poyang Lake region of Jiangxi Province, PR China. The total target population was 3148 from four separate administrative villages. Two thousand seven hundred and five of these comprised 400 families ranging in size from 3 to 188. After adjustments were made for gender, water contact and past history of having had schistosomiasis, the heritable component was estimated to account for as much as 58% of the phenotype variation under the polygenic model. Household was not shown to be an important environmental factor. Incorporating village effects indicated that the results were valid for the total population. We conclude that genetic heritability in this region is high and plays an important role in determining risk of infection with S. japonicum. (c) 2005 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.