Highlights d Six high-quality ixodid tick genomes and 678 re-sequenced tick specimens d Insights into the genetic basis of tick hematophagy and related phenotypes d Population structure and genetic diversity of six tick species d Tick-borne pathogen composition and distribution by metagenome analyses

Highlights d Chromosome-level assembly and methylome of the largest gymnosperm genome so far d Continuous expansion and slow removal of transposons cause conifer huge genome d Large genes with ultra-long introns tend to be expressed at higher levels d Distinctive reproductive evolutionary trajectory compared to angiosperms

The fall armyworm (Spodoptera frugiperda) is a lepidopteran insect pest that causes huge economic losses. This notorious insect pest has rapidly spread over the world in the past few years. However, the mechanisms of rapid dispersal are not well understood. Here, we report a chromosome-level assembled genome of the fall armyworm, named the ZJ-version, using PacBio and Hi-C technology. The sequenced individual was a female collected from the Zhejiang province of China and had high heterozygosity. The assembled genome size of ZJ-version was 486 Mb, containing 361 contigs with an N50 of 1.13 Mb. Hi-C scaffolding further assembled the genome into 31 chromosomes and a portion of W chromosome, representing 97.4% of all contigs and resulted in a chromosome-level genome with scaffold N50 of 16.3 Mb. The sex chromosomes were identified by genome resequencing of a single male pupa and a single female pupa. About 28% of the genome was annotated as repeat sequences, and 22,623 protein-coding genes were identified. Comparative genomics revealed the expansion of the detoxification-associated gene families, chemoreception-associated gene families, nutrition metabolism and transport system gene families in the fall armyworm. Transcriptomic and phylogenetic analyses focused on these gene families revealed the potential roles of the genes in polyphagia and invasion of fall armyworm. The high-quality of the fall armyworm genome provides an important genomic resource for further explorations of the mechanisms of polyphagia and insecticide resistance, as well as for pest management of fall armyworm. K E Y W O R D S chromosome-level genome, comparative genomics, fall armyworm, insecticide resistance, polyphagia | 1051 XIAO et Al.

Cannabis sativa is a well-known plant species that has great economic and ecological significance. An incomplete genome of cloned C. sativa was obtained by using SOAPdenovo software in 2011. To further explore the utilization of this plant resource, we generated an updated draft genome sequence for wild-type varieties of C. sativa in China using PacBio single-molecule sequencing and Hi-C technology. Our assembled genome is approximately 808 Mb, with scaffold and contig N50 sizes of 83.00 Mb and 513.57 kb, respectively. Repetitive elements account for 74.75% of the genome. A total of 38,828 protein-coding genes were annotated, 98.20% of which were functionally annotated. We provide the first comprehensive de novo genome of wild-type varieties of C. sativa distributed in Tibet, China. Due to long-term growth in the wild environment, these varieties exhibit higher heterozygosity and contain more genetic information. This genetic resource is of great value for future investigations of cannabinoid metabolic pathways and will aid in promoting the commercial production of C. sativa and the effective utilization of cannabinoids. The assembled genome is also a valuable resource for intensively and effectively investigating the C. sativa genome further in the future.

Panax notoginseng , a perennial herb of the genus Panax in the family Araliaceae, has played an important role in clinical treatment in China for thousands of years because of its extensive pharmacological effects. Here, we report a high-quality reference genome of P. notoginseng , with a genome size up to 2.66 Gb and a contig N50 of 1.12 Mb, produced with third-generation PacBio sequencing technology. This is the first chromosome-level genome assembly for the genus Panax . Through genome evolution analysis, we explored phylogenetic and whole-genome duplication events and examined their impact on saponin biosynthesis. We performed a detailed transcriptional analysis of P. notoginseng and explored gene-level mechanisms that regulate the formation of characteristic tubercles. Next, we studied the biosynthesis and regulation of saponins at temporal and spatial levels. We combined multi-omics data to identify genes that encode key enzymes in the P. notoginseng terpenoid biosynthetic pathway. Finally, we identified five glycosyltransferase genes whose products catalyzed the formation of different ginsenosides in P. notoginseng . The genetic information obtained in this study provides a resource for further exploration of the growth characteristics, cultivation, breeding, and saponin biosynthesis of P. notoginseng .

The fruits of Physalis (Solanaceae) have a unique structure, a lantern-like fruiting calyx known as inflated calyx syndrome (ICS) or the Chinese lantern, and are rich in steroid-related compounds. However, the genetic variations underlying the origin of these characteristic traits and diversity in Physalis remain largely unknown. Here, we present a high-quality chromosome-level reference genome assembly of Physalis floridana (~1.40 Gb in size) with a contig N50 of ~4.87 Mb. Through evolutionary genomics and experimental approaches, we found that the loss of the SEP-like MADS-box gene MBP21 subclade is likely a key mutation that, together with the previously revealed mutation affecting floral MPF2 expression, might have contributed to the origination of ICS in Physaleae, suggesting that the origination of a morphological novelty may have resulted from an evolutionary scenario in which one mutation compensated for another deleterious mutation. Moreover, the significant expansion of squalene epoxidase genes is potentially associated with the natural variation of steroid-related compounds in Physalis fruits. The results reveal the importance of gene gains (duplication) and/or subsequent losses as genetic bases of the evolution of distinct fruit traits, and the data serve as a valuable resource for the evolutionary genetics and breeding of solanaceous crops.

30Rhizobium nitrogen-fixing nodules are a well-known trait of legumes, but nodules also occur 31 in other plant lineages either with rhizobium or the actinomycete Frankia as microsymbiont. 32The widely accepted hypothesis is that nodulation evolved independently multiple times, with 33 only a few losses. However, insight in the evolutionary trajectory of nodulation is lacking. We 34 conducted comparative studies using Parasponia (Cannabaceae), the only non-legume able This finding challenges a long-standing hypothesis on evolution of nitrogen-fixing symbioses, 42and has profound implications for translational approaches aimed at engineering nitrogen-43 fixing nodules in crop plants.

The importance of structural variants (SVs) on phenotypes and human diseases is now recognized. Although a variety of SV detection platforms and strategies that vary in sensitivity and specificity have been developed, few benchmarking procedures are available to confidently assess their performances in biological and clinical research. To facilitate the validation and application of those approaches, our work established an Asian reference material comprising identified benchmark regions and high-confidence SV calls. We established a high-confidence SV callset with 8,938 SVs in an EBV immortalized B lymphocyte line, by integrating four alignment-based SV callers [from 109× PacBio continuous long read (CLR), 22× PacBio circular consensus sequencing (CCS) reads, 104× Oxford Nanopore long reads, and 114× optical mapping platform (Bionano)] and one de novo assembly-based SV caller using CCS reads. A total of 544 randomly selected SVs were validated by PCR and Sanger sequencing, proofing the robustness of our SV calls. Combining trio-binning based haplotype assemblies, we established an SV benchmark for identification of false negatives and false positives by constructing the continuous high confident regions (CHCRs), which cover 1.46Gb and 6,882 SVs supported by at least one diploid haplotype assembly. Establishing high-confidence SV calls for a benchmark sample that has been characterized by multiple technologies provides a valuable resource for investigating SVs in human biology, disease, and clinical diagnosis.