Sequence locally, think globally: The Darwin Tree of Life Project

Genomics of insect natural enemies in agroecosystems

Currently, a wealth of genomic data are now accessible for numerous insect natural enemies, serving as valuable resources that deepen our understanding of the genetic basis of biocontrol traits in these organisms. We summarize the current state of genome sequencing and highlight candidate genes related to biocontrol traits that hold promise for genetic improvement. We also review the recent population genomic studies in biological control and the discovery of potential insecticidal genes in parasitoid wasps. Collectively, current genomic works have shown the powerful ability to identify candidate genes responsible for desirable traits or promising effectors. However, further functional study is necessary to gain a mechanistic understanding of these genes, and future efforts are also needed to develop suitable approaches to translate genomic insights into field applications.

Advances in euglenoid genomics: unravelling the fascinating biology of a complex clade

Euglenids have long been studied due to their unique physiology and versatile metabolism, providing underpinnings for much of our understanding of photosynthesis and biochemistry, and a growing opportunity in biotechnology. Until recently there has been a lack of genetic studies due to their large and complex genomes, but recently new technologies have begun to unveil their genetic capabilities. Whilst much research has focused on the model organism Euglena gracilis, other members of the euglenids have now started to receive due attention. Currently only poor nuclear genome assemblies of E. gracilis and Rhabdomonas costata are available, but there are many more plastid genome sequences and an increasing number of transcriptomes. As more assemblies become available, there are great opportunities to understand the fundamental biology of these organisms and to exploit them for biotechnology.

Gene novelty and gene family expansion in the early evolution of Lepidoptera

BACKGROUND

Almost 10% of all known animal species belong to Lepidoptera: moths and butterflies. To understand how this incredible diversity evolved we assess the role of gene gain in driving early lepidopteran evolution. Here, we compared the complete genomes of 115 insect species, including 99 Lepidoptera, to search for novel genes coincident with the emergence of Lepidoptera.

RESULTS

We find 217 orthogroups or gene families which emerged on the branch leading to Lepidoptera; of these 177 likely arose by gene duplication followed by extensive sequence divergence, 2 are candidates for origin by horizontal gene transfer, and 38 have no known homology outside of Lepidoptera and possibly arose via de novo gene genesis. We focus on two new gene families that are conserved across all lepidopteran species and underwent extensive duplication, suggesting important roles in lepidopteran biology. One encodes a family of sugar and ion transporter molecules, potentially involved in the evolution of diverse feeding behaviours in early Lepidoptera. The second encodes a family of unusual propeller-shaped proteins that likely originated by horizontal gene transfer from Spiroplasma bacteria; we name these the Lepidoptera propellin genes.

CONCLUSION

We provide the first insights into the role of genetic novelty in the early evolution of Lepidoptera. This gives new insight into the rate of gene gain during the evolution of the order as well as providing context on the likely mechanisms of origin. We describe examples of new genes which were retained and duplicated further in all lepidopteran species, suggesting their importance in Lepidoptera evolution.

Accurate de novo transcription unit annotation from run-on and sequencing data

Functional element annotations are critical tools used to provide insight into the molecular processes governing cell development, differentiation, and disease. Run-on and sequencing assays measure the production of nascent RNAs and can provide an effective data source for discovering functional elements. However, the accurate inference of functional elements from run-on sequencing data remains an open problem because the signal is noisy and challenging to model. Here we investigated computational approaches that convert run-on and sequencing data into annotations representing transcription units, including genes and non-coding RNAs. We developed a convolutional neural network, called c onvolutional discovery of g ene a natomy using P RO-seq (CGAP), trained to identify different anatomical features of a transcription unit, which were then stitched together into transcript annotations using a hidden Markov model (HMM). Comparison with existing methods showed a significant performance improvement using our novel CGAP-HMM approach. We developed a voting system that ensembles the top three annotation strategies, resulting in large and significant improvements in transcription unit annotation accuracy over the best performing individual method. Finally, we also report a conditional generative adversarial network (cGAN) as a generative approach to transcription unit annotation that shows promise for further development. Collectively our work provides novel tools for de novo transcription unit annotation from run-on and sequencing data that are accurate enough to be useful in many applications.

Identification of the shortest species-specific oligonucleotide sequences

Despite the exponential increase in sequencing information driven by massively parallel DNA sequencing technologies, universal and succinct genomic fingerprints for each organism are still missing. Identifying the shortest species-specific nucleotide sequences offers insights into species evolution and holds potential practical applications in agriculture, wildlife conservation, and healthcare. We propose a new method for sequence analysis termed nucleic "quasi-primes," the shortest occurring sequences in each of 45,076 organismal reference genomes, present in one genome and absent from every other examined genome. In the human genome, we find that the genomic loci of nucleic quasi-primes are most enriched for genes associated with brain development and cognitive function. In a single-cell case study focusing on the human primary motor cortex, nucleic quasi-prime genes account for a significantly larger proportion of the variation based on average gene expression. Nonneuronal cell types, including astrocytes, endothelial cells, microglia perivascular-macrophages, oligodendrocytes, and vascular and leptomeningeal cells, exhibit significant activation of quasi-prime-containing gene associations related to cancer, whereas simultaneously suppressing quasi-prime-containing genes are associated with cognitive, mental, and developmental disorders. We also show that human disease-causing variants, eQTLs, mQTLs, and sQTLs are 4.43-fold, 4.34-fold, 4.29-fold, and 4.21-fold enriched at human quasi-prime loci, respectively. These findings indicate that nucleic quasi-primes are genomic loci linked to the evolution of species-specific traits, and in humans, they provide insights in the development of cognitive traits and human diseases, including neurodevelopmental disorders.

k-mer approaches for biodiversity genomics

The wide array of currently available genomes displays a wonderful diversity in size, composition, and structure and is quickly expanding thanks to several global biodiversity genomics initiatives. However, sequencing of genomes, even with the latest technologies, can still be challenging for both technical (e.g., small physical size, contaminated samples, or access to appropriate sequencing platforms) and biological reasons (e.g., germline-restricted DNA, variable ploidy levels, sex chromosomes, or very large genomes). In recent years, k-mer-based techniques have become popular to overcome some of these challenges. They are based on the simple process of dividing the analyzed sequences (e.g., raw reads or genomes) into a set of subsequences of length k, called k-mers, and then analyzing the frequency or sequences of those k-mers. Analyses based on k-mers allow for a rapid and intuitive assessment of complex sequencing data sets. Here, we provide a comprehensive review to the theoretical properties and practical applications of k-mers in biodiversity genomics with a special focus on genome modeling.

Telomere DNA in the insect order Dermaptera and the first evidence for the non-canonical telomeric motif TTCGG in Arthropoda

Despite recent advances in telomere research, the telomere DNA organization remains unknown for representatives of several insect orders. In this study, analysis of the chromosome-level genome assembly shows that the telomeric DNA of the earwig Labiaminor (Linnaeus, 1758) (Polyneoptera, Dermaptera, Spongiphoridae) consists of repeats of the 5 bp motif TTCGG/CCGAA. This is the first record describing the structure of telomeric DNA in the order Dermaptera. This record expands the spectrum of the known telomeric sequences, since the TTCGG motif has not been reported for insects previously.

tidk: a toolkit to rapidly identify telomeric repeats from genomic datasets

SUMMARY

"tidk" (short for telomere identification toolkit) uses a simple, fast algorithm to scan long DNA reads for the presence of short tandemly repeated DNA in runs, and to aggregate them based on canonical DNA string representation. These are telomeric repeat candidates. Our algorithm is shown to be accurate in genomes for which the telomeric repeat unit is known and is tested across a wide variety of newly assembled genomes to uncover new telomeric repeat units. Tools are provided to identify telomeric repeats de novo, scan genomes for known telomeric repeats, and to visualize telomeric repeats on the assembly. "tidk" is implemented in Rust and is available as a command line tool which can be compiled using the Rust toolchain or downloaded as a binary from bioconda.

AVAILABILITY AND IMPLEMENTATION

The "tidk" Rust crate is freely available under the MIT license (https://crates.io/crates/tidk), and the source code is available at https://github.com/tolkit/telomeric-identifier.

Comparative Genomics Reveals LINE-1 Recombination with Diverse RNAs

Long interspersed element-1 (LINE-1, L1) retrotransposons are the most abundant protein-coding transposable elements (TE) in mammalian genomes, and have shaped genome content over 170 million years of evolution. LINE-1 is self-propagating and mobilizes other sequences, including Alu elements. Occasionally, LINE-1 forms chimeric insertions with non-coding RNAs and mRNAs. U6 spliceosomal small nuclear RNA/LINE-1 chimeras are best known, though there are no comprehensive catalogs of LINE-1 chimeras. To address this, we developed TiMEstamp, a computational pipeline that leverages multiple sequence alignments (MSA) to estimate the age of LINE-1 insertions and identify candidate chimeric insertions where an adjacent sequence arrives contemporaneously. Candidates were refined by detecting hallmark features of L1 retrotransposition, such as target site duplication (TSD). Applying this pipeline to the human genome, we recovered all known species of LINE-1 chimeras and discovered new chimeric insertions involving small RNAs, Alu elements, and mRNA fragments. Some insertions are compatible with known mechanisms, such as RNA ligation. Other structures nominate novel mechanisms, such as trans-splicing. We also see evidence that LINE-1 loci with defunct promoters can acquire regulatory elements from nearby genes to restore retrotransposition activity. These discoveries highlight the recombinatory potential of LINE-1 RNA with implications for genome evolution and TE domestication.

High-quality sika deer omics data and integrative analysis reveal genic and cellular regulation of antler regeneration

The antler is the only organ that can fully regenerate annually in mammals. However, the regulatory pattern and mechanism of gene expression and cell differentiation during this process remain largely unknown. Here, we obtain comprehensive assembly and gene annotation of the sika deer (Cervus nippon) genome. We construct, together with large-scale chromatin accessibility and gene expression data, gene regulatory networks involved in antler regeneration, identifying four transcription factors, MYC, KLF4, NFE2L2, and JDP2, with high regulatory activity across the whole regeneration process. Comparative studies and luciferase reporter assay suggest the MYC expression driven by a cervid-specific regulatory element might be important for antler regenerative ability. We further develop a model called combinatorial TF Oriented Program (cTOP), which integrates single-cell data with bulk regulatory networks and find PRDM1, FOSL1, BACH1, and NFATC1 as potential pivotal factors in antler stem cell activation and osteogenic differentiation. Additionally, we uncover interactions within and between cell programs and pathways during the regeneration process. These findings provide insights into the gene and cell regulatory mechanisms of antler regeneration, particularly in stem cell activation and differentiation.

Fusion, fission, and scrambling of the bilaterian genome in Bryozoa

Groups of orthologous genes are commonly found together on the same chromosome over vast evolutionary distances. This extensive physical gene linkage, known as macrosynteny, is seen between bilaterian phyla as divergent as Chordata, Echinodermata, Mollusca, and Nemertea. Here, we report a unique pattern of genome evolution in Bryozoa, an understudied phylum of colonial invertebrates. Using comparative genomics, we reconstruct the chromosomal evolutionary history of five bryozoans. Multiple ancient chromosome fusions followed by gene mixing led to the near-complete loss of bilaterian linkage groups in the ancestor of extant bryozoans. A second wave of rearrangements, including chromosome fission, then occurred independently in two bryozoan classes, further scrambling bryozoan genomes. We also discover at least five derived chromosomal fusion events shared between bryozoans and brachiopods, supporting the traditional but highly debated Lophophorata hypothesis and suggesting macrosynteny to be a potentially powerful source of phylogenetic information. Finally, we show that genome rearrangements led to the dispersion of genes from bryozoan Hox clusters onto multiple chromosomes. Our findings demonstrate that the canonical bilaterian genome structure has been lost across all studied representatives of an entire phylum, and reveal that linkage group fission can occur very frequently in specific lineages.

The genome sequence of the Ruddy Darter, Sympetrum sanguineum (Müller, 1764)

We present a genome assembly from a male specimen of Sympetrum sanguineum (Ruddy Darter; Arthropoda; Insecta; Odonata; Libellulidae). The haplotype-resolved assembly contains two haplotypes with total lengths of 1,500.53 megabases and 1,304.05 megabases. Most of haplotype 1 is scaffolded into 13 chromosomal pseudomolecules, including the X sex chromosome, while haplotype 2 is scaffolded into 12 autosomes.

The genome sequence of the tawny cockroach, Ectobius (Ectobius) pallidus (Olivier, 1789)

We present a genome assembly from a specimen of Ectobius pallidus (tawny cockroach; Arthropoda; Insecta; Blattodea; Ectobiidae). The assembly contains two haplotypes with total lengths of 2,087.55 megabases and 2,124.67 megabases, respectively. Most of haplotype 1 (98.55%) is scaffolded into 11 chromosomal pseudomolecules, while haplotype 2 is assembled to scaffold level. The mitochondrial genome has also been assembled and is 15.75 kilobases in length.

The genome sequence of a sawfly, Athalia cordata Serville, 1823

We present a genome assembly from a female specimen of Athalia cordata (sawfly; Arthropoda; Insecta; Hymenoptera; Athaliidae). The genome sequence has a total length of 169.00 megabases. Most of the assembly (99.98%) is scaffolded into 4 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 27.47 kilobases in length.

The genome sequence of a cranefly, Diogma glabrata (Meigen, 1818)

We present a genome assembly from a specimen of Diogma glabrata (cranefly; Arthropoda; Insecta; Diptera; Cylindrotomidae). The genome sequence has a total length of 1,328.70 megabases. Most of the assembly (90.7%) is scaffolded into 4 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 17.5 kilobases in length.

The genome sequence of a leafhopper, Allygus modestus Scott, 1876

We present a genome assembly from an individual male specimen of Allygus modestus (leafhopper; Arthropoda; Insecta; Hemiptera; Cicadellidae). The genome sequence has a total length of 1,819.90 megabases. Most of the assembly (99.86%) is scaffolded into 7 chromosomal pseudomolecules, including the X sex chromosome. The mitochondrial genome has also been assembled and is 16.69 kilobases in length.

The genome sequence of a sawfly, Abia candens Konow, 1887

We present a genome assembly from an individual male specimen of Abia candens (sawfly; Arthropoda; Insecta; Hymenoptera; Cimbicidae). The genome sequence has a total length of 261.00 megabases. Most of the assembly (82.7%) is scaffolded into 16 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 19.72 kilobases in length.

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

TIPPo: A User-Friendly Tool for De Novo Assembly of Organellar Genomes with High-Fidelity Data

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

The 2025 Nucleic Acids Research database issue and the online molecular biology database collection

The 2025 Nucleic Acids Research database issue contains 185 papers spanning biology and related areas. Seventy three new databases are covered, while resources previously described in the issue account for 101 update articles. Databases most recently published elsewhere account for a further 11 papers. Nucleic acid databases include EXPRESSO for multi-omics of 3D genome structure (this issue's chosen Breakthrough Resource and Article) and NAIRDB for Fourier transform infrared data. New protein databases include structure predictions for human isoforms at ASpdb and for viral proteins at BFVD. UniProt, Pfam and InterPro have all provided updates: metabolism and signalling are covered by new descriptions of STRING, KEGG and CAZy, while updated microbe-oriented databases include Enterobase, VFDB and PHI-base. Biomedical research is supported, among others, by ClinVar, PubChem and DrugMAP. Genomics-related resources include Ensembl, UCSC Genome Browser and dbSNP. New plant databases cover the Solanaceae (SolR) and Asteraceae (AMIR) families while an update from NCBI Taxonomy also features. The Database Issue is freely available on the Nucleic Acids Research website (https://academic.oup.com/nar). At the NAR online Molecular Biology Database Collection (http://www.oxfordjournals.org/nar/database/c/), 932 entries have been reviewed in the last year, 74 new resources added and 226 discontinued URLs eliminated bringing the current total to 2236 databases.

Evaluation of controls, quality control assays, and protocol optimisations for PacBio HiFi sequencing on diverse and challenging samples

The Darwin Tree of Life (DToL) project aims to generate high-quality reference genomes for all eukaryotic organisms in Britain and Ireland. At the time of writing, PacBio HiFi reads are generated for all samples using the Sequel IIe systems by the Wellcome Sanger Institute's Scientific Operations teams, however we expect lessons from this work to apply directly to the Revio system too, as core principles of SMRT sequencing remain the same. We observed that HiFi yield is highly variable for DToL samples. We have investigated what drives this variation, and potential mitigations. To support these investigations a number of controls were evaluated to ensure that the library and sequencing preparation procedures, reagents, consumables, and Sequel IIe instruments, were performing as expected. Our findings support that a primary factor driving variability in HiFi yield is the quality of the DNA prior to library construction, e.g., purity, size, and damage. We investigated whether quality assessment assays could link measurable DNA damage or purity to sequencing yield. Some correlation could be established, however no assay was predictive of sequencing yield for all samples, indicating that the variability is driven by multiple factors that may interact. We demonstrate that contaminants present in some samples are the cause of very low HiFi yield, and show that these contaminants can negatively affect the PacBio internal sequencing control and samples multiplexed on the same SMRT Cell. We found that consistently high yields could be obtained if an amplification workflow was utilised, namely PacBio's ultra-low input library preparation protocol.

Ensembl 2025

Ensembl (www.ensembl.org) is an open platform integrating publicly available genomics data across the tree of life with a focus on eukaryotic species related to human health, agriculture and biodiversity. This year has seen a continued expansion in the number of species represented, with >4800 eukaryotic and >31 300 prokaryotic genomes available. The new Ensembl site, currently in beta, has continued to develop, currently holding >2700 eukaryotic genome assemblies. The new site provides genome, gene, transcript, homology and variation views, and will replace the current Rapid Release site; this represents a key step towards provision of a single integrated Ensembl site. Additional activities have included developing improved regulatory annotation for human, mouse and agricultural species, and expanding the Ensembl Variant Effect Predictor tool. To learn more about Ensembl, help and documentation are available along with an extensive training program that can be accessed via our training pages.

Eimeria of chickens: the changing face of an old foe

RESEARCH HIGHLIGHTS

The cost of coccidiosis in chickens fluctuates considerably, peaking in 2022.Three new Eimeria species can infect chickens and escape current vaccines.Eimeria infection exerts wide-ranging effects on enteric microbiota.

Pangenome graphs and their applications in biodiversity genomics

Complete datasets of genetic variants are key to biodiversity genomic studies. Long-read sequencing technologies allow the routine assembly of highly contiguous, haplotype-resolved reference genomes. However, even when complete, reference genomes from a single individual may bias downstream analyses and fail to adequately represent genetic diversity within a population or species. Pangenome graphs assembled from aligned collections of high-quality genomes can overcome representation bias by integrating sequence information from multiple genomes from the same population, species or genus into a single reference. Here, we review the available tools and data structures to build, visualize and manipulate pangenome graphs while providing practical examples and discussing their applications in biodiversity and conservation genomics across the tree of life.

De Novo Genome Assembly for an Endangered Lemur Using Portable Nanopore Sequencing in Rural Madagascar

As one of the most threatened mammalian taxa, lemurs of Madagascar are facing unprecedented anthropogenic pressures. To address conservation imperatives such as this, researchers have increasingly relied on conservation genomics to identify populations of particular concern. However, many of these genomic approaches necessitate high-quality genomes. While the advent of next-generation sequencing technologies and the resulting reduction in associated costs have led to the proliferation of genomic data and high-quality reference genomes, global discrepancies in genomic sequencing capabilities often result in biological samples from biodiverse host countries being exported to facilities in the Global North, creating inequalities in access and training within genomic research. Here, we present the first published reference genome for the endangered red-fronted brown lemur (Eulemur rufifrons) from sequencing efforts conducted entirely within the host country using portable Oxford Nanopore sequencing. Using an archived E. rufifrons specimen, we conducted long-read, nanopore sequencing at the Centre ValBio Research Station near Ranomafana National Park, in rural Madagascar, generating over 750 Gb of sequencing data from 10 MinION flow cells. Exclusively using this long-read data, we assembled 2.157 gigabase, 2980-contig nuclear assembly with an N50 of 101.6 Mb and a 17,108 bp mitogenome. The nuclear assembly had 30× average coverage and was comparable in completeness to other primate reference genomes, with a 96.1% BUSCO completeness score for primate-specific genes. As the first published reference genome for E. rufifrons and the only annotated genome available for the speciose Eulemur genus, this resource will prove vital for conservation genomic studies while our efforts exhibit the potential of this protocol to address research inequalities and build genomic capacity.

The Ribosomal Operon Database: A Full-Length rDNA Operon Database Derived From Genome Assemblies

Current rDNA reference sequence databases are tailored towards shorter DNA markers, such as parts of the 16/18S marker or the internally transcribed spacer (ITS) region. However, due to advances in long-read DNA sequencing technologies, longer stretches of the rDNA operon are increasingly used in environmental sequencing studies to increase the phylogenetic resolution. There is, therefore, a growing need for longer rDNA reference sequences. Here, we present the ribosomal operon database (ROD), which includes eukaryotic full-length rDNA operons fished from publicly available genome assemblies. Full-length operons were detected in 34.1% of the 34,701 examined eukaryotic genome assemblies from NCBI. In most cases (53.1%), more than one operon variant was detected, which can be due to intragenomic operon copy variability, allelic variation in non-haploid genomes, or technical errors from the sequencing and assembly process. The highest copy number found was 5947 in Zea mays. In total, 453,697 unique operons were detected, with 69,480 operon variant clusters remaining after intragenomic clustering at 99% sequence identity. The operon length varied extensively across eukaryotes, ranging from 4136 to 16,463 bp, which will lead to considerable polymerase chain reaction (PCR) bias during amplification of the entire operon. Clustering the full-length operons revealed that the different parts (i.e., 18S, 28S, and the hypervariable regions V4 and V9 of 18S) provide divergent taxonomic resolution, with 18S, the V4 and V9 regions being the most conserved. The ROD will be updated regularly to provide an increasing number of full-length rDNA operons to the scientific community.

MCHelper automatically curates transposable element libraries across eukaryotic species

The number of species with high-quality genome sequences continues to increase, in part due to the scaling up of multiple large-scale biodiversity sequencing projects. While the need to annotate genic sequences in these genomes is widely acknowledged, the parallel need to annotate transposable element (TE) sequences that have been shown to alter genome architecture, rewire gene regulatory networks, and contribute to the evolution of host traits is becoming ever more evident. However, accurate genome-wide annotation of TE sequences is still technically challenging. Several de novo TE identification tools are now available, but manual curation of the libraries produced by these tools is needed to generate high-quality genome annotations. Manual curation is time-consuming, and thus impractical for large-scale genomic studies, and lacks reproducibility. In this work, we present the Manual Curator Helper tool MCHelper, which automates the TE library curation process. By leveraging MCHelper's fully automated mode with the outputs from three de novo TE identification tools, RepeatModeler2, EDTA, and REPET, in the fruit fly, rice, hooded crow, zebrafish, maize, and human, we show a substantial improvement in the quality of the TE libraries and genome annotations. MCHelper libraries are less redundant, with up to 65% reduction in the number of consensus sequences, have up to 11.4% fewer false positive sequences, and up to ∼48% fewer "unclassified/unknown" TE consensus sequences. Genome-wide TE annotations are also improved, including larger unfragmented insertions. Moreover, MCHelper is an easy-to-install and easy-to-use tool.

Long-read RNA sequencing can probe organelle genome pervasive transcription

40 years ago, organelle genomes were assumed to be streamlined and, perhaps, unexciting remnants of their prokaryotic past. However, the field of organelle genomics has exposed an unparallel diversity in genome architecture (i.e. genome size, structure, and content). The transcription of these eccentric genomes can be just as elaborate - organelle genomes are pervasively transcribed into a plethora of RNA types. However, while organelle protein-coding genes are known to produce polycistronic transcripts that undergo heavy posttranscriptional processing, the nature of organelle noncoding transcriptomes is still poorly resolved. Here, we review how wet-lab experiments and second-generation sequencing data (i.e. short reads) have been useful to determine certain types of organelle RNAs, particularly noncoding RNAs. We then explain how third-generation (long-read) RNA-Seq data represent the new frontier in organelle transcriptomics. We show that public repositories (e.g. NCBI SRA) already contain enough data for inter-phyla comparative studies and argue that organelle biologists can benefit from such data. We discuss the prospects of using publicly available sequencing data for organelle-focused studies and examine the challenges of such an approach. We highlight that the lack of a comprehensive database dedicated to organelle genomics/transcriptomics is a major impediment to the development of a field with implications in basic and applied science.

A Faroese perspective on decoding life for sustainable use of nature and protection of biodiversity

Biodiversity is under pressure, mainly due to human activities and climate change. At the international policy level, it is now recognised that genetic diversity is an important part of biodiversity. The availability of high-quality reference genomes gives the best basis for using genetics and genetic diversity towards the global aims of (1) the protection of species, biodiversity, and nature, and (2) the management of biodiversity for achieving sustainable harvesting of nature. Protecting biodiversity is a global responsibility, also resting on small nations, like the Faroe Islands. Being in the middle of the North Atlantic Ocean and having large fisheries activity, the nation has a particular responsibility towards maritime matters. We here provide the reasoning behind the Genome Atlas of Faroese Ecology (Gen@FarE), a project based on our participation in the European Reference Genome Atlas consortium (ERGA). Gen@FarE has three major aims: (1) To acquire high-quality genomes of all eukaryotic species in the Faroe Islands and Faroese waters. (2) To establish population genetics for species of commercial or ecological interest. (3) To establish an information databank for all Faroese species, combined with a citizen science registration database, making it possible for the public to participate in acquiring and maintaining the overview of Faroese species in both terrestrial and marine environments. Altogether, we believe that this will enhance the society's interest in and awareness of biodiversity, thereby protecting the foundations of our lives. Furthermore, the combination of a wide and highly competent ERGA umbrella and more targeted national projects will help fulfil the formal and moral responsibilities that all nations, also those with limited resources, have in protecting biodiversity and achieving sustainability in harvesting from nature.

The genome sequence of the Gold Spot moth, Plusia festucae (Linnaeus, 1758)

We present a genome assembly from an individual male specimen of Plusia festucae (Gold Spot; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence has a total length of 422.50 megabases. Most of the assembly (99.92%) is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.3 kilobases in length. Gene annotation of this assembly on Ensembl identified 19,273 protein-coding genes.

Microsatellites explorer: A database of short tandem repeats across genomes

Short tandem repeats (STRs) are widespread, repetitive elements, with a number of biological functions and are among the most rapidly mutating regions in the genome. Their distribution varies significantly between taxonomic groups in the tree of life and are highly polymorphic within the human population. Advances in sequencing technologies coupled with decreasing costs have enabled the generation of an ever-growing number of complete genomes. Additionally, the arrival of accurate long reads has facilitated the generation of Telomere-to-Telomere (T2T) assemblies of complete genomes. Nevertheless, there is no comprehensive database that encompasses the STRs found per genome across different organisms and for different human genomes across diverse ancestries. Here we introduce Microsatellites Explorer, a database of STRs found in the genomes of 117,253 organisms across all major taxonomic groups, 15 T2T genome assemblies of different organisms, and 94 human haplotypes from the human pangenome. The database currently hosts 406,758,798 STR sequences, serving as a centralized user-friendly repository to perform searches, interactive visualizations, and download existing STR data for independent analysis. Microsatellites Explorer is implemented as a web-portal for browsing, analyzing and downloading STR data. Microsatellites Explorer is publicly available at https://www.microsatellitesexplorer.com.

Three decades of rat genomics: approaching the finish(ed) line

The rat, Rattus norvegicus, has provided an important model for investigation of a range of characteristics of biomedical importance. Here we survey the origins of this species, its introduction into laboratory research, and the emergence of genetic and genomic methods that utilize this model organism. Genomic studies have yielded important progress and provided new insight into several biologically important traits. However, some studies have been impeded by the lack of a complete and accurate reference genome for this species. New sequencing and genome assembly methods applied to the rat have resulted in a new reference genome assembly, GRCr8, which is a near telomere-to-telomere assembly of high base-level accuracy that incorporates several elements not captured in prior assemblies. As genome assembly methods continue to advance and production costs become a less significant obstacle, genome assemblies for multiple inbred rat strains are emerging. These assemblies will allow a rat pangenome assembly to be constructed that captures all the genetic variations in strains selected for their utility in research and will overcome reference bias, a limitation associated with reliance on a single reference assembly. By this means, the full utility of this model organism to genomic studies will begin to be revealed.

A map of canine sequence variation relative to a Greenland wolf outgroup

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

MAFcounter: An efficient tool for counting the occurrences of k-mers in MAF files

Motivation

With the rapid expansion of large-scale biological datasets, DNA and protein sequence alignments have become essential for comparative genomics and proteomics. These alignments facilitate the exploration of sequence similarity patterns, providing valuable insights into sequence conservation, evolutionary relationships and for functional analyses. Typically, sequence alignments are stored in formats such as the Multiple Alignment Format (MAF). Counting k-mer occurrences is a crucial task in many computational biology applications, but currently, there is no algorithm designed for k-mer counting in alignment files.

Results

We have developed MAFcounter, the first k-mer counter dedicated to alignment files. MAFcounter is multithreaded, fast, and memory efficient, enabling k-mer counting in DNA and protein sequence alignment files.

Availability

The MAFcounter package and its Python bindings are released under GPL license as a multi-platform application and are available at: https://github.com/Georgakopoulos-Soares-lab/MAFcounter.

The genome sequence of the Poplar Grey moth, Subacronicta megacephala (Denis & Schiffermüller, 1775)

We present a genome assembly from an individual male Subacronicta megacephala (Poplar Grey moth; Arthropoda; Insecta; Lepidoptera; Noctuidae). The genome sequence has a total length of 424.20 megabases. Most of the assembly (99.02%) is scaffolded into 31 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genome has also been assembled and is 15.35 kilobases in length. Gene annotation of this assembly on Ensembl identified 18,189 protein-coding genes.

Extensive longevity and DNA virus-driven adaptation in nearctic Myotis bats

The genus Myotis is one of the largest clades of bats, and exhibits some of the most extreme variation in lifespans among mammals alongside unique adaptations to viral tolerance and immune defense. To study the evolution of longevity-associated traits and infectious disease, we generated near-complete genome assemblies and cell lines for 8 closely related species of Myotis. Using genome-wide screens of positive selection, analyses of structural variation, and functional experiments in primary cell lines, we identify new patterns of adaptation contributing to longevity, cancer resistance, and viral interactions in bats. We find that Myotis bats have some of the most significant variation in cancer risk across mammals and demonstrate a unique DNA damage response in primary cells of the long-lived M. lucifugus. We also find evidence of abundant adaptation in response to DNA viruses - but not RNA viruses - in Myotis and other bats in sharp contrast with other mammals, potentially contributing to the role of bats as reservoirs of zoonoses. Together, our results demonstrate how genomics and primary cells derived from diverse taxa uncover the molecular bases of extreme adaptations in non-model organisms.

The genome sequence of the whiskered bat, Myotis mystacinus (Kuhl, 1817)

We present a genome assembly from an individual male Myotis mystacinus (whiskered bat; Chordata; Mammalia; Chiroptera; Vespertilionidae). The genome sequence has a total length of 2,081.20 megabases. Most of the assembly (97.52%) is scaffolded into 23 chromosomal pseudomolecules, including the X and Y sex chromosomes. The mitochondrial genome has also been assembled and is 16.93 kilobases in length.

Harnessing genomic technologies for one health solutions in the tropics

BACKGROUND

The targeted application of cutting-edge high-throughput molecular data technologies provides an enormous opportunity to address key health, economic and environmental issues in the tropics within the One Health framework. The Earth's tropical regions are projected to contain > 50% of the world's population by 2050 coupled with 80% of its biodiversity however these regions are relatively less developed economically, with agricultural productivity substantially lower than temperate zones, a large percentage of its population having limited health care options and much of its biodiversity understudied and undescribed. The generation of high-throughput molecular data and bespoke bioinformatics capability to address these unique challenges offers an enormous opportunity for people living in the tropics. MAIN: In this review we discuss in depth solutions to challenges to populations living in tropical zones across three critical One Health areas: human health, biodiversity and food production. This review will examine how some of the challenges in the tropics can be addressed through the targeted application of advanced omics and bioinformatics and will discuss how local populations can embrace these technologies through strategic outreach and education ensuring the benefits of the One Health approach is fully realised through local engagement.

CONCLUSION

Within the context of the One Health framework, we will demonstrate how genomic technologies can be utilised to improve the overall quality of life for half the world's population.

The genome sequence of the long-horned nomad bee, Nomada hirtipes Pérez, 1884

We present a genome assembly from an individual female Nomada hirtipes (the long-horned nomad bee; Arthropoda; Insecta; Hymenoptera; Apidae). The genome sequence has a total length of 316.5 megabases. Most of the assembly (90.79%) is scaffolded into 16 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 29.88 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,693 protein-coding genes.

The genome sequence of the jet ant, Lasius fuliginosus (Latreille, 1798)

We present a genome assembly from an individual female Lasius fuliginosus (the jet ant; Arthropoda; Insecta; Hymenoptera; Formicidae). The genome sequence is 256.2 megabases in span. Most of the assembly is scaffolded into 14 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 18.75 kilobases in length.

The genome sequence of the Scalloped Hook-tip moth, Falcaria lacertinaria (Linnaeus, 1758)

We present a genome assembly from an individual female Falcaria lacertinaria (the Scalloped Hook-tip; Arthropoda; Insecta; Lepidoptera; Drepanidae). The genome sequence has a total length of 300.20 megabases. Most of the assembly is scaffolded into 32 chromosomal pseudomolecules, including the W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 16.07 kilobases in length. Gene annotation of this assembly on Ensembl identified 11,709 protein-coding genes.

Disentangling cobionts and contamination in long-read genomic data using sequence composition

The recent acceleration in genome sequencing targeting previously unexplored parts of the tree of life presents computational challenges. Samples collected from the wild often contain sequences from several organisms, including the target, its cobionts, and contaminants. Effective methods are therefore needed to separate sequences. Though advances in sequencing technology make this task easier, it remains difficult to taxonomically assign sequences from eukaryotic taxa that are not well represented in databases. Therefore, reference-based methods alone are insufficient. Here, I examine how we can take advantage of differences in sequence composition between organisms to identify symbionts, parasites, and contaminants in samples, with minimal reliance on reference data. To this end, I explore data from the Darwin Tree of Life project, including hundreds of high-quality HiFi read sets from insects. Visualizing two-dimensional representations of read tetranucleotide composition learned by a variational autoencoder can reveal distinct components of a sample. Annotating the embeddings with additional information, such as coding density, estimated coverage, or taxonomic labels allows rapid assessment of the contents of a dataset. The approach scales to millions of sequences, making it possible to explore unassembled read sets, even for large genomes. Combined with interactive visualization tools, it allows a large fraction of cobionts reported by reference-based screening to be identified. Crucially, it also facilitates retrieving genomes for which suitable reference data are absent.

Pictures at an exhibition: How to share your imaging data

Open access to data underpinning published results is a key pillar of scientific reproducibility. Making data available at scale also provides opportunities for data reuse, encouraging the development of new analysis approaches. In this poster article, accompanying a recorded talk, we will explain the benefits of publicly archiving your image data alongside your published manuscripts, as well as highlight what resources are available to do this. This will include the BioImage Archive, EMBL-EBI's new resource for biological image data, https://www.ebi.ac.uk/bioimage-archive/. We will look at how image data submission works, how to prepare in advance for archiving your data and upcoming developments.

TeloSearchLR: an algorithm to detect novel telomere repeat motifs using long sequencing reads

Telomeres are eukaryotic chromosome end structures that guard against sequence loss and aberrant chromosome fusions. Telomeric repeat motifs (TRMs), the minimal repeating unit of a telomere, vary from species to species, with some evolutionary clades experiencing a rapid sequence divergence. To explore the full scope of this evolutionary divergence, many bioinformatic tools have been developed to infer novel TRMs using repetitive sequence search on short sequencing reads. However, novel telomeric motifs remain unidentified in up to half of the sequencing libraries assayed with these tools. A possible reason may be that short reads, derived from extensively sheared DNA, preserve little to no positional context of the repetitive sequences assayed. On the other hand, if a sequencing read is sufficiently long, telomeric sequences must appear at either end rather than in the middle. The TeloSearchLR algorithm relies on this to help identify novel TRMs on long reads, in many cases where short-read search tools have failed. In addition, we demonstrate that TeloSearchLR can reveal unusually long telomeric motifs not maintained by telomerase, and it can also be used to anchor terminal scaffolds in new genome assemblies.

GENCODE: massively expanding the lncRNA catalog through capture long-read RNA sequencing

Accurate and complete gene annotations are indispensable for understanding how genome sequences encode biological functions. For twenty years, the GENCODE consortium has developed reference annotations for the human and mouse genomes, becoming a foundation for biomedical and genomics communities worldwide. Nevertheless, collections of important yet poorly-understood gene classes like long non-coding RNAs (lncRNAs) remain incomplete and scattered across multiple, uncoordinated catalogs, slowing down progress in the field. To address these issues, GENCODE has undertaken the most comprehensive lncRNAs annotation effort to date. This is founded on the manual annotation of full-length targeted long-read sequencing, on matched embryonic and adult tissues, of orthologous regions in human and mouse. Altogether 17,931 novel human genes (140,268 novel transcripts) and 22,784 novel mouse genes (136,169 novel transcripts) have been added to the GENCODE catalog representing a 2-fold and 6-fold increase in transcripts, respectively - the greatest increase since the sequencing of the human genome. Novel gene annotations display evolutionary constraints, have well-formed promoter regions, and link to phenotype-associated genetic variants. They greatly enhance the functional interpretability of the human genome, as they help explain millions of previously-mapped "orphan" omics measurements corresponding to transcription start sites, chromatin modifications and transcription factor binding sites. Crucially, our targeted design assigned human-mouse orthologs at a rate beyond previous studies, tripling the number of human disease-associated lncRNAs with mouse orthologs. The expanded and enhanced GENCODE lncRNA annotations mark a critical step towards deciphering the human and mouse genomes.

The genome sequence of the John Dory, Zeus faber Linnaeus, 1758

We present a genome assembly from an individual Zeus faber (the John Dory; Chordata; Actinopteri; Zeiformes; Zeidae). The genome sequence is 804.7 megabases in span. Most of the assembly is scaffolded into 22 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 16.72 kilobases in length.

The cryptonephridial/rectal complex: an evolutionary adaptation for water and ion conservation

Arthropods have integrated digestive and renal systems, which function to acquire and maintain homeostatically the substances they require for survival. The cryptonephridial complex (CNC) is an evolutionary novelty in which the renal organs and gut have been dramatically reorganised. Parts of the renal or Malpighian tubules (MpTs) form a close association with the surface of the rectum, and are surrounded by a novel tissue, the perinephric membrane, which acts to insulate the system from the haemolymph and thus allows tight regulation of ions and water into and out of the CNC. The CNC can reclaim water and solutes from the rectal contents and recycle these back into the haemolymph. Fluid flow in the MpTs runs counter to flow within the rectum. It is this countercurrent arrangement that underpins its powerful recycling capabilities, and represents one of the most efficient water conservation mechanisms in nature. CNCs appear to have evolved multiple times, and are present in some of the largest and most evolutionarily successful insect groups including the larvae of most Lepidoptera and in a major beetle lineage (Cucujiformia + Bostrichoidea), suggesting that the CNC is an important adaptation. Here we review the knowledge of this remarkable organ system gained over the past 200 years. We first focus on the CNCs of tenebrionid beetles, for which we have an in-depth understanding from physiological, structural and ultrastructural studies (primarily in Tenebrio molitor), which are now being extended by studies in Tribolium castaneum enabled by advances in molecular and microscopy approaches established for this species. These recent studies are beginning to illuminate CNC development, physiology and endocrine control. We then take a broader view of arthropod CNCs, phylogenetically mapping their reported occurrence to assess their distribution and likely evolutionary origins. We explore CNCs from an ecological viewpoint, put forward evidence that CNCs may primarily be adaptations for facing the challenges of larval life, and argue that their loss in many aquatic species could point to a primary function in conserving water in terrestrial species. Finally, by considering the functions of renal and digestive epithelia in insects lacking CNCs, as well as the typical architecture of these organs in relation to one another, we propose that ancestral features of these organs predispose them for the evolution of CNCs.

Taxonomy Identifiers (TaxId) for Biodiversity Genomics: a guide to getting TaxId for submission of data to public databases

Biodiversity genomics critically depends on correct taxonomic identification of the sample from which data are derived. Tracking of that taxonomic information through systems that archive data and report on genome sequencing efforts. For submission of data to the International Nucleotide Sequence Database Collaboration (INSDC) databases (DNA DataBank of Japan [DDBJ], European Nucleotide Archive [ENA] and National Center for Biotechnology Information [NCBI]), samples and data derived from them must be assigned a species-level NCBI Taxonomy taxonomic identifier (TaxId, sometimes referred to as taxId or txid). We thus need to be able to identify the TaxId for a target species efficiently. Because the NCBI Taxonomy does not include all known species and cannot preemptively represent unknown taxa, we also need an efficient process for generating new TaxIds for species not yet listed. This document provides workflows for different kinds of TaxId acquisition scenarios and was created to guide users in these processes. Although developed for European projects such as Darwin Tree of Life and the European Reference Genome Atlas, the workflows are universally applicable and describe the use of ENA in resolving taxonomic issues. Too Long: Didn't Read (TL;DR): Use the ENA REST API programmatically to retrieve TaxIds for target species and confirm that sequence data can be submitted to those TaxIds.Use the NCBI Web interface to NCBI Taxonomy to identify potential homotypic synonyms.Request a new TaxId from ENA for a species not yet in NCBI Taxonomy, and for species-like entries for which the full Linnaean binomen is not determined (see https://ena-docs.readthedocs.io/en/latest/faq/taxonomy_requests.html#creating-taxon-requests).Discuss directly with the NCBI Taxonomy curators or the curators at ENA and NCBI whenever you think there is an opportunity to improve their database.

MAFin: Motif Detection in Multiple Alignment Files

Motivation

Genome and Proteome Alignments, represented by the Multiple Alignment File (MAF) format, have become a standard approach in the field of comparative genomics and proteomics. However, current approaches lack a direct method for motif detection within MAF files. To address this gap, we present MAFin, a novel tool that enables efficient motif detection and conservation analysis in MAF files, streamlining genomic and proteomic research.

Results

We developed MAFin, the first motif detection tool for Multiple Alignment Format files. MAFin enables the multithreaded search of conserved motifs using three approaches: 1) by using user-specified k-mers to search the sequences. 2) with regular expressions, in which case one or more patterns are searched, and 3) with predefined Position Weight Matrices. Once the motif has been found, MAFin detects the motif instances and calculates the conservation across the aligned sequences. MAFin also calculates a conservation percentage, which provides information about the conservation levels of each motif across the aligned sequences, based on the number of matches relative to the length of the motif. A set of statistics enable the interpretation of each motif's conservation level, and the detected motifs are exported in JSON and CSV files for downstream analyses.

Availability

MAFin is released as a Python package under the GPL license as a multi-platform application and is available at: https://github.com/Georgakopoulos-Soares-lab/MAFin.