Refined human artificial chromosome vectors for gene therapy and animal transgenesis.

Human artificial chromosomes (HACs) have several advantages as gene therapy vectors, including stable episomal maintenance, and the ability to carry large gene inserts. We previously developed HAC vectors from the normal human chromosomes using a chromosome engineering technique. However, endogenous genes were remained in these HACs, limiting their therapeutic applications. In this study, we refined a HAC vector without endogenous genes from human chromosome 21 in homologous recombination-proficient chicken DT40 cells. The HAC was physically characterized using a transformation-associated recombination (TAR) cloning strategy followed by sequencing of TAR-bacterial artificial chromosome clones. No endogenous genes were remained in the HAC. We demonstrated that any desired gene can be cloned into the HAC using the Cre-loxP system in Chinese hamster ovary cells, or a homologous recombination system in DT40 cells. The HAC can be efficiently transferred to other type of cells including mouse ES cells via microcell-mediated chromosome transfer. The transferred HAC was stably maintained in vitro and in vivo. Furthermore, tumor cells containing a HAC carrying the suicide gene, herpes simplex virus thymidine kinase (HSV-TK), were selectively killed by ganciclovir in vitro and in vivo. Thus, this novel HAC vector may be useful not only for gene and cell therapy, but also for animal transgenesis.

Estimation of the species-specific mutation rates at the DRB1 locus in humans and chimpanzee.

To estimate the species-specific mutation rates at the DRB1 locus in humans and chimpanzee, we analyzed the nucleotide sequence of a 37.6-kb chimpanzee chromosomal segment containing the entire Patr-DRB1*0701 allele and the flanking nongenic region and we compared it with two corresponding human sequences containing the HLA-DRB1*070101 allele using the sequence of HLA-DRB1*04011 as an outgroup. Because the allelic pair of HLA-DRB1*070101 and Patr-DRB1*0701 shows the lowest number of substitutions between the two species, it appears that these sequences diverged close to the time of the humans-chimpanzee divergence (6 million years ago). Alignment of the nucleotide sequences for HLA-DRB1*070101 and Patr-DRB1*0701 alleles showed that they share a high degree of similarity, suggesting that the studied chromosomal segments with these sequences have not been subjected to recombination since the humans-chimpanzee divergence. Comparison of the flanking 10.6 kb of nongenic sequences revealed an average of 41.5 and 83 single nucleotide substitutions in humans and chimpanzee, respectively. Thus, the species-specific nucleotide substitution rates in the flanking nongenic region were estimated to be 6.53 x 10(-10) and 1.31 x 10(-9) per site per year in humans and chimpanzee, respectively. Unexpectedly, the estimated rate in humans was twofold lower than in chimpanzee (P < 10(-3), Tajima's relative rate test) and lower than the average substitution rate in the human genome. Because the nucleotide substitution rate in nongenic regions free from selection is expected to be equal to the mutation rate, the estimated substitution rate should correspond to the species-specific mutation rate at the DRB1 locus. Our results strongly suggest that the mutation rate at DRB1 locus differs among species.

Evaluation of RANK/RANKL/OPG gene polymorphisms in aggressive periodontitis.

BACKGROUND AND OBJECTIVE: Aggressive periodontitis (AgP) is a specific type of periodontal disease that is characterized by rapid attachment loss and bone destruction. While attempting to identify genetic polymorphisms associated with AgP, previous research has focused on candidate genes that may be involved in immune responses to microbial infections. In this study, the focus was on single nucleotide polymorphisms (SNPs) in the key mediators of osteoclast differentiation and activation, which involve receptor activator of nuclear factor-kappaB (RANK), RANK ligand (RANKL) and osteoprotegrin (OPG), in the Japanese population. The aim of this study was to evaluate the association of RANK/RANKL/OPG gene polymorphisms with AgP in the Japanese population. MATERIAL AND METHODS: We examined 99 patients with AgP and 89 controls from the Japanese population to explore the possibility of RANK/RANKL/OPG loci as candidate regions associated with the disease. All exons and relevant exon-intron boundaries of these three candidate genes were amplified by polymerase chain reaction (PCR) using 19 primers, followed by direct sequencing. The polymorphisms were identified by comparing the sequences obtained from 48 subjects. RESULTS: We identified 27 SNPs in RANK, including 10 novel SNPs and seven SNPs each in both RANKL and OPG. A pairwise linkage disequilibrium analysis using the r2 statistic showed that some SNP pairs from the three loci are in tight linkage disequilibrium. CONCLUSION: An association analysis with allelotypes showed that SNPs identified in the RANK/RANKL/OPG genes have no significant association with AgP in the Japanese population.

CSF1 gene associated with aggressive periodontitis in the Japanese population.

Aggressive periodontitis (AgP) is characterized by the early onset of the rapid and progressive destruction of the alveolar bone. We investigated the correlation of single nucleotide polymorphisms (SNPs) in candidate genes with AgP in the Japanese population in order to determine the genetic risk factors for this complex disease. Among 11 genes related to bone formation and resorption, 43 known SNPs were tested in 98 case and 88 control samples for association with AgP by using SNP genotyping techniques. Among these, three polymorphisms located in the colony stimulating factor 1 (CSF1) gene showed a positive association with AgP. This is the first case of an association between a CSF1 polymorphism and a human disease.

FLT3 is fused to ETV6 in a myeloproliferative disorder with hypereosinophilia and a t(12;13)(p13;q12) translocation.

The FMS-like tyrosine kinase 3 (FLT3) gene, belonging to the receptor tyrosine kinase (TK) subclass III family, plays an important role in normal hematopoiesis and is one of the most frequently mutated genes in hematologic malignancies as well as an attractive target for directed inhibition. Activating mutations of this gene, including internal tandem duplication in the juxtamembrane (JM) domain and point mutations in the TK domain, are found in approximately one-third of patients with acute myeloid leukemia and in a smaller subset of patients with acute lymphoblastic leukemia. We report here that FLT3 may contribute to leukemogenesis in a patient with myeloproliferative disorder and a t(12;13)(p13;q12) translocation through generating a fusion gene with the ETS variant gene 6 (ETV6) gene. ETV6 has been reported to fuse to various partner genes, including TK and transcription factors. Both ETV6/FLT3 and reciprocal FLT3/ETV6 transcripts were detected in the patient mRNA by reverse transcriptase-polymerase chain reaction. At the protein level, however, only ETV6/FLT3 products were expressed. Among them, one retains the helix-loop-helix (HLH) oligomerization domain of ETV6 and the JM as well as TK domain of FLT3. FLT3 receptor in leukemic cells might be inappropriately activated through dimerization by HLH domain of ETV6, which consequently interfered with proliferation and differentiation of hematopoietic cells.

DNA sequence and comparative analysis of chimpanzee chromosome 22.

Human-chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.

DIGIT: a novel gene finding program by combining gene-finders.

We have developed a general purpose algorithm which finds genes by combining plural existing gene-finders. The algorithm has been implemented into a novel gene-finder named DIGIT. An outline of the algorithm is as follows. First, existing gene-finders are applied to an uncharacterized genomic sequence (input sequence). Next, DIGIT produces all possible exons from the results of gene-finders, and assigns them their exon types, reading frames and exon scores. Finally, DIGIT searches a set of exons whose additive score is maximized under their reading frame constraints. Bayesian procedure and a hidden Markov model are used to infer exon scores and search the exon set, respectively. We have designed DIGIT so as to combine the results of FGENESH, GENSCAN and HMMgene, and have assessed its prediction accuracy by using recently compiled benchmark data sets. For all data sets, DIGIT successfully discarded many false-positive exons predicted by individual gene-finders and yielded remarkable improvements in sensitivity and specificity at the gene level compared with the best gene level accuracies achieved by any single gene-finder.

Accumulation of species-specific amino acid replacements that cause loss of particular protein functions in Buchnera, an endocellular bacterial symbiont.

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected.

Genome sequence of an industrial microorganism Streptomyces avermitilis: deducing the ability of producing secondary metabolites.

Streptomyces avermitilis is a soil bacterium that carries out not only a complex morphological differentiation but also the production of secondary metabolites, one of which, avermectin, is commercially important in human and veterinary medicine. The major interest in this genus Streptomyces is the diversity of its production of secondary metabolites as an industrial microorganism. A major factor in its prominence as a producer of the variety of secondary metabolites is its possession of several metabolic pathways for biosynthesis. Here we report sequence analysis of S. avermitilis, covering 99% of its genome. At least 8.7 million base pairs exist in the linear chromosome; this is the largest bacterial genome sequence, and it provides insights into the intrinsic diversity of the production of the secondary metabolites of Streptomyces. Twenty-five kinds of secondary metabolite gene clusters were found in the genome of S. avermitilis. Four of them are concerned with the biosyntheses of melanin pigments, in which two clusters encode tyrosinase and its cofactor, another two encode an ochronotic pigment derived from homogentiginic acid, and another polyketide-derived melanin. The gene clusters for carotenoid and siderophore biosyntheses are composed of seven and five genes, respectively. There are eight kinds of gene clusters for type-I polyketide compound biosyntheses, and two clusters are involved in the biosyntheses of type-II polyketide-derived compounds. Furthermore, a polyketide synthase that resembles phloroglucinol synthase was detected. Eight clusters are involved in the biosyntheses of peptide compounds that are synthesized by nonribosomal peptide synthetases. These secondary metabolite clusters are widely located in the genome but half of them are near both ends of the genome. The total length of these clusters occupies about 6.4% of the genome.

Differential display analysis of mutants for the transcription factor Pdr1p regulating multidrug resistance in the budding yeast.

The transcription factor Pdr1p recognizes Pdr1p/Pdr3p-response element (PDRE) to activate genes involved in multidrug resistance of the budding yeast. To identify novel targets of Pdr1p, we compared transcriptomes among the yeast cells bearing wild, disrupted and gain-of-function alleles of PDR1 using a high-throughput fluorescent differential display PCR. Consequently, we identified 20 transcripts apparently regulated by Pdr1p, which are derived from well-known target genes as well as those that have never been described in the context of drug resistance. Intriguingly, among the latter, a previously unrecognized gene bearing a small putative open reading frame preceded by a functional PDRE was found.

Multiplex polymerase chain reaction (PCR) with color-tagged module-shuffling primers for comparing gene expression levels in various cells.

A method based on the multiplex polymerase chain reaction (PCR) and gel electrophoresis for the comparative analysis of gene expression levels was developed. Using the method many cDNA fragments from different sources can be compared simultaneously. Competitive PCR amplification of expressed genes from different sources was performed by using 'module-shuffling primers' (MPSs). The MPSs (labeled with different fluorophores) consist of sequence modules of 3 or 4 nt. The modules are arranged in different orders in each primer; therefore, the base sequences of the primers are different but their melting temperatures are identical. The genes expressed in different sources are ligated with tags complementary with the MPSs. Tag-ligated fragments are mixed in one tube and amplified at the same amplification efficiency by the MPSs. Amplified fragments are detected separately by multiple-color gel electrophoresis. This method can detect different amounts of each expressed gene, up to a difference in amounts of 30%, and its detection limit is 0.1 amol per assay.

Discovery of a novel compound: insight into mechanisms for acrylamide-induced axonopathy and colchicine-induced apoptotic neuronal cell death.

The exposure of humans and experimental animals to certain industrial toxins such as acrylamide is known to cause nerve damage classified as axonopathy, but the mechanisms involved are poorly understood. Here we show that acrylamide induces morphological changes and tyrosine phosphorylation of focal adhesion kinase (FAK) and proline-rich tyrosine kinase 2 (Pyk2), a member of the FAK subfamily, in human differentiating neuroblastoma SH-SY5Y cells. Furthermore, we identified a novel molecule designated 'compound-1' that inhibits the morphological and biochemical events. Daily oral administrations of the compound also effectively alleviated behavioral deficits in animals elicited by acrylamide in inclined plane testing, landing foot spread testing and rota-rod performance testing. The compound also effectively inhibited the biological and biochemical responses caused by another axonopathy inducer, colchicine, including tyrosine phosphorylation of Pyk2, formation of an 85-kDa poly(ADP-ribose)polymerase (PARP) fragment and apoptosis-associated induction of the NAPOR gene as well as neuronal cell death. Our findings not only provide insight into FAK and Pyk2 functions in neuronal cells, but may also be important in the development of therapeutic agents for peripheral neuropathy and neurodegeneration.

The genomic structure and expression of MJD, the Machado-Joseph disease gene.

Machado-Joseph disease (MJD) is an autosomal dominant neurodegenerative disorder that is clinically characterized by cerebellar ataxia and various associated symptoms. The disease is caused by an unstable expansion of the CAG repeat in the MJD gene. This gene is mapped to chromosome 14q32.1. To determine its genomic structure, we constructed a contig composed of six cosmid clones and eight bacterial artificial chromosome (BAC) clones. It spans approximately 300kb and includes MJD. We also determined the complete sequence (175,330bp) of B445M7, a human BAC clone that contains MJD. The MJD gene was found to span 48,240bp and to contain 11 exons. Northern blot analysis showed that MJD mRNA is ubiquitously expressed in human tissues, and in at least four different sizes; namely, 1.4, 1.8, 4.5, and 7.5kb. These different mRNA species probably result from differential splicing and polyadenylation, as shown by sequences of the 21 independent cDNA clones isolated after the screening of four human cDNA libraries prepared from whole brain, caudate, retina, and testis. The sequences of these latter clones relative to the MJD gene in B445M7 indicate that there are three alternative splicing sites and eight polyadenylation signals in MJD that are used to generate the differently sized transcripts.

Isolation of hNap1BP which interacts with human Nap1 (NCKAP1) whose expression is down-regulated in Alzheimer's disease.

We previously reported the isolation of a novel apoptosis-related gene, human Nap1 (HGMW-approved symbol NCKAP1), the expression of which was strongly down-regulated in sporadic Alzheimer's disease (AD). Human Nap1 proved to be an orthologue of rat Nap1 which binds to the adaptor molecule Nck in signal transduction. In order to further elucidate the function of human Nap1, we performed yeast two-hybrid screening. As a result of screening, we discovered a protein designated hNap1BP (human Nap1 binding protein) which is a member of the tyrosine kinase-binding protein family. In addition, hNap1BP bound to the SH3 domain of c-Abl and Nck. hNap1BP is expressed ubiquitously in various tissues like human Nap1, and intriguingly these genes are co-expressed in hippocampus and cerebral cortex in mouse brain where AD pathological features are strongly evident. Further functional analysis of hNap1BP may clarify its contribution to AD pathology.

Diverse transcriptional initiation revealed by fine, large-scale mapping of mRNA start sites.

Determination of the mRNA start site is the first step in identifying the promoter region, which is of key importance for transcriptional regulation of gene expression. The 'oligo-capping' method enabled us to introduce a sequence tag to the first base of an mRNA by replacing the cap structure of the mRNA. Using cDNA libraries made from oligo-capped mRNAs, we could identify the transcriptional start site of an individual mRNA just by sequencing the 5'-end of the cDNA. The fine mapping of transcriptional start sites was performed for 5880 mRNAs in 276 human genes. Contrary to our expectations, the majority of the genes showed a diverse distribution of transcriptional start sites. They were distributed over 61.7 bp with a standard deviation of 19.5. Our finding may reflect the dynamic nature of transcriptional initiation events of human genes in vivo.

Puromycin-sensitive aminopeptidase is essential for the maternal recognition of pregnancy in mice.

Maternal recognition of pregnancy in rodents requires semicircadian surges of hypophyseal PRL secretion during early gestation, which are required for the formation of the corpus luteum of pregnancy (CLP). Here we show that puromycin-sensitive aminopeptidase (Psa)-deficient mice display female infertility that results from impaired formation of CLP. Transplantation of mutant ovaries into normal females restored fertility but not vice versa. Psa-deficient females revealed no semicircadian surges of PRL induced after mating stimuli. Pregnancy in the mutant females was restored by grafting intact pituitaries to elevate circulating levels of PRL. Psa is thus required for the appearance of the semicircadian surges of PRL secretion that are crucial for maintaining pregnancy in rodents.

Male reproductive defects caused by puromycin-sensitive aminopeptidase deficiency in mice.

Male reproductive performance is composed of two principal elements, copulation and spermatogenesis. A wealth of literature has described the intricate web of endocrine events underlying these biological processes. In the present study we show that puromycin-sensitive aminopeptidase (Psa)-deficient mice are infertile, lack copulatory behavior, and have impaired spermatogenesis. The reproductive deficits of the mutants are not restored by androgen administration, although no aberrant localization of the sex steroid receptors was detectable in their brains and testes. Considering the strong expression of the Psa gene in the brain and Sertoli cells and the degenerative morphology of Sertoli cells in Psa-deficient mice, Psa may participate in testosterone-mediated reproductive signal pathways in the brain and testis.

Identification and characterization of the potential promoter regions of 1031 kinds of human genes.

To understand the mechanism of transcriptional regulation, it is essential to identify and characterize the promoter, which is located proximal to the mRNA start site. To identify the promoters from the large volumes of genomic sequences, we used mRNA start sites determined by a large-scale sequencing of the cDNA libraries constructed by the "oligo-capping" method. We aligned the mRNA start sites with the genomic sequences and retrieved adjacent sequences as potential promoter regions (PPRs) for 1031 genes. The PPR sequences were searched to determine the frequencies of major promoter elements. Among 1031 PPRs, 329 (32%) contained TATA boxes, 872 (85%) contained initiators, 999 (97%) contained GC box, and 663 (64%) contained CAAT box. Furthermore, 493 (48%) PPRs were located in CpG islands. This frequency of CpG islands was reduced in TATA(+)/Inr(+) PPRs and in the PPRs of ubiquitously expressed genes. In the PPRs of the CGM2 gene, the DRA gene, and the TM30pl genes, which showed highly colon specific expression patterns, the consensus sequences of E boxes were commonly observed. The PPRs were also useful for exploring promoter SNPs.

Molecular cloning and characterization of a gene expressed in mouse developing tongue, mDscr5 gene, a homolog of human DSCR5 (Down syndrome Critical Region gene 5).

For understanding the pathogenesis of Down syndrome (DS), it is important to identify and characterize the genes on Chromosome (Chr) 21, especially those in the Down syndrome critical region (DSCR) on Chr 21q22.2. Recently we have determined 33.5 Mb (more than 99%) of DNA sequence of Chr 21 and, from these sequence data, we identified a novel gene, DSCR5 (transcript = 0.8 kb), from DSCR by combination of computational gene prediction and cDNA screening. For functional analysis of DSCR5, we identified a mouse homolog of the DSCR5 cDNA, and termed it mDscr5 (transcript length = 0.8 kb). The gene was mapped to mouse Chr 16 C3-C4, the syntenic region of human Chr 21, and encodes an amino acid of 132 residues with 90% identity to DSCR5. In situ hybridization showed that mDscr5 is predominantly expressed in the developing tongue. To our best knowledge, no other gene in DSCR is reported to be expressed in tongue, so that DSCR5 may be the first candidate to elucidate the pathophysiology of tongue malformation observed in DS.