AUGUSTUS: ab initio prediction of alternative transcripts

Several secondary metabolite gene clusters in the genomes of ten Penicillium spp. raise the risk of multiple mycotoxin occurrence in chestnuts

Penicillium spp. produce a great variety of secondary metabolites, including several mycotoxins, on food substrates. Chestnuts represent a favorable substrate for Penicillium spp. development. In this study, the genomes of ten Penicillium species, virulent on chestnuts, were sequenced and annotated: P. bialowiezense. P. pancosmium, P. manginii, P. discolor, P. crustosum, P. palitans, P. viridicatum, P. glandicola, P. taurinense and P. terrarumae. Assembly size ranges from 27.5 to 36.8 Mb and the number of encoded genes ranges from 9,867 to 12,520. The total number of predicted biosynthetic gene clusters (BGCs) in the ten species is 551. The most represented families of BGCs are non ribosomal peptide synthase (191) and polyketide synthase (175), followed by terpene synthases (87). Genome-wide collections of gene phylogenies (phylomes) were reconstructed for each of the newly sequenced Penicillium species allowing for the prediction of orthologous relationships among our species, as well as other 20 annotated Penicillium species available in the public domain. We investigated in silico the presence of BGCs for 10 secondary metabolites, including 5 mycotoxins, whose production was validated in vivo through chemical analyses. Among the clusters present in this set of species we found andrastin A and its related cluster atlantinone A, mycophenolic acid, patulin, penitrem A and the cluster responsible for the synthesis of roquefortine C/glandicoline A/glandicoline B/meleagrin. We confirmed the presence of these clusters in several of the Penicillium species conforming our dataset and verified their capacity to synthesize them in a chestnut-based medium with chemical analysis. Interestingly, we identified mycotoxin clusters in some species for the first time, such as the andrastin A cluster in P. flavigenum and P. taurinense, and the roquefortine C cluster in P. nalgiovense and P. taurinense. Chestnuts proved to be an optimal substrate for species of Penicillium with different mycotoxigenic potential, opening the door to risks related to the occurrence of multiple mycotoxins in the same food matrix.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

SynGAP: a synteny-based toolkit for gene structure annotation polishing

Genome sequencing has become a routine task for biologists, but the challenge of gene structure annotation persists, impeding accurate genomic and genetic research. Here, we present a bioinformatics toolkit, SynGAP (Synteny-based Gene structure Annotation Polisher), which uses gene synteny information to accomplish precise and automated polishing of gene structure annotation of genomes. SynGAP offers exceptional capabilities in the improvement of gene structure annotation quality and the profiling of integrative gene synteny between species. Furthermore, an expression variation index is designed for comparative transcriptomics analysis to explore candidate genes responsible for the development of distinct traits observed in phylogenetically related species.

Genomic variation, environmental adaptation, and feralization in ramie, an ancient fiber crop

Feralization is an important evolutionary process, but the mechanisms behind it remain poorly understood. Here, we use the ancient fiber crop ramie (Boehmeria nivea (L.) Gaudich.) as a model to investigate genomic changes associated with both domestication and feralization. We first produced a chromosome-scale de novo genome assembly of feral ramie and investigated structural variations between feral and domesticated ramie genomes. Next, we gathered 915 accessions from 23 countries, comprising cultivars, major landraces, feral populations, and the wild progenitor. Based on whole-genome resequencing of these accessions, we constructed the most comprehensive ramie genomic variation map to date. Phylogenetic, demographic, and admixture signal detection analyses indicated that feral ramie is of exoferal or exo-endo origin, i.e., descended from hybridization between domesticated ramie and the wild progenitor or ancient landraces. Feral ramie has higher genetic diversity than wild or domesticated ramie, and genomic regions affected by natural selection during feralization differ from those under selection during domestication. Ecological analyses showed that feral and domesticated ramie have similar ecological niches that differ substantially from the niche of the wild progenitor, and three environmental variables are associated with habitat-specific adaptation in feral ramie. These findings advance our understanding of feralization, providing a scientific basis for the excavation of new crop germplasm resources and offering novel insights into the evolution of feralization in nature.

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.

A chromosome-scale Mytilus edulis genome assembly for aquaculture, marine ecology, and evolution

The smooth-shelled blue mussel, Mytilus edulis is part of the Mytilus species complex, encompassing at least three putative species: M. edulis, Mytilus galloprovincialis, and Mytilus trossulus. These three species occur on both sides of the Atlantic and hybridize in nature, and both M. edulis and M. galloprovincialis are important aquaculture species. They are also invasive species in many parts of the world. Here, we present a chromosome-level assembly of M. edulis. We used a combination of PacBio sequencing and Dovetail's Omni-C technology to generate an assembly with 14 long scaffolds containing 94% of the predicted length of the M. edulis genome (1.6 out of 1.7 Gb). Assembly statistics were as follows: total length = 1.65 Gb, N50 = 116 Mb, L50 = 7, and L90 = 13. BUSCO analysis showed 92.55% eukaryote BUSCOs identified. AB-Initio annotation using RNA-seq from mantle, gills, muscle, and foot predicted 47,128 genes. These gene models were combined with IsoSeq validation resulting in 45,379 full CDS protein sequences and 129,708 isoforms. Using GBS and shotgun sequencing, we also sequenced several eastern Canadian populations of Mytilus to characterize single-nucleotide as well as structural variance. This high-quality genome for M. edulis provides a platform to develop tools that can be used in breeding, molecular ecology and evolution to address questions of both commercial and environmental perspectives.

Improved assembly of the Pungitius pungitius reference genome

The nine-spined stickleback (Pungitius pungitius) has been increasingly used as a model system in studies of local adaptation and sex chromosome evolution but its current reference genome assembly is far from perfect, lacking distinct sex chromosomes. We generated an improved assembly of the nine-spined stickleback reference genome (98.3% BUSCO completeness) with the aid of linked-read mapping. While the new assembly (v8) was of similar size as the earlier version (v7), we were able to assign 4.4 times more contigs to the linkage groups and improve the contiguity of the genome. Moreover, the new assembly contains a ∼22.8 Mb Y-linked scaffold (LG22) consisting mainly of previously assigned X-contigs, putative Y-contigs, putative centromere contigs, and highly repetitive elements. The male individual showed an even mapping depth on LG12 (pseudo X chromosome) and LG22 (Y-linked scaffold) in the segregating sites, suggesting near-pure X and Y representation in the v8 assembly. A total of 26,803 genes were annotated, and about 33% of the assembly was found to consist of repetitive elements. The high proportion of repetitive elements in LG22 (53.10%) suggests it can be difficult to assemble the complete sequence of the species' Y chromosome. Nevertheless, the new assembly is a significant improvement over the previous version and should provide a valuable resource for genomic studies of stickleback fishes.

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.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Hi-C metagenomics facilitate comparative genome analysis of bacteria and yeast from spontaneous beer and cider

Sequence-based analysis of fermented foods and beverages' microbiomes offers insights into their impact on taste and consumer health. High-throughput metagenomics provide detailed taxonomic and functional community profiling, but bacterial and yeast genome reconstruction and mobile genetic elements tracking are to be improved. We established a pipeline for exploring fermented foods microbiomes using metagenomics coupled with chromosome conformation capture (Hi-C metagenomics). The approach was applied to analyze a collection of spontaneously fermented beers and ciders (n = 12). The Hi-C reads were used to reconstruct the metagenome-assembled genomes (MAGs) of bacteria and yeasts facilitating subsequent comparative genomic analysis, assembly scaffolding and exploration of "plasmid-bacteria" links. For a subset of beverages, yeasts were isolated and characterized phenotypically. The reconstructed Hi-C MAGs primarily belonged to the Lactobacillaceae family in beers, along with Acetobacteraceae and Enterobacteriaceae in ciders, exhibiting improved quality compared to conventional metagenomic MAGs. Comparative genomic analysis of Lactobacillaceae Hi-C MAGs revealed clustering by niche and suggested genetic determinants of survival and probiotic potential. For Pediococcus damnosus, Hi-C-based networks of contigs enabled linking bacteria with plasmids. Analyzing phylogeny and accessory genes in the context of known reference genomes offered insights into the niche specialization of beer lactobacilli. The subspecies-level diversity of cider Tatumella spp. was disentangled using a Hi-C-based graph. We obtained highly complete yeast Hi-C MAGs primarily represented by Brettanomyces and Saccharomyces, with Hi-C-facilitated chromosome-level genome assembly for the former. Utilizing Hi-C metagenomics to unravel the genomic content of individual species can provide a deeper understanding of the ecological interactions within the food microbiome, aid in bioprospecting beneficial microorganisms, improving quality control and improving innovative fermented products.

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.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Genome sequence data of Saccharomyces cerevisiae CBS 493.94

Whole genome sequencing (WGS) and data concerning identity and safety for Saccharomyces cerevisiae CBS 493.94 are reported. This strain was isolated from a British brewery in 1958 and deposited at the CBS culture collection Westerdijk Fungal Biodiversity Institute under the accession number CBS 493.94. The long-reads sequencing data, obtained via PacBio Sequel, and short-reads data, via Illumina NovaSeq 6000, were deposited at NCBI under accession number PRJNA1044661. The hybrid assembly was made publicly available via Zenodo and NCBI. For strain identification, data from 18S rRNA, ANI dendrogram and Core Genome single nucleotide polymorphism (SNP) Tree showed that the present isolate belongs to the genus Saccharomyces, species cerevisiae. The potential genes of concern, e.g. antimycotic resestance genes, were not detected. This strain is commonly used as a feed additive for animal health improvement and the present data summarise the unambiguous identity and strain's FKS1 gene does not code for any amino acid variants of concern.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Phylogenomic data exploration with increased sampling provides new insights into the higher-level relationships of butterflies and moths (Lepidoptera)

A robust and stable phylogenetic framework is a fundamental goal of evolutionary biology. As the third largest insect order in the world following Coleoptera and Diptera, Lepidoptera (butterflies and moths) play a central role in almost every terrestrial ecosystem as indicators of environmental change and serve as important models for biologists exploring questions related to ecology and evolutionary biology. However, for such a charismatic insect group, the higher-level phylogenetic relationships among its superfamilies are still poorly resolved. Compared to earlier phylogenomic studies, we increased taxon sampling among Lepidoptera (37 superfamilies and 68 families containing 263 taxa) and acquired a series of large amino-acid datasets from 69,680 to 400,330 for phylogenomic reconstructions. Using these datasets, we explored the effect of different taxon sampling with significant increases in the number of included genes on tree topology by considering a series of systematic errors using maximum-likelihood (ML) and Bayesian inference (BI) methods. Moreover, we also tested the effectiveness in topology robustness among the three ML-based models. The results showed that taxon sampling is an important determinant in tree robustness of accurate lepidopteran phylogenetic estimation. Long-branch attraction (LBA) caused by site-wise heterogeneity is a significant source of bias giving rise to unstable positions of ditrysian groups in phylogenomic reconstruction. Phylogenetic inference showed the most comprehensive framework to reveal the relationships among lepidopteran superfamilies, and presented some newly relationships with strong supports (Papilionoidea was sister to Gelechioidea and Immoidea was sister to Galacticoidea, respectively), but limited by taxon sampling, the relationships within the species-rich and relatively rapid radiation Ditrysia and especially Apoditrysia remain poorly resolved, which need to increase taxon sampling for further phylogenomic reconstruction. The present study demonstrates that taxon sampling is an important determinant for an accurate lepidopteran tree of life and provides some essential insights for future lepidopteran phylogenomic studies.

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

A chromosome-level genome assembly provides insights into the local adaptation of Tamarix austromongolica in the Yellow River Basin, China

Tamarix austromongolica is endemic to the Yellow River Basin and has adapted to diverse ecological settings in the region, including the arid areas of northwestern China and the saline soil regions of the Yellow River Delta. However, the genetic basis of its local adaptation remains unclear. We report a chromosome-level assembly of the T. austromongolica genome based on PacBio high-fidelity sequencing and Hi-C technology. The 12 pseudochromosomes cover 98.44% of the 1.32 Gb assembly, with a contig N50 of 52.57 Mb and a BUSCO score of 98.2%. The genome comprises 913.6 Mb (68.83%) of repetitive sequences and 22,374 protein-coding genes. Genome evolution analyses suggest that genes under positive selection and significantly expanded gene families have facilitated T. austromongolica's adaptability to diverse environmental factors and high resistance to diseases. Using genotyping-by-sequencing, we conducted population structure and selection analyses of 114 samples from 15 sites. Two genetic groups were identified, and 114 and 289 candidate genes were assigned to the populations of the northwestern and eastern parts of the Yellow River, respectively. Furthermore, we discovered numerous candidate genes associated with high-altitude adaptability and salt tolerance. This research provides valuable genomic resources for the evolutionary study and genetic breeding of tamarisk.

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.

Comparative genomics of Cryptosporidium parvum reveals the emergence of an outbreak-associated population in Europe and its spread to the United States

The zoonotic parasite Cryptosporidium parvum is a global cause of gastrointestinal disease in humans and ruminants. Sequence analysis of the highly polymorphic gp60 gene enabled the classification of C. parvum isolates into multiple groups (e.g., IIa, IIc, Id) and a large number of subtypes. In Europe, subtype IIaA15G2R1 is largely predominant and has been associated with many water- and food-borne outbreaks. In this study, we generated new whole-genome sequence (WGS) data from 123 human- and ruminant-derived isolates collected in 13 European countries and included other available WGS data from Europe, Egypt, China, and the United States (n = 72) in the largest comparative genomics study to date. We applied rigorous filters to exclude mixed infections and analyzed a data set from 141 isolates from the zoonotic groups IIa (n = 119) and IId (n = 22). Based on 28,047 high-quality, biallelic genomic SNPs, we identified three distinct and strongly supported populations: Isolates from China (IId) and Egypt (IIa and IId) formed population 1; a minority of European isolates (IIa and IId) formed population 2; and the majority of European (IIa, including all IIaA15G2R1 isolates) and all isolates from the United States (IIa) clustered in population 3. Based on analyses of the population structure, population genetics, and recombination, we show that population 3 has recently emerged and expanded throughout Europe to then, possibly from the United Kingdom, reach the United States, where it also expanded. The reason(s) for the successful spread of population 3 remain elusive, although genes under selective pressure uniquely in this population were identified.

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.

Genome assembly, Full-length transcriptome, and isoform diversity of Red Snapper, Lutjanus argentimaculatus

The mangrove red snapper, Lutjanus argentimaculatus, is a marine food fish of economic and aquaculture importance. The application of genomic selection-based breeding programs for this species is limited by the absence of a reference genome and transcriptome profiles. The current study attempted to fill this void by generating genomic and transcriptomic resources for red snapper. Using PacBio long reads, and Arima Hi-C linked reads, a scaffold-level genome assembly was generated for L. argentimaculatus. The assembly is of 1.03 Gb comprising of 400 scaffolds with N50 of 33.8 Mb and was assessed to be 97.2% complete upon benchmarking with BUSCO. Full-length transcriptome generated with PacBio Iso-Sequencing strategy using six tissues (muscle, gills, liver, kidney, stomach, and gonad) contained 56,515 isoforms belonging to 18,108 unique genes with N50 length of 3,973 bp. The resources generated will have potential applications in the functional studies, conservation, broodstock management and selective breeding programmes of L. argentimaculatus.

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.

A chromosome-level genome assembly for the paramylon-producing microalga Euglena gracilis

Euglena gracilis (E. gracilis), pivotal in the study of photosynthesis, endosymbiosis, and chloroplast development, is also an industrial microalga for paramylon production. Despite its importance, E. gracilis genome exploration faces challenges due to its intricate nature. In this study, we achieved a chromosome-level de novo assembly (2.37 Gb) using Illumina, PacBio, Bionano, and Hi-C data. The assembly exhibited a contig N50 of 619 Kb and scaffold N50 of 1.12 Mb, indicating superior continuity. Approximately 99.83% of the genome was anchored to 46 chromosomes, revealing structural insights. Repetitive elements constituted 58.84% of the sequences. Functional annotations were assigned to 39,362 proteins, enhancing interpretative power. BUSCO analysis confirmed assembly completeness at 80.39%. This first high-quality E. gracilis genome offers insights for genetics and genomics studies, overcoming previous limitations. The impact extends to academic and industrial research, providing a foundational resource.

Genomic and single-cell analyses reveal genetic signatures of swimming pattern and diapause strategy in jellyfish

Jellyfish exhibit innovative swimming patterns that contribute to exploring the origins of animal locomotion. However, the genetic and cellular basis of these patterns remains unclear. Herein, we generated chromosome-level genome assemblies of two jellyfish species, Turritopsis rubra and Aurelia coerulea, which exhibit straight and free-swimming patterns, respectively. We observe positive selection of numerous genes involved in statolith formation, hair cell ciliogenesis, ciliary motility, and motor neuron function. The lineage-specific absence of otolith morphogenesis- and ciliary movement-related genes in T. rubra may be associated with homeostatic structural statocyst loss and straight swimming pattern. Notably, single-cell transcriptomic analyses covering key developmental stages reveal the enrichment of diapause-related genes in the cyst during reverse development, suggesting that the sustained diapause state favours the development of new polyps under favourable conditions. This study highlights the complex relationship between genetics, locomotion patterns and survival strategies in jellyfish, thereby providing valuable insights into the evolutionary lineages of movement and adaptation in the animal kingdom.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

Combined genome and transcriptome provides insight into the genetic evolution of an edible halophyte Suaeda salsa adaptation to high salinity

Suaeda salsa L. is a typical halophyte with high value as a vegetable. Here, we report a 447.98 Mb, chromosomal-level genome of S. salsa, assembled into nine pseudomolecules (contig N50 = 1.36 Mb) and annotated with 27,927 annotated protein-coding genes. Most of the assembled S. salsa genome, 58.03%, consists of transposable elements. Some gene families including HKT1, NHX, SOS and CASP related to salt resistance were significantly amplified. We also observed expansion of genes encoding protein that bind the trace elements Zn, Fe, Cu and Mn, and genes related to flavonoid and α-linolenic acid metabolism. Many expanded genes were significantly up-regulated under salinity, which might have contributed to the acquisition of salt tolerance in S. salsa. Transcriptomic data showed that high salinity markedly up-regulated salt-resistance related genes, compared to low salinity. Abundant metabolic pathways of secondary metabolites including flavonoid, unsaturated fatty acids and selenocompound were enriched, which indicates that the species is a nutrient-rich vegetable. Particularly worth mentioning is that there was no significant difference in the numbers of cis-elements in the promoters of salt-related and randomly selected genes in S. salsa when compared with Arabidopsis thaliana, which may affirm that plant salt tolerance is a quantitative rather than a qualitative trait in terms of promoter evolution. Our findings provide deep insight into the adaptation of halophytes to salinity from a genetic evolution perspective.

EquCab_Finn: A new reference genome assembly for the domestic horse, Finnhorse

Finnhorse is Finland's native and national horse breed and it has genetic affinities to northern European and Asian horses. It has historical importance for agriculture, forest work and transport and as a war horse. Finnhorse has four breeding sections in the studbook and is under conservation and characterisation efforts. We sequenced and annotated the genome of a Finnhorse mare from the working horse section using PacBio and Omni-C data. This genome can complement the existing Thoroughbred reference genome (EquCab 3.0) and facilitate genetic studies of horses from northern Eurasia. We assembled 2.4 Gb of the genome with an N50 scaffold length of 83.8 Mb and the genome annotation resulted in a total of 19 748 protein coding genes of which 1200 were Finnhorse specific. The assembly has high quality and synteny with the current horse reference genome. We manually curated five genes of interest and deposited the final assembly in the European Nucleotide Archive under the accession no. PRJEB71364.

The miniature genome of broad mite, Polyphagotarsonemus latus (Tarsonemidae: Acari)

The broad mite, Polyphagotarsonemus latus (Tarsonemidae: Acari) is a highly polyphagous species that damage plant species spread across 57 different families. This pest has developed high levels of resistance to some commonly used acaricides. In the present investigation, we deciphered the genome information of P. latus by PacBio HiFi sequencing. P. latus is the third smallest arthropod genome sequenced so far with a size of 49.1 Mb. The entire genome was assembled into two contigs. A set of 9,286 protein-coding genes were annotated. Its compact genome size could be credited with multiple features such as very low repeat content (5.1%) due to the lack of proliferation of transposable elements, high gene density (189.1/Mb), more intronless genes (20.3%) and low microsatellite density (0.63%).

Unraveling the genomic secrets of Tritonibacter mobilis AK171: a plant growth-promoting bacterium isolated from Avicennia marina

The scarcity of freshwater resources resulting in a significant yield loss presents a pressing challenge in agriculture. To address this issue, utilizing abundantly available saline water could offer a smart solution. In this study, we demonstrate that the genome sequence rhizosphere bacterium Tritonibacter mobilis AK171, a halophilic marine bacterium recognized for its ability to thrive in saline and waterlogged environments, isolated from mangroves, has the remarkable ability to enable plant growth using saline irrigation. AK171 is characterized as rod-shaped cells, displays agile movement in free-living conditions, and adopts a rosette arrangement in static media. Moreover, The qualitative evaluation of PGP traits showed that AK171 could produce siderophores and IAA but could not solubilize phosphate nor produce hydrolytic enzymes it exhibits a remarkable tolerance to high temperatures and salinity. In this study, we conducted a comprehensive genome sequence analysis of T. mobilis AK171 to unravel the genetic mechanisms underlying its plant growth-promoting abilities in such challenging conditions. Our analysis revealed diverse genes and pathways involved in the bacterium's adaptation to salinity and waterlogging stress. Notably, T. mobilis AK171 exhibited a high level of tolerance to salinity and waterlogging through the activation of stress-responsive genes and the production of specific enzymes and metabolites. Additionally, we identified genes associated with biofilm formation, indicating its potential role in establishing symbiotic relationships with host plants. Furthermore, our analysis unveiled the presence of genes responsible for synthesizing antimicrobial compounds, including tropodithietic acid (TDA), which can effectively control phytopathogens. This genomic insight into T. mobilis AK171 provides valuable information for understanding the molecular basis of plant-microbial interactions in saline and waterlogged environments. It offers potential applications for sustainable agriculture in challenging conditions.

De Novo Long-Read Genome Assembly and Annotation of the Luna Moth (Actias luna) Fully Resolves Repeat-Rich Silk Genes

We present the first long-read de novo assembly and annotation of the luna moth (Actias luna) and provide the full characterization of heavy chain fibroin (h-fibroin), a long and highly repetitive gene (>20 kb) essential in silk fiber production. There are >160,000 described species of moths and butterflies (Lepidoptera), but only within the last 5 years have we begun to recover high-quality annotated whole genomes across the order that capture h-fibroin. Using PacBio HiFi reads, we produce the first high-quality long-read reference genome for this species. The assembled genome has a length of 532 Mb, a contig N50 of 16.8 Mb, an L50 of 14 contigs, and 99.4% completeness (BUSCO). Our annotation using Bombyx mori protein and A. luna RNAseq evidence captured a total of 20,866 genes at 98.9% completeness with 10,267 functionally annotated proteins and a full-length h-fibroin annotation of 2,679 amino acid residues.

Chromosome-level genome assembly of Helwingia omeiensis: the first genome in the family Helwingiaceae

Helwingia, a shrub of the monotypic cosmopolitan family Helwingiaceae, is distinguished by its inflorescence, in which flowers are borne on the midrib of the leaf-a trait not commonly observed in related plant families. Previous studies have investigated the development of this unusual structure using comparative anatomical methods. However, the scarcity of genomic data has hindered our understanding of the origins and evolutionary history of this uncommon trait at the molecular level. Here, we report the first high-quality genome of the family Helwingiaceae. Assembled using HiFi sequencing and Hi-C technologies, the genome of H. omeiensis is anchored to 19 chromosomes, with a total length of 2.75 Gb and a contig N50 length of 6.78 Mb. The BUSCO completeness score of the assembled genome was 98.2%. 53,951 genes were identified, of which 99.7% were annotated in at least one protein database. The high-quality reference genome of H. omeiensis provides an essential genetic resource and sheds light on the phylogeny and evolution of specific traits in the family Helwingiaceae.

Sirenian genomes illuminate the evolution of fully aquatic species within the mammalian superorder afrotheria

Sirenians of the superorder Afrotheria were the first mammals to transition from land to water and are the only herbivorous marine mammals. Here, we generated a chromosome-level dugong (Dugong dugon) genome. A comparison of our assembly with other afrotherian genomes reveals possible molecular adaptations to aquatic life by sirenians, including a shift in daily activity patterns (circadian clock) and tolerance to a high-iodine plant diet mediated through changes in the iodide transporter NIS (SLC5A5) and its co-transporters. Functional in vitro assays confirm that sirenian amino acid substitutions alter the properties of the circadian clock protein PER2 and NIS. Sirenians show evidence of convergent regression of integumentary system (skin and its appendages) genes with cetaceans. Our analysis also uncovers gene losses that may be maladaptive in a modern environment, including a candidate gene (KCNK18) for sirenian cold stress syndrome likely lost during their evolutionary shift in daily activity patterns. Genomes from nine Australian locations and the functionally extinct Okinawan population confirm and date a genetic break ~10.7 thousand years ago on the Australian east coast and provide evidence of an associated ecotype, and highlight the need for whole-genome resequencing data from dugong populations worldwide for conservation and genetic management.

Convergent reductive evolution in bee-associated lactic acid bacteria

Distantly related organisms may evolve similar traits when exposed to similar environments or engaging in certain lifestyles. Several members of the Lactobacillaceae (LAB) family are frequently isolated from the floral niche, mostly from bees and flowers. In some floral LAB species (henceforth referred to as bee-associated), distinctive genomic (e.g., genome reduction) and phenotypic (e.g., preference for fructose over glucose or fructophily) features were recently documented. These features are found across distantly related species, raising the hypothesis that specific genomic and phenotypic traits evolved convergently during adaptation to the floral environment. To test this hypothesis, we examined representative genomes of 369 species of bee-associated and non-bee-associated LAB. Phylogenomic analysis unveiled seven independent ecological shifts towards the floral niche in LAB. In these bee-associated LAB, we observed pervasive, significant reductions of genome size, gene repertoire, and GC content. Using machine leaning, we could distinguish bee-associated from non-bee-associated species with 94% accuracy, based on the absence of genes involved in metabolism, osmotic stress, or DNA repair. Moreover, we found that the most important genes for the machine learning classifier were seemingly lost, independently, in multiple bee-associated lineages. One of these genes, adhE, encodes a bifunctional aldehyde-alcohol dehydrogenase associated with the evolution of fructophily, a rare phenotypic trait that was recently identified in many floral LAB species. These results suggest that the independent evolution of distinctive phenotypes in bee-associated LAB has been largely driven by independent loss of the same set of genes.

Importance

Several lactic acid bacteria (LAB) species are intimately associated with bees and exhibit unique biochemical properties with potential for food applications and honeybee health. Using a machine-learning based approach, our study shows that adaptation of LAB to the bee environment was accompanied by a distinctive genomic trajectory deeply shaped by gene loss. Several of these gene losses occurred independently in distantly related species and are linked to some of their unique biotechnologically relevant traits, such as the preference of fructose over glucose (fructophily). This study underscores the potential of machine learning in identifying fingerprints of adaptation and detecting instances of convergent evolution. Furthermore, it sheds light onto the genomic and phenotypic particularities of bee-associated bacteria, thereby deepening the understanding of their positive impact on honeybee health.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

Haplotype-resolved assembly of the mule duck genome using high-fidelity sequencing technology

Mule duck is vitally important to the production of global duck meat. Here, we present two high-quality haplotypes of a female mule duck (haplotype 1 (H1):1.28 Gb, haplotype 2 (H2): 1.40 Gb). The continuity (H1: contig N50 = 14.90 Mb, H2: contig N50 = 15.70 Mb) and completeness (BUSCO: H1 = 96.9%, H2 = 97.3%) are substantially better than those of other duck genomes. We detected the structural variations (SVs) in H1 and H2. We observed a positive correlation between autosome length and the number of SVs. Z chromosome was some deficient in deletions and insertions, but W chromosome was some excessive. A total of 1,451 genes were haplotype specific expression (HSEs). Among them, 737 specifically expressed in H1, and 714 specifically expressed in H2. We found that H1 and H2 HSEs tended to be involved in similar biological processes, such as myometrial relaxation and contraction pathways, muscle structure development and phosphorylation. Our haplotype-resolved genome assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in mule duck.

Population genomics highlights structural variations in local adaptation to saline coastal environments in woolly grape

Structural variations (SVs) are a feature of plant genomes that has been largely unexplored despite their significant impact on plant phenotypic traits and local adaptation to abiotic and biotic stress. In this study, we employed woolly grape (Vitis retordii), a species native to the tropical and subtropical regions of East Asia with both coastal and inland habitats, as a valuable model for examining the impact of SVs on local adaptation. We assembled a haplotype-resolved chromosomal reference genome for woolly grape, and conducted population genetic analyses based on whole-genome sequencing (WGS) data from coastal and inland populations. The demographic analyses revealed recent bottlenecks in all populations and asymmetric gene flow from the inland to the coastal population. In total, 1,035 genes associated with plant adaptive regulation for salt stress, radiation, and environmental adaptation were detected underlying local selection by SVs and SNPs in the coastal population, of which 37.29% and 65.26% were detected by SVs and SNPs, respectively. Candidate genes such as FSD2, RGA1, and AAP8 associated with salt tolerance were found to be highly differentiated and selected during the process of local adaptation to coastal habitats in SV regions. Our study highlights the importance of SVs in local adaptation; candidate genes related to salt stress and climatic adaptation to tropical and subtropical environments are important genomic resources for future breeding programs of grapevine and its rootstocks.

Chromosome-level genome assembly of Guide Black-Fur sheep (Ovis aries)

Guide Black-Fur sheep (GD) is a breed of Tibetan sheep (Ovis aries) that lives in the Qinghai-Tibetan plateau region at an altitude of over 4,000 m. However, a lack of genomic information has made it difficult to understand the high-altitude adaptation of these sheep. We sequenced and assembled the GD reference genome using PacBio, Hi-C, and Illumina sequencing technologies. The final assembled genome size was 2.73 Gb, with a contig N50 of 20.30 Mb and a scaffold N50 of 107.63 Mb. The genome is predicted to contain 20,759 protein-coding genes, of which 98.42 have functional annotations. Repeat elements account for approximately 52.2% of the genomic landscape. The completeness of the GD genome assembly is highlighted by a BUSCO score of 93.1%. This high-quality genome assembly provides a critical resource for future molecular breeding and genetic improvement of Tibetan sheep.

Identifying genes within pathways in unannotated genomes with PaGeSearch

In biological research, the identification and comparison of genes within specific pathways across the genomes of various species are invaluable. However, annotating the entire genome is resource intensive, and sequence similarity searches often yield results that are not actually genes. To address these limitations, we introduce Pathway Gene Search (PaGeSearch), a tool designed to identify genes from predefined lists, especially those in specific pathways, within genomes. The tool uses an initial sequence similarity search to identify relevant genomic regions, followed by targeted gene prediction and neural network-based result filtering. PaGeSearch suggests the regions that are most likely the orthologs of the genes in the query and is designed to be applicable for species within five classes: mammals, fish, birds, eudicotyledons, and Liliopsida. Compared with GeMoMa and miniprot, PaGeSearch generally outperforms in terms of sensitivity and positive predictive value, as well as negative predictive value. Also, the exon coverage of gene models from PaGeSearch is higher compared with those in GeMoMa and miniprot. Although its performance shows increased variability when applied to actual biological pathways, it nonetheless maintains an acceptable level of accuracy. Evaluating PaGeSearch across different assembly levels, chromosome, scaffold, and contig shows minimal variation in outcomes, indicating that PaGeSearch is resilient to variations in assembly quality.

Chromosome-level genome assembly of a cliff plant Taihangia rupestris var. ciliata provides insights into its adaptation and demographic history

BACKGROUND

Cliffs are recognized as one of the most challenging environments for plants, characterized by harsh conditions such as drought, infertile soil, and steep terrain. However, they surprisingly host ancient and diverse plant communities and play a crucial role in protecting biodiversity. The Taihang Mountains, which act as a natural boundary in eastern China, support a rich variety of plant species, including many unique to cliff habitats. However, it is little known how cliff plants adapt to harsh habitats and the demographic history in this region.

RESULTS

To better understand the demographic history and adaptation of cliff plants in this area, we analyzed the chromosome-level genome of a representative cliff plant, T. rupestris var. ciliata, which has a genome size of 769.5 Mb, with a scaffold N50 of 104.92 Mb. The rapid expansion of transposable elements may have contributed to the increasing genome and its ability to adapt to unique and challenging cliff habitats. Comparative analysis of the genome evolution between Taihangia and non-cliff plants in Rosaceae revealed a significant expansion of gene families associated with oxidative phosphorylation, which is likely a response to the abiotic stresses faced by cliff plants. This expansion may explain the long-term adaptation of Taihangia to harsh cliff environments. The effective population size of the two varieties has continuously decreased due to climatic fluctuations during the Quaternary period. Furthermore, significant differences in gene expression between the two varieties may explain the varied leaf phenotypes and adaptations to harsh conditions in different natural distributions.

CONCLUSION

Our study highlights the extraordinary adaptation of T. rupestris var. ciliata, shedding light on the evolution of cliff plants worldwide.

BRAKER3: Fully automated genome annotation using RNA-seq and protein evidence with GeneMark-ETP, AUGUSTUS, and TSEBRA

Gene prediction has remained an active area of bioinformatics research for a long time. Still, gene prediction in large eukaryotic genomes presents a challenge that must be addressed by new algorithms. The amount and significance of the evidence available from transcriptomes and proteomes vary across genomes, between genes, and even along a single gene. User-friendly and accurate annotation pipelines that can cope with such data heterogeneity are needed. The previously developed annotation pipelines BRAKER1 and BRAKER2 use RNA-seq or protein data, respectively, but not both. A further significant performance improvement integrating all three data types was made by the recently released GeneMark-ETP. We here present the BRAKER3 pipeline that builds on GeneMark-ETP and AUGUSTUS, and further improves accuracy using the TSEBRA combiner. BRAKER3 annotates protein-coding genes in eukaryotic genomes using both short-read RNA-seq and a large protein database, along with statistical models learned iteratively and specifically for the target genome. We benchmarked the new pipeline on genomes of 11 species under an assumed level of relatedness of the target species proteome to available proteomes. BRAKER3 outperforms BRAKER1 and BRAKER2. The average transcript-level F1-score is increased by about 20 percentage points on average, whereas the difference is most pronounced for species with large and complex genomes. BRAKER3 also outperforms other existing tools, MAKER2, Funannotate, and FINDER. The code of BRAKER3 is available on GitHub and as a ready-to-run Docker container for execution with Docker or Singularity. Overall, BRAKER3 is an accurate, easy-to-use tool for eukaryotic genome annotation.

Chromosome-level genome of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus

Recent conservation efforts to protect rare and endangered aquatic species have intensified. Nevertheless, the ornate spiny lobster (Panulirus ornatus), which is prevalent in the Indo-Pacific waters, has been largely ignored. In the absence of a detailed genomic reference, the conservation and population genetics of this crustacean are poorly understood. Here, We assembled a comprehensive chromosome-level genome for P. ornatus. This genome-among the most detailed for lobsters-spans 2.65 Gb with a contig N50 of 51.05 Mb, and 99.11% of the sequences with incorporated to 73 chromosomes. The ornate spiny lobster genome comprises 65.67% repeat sequences and 22,752 protein-coding genes with 99.20% of the genes functionally annotated. The assembly of the P. ornatus genome provides valuable insights into comparative crustacean genomics and endangered species conservation, and lays the groundwork for future research on the speciation, ecology, and evolution of the ornate spiny lobster.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

Chromosome-level genome assembly of Euphorbia tirucalli (Euphorbiaceae), a highly stress-tolerant oil plant

Euphorbia, one of the largest genera of flowering plants, is well-known for containing many biofuel crops. Euphorbia tirucalli, an evergreen succulent mainly native to the Africa continent but cultivated worldwide, is a promising petroleum plant with high tolerance to drought and salt stress. However, the exploration of such an important plant resource is severely hampered by the lack of a reference genome. Here, we present the chromosome-level genome assembly of E. tirucalli using PacBio HiFi sequencing and Hi-C technology. Its genome size was approximately 745.62 Mb, with a contig N50 of 74.16 Mb. A total of 743.63 Mb (99.73%) of the assembled sequences were anchored to 10 chromosomes with a complete BUSCO score of 97.80%. Genome annotation revealed 26,304 protein-coding genes, and 76.37% of the genome was identified as repeat elements. The high-quality genome provides valuable genetic resources that would be useful for unraveling the genetic mechanisms of biofuel synthesis and evolutionary adaptation of E. tirucalli.

Description of new micro-colonial fungi species Neophaeococcomyces mojavensis, Coniosporium tulheliwenetii, and Taxawa tesnikishii cultured from biological soil crusts

Black yeasts and relatives comprise Micro-Colonial Fungi (MCFs) which are slow-growing stress-tolerant micro-eukaryotes that specialize in extreme environments. MCFs are paraphyletic and found in the Orders Chaetothyriales (Eurotiomycetes) and Dothideales (Dothidiomycetes). We have isolated and described three new MCFs species from desert biological soil crusts (BSCs) collected from two arid land regions: Joshua Tree National Park (Mojave Desert) and UC Natural Reserve at Boyd Deep Canyon (confluence of Mojave and Sonoran Deserts). BSCs are composite assemblages of cyanobacteria, eukaryotic algae, fungi, lichens, and bryophytes embedded into the surface of desert soils, providing a protective buffer against the harsh desert environment. Our work focused on one type of desert BSC, the cyanolichen crust dominated by Collema sp. Using culture-dependent protocols, three MCFs were axenically isolated from their respective samples along with the extracted DNA. Their genomes were sequenced using Illumina and Nanopore, and finally assembled and annotated using hybrid assembly approaches and established bioinformatics pipelines to conduct final taxonomic phylogenetic analysis and placement. The three species described here are unique specimen from desert BSCs, here we introduce, Neophaeococcomyces mojavensis (Chaetothyriales), Cladosporium tulheliwenetii (Dothideales), and Taxawa tesnikishii (Dothideales).

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

Genetic Diversity and Population Structure of Rhodeus uyekii in the Republic of Korea Revealed by Microsatellite Markers from Whole Genome Assembly

This study is the first report to characterize the Rhodus uyekii genome and study the development of microsatellite markers and their markers applied to the genetic structure of the wild population. Genome assembly was based on PacBio HiFi and Illumina HiSeq paired-end sequencing, resulting in a draft genome assembly of R. uyekii. The draft genome was assembled into 2652 contigs. The integrity assessment of the assemblies indicates that the quality of the draft assemblies is high, with 3259 complete BUSCOs (97.2%) in the database of Verbrata. A total of 31,166 predicted protein-coding genes were annotated in the protein database. The phylogenetic tree showed that R. uyekii is a close but distinct relative of Onychostoma macrolepis. Among the 10 fish genomes, there were significant gene family expansions (8-2387) and contractions (16-2886). The average number of alleles amplified by the 21 polymorphic markers ranged from 6 to 23, and the average PIC value was 0.753, which will be useful for evolutionary and genetic analysis. Using population genetic analysis, we analyzed genetic diversity and the genetic structures of 120 individuals from 6 populations. The average number of alleles per population ranged from 7.6 to 9.9, observed heterozygosity ranged from 0.496 to 0.642, and expected heterozygosity ranged from 0.587 to 0.783. Discriminant analysis of principal components According to the analysis method, the population was divided into three populations (BS vs. DC vs. GG, GC, MS, DC). In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of R. uyekii.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Chromosome-level genome assembly of the cereal cyst nematode Heterodera flipjevi

As an economically important plant parasitic nematode (PPN), Heterodera filipjevi causes great damage on wheat, and now it was widely recorded in many countries. While multiple genomes of PPNs have been published, high-quality genome assembly and annotation on H. filipjevi have yet to be performed. This study presents a chromosome-scale genome assembly and annotation for H. filipjevi, utilizing a combination of Illumina short-read, PacBio long-read, and Hi-C sequencing technologies. The genome consists of 9 pseudo-chromosomes that contain 134.19 Mb of sequence, with a scaffold N50 length of 11.88 Mb. In total, 10,036 genes were annotated, representing 75.20% of the total predicted protein-coding genes. Our study provides the first chromosome-scale genome for H. filipjevi, which is also the inaugural high-quality genome of cereal cyst nematodes (CCNs). It provides a valuable genomic resource for further biological research and pest management of cereal cyst nematodes disease.