Search

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase ( http://silkbase.ab.a.u-tokyo.ac.jp ).

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister–taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault’s rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Heterochromatin is marked by methylation of lysine 9 on histone H3 (H3K9me). A puzzling feature of H3K9me is that this modification localizes not only in promoters but also in internal regions (bodies) of silent transcription units. Despite its prevalence, the biological significance of gene‐body H3K9me remains enigmatic. Here we show that H3K9me‐associated removal of H3K4 monomethylation (H3K4me1) in gene bodies mediates transcriptional silencing. Mutations in an Arabidopsis H3K9 demethylase gene IBM1 induce ectopic H3K9me2 accumulation in gene bodies, with accompanying severe developmental defects. Through suppressor screening of the ibm1 ‐induced developmental defects, we identified the LDL2 gene, which encodes a homolog of conserved H3K4 demethylases. The ldl2 mutation suppressed the developmental defects, without suppressing the ibm1 ‐induced ectopic H3K9me2. The ectopic H3K9me2 mark directed removal of gene‐body H3K4me1 and caused transcriptional repression in an LDL2‐dependent manner. Furthermore, mutations of H3K9 methylases increased the level of H3K4me1 in the gene bodies of various transposable elements, and this H3K4me1 increase is a prerequisite for their transcriptional derepression. Our results uncover an unexpected role of gene‐body H3K9me2/H3K4me1 dynamics as a mediator of heterochromatin silencing and epigenome differentiation. ![][1] A genetic interaction between histone demethylases IBM1 and LDL2 reveals that H3K9 methylation in gene bodies induces transcriptional silencing by triggering the loss of H3K4 monomethylation. [1]: /embed/graphic-1.gif

Sex-determining regions (SDRs) or mating-type (MT) loci in two sequenced volvocine algal species, Chlamydomonas reinhardtii and Volvox carteri, exhibit major differences in size, structure, gene content, and gametolog differentiation. Understanding the origin of these differences requires investigation of MT loci from related species. Here, we determined the sequences of the minus and plus MT haplotypes of the isogamous 16-celled volvocine alga, Gonium pectorale, which is more closely related to the multicellular V. carteri than to C. reinhardtii. Compared to C. reinhardtii MT, G. pectorale MT is moderately larger in size, and has a less complex structure, with only two major syntenic blocs of collinear gametologs. However, the gametolog content of G. pectorale MT has more overlap with that of V. carteri MT than with C. reinhardtii MT, while the allelic divergence between gametologs in G. pectorale is even lower than that in C. reinhardtii. Three key sex-related genes are conserved in G. pectorale MT: GpMID and GpMTD1 in MT–, and GpFUS1 in MT+. GpFUS1 protein exhibited specific localization at the plus-gametic mating structure, indicating a conserved function in fertilization. Our results suggest that the G. pectorale–V. carteri common ancestral MT experienced at least one major reformation after the split from C. reinhardtii, and that the V. carteri ancestral MT underwent a subsequent expansion and loss of recombination after the divergence from G. pectorale. These data begin to polarize important changes that occurred in volvocine MT loci, and highlight the potential for discontinuous and dynamic evolution in SDRs.

Environmental adaptation is one of the most fundamental features of organisms. Modern genome science has identified some genes associated with adaptive traits of organisms, and has provided insights into environmental adaptation and evolution. However, how genes contribute to adaptive traits and how traits are selected under an environment in the course of evolution remain mostly unclear. To approach these issues, we utilize "Dark-fly", a Drosophila melanogaster line maintained in constant dark conditions for more than 60 years. Our previous analysis identified 220, 000 single nucleotide polymorphisms (SNPs) in the Dark-fly genome, but did not clarify which SNPs of Dark-fly are truly adaptive for living in the dark. We found here that Dark-fly dominated over the wild-type fly in a mixed population under dark conditions, and based on this domination we designed an experiment for genome reselection to identify adaptive genes of Dark-fly. For this experiment, large mixed populations of Dark-fly and the wild-type fly were maintained in light conditions or in dark conditions, and the frequencies of Dark-fly SNPs were compared between these populations across the whole genome. We thereby detected condition-dependent selections toward approximately 6% of the genome. In addition, we observed the time-course trajectory of SNP frequency in the mixed populations through generations 0, 22, and 49, which resulted in notable categorization of the selected SNPs into three types with different combinations of positive and negative selections. Our data provided a list of about 100 strong candidate genes associated with the adaptive traits of Dark-fly., 「暗黒バエ」の暗闇適応に関わるゲノム配列の解析 -ゲノムの視点から環境適応のメカニズムに迫る-. 京都大学プレスリリース. 2016-03-09.

There is increasing evidence that long non-coding RNAs (lncRNAs) are important for normal reproductive development, yet very few lncRNAs have been identified in phalluses so far. Unlike eutherians, phallus development in the marsupial tammar wallaby occurs post-natally, enabling manipulation not possible in eutherians in which differentiation occurs in utero. We treated with sex steroids to determine the effects of androgen and oestrogen on lncRNA expression during phallus development. Hormonal manipulations altered the coding and non-coding gene expression profile of phalluses. We identified several predicted co-regulatory lncRNAs that appear to be co-expressed with the hormone-responsive candidate genes regulating urethral closure and phallus growth, namely IGF1, AR and ESR1. Interestingly, more than 50% of AR-associated coding genes and lncRNAs were also associated with ESR1. In addition, we identified and validated three novel co-regulatory and hormone-responsive lncRNAs: lnc-BMP5, lnc-ZBTB16 and lncRSPO4. Lnc-BMP5 was detected in the urethral epithelium of male phalluses and was downregulated by oestrogen in males. Lnc-ZBTB16 was downregulated by oestrogen treatment in male phalluses at day 50 post-partum (pp). LncRSPO4 was downregulated by adiol treatment in female phalluses but increased in male phalluses after castration. Thus, the expression pattern and hormone responsiveness of these lncRNAs suggests a physiological role in the development of the phallus.

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68–1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow., Author summary When species evolve, reproductive isolation leads to a build-up of differentiation in the genome where genes involved in the process occur. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence accumulates during speciation. However, much of our understanding of stickleback speciation comes from early stage divergence, with relatively few examples from more divergent species pairs that still exchange genes. To address this, we focused on Pacific Ocean and Japan Sea sticklebacks, which co-occur in the Japanese islands. We established that they are the oldest and most divergent known stickleback species pair, that they evolved in the face of gene flow and that this gene flow is still on going. We found introgression is confined to small, localised genomic regions where recombination rate is high. Our results show high divergence can be maintained between species, despite extensive gene flow.

HJURP exhibits a centromere expansion activity that is separable from its CENP-A–binding activity. It is proposed that this centromere expansion activity reflects the functional properties of HJURP, which uses this activity to contribute to the plastic establishment of a centromeric chromatin structure., The CENP-A–specific chaperone HJURP mediates CENP-A deposition at centromeres. The N-terminal region of HJURP is responsible for binding to soluble CENP-A. However, it is unclear whether other regions of HJURP have additional functions for centromere formation and maintenance. In this study, we generated chicken DT40 knockout cell lines and gene replacement constructs for HJURP to assess the additional functions of HJURP in vivo. Our analysis revealed that the middle region of HJURP associates with the Mis18 complex protein M18BP1/KNL2 and that the HJURP-M18BP1 association is required for HJURP function. In addition, on the basis of the analysis of artificial centromeres induced by ectopic HJURP localization, we demonstrate that HJURP exhibits a centromere expansion activity that is separable from its CENP-A–binding activity. We also observed centromere expansion surrounding natural centromeres after HJURP overexpression. We propose that this centromere expansion activity reflects the functional properties of HJURP, which uses this activity to contribute to the plastic establishment of a centromeric chromatin structure.

Defensive and reproductive protein genes undergo rapid evolution. Small, cysteine-rich secreted peptides (CRPs) act as antimicrobial agents and function in plant intercellular signaling and are over-represented among reproductively expressed proteins. Because of their roles in defense, reproduction and development and their presence in multigene families, CRP variation can have major consequences for plant phenotypic and functional diversification. We surveyed the CRP genes of six closely related Oryza genomes comprising Oryza sativa ssp. japonica and ssp. indica, Oryza glaberrima and three accessions of Oryza rufipogon to observe patterns of evolution in these gene families and the effects of variation on their gene expression. These Oryza genomes, like other plant genomes, have accumulated large reservoirs of CRP sequences, comprising 26 groups totaling between 676 and 843 genes, in contrast to antimicrobial CRPs in animal genomes. Despite the close evolutionary relationships between the genomes, we observed rapid changes in number and structure among CRP gene families. Many CRP sequences are in gene clusters generated by local duplications, have undergone rapid turnover and are more likely to be silent or specifically expressed. By contrast, conserved CRP genes are more likely to be highly and broadly expressed. Variable CRP genes created by repeated duplication, gene modification and inactivation can gain new functions and expression patterns in newly evolved gene copies. For the CRP proteins, the process of gain/loss by deletion or duplication at gene clusters seems to be an important mechanism in evolution of the gene families, which also contributes to their expression evolution.

Genetic transformation systems using reporter genes in whole plants have a wide variety of applications for molecular biological study including the visualization of expression patterns of particular genes and intracellular biological phenomena as well as the identification of novel genes. In this study, we assessed co-expression of each three codon-optimized reporter genes and a selectable marker in the nuclear transformation system of whole Pyropia yezoensis, a red marine alga. With the use of an endogenous promoter, both the codon-optimized hygromycin resistance gene and s-glucuronidase gene (PyGUS) were co-expressed in P. yezoensis cells. A high level of GUS activity was observed in 60 % of the individuals in hygromycin-resistant lines. A histochemical GUS assay revealed that the PyGUS reporter gene was stably introduced and expressed throughout the algae's life cycle. In addition, two live cell reporters, humanized cyan fluorescent protein from Anemonia majano and luciferase from Gaussia princeps, were successfully expressed in whole P. yezoensis. The development of this transformation system involving three types of reporter genes provides opportunities for monitoring temporal changes in gene expression and for genetic screening in red marine algae.

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.

Transposable elements (TEs) have a major impact on genome evolution, but they are potentially deleterious, and most of them are silenced by epigenetic mechanisms, such as DNA methylation. Here, we report the characterization of a TE encoding an activity to counteract epigenetic silencing by the host. In Arabidopsis thaliana, we identified a mobile copy of the Mutator-like element (MULE) with degenerated terminal inverted repeats (TIRs). This TE, named Hiun (Hi), is silent in wild-type plants, but it transposes when DNA methylation is abolished. When a Hi transgene was introduced into the wild-type background, it induced excision of the endogenous Hi copy, suggesting that Hi is the autonomously mobile copy. In addition, the transgene induced loss of DNA methylation and transcriptional activation of the endogenous Hi. Most importantly, the trans-activation of Hi depends on a Hi-encoded protein different from the conserved transposase. Proteins related to this anti-silencing factor, which we named VANC, are widespread in the non-TIR MULEs and may have contributed to the recent success of these TEs in natural Arabidopsis populations.

Most indigenous citrus varieties are assumed to be natural hybrids, but their parentage has so far been determined in only a few cases because of their wide genetic diversity and the low transferability of DNA markers. Here we infer the parentage of indigenous citrus varieties using simple sequence repeat and indel markers developed from various citrus genome sequence resources. Parentage tests with 122 known hybrids using the selected DNA markers certify their transferability among those hybrids. Identity tests confirm that most variant strains are selected mutants, but we find four types of kunenbo (Citrus nobilis) and three types of tachibana (Citrus tachibana) for which we suggest different origins. Structure analysis with DNA markers that are in Hardy-Weinberg equilibrium deduce three basic taxa coinciding with the current understanding of citrus ancestors. Genotyping analysis of 101 indigenous citrus varieties with 123 selected DNA markers infers the parentages of 22 indigenous citrus varieties including Satsuma, Temple, and iyo, and single parents of 45 indigenous citrus varieties, including kunenbo, C. ichangensis, and Ichang lemon by allele-sharing and parentage tests. Genotyping analysis of chloroplast and mitochondrial genomes using 11 DNA markers classifies their cytoplasmic genotypes into 18 categories and deduces the combination of seed and pollen parents. Likelihood ratio analysis verifies the inferred parentages with significant scores. The reconstructed genealogy identifies 12 types of varieties consisting of Kishu, kunenbo, yuzu, koji, sour orange, dancy, kobeni mikan, sweet orange, tachibana, Cleopatra, willowleaf mandarin, and pummelo, which have played pivotal roles in the occurrence of these indigenous varieties. The inferred parentage of the indigenous varieties confirms their hybrid origins, as found by recent studies., ミカンの親はどの品種? : 遺伝解析により60種以上のカンキツ類の親子関係が明らかに. 京都大学プレスリリース. 2017-01-13.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Aurantimicrobium minutum type strain KNC T is a planktonic ultramicrobacterium isolated from river water in western Japan. Strain KNC T has an extremely small, streamlined genome of 1,622,386 bp comprising 1,575 protein-coding sequences. The genome annotation suggests that strain KNC T has an actinorhodopsin-based photometabolism.

In humans and mice, there are 11 genes derived from sushi-ichi related retrotransposons, some of which are known to play essential roles in placental development. Interestingly, this family of retrotransposons was thought to exist only in eutherian mammals, indicating their significant contributions to the eutherian evolution, but at least one, PEG10, is conserved between marsupials and eutherians. Here we report a novel sushi-ichi retrotransposon-derived gene, SIRH12, in the tammar wallaby, an Australian marsupial species of the kangaroo family. SIRH12 encodes a protein highly homologous to the sushi-ichi retrotransposon Gag protein in the tammar wallaby, while SIRH12 in the South American short-tailed grey opossum is a pseudogene degenerated by accumulation of multiple nonsense mutations. This suggests that SIRH12 retrotransposition occurred only in the marsupial lineage but acquired and retained some as yet unidentified novel function, at least in the lineage of the tammar wallaby.

MITOPLD is a member of the phospholipase D superfamily proteins conserved among diverse species. Zucchini (Zuc), the Drosophila homolog of MITOPLD, has been implicated in primary biogenesis of Piwi-interacting RNAs (piRNAs). By contrast, MITOPLD has been shown to hydrolyze cardiolipin in the outer membrane of mitochondria to generate phosphatidic acid, which is a signaling molecule. To assess whether the mammalian MITOPLD is involved in piRNA biogenesis, we generated Mitopld mutant mice. The mice display meiotic arrest during spermatogenesis, demethylation and derepression of retrotransposons, and defects in primary piRNA biogenesis. Furthermore, in mutant germ cells, mitochondria and the components of the nuage, a perinuclear structure involved in piRNA biogenesis/function, are mislocalized to regions around the centrosome, suggesting that MITOPLD may be involved in microtubule-dependent localization of mitochondria and these proteins. Our results indicate a conserved role for MITOPLD/Zuc in the piRNA pathway and link mitochondrial membrane metabolism/signaling to small RNA biogenesis.

Shusei Sato, Hideki Hirakawa, Sachiko Isobe, Eigo Fukai, Akiko Watanabe, Midori Kato, Kumiko Kawashima, Chiharu Minami, Akiko Muraki, Naomi Nakazaki, Chika Takahashi, Shinobu Nakayama, Yoshie Kishida, Mitsuyo Kohara, Manabu Yamada, Hisano Tsuruoka, Shigemi Sasamoto, Satoshi Tabata, Tomoyuki Aizu, Atsushi Toyoda, Tadasu Shin-i, Yohei Minakuchi, Yuji Kohara, Asao Fujiyama, Suguru Tsuchimoto, Shin'ichiro Kajiyama, Eri Makigano, Nobuko Ohmido, Nakako Shibagaki, Joyce A. Cartagena, Naoki Wada, Tsutomu Kohinata, Alipour Atefeh, Shota Yuasa, Sachihiro Matsunaga, Kiichi Fukui, Sequence Analysis of the Genome of an Oil-Bearing Tree, Jatropha curcas L., DNA Research, Volume 18, Issue 1, February 2011, Pages 65–76, https://doi.org/10.1093/dnares/dsq030., The whole genome of Jatropha curcas was sequenced, using a combination of the conventional Sanger method and new-generation multiplex sequencing methods. Total length of the non-redundant sequences thus obtained was 285 858 490 bp consisting of 120 586 contigs and 29 831 singlets. They accounted for ∼95% of the gene-containing regions with the average G + C content was 34.3%. A total of 40 929 complete and partial structures of protein encoding genes have been deduced. Comparison with genes of other plant species indicated that 1529 (4%) of the putative protein-encoding genes are specific to the Euphorbiaceae family. A high degree of microsynteny was observed with the genome of castor bean and, to a lesser extent, with those of soybean and Arabidopsis thaliana. In parallel with genome sequencing, cDNAs derived from leaf and callus tissues were subjected to pyrosequencing, and a total of 21 225 unigene data have been generated. Polymorphism analysis using microsatellite markers developed from the genomic sequence data obtained was performed with 12 J. curcas lines collected from various parts of the world to estimate their genetic diversity. The genomic sequence and accompanying information presented here are expected to serve as valuable resources for the acceleration of fundamental and applied research with J. curcas, especially in the fields of environment-related research such as biofuel production. Further information on the genomic sequences and DNA markers is available at http://www.kazusa.or.jp/jatropha/.

We constructed and characterized arrayed bacterial artificial chromosome (BAC) libraries of five Drosophila species (D. melanogaster, D. simulans, D. sechellia, D. auraria, and D. ananassae), which are genetically well characterized in the studies of meiosis, evolution, population genetics, and developmental biology. The BAC libraries comprise 8,000 to 12,500 clones for each species, estimated to cover the most of the genomes. We sequenced both ends of most of these BAC clones with a success rate of 91%. Of these, 53,701 clones consisting of non-repetitive BAC end sequences (BESs) were mapped with reference of the public D. melanogaster genome sequences. The BES mapping estimated that the BAC libraries of D. auraria and D. ananassae covered 47% and 57% of the D. melanogaster genome, respectively, and those of D. melanogaster, D. sechellia, and D. simulans covered 94-97%. The low coverage by BESs of D. auraria and D. ananassae may be due to the high sequence divergence with D. melanogaster. From the comparative BES mapping, 111 possible breakpoints of chromosomal rearrangements were identified in these four species. The breakpoints of the major chromosome rearrangement between D. simulans and D. melanogaster on the third chromosome were determined within 20 kb in 84E and 30 kb in 93E/F. Corresponding breakpoints were also identified in D. sechellia. The BAC clones described here will be an important addition to the Drosophila genomic resources.

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

Soil microbial communities have great potential for bioremediation of recalcitrant aromatic compounds. However, it is unclear which taxa and genes in the communities, and how they contribute to the bioremediation in the polluted soils. To get clues about this fundamental question here, time-course (up to 24 weeks) metagenomic analysis of microbial community in a closed soil microcosm artificially polluted with four aromatic compounds, including phenanthrene, was conducted to investigate the changes in the community structures and gene pools. The pollution led to drastic changes in the community structures and the gene sets for pollutant degradation. Complete degradation of phenanthrene was strongly suggested to occur by the syntrophic metabolism by Mycobacterium and the most proliferating genus, Burkholderia. The community structure at Week 24 (-12 weeks after disappearance of the pollutants) returned to the structure similar to that before pollution. Our time-course metagenomic analysis of phage genes strongly suggested the involvement of the 'kill-the-winner' phenomenon (i.e. phage predation of Burkholderia cells) for the returning of the microbial community structure. The pollution resulted in a decrease in taxonomic diversity and a drastic increase in diversity of gene pools in the communities, showing the functional redundancy and robustness of the communities against chemical disturbance.

A phenanthrene degrader, Mycobacterium sp. EPa45, was isolated from a phenanthrene-degrading consortium. Here, we report the complete genome sequence of EPa45, which has a 6.2-Mb single circular chromosome. We propose a phenanthrene degradation pathway in EPa45 based on the complete genome sequence.

Epigenetic variations of phenotypes, especially those associated with DNA methylation, are often inherited over multiple generations in plants. The active and inactive chromatin states are heritable and can be maintained or even be amplified by positive feedback in a transgenerational manner. However, mechanisms controlling the transgenerational DNA methylation dynamics are largely unknown. As an approach to understand the transgenerational dynamics, we examined long-term effect of impaired DNA methylation in Arabidopsis mutants of the chromatin remodeler gene DDM1 (Decrease in DNA Methylation 1) through whole genome DNA methylation sequencing. The ddm1 mutation induces a drastic decrease in DNA methylation of transposable elements (TEs) and repeats in the initial generation, while also inducing ectopic DNA methylation at hundreds of loci. Unexpectedly, this ectopic methylation can only be seen after repeated self-pollination. The ectopic cytosine methylation is found primarily in the non-CG context and starts from 3’ regions within transcription units and spreads upstream. Remarkably, when chromosomes with reduced DNA methylation were introduced from a ddm1 mutant into a DDM1 wild-type background, the ddm1-derived chromosomes also induced analogous de novo accumulation of DNA methylation in trans. These results lead us to propose a model to explain the transgenerational DNA methylation redistribution by genome-wide negative feedback. The global negative feedback, together with local positive feedback, would ensure robust and balanced differentiation of chromatin states within the genome., Author Summary DNA methylation is important for controlling activity of transposable elements and genes. An intriguing feature of DNA methylation in plants is that its pattern can be inherited over multiple generations at high fidelity in a Mendelian manner. However, mechanisms controlling the trans-generational DNA methylation dynamics are largely unknown. Arabidopsis mutants of a chromatin remodeler gene DDM1 (Decrease in DNA Methylation 1) show drastic reduction of DNA methylation in transposons and repeats, and also show progressive changes in developmental phenotypes during propagation through self-pollination. We now show using whole genome DNA methylation sequencing that upon repeated selfing, the ddm1 mutation induces an ectopic accumulation of DNA methylation at hundreds of loci. Remarkably, even in the wild type background, the analogous de novo increase of DNA methylation can be induced in trans by chromosomes with reduced DNA methylation. Collectively, our findings support a model to explain the transgenerational DNA methylation redistribution by genome-wide negative feedback, which should be important for balanced differentiation of DNA methylation states within the genome.

Despite the recent discovery of interchromosomal telomere length variation, a role for heterogeneity in telomere maintenance has yet to be established. This study aimed to clarify relative telomere length differences between chromosomes. Combined chromosomal sorting and telomeric repeat content analysis in GM130B cells, the relative telomeric repeat content in each chromosome, were calculated. Each chromosome could be isolated except for chromosomes 1 and 2 and chromosomes 9 to 12, which were isolated in a group. Telomere length was correlated with the size of the corresponding chromosome. Concomitant relative telomeric repeat content analysis in each chromosome and terminal restriction fragment analysis using the whole genome revealed that the terminal restriction fragments of each chromosome were heterogenously distributed through the smear of the fragments. This is the first description of an association between telomere length and chromosome size.

De novo microbial genome sequencing reached a turning point with third-generation sequencing (TGS) platforms, and several microbial genomes have been improved by TGS long reads. Bacillus subtilis natto is closely related to the laboratory standard strain B. subtilis Marburg 168, and it has a function in the production of the traditional Japanese fermented food "natto." The B. subtilis natto BEST195 genome was previously sequenced with short reads, but it included some incomplete regions. We resequenced the BEST195 genome using a PacBio RS sequencer, and we successfully obtained a complete genome sequence from one scaffold without any gaps, and we also applied Illumina MiSeq short reads to enhance quality. Compared with the previous BEST195 draft genome and Marburg 168 genome, we found that incomplete regions in the previous genome sequence were attributed to GC-bias and repetitive sequences, and we also identified some novel genes that are found only in the new genome.

Male and female, generally defined based on differences in gamete size and motility, likely have multiple independent origins, appearing to have evolved from isogamous organisms in various eukaryotic lineages. Recent studies of the gamete fusogen GCS1/HAP2 indicate that this protein is deeply conserved across eukaryotes, and its exclusive and/or functional expression generally resides in males or in male homologues. However, little is known regarding the conserved or primitive molecular traits of males and females within eukaryotes. Here, using morphologically indistinguishable isogametes of the colonial volvocine Gonium pectorale , we demonstrated that GCS1 is differently regulated between the sexes. G. pectorale GCS1 molecules in one sex (homologous to male) are transported from the gamete cytoplasm to the protruded fusion site, whereas those of the other sex (females) are quickly degraded within the cytoplasm upon gamete activation. This molecular trait difference might be conserved across various eukaryotic lineages and may represent male and female prototypes originating from a common eukaryotic ancestor.

Oryzias latipes (medaka) has been established as a vertebrate genetic model for more than a century and recently has been rediscovered outside its native Japan. The power of new sequencing methods now makes it possible to reinvigorate medaka genetics, in particular by establishing a near-isogenic panel derived from a single wild population. Here we characterize the genomes of wild medaka catches obtained from a single Southern Japanese population in Kiyosu as a precursor for the establishment of a near-isogenic panel of wild lines. The population is free of significant detrimental population structure and has advantageous linkage disequilibrium properties suitable for the establishment of the proposed panel. Analysis of morphometric traits in five representative inbred strains suggests phenotypic mapping will be feasible in the panel. In addition, high-throughput genome sequencing of these medaka strains confirms their evolutionary relationships on lines of geographic separation and provides further evidence that there has been little significant interbreeding between the Southern and Northern medaka population since the Southern/Northern population split. The sequence data suggest that the Southern Japanese medaka existed as a larger older population that went through a relatively recent bottleneck approximately 10,000 years ago. In addition, we detect patterns of recent positive selection in the Southern population. These data indicate that the genetic structure of the Kiyosu medaka samples is suitable for the establishment of a vertebrate near-isogenic panel and therefore inbreeding of 200 lines based on this population has commenced. Progress of this project can be tracked at http://www.ebi.ac.uk/birney-srv/medaka-ref-panel.

Using a two-dimensional DNA mapping method, we detected four NotI restriction sites that escape chromosome X-specific methylation in humans. Two genes corresponding to two of these sites that lie in the region of Xp11.2-p11.4 were cloned and their properties studied. One of the genes matched a known gene, but the other, termed EXLM1, is novel and is predominantly expressed in cultured lymphocytes and skeletal muscle.

DNA aberrations in human hepatocellular carcinoma (HCC) were studied by two-dimensional DNA electrophoresis analysis. Five intensified and 60 dwindling spots were detected recurrently in the two-dimensional profile which showed about 3000 restriction DNA fragments as distinctive spots. We assigned these aberrant spots to chromosomes, using the chromosome-assigned two-dimensional profile. Four of the five intensified, and 53 of the 60 dwindling spots were given chromosome assignments. Intensified spots were assigned to chromosomes 5, 6, 9 through 12, 16 and 18. Among the dwindling spots, the highest incidence of aberrations was found on chromosome 16, followed by 9 through 12 and chromosome 2. No aberrations were detected in chromosomes 7, 21, 22 or Y.

The human genome is composed of long-range G+C% (GC%) mosaic structures thought to be related to chromosome bands. We previously reported a boundary of megabase-sized GC% mosaic domains at the junction area between major histocompatibility complex (MHC) classes II and III, proposing it as a possible chromosome band boundary. DNA replication timing during the S phase is known to be correlated cytogenetically with chromosome band zones, and thus the band boundaries have been predicted to contain a switch point for DNA replication timing. In this study, to identify to the nucleotide sequence level the replication switch point during the S phase, we determined the precise DNA replication timing for MHC classes II and III, focusing on the junction area. To do this, we used PCR-based quantitation of nascent DNA obtained from synchronized human myeloid leukemia HL60 cells. The replication timing changed precisely in the boundary region with a 2-h difference between the two sides, supporting the prediction that this region may be a chromosome band boundary. We supposed that replication fork movement terminates (pauses) or significantly slows in the switch region, which contains dense Alu clusters; polypurine/polypyrimidine tracts; di-, tri-, or tetranucleotide repeats; and medium-reiteration-frequency sequences. Because the nascent DNA in the switch region was recovered at low efficiency, we investigated whether this region is associated with the nuclear scaffold and found three scaffold-associated regions in and around the switch region.

We have cloned a new gene, pdh1, from genomic DNA of fission yeast Schizosaccharomyces pombe. pdh1 is actively transcribed as 1400-nucleotide mRNA in vegetatively growing cells and can code for a 226 amino acid polypeptide (pdh1p). Computational structural prediction has revealed that the pdh1p is a highly hydrophobic protein with seven transmembrane domains. The prediction has also detected a possible C-kinase phosphorylation site within the longest hydrophilic loop.

The jawless fish occupy an important phylogenetic position for understanding the evolution of body plans, the origin of adaptive immunity and genome evolution in chordates. We describe here the construction of a large-insert bacterial artificial chromosome (BAC) library from the inshore hagfish, Eptatretus burgeri. The BAC library contains 93,978 clones with an average insert size of 100 kb and is estimated to represent threefold genome-equivalent coverage. The library was organized in three-dimensional pools to facilitate screening by PCR. We have screened this library by PCR and isolated several BAC clones; the average number of positive clones was compatible with the estimated genome coverage of the library. This BAC library, constructed for the first time from the jawless fish, should serve as a useful resource for the scientific community.

RodZ interacts with MreB and both factors are required to maintain the rod shape of Escherichia coli. The assembly of MreB into filaments regulates the subcellular arrangement of a group of enzymes that synthesizes the peptidoglycan (PG) layer. However, it is still unknown how polymerization of MreB determines the rod shape of bacterial cells. Regulatory factor(s) are likely to be involved in controlling the function and dynamics of MreB. We isolated suppressor mutations to partially recover the rod shape in rodZ deletion mutants and found that some of the suppressor mutations occurred in mreB. All of the mreB mutations were in or in the vicinity of domain IA of MreB. Those mreB mutations changed the property of MreB filaments in vivo. In addition, suppressor mutations were found in the periplasmic regions in PBP2 and RodA, encoded by mrdA and mrdB genes. Similar to MreB and RodZ, PBP2 and RodA are pivotal to the cell wall elongation process. Thus, we found that mutations in domain IA of MreB and in the periplasmic domain of PBP2 and RodA can restore growth and rod shape to ΔrodZ cells, possibly by changing the requirements of MreB in the process.

Tardigrades are able to tolerate almost complete dehydration by reversibly switching to an ametabolic state. This ability is called anhydrobiosis. In the anhydrobiotic state, tardigrades can withstand various extreme environments including space, but their molecular basis remains largely unknown. Late embryogenesis abundant (LEA) proteins are heat-soluble proteins and can prevent protein-aggregation in dehydrated conditions in other anhydrobiotic organisms, but their relevance to tardigrade anhydrobiosis is not clarified. In this study, we focused on the heat-soluble property characteristic of LEA proteins and conducted heat-soluble proteomics using an anhydrobiotic tardigrade. Our heat-soluble proteomics identified five abundant heat-soluble proteins. All of them showed no sequence similarity with LEA proteins and formed two novel protein families with distinct subcellular localizations. We named them Cytoplasmic Abundant Heat Soluble (CAHS) and Secretory Abundant Heat Soluble (SAHS) protein families, according to their localization. Both protein families were conserved among tardigrades, but not found in other phyla. Although CAHS protein was intrinsically unstructured and SAHS protein was rich in β-structure in the hydrated condition, proteins in both families changed their conformation to an α-helical structure in water-deficient conditions as LEA proteins do. Two conserved repeats of 19-mer motifs in CAHS proteins were capable to form amphiphilic stripes in α-helices, suggesting their roles as molecular shield in water-deficient condition, though charge distribution pattern in α-helices were different between CAHS and LEA proteins. Tardigrades might have evolved novel protein families with a heat-soluble property and this study revealed a novel repertoire of major heat-soluble proteins in these anhydrobiotic animals.

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.

Posttranslational processing/modification is required for membrane localization and activation of ras proteins. In the case of yeast RAS2 protein, we have reported that the process starts with the removal of the initiator methionine followed by polyisoprenylation, removal of 3 amino acid residues from the C terminus, methyl esterification, and fatty acid acylation (Fujiyama, A., and Tamanoi, F. (1990) J. Biol. Chem. 265, 3362-3368). In this study, we demonstrate that polyisoprenylation and methyl esterification of the cysteine residue in the C-terminal domain of the RAS2 protein are involved in the conversion process from precursor form to intermediate form. The polyisoprenoid moiety attached to the RAS2 protein was identified as a 15-carbon farnesyl group through two independent experiments: the release of S-farnesylcysteine with carboxypeptidase Y from the RAS2 protein, and the recovery of radioactive farnesol through methyliodide treatment of the RAS2 protein purified from yeast cells labeled with [3H]mevalonic acid. The farnesyl group attached to the RAS2 protein was detected predominantly in the C-terminal peptide, SGSGGCC, both in the intermediate and in the fatty acid acylated RAS2 protein. The C-terminal cysteine of the intermediate protein is also modified by methyl esterification in a nearly stoichiometric manner.

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation., Author Summary Where does a local adaptation take place? In general, an adaptive divergence is predicted to occur between isolated populations because gene flow will erode and prevent the divergence. Therefore, previous genome-wide studies that aim to find the adaptive genes have compared populations that are usually tens of hundreds of kilometers apart. However, because nearby populations are likely to be genetically connected or connected until recently, most of the genome should be undifferentiated, leaving the genetic footprints of natural selections more pronounced. Thus, if an adaptive divergence is to be found within a small spatial scale, such case may favor the screening for the adaptive genes. Here, we took advantage of a unique small-scale local adaptation in Arabidopsis halleri subsp. gemmifera, where similar phenotypic differentiation is found across an altitudinal cline on two distinct mountains. By scanning the genome with a focus on the presence of unidirectional allele frequency shift along the altitudes, we successfully obtained genes with functions that were in line with the known phenotypic and environmental difference between altitudes. Our approach is applicable to any species that show microgeographic divergence and should help understand the genetic basis of small-scale evolution.

Cichlid fishes in Lake Victoria are model organisms for studying rapid radiation and speciation. On the way to examine the molecular basis of how these cichlid fishes achieved such a remarkable morphological diversification, we constructed a bacterial artificial chromosome (BAC) library derived from a cichlid species, Haplochromis chilotes, from Lake Victoria. The library includes 157,056 clones with the average insert size of 128 kb, corresponding to a 10-fold coverage of the H. chilotes genome. Given that the cichlid fishes endemic to Lake Victoria are closely related to one another phylogenetically and their genomes are nearly identical, this BAC library can be utilized to isolate genes from the more than 200 Haplochromine cichlid species in Lake Victoria.

We recently helped to complete the sequence of human chromosome 21 at a very high level of accuracy. Using this sequence we identified two novel genes, designated DSCR9 and DSCR10, in the so-called Down Syndrome Critical Region (DSCR) by computational gene prediction and subsequent cDNA cloning. Both DSCR9 and DSCR10 are expressed preferentially in testis and encode functionally unknown proteins with 149 and 87 amino acid residues, respectively. Zoo blot analysis suggested that both genes are exclusive to primate genomes such as chimpanzee, gorilla, orangutan, crab-eating monkey and African green monkey but are not present in other non-primate mammals including mouse, dog, cat, and chicken. Comparative genomic sequence analysis of DSCR9 and DSCR10 with the corresponding mouse syntenic region confirmed the lack of these genes in the mouse. These results strongly suggest that DSCR9 and DSCR10 have emerged as a new class of gene in the primate lineage during evolution.

The completion of the human genome sequence will greatly accelerate development of a new branch of bioscience and provide fundamental knowledge to biomedical research. We used the sequence information to measure replication timing of the entire lengths of human chromosomes 11q and 21q. Megabase-sized zones that replicate early or late in S phase (thus early/late transition) were defined at the sequence level. Early zones were more GC-rich and gene-rich than were late zones, and early/late transitions occurred primarily at positions identical to or near GC% transitions. We also found the single nucleotide polymorphism (SNP) frequency was high in the late-replicating and replication-transition regions. In the early/late transition regions, concentrated occurrence of cancer-related genes that include CCND1 encoding cyclin D1 (BCL1), FGF4 (KFGF), TIAM1 and FLI1, was observed. The transition regions contained other disease-related genes including APP associated with familial Alzheimer's disease (AD1), SOD1 associated with familial amyotrophic lateral sclerosis (ALS1) and PTS associated with phenylketonuria. These findings are discussed with respect to the prediction that increased DNA damage occurs in replication-transition regions. We propose that genome-wide assessment of replication timing serves as an efficient strategy for identifying disease-related genes.

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration., Author Summary Sex is ubiquitous, but sex-determining mechanisms are evolutionarily labile. Even in taxa with genetic sex determination mechanisms, closely related species often have different sex chromosomes. Although there are many examples of such sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. Using the Japan Sea species of threespine stickleback (Gasterosteus aculeatus), which has a neo-sex chromosome derived from a Y-autosome fusion, we found that incipient neo-X and neo-Y chromosomes have started to diverge in both coding and non-coding regions and to accumulate genes with sex-biased expression. We further found that non-synonymous mutation rates are elevated on the neo-X chromosome in the Japan Sea lineage. Finally, the neo-X chromosome contributes to phenotypic divergence between the Japan Sea species and the ancestral Pacific Ocean species, but the causative genes were distributed broadly on the neo-X chromosome. These results suggest that neo-sex chromosome formation may create a genomic location where mutations causing inter-species differences can accumulate, thereby promoting phenotypic divergence between incipient species.

Using DNA from sorted human chromosomes and two-dimensional gel electrophoresis, we assigned 2295 NotI sites, 43% of the total, to specific chromosomes and designated the procedure CA-RLGS (chromosome-assigned restriction landmark genomic scanning). Although the NotI enzyme is sensitive to DNA methylation, our results suggested that the majority of the spots did not seem to be affected by this modification. The NotI sites were distributed at higher levels in chromosomes 17, 19, and 22, suggesting higher gene content in these chromosomes. Most spots were assigned to unique chromosomes, but some spots were found on two or more chromosomes. Quantitative analysis revealed the intensity of the DNA spots on the sex chromosomes to be haploid and that of the chromosome 21 spots in DNA from a male with Down syndrome to be trisomic, although there were exceptions. We report here the first-generation CA-RLGS map of the human genome.

Background Moyamoya disease is an idiopathic vascular disorder of intracranial arteries. Its susceptibility locus has been mapped to 17q25.3 in Japanese families, but the susceptibility gene is unknown. Methodology/Principal Findings Genome-wide linkage analysis in eight three-generation families with moyamoya disease revealed linkage to 17q25.3 (P, もやもや病感受性遺伝子の特定とその機能についての発見. 京都大学プレスリリース. 2011-7-21.

RAS2 protein of Saccharomyces cerevisiae undergoes post-translational modifications involving methyl esterification and palmitic acid addition, resulting in their association with the plasma membrane. In this paper, we provide evidence that two kinds of proteolytic events accompany the biosynthesis. This is shown by separating and characterizing three intracellular forms of RAS2 protein: precursor, intermediate, and mature (fatty acid-acylated) forms. N-Terminal sequencing has revealed that all three forms start with proline, which is the second amino acid expected from the RAS2 gene sequence. Thus, the first methionine is removed very early during the biosynthesis. Isolation and sequencing of C-terminal peptides indicate that three C-terminal amino acids present in the precursor form are removed in the intermediate and in the fatty acid acylated forms. C-Terminal proteolysis appears to accompany methyl esterification, since the methylation occurs with the intermediate and the fatty acid-acylated forms, but not with the precursor. Palmitic acid is identified as the major fatty acid attached to the fatty acid-acylated form.

A transforming DNA, named lca (for liver cancer), was obtained from a primary human hepatocellular carcinoma (HCC) in transformation assays using NIH 3T3 cells and a calcium phosphate coprecipitation method. High molecular weight DNA obtained from the HCC tissue was employed for this purpose. This transforming DNA had a linkage to the Alu sequence and was cloned in lambda phage for further studies. Restriction enzyme analyses showed that the minimal size of the lca transforming DNA is about 10 kilobase pairs and that its cleavage profiles are different from those of any one of the previously reported human transforming genes or retroviral oncogenes. No cross-hybridization was observed between these genes and the lca DNA. Southern blot analyses of DNAs from flow-sorted human chromosomes and human-mouse somatic cell hybrids indicated that the lca DNA is located on human chromosome 2. An independently obtained transforming DNA from another HCC exhibited identical restriction enzyme cleavage profiles. Thus, lca DNA is likely to represent a commonly encountered transforming DNA in HCC.

Entire genomes of hepatitis B virus (subtype adr) have been cloned. The nucleotide sequence data were compared with other sequences of HBV genome including: adw [Valenzuela et al. (1981) in Animal Virus Genetics. Fields et al. eds. Academic Press, Inc., NY. pp. 57-70], ayw [Galibert et al. (1979) Nature, 281, 646-650], and adyw [Pasek et al. (1979) Nature 282, 575-579]. Four open coding frames for polypeptides larger than 6,000 dalton were found to be conserved and were highly compressed by overlapping with each other in one strand (L-strand). Sites of initiation of the S gene and termination of the P gene were not conserved. No conserved coding frame was found on the opposite strand (S strand). Amino acid sequences of six surface antigen (HBsAg) peptides, including subtypes adr, adw, and ayw, are deduced from the DNA sequences, and the substitution of amino acid residues which are consistent with the change of subtypes are demonstrated.

Products of ras genes are synthesized as precursors in the cytosol and transported to the plasma membrane by a process which involves posttranslational modification by fatty acid. In this paper, we present evidence for the occurrence in the cytosol of an intermediate modification of ras proteins prior to the fatty acid acylation. The modification is detected by a slight shift in the mobility of the protein on SDS polyacrylamide gel. The fatty acid acylation does not contribute to this mobility shift. This modification is affected by the dpr1 mutation which has recently been shown to affect the processing of yeast RAS proteins. To further characterize the nature of the modification event, we have cloned DPR1 gene from the DNA of Saccharomyces cerevisiae. The gene is actively transcribed in yeast cells producing mRNA of approximately 1.6 kb. Genes related to the DRP1 appear to be present in a distantly related yeast, Schizosaccharomyces pombe as well as in guinea pig and human cells.

The initiation sites for DNA synthesis, that mark the transition sites from primer RNA to DNA synthesis, have been located precisely for the discontinuous DNA replication of bacteriophage T7 and for the replication at the origin of E. coli chromosome. On a 339-nucleotide segment of the H-strand of T7 genome (H3 fragment of a Hpa II digest), five initiation sites of DNA synthesis were located. All shared the common sequence of 5′ GAC A C NN′N″ 3′, in which DNA synthesis initiated at N′ or N″. The complementary sequence, 3′ CTG T G N(N′) 5′ in the template strand is very likely the signal for the synthesis of primer RNA which we have shown to have the sequence of pppApCp A p N(pN′). Since the observed common sequence would be present once 100 to 200 nucleotides assuming a random base sequence of T7 DNA, it is used only once in ten rounds of replication so that a 1000 to 2000 nucleotides-long nascent fragment results. On the 245-nucleotide E. coli ori DNA fragment, two transition sites to DNA synthesis were detected but both were on the same strand; the strand which runs from 5′ to 3′ polarity counterclockwise on the E. coli genetic map. The transition sites were separated by 37 nucleotides. No common sequence was found adjacent to the two transition sites but two kinds of common sequences were upstream to each. The results suggest that DNA replication starts uni-directionally within the origin and becomes bidirectional at outside the origin.

An autonomously replicating 2,248-base-pair DNA segment of the mini-F plasmid carries nine 19-base-pair repeating sequences. Five of the repeats are arranged in one direction and form the right cluster, whereas the remaining four repeats are arranged in the opposite direction and form the left cluster (Murotsu et al., Gene 15: 257-271, 1981). Each cluster, cloned separately into the multicopy plasmid vector pBR322, exhibited a strong F-specific incompatibility phenotype (FIP). These clusters were thought to be responsible for the expression of IncB and IncC phenotypes, causing a switchoff function on mini-F replication. Mini-F DNA fragments containing two, three, or more than four repeats were inserted into pBR322. Cells carrying these recombinant plasmids exhibited, respectively, no, intermediate, and strong FIP intensity. Cloning of five repeats into pSC101, whose copy number is about 6 in contrast to 20 for pBR322, resulted in an FIP of intermediate intensity. Thus, the intensity of FIP reflects the dosage of repeats in a cell. The five repeats in the right cluster were eliminated from the mini-F derivative without impairing its autonomous-replicating ability (Bergquist et al., J. Bacteriol. 147: 888-889, 1981; Kline and Palchavdhuri, Plasmid 4: 281-289). Such deletion, however, caused a sixfold elevation of the copy number. When the eliminated cluster of repeats was reinserted in the derivative, the copy number was lowered to the original value, viz., 1 to 2. The position and orientation of this insertion was not important in the copy number control. Thus, the repeats are also related to copy number control. A model to account for the role of the repeating sequences in the control of copy number and FIP is discussed.