The novel imprinted carboxypeptidase A4 gene ( CPA4) in the 7q32 imprinting domain

By a search for novel human imprinted genes in the vicinity of the imprinted gene MEST, at chromosome 7q32, we identified the carboxypeptidase A4 gene ( CPA4) in a gene cluster of the carboxypeptidase family, 200 kb centromeric to MEST. Because CPA4 was originally identified as a protein induced in a prostate cancer cell line (PC-3) by histone deacetylase inhibitors, and was located at the putative prostate cancer-aggressiveness locus at 7q32, we investigated its imprinting status in fetal tissues and in adult benign hypertrophic prostate (BPH). RT-PCR using four intragenic polymorphisms as markers showed that CPA4 was expressed preferentially from the maternal allele in the fetal heart, lung, liver, intestine, kidney, adrenal gland, and spleen, but not in the fetal brain. It was also preferentially expressed in the BPH. These findings support that CPA4 is imprinted and may become a strong candidate gene for prostate cancer-aggressiveness. As a Silver-Russell syndrome (SRS) locus has been proposed to be located to a region near MEST and to be involved in imprinting, CPA4 would have been a candidate gene for SRS. However, analysis of ten SRS patients revealed no mutations in CPA4.     

Detection of closely linked multiple quantitative trait loci using a genetic algorithm

The existence of a quantitative trait locus (QTL) is usually tested using the likelihood of the quantitative trait on the basis of phenotypic character data plus the recombination fraction between QTL and flanking markers. When doing this, the likelihood is calculated for all possible locations on the linkage map. When multiple QTL are suspected close by, it is impractical to calculate the likelihood for all possible combinations of numbers and locations of QTL. Here, we propose a genetic algorithm (GA) for the heuristic solution of this problem. GA can globally search the optimum by improving the "genotype" with alterations called "recombination" and "mutation." The "genotype" of our GA is the number and location of QTL. The "fitness" is a function based on the likelihood plus Akaike's information criterion (AIC), which helps avoid false-positive QTL. A simulation study comparing the new method with existing QTL mapping packages shows the advantage of the new GA. The GA reliably distinguishes multiple QTL located in a single marker interval.     

Multidimensional vector space representation for convergent evolution and molecular phylogeny

With growing amounts of genome data and constant improvement of models of molecular evolution, phylogenetic reconstruction became more reliable. However, our knowledge of the real process of molecular evolution is still limited. When enough large-sized data sets are analyzed, any subtle biases in statistical models can support incorrect topologies significantly because of the high signal-to-noise ratio. We propose a procedure to locate sequences in a multidimensional vector space (MVS), in which the geometry of the space is uniquely determined in such a way that the vectors of sequence evolution are orthogonal among different branches. In this paper, the MVS approach is developed to detect and remove biases in models of molecular evolution caused by unrecognized convergent evolution among lineages or unexpected patterns of substitutions. Biases in the estimated pairwise distances are identified as deviations (outliers) of sequence spatial vectors from the expected orthogonality. Modifications to the estimated distances are made by minimizing an index to quantify the deviations. In this way, it becomes possible to reconstruct the phylogenetic tree, taking account of possible biases in the model of molecular evolution. The efficacy of the modification procedure was verified by simulating evolution on various topologies with rate heterogeneity and convergent change. The phylogeny of placental mammals in previous analyses of large data sets has varied according to the genes being analyzed. Systematic deviations caused by convergent evolution were detected by our procedure in all representative data sets and were found to strongly affect the tree structure. However, the bias correction yielded a consistent topology among data sets. The existence of strong biases was validated by examining the sites of convergent evolution between the hedgehog and other species in mitochondrial data set. This convergent evolution explains why it has been difficult to determine the phylogenetic placement of the hedgehog in previous studies.     

Divergence pattern of duplicate genes in protein-protein interactions follows the power law

The impact of the biological network structures on the divergence between the two copies of one duplicate gene pair involved in the networks has not been documented on a genome scale. Having analyzed the most recently updated Database of Interacting Proteins (DIP) by incorporating the information for duplicate genes of the same age in yeast, we find that there was a highly significantly positive correlation between the level of connectivity of ancient genes and the number of shared partners of their duplicates in the protein-protein interaction networks. This suggests that duplicate genes with a low ancestral connectivity tend to provide raw materials for functional novelty, whereas those duplicate genes with a high ancestral connectivity tend to create functional redundancy for a genome during the same evolutionary period. Moreover, the difference in the number of partners between two copies of a duplicate pair was found to follow a power-law distribution. This suggests that loss and gain of interacting partners for most duplicate genes with a lower level of ancestral connectivity is largely symmetrical, whereas the "hub duplicate genes" with a higher level of ancient connectivity display an asymmetrical divergence pattern in protein-protein interactions. Thus, it is clear that the protein-protein interaction network structures affect the divergence pattern of duplicate genes. Our findings also provide insights into the origin and development of biological networks.     

Estimating absolute rates of synonymous and nonsynonymous nucleotide substitution in order to characterize natural selection and date species divergences

The rate of molecular evolution can vary among lineages. Sources of this variation have differential effects on synonymous and nonsynonymous substitution rates. Changes in effective population size or patterns of natural selection will mainly alter nonsynonymous substitution rates. Changes in generation length or mutation rates are likely to have an impact on both synonymous and nonsynonymous substitution rates. By comparing changes in synonymous and nonsynonymous rates, the relative contributions of the driving forces of evolution can be better characterized. Here, we introduce a procedure for estimating the chronological rates of synonymous and nonsynonymous substitutions on the branches of an evolutionary tree. Because the widely used ratio of nonsynonymous and synonymous rates is not designed to detect simultaneous increases or simultaneous decreases in synonymous and nonsynonymous rates, the estimation of these rates rather than their ratio can improve characterization of the evolutionary process. With our Bayesian approach, we analyze cytochrome oxidase subunit I evolution in primates and infer that nonsynonymous rates have a greater tendency to change over time than do synonymous rates. Our analysis of these data also suggests that rates have been positively correlated.     

Simultaneous detection of linkage disequilibrium and genetic differentiation of subdivided populations

We propose a new method for simultaneously detecting linkage disequilibrium and genetic structure in subdivided populations. Taking subpopulation structure into account with a hierarchical model, we estimate the magnitude of genetic differentiation and linkage disequilibrium in a metapopulation on the basis of geographical samples, rather than decompose a population into a finite number of random-mating subpopulations. We assume that Hardy-Weinberg equilibrium is satisfied in each locality, but do not assume independence between marker loci. Linkage states remain unknown. Genetic differentiation and linkage disequilibrium are expressed as hyperparameters describing the prior distribution of genotypes or haplotypes. We estimate related parameters by maximizing marginal-likelihood functions and detect linkage equilibrium or disequilibrium by the Akaike information criterion. Our empirical Bayesian model analyzes genotype and haplotype frequencies regardless of haploid or diploid data, so it can be applied to most commonly used genetic markers. The performance of our procedure is examined via numerical simulations in comparison with classical procedures. Finally, we analyze isozyme data of ayu, a severely exploited fish species, and single-nucleotide polymorphisms in human ALDH2.     

Ricinoleic acid and castor oil as substrates for conjugated linoleic acid production by washed cells of Lactobacillus plantarum

Ricinoleic acid (12-hydroxy-cis-9-octadecaenoic acid) was an effective substrate for conjugated linoleic acid (CLA) production by washed cells of Lactobacillus plantarum AKU 1009a. The CLA produced was a mixture of cis-9,trans-11- and trans-9,trans-11-octadecadienoic acids. Addition of alpha-linolenic acid to the culture medium increased the CLA productivity of the washed cells. In the presence of lipase, castor oil, in which the main fatty acid component is ricinoleic acid, also was a substrate for CLA.     

Characteristics of45Ca uptake stimulated by high KCl of differentiated and undifferentiated NG108-15 and PC12h cells

The characteristics of KCl-stimulated⁴⁵Ca uptake by neuroblastoma x glioma hybrid NG108-15 cells induced to differentiate with dibutyryl cAMP (Bt₂cAMP) and of PC12h pheochromocytoma cells induced to differentiate with nerve growth factor (NGF) were studied. The extent and rate of KCl-stimulated⁴⁵Ca uptake by differentiated NG108-15 cells induced with Bt₂cAMP were significantly higher than those of the undifferentiated cells. However, differentiation of PC12h cells induced with NGF did not enhance their extent or rate of KCl-stimulated⁴⁵Ca uptake. The effects of Ca agonist and antagonists indicated that the characteristics of KCl-stimulated⁴⁵Ca uptake by Bt₂cAMP-treated NG108-15 cells and NGF-treated PC12h cells mainly reflected those of peripheral L-type voltage-sensitive calcium channels activated by high KCl. These results suggest that differentiated neural cells did not all show an enhanced capacity for KCl-stimulated⁴⁵Ca uptake, although the characteristic patterns of differentiation (extension of neurite-like processes, etc.) and that of effect by Ca agonist or antagonists on NG108-15 cells and PC12h cells were similar.