A Simple Clinical Scoring System to Determine the Risk of Pancreatic Cancer in the General Population.

This study aimed to develop and validate a simple scoring system to determine the high-risk group for pancreatic cancer (PC) in the asymptomatic general population. The scoring system was developed using data from PC cases and randomly selected non-PC cases undergoing annual medical checkups between 2008 and 2013. The performance of this score was validated for participants with medical checkups between 2014 and 2016. In the development set, 45 PC cases were diagnosed and 450 non-PC cases were identified. Multivariate analysis showed three changes in clinical data from 1 year before diagnosis as independent risk factors: ΔHbA1c ≥ 0.3%, ΔBMI ≤ -0.5, and ΔLDL ≤ -20 mg/dL. A simple scoring system, incorporating variables and abdominal ultrasound findings, was developed. In the validation set, 36 PC cases were diagnosed over a 3-year period from 32,877 participants. The AUROC curve of the scoring system was 0.925 (95%CI 0.877-0.973). The positive score of early-stage PC cases, including Stage 0 and I cases, was significantly higher than that of non-PC cases (80% vs. 6%, p = 0.001). The simple scoring system effectively narrows down high-risk PC cases in the general population and provides a reasonable approach for early detection of PC.

A Novel Method of Calculating the Drained Liver Volume Using a 3D Volume Analyzer for Biliary Drainage of Unresectable Malignant Hilar Biliary Obstruction.

OBJECTIVES: In patients with unresectable malignant hilar biliary obstruction (UMHBO), drainage of ≥ 50% liver volume correlates with better clinical outcomes. Accurately measuring the liver volume to be drained by biliary stents is required. We aimed to develop a novel method for calculating the drained liver volume (DLV) using a 3D volume analyzer (3D volumetry), and assess the usefulness for drainage in patients with UMHBO. METHODS: Three-dimensional volumetry comprises the following steps: (1) manual tracing of bile duct using 3D imaging system; (2) 3D reconstruction of bile duct and liver parenchyma; and (3) calculating DLV according to the 3D distribution of bile ducts. Using 3D volumetry, we reviewed data of patients who underwent biliary drainage for UMHBO, calculated the DLV, and determined the association between DLV and biliary drainage outcome. RESULTS: There were 104 eligible cases. The mean DLV was 708 ± 393 ml (53% ± 21%). and 65 patients (63%) underwent drainage of ≥50% liver volume. The clinical success rate was significantly higher in patients with DLV ≥ 50% than in patients with DLV < 50% (89% vs. 28%, P < 0.001). The median time to recurrence of biliary obstruction (TRBO) and survival time were significantly longer in patients with DLV ≥ 50% than in patients with DLV < 50% (TRBO, 292 vs. 119 days, P = 0.03; survival, 285 vs. 65days, P = 0.004, log-rank test, respectively). CONCLUSIONS: Three-dimensional volumetry, a novel method to calculate DLV accurately according to bile duct distribution was useful for drainage in UMHBO patients.

Complemental Diagnosis of IgG4-Related Pancreaticobiliary Diseases by Multiple Hypoechoic Lesions in the Submandibular Glands.

The diagnosis of autoimmune pancreatitis (AIP) and immunoglobulin G4-related sclerosing cholangitis (IgG4-SC) may require a somewhat invasive pathological examination and steroid responsiveness. This retrospective study assessed the complemental diagnosis of AIP and IgG4-SC using submandibular gland (SG) ultrasonography (US) in 69 patients, including 54 patients with AIP, 2 patients with IgG4-SC, and 13 patients with both AIP and IgG4-SC. The data from the physical examination and US of SGs to diagnose AIP (n = 67) and IgG4-SC (n = 15) were analyzed. The steroid therapy efficacy in resolving hypoechoic lesions in SGs was evaluated in 36 cases. The presence of IgG4-related pancreaticobiliary disease with multiple hypoechoic lesions in SGs was reduced from 31 to 11 cases after steroid therapy, suggesting that multiple hypoechoic lesions in SGs are strongly associated with IgG4-positive cell infiltrations. Multiple hypoechoic lesions in SGs were observed in 53 cases, whereas submandibular swelling on palpation was observed in 21 cases of IgG4-related pancreaticobiliary diseases. A complemental diagnosis of IgG4-related pancreaticobiliary diseases without a histological diagnosis and steroid therapy was achieved in 57 and 68 cases without and with multiple hypoechoic lesions in SGs, respectively. In conclusion, multiple hypoechoic lesions in SGs are useful for the complemental diagnosis of IgG4-related pancreaticobiliary diseases.

Global population structure of the Serratia marcescens complex and identification of hospital-adapted lineages in the complex.

Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the 'S. marcescens complex'), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila, S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum ß-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features.

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.

A ubiquitous subcuticular bacterial symbiont of a coral predator, the crown-of-thorns starfish, in the Indo-Pacific.

BACKGROUND: Population outbreaks of the crown-of-thorns starfish (Acanthaster planci sensu lato; COTS), a primary predator of reef-building corals in the Indo-Pacific Ocean, are a major threat to coral reefs. While biological and ecological knowledge of COTS has been accumulating since the 1960s, little is known about its associated bacteria. The aim of this study was to provide fundamental information on the dominant COTS-associated bacteria through a multifaceted molecular approach. METHODS: A total of 205 COTS individuals from 17 locations throughout the Indo-Pacific Ocean were examined for the presence of COTS-associated bacteria. We conducted 16S rRNA metabarcoding of COTS to determine the bacterial profiles of different parts of the body and generated a full-length 16S rRNA gene sequence from a single dominant bacterium, which we designated COTS27. We performed phylogenetic analysis to determine the taxonomy, screening of COTS27 across the Indo-Pacific, FISH to visualize it within the COTS tissues, and reconstruction of the bacterial genome from the hologenome sequence data. RESULTS: We discovered that a single bacterium exists at high densities in the subcuticular space in COTS forming a biofilm-like structure between the cuticle and the epidermis. COTS27 belongs to a clade that presumably represents a distinct order (so-called marine spirochetes) in the phylum Spirochaetes and is universally present in COTS throughout the Indo-Pacific Ocean. The reconstructed genome of COTS27 includes some genetic traits that are probably linked to adaptation to marine environments and evolution as an extracellular endosymbiont in subcuticular spaces. CONCLUSIONS: COTS27 can be found in three allopatric COTS species, ranging from the northern Red Sea to the Pacific, implying that the symbiotic relationship arose before the speciation events (approximately 2 million years ago). The universal association of COTS27 with COTS and nearly mono-specific association at least with the Indo-Pacific COTS provides a useful model system for studying symbiont-host interactions in marine invertebrates and may have applications for coral reef conservation. Video Abstract.

Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process.

De novo assembly of short DNA reads remains an essential technology, especially for large-scale projects and high-resolution variant analyses in epidemiology. However, the existing tools often lack sufficient accuracy required to compare closely related strains. To facilitate such studies on bacterial genomes, we developed Platanus_B, a de novo assembler that employs iterations of multiple error-removal algorithms. The benchmarks demonstrated the superior accuracy and high contiguity of Platanus_B, in addition to its ability to enhance the hybrid assembly of both short and nanopore long reads. Although the hybrid strategies for short and long reads were effective in achieving near full-length genomes, we found that short-read-only assemblies generated with Platanus_B were sufficient to obtain ≥90% of exact coding sequences in most cases. In addition, while nanopore long-read-only assemblies lacked fine-scale accuracies, inclusion of short reads was effective in improving the accuracies. Platanus_B can, therefore, be used for comprehensive genomic surveillances of bacterial pathogens and high-resolution phylogenomic analyses of a wide range of bacteria.

Rapidly declining trend of signet ring cell cancer of the stomach may parallel the infection rate of Helicobacter pylori.

BACKGROUND: Studies indicate that gastric cancer (GC) incidence has decreased, whereas signet ring cell carcinoma (SRC) incidence has increased. However, recent trends in GC incidence are unclear. We used our hospital cancer registry to evaluate the changes in the incidence of GC, SRC, and non-SRC (NSRC) over time in comparison to changes in the H. pylori infection rates over time. METHODS: We identified 2532 patients with GC enrolled in our registry between January 2007 and December 2018 and statistically analyzed SRC and NSRC incidence. The H. pylori infection rate in patients with SRC was determined by serum anti-H. pylori antibody testing, urea breath test, biopsy specimen culture, and immunohistochemical analysis (IHC) of gastric tissue. Additionally, genomic detection of H. pylori was performed in SRCs by extracting DNA from formalin-fixed paraffin-embedded gastric tissue and targeting 16S ribosomal RNA of H. pylori. RESULTS: Overall, 211 patients had SRC (8.3%). Compared with patients with NSRC, those with SRC were younger (P <  0.001) and more likely to be female (P <  0.001). Time series analysis using an autoregressive integrated moving average model revealed a significant decrease in SRC (P <  0.001) incidence; NSRC incidence showed no decline. There was no difference in H. pylori infection prevalence between the SRC and NSRC groups. IHC and genomic methods detected H. pylori in 30 of 37 (81.1%) SRCs. CONCLUSIONS: Reduction in H. pylori infection prevalence may be associated with the decrease in the incidence of SRC, which was higher than that of NSRC.

Genome Sequencing Verifies Relapsed Infection of Helicobacter cinaedi.

BACKGROUND: Recurrent infections of Helicobacter cinaedi are often reported, and long-term antimicrobial treatment is empirically recommended to prevent such infections. However, there have been no studies examining whether recurrent infections are relapses of former infections or reinfections with different clones. METHODS: A 69-year-old woman presented with recurrent H cinaedi bacteremia-associated cellulitis after a 51-day interval. We isolated 10 colonies from the blood cultures obtained during each of the 2 episodes and subjected them to whole-genome sequencing (WGS). High-confidence single-nucleotide polymorphisms (SNPs) were identified by an assembly based method. Heterogeneous SNP sites were identified by read mapping. The susceptibility of a representative isolate to 14 antimicrobials was also examined. RESULTS: Whole-genome sequence analysis revealed only 6 SNP sites among the 20 isolates at the whole-genome level. Based on the 6 SNPs, 5 within-host variants (referred to as genotypes) were identified. All 5 genotypes were detected in the first infection; however, only 2 genotypes were detected in the second infection. Although the H cinaedi clone showed a higher minimum inhibitory concentration to fluoroquinolones and macrolides and responsible mutations were identified, none of the 6 SNPs appeared related to additional resistance. CONCLUSIONS: The second infection analyzed here was a relapse of the first infection. A certain level of within-host genomic heterogeneity of the H cinaedi clone was already present in the first infection. Our results suggest the importance of longer treatment courses to eradicate H cinaedi for preventing the relapse of its infection.

Evaluation of SNP calling methods for closely related bacterial isolates and a novel high-accuracy pipeline: BactSNP.

Bacteria are highly diverse, even within a species; thus, there have been many studies which classify a single species into multiple types and analyze the genetic differences between them. Recently, the use of whole-genome sequencing (WGS) has been popular for these analyses, and the identification of single-nucleotide polymorphisms (SNPs) between isolates is the most basic analysis performed following WGS. The performance of SNP-calling methods therefore has a significant effect on the accuracy of downstream analyses, such as phylogenetic tree inference. In particular, when closely related isolates are analyzed, e.g. in outbreak investigations, some SNP callers tend to detect a high number of false-positive SNPs compared with the limited number of true SNPs among isolates. However, the performances of various SNP callers in such a situation have not been validated sufficiently. Here, we show the results of realistic benchmarks of commonly used SNP callers, revealing that some of them exhibit markedly low accuracy when target isolates are closely related. As an alternative, we developed a novel pipeline BactSNP, which utilizes both assembly and mapping information and is capable of highly accurate and sensitive SNP calling in a single step. BactSNP is also able to call SNPs among isolates when the reference genome is a draft one or even when the user does not input the reference genome. BactSNP is available at https://github.com/IEkAdN/BactSNP.

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.

Validation of Tokyo Guideline 2013 as Treatment of Acute Cholecystitis by Real World Data.

BACKGROUND: The Tokyo Guidelines (TG; 2013) indicated that emergency cholecystectomy is an important early treatment option for acute cholecystitis; however, surgical intervention is not necessarily indicated in patients with advanced age. We evaluated percutaneous transhepatic gallbladder aspiration (PTGBA), percutaneous transhepatic gallbladder drainage (PTGBD), and the administration of antibiotics alone as treatment options for acute -cholecystitis. METHODS: From January 2010 to December 2017, 159 patients with acute cholecystitis were treated at our institution. The data from these patients were retrospectively analyzed. RESULTS: Of these 159 cases, 109 underwent PTGBA, 28 underwent PTGBD, and 22 were administered antibiotics alone. None of the 159 patients needed urgent (early) cholecystectomy, and all patients were discharged without mortality. PTGBA was unsuccessful in only 6 of 109 patients; PTGBD was performed in these 6 cases. Long-term follow-up was conducted in all cases. Of the 159 patients, 146 had gallbladder stones initially, while 13 had none at the time of presentation. Of these 146 patients with gallbladder stones, 84 underwent elective cholecystectomy, while 62 did not. Of the 84 patients who underwent elective cholecystectomy, 2 developed choledocholithiasis; of the 62 patients who did not undergo elective cholecystectomy, 5 developed choledocholithiasis and 2 developed acute cholecystitis. The incidences of choledocholithiasis and acute cholecystitis did not significantly differ between the 2 groups (p = 0.06). CONCLUSIONS: Despite the recommendations in the TG (2013), emergency cholecystectomy was not needed in any of the present patients with acute cholecystitis. Acute cholecystitis can be successfully treated with -PTGBA or PTGBD, which are simple procedures with good short- and long-term safety. These procedures are highly recommended for patients with acute cholecystitis, especially in the elderly population.

Deficiency of mismatch repair genes is less frequently observed in signet ring cell compared with non-signet ring cell gastric cancer.

Signet ring cell (SRC) gastric cancer at advanced stage has poor prognosis. While a recent study reported nearly one-third of SRC cases contain tumors with deficient mismatch repair (MMR) genes, other studies in SRC have been inconclusive. To re-analyze the results, we performed immunohistochemical staining of MLH1, MSH2, MSH6 and PMS2 proteins in 38 SRC gastric tumors compared with 109 non-SRC (NSRC) tumors from 94 patients. In contrast to the previous study, all SRC gastric tumors normally expressed MMR proteins, whereas 22 of 109 of NSRC (20%) showed deficient MMR proteins. To reinforce our results, we referred to the Cancer Genome Atlas (TCGA) genomic database and found that only 6 (6%) of 99 samples with diffuse gastric tumors showed deficient MMR, whereas 64 (21%) of 304 in intestinal gastric tumors showed deficient MMR. Our results as well as the TCGA database indicated that MMR genes are infrequently inactivated in SRC gastric cancer. These findings indicate that SRC patients may not be the best candidates for immuno-oncology therapy.

Molecular subtype switching in early-stage gastric cancers with multiple occurrences.

BACKGROUND: Multiple gastric cancers at the same time (synchronous) or recurrence after 1 year (metachronous) are frequently encountered. Since their genetic profiles were not well elucidated, we molecularly subtyped the genetic events of synchronous and metachronous early-stage gastric cancers. METHODS: We studied mismatch repair (MMR) genes in 84 tumors from 31 patients (15 synchronous and 16 metachronous) by immunohistochemistry. We performed microsatellite instability analysis and targeted sequencing of 58 significantly mutated genes (SMGs) in 35 tumors from thirteen patients. Genomic data from TCGA were used for comparisons with advanced-stage cancers. RESULTS: Among the 31 patients, at least one deficient-MMR (dMMR) tumor was observed in eight (26%). Of eight patients, seven showed a mixture of proficient-MMR (pMMR) and dMMR tumors. The one case with only dMMR had six recurrent tumors within 2 years. To further subtype, we sequenced 58 SMGs in 35 samples (25 pMMR and 10 dMMR) from thirteen patients. In 35 samples, 163 mutations were identified, but none matched in almost cases, strongly indicating different clonal origins, whether synchronous or metachronous occurrences. Of the 25 pMMR cases, 1 belonged to Epstein-Barr virus (EBV), 24 belonged to chromosomal instability (CIN) subtypes. Of the thirteen cases, repetitive CIN, a mixture of CIN and MSI, a mixture of CIN and EBV, and repetitive MSI were observed in nine (70%), two (15%), one (8%) and one (8%), respectively. CONCLUSIONS: Despite multiple tumors occurring in the same patient simultaneously or several years apart, clonal origin was totally different. 'Switching' or 'mixing' of dMMR and pMMR, EBV or CIN occurred, which had clinical relevance with regard to immunotherapy.

Influence of NUDT15 variants on hematological pictures of patients with inflammatory bowel disease treated with thiopurines.

AIM: The single nucleotide polymorphism (SNP) c.415C>T in exon 3 of NUDT15 affects thiopurine-induced leukopenia in Asian patients with Crohn's disease. Meanwhile, three additional genetic variants of NUDT15 were reported in patients with acute lymphoblastic leukemia. We evaluated the effects of these additional genetic variants of NUDT15 in patients with inflammatory bowel disease (IBD) treated with thiopurines. METHODS: Ninety-six Japanese patients with IBD were enrolled. Genotyping for the NUDT15 and TPMT genes was performed using Custom TaqMan SNP genotyping assays or Sanger sequencing. The changes in white blood cell (WBC) count, mean corpuscular volume (MCV), platelet count, hemoglobin, CRP, amylase, albumin, AST, ALT, and ESR were evaluated. RESULTS: Genetic variants of exon 1 and exon 3 of NUDT15 were identified in 24 of 96 patients (25.0%). C.52G > A and c.36_37insGGAGTC in exon 1 were found in three patients each. All three patients with c.36_37insGGAGTC in exon 1 were heterozygotes of p.Arg139Cys in exon 3. Eighteen patients had p.Arg139Cys in exon 3 alone. The WBC count gradually decreased after initiation of thiopurine treatment in the mutated cases (n = 24), and was significantly lower at 6, 8, 10, and 16 wk (P = 0.0271, 0.0037, 0.0051, and 0.0185, respectively). The WBC counts were also evaluated in patients with and without prednisolone treatment. In the patients with prednisolone treatment, the WBC count tended to show a greater decrease in the mutated cases, with significant differences at 8 and 10 wk (P = 0.012 and 0.029, respectively). In the patients without prednisolone treatment, the WBC count was significantly lower at 2, 4, 8, and 14 wk in mutated cases (P = 0.0196, 0.0182, 0.0237 and 0.0241, respectively). MCV increased after starting thiopurine treatment in the mutated cases, and was significantly higher at 10 wk (P = 0.0085). Platelet count, hemoglobin, CRP, amylase, albumin, AST, ALT and ESR did not differ significantly between the wild-type and mutated cases. TPMT mutations were not found in any of the patients. CONCLUSION: Mutations in exon 1 of NUDT15 also affect thiopurine-induced leukopenia in patients with IBD. To discuss thiopurine-induced leukopenia in more detail, investigation of SNPs in both exon 1 and exon 3 of NUDT15 is needed.

Multi-step genomic dissection of a suspected intra-hospital Helicobacter cinaedi outbreak.

Helicobacter cinaedi is an emerging pathogen causing bacteraemia and cellulitis. Nosocomial transmission of this microbe has been described, but detailed molecular-epidemiological analyses have not been performed. Here, we describe the results of a multi-step genome-wide phylogenetic analysis of a suspected intra-hospital outbreak of H. cinaedi that occurred in a hospital in Japan. The outbreak was recognized by the infectious control team (ICT) of the hospital as a sudden increase in H. cinaedi bacteraemia. ICT defined this outbreak case based on 16S rRNA sequence data and epidemiological information, but were unable to determine the source and route of the infections. We therefore re-investigated this case using whole-genome sequencing (WGS). We first performed a species-wide analysis using publicly available genome sequences to understand the level of genomic diversity of this under-studied species. The clusters identified were then separately analysed using the genome sequence of a representative strain in each cluster as a reference. These analyses provided a high-level phylogenetic resolution of each cluster, identified a confident set of outbreak isolates, and discriminated them from other closely related but distinct clones, which were locally circulating and invaded the hospital during the same period. By considering the epidemiological data, possible strain transmission chains were inferred, which highlighted the role of asymptomatic carriers or environmental contamination. The emergence of a subclone with increased resistance to fluoroquinolones in the outbreak was also recognized. Our results demonstrate the impact of the use of a closely related genome as a reference to maximize the power of WGS.

Extremely Low Genomic Diversity of Rickettsia japonica Distributed in Japan.

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.