Universal Single-Copy Orthologs Research Articles

The gyrfalcon (Falco rusticolus) genome.

High-quality genome assemblies are characterized by high-sequence contiguity, completeness, and a low error rate, thus providing the basis for a wide array of studies focusing on natural species ecology, conservation, evolution, and population genomics. To provide this valuable resource for conservation projects and comparative genomics studies on gyrfalcon (Falco rusticolus), we sequenced and assembled the genome of this species using third-generation sequencing strategies and optical maps. Here, we describe a highly contiguous and complete genome assembly comprising 20 scaffolds and 13 contigs with a total size of 1.193 Gbp, including 8,064 complete Benchmarking Universal Single-Copy Orthologs (BUSCOs) of the total 8,338 BUSCO groups present in the library aves_odb10. Of these BUSCO genes, 96.7% were complete, 96.1% were present as a single copy, and 0.6% were duplicated. Furthermore, 0.8% of BUSCO genes were fragmented and 2.5% (210) were missing. A de novo search for transposable elements (TEs) identified 5,716 TEs that masked 7.61% of the F. rusticolus genome assembly when combined with publicly available TE collections. Long interspersed nuclear elements, in particular, the element Chicken-repeat 1 (CR1), were the most abundant TEs in the F. rusticolus genome. A de novo first-pass gene annotation was performed using 293,349 PacBio Iso-Seq transcripts and 496,195 transcripts derived from the assembly of 42,429,525 Illumina PE RNA-seq reads. In all, 19,602 putative genes, of which 59.31% were functionally characterized and associated with Gene Ontology terms, were annotated. A comparison of the gyrfalcon genome assembly with the publicly available assemblies of the domestic chicken (Gallus gallus), zebra finch (Taeniopygia guttata), and hummingbird (Calypte anna) revealed several genome rearrangements. In particular, nine putative chromosome fusions were identified in the gyrfalcon genome assembly compared with those in the G. gallus genome assembly. This genome assembly, its annotation for TEs and genes, and the comparative analyses presented, complement and strength the base of high-quality genome assemblies and associated resources available for comparative studies focusing on the evolution, ecology, and conservation of Aves.

Genome Sequence Resource of Ustilago crameri, a Fungal Pathogen Causing Millet Smut Disease of Foxtail Millet.

Genome Sequence Resource of Ustilago crameri, a Fungal Pathogen Causing Millet Smut Disease of Foxtail Millet.

Long-Read Genome Sequencing of Abscondita cerata (Coleoptera: Lampyridae), the Endemic Firefly of Taiwan.

Abscondita cerata is the most abundant and widely distributed endemic firefly species in Taiwan and is considered a key environmental and ecological indicator organism. In this study, we report the first long-read genome sequencing of Abs. cerata sequenced by Nanopore technology. The draft genome size, 967 Mb, was measured through a hybrid approach that consisted of assembling using 11.25-Gb Nanopore long reads and polishing using 9.47-Gb BGI PE100 short reads. The drafted genome was assembled into 4,855 contigs, with the N50 reaching 325.269 kb length. The assembled genome was predicted to possess 55,206 protein-coding genes, of which 20,862 (37.78%) were functionally annotated with public databases. 47.11% of the genome sequences consisted of repeat elements; among them DNA transposons accounted for the largest proportion (26.79%). A BUSCO (Benchmarking Universal Single Copy Orthologs) evaluation demonstrated that the genome and gene completeness were 84.8% and 79%, respectively. The phylogeny constructed using 1,792 single copy genes was consistent with previous studies. The comparative transcriptome between adult male head and lantern tissues revealed (1) the vision of Abs. cerata is primarily UV-sensitive to environmental twilight, which determines when it begins its nocturnal activity, (2) the major expressed OR56d receptor may be correlated to suitable humidity sensing, and (3) Luc1-type luciferase is responsible for Abs. cerata's luminescent spectrum.

Whole-Genome Resource of Lasiodiplodia pseudotheobromae BaA, the Causative Agent of Black Root Rot Morinda officinalis.

Whole-Genome Resource of Lasiodiplodia pseudotheobromae BaA, the Causative Agent of Black Root Rot Morinda officinalis.

Draft genome of the medicinal tea tree Melaleuca alternifolia

Melaleuca alternifolia is a commercially important medicinal tea tree native to Australia. Tea tree oil, the essential oil distilled from its branches and leaves, has broad-spectrum germicidal activity and is highly valued in the pharmaceutical and cosmetic industries. Thus, the study of genome, which can provide reference for the investigation of genes involved in terpinen-4-ol biosynthesis, is quite crucial for improving the productivity of Tea tree oil. In our study, the next-generation sequencing was used to investigate the whole genome of Melaleuca alternifolia. About 114Gb high quality sequence data were obtained and assembled into 1,838,159 scafolds with an N50 length of 1021bp. The assembled genome size is about 595Mb, twice of that predicted by flow cytometer (300Mb) and k-mer analysis (345Mb). Benchmarking Universal Single-Copy Orthologs analyses indicated that only 11.3% of the conserved single-copy genes were miss. Repetitive regions cover over 40.43% of the genome. A total of 44,369 protein-coding genes were predicted and annotated against Nr, Swissprot, Refseq, COG, KOG, and KEGG database. Among these genes, 32,909 and 16,241 genes were functionally annotated in Nr and KEGG, respectively. Moreover, 29,411 and 14,435 genes were functionally annotated in COG and KOG. Additionally, 457,661 simple sequence repeats and 1109 transcription factors (TFs) form 67 TF families were identified in the assembled genome. Our findings provide a draft genome sequencing of M. alternifolia which can act as a reference for the deep sequencing strategies, and are useful for future functional and comparative genomics analyses.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus).

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

A draft genome assembly of ectoparasitic mite Pyemotes zhonghuajia Yu, Zhang & He (Prostigmata: Pyemotidae)

Pyemotes zhonghuajia Yu, Zhang & He (Prostigmata: Pyemotidae) is an ectoparasitic mite species with a variety of insects-specific protein toxins, which are harmless to mammals and human. This mite can cause the death of many agricultural and forestry pests which can be 680,000 times greater than its own body weight (Tian et al., 2020; Chen et al., 2021). Exploring its genome will help us understand its parasitic mechanism. PacBio technology was used to sequence the genome of P. zhonghuajia. The genome size of P. zhonghuajia was estimated to be 71.943 Mb, which captured 90.6% of the Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness (Song et al., 2022). A total of 11,183 protein coding genes were annotated with 91.2% BUSCO completeness. Maximum Likelihood phylogenetic tree was constructed using 473 single-copy genes of 12 species including two species of spiders, two species of scorpions and 10 species of mites, of which spiders and scorpions were outgroups. The result shows that Pyemotes zhonghuajia is sister to Tetranychus urticae. An estimate of species divergence time showed that P. zhonghuajia and T. urticae emerged during Triassic (223.43~247.66 Mya). In P. zhonghuajia, there were 621 gene families (containing 1,224 genes) expanded, and 3664 gene families (containing 3,927 genes) contracted. Most importantly, we found that P. zhonghuajia contained three types of toxin proteins, which were neurotoxins, dermonecrotic toxins and agatoxins. Neurotoxin and dermonecrotic toxin-related genes are significantly expanded. The toxin-related gene expression level of one-day-old mites was higher than that of seven-day-pregnant mites, which might explain why one one-day-old mite had such strongly toxin (Chen et al., 2021). We also verified the function of the toxin proteins and explored its effect on insect cells to lay a theoretical foundation for the study and development of biopesticides with multi-target insecticidal toxins (Song et al., 2022).

The Assembled Genome of the Stroke-Prone Spontaneously Hypertensive Rat.

We report the creation and evaluation of a de novo assembly of the genome of the spontaneously hypertensive rat, the most widely used model of human cardiovascular disease. The genome is assembled from long read sequencing (PacBio HiFi and continuous long read data [CLR]) and scaffolded with long-range structural information obtained from Bionano optical maps and proximity ligation sequencing proximity analysis of the genome. The genome assembly was polished with Illumina short reads. Completeness of the assembly was investigated using Benchmarking Universal Single Copy Orthologs analysis. The genome assembly was also evaluated with the rat reference gene set, using NCBI automated protocols. We also generated orthogonal single molecule transcript sequence reads (Iso-Seq) from 8 tissues and used them to validate the coding assembly, to annotate the assembly with RNA transcripts representing unique full length transcript isoforms for each gene and to determine whether divergences between RefSeq sequences and the assembly were attributable to assembly errors or polymorphisms. The assembly analysis indicates that this assembly is comparable in contiguity and completeness to the current rat reference assembly, while the use of HiFi sequencing yields an assembly that is more correct at the single base level. Synteny analysis was performed to uncover the extent of synteny and the presence and distribution of chromosomal rearrangements between the reference and this assembly. The resulting genome assembly is reference quality and captures significant structural variation.

Draft genome of six Cuban Anolis lizards and insights into genetic changes during their diversification

BackgroundDetecting genomic variants and their accumulation processes during species diversification and adaptive radiation is important for understanding the molecular and genetic basis of evolution. Anolis lizards in the West Indies are good models for studying evolutionary mechanisms because of the repeated evolution of their morphology and the ecology. We performed de novo genome assembly of six Cuban Anolis lizards with different ecomorphs and thermal habitats (Anolis isolepis, Anolis allisoni, Anolis porcatus, Anolis allogus, Anolis homolechis, and Anolis sagrei). We carried out a comparative analysis of these genome assemblies to investigate the genetic changes that occurred during their diversification.ResultsWe reconstructed novel draft genomes with relatively long scaffolds and high gene completeness, with the scaffold N50 ranging from 5.56 to 39.79 Mb and vertebrate Benchmarking Universal Single-Copy Orthologs completeness ranging from 77.5% to 86.9%. Comparing the repeat element compositions and landscapes revealed differences in the accumulation process between Cuban trunk-crown and trunk-ground species and separate expansions of several families of LINE in each Cuban trunk-ground species. Duplicated gene analysis suggested that the proportional differences in duplicated gene numbers among Cuban Anolis lizards may be associated with differences in their habitat ranges. Additionally, Pairwise Sequentially Markovian Coalescent analysis suggested that the effective population sizes of each species may have been affected by Cuba’s geohistory.ConclusionsWe provide draft genomes of six Cuban Anolis lizards and detected species and lineage-specific transposon accumulation and gene copy number changes that may be involved in adaptive evolution. The change processes in the past effective population size was also estimated, and the factors involved were inferred. These results provide new insights into the genetic basis of Anolis lizard diversification and are expected to serve as a stepping stone for the further elucidation of their diversification mechanisms.

Combination of long-read and short-read sequencing provides comprehensive transcriptome and new insight for Chrysanthemum morifolium ray-floret colorization

Chrysanthemum morifolium is one of the most popular ornamental plants globally. Owing to its large and complex genome (around 10 Gb, segmental hexaploid), it has been difficult to obtain comprehensive transcriptome, which will promote to perform new breeding technique, such as genome editing, in C. morifolium. In this study, we used single-molecule real-time (SMRT) sequencing and RNA-seq technologies, combined them with an error-correcting process, and obtained high-coverage ray-floret transcriptome. The SMRT-seq data increased the ratio of long mRNAs containing complete open-reading frames, and the combined dataset provided a more complete transcriptomic data than those produced from either SMRT-seq or RNA-seq-derived transcripts. We finally obtained ‘Sei Arabella’ transcripts containing 928,645 non-redundant mRNA, which showed 96.6% Benchmarking Universal Single-Copy Orthologs (BUSCO) score. We also validated the reliability of the dataset by analyzing a mapping rate, annotation and transcript expression. Using the dataset, we searched anthocyanin biosynthesis gene orthologs and performed a qRT-PCR experiment to assess the usability of the dataset. The assessment of the dataset and the following analysis indicated that our dataset is reliable and useful for molecular biology. The combination of sequencing methods provided genetic information and a way to analyze the complicated C. morifolium transcriptome.

Draft genome of the bluefin tuna blood fluke, Cardicola forsteri.

The blood fluke Cardicola forsteri (Trematoda: Aporocotylidae) is a pathogen of ranched bluefin tuna in Japan and Australia. Genomics of Cardicola spp. have thus far been limited to molecular phylogenetics of select gene sequences. In this study, sequencing of the C. forsteri genome was performed using Illumina short-read and Oxford Nanopore long-read technologies. The sequences were assembled de novo using a hybrid of short and long reads, which produced a high-quality contig-level assembly (N50 > 430 kb and L50 = 138). The assembly was also relatively complete and unfragmented, comprising 66% and 7.2% complete and fragmented metazoan Benchmarking Universal Single-Copy Orthologs (BUSCOs), respectively. A large portion (> 55%) of the genome was made up of intergenic repetitive elements, primarily long interspersed nuclear elements (LINEs), while protein-coding regions cover > 6%. Gene prediction identified 8,564 hypothetical polypeptides, > 77% of which are homologous to published sequences of other species. The identification of select putative proteins, including cathepsins, calpains, tetraspanins, and glycosyltransferases is discussed. This is the first genome assembly of any aporocotylid, a major step toward understanding of the biology of this family of fish blood flukes and their interactions within hosts.

Genomic novelty within a "great speciator" revealed by a high-quality reference genome of the collared kingfisher (Todiramphus chloris collaris).

Islands are natural laboratories for studying patterns and processes of evolution. Research on island endemic birds has revealed elevated speciation rates and rapid phenotypic evolution in several groups (e.g. white-eyes, Darwin's finches). However, understanding the evolutionary processes behind these patterns requires an understanding of how genotypes map to novel phenotypes. To date, there are few high-quality reference genomes for species found on islands. Here, we sequence the genome of one of Ernst Mayr's "great speciators," the collared kingfisher (Todiramphus chloris collaris). Utilizing high molecular weight DNA and linked-read sequencing technology, we assembled a draft high-quality genome with highly contiguous scaffolds (scaffold N50 = 19 Mb). Based on universal single-copy orthologs, we estimated a gene space completeness of 96.6% for the draft genome assembly. The population demographic history analyses reveal a distinct pattern of contraction and expansion in population size throughout the Pleistocene. Comparative genomic analysis of gene family evolution revealed that species-specific and rapidly expanding gene families in the collared kingfisher (relative to other Coraciiformes) are mainly involved in the ErbB signaling pathway and focal adhesion. Todiramphus kingfishers are a species-rich group that has become a focus of speciation research. This draft genome will be a platform for future taxonomic, phylogeographic, and speciation research in the group. For example, target genes will enable testing of changes in sensory structures associated with changes in vision and taste genes across kingfishers.

Draft genome of the lowland anoa (Bubalus depressicornis) and comparison with buffalo genome assemblies (Bovidae, Bubalina).

Genomic data for wild species of the genus Bubalus (Asian buffaloes) are still lacking while several whole genomes are currently available for domestic water buffaloes. To address this, we sequenced the genome of a wild endangered dwarf buffalo, the lowland anoa (Bubalus depressicornis), produced a draft genome assembly and made comparison to published buffalo genomes. The lowland anoa genome assembly was 2.56 Gbp long and contained 103,135 contigs, the longest contig being 337.39 kbp long. N50 and L50 values were 38.73 and 19.83 kbp, respectively, mean coverage was 44× and GC content was 41.74%. Two strategies were adopted to evaluate genome completeness: (1) determination of genomic features with de novo and homology-based predictions using annotations of chromosome-level genome assembly of the river buffalo and (2) employment of benchmarking against universal single-copy orthologs (BUSCO). Homology-based predictions identified 94.51% complete and 3.65% partial genomic features. De novo gene predictions identified 32,393 genes, representing 97.14% of the reference's annotated genes, whilst BUSCO search against the mammalian orthologs database identified 71.1% complete, 11.7% fragmented, and 17.2% missing orthologs, indicating a good level of completeness for downstream analyses. Repeat analyses indicated that the lowland anoa genome contains 42.12% of repetitive regions. The genome assembly of the lowland anoa is expected to contribute to comparative genome analyses among bovid species.

Reference quality genome sequence of Indian pomegranate cv. 'Bhagawa' (Punica granatum L.).

Pomegranate is an important fruit crop for ensuring livelihood and nutrition security in fragile semi-arid regions of the globe having limited irrigation resources. This is a high-value, nutritionally rich, and export-oriented agri-commodity that ensures high returns on investment to growers across the world. Although it is a valuable fruit crop, it has received only a limited genomics research outcome. To fast-track the pomegranate improvement program, de novo whole-genome sequencing of the main Indian cultivar ‘Bhagawa’ was initiated by the Indian Council of Agricultural Research–National Research Center on Pomegranate (ICAR–NRCP). We have demonstrated that a combination of commercially available technologies from Illumina, PacBio, 10X Genomics, and BioNano Genomics could be used efficiently for sequencing and reference-grade de novo assembly of the pomegranate genome. The research led to a final reference-quality genome assembly for ‘Bhagawa’ of 346.08 Mb in 342 scaffolds and an average N50 of 16.12 Mb and N90 of 1088.62 Kb. This assembly covered more than 98% of the estimated pomegranate genome size, 352.54 Mb. The LTR assembly index (LAI) value of 10 and 93.68% Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score over the 1,440 ortholog genes of the completed pomegranate genome indicates the quality of the assembled pomegranate genome. Furthermore, 29,435 gene models were discovered with a mean transcript length of 2,954 bp and a mean coding sequence length 1,090 bp. Four transcript data samples of pomegranate tissues were mapped over the assembled ‘Bhagawa’ genome up to 95% significant matches, indicating the high quality of the assembled genome. We have compared the ‘Bhagawa’ genome with the genomes of the pomegranate cultivars ‘Dabenzi’ and ‘Taishanhong.’ We have also performed whole-genome phylogenetic analysis using Computational Analysis of Gene Family Evolution (CAFE) and found that Eucalyptus grandis and pomegranate diverged 64 (60–70) million years ago. About 1,573 protein-coding resistance genes identified in the ‘Bhagawa’ genome were classified into 32 domains. In all, 314 copies of miRNA belonging to 26 different families were identified in the ‘Bhagawa’ genome. The reference-quality genome assembly of ‘Bhagawa’ is certainly a significant genomic resource for accelerated pomegranate improvement.

Genomic diversity and molecular epidemiology of a multidrug-resistant Pseudomonas aeruginosa DMC30b isolated from a hospitalized burn patient in Bangladesh

Objectives: Pseudomonas aeruginosa is a key opportunistic pathogen causing a wide range of community- and hospital-acquired infections in immunocompromised or catheterized patients. Here, we report the complete genome sequence of a multidrug-resistant (MDR) P. aeruginosa DMC30b to elucidate the genetic diversity, molecular epidemiology, and underlying mechanisms for antimicrobial resistance and virulence. Methods: P. aeruginosa DMC30b was isolated from septic wound swab of a severe burn patient. Whole-genome sequencing was performed under Ion Torrent platform. The genome was assembled using the SPAdes v. 3.12.01 in an integrated Genome Analysis Platform for Ion Torrent sequence data. The genome was annotated using the National Center for Biotechnology Information (NCBI) Prokaryotic Genome Annotation Pipeline. In-silico predictions of antimicrobial resistance genes, virulence factor genes, and metabolic functional potentials were performed using different curated bioinformatics tools. Results: P. aeruginosa DMC30b was found as a MDR strain and belonged to sequence type 244 (ST244). The complete genome size is 6 994 756 bp with a coverage of 76.76x, guanine-cytosine content of 65.7% and a Benchmarking Universal Single-Copy Orthologs score of 100. The genome of P. aeruginosa DMC30b harboured two predicted plasmid replicons (e,g. IncP-6; 78 007 bp and ColRNAI; 9359 bp), 35 resistomes (antimicrobial resistance genes) conferring resistance to 18 different antibiotics (including four beta-lactam classes), and 214 virulence factor genes. It was identified as the 167th ST244 strain among ∼ 5800 whole-genome sequences of P. aeruginosa available in the NCBI database. Conclusion: The MDR P. aeruginosa DMC30b was identified as the 167th ST244 complete genome to be submitted to the NCBI, and the first ST244 isolate sequenced from Bangladesh. The complete genome data with high genetic diversity and underlying mechanisms for antimicrobial resistance and virulence of P. aeruginosa DMC30b will aid in understanding the evolution and phylogeny of such high-risk clones and provide a solid basis for further research on MDR or extensively drug resistant strains.

Comparative Genome and Evolution Analyses of an Endangered Stony Coral Species Dendrophyllia cribrosa Near Dokdo Islands in the East Sea.

Stony corals often harbor intracellular photosynthetic dinoflagellate algae that receive dissolved inorganic nutrients. However, Dendrophyllia cribrosa is a nonsymbiotic stony coral distributed in the western Pacific. We assembled a chromosome-level D. cribrosa genome using PacBio and Hi-C technologies. The final assembly was 625 Mb, distributed on 14 chromosomes, and contained 30,493 protein-coding genes. The Benchmarking Universal Single-Copy Orthologs analysis revealed a percentage of 96.8 of the metazoan genome. A comparative phylogenetic analysis revealed that D. cribrosa, which lacks symbionts, evolved to acquire cellular energy by expanding genes related to acyl-CoA metabolism and carbohydrate transporters. This species also has expanded immune-related genes involved in the receptor protein tyrosine kinase signaling pathway. In addition, we observed a specific expansion of calcification genes, such as coral acid-rich proteins and carbonic anhydrase, in D. cribrosa. This high-quality reference genome and comparative analysis provides insights into the ecology and evolution of nonsymbiotic stony corals.

Genome assembly and annotation of the European earwig Forficula auricularia (subspecies B)

The European earwig Forficula auricularia is an important model for studies of maternal care, sexual selection, sociality, and host–parasite interactions. However, detailed genetic investigations of this species are hindered by a lack of genomic resources. Here, we present a high-quality hybrid genome assembly for Forficula auricularia using Nanopore long-reads and 10× linked-reads. The final assembly is 1.06 Gb in length with 31.03% GC content. It consists of 919 scaffolds with an N50 of 12.55 Mb. Half of the genome is present in only 20 scaffolds. Benchmarking Universal Single-Copy Orthologs scores are ∼90% from 3 sets of single-copy orthologs (eukaryotic, insect, and arthropod). The total repeat elements in the genome are 64.62%. The MAKER2 pipeline annotated 12,876 protein-coding genes and 21,031 mRNAs. Phylogenetic analysis revealed the assembled genome as that of species B, one of the 2 known genetic subspecies of Forficula auricularia. The genome assembly, annotation, and associated resources will be of high value to a large and diverse group of researchers working on dermapterans.

The genome of the rice planthopper egg parasitoid wasps Anagrus nilaparvatae casts light on the chemo- and mechanosensation in parasitism

BackgroundMymaridae is an ancient insect group and is a basal lineage of the superfamily Chalcidoidea. Species of Mymaridae have great potential for biological control. Anagrus nilaparvatae, a representative species of Mymaridae, is ideal for controlling rice planthopper due to its high rate of parasitism and ability to find hosts efficiently in paddy ridges and fields.ResultsUsing both PacBio single-molecule real-time and Illumina sequencing, we sequenced and assembled the whole genome of A. nilaparvatae, a first for the family Mymaridae. The assembly consists of 394 scaffolds, totaling 488.8 Mb. The assembly is of high continuity and completeness, indicated by the N50 value of 25.4 Mb and 98.2% mapping rate of Benchmarking Universal Single-Copy Orthologs. In total, 16,894 protein-coding genes in the genome were annotated. A phylogenomic tree constructed for A. nilaparvatae and other 12 species of Hymenoptera confirmed that the family Mymaridae is sister to all remaining chalcidoids. The divergence time between A. nilaparvatae and the other seven Chalcidoidea species was dated at ~ 126.9 Mya. Chemoreceptor and mechanoreceptor genes are important in explaining parasitic behavior. We identified 17 odorant binding proteins, 11 chemosensory proteins, four Niemann-Pick type C2 proteins, 88 olfactory receptors, 12 gustatory receptors, 22 ionotropic receptors and 13 sensory neuron membrane proteins in the genome of A. nilaparvatae, which are associated with the chemosensory functions. Strikingly, there is only one pickpocket receptors and nine transient receptor potential genes in the genome that have a mechanosensory function.ConclusionsWe obtained a high-quality genome assembly for A. nilaparvatae using PacBio single-molecule real-time sequencing, which provides phylogenomic insights for its evolutionary history. The small numbers of chemo- and mechanosensory genes in A. nilaparvatae indicate the species-specific host detection and oviposition behavior of A. nilaparvatae might be regulated by relatively simple molecular pathways.

A chromosome-level reference genome of a Convolvulaceae species Ipomoea cairica.

Ipomoea cairica is a perennial creeper that has been widely introduced as a garden ornamental across tropical, subtropical, and temperate regions. Because it grows extremely fast and spreads easily, it has been listed as an invasive species in many countries. Here, we constructed the chromosome-level reference genome of Ipomoea cairica by Pacific Biosciences HiFi and Hi-C sequencing, with the assembly size of 733.0 Mb, the contig N50 of 43.8 Mb, the scaffold N50 of 45.7 Mb, and the Benchmarking Universal Single-Copy Orthologs complete rate of 98.0%. Hi-C scaffolding assigned 97.9% of the contigs to 15 pseudo-chromosomes. Telomeric repeat analysis reveals that 7 of the 15 pseudo-chromosomes are gapless and telomere to telomere. The transposable element content of Ipomoea cairica is 73.4%, obviously higher than that of other Ipomoea species. A total of 38,115 protein-coding genes were predicted, with the Benchmarking Universal Single-Copy Orthologs complete rate of 98.5%, comparable to that of the genome assembly, and 92.6% of genes were functional annotated. In addition, we identified 3,039 tRNA genes and 2,403 rRNA genes in the assembled genome. Phylogenetic analysis showed that Ipomoea cairica formed a clade with Ipomoea aquatica, and they diverged from each other 8.1 million years ago. Through comparative genome analysis, we reconfirmed that a whole genome triplication event occurred specific to Convolvulaceae family and in the ancestor of the genus Ipomoea and Cuscuta. This high-quality reference genome of Ipomoea cairica will greatly facilitate the studies on the molecular mechanisms of its rapid growth and invasiveness.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus.

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.