Rickettsiae are obligate intracellular bacteria that have small genomes as a result of reductive evolution. Many Rickettsia species of the spotted fever group (SFG) cause tick-borne diseases known as "spotted fevers". The life cycle of SFG rickettsiae is closely associated with that of the tick, which is generally thought to act as a bacterial vector and reservoir that maintains the bacterium through transstadial and transovarial transmission. Each SFG member is thought to have adapted to a specific tick species, thus restricting the bacterial distribution to a relatively limited geographic region. These unique features of SFG rickettsiae allow investigation of how the genomes of such biologically and ecologically specialized bacteria evolve after genome reduction and the types of population structures that are generated. Here, we performed a nationwide, high-resolution phylogenetic analysis of Rickettsia japonica, an etiological agent of Japanese spotted fever that is distributed in Japan and Korea. The comparison of complete or nearly complete sequences obtained from 31 R. japonica strains isolated from various sources in Japan over the past 30 years demonstrated an extremely low level of genomic diversity. In particular, only 34 single nucleotide polymorphisms were identified among the 27 strains of the major lineage containing all clinical isolates and tick isolates from the three tick species. Our data provide novel insights into the biology and genome evolution of R. japonica, including the possibilities of recent clonal expansion and a long generation time in nature due to the long dormant phase associated with tick life cycles.

Recent studies on paleogenomics have reported some Paleolithic and Neolithic genomes that have provided new insights into the human population history in East and Northeast Asia. However, there remain some cases where more recent migration events need to be examined to elucidate the detailed formation process of local populations. Although the area around northern Japan is one of the regions archaeologically suggested to have been affected by migration waves after the Neolithic period, the genetic source of these migrations are still unclear. Thus, genomic data from such past migrant populations would be highly informative to clarify the detailed formation process of local populations in this region. Here, we report the genome sequence of a 900-year-old adult female (NAT002) belonging to the prehistoric Okhotsk people, who have been considered to be the past migrants to northern Japan after the Neolithic period. We found a close relationship between NAT002 and modern Lower Amur populations and past admixture events between the Amur, Jomon, and Kamchatka ancestries. The admixture dating suggested migration of Amur-related ancestry at approximately 1,600 BP, which is compatible with the archaeological evidence regarding the settlement of the Okhotsk people. Our results also imply migration of Kamchatka-related ancestry at approximately 2,000 BP. In addition, human leukocyte antigen (HLA) typing detected the HLA-B*40 allele, which is reported to increase the risk of arthritis, suggesting the genetic vulnerability of NAT002 to hyperostosis, which was observed around her chest clavicle.

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.

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than 2-fold heritable genome size variation among the unicellular Zygnematophycean alga, Closterium peracerosum-strigosum-littorale (C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1 segregation analysis. Six de novo assembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in the C. psl. complex, despite its possible detrimental effects.

The paper nautilus or greater argonaut, Argonauta argo, is a species of octopods which is characterized by its pelagic lifestyle and by the presence of a protective spiral-shaped shell-like eggcase in females. To reveal the genomic background of how the species adapted to the pelagic lifestyle and acquired its shell-like eggcase, we sequenced the draft genome of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters, and some gene clusters that could probably be related to the pelagic lifestyle, such as reflectin, tyrosinase, and opsin. The gene models also revealed several homologous genes related to calcified shell formation in Conchiferan mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the shell-less octopus, as well as Nautilus, which has a true outer shell. Therefore, the draft genome sequence of Arg. argo presented here has helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes to form an important, converging extended phenotypic structure such as the shell and the shell-like eggcase. Additionally, it allows us to explore the evolution of from benthic to pelagic lifestyles in cephalopods and octopods.

Plastid genomes are not normally celebrated for being large. But researchers are steadily uncovering algal lineages with big and, in rare cases, enormous plastid DNAs (ptDNAs), such as volvocine green algae. Plastome sequencing of five different volvocine species has revealed some of the largest, most repeat-dense plastomes on record, including that of Volvox carteri (∼525 kb). Volvocine algae have also been used as models for testing leading hypotheses on organelle genome evolution (e.g., the mutational hazard hypothesis), and it has been suggested that ptDNA inflation within this group might be a consequence of low mutation rates and/or the transition from a unicellular to multicellular existence. Here, we further our understanding of plastome size variation in the volvocine line by examining the ptDNA sequences of the colonial species Yamagishiella unicocca and Eudorina sp. NIES-3984 and the multicellular Volvox africanus, which are phylogenetically situated between species with known ptDNA sizes. Although V. africanus is closely related and similar in multicellular organization to V. carteri, its ptDNA was much less inflated than that of V. carteri. Synonymous- and noncoding-site nucleotide substitution rate analyses of these two Volvox ptDNAs suggest that there are drastically different plastid mutation rates operating in the coding versus intergenic regions, supporting the idea that error-prone DNA repair in repeat-rich intergenic spacers is contributing to genome expansion. Our results reinforce the idea that the volvocine line harbors extremes in plastome size but ultimately shed doubt on some of the previously proposed hypotheses for ptDNA inflation within the lineage.

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

Semisulcospira habei is a freshwater snail species endemic to the Lake Biwa drainage and belongs to a species group radiated within the lake system. We report the chromosome-scale genome assembly of S. habei, including eight megascaffolds larger than 150 Mb. The genome assembly size is about 2.0 Gb with an N50 of 237 Mb. There are 41,547 protein-coding genes modeled by ab initio gene prediction based on the transcriptome data set, and the BUSCO completeness of the annotated genes was 92.2%. The repeat elements comprise approximately 76% of the genome assembly. The Hi-C contact map showed seven well-resolved scaffolds that correspond to the basic haploid chromosome number of S. habei inferred from the preceding karyotypic study, while it also exhibited one scaffold with a complicated mosaic pattern that is likely to represent the complex of multiple supernumerary chromosomes. The genome assembly reported here represents a high-quality genome resource in disentangling the genomic background of the adaptive radiation of Semisulcospira and also facilitates evolutionary studies in the superfamily Cerithioidea.

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Fish use olfaction to sense a variety of nonvolatile chemical signals in water. However, the evolutionary importance of olfaction in species-rich cichlids is controversial. Here, we determined an almost complete sequence of the vomeronasal type 2 receptor-like (OlfC: putative amino acids receptor in teleosts) gene cluster using the bacterial artificial chromosome library of the Lake Victoria cichlid, Haplochromis chilotes. In the cluster region, we found 61 intact OlfC genes, which is the largest number of OlfC genes identified among the seven teleost fish investigated to date. Data mining of the Oreochromis niloticus (Nile tilapia) draft genome sequence, and genomic Southern hybridization analysis revealed that the ancestor of all modern cichlids had already developed almost the same OlfC gene repertoire, which was accomplished by lineage-specific gene expansions. Furthermore, comparison of receptor sequences showed that recently duplicated paralogs are more variable than orthologs of different species at particular sites that were predicted to be involved in amino acid selectivity. Thus, the increase of paralogs through gene expansion may lead to functional diversification in detection of amino acids. This study implies that cichlids have developed a potent capacity to detect a variety of amino acids (and their derivatives) through OlfCs, which may have contributed to the extraordinary diversity of their feeding habitats.