The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution.

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Rapid seasonal changes in phenotypes in a wild Drosophila population.

Seasonal environmental change is one of the most rapid and striking environmental variables. Although relatively rapid adaptation to environmental changes over several years or several decades has been described in many taxa, rapid responses to seasonal environments are delicate, and therefore, the detection of the evolutionary responses requires sensitive methods. In this study, we examined seasonal changes in phenotypes related to thermal tolerance and morphological traits of Drosophila lutescens collected at the spring and autumn periods from a single location. We first demonstrated that flies in the two seasonal periods were almost genetically identical using double-digest restriction site-associated DNA sequencing and analysis. Using an experimental design to eliminate the effect of possible confounding factors that influence phenotypes (i.e., maternal effects and the environmental conditions in which each phenotype was analyzed), we showed that the heat tolerance of D. lutescens was significantly higher in the autumn population than in the spring population. Furthermore, cold tolerance was slightly higher in the spring population than in the autumn one. Although wing length and thorax length did not change significantly between seasons, the ratio of wing length to thorax length changed significantly between them. These results suggest that seasonal environmental heterogeneity induces rapid phenotypic changes within a year. Finally, we discuss the possibility of rapid evolutionary responses to seasonal changes.

Telomere-to-telomere genome assembly of an allotetraploid pernicious weed, Echinochloa phyllopogon.

Echinochloa phyllopogon is an allotetraploid pernicious weed species found in rice fields worldwide that often exhibit resistance to multiple herbicides. An accurate genome sequence is essential to comprehensively understand the genetic basis underlying the traits of this species. Here, the telomere-to-telomere genome sequence of E. phyllopogon was presented. Eighteen chromosome sequences spanning 1.0 Gb were constructed using the PacBio highly fidelity long technology. Of the 18 chromosomes, 12 sequences were entirely assembled into telomere-to-telomere and gap-free contigs, whereas the remaining six sequences were constructed at the chromosomal level with only eight gaps. The sequences were assigned to the A and B genome with total lengths of 453 and 520 Mb, respectively. Repetitive sequences occupied 42.93% of the A genome and 48.47% of the B genome, although 32,337, and 30,889 high-confidence genes were predicted in the A and B genomes, respectively. This suggested that genome extensions and gene disruptions caused by repeated sequence accumulation often occur in the B genome before polyploidization to establish a tetraploid genome. The highly accurate and comprehensive genome sequence could be a milestone in understanding the molecular mechanisms of the pernicious traits and in developing effective weed control strategies to avoid yield loss in rice production.

From green to red: Urban heat stress drives leaf color evolution.

Prevalence of impervious surface and resulting higher temperatures in urban areas, known as urban heat islands, comprises prominent characteristics in global cities. However, it is not known whether and how urban plants adapt to such heat stress. This study focused on Oxalis corniculata, which has intraspecific polymorphism in leaf color (green and red) and examined whether the leaf color variation is associated with urban heat stress. Field observations revealed that green-leaved plants were dominant in green habitats, and red-leaved individuals were dominant in urban habitats, at local (<500 meters), landscape (<50 kilometers), and global scales. Growth and photosynthesis experiments demonstrated that red-leaved individuals performed better under heat stress, while green-leaved individuals performed better under nonstressful conditions. Genome-wide SNP analysis suggests that the red leaf may have evolved multiple times from the ancestral green leaf. Overall, the results suggest that the red leaves of O. corniculata observed in cities worldwide are evidence of plant adaptive evolution due to urban heat islands.

Discrimination of Camellia cultivars using iD-NA analysis.

Recently, many new cultivars have been taken abroad illegally, which is now considered an international issue. Botanical evidence found at a crime scene provides valuable information about the origin of the sample. However, botanical resources for forensic evidence remain underutilized because molecular markers, such as microsatellites, are not available without a limited set of species. Multiplexed intersimple sequence repeat (ISSR) genotyping by sequencing (MIG-seq) and its analysis method, identification of not applicable (iD-NA), have been used to determine several genome-wide genetic markers, making them applicable to all plant species, including those with limited available genetic information. Camellia cultivars are popular worldwide and are often planted in many gardens and bred to make new cultivars. In this study, we aimed to analyze Camellia cultivars/species through MIG-seq. MIG-seq could discriminate similar samples, such as bud mutants and closely related samples that could not be distinguished based on morphological features. This discrimination was consistent with that of a previous study that classified cultivars based on short tandem repeat (STR) markers, indicating that MIG-seq has the same or higher discrimination ability as STR markers. Furthermore, we observed unknown phylogenetic relationships. Because MIG-seq can be applied to unlimited species and low-quality DNA, it may be useful in various scientific fields.

Potential contribution of floral thermogenesis to cold adaptation, distribution pattern, and population structure of thermogenic and non/slightly thermogenic Symplocarpus species.

The genus Symplocarpus in basal Araceae includes both thermogenic and non/slightly thermogenic species that prefer cold environments. If floral thermogenesis of Symplocarpus contributes to cold adaptation, it would be expected that thermogenic species have a larger habitat than non/slightly thermogenic species during an ice age, leading to increased genetic diversity in the current population. To address this question, potential distribution in past environment predicted by ecological niche modeling (ENM), genetic diversity, and population structure of chloroplast and genome-wide single nucleotide polymorphisms were compared between thermogenic Symplocarpus renifolius and non/slightly thermogenic Symplocarpus nipponicus. ENM revealed that the distribution of S. nipponicus decreased, whereas that of S. renifolius expanded in the Last Glacial Maximum. Phylogeographic analyses have shown that the population structures of the two species were genetically segmented and that the genetic diversity of S. renifolius was higher than that of S. nipponicus. The phylogenetic relationship between chloroplast and nuclear DNA is topologically different in the two species, which may be due to the asymmetric gene flow ubiquitously observed in plants. The results of this study imply that floral thermogenesis of Symplocarpus contributes to expanding the distribution during an ice age, resulting in increased genetic diversity due to cold adaptation.

Molecular characterization of SrSTP14, a sugar transporter from thermogenic skunk cabbage, and its possible role in developing pollen.

In floral thermogenesis, sugars play an important role not only as energy providers but also as growth and development facilitators. Yet, the mechanisms underlying sugar translocation and transport in thermogenic plants remain to be studied. Asian skunk cabbage (Symplocarpus renifolius) is a species that can produce durable and intense heat in its reproductive organ, the spadix. Significant morphological and developmental changes in the stamen are well-characterized in this plant. In this study, we focused on the sugar transporters (STPs), SrSTP1 and SrSTP14, whose genes were identified by RNA-seq as the upregulated STPs during thermogenesis. Real-time PCR confirmed that mRNA expression of both STP genes was increased from the pre-thermogenic to the thermogenic stage in the spadix, where it is predominantly expressed in the stamen. SrSTP1 and SrSTP14 complemented the growth defects of a hexose transporter-deficient yeast strain, EBY4000, on media containing 0.02, 0.2, and 2% (w/v) glucose and galactose. Using a recently developed transient expression system in skunk cabbage leaf protoplasts, we revealed that SrSTP1 and SrSTP14-GFP fusion proteins were mainly localized to the plasma membrane. To dig further into the functional analysis of SrSTPs, tissue-specific localization of SrSTPs was investigated by in situ hybridization. Using probes for SrSTP14, mRNA expression was observed in the microspores within the developing anther at the thermogenic female stage. These results indicate that SrSTP1 and SrSTP14 transport hexoses (e.g., glucose and galactose) at the plasma membrane and suggest that SrSTP14 may play a role in pollen development through the uptake of hexoses into pollen precursor cells.

Telomere-to-telomere genome assembly of matsutake (Tricholoma matsutake).

Here, we report the first telomere-to-telomere genome assembly of matsutake (Tricholoma matsutake), which consists of 13 sequences (spanning 161.0 Mb) and a 76 kb circular mitochondrial genome. All the 13 sequences were supported with telomeric repeats at the ends. GC-rich regions are located at the middle of the sequences and are enriched with long interspersed nuclear elements (LINEs). Repetitive sequences including long-terminal repeats (LTRs) and LINEs occupy 71.6% of the genome. A total of 21,887 potential protein-coding genes were predicted. The genomic data reported in this study served not only matsutake gene sequences but also genome structures and intergenic sequences. The information gained would be a great reference for exploring the genetics, genomics, and evolutionary study of matsutake in the future, and ultimately facilitate the conservation of this vulnerable genetic resource.

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.

Establishing an efficient protoplast transient expression system for investigation of floral thermogenesis in aroids.

KEY MESSAGE: Floral thermogenesis is an important reproductive strategy for attracting pollinators. We developed essential biological tools for studying floral thermogenesis using two species of thermogenic aroids, Symplocarpus renifolius and Alocasia odora. Aroids contain many species with intense heat-producing abilities in their inflorescences. Several genes have been proposed to be involved in thermogenesis of these species, but biological tools for gene functional analyses are lacking. In this study, we aimed to develop a protoplast-based transient expression (PTE) system for the study of thermogenic aroids. Initially, we focused on skunk cabbage (Symplocarpus renifolius) because of its ability to produce intense as well as durable heat. In this plant, leaf protoplasts were isolated from potted and shoot tip-cultured plants with high efficiency (ca. 1.0 × 105/g fresh weight), and more than half of these protoplasts were successfully transfected. Using this PTE system, we determined the protein localization of three mitochondrial energy-dissipating proteins, SrAOX, SrUCPA, and SrNDA1, fused to green fluorescent protein (GFP). These three GFP-fused proteins were localized in MitoTracker-stained mitochondria in leaf protoplasts, although the green fluorescent particles in protoplasts expressing SrUCPA-GFP were significantly enlarged. Finally, to assess whether the PTE system established in the leaves of S. renifolius is applicable for floral tissues of thermogenic aroids, inflorescences of S. renifolius and another thermogenic aroid (Alocasia odora) were used. Although protoplasts were successfully isolated from several tissues of the inflorescences, PTE systems worked well only for the protoplasts isolated from the female parts (slightly thermogenic or nonthermogenic) of A. odora inflorescences. Our developed system has a potential to be widely used in inflorescences as well as leaves in thermogenic aroids and therefore may be a useful biological tool for investigating floral thermogenesis.

Positive selective sweeps of epigenetic mutations regulating specialized metabolites in plants.

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana In comparison with single nucleotide polymorphisms, three types of methylation-methylated CGs (mCGs), mCHGs, and mCHHs-contributed highly to variable gene expression levels among an A thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

Draft Genome Sequences of Bifidobacterium animalis Consecutively Isolated from Healthy Japanese Individuals.

Bifidobacterium species are well recognized as probiotics and colonized in various parts of the human body. Here, we report the draft genome sequences of Bifidobacterium animalis isolated from two healthy Japanese volunteers, one of which was sampled twice before and after a 10-year interval. A core genome phylogeny analysis indicated that the strains isolated from the same volunteer were closely related. This paper is the first report of multiple draft genome sequences of B. animalis independently isolated from the same individual and provides insight into the probiotic potential of a member of this species.

Differential dynamics and impacts of prophages and plasmids on the pangenome and virulence factor repertoires of Shiga toxin-producing Escherichia coli O145:H28.

Phages and plasmids play important roles in bacterial evolution and diversification. Although many draft genomes have been generated, phage and plasmid genomes are usually fragmented, limiting our understanding of their dynamics. Here, we performed a systematic analysis of 239 draft genomes and 7 complete genomes of Shiga toxin (Stx)-producing Escherichia coli O145:H28, the major virulence factors of which are encoded by prophages (PPs) or plasmids. The results indicated that PPs are more stably maintained than plasmids. A set of ancestrally acquired PPs was well conserved, while various PPs, including Stx phages, were acquired by multiple sublineages. In contrast, gains and losses of a wide range of plasmids have frequently occurred across the O145:H28 lineage, and only the virulence plasmid was well conserved. The different dynamics of PPs and plasmids have differentially impacted the pangenome of O145:H28, with high proportions of PP- and plasmid-associated genes in the variably present and rare gene fractions, respectively. The dynamics of PPs and plasmids have also strongly impacted virulence gene repertoires, such as the highly variable distribution of stx genes and the high conservation of a set of type III secretion effectors, which probably represents the core effectors of O145:H28 and the genes on the virulence plasmid in the entire O145:H28 population. These results provide detailed insights into the dynamics of PPs and plasmids, and show the application of genomic analyses using a large set of draft genomes and appropriately selected complete genomes.

Draft Genome Sequences of Enterococcus faecalis Strains Isolated from Healthy Japanese Individuals.

Enterococcus faecalis is a common commensal of the intestines of humans and other mammals but is also a frequent cause of serious ailments. Here, we report 14 draft genome sequences of strains of Enterococcus faecalis, a normal inhabitant and Gram-positive bacterium that was isolated from 7 healthy Japanese volunteers.

Comparison of the sequencing bias of currently available library preparation kits for Illumina sequencing of bacterial genomes and metagenomes.

In bacterial genome and metagenome sequencing, Illumina sequencers are most frequently used due to their high throughput capacity, and multiple library preparation kits have been developed for Illumina platforms. Here, we systematically analysed and compared the sequencing bias generated by currently available library preparation kits for Illumina sequencing. Our analyses revealed that a strong sequencing bias is introduced in low-GC regions by the Nextera XT kit. The level of bias introduced is dependent on the level of GC content; stronger bias is generated as the GC content decreases. Other analysed kits did not introduce this strong sequencing bias. The GC content-associated sequencing bias introduced by Nextera XT was more remarkable in metagenome sequencing of a mock bacterial community and seriously affected estimation of the relative abundance of low-GC species. The results of our analyses highlight the importance of selecting proper library preparation kits according to the purposes and targets of sequencing, particularly in metagenome sequencing, where a wide range of microbial species with various degrees of GC content is present. Our data also indicate that special attention should be paid to which library preparation kit was used when analysing and interpreting publicly available metagenomic data.

Large-scale genome analysis of bovine commensal Escherichia coli reveals that bovine-adapted E. coli lineages are serving as evolutionary sources of the emergence of human intestinal pathogenic strains.

How pathogens evolve their virulence to humans in nature is a scientific issue of great medical and biological importance. Shiga toxin (Stx)-producing Escherichia coli (STEC) and enteropathogenic E. coli (EPEC) are the major foodborne pathogens that can cause hemolytic uremic syndrome and infantile diarrhea, respectively. The locus of enterocyte effacement (LEE)-encoded type 3 secretion system (T3SS) is the major virulence determinant of EPEC and is also possessed by major STEC lineages. Cattle are thought to be the primary reservoir of STEC and EPEC. However, genome sequences of bovine commensal E. coli are limited, and the emerging process of STEC and EPEC is largely unknown. Here, we performed a large-scale genomic comparison of bovine commensal E. coli with human commensal and clinical strains, including EPEC and STEC, at a global level. The analyses identified two distinct lineages, in which bovine and human commensal strains are enriched, respectively, and revealed that STEC and EPEC strains have emerged in multiple sublineages of the bovine-associated lineage. In addition to the bovine-associated lineage-specific genes, including fimbriae, capsule, and nutrition utilization genes, specific virulence gene communities have been accumulated in stx- and LEE-positive strains, respectively, with notable overlaps of community members. Functional associations of these genes probably confer benefits to these E. coli strains in inhabiting and/or adapting to the bovine intestinal environment and drive their evolution to highly virulent human pathogens under the bovine-adapted genetic background. Our data highlight the importance of large-scale genome sequencing of animal strains in the studies of zoonotic pathogens.

Alternative Oxidase Capacity of Mitochondria in Microsporophylls May Function in Cycad Thermogenesis.

Cone thermogenesis is a widespread phenomenon in cycads and may function to promote volatile emissions that affect pollinator behavior. Given their large population size and intense and durable heat-producing effects, cycads are important organisms for comprehensive studies of plant thermogenesis. However, knowledge of mitochondrial morphology and function in cone thermogenesis is limited. Therefore, we investigated these mitochondrial properties in the thermogenic cycad species Cycas revoluta Male cones generated heat even in cool weather conditions. Female cones produced heat, but to a lesser extent than male cones. Ultrastructural analyses of the two major tissues of male cones, microsporophylls and microsporangia, revealed the existence of a population of mitochondria with a distinct morphology in the microsporophylls. In these cells, we observed large mitochondria (cross-sectional area of 2 µm2 or more) with a uniform matrix density that occupied >10% of the total mitochondrial volume. Despite the size difference, many nonlarge mitochondria (cross-sectional area <2 µm2) also exhibited a shape and a matrix density similar to those of large mitochondria. Alternative oxidase (AOX) capacity and expression levels in microsporophylls were much higher than those in microsporangia. The AOX genes expressed in male cones revealed two different AOX complementary DNA sequences: CrAOX1 and CrAOX2 The expression level of CrAOX1 mRNA in the microsporophylls was 100 times greater than that of CrAOX2 mRNA. Collectively, these results suggest that distinctive mitochondrial morphology and CrAOX1-mediated respiration in microsporophylls might play a role in cycad cone thermogenesis.

Genomic Characterization of ß-Glucuronidase-Positive Escherichia coli O157:H7 Producing Stx2a.

Among Shiga toxin (Stx)-producing Escherichia coli (STEC) O157:H7 strains, those producing Stx2a cause more severe diseases. Atypical STEC O157:H7 strains showing a ß-glucuronidase-positive phenotype (GP STEC O157:H7) have rarely been isolated from humans, mostly from persons with asymptomatic or mild infections; Stx2a-producing strains have not been reported. We isolated, from a patient with bloody diarrhea, a GP STEC O157:H7 strain (PV15-279) that produces Stx2a in addition to Stx1a and Stx2c. Genomic comparison with other STEC O157 strains revealed that PV15-279 recently emerged from the stx1a/stx2c-positive GP STEC O157:H7 clone circulating in Japan. Major virulence genes are shared between typical (ß-glucuronidase-negative) and GP STEC O157:H7 strains, and the Stx2-producing ability of PV15-279 is comparable to that of typical STEC O157:H7 strains; therefore, PV15-279 presents a virulence potential similar to that of typical STEC O157:H7. This study reveals the importance of GP O157:H7 as a source of highly pathogenic STEC clones.

Population structure of Escherichia coli O26 : H11 with recent and repeated stx2 acquisition in multiple lineages.

A key virulence factor of enterohaemorrhagic Escherichia coli (EHEC) is the bacteriophage-encoded Shiga toxin (Stx). Stxs are classified into two types, Stx1 and Stx2, and Stx2-producing strains are thought to cause more severe infections than strains producing only Stx1. Although O26 : H11 is the second most prevalent EHEC following O157 : H7, the majority of O26 : H11 strains produce Stx1 alone. However, Stx2-producing O26 strains have increasingly been detected worldwide. Through a large-scale genome analysis, we present a global phylogenetic overview and evolutionary timescale for E. coli O26 : H11. The origin of O26 has been estimated to be 415 years ago. Sequence type 21C1 (ST21C1), one of the two sublineages of ST21, the most predominant O26 : H11 lineage worldwide, emerged 213 years ago from one of the three ST29 sublineages (ST29C2). The other ST21 lineage (ST21C2) emerged 95 years ago from ST21C1. Increases in population size occurred in the late 20th century for all of the O26 lineages, but most remarkably for ST21C2. Analysis of the distribution of stx2-positive strains revealed the recent and repeated acquisition of the stx2 gene in multiple lineages of O26, both in ST21 and ST29. Other major EHEC virulence genes, such as type III secretion system effector genes and plasmid-encoded virulence genes, were well conserved in ST21 compared to ST29. In addition, more antimicrobial-resistance genes have accumulated in the ST21C1 lineage. Although current attention is focused on several highly virulent ST29 clones that have acquired the stx2 gene, there is also a considerable risk that the ST21 lineage could yield highly virulent clones.

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.