GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes

Assessing de novo parasite genomes assembled using only Oxford Nanopore Technologies MinION data

In this study, we assessed the quality of de novo genome assemblies for three species of parasitic nematodes (Brugia malayi, Trichuris trichiura, and Ancylostoma caninum) generated using only Oxford Nanopore Technologies MinION data. Assemblies were compared to current reference genomes and against additional assemblies that were supplemented with short-read Illumina data through polishing or hybrid assembly approaches. For each species, assemblies generated using only MinION data had similar or superior measures of contiguity, completeness, and gene content. In terms of gene composition, depending on the species, between 88.9 and 97.6% of complete coding sequences predicted in MinION data only assemblies were identical to those predicted in assemblies polished with Illumina data. Polishing MinION data only assemblies with Illumina data therefore improved gene-level accuracy to a degree. Furthermore, modified DNA extraction and library preparation protocols produced sufficient genomic DNA from B. malayi and T. trichiura to generate de novo assemblies from individual specimens.

Genomic resources of aquatic Lepidoptera, Elophila obliteralis and Hyposmocoma kahamanoa, reveal similarities with Trichoptera in amino acid composition of major silk genes

While most species of butterflies and moths (Lepidoptera) have entirely terrestrial life histories, ∼0.5% of the described species are known to have an aquatic larval stage. Larvae of aquatic Lepidoptera are similar to caddisflies (Trichoptera) in that they use silk to anchor themselves to underwater substrates or to build protective cases. However, the physical properties and genetic elements of silks in aquatic Lepidoptera remain unstudied, as most research on lepidopteran silk has focused on the commercially important silkworm, Bombyx mori. Here, we provide high-quality PacBio HiFi genome assemblies of 2 distantly-related aquatic Lepidoptera species [Elophila obliteralis (Pyraloidea: Crambidae) and Hyposmocoma kahamanoa (Gelechioidea: Cosmopterigidae)]. As a step toward understanding the evolution of underwater silk in aquatic Lepidoptera, we used the genome assemblies and compared them to published genetic data of aquatic and terrestrial Lepidoptera. Sequences of the primary silk protein, h-fibroin, in aquatic moths have conserved termini and share a basic motif structure with terrestrial Lepidoptera. However, these sequences were similar to aquatic Trichoptera in that the percentage of positively and negatively charged amino acids was much higher than in terrestrial Lepidoptera, indicating a possible adaptation of silks to aquatic environments.

Chromosome-level genome assembly of predatory Arma chinensis

Arma chinensis is a natural enemy that preys on various species and can suppress agricultural and forest pests in the orders Lepidoptera and Coleoptera. Here, we aimed to determine the genome of A. chinensis assembled at the chromosome-level using PacBio and Hi-C technologies. The assembled genome was 986 Mb, with a contig N50 of 2.40 Mb, scaffold N50 of 134.98 Mb, and BUSCO completeness of 96.10%. Hi-C data aided in anchoring the assembly onto seven chromosomes. A sequence of ~ 496.2 Mb was annotated as a repeat element, constituting 51.15% of the genome. We functionally annotated 84.79% of 20,853 predicted protein-encoding genes. This high-quality A. chinensis genome provides a novel genomic resource for future research on Pentatomidae insects.

Chromosome-scale Elaeis guineensis and E. oleifera assemblies: comparative genomics of oil palm and other Arecaceae

Elaeis guineensis and E. oleifera are the two species of oil palm. E. guineensis is the most widely cultivated commercial species, and introgression of desirable traits from E. oleifera is ongoing. We report an improved E. guineensis genome assembly with substantially increased continuity and completeness, as well as the first chromosome-scale E. oleifera genome assembly. Each assembly was obtained by integration of long-read sequencing, proximity ligation sequencing, optical mapping, and genetic mapping. High interspecific genome conservation is observed between the two species. The study provides the most extensive gene annotation to date, including 46,697 E. guineensis and 38,658 E. oleifera gene predictions. Analyses of repetitive element families further resolve the DNA repeat architecture of both genomes. Comparative genomic analyses identified experimentally validated small structural variants between the oil palm species and resolved the mechanism of chromosomal fusions responsible for the evolutionary descending dysploidy from 18 to 16 chromosomes.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

Comparative genomics provides insights into molecular adaptation to hypermetamorphosis and cantharidin metabolism in blister beetles (Coleoptera: Meloidae)

Blister beetles (Coleoptera: Meloidae) are currently subdivided into three subfamilies: Eleticinae (a basal group), Nemognathinae, and Meloinae. These are all characterized by the endogenous production of the defensive terpene cantharidin (CA), whereas the two most derived subfamilies show a hypermetamorphic larval development. Here, we provide novel draft genome assemblies of five species sampled across the three blister beetle subfamilies (Iselma pallidipennis, Stenodera caucasica, Zonitis immaculata, Lydus trimaculatus, and Mylabris variabilis) and performed a comparative analysis with other available Meloidae genomes and the closely-related canthariphilous species (Pyrochroa serraticornis) to disclose adaptations at a molecular level. Our results highlighted the expansion and selection of genes potentially responsible for CA production and metabolism, as well as its mobilization and vesicular compartmentalization. Furthermore, we observed adaptive selection patterns and gain of genes devoted to epigenetic regulation, development, and morphogenesis, possibly related to hypermetamorphosis. We hypothesize that most genetic adaptations occurred to support both CA biosynthesis and hypermetamorphosis, two crucial aspects of Meloidae biology that likely contributed to their evolutionary success.

Highly contiguous genome assembly and gene annotation of the short-finned eel (Anguilla bicolor pacifica)

In East Asia, anguillid eels are commercially important. However, unlike other species, they have not been successfully cultivated throughout their lifecycle. Facing population decline due to overharvesting and environmental pressures, the industry is turning to alternatives, such as Anguilla bicolor pacifica (short-finned eel). However, genomic data for short-finned eels are unavailable. Here, we present in-depth whole-genome sequencing results for short-finned eel obtained using two sequencing platforms (PacBio Revio, and Illumina). In this study, we achieved a highly contiguous genome assembly of the short-finned eel, comprising 19 pseudochromosomes encompassing 99.76% of the 1.087 Gb genome sequence with an N50 of 16.88 and 61.07 Mb from contig and scaffold, respectively. Transcripts from four different tissues led to the annotation of 23,095 protein-coding genes in the eel genome, 98.66% of which were functionally annotated. This high-quality genome assembly, along with the annotation data, provides a foundation for future functional genomic studies of short-finned eels.

A high-quality genome assembly of the shrubby cinquefoil (Dasiphora fruticosa)

Dasiphora fruticosa (Rosaceae), commonly known as shrubby cinquefoil, is a flowering shrub of high ornamental value yet underutilized in East Asian landscapes. Given its broad elevational distribution range, D. fruticosa serves as an ideal model for studying genetic adaptations and speciation along elevation gradients. Here, we present a high-quality chromosome-scale assembly of D. fruticosa with a genome size of 249.23 Mb and a contig N50 length of 14.01 Mb. The genome sequence contains 32,613 protein-coding genes, of which 30,643 (93.96%) were functionally annotated. Compared to the published D. fruticosa genome sequence, our assembly demonstrates higher completeness and continuity. Furthermore, comparative genomic analyses provide insights into the phylogenetic relationship and high-altitude adaptation of D. fruticosa. Overall, our study offers a valuable genetic resource for both molecular and evolutionary research on shrubby cinquefoil.

First genome assembly of the order Strepsiptera using PacBio HiFi reads reveals a miniature genome

Twisted-wing insects (Strepsiptera) are an enigmatic order of parasites with unusual life histories and striking sexual dimorphism. Males emerge from hosts as free-living winged adults, while females from most species remain as endoparasites that retain larval traits. Due to scarce genomic data and phylogenetic controversies, Strepsiptera was only recently placed as the closest living relative to beetles (Coleoptera). Here, we report the first PacBio HiFi genome assembly of the strepsipteran Xenos peckii (Xenidae). This de novo assembly size is 72.1 Mb, with a BUSCO score of 87.4%, N50 of 7.3 Mb, 23.4% GC content, and 38.41% repeat content. We identified 8 contigs that contain >75% of the assembly and reflect the haploid chromosome number reported from karyotypic data, and 3 contigs that exhibit sex chromosome coverage patterns. Additionally, the mitochondrial genome is 16,111 bp long and has 37 genes. This long-read assembly for Strepsiptera reveals a miniature genome and provides a unique tool to understand complex genome evolution associated with a parasitic lifestyle and extreme sexual dimorphism.

The chromosome level genome assembly of the Asian green mussel, Perna viridis

The Asian green mussel, Perna viridis is an important aquaculture species in the family Mytilidae contributing substantially to molluscan aquaculture. We generated a high-quality chromosome level assembly of this species by combining PacBio single molecule sequencing technique (SMRT), Illumina paired-end sequencing, high-throughput chromosome conformation capture technique (Hi-C) and Bionano mapping. The final assembly resulted in a genome of 723.49 Mb in size with a scaffold N50 of 49.74 Mb with 99% anchored into 15 chromosomes. A total of 49654 protein-coding genes were predicted from the genome. The presence of 634 genes associated with the cancer pathway and 408 genes associated with viral carcinogenesis indicates the potential of this species to be used as a model for cancer studies. The chromosome-level assembly of this species is also a valuable resource for further genomic selection and selective breeding for improving economically important aquaculture traits and augmenting aquaculture productivity.

Kinetochore and ionomic adaptation to whole-genome duplication in Cochlearia shows evolutionary convergence in three autopolyploids

Whole-genome duplication (WGD) occurs in all kingdoms and impacts speciation, domestication, and cancer outcome. However, doubled DNA management can be challenging for nascent polyploids. The study of within-species polyploidy (autopolyploidy) permits focus on this DNA management aspect, decoupling it from the confounding effects of hybridization (in allopolyploid hybrids). How is autopolyploidy tolerated, and how do young polyploids stabilize? Here, we introduce a powerful model to address this: the genus Cochlearia, which has experienced many polyploidization events. We assess meiosis and other polyploid-relevant phenotypes, generate a chromosome-scale genome, and sequence 113 individuals from 33 ploidy-contrasting populations. We detect an obvious autopolyploidy-associated selection signal at kinetochore components and ion transporters. Modeling the selected alleles, we detail evidence of the kinetochore complex mediating adaptation to polyploidy. We compare candidates in independent autopolyploids across three genera separated by 40 million years, highlighting a common function at the process and gene levels, indicating evolutionary flexibility in response to polyploidy.

Two outbreaks and sporadic occurrences of Candida auris from one hospital in China: an epidemiological, genomic retrospective study

OBJECTIVES

To investigate the clinical relevance, origin, transmission, and resistance of Candida auris (C. auris) isolates from two outbreaks and sporadic occurrences from one hospital in China.

METHODS

A total of 135 C. auris isolates were collected. Clinical characteristics were obtained and antifungal susceptibility testing (AFST) was performed using the method of broth microdilution. Phylogenetic tree, WGS analysis, and single nucleotide polymorphisms (SNPs) were used to determine the origin, transmission, and resistance mechanisms.

RESULTS

A total of 31 patients (91.2%, 31/34) received invasive medical procedures and 13 patients (38.2%, 13/34) had antifungal agents before C. auris infection/colonization, except one patient whose clinical information was missing. Only 4 cases of C. auris candidemia were observed. 18 patients died, 13 patients recovered, and the outcomes of 3 patients were not available. A total of 35 C. auris isolates, which were successfully cultivated and the first isolated or harbored specific drug-resistant phenotype from each patient, were selected to be sequenced and further analyzed. C. auris isolates presented low genetic variability and belonged to clade I, possibly originating from BJ004-H7 in Beijing. All 35 isolates were resistant to Fluconazole (FCZ) and amphotericin B (AMB), and 3 isolates were resistant to caspofungin (CAS). Mutations in ERG11 and FKS1 were linked to reduced azole and echinocandin susceptibility, respectively.

CONCLUSIONS

Two outbreaks of highly clonal, multidrug-resistant C. auris isolates within the medical facility were reported. The intensive performance of disinfection measures helped block in-hospital transmission. Understanding the epidemiology, drug resistance and management of C. auris will be helpful for implementing effective infection control and treatment strategies.

Improved pokeweed genome assembly and early gene expression changes in response to jasmonic acid

BACKGROUND

Jasmonic acid (JA) is a phytohormone involved in regulating responses to biotic and abiotic stress. Although the JA pathway is well characterized in model plants such as Arabidopsis thaliana, less is known about many non-model plants. Phytolacca americana (pokeweed) is native to eastern North Americana and is resilient to environmental stress. The goal of this study was to produce a publicly available pokeweed genome assembly and annotations and use this resource to determine how early response to JA changes gene expression, with particular focus on genes involved in defense.

RESULTS

We assembled the pokeweed genome de novo from approximately 30 Gb of PacBio Hifi long reads and achieved an NG50 of ~ 13.2 Mb and a minimum 93.9% complete BUSCO score for gene annotations. With this reference, we investigated the early changes in pokeweed gene expression following JA treatment. Approximately 5,100 genes were differentially expressed during the 0-6 h time course with almost equal number of genes with increased and decreased transcript levels. Cluster and gene ontology analyses indicated the downregulation of genes associated with photosynthesis and upregulation of genes involved in hormone signaling and defense. We identified orthologues of key transcription factors and constructed the first JA gene response network integrated with our transcriptomic data from orthologues of Arabidopsis genes. We discovered that pokeweed did not use leaf senescence as a means of reallocating resources during stress; rather, most secondary metabolite synthesis genes were constitutively expressed, suggesting that pokeweed directs its resources for survival over the long term. In addition, pokeweed synthesizes several RNA N-glycosylases hypothesized to function in defense, each with unique expression profiles in response to JA.

CONCLUSIONS

Our investigation of the early response of pokeweed to JA illustrates patterns of gene expression involved in defence and stress tolerance. Pokeweed provides insight into the defense mechanisms of plants beyond those observed in research models and crops, and further study may yield novel approaches to improving the resilience of plants to environmental changes. Our assembled pokeweed genome is the first within the taxonomic family Phytolaccaceae to be publicly available for continued research.

Pan-genome and phylogenomic analyses highlight Hevea species delineation and rubber trait evolution

The para rubber tree (Hevea brasiliensis) is the world's sole commercial source of natural rubber, a vital industrial raw material. However, the narrow genetic diversity of this crop poses challenges for rubber breeding. Here, we generate high-quality de novo genome assemblies for three H. brasiliensis cultivars, two H. brasiliensis wild accessions, and three other Hevea species (H. nitida, H. pauciflora, and H. benthamiana). Through analyzing genomes of 94 Hevea accessions, we identify five distinct lineages that do not align with their previous species delineations. We discover multiple accessions with hybrid origins between these lineages, indicating incomplete reproductive isolation between them. Only two out of four wild lineages have been introduced to commercial rubber cultivars. Furthermore, we reveal that the rubber production traits emerged following the development of a large REF/SRPP gene cluster and its functional specialization in rubber-producing laticifers within this genus. These findings would enhance rubber breeding and benefit research communities.

Hybrid assembly and comparative genomics unveil insights into the evolution and biology of the red-legged partridge

The red-legged partridge Alectoris rufa plays a crucial role in the ecosystem of southwestern Europe, and understanding its genetics is vital for conservation and management. Here we sequence, assemble, and annotate a highly contiguous and nearly complete version of its genome. This assembly encompasses 96.9% of the avian genes flagged as essential in the BUSCO aves_odb10 dataset. Moreover, we pinpointed RNA and protein-coding genes, 95% of which had functional annotations. Notably, we observed significant chromosome rearrangements in comparison to quail (Coturnix japonica) and chicken (Gallus gallus). In addition, a comparative phylogenetic analysis of these genomes suggests that A. rufa and C. japonica diverged roughly 20 million years ago and that their common ancestor diverged from G. gallus 35 million years ago. Our assembly represents a significant advancement towards a complete reference genome for A. rufa, facilitating comparative avian genomics, and providing a valuable resource for future research and conservation efforts for the red-legged partridge.

The genome assembly of Carex breviculmis provides evidence for its phylogenetic localization and environmental adaptation

BACKGROUND AND AIMS

Carex breviculmis is a perennial herb with good resistance and is widely used for forage production and turf management. It is important in ecology, environmental protection and biodiversity conservation, but faces several challenges due to human activities. However, the absence of genome sequences has limited basic research and the improvement of wild plants.

METHODS

We annotated the genome of C. breviculmis and conducted a systematic analysis to explore its resistance to harsh environments. We also conducted a comparative analysis of Achnatherum splendens, which is similarly tolerant to harsh environments.

KEY RESULTS

The assembled the genome comprises 469.01 Mb, revealing 37 372 genes with a BUSCO completeness score of 99.0 %. The genome has 52.03 % repetitive sequences, primarily influenced by recent LTR insertions that have contributed to its expansion. Phylogenetic analysis suggested that C. breviculmis diverged from C. littledalei ~6.61 million years ago. Investigation of repetitive sequences and expanded gene families highlighted a rapid expansion of tandem duplicate genes, particularly in areas related to sugar metabolism, synthesis of various amino acids, and phenylpropanoid biosynthesis. Additionally, our analysis identified crucial genes involved in secondary metabolic pathways, such as glycolysis, phenylpropanoid biosynthesis and amino acid metabolism, which have undergone positive selection. We reconstructed the sucrose metabolic pathway and identified significant gene expansions, including 16 invertase, 9 sucrose phosphate synthase and 12 sucrose synthase genes associated with sucrose metabolism, which showed varying levels of expansion.

CONCLUSIONS

The expansion of these genes, coupled with subsequent positive selection, contributed to the ability of C. breviculmis to adapt to environmental stressors. This study lays the foundation for future research on the evolution of Carex plants, their environmental adaptations, and potential genetic breeding.

A draft reference genome assembly of California Pipevine, Aristolochia californica Torr

The California Pipevine, Aristolochia californica Torr., is the only endemic California species within the cosmopolitan birthwort family Aristolochiaceae. It occurs as an understory vine in riparian and chaparral areas and in forest edges and windrows. The geographic range of this plant species almost entirely overlaps with that of its major specialized herbivore, the California Pipevine Swallowtail Butterfly Battus philenor hirsuta. While this species pair is a useful, ecologically well-understood system to study co-evolution, until recently, genomic resources for both have been lacking. Here, we report a new, chromosome-level assembly of A. californica as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 531 scaffolds spanning 661 megabase (Mb) pairs, with a contig N50 of 6.53 Mb, a scaffold N50 of 42.2 Mb, and BUSCO complete score of 98%. In combination with the recently published B. philenor hirsuta reference genome assembly, the A. californica reference genome assembly will be a powerful tool for studying co-evolution in a rapidly changing California landscape.

Chromosomal-level genome assembly of Hylurgus ligniperda: insights into host adaptation and environmental tolerance

BACKGROUND

Hylurgus ligniperda (Coleoptera: Curculionidae) is a worldwide forest quarantine pest. It is widely distributed, has many host tree species, and possesses strong adaptability. To explore its environmental adaptability and the related molecular mechanisms, we conducted chromosome-level genome sequencing and analyzed the transcriptome under different environmental factors, identifying key expressed genes.

RESULTS

We employed PacBio, Illumina, and Hi-C sequencing techniques to assemble a 520 Mb chromosomal-level genome of H. ligniperda, obtaining an N50 of 39.97 Mb across 138 scaffolds. A total of 10,765 protein-coding genes were annotated after repeat masking. Fourteen chromosomes were identified, among which Hyli14 was determined to be the sex chromosome. Survival statistics were tested over various growth periods under high temperature and low humidity conditions. The maximum survival period of adults reached 292 days at 25 °C, 65% relative humidity. In comparison, the maximum survival period was 14 days under 35 °C, 65% relative humidity, and 106 days under 25°C, 40% relative humidity. This indicated that environmental stress conditions significantly reduced adults' survival period. We further conducted transcriptome analysis to screen for potentially influential differentially expressed genes, such as CYP450 and Histone. Subsequently, we performed gene family analysis to gain insights into their functions and interactions, such as CYP450 and Histone. CYP450 genes affected the detoxification metabolism of enzymes in the Cytochrome P450 pathway to adapt to different environments. Histone genes are involved in insect hormone biosynthesis and longevity-regulating pathways in H. ligniperda to adapt to environmental stress.

CONCLUSIONS

The genome at the chromosome level of H. ligniperda was assembled for the first time. The mortality of H. ligniperda increased significantly at 35 ℃, 65% RH, and 25 ℃, 40% RH. CYP450 and Histone genes played an important role in response to environmental stress. This genome offers a substantial genetic resource for investigating the molecular mechanisms behind beetle invasion and spread.

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

A genome assembly of the American black bear, Ursus americanus, from California

The American black bear, Ursus americanus, is a widespread and ecologically important species in North America. In California, the black bear plays an important role in a variety of ecosystems and serves as an important species for recreational hunting. While research suggests that the populations in California are currently healthy, continued monitoring is critical, with genomic analyses providing an important surveillance tool. Here we report a high-quality, near chromosome-level genome assembly from a U. americanus sample from California. The primary assembly has a total length of 2.5 Gb contained in 316 scaffolds, a contig N50 of 58.9 Mb, a scaffold N50 of 67.6 Mb, and a BUSCO completeness score of 96%. This U. americanus genome assembly will provide an important resource for the targeted management of black bear populations in California, with the goal of achieving an appropriate balance between the recreational value of black bears and the maintenance of viable populations. The high quality of this genome assembly will also make it a valuable resource for comparative genomic analyses among black bear populations and among bear species.

Integrative genomics reveals the polygenic basis of seedlessness in grapevine

Seedlessness is a crucial quality trait in table grape (Vitis vinifera L.) breeding. However, the development of seeds involved intricate regulations, and the polygenic basis of seed abortion remains unclear. Here, we combine comparative genomics, population genetics, quantitative genetics, and integrative genomics to unravel the evolution and polygenic basis of seedlessness in grapes. We generated the haplotype-resolved genomes for two seedless grape cultivars, "Thompson Seedless" (TS, syn. "Sultania") and "Black Monukka" (BM). Comparative genomics identified a ∼4.25 Mb hemizygous inversion on Chr10 specific in seedless cultivars, with seedless-associated genes VvTT16 and VvSUS2 located at breakpoints. Population genomic analyses of 548 grapevine accessions revealed two distinct clusters of seedless cultivars, and the identity-by-descent (IBD) results indicated that the origin of the seedlessness trait could be traced back to "Sultania." Introgression, rather than convergent selection, shaped the evolutionary history of seedlessness in grape improvement. Genome-wide association study (GWAS) analysis identified 110 quantitative trait loci (QTLs) associated with 634 candidate genes, including previously unidentified candidate genes, such as three 11S GLOBULIN SEED STORAGE PROTEIN and two CYTOCHROME P450 genes, and well-known genes like VviAGL11. Integrative genomic analyses resulted in 339 core candidate genes categorized into 13 functional categories related to seed development. Machine learning-based genomic selection achieved a remarkable prediction accuracy of 97% for seedlessness in grapevines. Our findings highlight the polygenic nature of seedlessness and provide candidate genes for molecular genetics and an effective prediction for seedlessness in grape genomic breeding.

Chromosome-contiguous genome for the Haecon-5 strain of Haemonchus contortus reveals marked genetic variability and enables the discovery of essential gene candidates

Millions of livestock animals worldwide are infected with the haematophagous barber's pole worm, Haemonchus contortus, the aetiological agent of haemonchosis. Despite the major significance of this parasite worldwide and its widespread resistance to current treatments, the lack of a high-quality genome for the well-defined strain of this parasite from Australia, called Haecon-5, has constrained research in a number of areas including host-parasite interactions, drug discovery and population genetics. To enable research in these areas, we report here a chromosome-contiguous genome (∼280 Mb) for Haecon-5 with high-quality models for 19,234 protein-coding genes. Comparative genomic analyses show significant genomic similarity (synteny) with a UK strain of H. contortus, called MHco3(ISE).N1 (abbreviated as "ISE"), but we also discover marked differences in genomic structure/gene arrangements, distribution of nucleotide variability (single nucleotide polymorphisms (SNPs) and indels) and orthology between Haecon-5 and ISE. We used the genome and extensive transcriptomic resources for Haecon-5 to predict a subset of essential single-copy genes employing a "cross-species" machine learning (ML) approach using a range of features from nucleotide/protein sequences, protein orthology, subcellular localisation, single-cell RNA-seq and/or histone methylation data available for the model organisms Caenorhabditis elegans and Drosophila melanogaster. From a set of 1,464 conserved single copy genes, transcribed in key life-cycle stages of H. contortus, we identified 232 genes whose homologs have critical functions in C. elegans and/or D. melanogaster, and prioritised 10 of them for further characterisation; nine of the 10 genes likely play roles in neurophysiological processes, germline, hypodermis and/or respiration, and one is an unknown (orphan) gene for which no detailed functional information exists. Future studies of these genes/gene products are warranted to elucidate their roles in parasite biology, host-parasite interplay and/or disease. Clearly, the present Haecon-5 reference genome and associated resources now underpin a broad range of fundamental investigations of H. contortus and could assist in accelerating the discovery of novel intervention targets and drug candidates to combat haemonchosis.

Origin, structure, and composition of the spider major ampullate silk fiber revealed by genomics, proteomics, and single-cell and spatial transcriptomics

Spiders produce nature's toughest fiber using renewable components at ambient temperatures and with water as solvent, making it highly interesting to replicate for the materials industry. Despite this, much remains to be understood about the bioprocessing and composition of spider silk fibers. Here, we identify 18 proteins that make up the spiders' strongest silk type, the major ampullate fiber. Single-cell RNA sequencing and spatial transcriptomics revealed that the secretory epithelium of the gland harbors six cell types. These cell types are confined to three distinct glandular zones that produce specific combinations of silk proteins. Image analysis of histological sections showed that the secretions from the three zones do not mix, and proteomics analysis revealed that these secretions form layers in the final fiber. Using a multi-omics approach, we provide substantial advancements in the understanding of the structure and function of the major ampullate silk gland as well as of the architecture and composition of the fiber it produces.

Chromosome-level genome assembly and annotation of Flueggea virosa (Phyllanthaceae)

Flueggea virosa (Roxb. ex Willd.) Royle, an evergreen shrub and small tree in the Phyllanthaceae family, holds significant potential in garden landscaping and pharmacological applications. However, the lack of genomic data has hindered further scientific understanding of its horticultural and medicinal values. In this study, we have assembled a haplotype-resolved genome of F. virosa for the first time. The two haploid genomes, named haplotype A genome and haplotype B genome, are 487.33 Mb and 477.53 Mb in size, respectively, with contig N50 lengths of 31.45 Mb and 32.81 Mb. More than 99% of the assembled sequences were anchored to 13 pairs of pseudo-chromosomes. Furthermore, 21,587 and 21,533 protein-coding genes were predicted in haplotype A and haplotype B genomes, respectively. The availability of this chromosome-level genome fills the gap in genomic data for F. virosa and provides valuable resources for molecular studies of this species, supporting future research on speciation, functional genomics, and comparative genomics within the Phyllanthaceae family.

Chromosome-level and haplotype-resolved genome assembly of Dracaena cambodiana (Asparagaceae)

Dracaena cambodiana Pierre ex Gagn. (Asparagaceae) is the source plant of Dragon's blood and has high ornamental values in gardening. Currently, this species is classified as the second-class state-protected species in the National Key Protected Wild Plants (NKPWP) of China. However, limited genomic data has hindered a more comprehensive scientific understanding of the processes involved in the production of Dragon's blood and the related conservation genomics research. In this study, we assembled a haplotype-resolved genome of D. cambodiana. The haploid genomes, haplotype A and haplotype B, are 1,015.22 Mb and 1,003.13 Mb in size, respectively. The completeness of haplotype A and haplotype B genomes was 98.60% and 98.20%, respectively, using the "embryophyta_10" dataset. Haplotype A and haplotype B genomes contained 27,361 and 27,066 protein-coding genes, respectively, with nearly all being functionally annotated. These findings provide new insights into the genomic characteristics of D. cambodiana and will offer additional genomic resources for studying the biosynthesis mechanism of Dragon's blood and the horticultural application of Dragon trees.

Chromosome-Level Assembly Reveals a Fifteen-Chromosome Aneuploid Genome and Environmental Adaptation Strategy of Chinese Traditional Medical Fungus Wolfiporia hoelen

The sclerotia of Wolfiporia hoelen are one of the most important traditional Chinese medicines and foods commonly used in China, Japan, Korea, and other Asian countries. To provide a high-quality reference genome and deepen our understanding of the genome of W. hoelen to elucidate various biological phenomena. In this study, we assembled three genomes of W. hoelen using a combination of Nanopore and Illumina sequencing strategies. The fifteen-chromosome genome L7 of W. hoelen was assembled with two-sided telomere and rDNA sequences for the first time. The chromosome count was subsequently confirmed through collinearity analysis, correcting the previous belief that W. hoelen had only fourteen chromosomes. Moreover, the aneuploid genome was discovered in W. hoelen for the first time through sequencing depth analysis of different chromosomes, and only some strains of W. hoelen exhibit aneuploid genomes. According to the genome analysis of homokaryotic offspring and protoplast-isolated strains, a potential variation in chromosome allocation patterns was revealed. Moreover, the gene function enrichment analysis of genes on reduplicated chromosomes demonstrated that aneuploidy in the genome may be the result of environmental adaptation for W. hoelen. The discovery of an aneuploid genome also provides new ideas for genetic improvement of W. hoelen.

Allopolyploidization from two dioecious ancestors leads to recurrent evolution of sex chromosomes

Polyploidization presents an unusual challenge for species with sex chromosomes, as it can lead to complex combinations of sex chromosomes that disrupt reproductive development. This is particularly true for allopolyploidization between species with different sex chromosome systems. Here, we assemble haplotype-resolved chromosome-level genomes of a female allotetraploid weeping willow (Salix babylonica) and a male diploid S. dunnii. We show that weeping willow arose from crosses between a female ancestor from the Salix-clade, which has XY sex chromosomes on chromosome 7, and a male ancestor from the Vetrix-clade, which has ancestral XY sex chromosomes on chromosome 15. We find that weeping willow has one pair of sex chromosomes, ZW on chromosome 15, that derived from the ancestral XY sex chromosomes in the male ancestor of the Vetrix-clade. Moreover, the ancestral 7X chromosomes from the female ancestor of the Salix-clade have reverted to autosomal inheritance. Duplicated intact ARR17-like genes on the four homologous chromosomes 19 likely have contributed to the maintenance of dioecy during polyploidization and sex chromosome turnover. Taken together, our results suggest the rapid evolution and reversion of sex chromosomes following allopolyploidization in weeping willow.

Genomic and transcriptomic analyses of the elite rice variety Huizhan provide insight into disease resistance and heat tolerance

The indica rice variety Huizhan shows elite traits of disease resistance and heat tolerance. However, the underlying genetic basis of these traits is not fully understood due to limited genomic resources. Here, we used Nanopore long-read and next-generation sequencing technologies to generate a chromosome-scale genome assembly of Huizhan. Comparative genomics analysis uncovered a large chromosomal inversion and expanded gene families that are associated with plant growth, development and stress responses. Functional rice blast resistance genes, including Pi2, Pib and Ptr, and bacterial blight resistance gene Xa27, contribute to disease resistance of Huizhan. Furthermore, integrated genomics and transcriptomics analyses showed that OsHIRP1, OsbZIP60, the SOD gene family, and various transcription factors are involved in heat tolerance of Huizhan. The high-quality genome assembly and comparative genomics results presented in this study facilitate the use of Huizhan as an elite parental line in developing rice varieties adapted to disease pressure and climate challenges.

otb: an automated HiC/HiFi pipeline assembles the Prosapia bicincta Genome

The implementation of a new genomic assembly pipeline named only the best (otb) has effectively addressed various challenges associated with data management during the development and storage of genome assemblies. otb, which incorporates a comprehensive pipeline involving a setup layer, quality checks, templating, and the integration of Nextflow and Singularity. The primary objective of otb is to streamline the process of creating a HiFi/HiC genome, aiming to minimize the manual intervention required in the genome assembly process. The 2-lined spittlebug, (Prosapia bicincta, Hemiptera: Cercopidae), a true bug insect herbivore, serves as a practical test case for evaluating otb. The 2-lined spittlebug is both a crucial agricultural pest and a genomically understudied insect belonging to the order Hemiptera. This insect is a significant threat to grasslands and pastures, leading to plant wilting and phytotoxemia when infested. Its presence in tropical and subtropical regions around the world poses a long-term threat to the composition of plant communities in grassland landscapes, impacting rangelands, and posing a substantial risk to cattle production.

The reference genome of an endangered Asteraceae, Deinandra increscens subsp. villosa, endemic to the Central Coast of California

We present a reference genome for the federally endangered Gaviota tarplant, Deinandra increscens subsp. villosa (Madiinae, Asteraceae), an annual herb endemic to the Central California coast. Generating PacBio HiFi, Oxford Nanopore Technologies, and Dovetail Omni-C data, we assembled a haploid consensus genome of 1.67 Gb as 28.7 K scaffolds with a scaffold N50 of 74.9 Mb. We annotated repeat content in 74.8% of the genome. Long terminal repeats (LTRs) covered 44.0% of the genome with Copia families predominant at 22.9% followed by Gypsy at 14.2%. Both Gypsy and Copia elements were common in ancestral peaks of LTRs, and the most abundant element was a Gypsy element containing nested Copia/Angela sequence similarity, reflecting a complex evolutionary history of repeat activity. Gene annotation produced 33,257 genes and 68,942 transcripts, of which 99% were functionally annotated. BUSCO scores for the annotated proteins were 96.0% complete of which 77.6% was single copy and 18.4% duplicates. Whole genome duplication synonymous mutation rates of Gaviota tarplant and sunflower (Helianthus annuus) shared peaks that correspond to the last Asteraceae polyploidization event and subsequent divergence from a common ancestor at ∼27 MYA. Regions of high-density tandem genes were identified, pointing to potentially important loci of environmental adaptation in this species.

Chromosome-Scale Assembly and Annotation of Eight Arabidopsis thaliana Ecotypes

The plant Arabidopsis thaliana is a model system used by researchers through much of plant research. Recent efforts have focused on discovering the genomic variation found in naturally occurring ecotypes isolated from around the world. These ecotypes have come from diverse climates and therefore have faced and adapted to a variety of abiotic and biotic stressors. The sequencing and comparative analysis of these genomes can offer insight into the adaptive strategies of plants. While there are a large number of ecotype genome sequences available, the majority were created using short-read technology. Mapping of short-reads containing structural variation to a reference genome bereft of that variation leads to incorrect mapping of those reads, resulting in a loss of genetic information and introduction of false heterozygosity. For this reason, long-read de novo sequencing of genomes is required to resolve structural variation events. In this article, we sequenced the genomes of eight natural variants of A. thaliana using nanopore sequencing. This resulted in highly contiguous assemblies with >95% of the genome contained within five contigs. The sequencing results from this study include five ecotypes from relict and African populations, an area of untapped genetic diversity. With this study, we increase the knowledge of diversity we have across A. thaliana ecotypes and contribute to ongoing production of an A. thaliana pan-genome.

Genome sequencing and CAZymes repertoire analysis of Diaporthe eres P3-1W causing postharvest fruit rot of 'Hongyang' kiwifruit in China

Postharvest rot caused by various fungal pathogens is a damaging disease affecting kiwifruit production and quality, resulting in significant annual economic losses. This study focused on isolating the strain P3-1W, identified as Diaporthe eres, as the causal agent of 'Hongyang' postharvest rot disease in China. The investigation highlighted cell wall degrading enzymes (CWDEs) as crucial pathogenic factors. Specially, the enzymatic activities of cellulase, β-galactosidase, polygalacturonase, and pectin methylesterases peaked significantly on the second day after infection of D. eres P3-1W. To gain a comprehensive understanding of these CWDEs, the genome of this strain was sequenced using PacBio and Illumina sequencing technologies. The analysis revealed that the genome of D. eres P3-1W spans 58,489,835 bp, with an N50 of 5,939,879 bp and a GC content of 50.7%. A total of 15,407 total protein-coding genes (PCGs) were predicted and functionally annotated. Notably, 857 carbohydrate-active enzymes (CAZymes) were identified in D. eres P3-1W, with 521 CWDEs consisting of 374 glycoside hydrolases (GHs), 108 carbohydrate esterase (CEs) and 91 polysaccharide lyases (PLs). Additionally, 221 auxiliary activities (AAs), 91 glycosyltransferases (GTs), and 108 carbohydrate binding modules (CBMs) were detected. These findings offer valuable insights into the CAZymes of D. eres P3-1W.

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages.

Chromosome-level genome assembly of marmalade hoverfly Episyrphus balteatus (Diptera: Syrphidae)

Episyrphus balteatus can provide dual ecosystem services including pest control and pollination, which the larvae are excellent predators of aphid pest whereas adults are efficient pollinator. In this study, we assembled a high-quality genome of E. balteatus from northern China geographical population at the chromosome level by using Illumina, PacBio long reads, and Hi-C technologies. The 467.42 Mb genome was obtained from 723 contigs, with a contig N50 of 9.16 Mb and Scaffold N50 of 118.85 Mb, and 90.25% (431.75 Mb) of the assembly was anchored to 4 pseudo-autosomes and one pseudo-heterosome. In total, 14,848 protein-coding genes were annotated, and 95.14% of genes were fully represented in NR, GO, KEGG databases. Besides, we also obtained the mitochondrial genome of E. balteatus of 16, 837 bp in length with 37 typical mitochondrial genes. Overall, this high-quality genome is valuable for evolutionary and genetic studies of E. balteatus and other Syrphidae hoverfly species.

Convergent evolution of plant prickles by repeated gene co-option over deep time

An enduring question in evolutionary biology concerns the degree to which episodes of convergent trait evolution depend on the same genetic programs, particularly over long timescales. In this work, we genetically dissected repeated origins and losses of prickles-sharp epidermal projections-that convergently evolved in numerous plant lineages. Mutations in a cytokinin hormone biosynthetic gene caused at least 16 independent losses of prickles in eggplants and wild relatives in the genus Solanum. Homologs underlie prickle formation across angiosperms that collectively diverged more than 150 million years ago, including rice and roses. By developing new Solanum genetic systems, we leveraged this discovery to eliminate prickles in a wild species and an indigenously foraged berry. Our findings implicate a shared hormone activation genetic program underlying evolutionarily widespread and recurrent instances of plant morphological innovation.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

Chromosome-level genome assembly and annotation of the black sea urchin Arbacia lixula (Linnaeus, 1758)

The black sea urchin (Arbacia lixula) is a keystone species inhabiting the coastal shallow waters of the Mediterranean Sea, which is a key driver of littoral communities' structure. Here, we present the first genome assembly and annotation of this species, standing as the first Arbacioida genome, including both nuclear and mitochondrial genomes. To obtain a chromosome-level assembly, we used a combination of PacBio high fidelity (HiFi) reads and chromatin capture reads (Omni-C). In addition, we generated a high-quality nuclear annotation of both coding and non-coding genes, by using published RNA-Seq data from several individuals of A. lixula and gene models from closely related species. The nuclear genome assembly has a total span of 607.91 Mb, being consistent with its experimentally estimated genome size. The assembly contains 22 chromosome-scale scaffolds (96.52% of the total length), which coincides with its known karyotype. A total of 72,767 transcripts were predicted from the nuclear genome, 24,171 coding, and 48,596 non-coding that included lncRNA, snoRNA, and tRNAs. The circularized mitochondrial genome had 15,740 bp comprising 13 protein-coding genes, 2 rRNA, and 22 tRNA. This reference genome will enhance ongoing A. lixula studies and benefit the wider sea urchin scientific community.

JCVI: A versatile toolkit for comparative genomics analysis

The life cycle of genome builds spans interlocking pillars of assembly, annotation, and comparative genomics to drive biological insights. While tools exist to address each pillar separately, there is a growing need for tools to integrate different pillars of a genome project holistically. For example, comparative approaches can provide quality control of assembly or annotation; genome assembly, in turn, can help to identify artifacts that may complicate the interpretation of genome comparisons. The JCVI library is a versatile Python-based library that offers a suite of tools that excel across these pillars. Featuring a modular design, the JCVI library provides high-level utilities for tasks such as format parsing, graphics generation, and manipulation of genome assemblies and annotations. Supporting genomics algorithms like MCscan and ALLMAPS are widely employed in building genome releases, producing publication-ready figures for quality assessment and evolutionary inference. Developed and maintained collaboratively, the JCVI library emphasizes quality and reusability.

In-depth exploration of the genomic diversity in tea varieties based on a newly constructed pangenome of Camellia sinensis

Tea, one of the most widely consumed beverages globally, exhibits remarkable genomic diversity in its underlying flavour and health-related compounds. In this study, we present the construction and analysis of a tea pangenome comprising a total of 11 genomes, with a focus on three newly sequenced genomes comprising the purple-leaved assamica cultivar "Zijuan", the temperature-sensitive sinensis cultivar "Anjibaicha" and the wild accession "L618" whose assemblies exhibited excellent quality scores as they profited from latest sequencing technologies. Our analysis incorporates a detailed investigation of transposon complement across the tea pangenome, revealing shared patterns of transposon distribution among the studied genomes and improved transposon resolution with long read technologies, as shown by long terminal repeat (LTR) Assembly Index analysis. Furthermore, our study encompasses a gene-centric exploration of the pangenome, exploring the genomic landscape of the catechin pathway with our study, providing insights on copy number alterations and gene-centric variants, especially for Anthocyanidin synthases. We constructed a gene-centric pangenome by structurally and functionally annotating all available genomes using an identical pipeline, which both increased gene completeness and allowed for a high functional annotation rate. This improved and consistently annotated gene set will allow for a better comparison between tea genomes. We used this improved pangenome to capture the core and dispensable gene repertoire, elucidating the functional diversity present within the tea species. This pangenome resource might serve as a valuable resource for understanding the fundamental genetic basis of traits such as flavour, stress tolerance, and disease resistance, with implications for tea breeding programmes.

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

Diatom abundance in the polar oceans is predicted by genome size

A principal goal in ecology is to identify the determinants of species abundances in nature. Body size has emerged as a fundamental and repeatable predictor of abundance, with smaller organisms occurring in greater numbers than larger ones. A biogeographic component, known as Bergmann's rule, describes the preponderance, across taxonomic groups, of larger-bodied organisms in colder areas. Although undeniably important, the extent to which body size is the key trait underlying these patterns is unclear. We explored these questions in diatoms, unicellular algae of global importance for their roles in carbon fixation and energy flow through marine food webs. Using a phylogenomic dataset from a single lineage with worldwide distribution, we found that body size (cell volume) was strongly correlated with genome size, which varied by 50-fold across species and was driven by differences in the amount of repetitive DNA. However, directional models identified temperature and genome size, not cell size, as having the greatest influence on maximum population growth rate. A global metabarcoding dataset further identified genome size as a strong predictor of species abundance in the ocean, but only in colder regions at high and low latitudes where diatoms with large genomes dominated, a pattern consistent with Bergmann's rule. Although species abundances are shaped by myriad interacting abiotic and biotic factors, genome size alone was a remarkably strong predictor of abundance. Taken together, these results highlight the cascading cellular and ecological consequences of macroevolutionary changes in an emergent trait, genome size, one of the most fundamental and irreducible properties of an organism.

Water wisteria genome reveals environmental adaptation and heterophylly regulation in amphibious plants

Heterophylly is a phenomenon whereby an individual plant dramatically changes leaf shape in response to the surroundings. Hygrophila difformis (Acanthaceae; water wisteria), has recently emerged as a model plant to study heterophylly because of its striking leaf shape variation in response to various environmental factors. When submerged, H. difformis often develops complex leaves, but on land it develops simple leaves. Leaf complexity is also influenced by other factors, such as light density, humidity, and temperature. Here, we sequenced and assembled the H. difformis chromosome-level genome (scaffold N50: 60.43 Mb, genome size: 871.92 Mb), which revealed 36 099 predicted protein-coding genes distributed over 15 pseudochromosomes. H. difformis diverged from its relatives during the Oligocene climate-change period and expanded gene families related to its amphibious habit. Genes related to environmental stimuli, leaf development, and other pathways were differentially expressed in submerged and terrestrial conditions, possibly modulating morphological and physiological acclimation to changing environments. We also found that auxin plays a role in H. difformis heterophylly. Finally, we discovered candidate genes that respond to different environmental conditions and elucidated the role of LATE MERISTEM IDENTITY 1 (LMI1) in heterophylly. We established H. difformis as a model for studying interconnections between environmental adaptation and morphogenesis.

Polyploids of Brassicaceae: Genomic Insights and Assembly Strategies

The Brassicaceae family is distinguished by its inclusion of high-value crops such as cabbage, broccoli, mustard, and wasabi, all noted for their glucosinolates. In this family, many polyploidy species are distributed and shaped by numerous whole-genome duplications, independent genome doublings, and hybridization events. The evolutionary trajectory of the family is marked by enhanced diversification and lineage splitting after paleo- and meso-polyploidization, with discernible remnants of whole-genome duplications within their genomes. The recent neopolyploidization events notably increased the proportion of polyploid species within the family. Although sequencing efforts for the Brassicaceae genome have been robust, accurately distinguishing sub-genomes remains a significant challenge, frequently complicating the assembly process. Assembly strategies include comparative analyses with ancestral species and examining k-mers, long terminal repeat retrotransposons, and pollen sequencing. This review comprehensively explores the unique genomic characteristics of the Brassicaceae family, with a particular emphasis on polyploidization events and the latest strategies for sequencing and assembly. This review will significantly improve our understanding of polyploidy in the Brassicaceae family and assist in future genome assembly methods.

The genome of Vicia sativa ssp. amphicarpa provides insights into the role of terpenoids in antimicrobial resistance within subterranean fruits

Vicia sativa ssp. amphicarpa is a unique forage crop capable of simultaneously producing fruits above and below ground, representing a typical amphicarpic plant. In this study, we sequenced and assembled seven pseudo-chromosomes of the genome of V. sativa ssp. amphicarpa (n = 7) yielding a genome size of 1.59 Gb, with a total annotation of 48 932 protein-coding genes. Long terminal repeat (LTR) elements constituted 62.28% of the genome, significantly contributing to the expansion of genome size. Phylogenetic analysis revealed that the divergence between V. sativa ssp. amphicarpa and V. sativa was around 0.88 million years ago (MYA). Comparative transcriptomic and metabolomic analysis of aerial and subterranean pod shells showed biosynthesis of terpenoids in the subterranean pod shells indicating a correlation between the antimicrobial activity of subterranean pod shells and the biosynthesis of terpenoids. Furthermore, functional validation indicates that overexpression of VsTPS5 and VsTPS16 enhances terpenoid biosynthesis for antibacterial activity. Metabolomic analysis suggests the involvement of terpenoids in the antimicrobial properties of subterranean pod shells. Deciphering the genome of V. sativa ssp. amphicarpa elucidated the molecular mechanisms behind the antimicrobial properties of subterranean fruits in amphicarpic plants, providing valuable insights for the study of amphicarpic plant biology.

Chromosome-level genome assembly of Plagiognathops microlepis based on PacBio HiFi and Hi-C sequencing

Plagiognathops microlepis is an economic freshwater fish in the subfamily Xenocyprinae of Cyprinidae. It is widely distributed in the freshwater ecosystem of China, with moderate economic value and broad development prospects. However, the lack of genomic resources has limited our understanding on the genetic basis, phylogenetic status and adaptive evolution strategies of this fish. Here, we assembled a chromosome-level reference genome of P. microlepis by integrating Pacbio HiFi long-reads, Illumina short-reads and Hi-C sequencing data. The size of this genome is 1004.34 Mb with a contig N50 of 38.80 Mb. Using Hi-C sequencing data, 99.59% of the assembled sequences were further anchored to 24 chromosomes. A total of 578.91 Mb repeat sequences and 28,337 protein-coding genes were predicted in the current genome, of which, 26,929 genes were functionally annotated. This genome provides valuable information for investigating the phylogeny and evolutionary history of cyprinid fishes, as well as the genetic basis of adaptive strategies and special traits in P. microlepis.

Whole-genome duplication in an algal symbiont bolsters coral heat tolerance

The algal endosymbiont Durusdinium trenchii enhances the resilience of coral reefs under thermal stress. D. trenchii can live freely or in endosymbiosis, and the analysis of genetic markers suggests that this species has undergone whole-genome duplication (WGD). However, the evolutionary mechanisms that underpin the thermotolerance of this species are largely unknown. Here, we present genome assemblies for two D. trenchii isolates, confirm WGD in these taxa, and examine how selection has shaped the duplicated genome regions using gene expression data. We assess how the free-living versus endosymbiotic lifestyles have contributed to the retention and divergence of duplicated genes, and how these processes have enhanced the thermotolerance of D. trenchii. Our combined results suggest that lifestyle is the driver of post-WGD evolution in D. trenchii, with the free-living phase being the most important, followed by endosymbiosis. Adaptations to both lifestyles likely enabled D. trenchii to provide enhanced thermal stress protection to the host coral.

Chromosomal-level genome assembly and single-nucleotide polymorphism sites of black-faced spoonbill Platalea minor

Platalea minor, or black-faced spoonbill (Threskiornithidae), is a wading bird confined to coastal areas in East Asia. Due to habitat destruction, it was classified as globally endangered by the International Union for Conservation of Nature. However, the lack of genomic resources for this species hinders the understanding of its biology and diversity, and the development of conservation measures. Here, we report the first chromosomal-level genome assembly of P. minor using a combination of PacBio SMRT and Omni-C scaffolding technologies. The assembled genome (1.24 Gb) contains 95.33% of the sequences anchored to 31 pseudomolecules. The genome assembly has high sequence continuity with scaffold length N50 = 53 Mb. We predicted 18,780 protein-coding genes and measured high BUSCO score completeness (97.3%). Finally, we revealed 6,155,417 bi-allelic single nucleotide polymorphisms, accounting for ∼5% of the genome. This resource offers new opportunities for studying the black-faced spoonbill and developing conservation measures for this species.

Whole-genome sequencing of 13 Arctic plants and draft genomes of Oxyria digyna and Cochlearia groenlandica

To understand the genomic characteristics of Arctic plants, we generated 28-44 Gb of short-read sequencing data from 13 Arctic plants collected from the High Arctic Svalbard. We successfully estimated the genome sizes of eight species by using the k-mer-based method (180-894 Mb). Among these plants, the mountain sorrel (Oxyria digyna) and Greenland scurvy grass (Cochlearia groenlandica) had relatively small genome sizes and chromosome numbers. We obtained 45 × and 121 × high-fidelity long-read sequencing data. We assembled their reads into high-quality draft genomes (genome size: 561 and 250 Mb; contig N50 length: 36.9 and 14.8 Mb, respectively), and correspondingly annotated 43,105 and 29,675 genes using ~46 and ~85 million RNA sequencing reads. We identified 765,012 and 88,959 single-nucleotide variants, and 18,082 and 7,698 structural variants (variant size ≥ 50 bp). This study provided high-quality genome assemblies of O. digyna and C. groenlandica, which are valuable resources for the population and molecular genetic studies of these plants.

Chromosome-Scale Genome Assembly of the Sheep-Biting Louse Bovicola ovis Using Nanopore Sequencing Data and Pore-C Analysis

Bovicola ovis, commonly known as the sheep-biting louse, is an ectoparasite that adversely affects the sheep industry. Sheep louse infestation lowers the quality of products, including wool and leather, causing a loss of approximately AUD 123M per annum in Australia alone. The lack of a high-quality genome assembly for the sheep-biting louse, as well as any closely related livestock lice, has hindered the development of louse research and management control tools. In this study, we present the assembly of B. ovis with a genome size of ~123 Mbp based on a nanopore long-read sequencing library and Illumina RNA sequencing, complemented with a chromosome-level scaffolding using the Pore-C multiway chromatin contact dataset. Combining multiple alignment and gene prediction tools, a comprehensive annotation on the assembled B. ovis genome was conducted and recalled 11,810 genes as well as other genomic features including orf, ssr, rRNA and tRNA. A manual curation using alignment with the available closely related louse species, Pediculus humanus, increased the number of annotated genes to 16,024. Overall, this study reported critical genetic resources and biological insights for the advancement of sheep louse research and the development of sustainable control strategies in the sheep industry.