Genome assembly of Genji firefly (Nipponoluciola cruciata) reveals novel luciferase-like luminescent proteins without peroxisome targeting signal.

The Genji firefly, Nipponoluciola cruciata, is an aquatic firefly endemic to Japan, inhabiting a wide area of the Japanese archipelago. The luminescence of fireflies is a scientifically interesting phenomenon, and many studies have evaluated this species in Japan. In this study, we sequenced the whole genome of male N. cruciata and constructed a high-quality genome assembly of 662 Mb with a BUSCO completeness of 99.1% in the genome mode. Using the detected set of 15,169 protein-coding genes, the genomic structures and genetic background of luminescence-related genes were also investigated. We found four new firefly luciferase-like genes in the genome. The highest bioluminescent activity was observed for LLa2, which originated from ancestral PDGY, a mitochondrial acyl-CoA synthetase. A thioesterase candidate, NcruACOT1, which is involved in d-luciferin biosynthesis, was expressed in the lantern. Two opsins were also detected and the absorption wavelength of the UV-type opsin candidate shifted from UV to blue. These findings provide an important resource for unravelling the adaptive evolution of fireflies in terms of luminescence and vision.

Arms race between anti-silencing and RdDM in noncoding regions of transposable elements.

Transposable elements (TEs) are among the most dynamic parts of genomes. Since TEs are potentially deleterious, eukaryotes silence them through epigenetic mechanisms such as repressive histone modifications and DNA methylation. We previously reported that Arabidopsis TEs, called VANDALs, counteract epigenetic silencing through a group of sequence-specific anti-silencing proteins, VANCs. VANC proteins bind to noncoding regions of specific VANDAL copies and induce loss of silent chromatin marks. The VANC-target regions form tandem repeats, which diverge rapidly. Sequence-specific anti-silencing allows these TEs to proliferate with minimum host damage. Here, we show that RNA-directed DNA methylation (RdDM) efficiently targets noncoding regions of VANDAL TEs to silence them de novo. Thus, escape from RdDM could be a primary event leading to the rapid evolution and diversification of sequence-specific anti-silencing systems. We propose that this selfish behavior of TEs paradoxically could make them diverse and less harmful to the host.

A divergent RWP-RK transcription factor determines mating type in heterothallic Closterium.

The Closterium peracerosum-strigosum-littorale complex (Closterium, Zygnematophyceae) has an isogamous mating system. Members of the Zygnematophyceae are the closest relatives to extant land plants and are distantly related to chlorophytic models, for which a genetic basis of mating type (MT) determination has been reported. We thus investigated MT determination in Closterium. We sequenced genomes representing the two MTs, mt+ and mt-, in Closterium and identified CpMinus1, a gene linked to the mt- phenotype. We analyzed its function using reverse genetics methods. CpMinus1 encodes a divergent RWP-RK domain-containing-like transcription factor and is specifically expressed during gamete differentiation. Introduction of CpMinus1 into an mt+ strain was sufficient to convert it to a phenotypically mt- strain, while CpMinus1-knockout mt- strains were phenotypically mt+. We propose that CpMinus1 is the major MT determinant that acts by evoking the mt- phenotype and suppressing the mt+ phenotype in heterothallic Closterium. CpMinus1 likely evolved independently in the Zygnematophyceae lineage, which lost an egg-sperm anisogamous mating system. mt- specific regions possibly constitute an MT locus flanked by common sequences that undergo some recombination.

MoG+: a database of genomic variations across three mouse subspecies for biomedical research.

Laboratory mouse strains have mosaic genomes derived from at least three major subspecies that are distributed in Eurasia. Here, we describe genomic variations in ten inbred strains: Mus musculus musculus-derived BLG2/Ms, NJL/Ms, CHD/Ms, SWN/Ms, and KJR/Ms; M. m. domesticus-derived PGN2/Ms and BFM/Ms; M. m. castaneus-derived HMI/Ms; and JF1/Ms and MSM/Ms, which were derived from a hybrid between M. m. musculus and M. m. castaneus. These strains were established by Prof. Moriwaki in the 1980s and are collectively named the "Mishima Battery". These strains show large phenotypic variations in body size and in many physiological traits. We resequenced the genomes of the Mishima Battery strains and performed a comparative genomic analysis with dbSNP data. More than 81 million nucleotide coordinates were identified as variant sites due to the large genetic distances among the mouse subspecies; 8,062,070 new SNP sites were detected in this study, and these may underlie the large phenotypic diversity observed in the Mishima Battery. The new information was collected in a reconstructed genome database, termed MoG+ that includes new application software and viewers. MoG+ intuitively visualizes nucleotide variants in genes and intergenic regions, and amino acid substitutions across the three mouse subspecies. We report statistical data from the resequencing and comparative genomic analyses and newly collected phenotype data of the Mishima Battery, and provide a brief description of the functions of MoG+, which provides a searchable and unique data resource of the numerous genomic variations across the three mouse subspecies. The data in MoG+ will be invaluable for research into phenotype-genotype links in diverse mouse strains.

OryzaGenome2.1: Database of Diverse Genotypes in Wild Oryza Species.

BACKGROUND: OryzaGenome ( http://viewer.shigen.info/oryzagenome21detail/index.xhtml ), a feature within Oryzabase ( https://shigen.nig.ac.jp/rice/oryzabase/ ), is a genomic database for wild Oryza species that provides comparative and evolutionary genomics approaches for the rice research community. RESULTS: Here we release OryzaGenome2.1, the first major update of OryzaGenome. The main feature in this version is the inclusion of newly sequenced genotypes and their meta-information, giving a total of 217 accessions of 19 wild Oryza species (O. rufipogon, O. barthii, O. longistaminata, O. meridionalis, O. glumaepatula, O. punctata, O. minuta, O. officinalis, O. rhizomatis, O. eichingeri, O. latifolia, O. alta, O. grandiglumis, O. australiensis, O. brachyantha, O. granulata, O. meyeriana, O. ridleyi, and O. longiglumis). These 19 wild species belong to 9 genome types (AA, BB, CC, BBCC, CCDD, EE, FF, GG, and HHJJ), representing wide genomic diversity in the genus. Using the genotype information, we analyzed the genome diversity of Oryza species. Other features of OryzaGenome facilitate the use of information on single nucleotide polymorphisms (SNPs) between O. sativa and its wild progenitor O. rufipogon in rice research, including breeding as well as basic science. For example, we provide Variant Call Format (VCF) files for genome-wide SNPs of 33 O. rufipogon accessions against the O. sativa reference genome, IRGSP1.0. In addition, we provide a new SNP Effect Table function, allowing users to identify SNPs or small insertion/deletion polymorphisms in the 33 O. rufipogon accessions and to search for the effect of these polymorphisms on protein function if they reside in the coding region (e.g., are missense or nonsense mutations). Furthermore, the SNP Viewer for 446 O. rufipogon accessions was updated by implementing new tracks for possible selective sweep regions and highly mutated regions that were potentially exposed to selective pressures during the process of domestication. CONCLUSION: OryzaGenome2.1 focuses on comparative genomic analysis of diverse wild Oryza accessions collected around the world and on the development of resources to speed up the identification of critical trait-related genes, especially from O. rufipogon. It aims to promote the use of genotype information from wild accessions in rice breeding and potential future crop improvements. Diverse genotypes will be a key resource for evolutionary studies in Oryza, including polyploid biology.

RNA interference-independent reprogramming of DNA methylation in Arabidopsis.

DNA methylation is important for silencing transposable elements (TEs) in diverse eukaryotes, including plants. In plant genomes, TEs are silenced by methylation of histone H3 lysine 9 (H3K9) and cytosines in both CG and non-CG contexts. The role of RNA interference (RNAi) in establishing TE-specific silent marks has been extensively studied, but the importance of RNAi-independent pathways remains largely unexplored. Here, we directly investigated transgenerational de novo DNA methylation of TEs after the loss of silent marks. Our analyses uncovered potent and precise RNAi-independent pathways for recovering non-CG methylation and H3K9 methylation in most TE genes (that is, coding regions within TEs). Characterization of a subset of TE genes without the recovery revealed the effects of H3K9 demethylation, replacement of histone H2A variants and their interaction with CG methylation, together with feedback from transcription. These chromatin components are conserved among eukaryotes and may contribute to chromatin reprogramming in a conserved manner.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera.

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.

Evolution and Diversity of the Wild Rice Oryza officinalis Complex, across Continents, Genome Types, and Ploidy Levels.

The Oryza officinalis complex is the largest species group in Oryza, with more than nine species from four continents, and is a tertiary gene pool that can be exploited in breeding programs for the improvement of cultivated rice. Most diploid and tetraploid members of this group have a C genome. Using a new reference C genome for the diploid species O. officinalis, and draft genomes for two other C genome diploid species Oryza eichingeri and Oryza rhizomatis, we examine the influence of transposable elements on genome structure and provide a detailed phylogeny and evolutionary history of the Oryza C genomes. The O. officinalis genome is 1.6 times larger than the A genome of cultivated Oryza sativa, mostly due to proliferation of Gypsy type long-terminal repeat transposable elements, but overall syntenic relationships are maintained with other Oryza genomes (A, B, and F). Draft genome assemblies of the two other C genome diploid species, Oryza eichingeri and Oryza rhizomatis, and short-read resequencing of a series of other C genome species and accessions reveal that after the divergence of the C genome progenitor, there was still a substantial degree of variation within the C genome species through proliferation and loss of both DNA and long-terminal repeat transposable elements. We provide a detailed phylogeny and evolutionary history of the Oryza C genomes and a genomic resource for the exploitation of the Oryza tertiary gene pool.

Seasonal Stability and Dynamics of DNA Methylation in Plants in a Natural Environment.

: DNA methylation has been considered a stable epigenetic mark but may respond to fluctuating environments. However, it is unclear how they behave in natural environments. Here, we analyzed seasonal patterns of genome-wide DNA methylation in a single clone from a natural population of the perennial Arabidopsishalleri. The genome-wide pattern of DNA methylation was primarily stable, and most of the repetitive regions were methylated across the year. Although the proportion was small, we detected seasonally methylated cytosines (SeMCs) in the genome. SeMCs in the CHH context were detected predominantly at repetitive sequences in intergenic regions. In contrast, gene-body CG methylation (gbM) itself was generally stable across seasons, but the levels of gbM were positively associated with seasonal stability of RNA expression of the genes. These results suggest the existence of two distinct aspects of DNA methylation in natural environments: sources of epigenetic variation and epigenetic marks for stable gene expression.

De novo assembly of the goldfish (Carassius auratus) genome and the evolution of genes after whole-genome duplication.

For over a thousand years, the common goldfish (Carassius auratus) was raised throughout Asia for food and as an ornamental pet. As a very close relative of the common carp (Cyprinus carpio), goldfish share the recent genome duplication that occurred approximately 14 million years ago in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies provides an exciting opportunity to link genotype to phenotype and to understand the dynamics of genome evolution and speciation. We generated a high-quality draft sequence and gene annotations of a "Wakin" goldfish using 71X PacBio long reads. The two subgenomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged substantially across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression levels across all tissues, while loss of conserved noncoding elements determined expression variance between different tissues. This assembly provides an important resource for comparative genomics and understanding the causes of goldfish variants.

Platanus-allee is a de novo haplotype assembler enabling a comprehensive access to divergent heterozygous regions.

The ultimate goal for diploid genome determination is to completely decode homologous chromosomes independently, and several phasing programs from consensus sequences have been developed. These methods work well for lowly heterozygous genomes, but the manifold species have high heterozygosity. Additionally, there are highly divergent regions (HDRs), where the haplotype sequences differ considerably. Because HDRs are likely to direct various interesting biological phenomena, many genomic analysis targets fall within these regions. However, they cannot be accessed by existing phasing methods, and we have to adopt costly traditional methods. Here, we develop a de novo haplotype assembler, Platanus-allee ( http://platanus.bio.titech.ac.jp/platanus2 ), which initially constructs each haplotype sequence and then untangles the assembly graphs utilizing sequence links and synteny information. A comprehensive benchmark analysis reveals that Platanus-allee exhibits high recall and precision, particularly for HDRs. Using this approach, previously unknown HDRs are detected in the human genome, which may uncover novel aspects of genome variability.

Note to: Hox gene cluster of the ascidian, Halocynthia roretzi, reveals multiple ancient steps of cluster disintegration during ascidian evolution.

BACKGROUND: In the previous paper published in 2017, we described the structure of Hox gene cluster of the ascidian, Halocynthia roretzi, and discussed the scenario for the disintegration of Hox gene clusters during evolution of ascidians. The description about the Hox gene cluster structure still represents the latest information, hence it has been left unchanged. In contrast, some points in Discussion, the description on the phylogenetic relationships among tunicates and the theoretical scenario for the disintegration of Hox gene cluster during evolution of ascidians, should be changed because the phylogenetic relationships among tunicates have recently been updated. The above mentioned points were made in accordance with the phylogenetic tree for tunicates based on the mitochondrial DNA sequences, which was the latest at the time of publication. In 2018, however, Kocot et al. and Delsuc et al. proposed new phylogenetic trees for tunicates based on a large number of nuclear gene sequences. The trees obtained by the two groups are essentially the same and different from the previous one in the phylogenetic positions of Appendicularia and Thaliacea, which leads to a change in the order of the emergence of ascidians and the Hox gene cluster disintegration during evolution of ascidians or tunicates. RESULTS: We add here a note to update the previous description on the phylogenetic relationships among tunicates and the theoretical scenario, including one Figure, so as to coincide with the new phylogenetic relationships among tunicates based on the nuclear gene sequences. CONCLUSION: The previous summarized conclusion remains unchanged: we suggest that the Hox gene cluster of the ancestral ascidian experienced extensive genome shuffling during the course of evolution to Hr and Ci. Nevertheless, some features are shared in Hox gene components and gene organization on the chromosomes, suggesting that Hox gene cluster disintegration in ascidians involved early events common to all ascidians and later lineage-specific events.

High-quality genome assembly of the silkworm, Bombyx mori.

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).

Comparative analysis of seven types of superoxide dismutases for their ability to respond to oxidative stress in Bombyx mori.

Insects are well adapted to changing environmental conditions. They have unique systems for eliminating reactive oxygen species (ROS). Superoxide dismutase (SOD) is a key enzyme that plays a primary role in removing ROS. Bombyx mori is a lepidopteran insect, whose body size is larger than the model insect Drosophila melanogaster, which enabled us to more easily examine gene expression at the tissue level. We searched B. mori SOD (BmSOD) genes using genome database, and we analyzed their function under different type of oxidative stress. Consequently, we identified four new types of BmSODs in addition to the three types already known. Two of the seven types had a unique domain architecture that has not been discovered previously in the SOD family, and they were expressed in different tissues and developmental stages. Furthermore, these BmSODs responded differently to several kinds of stressors. Our results showed that the seven types of BmSODs are likely to play different roles in B. mori; therefore, B. mori could be used to distinguish the functions of each SOD for resistance to oxidative stress that changes with the environmental conditions.

Effects of androgen and oestrogen on IGF pathways controlling phallus growth.

The development of the mammalian phallus involves hormone-dependent mesenchymal-epithelial signalling mechanisms that contribute to urethral closure and regulation of phallus elongation and growth. In marsupials, most differentiation and growth of the phallus occurs post-natally, making them amenable to direct hormone treatment. Expression of IGFs, FGFs, EFNB2, MAFB, DLX5 and AP-1 mRNAs in the phallus at day 50 post-partum (pp) were altered after treatment of tammar wallaby young from day 20 to 40 pp with androgen, oestrogen or after castration at day 25 pp. However, the most interesting changes occurred in the IGF pathway genes. Androgen treatment upregulated IGF1 in female phalluses and oestrogen treatment upregulated IGF1 in male phalluses, but it was downregulated by castration. IGFBP3 was higher in female phalluses and downregulated by androgen. IGF1 expression was higher in all untreated male than in female phalluses from day 50 to 150 pp, but IGFBP3 had the reverse pattern. At day 90 pp, when urethral closure in males is progressing and male phallus growth is accelerating. IGF1 and PCNA protein were only detected in the male urorectal septum, suggesting for the first time that closure and elongation may involve IGF1 activation of cell proliferation specifically in male phalluses. These effects of sex steroids on gene expression and on the IGF1 signalling pathway in particular, suggest that the developing phallus may be especially susceptible to perturbation by exogenous hormones.

Functional divergence of a heterochromatin-binding protein during stickleback speciation.

Intragenomic conflict, the conflict of interest between different genomic regions within an individual, is proposed as a mechanism driving both the rapid evolution of heterochromatin-related proteins and the establishment of intrinsic genomic incompatibility between species. Although molecular studies of laboratory model organisms have demonstrated the link between heterochromatin evolution and hybrid abnormalities, we know little about their link in natural systems. Previously, we showed that F1 hybrids between the Japan Sea stickleback and the Pacific Ocean stickleback show hybrid male sterility and found a region responsible for hybrid male sterility on the X chromosome, but did not identify any candidate genes. In this study, we first screened for genes rapidly evolving under positive selection during the speciation of Japanese sticklebacks to find genes possibly involved in intragenomic conflict. We found that the region responsible for hybrid male sterility contains a rapidly evolving gene encoding a heterochromatin-binding protein TRIM24B. We conducted biochemical experiments and showed that the binding affinity of TRIM24B to a heterochromatin mark found at centromeres and transposons, histone H4 lysine 20 trimethylation (H4K20me3), is reduced in the Japan Sea stickleback. In addition, mRNA expression levels of Trim24b were different between the Japan Sea and the Pacific Ocean testes. Further expression analysis of genes possibly in the TRIM24B-regulated pathway showed that some gypsy retrotransposons are overexpressed in the F1 hybrid testes. We, therefore, demonstrate that a heterochromatin-binding protein can evolve rapidly under positive selection and functionally diverge during stickleback speciation.

Assembling the genome of the African wild rice Oryza longistaminata by exploiting synteny in closely related Oryza species.

The African wild rice species Oryza longistaminata has several beneficial traits compared to cultivated rice species, such as resistance to biotic stresses, clonal propagation via rhizomes, and increased biomass production. To facilitate breeding efforts and functional genomics studies, we de-novo assembled a high-quality, haploid-phased genome. Here, we present our assembly, with a total length of 351 Mb, of which 92.2% was anchored onto 12 chromosomes. We detected 34,389 genes and 38.1% of the genome consisted of repetitive content. We validated our assembly by a comparative linkage analysis and by examining well-characterized gene families. This genome assembly will be a useful resource to exploit beneficial alleles found in O. longistaminata. Our results also show that it is possible to generate a high-quality, functionally complete rice genome assembly from moderate SMRT read coverage by exploiting synteny in a closely related Oryza species.

Anisogamy evolved with a reduced sex-determining region in volvocine green algae.

Male and female gametes differing in size-anisogamy-emerged independently from isogamous ancestors in various eukaryotic lineages, although genetic bases of this emergence are still unknown. Volvocine green algae are a model lineage for investigating the transition from isogamy to anisogamy. Here we focus on two closely related volvocine genera that bracket this transition-isogamous Yamagishiella and anisogamous Eudorina. We generated de novo nuclear genome assemblies of both sexes of Yamagishiella and Eudorina to identify the dimorphic sex-determining chromosomal region or mating-type locus (MT) from each. In contrast to the large (>1 Mb) and complex MT of oogamous Volvox, Yamagishiella and Eudorina MT are smaller (7-268 kb) and simpler with only two sex-limited genes-the minus/male-limited MID and the plus/female-limited FUS1. No prominently dimorphic gametologs were identified in either species. Thus, the first step to anisogamy in volvocine algae presumably occurred without an increase in MT size and complexity.

Draft genome of Dugesia japonica provides insights into conserved regulatory elements of the brain restriction gene nou-darake in planarians.

BACKGROUND: Planarians are non-parasitic Platyhelminthes (flatworms) famous for their regeneration ability and for having a well-organized brain. Dugesia japonica is a typical planarian species that is widely distributed in the East Asia. Extensive cellular and molecular experimental methods have been developed to identify the functions of thousands of genes in this species, making this planarian a good experimental model for regeneration biology and neurobiology. However, no genome-level information is available for D. japonica, and few gene regulatory networks have been identified thus far. RESULTS: To obtain whole-genome information on this species and to study its gene regulatory networks, we extracted genomic DNA from 200 planarians derived from a laboratory-bred asexual clonal strain, and sequenced 476 Gb of data by second-generation sequencing. Kmer frequency graphing and fosmid sequence analysis indicated a complex genome that would be difficult to assemble using second-generation sequencing short reads. To address this challenge, we developed a new assembly strategy and improved the de novo genome assembly, producing a 1.56 Gb genome sequence (DjGenome ver1.0, including 202,925 scaffolds and N50 length 27,741 bp) that covers 99.4% of all 19,543 genes in the assembled transcriptome, although the genome is fragmented as 80% of the genome consists of repeated sequences (genomic frequency ≥ 2). By genome comparison between two planarian genera, we identified conserved non-coding elements (CNEs), which are indicative of gene regulatory elements. Transgenic experiments using Xenopus laevis indicated that one of the CNEs in the Djndk gene may be a regulatory element, suggesting that the regulation of the ndk gene and the brain formation mechanism may be conserved between vertebrates and invertebrates. CONCLUSION: This draft genome and CNE analysis will contribute to resolving gene regulatory networks in planarians. The genome database is available at: http://www.planarian.jp.