Illumina Sequencing Data Research Articles

Comparative analysis of mitochondrial genomes of invasive weed Mikania micrantha and its indigenous congener Mikania cordata

Mikania micrantha and Mikania cordata are two distinct species in China. The former is notorious as one of the top 100 worst invasive species, whereas the latter is an indigenous species harmless to native plants or the environment. They form an ideal congener pair for comparative studies aimed at deeply understanding the invasion mechanisms of the exotic weed. In this study, we have assembled and annotated the mitogenomes of both species using Illumina and PacBio sequencing data and compared their characteristic differences. The complete mitogenome of M. micrantha is a double-stranded DNA with a length of 336,564 bp, while the mitogenome of M. cordata exhibits a branching structure, consisting of two small circular molecules and six linear molecules, with a combined length totaling 335,444 bp. Compared to M. cordata, M. micrantha has less SSRs, tandem repeats, dispersed repeats, mitochondrial protein coding genes (PCGs). The two plants show similar codon usage patterns. This comparative study has revealed the structure and function of the mitogenomes of the two species and laid a solid foundation for investigating the effects of gene loss and duplication on the development of invasive traits in M. micrantha.

De novo assembly and characterization of the complete mitochondrial genome of Phellodendron amurense reveals three repeat‐mediated recombination

Phellodendron amurense Rupr., a rare herb renowned for its medicinal and ecological significance, has remained genetically unexplored at the mitochondrial level until now. This study presents the first-ever systematic assembly and annotation of the complete mitochondrial genome of P. amurense, achieved through a hybrid strategy combining Illumina and Nanopore sequencing data. The mitochondrial genome spans 566,285 bp with a GC content of 45.51 %, structured into two circular molecules. Our comprehensive analysis identified 32 protein-coding genes (PCGs), 33 tRNA genes, and 3 rRNA genes, alongside 181 simple sequence repeats, 19 tandem repeats, and 310 dispersed repeats. Notably, multiple genome conformations were predicted due to repeat-mediated homologous recombination. Additionally, we assembled the chloroplast genome, identifying 21 mitochondrial plastid sequences that provide insights into organelle genome interactions. A total of 380 RNA-editing sites within the mitochondrial PCGs were predicted, enhancing our understanding of gene regulation and function. Phylogenetic analysis using mitochondrial PCGs from 30 species revealed evolutionary relationships, confirming the homology between P. amurense and Citrus species. This foundational study offers a valuable genetic resource for the Rutaceae family, facilitating further research into genetic evolution and molecular diversity in plant mitochondrial genomes.

The de novo assembly and characterization of the complete mitochondrial genome of bottle gourd (Lagenaria siceraria) reveals the presence of homologous conformations produced by repeat-mediated recombination.

Bottle gourd is an annual herbaceous plant that not only has high nutritional value and many medicinal applications but is also used as a rootstock for the grafting of cucurbit crops such as watermelon, cucumber and melon. Organellar genomes provide valuable resources for genetic breeding. A hybrid strategy with Illumina and Oxford Nanopore Technology sequencing data was used to assemble bottle gourd mitochondrial and chloroplast genomes. The length of the bottle gourd mitochondrial genome was 357547 bp, and that of the chloroplast genome was 157121 bp. These genomes had 27 homologous fragments, accounting for 6.50% of the total length of the bottle gourd mitochondrial genome. In the mitochondrial genome, 101 simple sequence repeats (SSRs) and 10 tandem repeats were identified. Moreover, 1 pair of repeats was shown to mediate homologous recombination into 1 major conformation and 1 minor conformation. The existence of these conformations was verified via PCR amplification and Sanger sequencing. Evolutionary analysis revealed that the mitochondrial genome sequence of bottle gourd was highly conserved. Furthermore, collinearity analysis revealed many rearrangements between the homologous fragments of Cucurbita and its relatives. The Ka/Ks values for most genes were between 0.3~0.9, which means that most of the genes in the bottle gourd mitochondrial genome are under purifying selection. We also identified a total of 589 potential RNA editing sites on 38 mitochondrial protein-coding genes (PCGs) on the basis of long noncoding RNA (lncRNA)-seq data. The RNA editing sites of nad1-2, nad4L-2, atp6-718, atp9-223 and rps10-391 were successfully verified via PCR amplification and Sanger sequencing. In conclusion, we assembled and annotated bottle gourd mitochondrial and chloroplast genomes to provide a theoretical basis for similar organelle genomic studies.

The complete chloroplast genome of Lindernia crustacea (L.) F. Muell 1882 (Linderniaceae) and its phylogenetic analysis

Lindernia crustacea (L.) F. Muell 1882, a species in the Linderniaceae family, holds traditional medicinal value in China. To investigate its genetic diversity, we assembled, annotated, and characterized the first complete chloroplast genome of L. crustacea using Illumina sequencing data and various bioinformatics tools. The genome is 153,647 bp in length, with a GC content of 37.6%. It exhibits a typical quadripartite structure, consisting of a large single-copy region (LSC) of 85,411 bp, a small single-copy region (SSC) of 18,724 bp, and two inverted repeat sequences (IRa and IRb) of 25,816 bp each. The genome was predicted to contain 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis indicated that L. crustacea is closely related to L. stricta. These findings provide a foundation for further research on the evolution and potential medicinal applications of the Linderniaceae family.

Complete mitochondrial genome of the edible mushroom Singerocybe alboinfundibuliformis (Clitocybaceae, Agaricales)

Singerocybe alboinfundibuliformis (Seok et al.) Yang, Qin & Takah 2014 is an edible mushroom distributed in several East or Southeast Asian countries. Herein, we report the mitochondrial genome of S. alboinfundibuliformis based on Illumina sequencing data. The overall length of the mitochondrial genome is 64,279 bp, with a GC content of 29.0%. It contains 14 typical protein-coding genes, 27 tRNA genes, two rRNA genes, and 13 intergenic ORFs. Most of these genes (39 out of 56) are transcribed at the forward strand, and few (17 out of 56) are transcribed at the reverse strand. Among these genes, only the rnl gene is invaded by an intron, and all other genes are intron-free. Phylogenetic analysis based on mitochondrial amino acid sequences supports the phylogenetic position of S. alboinfundibuliformis in Clitocybaceae, being close to Lepista sordida (Schumach.) Singer 1951. This study serves as a springboard for future investigation on fungal evolution in Clitocybaceae.

High-quality reference genome of cowpea beetle Callosobruchus maculatus

Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.

Graph-based mitochondrial genomes of three foundation species in the Saccharum genus.

We reported the graph-based mitochondrial genomes of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) for the first time. The results revealed pan-structural variation and evolutionary processes in the mitochondrial genomes within Saccharum. Saccharum belongs to the Andropogoneae, and cultivars species in Saccharum contribute nearly 80% of sugar production in the world. To explore the genomic studies in Saccharum, we assembled 15 complete mitochondrial genomes (mitogenome) of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) using Illumina and Oxford Nanopore Technologies sequencing data. The mitogenomes of the three species were divided into a total of eight types based on contig numbers and linkages. All mitogenomes in the three species encoded 51 unique genes, including 32 protein-coding, 3 ribosomal RNA (rRNA) and 16 transfer RNA (tRNA) genes. The existence of long and short-repeat-mediated recombinations in the mitogenome of S. officinarum and S. robustum was revealed and confirmed through PCR validation. Furthermore, employing comparative genomics and phylogenetic analyses of the organelle genomes, we unveiled the evolutionary relationships and history of the major interspecific lineages in Saccharum genus. Phylogenetic analyses of homologous fragments between S. officinarum and S. robustum showed that S. officinarum and S. robustum are phylogenetically distinct and that they were likely parallel rather than domesticated. The variations between ancient (S. sinense and S. barberi) and modern cultivated species (S. hybrid) possibly resulted from hybridization involving different S. officinarum accessions. Lastly, this project reported the first graph-based mitogenomes of three Saccharum species, and a systematic comparison of the structural organization, evolutionary processes, and pan-structural variation of the Saccharum mitogenomes revealed the differential features of the Saccharum mitogenomes.

Whole-genome assembly of a novel invertebrate herpesvirus from the gastropod Babylonia areolata.

Molluscan herpesviruses cause disease in species of major importance to aquaculture and are the only known herpesviruses to infect invertebrates, which lack an adaptive immune system. Understanding the evolution of malacoherpesviruses in relation to their hosts will likely require comparative genomic studies on multiple phylogenetic scales. Currently, only two malacoherpesvirus species have genomes that have been fully assembled, which limits the ability to perform comparative genomic studies on this family of viruses. In the present study, we fully assemble a herpesvirus from Illumina and Nanopore sequence data that were previously used to assemble the genome of the gastropod Babylonia areolata. We tentatively assign this novel herpesvirus to the genus Aurivirus within the family Malacoherpesviridae based on a phylogenetic analysis of DNA polymerase. While structurally similar to other malacoherpesvirus genomes, a synteny analysis of the novel herpesvirus with another Aurivirus species indicates that genomic rearrangements might be an important process in the evolution of this genus. We anticipate that future complete assemblies of malacoherpesviruses will be a valuable resource in comparative herpesvirus research.

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

The complete mitogenome of the pond wolf spider Pardosa pseudoannulata, with phylogenetic implications for the Lycosidae

The pond wolf spider Pardosa pseudoannulata Bösenberg & Strand, 1906 (Araneae: Lycosidae) is an important predator of agricultural pests in southern, eastern and southeastern Asia. Here, we report the complete mitogenome of this spider reconstructed from Illumina sequencing data. The circular mitogenome length is 14,533 bp with the nucleotide composition A (33.3%), C (8.2%), G (15.2%), and T (43.3%). The P. pseudoannulata mitogenome comprises 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region. Phylogenetic analyses of Lycosidae mitogenomes supported the monophyly of the subfamily Pardosinae and the two genera Pardosa and Alopecosa, and indicated the polyphyly of the subfamily Lycosinae and the paraphyly of its type genus Lycosa. In this study, P. pseudoannulata is the closest relative to P. pusiola. These results provide useful genetic information for future studies on the diversity, phylogeny, and evolution for wolf spiders.

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224–5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

Re-analysis of an outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with raw drinking milk using Nanopore sequencing

The aim of this study was to compare Illumina and Oxford Nanopore Technology (ONT) sequencing data to quantify genetic variation to assess within-outbreak strain relatedness and characterise microevolutionary events in the accessory genomes of a cluster of 23 genetically and epidemiologically linked isolates related to an outbreak of Shiga toxin-producing Escherichia coli O157:H7 caused by the consumption of raw drinking milk. There were seven discrepant variants called between the two technologies, five were false-negative or false-positive variants in the Illumina data and two were false-negative calls in ONT data. After masking horizontally acquired sequences such as prophages, analysis of both short and long-read sequences revealed the 20 isolates linked to the outbreak in 2017 had a maximum SNP distance of one SNP between each other, and a maximum of five SNPs when including three additional strains identified in 2019. Analysis of the ONT data revealed a 47 kbp deletion event in a terminal compound prophage within one sample relative to the remaining samples, and a 0.65 Mbp large chromosomal rearrangement (inversion), within one sample relative to the remaining samples. Furthermore, we detected two bacteriophages encoding the highly pathogenic Shiga toxin (Stx) subtype, Stx2a. One was typical of Stx2a-phage in this sub-lineage (Ic), the other was atypical and inserted into a site usually occupied by Stx2c-encoding phage. Finally, we observed an increase in the size of the pO157 IncFIB plasmid (1.6 kbp) in isolates from 2019 compared to those from 2017, due to the duplication of insertion elements within the plasmids from the more recently isolated strains. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the genome plasticity and virulence between strains of this zoonotic, foodborne pathogen.

Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum

BackgroundSchistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China, the Philippines and Indonesia. Oncomelania hupensis (Gastropoda: Pomatiopsidae) is the unique intermediate host of S. japonicum. A complete genome sequence of O. hupensis will enable the fundamental understanding of snail biology as well as its co-evolution with the S. japonicum parasite. Assembling a high-quality reference genome of O. hupehensis will provide data for further research on the snail biology and controlling the spread of S. japonicum.MethodsThe draft genome was de novo assembly using the long-read sequencing technology (PacBio Sequel II) and corrected with Illumina sequencing data. Then, using Hi-C sequencing data, the genome was assembled at the chromosomal level. CAFE was used to do analysis of contraction and expansion of the gene family and CodeML module in PAML was used for positive selection analysis in protein coding sequences.ResultsA total length of 1.46 Gb high-quality O. hupensis genome with 17 unique full-length chromosomes (2n = 34) of the individual including a contig N50 of 1.35 Mb and a scaffold N50 of 75.08 Mb. Additionally, 95.03% of these contig sequences were anchored in 17 chromosomes. After scanning the assembled genome, a total of 30,604 protein-coding genes were predicted. Among them, 86.67% were functionally annotated. Further phylogenetic analysis revealed that O. hupensis was separated from a common ancestor of Pomacea canaliculata and Bellamya purificata approximately 170 million years ago. Comparing the genome of O. hupensis with its most recent common ancestor, it showed 266 significantly expanded and 58 significantly contracted gene families (P < 0.05). Functional enrichment of the expanded gene families indicated that they were mainly involved with intracellular, DNA-mediated transposition, DNA integration and transposase activity.ConclusionsIntegrated use of multiple sequencing technologies, we have successfully constructed the genome at the chromosomal-level of O. hupensis. These data will not only provide the compressive genomic information, but also benefit future work on population genetics of this snail as well as evolutional studies between S. japonicum and the snail host.Graphical

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli.

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Do we still need Illumina sequencing data? Evaluating Oxford Nanopore Technologies R10.4.1 flow cells and the Rapid v14 library prep kit for Gram negative bacteria whole genome assemblies.

The best whole genome assemblies are currently built from a combination of highly accurate short-read sequencing data and long-read sequencing data that can bridge repetitive and problematic regions. Oxford Nanopore Technologies (ONT) produce long-read sequencing platforms and they are continually improving their technology to obtain higher quality read data that is approaching the quality obtained from short-read platforms such as Illumina. As these innovations continue, we evaluated how much ONT read coverage produced by the Rapid Barcoding Kit v14 (SQK-RBK114) is necessary to generate high-quality hybrid and long-read-only genome assemblies for a panel of carbapenemase-producing Enterobacterales bacterial isolates. We found that 30× long-read coverage is sufficient if Illumina data are available, and that more (at least 100× long-read coverage is recommended for long-read-only assemblies. Illumina polishing is still improving single nucleotide variants (SNVs) and INDELs in long-read-only assemblies. We also examined if antimicrobial resistance genes could be accurately identified in long-read-only data, and found that Flye assemblies regardless of ONT coverage detected>96% of resistance genes at 100% identity and length. Overall, the Rapid Barcoding Kit v14 and long-read-only assemblies can be an optimal sequencing strategy (i.e., plasmid characterization and AMR detection) but finer-scale analyses (i.e., SNV) still benefit from short-read data.

PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data.

Complete mitochondrial genomes (mitogenomes) of plants are valuable resources for nucleocytoplasmic interactions, plant evolution, and plant cytoplasmic male sterile line breeding. However, the complete assembly of plant mitogenomes is challenging due to frequent recombination events and horizontal gene transfers. Previous studies have adopted Illumina, PacBio, and Nanopore sequencing data to assemble plant mitogenomes, but the poor assembly completeness, low sequencing accuracy, and high cost limit the sampling capacity. Here, we present an efficient assembly toolkit (PMAT) for de novo assembly of plant mitogenomes using low-coverage HiFi sequencing data. PMAT has been applied to the de novo assembly of 13 broadly representative plant mitogenomes, outperforming existing organelle genome assemblers in terms of assembly accuracy and completeness. By evaluating the assembly of plant mitogenomes from different sequencing data, it was confirmed that PMAT only requires 1× HiFi sequencing data to obtain a complete plant mitogenome. The source code for PMAT is available at https://github.com/bichangwei/PMAT. The developed PMAT toolkit will indeed accelerate the understanding of evolutionary variation and breeding application of plant mitogenomes.

Comparative plastome assembly of the yellow ironweed (Verbesina alternifolia) using Nanopore and Illumina reads.

Chloroplast genomes (plastomes) represent a very important source of valuable information for phylogenetic and biogeographic reconstructions. The use of short reads (as those produced from Illumina sequencing), along with de novo read assembly, has been considered the "gold standard" for plastome reconstruction. However, short reads often cannot reconstruct long repetitive regions in chloroplast genomes. Long Nanopore (ONT) reads can help bridging long repetitive regions but are by far more error-prone than those produced by Illumina sequencing. Verbesina is the largest genus of tribe Heliantheae (Asteraceae) and includes species of economic importance as ornamental or as invasive weeds. However, no complete chloroplast genomes have been published yet for the genus. We utilized Illumina and Nanopore sequencing data and different assembly strategies to reconstruct the plastome of Verbesina alternifolia and evaluated the usefulness of the Nanopore assemblies. The two plastome sequence assemblages, one obtained with the Nanopore sequencing and the other inferred with Illumina reads, were identical, except for missing bases in homonucleotide regions. The best-assembled plastome of V. alternifolia was 152,050 bp in length and contained 80, 29, and four unique protein-coding genes, tRNAs, and rRNAs, respectively. When used as reference for mapping Illumina reads, all plastomes performed similarly. In a phylogenetic analysis including 28 other plastomes from closely related taxa (from the Heliantheae alliance), the two Verbesina chloroplast genomes grouped together and were nested among the other members of the tribe Heliantheae s.str. Our study highlights the usefulness of the Nanopore technology for assembling rapidly and cost-effectively chloroplast genomes, especially in taxonomic groups with paucity of publicly available plastomes.

Assembly and comparative genome analysis of four mitochondrial genomes from Saccharum complex species.

Saccharum complex includes genera Saccharum, Miscanthus, Erianthus, Narenga, and Tripidium. Since the Saccharum complex/Saccharinae constitutes the gene pool used by sugarcane breeders to introduce useful traits into sugarcane, studying the genomic characterization of the Saccharum complex has become particularly important. Here, we assembled graph-based mitochondrial genomes (mitogenomes) of four Saccharinae species (T. arundinaceum, E. rockii, M. sinensis, and N. porphyrocoma) using Illumina and PacBio sequencing data. The total lengths of the mitogenomes of T. arundinaceum, M. sinensis, E. rockii and N. porphyrocoma were 549,593 bp, 514,248 bp, 481,576 bp and 513,095 bp, respectively. Then, we performed a comparative mitogenomes analysis of Saccharinae species, including characterization, organelles transfer sequence, collinear sequence, phylogenetics analysis, and gene duplicated/loss. Our results provided the mitogenomes of four species closely related to sugarcane breeding, enriching the mitochondrial genomic resources of the Saccharinae. Additionally, our study offered new insights into the evolution of mitogenomes at the family and genus levels and enhanced our understanding of organelle evolution in the highly polyploid Saccharum genus.

High-Risk Sequence Prediction Model in DNA Storage: The LQSF Method.

Traditional DNA storage technologies rely on passive filtering methods for error correction during synthesis and sequencing, which result in redundancy and inadequate error correction. Addressing this, the Low Quality Sequence Filter (LQSF) was introduced, an innovative method employing deep learning models to predict high-risk sequences. The LQSF approach leverages a classification model trained on error-prone sequences, enabling efficient pre-sequencing filtration of low-quality sequences and reducing time and resources in subsequent stages. Analysis has demonstrated a clear distinction between high and low-quality sequences, confirming the efficacy of the LQSF method. Extensive training and testing were conducted across various neural networks and test sets. The results showed all models achieving an AUC value above 0.91 on ROC curves and over 0.95 on PR curves across different datasets. Notably, models such as Alexnet, VGG16, and VGG19 achieved a perfect AUC of 1.0 on the Original dataset, highlighting their precision in classification. Further validation using Illumina sequencing data substantiated a strong correlation between model scores and sequence error-proneness, emphasizing the model's applicability. The LQSF method marks a significant advancement in DNA storage technology, introducing active sequence filtering at the encoding stage. This pioneering approach holds substantial promise for future DNA storage research and applications.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.