Chromosome‐level genome assembly of a cyprinid fish Onychostoma macrolepis by integration of nanopore sequencing,Bionano and Hi‐C technology

Onychostoma macrolepis is an emerging commercial cyprinid fish species. It is a model system for studies of sexual dimorphism and genome evolution. Here, we report the chromosome‐level assembly of the O.macrolepis genome obtained from the integration of nanopore long‐read sequencing with physical maps produced using Bionano and Hi‐C technology. A total of 87.9 Gb of nanopore sequence provided approximately 100‐fold coverage of the genome. The preliminary genome assembly was 883.2 Mb in size with a contig N50 size of 11.2 Mb. The 969 corrected contigs obtained from Bionano optical mapping were assembled into 853 scaffolds and produced an assembly of 886.5 Mb with a scaffold N50 of 16.5 Mb. Finally, using the Hi‐C data, 881.3 Mb (99.4% of genome) in 526 scaffolds were anchored and oriented in 25 chromosomes ranging in size from 25.27 to 56.49 Mb. In total, 24,770 protein‐coding genes were predicted in the genome, and ~96.85% of the genes were functionally annotated. The annotated assembly contains 93.3% complete genes from the BUSCO reference set. In addition, we identified 409 Mb (46.23% of the genome) of repetitive sequence, and 11,213 non‐coding RNAs, in the genome. Evolutionary analysis revealed that O. macrolepis diverged from common carp approximately 24.25 million years ago. The chromosomes of O. macrolepis showed an unambiguous correspondence to the chromosomes of zebrafish. The high‐quality genome assembled in this work provides a valuable genomic resource for further biological and evolutionary studies of O. macrolepis.     

A chromosome‐level genome assembly of the woolly apple aphid,Eriosoma lanigerum Hausmann (Hemiptera: Aphididae)

Woolly apple aphid (WAA, Eriosoma lanigerum Hausmann) (Hemiptera: Aphididae) is a major pest of apple trees (Malus domestica, order Rosales) and is critical to the economics of the apple industry in most parts of the world. Here, we generated a chromosome‐level genome assembly of WAA—representing the first genome sequence from the aphid subfamily Eriosomatinae—using a combination of 10X Genomics linked‐reads and in vivo Hi‐C data. The final genome assembly is 327 Mb, with 91% of the assembled sequences anchored into six chromosomes. The contig and scaffold N50 values are 158 kb and 71 Mb, respectively, and we predicted a total of 28,186 protein‐coding genes. The assembly is highly complete, including 97% of conserved arthropod single‐copy orthologues based on Benchmarking Universal Single‐Copy Orthologs (busco ) analysis. Phylogenomic analysis of WAA and nine previously published aphid genomes, spanning four aphid tribes and three subfamilies, reveals that the tribe Eriosomatini (represented by WAA) is recovered as a sister group to Aphidini + Macrosiphini (subfamily Aphidinae). We identified syntenic blocks of genes between our WAA assembly and the genomes of other aphid species and find that two WAA chromosomes (El5 and El6) map to the conserved Macrosiphini and Aphidini X chromosome. Our high‐quality WAA genome assembly and annotation provides a valuable resource for research in a broad range of areas such as comparative and population genomics, insect–plant interactions and pest resistance management.     

Chromosome‐level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high‐quality chromosome‐level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single‐molecule real‐time sequencing, Illumina paired‐end sequencing, 10X Genomics linked‐reads and Hi‐C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single‐Copy Orthologs were determined in this genome assembly. And Hi‐C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein‐coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.     

A chromosome‐level genome assembly of the giant grouper (Epinephelus lanceolatus) provides insights into its innate immunity and rapid growth

The giant grouper (Epinephelus lanceolatus) is the largest coral reef teleost, with a native range that spans temperate and tropical waters in the Pacific and the Indian Oceans. It is cultured artificially and used as a breeding species in aquaculture due to its rapid growth rate. Here we report a giant grouper genome assembled at the chromosome scale from sequences generated using Illumina and high‐throughput chromatin conformation capture (Hi‐C) technology. The assembly comprised 1.086 Gb, with 98.4% of the scaffold sequences anchored into 24 chromosomes. The contig and scaffold N50 values were 119.9 kb and 46.2 Mb, respectively. The assembly is of high integrity, including 96.4% universal single‐copy orthologues based on BUSCO analysis. Through chromosome‐scale evolution analysis, we identified alignments of six giant grouper chromosomes to three stickleback chromosomes and some of the genes located within the breakpoints of reshuffling events may related to development and growth. From the 24,718 protein‐coding genes, we found that several gene families related to innate immunity and glycan biosynthesis were significantly expanded in the giant grouper genome compared to other teleost genomes. In addition, we identified several genes related to the hormone signalling pathway and innate immunity that have experienced positive selection or accelerated evolution, implicating their roles in immune defence and fast growth of the species. The high‐quality genome assembly will provide a valuable genomic resource for further biological and evolutionary studies, and useful genomic tools for breeding of the giant grouper.     

Finding Nemo’s Genes: A chromosome‐scale reference assembly of the genome of the orange clownfish Amphiprion percula

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome‐scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single‐molecule real‐time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi‐C‐based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein‐coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.     

High quality genome of Erigeron breviscapus provides a reference for herbal plants in Asteraceae

Erigeron breviscapus is an important medicinal plant in Compositae and the first species to realize the whole process from the decoding of the draft genome sequence to scutellarin biosynthesis in yeast. However, the previous low‐quality genome assembly has hindered the optimization of candidate genes involved in scutellarin synthesis and the development of molecular‐assisted breeding based on the genome. Here, the E. breviscapus genome was updated using PacBio RSII sequencing data and Hi‐C data, and increased in size from 1.2 Gb to 1.43 Gb, with a scaffold N50 of 156.82 Mb and contig N50 of 140.95 kb, and a total of 43,514 protein‐coding genes were obtained and oriented onto nine pseudo‐chromosomes, thus becoming the third plant species assembled to chromosome level after sunflower and lettuce in Compositae. Fourteen genes with evidence for positive selection were identified and found to be related to leaf morphology, flowering and secondary metabolism. The number of genes in some gene families involved in flavonoid biosynthesis in E. breviscapus have been significantly expanded. In particular, additional candidate genes involved in scutellarin biosynthesis, such as flavonoid‐7‐O‐glucuronosyltransferase genes (F7GATs) were identified using updated genome. In addition, three candidate genes encoding indole‐3‐pyruvate monooxygenase YUCCA2 (YUC2), serine carboxypeptidase‐like 18 (SCPL18), and F‐box protein (FBP), respectively, were identified to be probably related to leaf development and flowering by resequencing 99 individuals. These results provided a substantial genetic basis for improving agronomic and quality traits of E. breviscapus, and provided a platform for improving other draft genome assemblies to chromosome‐level.     

De novo assembly of a chromosome‐level reference genome of red‐spotted grouper (Epinephelus akaara) using nanopore sequencing and Hi‐C

The red‐spotted grouper Epinephelus akaara (E. akaara) is one of the most economically important marine fish in China, Japan and South‐East Asia and is a threatened species. The species is also considered a good model for studies of sex inversion, development, genetic diversity and immunity. Despite its importance, molecular resources for E. akaara remain limited and no reference genome has been published to date. In this study, we constructed a chromosome‐level reference genome of E. akaara by taking advantage of long‐read single‐molecule sequencing and de novo assembly by Oxford Nanopore Technology (ONT) and Hi‐C. A red‐spotted grouper genome of 1.135 Gb was assembled from a total of 106.29 Gb polished Nanopore sequence (GridION, ONT), equivalent to 96‐fold genome coverage. The assembled genome represents 96.8% completeness (BUSCO) with a contig N50 length of 5.25 Mb and a longest contig of 25.75 Mb. The contigs were clustered and ordered onto 24 pseudochromosomes covering approximately 95.55% of the genome assembly with Hi‐C data, with a scaffold N50 length of 46.03 Mb. The genome contained 43.02% repeat sequences and 5,480 noncoding RNAs. Furthermore, combined with several RNA‐seq data sets, 23,808 (99.5%) genes were functionally annotated from a total of 23,923 predicted protein‐coding sequences. The high‐quality chromosome‐level reference genome of E. akaara was assembled for the first time and will be a valuable resource for molecular breeding and functional genomics studies of red‐spotted grouper in the future.     

Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.     

Chromosome‐level reference genome of X12, a highly virulent race of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean that is spreading across major soybean production regions worldwide. Increased SCN virulence has recently been observed in both the United States and China. However, no study has reported a genome assembly for H. glycines at the chromosome scale. Herein, the first chromosome‐level reference genome of X12, an unusual SCN race with high infection ability, is presented. Using whole‐genome shotgun (WGS) sequencing, Pacific Biosciences (PacBio) sequencing, Illumina paired‐end sequencing, 10X Genomics linked reads and high‐throughput chromatin conformation capture (Hi‐C) genome scaffolding techniques, a 141.01‐megabase (Mb) assembled genome was obtained with scaffold and contig N50 sizes of 16.27 Mb and 330.54 kilobases (kb), respectively. The assembly showed high integrity and quality, with over 90% of Illumina reads mapped to the genome. The assembly quality was evaluated using Core Eukaryotic Genes Mapping Approach and Benchmarking Universal Single‐Copy Orthologs. A total of 11,882 genes were predicted using de novo, homolog and RNAseq data generated from eggs, second‐stage juveniles (J2), third‐stage juveniles (J3) and fourth‐stage juveniles (J4) of X12, and 79.0% of homologous sequences were annotated in the genome. These high‐quality X12 genome data will provide valuable resources for research in a broad range of areas, including fundamental nematode biology, SCN–plant interactions and co‐evolution, and also contribute to the development of technology for overall SCN management.     

Chromosome‐level genome assembly of the predator Propylea japonica to understand its tolerance to insecticides and high temperatures

The ladybird beetle Propylea japonica is an important natural enemy in agro‐ecological systems. Studies on the strong tolerance of P. japonica to high temperatures and insecticides, and its population and phenotype diversity have recently increased. However, abundant genome resources for obtaining insights into stress‐resistance mechanisms and genetic intra‐species diversity for P. japonica are lacking. Here, we constructed the P. japonica genome maps using Pacific Bioscience (PacBio) and Illumina sequencing technologies. The genome size was 850.90 Mb with a contig N50 of 813.13 kb. The Hi‐C sequence data were used to upgrade draft genome assemblies; 4,777 contigs were assembled to 10 chromosomes; and the final draft genome assembly was 803.93 Mb with a contig N50 of 813.98 kb and a scaffold N50 of 100.34 Mb. Approximately 495.38 Mb of repeated sequences was annotated. The 18,018 protein‐coding genes were predicted, of which 95.78% were functionally annotated, and 1,407 genes were species‐specific. The phylogenetic analysis showed that P. japonica diverged from the ancestor of Anoplophora glabripennis and Tribolium castaneum ~ 236.21 million years ago. We detected that some important gene families involved in detoxification of pesticides and tolerance to heat stress were expanded in P. japonica, especially cytochrome P450 and Hsp70 genes. Overall, the high‐quality draft genome sequence of P. japonica will provide invaluable resource for understanding the molecular mechanisms of stress resistance and will facilitate the research on population genetics, evolution and phylogeny of Coccinellidae. This genome will also provide new avenues for conserving the diversity of predator insects.     

Genome sequence of M6, a diploid inbred clone of the high‐glycoalkaloid‐producing tuber‐bearing potato species Solanum chacoense,reveals residual heterozygosity

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self‐compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high‐confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome‐enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.     

Chromosome‐level genome assembly of the greenhouse whitefly (Trialeurodes vaporariorum Westwood)

The greenhouse whitefly, Trialeurodes vaporariorum Westwood, is an agricultural pest of global importance. Here we report a 787‐Mb high‐quality draft genome sequence of T. vaporariorum assembled from PacBio long reads and Hi‐C chromatin interaction maps, which has scaffold and contig N50 lengths of 70 Mb and 500 kb, respectively, and contains 18,275 protein‐coding genes. About 98.8% of the assembled contigs were placed onto the 11 T. vaporariorum chromosomes. Comparative genomic analysis reveals significantly expanded gene families such as aspartyl proteases in T. vaporariorum compared to Bemisia tabaci Mediterranean (MED) and Middle East‐Asia Minor 1 (MEAM1). Furthermore, the cytochrome CYP6 subfamily shows significant expansion in T. vaporariorum and several genes in this subfamily display developmental stage‐specific expression patterns. The high‐quality T. vaporariorum genome provides a valuable resource for research in a broad range of areas such as fundamental molecular ecology, insect–plant/insect–microorganism or virus interactions and pest resistance management.     

Chromosome‐scale pseudomolecules refined by optical,physical and genetic maps in flax

Genomes of varying sizes have been sequenced with next‐generation sequencing platforms. However, most reference sequences include draft unordered scaffolds containing chimeras caused by mis‐scaffolding. A BioNano genome (BNG) optical map was constructed to improve the previously sequenced flax genome (Linum usitatissimum L., 2n = 30, about 373 Mb), which consisted of 3852 scaffolds larger than 1 kb and totalling 300.6 Mb. The high‐resolution BNG map of cv. CDC Bethune totalled 317 Mb and consisted of 251 BNG contigs with an N50 of 2.15 Mb. A total of 622 scaffolds (286.6 Mb, 94.9%) aligned to 211 BNG contigs (298.6 Mb, 94.2%). Of those, 99 scaffolds, diagnosed to contain assembly errors, were refined into 225 new scaffolds. Using the newly refined scaffold sequences and the validated bacterial artificial chromosome‐based physical map of CDC Bethune, the 211 BNG contigs were scaffolded into 94 super‐BNG contigs (N50 of 6.64 Mb) that were further assigned to the 15 flax chromosomes using the genetic map. The pseudomolecules total about 316 Mb, with individual chromosomes of 15.6 to 29.4 Mb, and cover 97% of the annotated genes. Evidence from the chromosome‐scale pseudomolecules suggests that flax has undergone palaeopolyploidization and mesopolyploidization events, followed by rearrangements and deletions or fusion of chromosome arms from an ancient progenitor with a haploid chromosome number of eight.     

A chromosome‐scale genome assembly of Antheraea pernyi (Saturniidae,Lepidoptera)

Antheraea pernyi is a semi‐domesticated lepidopteran insect species valuable to the silk industry, human health, and ecological tourism. Owing to its economic influence and developmental properties, it serves as an ideal model for investigating divergence of the Bombycoidea super family. However, studies on the karyotype evolution and functional genomics of A. pernyi are limited by scarce genomic resource. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first high‐quality A. pernyi genome from a single male individual. The genome is 720.67 Mb long with 49 chromosomes and a 13.77‐Mb scaffold N50. Approximately 441.75 Mb, accounting for 60.74% of the genome, was identified as repeats. The genome comprises 21,431 protein‐coding genes, 85.22% of which were functionally annotated. Comparative genomics analysis suggested that A. pernyi diverged from its common ancestor with A. yamamai ~30.3 million years ago, and that chromosome fission contributed to the increased chromosome number. The genome assembled in this work will not only facilitate future research on A. pernyi and related species but also help to progress comparative genomics analyses in Lepidoptera.     

Apolygus lucorum genome provides insights into omnivorousness and mesophyll feeding

Apolygus lucorum (Miridae) is an omnivorous pest that occurs worldwide and is notorious for the serious damage it causes to various crops and substantial economic losses. Although some studies have examined the biological characteristics of the mirid bug, no reference genome is available in Miridae, limiting in‐depth studies of this pest. Here, we present a chromosome‐scale reference genome of A. lucorum, the first sequenced Miridae species. The assembled genome size was 1.02 Gb with a contig N50 of 785 kb. With Hi‐C scaffolding, 1,016 Mb contig sequences were clustered, ordered and assembled into 17 large scaffolds with scaffold N50 length 68 Mb, each corresponding to a natural chromosome. Numerous transposable elements occur in this genome and contribute to the large genome size. Expansions of genes associated with omnivorousness and mesophyll feeding such as those related to digestion, chemosensory perception, and detoxification were observed in A. lucorum, suggesting that gene expansion contributed to its strong environmental adaptability and severe harm to crops. We clarified that a salivary enzyme polygalacturonase is unique in mirid bugs and has significantly expanded in A. lucorum, which may contribute to leaf damage from this pest. The reference genome of A. lucorum not only facilitates biological studies of Hemiptera as well as an understanding of the damage mechanism of mesophyll feeding, but also provides a basis on which to develop efficient control technologies for mirid bugs.     

Chromosome‐level genome assembly of the greenfin horse‐faced filefish (Thamnaconus septentrionalis) using Oxford Nanopore PromethION sequencing and Hi‐C technology

The greenfin horse‐faced filefish, Thamnaconus septentrionalis, is a valuable commercial fish species that is widely distributed in the Indo‐West Pacific Ocean. This fish has characteristic blue–green fins, rough skin and a spine‐like first dorsal fin. Thamnaconus septentrionalis is of conservation concern because its population has declined sharply, and it is an important marine aquaculture fish species in China. Genomic resources for the filefish are lacking, and no reference genome has been released. In this study, the first chromosome‐level genome of T. septentrionalis was constructed using nanopore sequencing and Hi‐C technology. A total of 50.95 Gb polished nanopore sequences were generated and were assembled into a 474.31‐Mb genome, accounting for 96.45% of the estimated genome size of this filefish. The assembled genome contained only 242 contigs, and the achieved contig N50 was 22.46 Mb, a surprisingly high value among all sequenced fish species. Hi‐C scaffolding of the genome resulted in 20 pseudochromosomes containing 99.44% of the total assembled sequences. The genome contained 67.35 Mb of repeat sequences, accounting for 14.2% of the assembly. A total of 22,067 protein‐coding genes were predicted, 94.82% of which were successfully annotated with putative functions. Furthermore, a phylogenetic tree was constructed using 1,872 single‐copy orthologous genes, and 67 unique gene families were identified in the filefish genome. This high‐quality assembled genome will be a valuable resource for a range of future genomic, conservation and breeding studies of T. septentrionalis.