OrthoDB v11: annotation of orthologs in the widest sampling of organismal diversity

Convergent evolution in nuclear and mitochondrial OXPHOS subunits underlies the phylogenetic discordance in deep lineages of Squamata

The order Squamata is a good candidate for detecting unusual patterns of mitochondrial evolution. The lineages leading to the snake and agamid clades likely experienced convergent evolution in mitochondrial OXidative PHOSphorylation (OXPHOS) genes, which provides strong support for the sister relationship of these two groups. The OXPHOS subunits are encoded by both the nuclear and mitochondrial genomes, which are subject to distinct evolutionary pressures. Nevertheless, the cooperation between OXPHOS subunits is essential for proper OXPHOS function, as incompatibilities between subunits can be highly deleterious. In the present study, we annotated OXPHOS genes of 56 Squamata species. The nuclear OXPHOS subunits that physically interact with mitochondrial proteins also support the clade sister relationship between snakes and agamids. Additionally, we found a significant number of convergent amino acid changes between agamids and snakes, not only in mitochondrial OXPHOS genes, but also in nuclear ones, with a higher rate of convergence in the nuclear OXPHOS subunits that play central roles in the OXPHOS complexes, like COX4 and NDUFA4. Overall, the common selective pressures in two distinct lineages can lead two sets of genes, encoded by two different genomes, to exhibit similar patterns of convergent evolution, as well as similar evolutionary rates. As a consequence, the coevolution of interdependent subunits and their adaptation to specific evolutionary pressures can heavily influence the molecular structure of cytonuclear enzyme complexes and blur phylogenetic signals.

Multi-omics insights into the response of Aspergillus parasiticus to long-chain alkanes in relation to polyethylene modification

Plastic pollution presents a global challenge, with polyethylene (PE) being among the most persistent plastics due to its durability and environmental resilience. Long-chain alkane (lcAlk) degrading microbes are a potential source of PE-degrading enzymes, as both lcAlk and PE are large hydrophobic compounds that consist exclusively of C-C and C-H bonds. In this work, we employed a multi-omics approach to study the ability of Aspergillus parasiticus MM36, an isolate derived from Tenebrio molitor intestines, to metabolize lcAlk and secrete enzymes that are potentially capable of modifying PE. The fungus was grown with hexadecane (C16) or a mixture of lcAlk (C24 to C36) as carbon sources and culture supernatants were tested daily for their ability to modify PE. Proteomic analysis identified induced oxidases hypothetically involved in lcAlk and PE functionalization. Key enzymes include multicopper oxidases, peroxidases, an unspecific peroxygenase and FAD-dependent monooxygenases. Surfactant proteins facilitating enzymatic and cellular interaction with hydrophobic substrates, such as one hydrophobin, three hydrophobic surface-binding proteins (HsbA) and one cerato platanin, were present in all secretomes. Transcriptomic analysis comparing lcAlk to C16 cultures highlighted the enrichment of oxidoreductase activities and carboxylic acid metabolism in both lcAlk incubation days, with transmembrane transporters and transferases predominating on day 2 and biosynthetic processes on day 3. In C16 cultures, hydrolytic enzymes, including esterases, were upregulated alongside Baeyer-Villiger monooxygenases, suggesting a shift toward sub-terminal hydroxylation. Integrating transcriptomic and secretomic data, we propose a mechanism for lcAlk assimilation by A. parasiticus MM36, involving extracellular oxyfunctionalization, hydrocarbon uptake via surface-modifying proteins and channeling through membrane transporters for energy consumption and biosynthetic processes. This study provides insights into fungal mechanisms for alkane metabolism and highlights their potential relevance to plastic biotransformation.

NMD-mediated posttranscriptional regulation fine-tunes the NLR-WRKY regulatory module to modulate bacterial defense response

Nonsense-mediated mRNA decay (NMD) is a conserved eukaryotic surveillance system that maintains transcriptome integrity by degrading aberrant RNA transcripts. NMD ensures proper growth and development by preventing autoimmunity through the direct regulation of nucleotide-binding, leucine-rich repeat (NLR) genes. Whether NMD directly regulates WRKY genes remains unclear, despite their upregulation in NMD-deficient plants, and potential feedback between NLRs and WRKYs is also poorly understood. In this study, we showed that NMD also directly regulates a subset of WRKY (WRKY15, 18, 25, 33, 46, 60, and 70) genes, particularly at lower temperatures (16°C). NMD signature-containing transcripts of WRKY46 and WRKY70, selected as representative NMD-regulated WRKY genes, showed increased half-lives in NMD-deficient mutants. Transcriptome analyses showed that these seven NMD-regulated WRKY genes are induced in response to bacterial infection. Potential homologues of these seven NMD-regulated WRKY genes in maize and rice showed similar induction in response to bacterial pathogen infection. Furthermore, these NMD-regulated WRKY genes are induced in plants overexpressing RESISTANT TO P. SYRINGAE 4 (RPS4) in a temperature-dependent manner. By using ChIP-seq and DAP-seq data of WRKY transcription factors, we showed that WRKYs potentially regulate a significant number of NLR genes by directly binding to the W-box in their promoter regions. Taken together, our findings revealed that in addition to the NLRs, the NMD machinery also regulates WRKY genes to keep the basal defense levels in check and the WRKY transcription factors directly regulate NLR genes to constitutes a positive feedback regulatory loop to optimize the plant response to invading pathogens.

3D Genome Constrains Breakpoints of Inversions That Can Act as Barriers to Gene Flow in the Stickleback

DNA within the nucleus is organised into a well-regulated three-dimensional (3D) structure. However, how such 3D genome structures influence speciation processes remains largely elusive. Recent studies have shown that 3D genome structures influence mutation rates, including the occurrence of chromosomal rearrangement. For example, breakpoints of chromosomal rearrangements tend to be located at topologically associating domain (TAD) boundaries. Here, we hypothesised that TAD structures may constrain the location of chromosomal inversions and thereby shape the genomic landscape of divergence between species with ongoing gene flow, given that inversions can act as barriers to gene flow. To test this hypothesis, we used a pair of Japanese stickleback species, Gasterosteus nipponicus (Japan Sea stickleback) and G. aculeatus (three-spined stickleback). We first constructed chromosome-scale genome assemblies of both species using high fidelity long reads and high-resolution proximity ligation data and identified several chromosomal inversions. Second, via population genomic analyses, we revealed higher genetic differentiation in inverted regions than in colinear regions and no gene flow within inversions, which contrasts with the significant gene flow in colinear regions. Third, using Hi-C data, we revealed 3D genome structures of sticklebacks, delineated by A/B compartments and TADs. Finally, we found that inversion breakpoints tend to be located at TAD boundaries. Thus, our study demonstrates that the 3D genome constrains breakpoints of inversions that can act as barriers to gene flow in the stickleback. Further integration of 3D genome analyses with population genomics could provide novel insights into how the 3D genome influences speciation.

Chromosome-scale whole genome assembly and annotation of the Jamaican field cricket Gryllus assimilis

Gryllus assimilis, commonly known as Jamaican field cricket, is an edible insect with significant economic value in sustainable food production. Despite its importance, a high-quality reference genome of G. assimilis has not yet been published. Here, we report a chromosome-level reference genome of G. assimilis based on Oxford Nanopore Technologies (ONT) sequencing, Illumina sequencing, and Hi-C technologies. The assembled genome has a total length of 1.60 Gbp with a scaffold N50 of 102 Mbp, and 96.80% of the nucleotides was assigned to 15 chromosome-scale scaffolds. The assembly completeness was validated using BUSCO, achieving 99.5% completeness against the arthropoda database. We predicted 27,645 protein-coding genes, and 825 Mb repetitive elements were annotated in the reference genome. This reference genome of G. assimilis can provide a basis for the subsequent development of genomic resources, offering insights for future functional genomic studies, comparative genomics, and DNA-informed breeding of this species.

Predicting Protein Function in the AI and Big Data Era

It is an exciting time for researchers working to link proteins to their functions. Most techniques for extracting functional information from genomic sequences were developed several years ago, with major progress driven by the availability of big data. Now, groundbreaking advances in deep-learning and AI-based methods have enriched protein databases with three-dimensional information and offer the potential to predict biochemical properties and biomolecular interactions, providing key functional insights. This progress is expected to increase the proportion of functionally bright proteins in databases and deepen our understanding of life at the molecular level.

DIGITtally-a new tool for streamlining and simplifying Drosophila melanogaster meta-analysis

Drosophila melanogaster has one of the deepest research bases within the life sciences, with a wealth of high-quality tissue- and cell type-specific transcriptomic data available. However, integrating large datasets derived from disparate sources is not trivial. We have designed a broadly applicable solution to this problem in the form of the Drosophila Interesting Genes in Individual Tissues-tally (DIGITtally) system. It is freely available online at www.digittally.org. DIGITtally is customizable and hypothesis-free, allowing meta-analysis across the Drosophila research space along with analysis of conservation in other species, querying 10 data sources for seven indicators of tissue-specific activity. We have applied DIGITtally to a pertinent question within entomology-that is, whether a specific pattern of gene expression underlies the transporting activity of epithelial tissues (an 'epitheliome'). By using DIGITtally to survey gene expression throughout the tissues comprising the D. melanogaster alimentary canal (salivary gland, midgut, Malpighian tubules, and hindgut), we have verified the existence of a specific 'epithelial' vacuolar-type ATPase configuration.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Current genome sequencing initiatives across a wide range of life forms offer significant potential to enhance our understanding of evolutionary relationships and support transformative biological and medical applications. Species trees play a central role in many of these applications; however, despite the widespread availability of genome assemblies, accurate inference of species trees remains challenging due to the limited automation, substantial domain expertise, and computational resources required by conventional methods. To address this limitation, we present ROADIES, a fully automated pipeline to infer species trees starting from raw genome assemblies. In contrast to the prominent approach, ROADIES incorporates a unique strategy of randomly sampling segments of the input genomes to generate gene trees. This eliminates the need for predefining a set of loci, limiting the analyses to a fixed number of genes, and performing the cumbersome gene annotation and/or whole genome alignment steps. ROADIES also eliminates the need to infer orthology by leveraging existing discordance-aware methods that allow multicopy genes. Using the genomic datasets from large-scale sequencing efforts across four diverse life forms (placental mammals, pomace flies, birds, and budding yeasts), we show that ROADIES infers species trees that are comparable in quality to the state-of-the-art studies but in a fraction of the time and effort, including on challenging datasets with rampant gene tree discordance and complex polyploidy. With its speed, accuracy, and automation, ROADIES has the potential to vastly simplify species tree inference, making it accessible to a broader range of scientists and applications.

Production of the light-activated elsinochrome phytotoxin in the soybean pathogen Coniothyrium glycines hints at virulence factor

The Dothideomycete pathogen Coniothyrium glycines causes red leaf blotch of soybean, a major disease in Africa. It is one of two fungal plant pathogens on the USDA PPQ Select Agents and Toxins list of pathogens important to the biosecurity of the United States, reflective of its potential to be highly destructive if introduced. Despite its importance, there are no published reports regarding the molecular basis of host infection. Examination of the C. glycines genome revealed a secondary metabolite gene cluster that is similar in gene content and organization to clusters that synthesize light-activated perylenequinone toxins, such as cercosporin. Perylenequinones are non-host specific toxins that, upon exposure to light, generate reactive oxygen species, which have near-universal toxicity to plant hosts. Coniothyrium glycines isolates from eastern and southern Africa were cultured axenically under light and dark conditions. Light-grown cultures produced red-pink pigmentation typical of perylenequinones. Differential gene expression analysis showed that six of the eight genes in the biosynthetic gene cluster, including the polyketide synthase gene, were significantly upregulated in light. Liquid chromatography-mass spectrometry confirmed production of the perylenequinone elsinochrome A, a known virulence factor in other fungal pathogens. On leaves incubated in the dark, significantly fewer lesions formed and symptoms were delayed, compared to leaves incubated in the light. In addition, we identified orthologous gene clusters in more distantly related Dothideomycete plant pathogens where their presence was previously unknown, indicating a broader importance of these toxins to agriculture and fungal ecology. This work provides the first evidence that elsinochrome A may contribute to the virulence of C. glycines.

The draft genome sequence of Diaporthe vaccinii, isolated from diseased cranberries

We report the assembly and annotation of the nuclear genome of Diaporthe vaccinii, a pathogenic fungus isolated from diseased cranberries in Quebec, Canada. The genome was sequenced with the Illumina paired-end sequencing technology, assembled into 67 Mbp across 588 contigs, with an N50 of 386 Kbp and 97.5% BUSCO completeness.

Draft genome of the endemic alpine ground beetle Carabus (Platycarabus) depressus (Coleoptera: Carabidae) from long-read sequencing of a frozen archived specimen

The rapid advancement of genomic technologies has enabled the production of highly contiguous reference genomes for nonmodel organisms. However, these methods often require exceptionally fresh material containing unfragmented high-molecular-weight nucleic acids. Researchers who preserve field-collected specimens in ethanol at ambient temperatures, prior to transferring them to long-term frozen archives, face challenges in applying advanced genomic approaches due to DNA and RNA fragmentation under suboptimal preservation conditions. To explore the potential of such preserved specimens as sources of reference genomes, we utilized Nanopore MinION technology to generate genomic data from a frozen archived specimen of the endemic alpine ground beetle Carabus (Platycarabus) depressus. Using a rapid in-house protocol for high-molecular-weight DNA extraction, followed by sequencing on a single flow cell, we produced 8.75 million raw reads with an N50 of 2.8 kb. The resulting assembly achieved remarkable completeness, recovering up to 98% of Benchmarking Universal Single-Copy Orthologs genes, despite a moderate N50 of 945 kb. This genome is only the second available for the taxonomically diverse genus Carabus, demonstrating the feasibility of using short-to-long-read sequencing on frozen archived specimens commonly housed in natural history collections. These findings open new avenues for advancing nonmodel organism genomics and its downstream applications.

A high-quality Oxford Nanopore assembly of the hourglass dolphin (Lagenorhynchus cruciger) genome

The hourglass dolphin (Lagenorhynchus cruciger) is a small cetacean species of the Southern Ocean, with significance to iwi Māori (Māori tribes) of Aotearoa New Zealand as taonga (treasured/valued). Due to the remoteness and difficulty of surveying Antarctic waters, it remains one of the least-studied dolphin species. A recent stranding of an hourglass dolphin represented a rare opportunity to generate a genome assembly as a resource for future study into the conservation and evolutionary biology of this species. In this study, we present a high-quality genome assembly of an hourglass dolphin individual using a single sequencing platform, Oxford Nanopore Technologies, coupled with computationally efficient assembly methods. Our assembly strategy yielded a genome of high contiguity (N50 of 8.07 Mbp) and quality (98.3% BUSCO completeness). Compared to other Delphinoidea reference genomes, this assembly has fewer missing BUSCOs than any except Orcinus orca, more single-copy complete BUSCOs than any except Phocoena sinus, and 20% fewer duplicated BUSCOs than the average Delphinoidea reference genome. This suggests that it is one of the most complete and accurate marine mammal genomes to date. This study showcases the feasibility of a cost-effective mammalian genome assembly method, allowing for genomic data generation outside the traditional confines of academia and/or resource-rich genome assembly hubs, and facilitating the ability to uphold Indigenous data sovereignty. In the future, the genome assembly presented here will allow valuable insights into the past population size changes, adaptation, vulnerability to future climate change of the hourglass dolphin and related species.

Chromatin loops are an ancestral hallmark of the animal regulatory genome

In bilaterian animals, gene regulation is shaped by a combination of linear and spatial regulatory information. Regulatory elements along the genome are integrated into gene regulatory landscapes through chromatin compartmentalization1,2, insulation of neighbouring genomic regions3,4 and chromatin looping that brings together distal cis-regulatory sequences5. However, the evolution of these regulatory features is unknown because the three-dimensional genome architecture of most animal lineages remains unexplored6,7. To trace the evolutionary origins of animal genome regulation, here we characterized the physical organization of the genome in non-bilaterian animals (sponges, ctenophores, placozoans and cnidarians)8,9 and their closest unicellular relatives (ichthyosporeans, filastereans and choanoflagellates)10 by combining high-resolution chromosome conformation capture11,12 with epigenomic marks and gene expression data. Our comparative analysis showed that chromatin looping is a conserved feature of genome architecture in ctenophores, placozoans and cnidarians. These sequence-determined distal contacts involve both promoter-enhancer and promoter-promoter interactions. By contrast, chromatin loops are absent in the unicellular relatives of animals. Our findings indicate that spatial genome regulation emerged early in animal evolution. This evolutionary innovation introduced regulatory complexity, ultimately facilitating the diversification of animal developmental programmes and cell type repertoires.

Redistribution of Ancestral Functions Underlies the Evolution of Venom Production in Marine Predatory Snails

Venom-secreting glands are highly specialized organs evolved throughout the animal kingdom to synthetize and secrete toxins for predation and defense. Venom is extensively studied for its toxin components and application potential; yet, how animals become venomous remains poorly understood. Venom systems therefore offer a unique opportunity to understand the molecular mechanisms underlying functional innovation. Here, we conducted a multispecies multi-tissue comparative transcriptomics analysis of 12 marine predatory gastropod species, including species with venom glands and species with homologous non-venom-producing glands, to examine how specialized functions evolve through gene expression changes. We found that while the venom gland specialized for the mass production of toxins, its homologous glands retained the ancestral digestive functions. The functional divergence and specialization of the venom gland were achieved through a redistribution of its ancestral digestive functions to other organs, specifically the esophagus. This entailed concerted expression changes and accelerated transcriptome evolution across the entire digestive system. The increase in venom gland secretory capacity was achieved through the modulation of an ancient secretory machinery, particularly genes involved in endoplasmic reticulum stress and unfolded protein response. This study shifts the focus from the well-explored evolution of toxins to the lesser-known evolution of the organ and mechanisms responsible for venom production. As such, it contributes to elucidating the molecular mechanisms underlying organ evolution at a fine evolutionary scale, highlighting the specific events that lead to functional divergence.

The First De Novo HiFi Genome Assemblies for Three Clownfish-hosting Sea Anemone Species (Anthozoa: Actiniaria)

The symbiosis between clownfish and giant tropical sea anemones (Order Actiniaria) is one of the most iconic on the planet. Distributed on tropical reefs, 28 species of clownfishes form obligate mutualistic relationships with 10 nominal species of venomous sea anemones. Our understanding of the symbiosis is limited by the fact that most research has been focused on the clownfishes. Chromosome-scale reference genomes are available for all clownfish species, yet only short reads-based reference genomes are available for five species of host sea anemones. Recent studies have shown that the clownfish-hosting sea anemones belong to three distinct clades of sea anemones that have evolved symbiosis with clownfishes independently. Here we present the first high-quality long-read assemblies for three species of clownfish-hosting sea anemones belonging to each of these clades: Entacmaea quadricolor, Stichodactyla haddoni, and Radianthus doreensis. PacBio HiFi sequencing yielded 1,597,562, 3,101,773, and 1,918,148 million reads for E. quadricolor, S. haddoni, and R. doreensis, respectively. All three assemblies were highly contiguous and complete with N50 values above 4 Mb and BUSCO completeness above 95% on the Metazoa dataset. Genome structural annotation with BRAKER3 predicted 20,454, 18,948, and 17,056 protein-coding genes in E. quadricolor, S. haddoni, and R. doreensis genome, respectively. These new resources will form the basis of comparative genomic analyses that will allow us to deepen our understanding of this mutualism from the host perspective.

Multiple Horizontal Transfers of Immune Genes Between Distantly Related Teleost Fishes

Horizontal gene transfer (HGT) is less frequent in eukaryotes than in prokaryotes, yet can have strong functional implications and was proposed as a causal factor for major adaptations in several eukaryotic lineages. Most cases of eukaryote HGT reported to date are inter-domain transfers, and few studies have investigated eukaryote-to-eukaryote HGTs. Here, we performed a large-scale survey of HGT among 242 species of ray-finned fishes. We found multiple lines of evidence supporting 19 teleost-to-teleost HGT events that involve 17 different genes in 11 teleost fish orders. The genes involved in these transfers show lower synonymous divergence than expected under vertical transmission, their phylogeny is inconsistent with that of teleost fishes, and they occur at non-syntenic positions in donor and recipient lineages. The distribution of HGT events in the teleost tree is heterogenous, with 8 of the 19 transfers occurring between the same two orders (Osmeriformes and Clupeiformes). Though we favor a scenario involving multiple HGT events, future work should evaluate whether hybridization between species belonging to different teleost orders may generate HGT-like patterns. Besides the previously reported transfer of an antifreeze protein, most transferred genes play roles in immunity or are pore-forming proteins, suggesting that such genes may be more likely than others to confer a strong selective advantage to the recipient species. Overall, our work shows that teleost-to-teleost HGT has occurred on multiple occasions, and it will be worth further quantifying these transfers and evaluating their impact on teleost evolution as more genomes are sequenced.

Chromosome-Level Genome Assembly of the Asian Tramp Snail Bradybaena similaris (Stylommatophora: Camaenidae)

While terrestrial land snails have long been used for evolutionary research, a lack of high-quality genomic resources has impeded recent progress. Bradybaena snails in particular have numerous intriguing traits that make them a good model for studying evolution, including shell pattern polymorphism and convergent evolution. They are also introduced and invasive across the world. In this study, we present a chromosome-level genome assembly of the Asian tramp snail Bradybaena similaris, utilizing 88-fold Illumina short-read sequences, 125-fold Nanopore long-read sequences, 63-fold PacBio HiFi sequences, and 47-fold Hi-C sequences. The assembled genome of 2.18 Gb is anchored to 28 chromosomes and exhibits high completeness (single copy, 91.7%; duplicates, 7.1%) and contiguity (N50 of 75.6 Mb). Additionally, we also obtained a high-quality transcriptome for annotation. This resource represents the first chromosome-level assembly for snails in the superfamily Helicoidea, which includes more than 5,000 species of terrestrial snails, and will facilitate genomic study in Bradybaena and, more broadly, in the superfamily Helicoidea.

Phylogenomic Inference Suggests Differential Deep Time Phylogenetic Signals from Nuclear and Organellar Genomes in Gymnosperms

The living gymnosperms include about 1200 species in five major groups: cycads, ginkgo, gnetophytes, Pinaceae (conifers I), and cupressophytes (conifers II). Molecular phylogenetic studies have yet to reach a unanimously agreed-upon relationship among them. Moreover, cytonuclear phylogenetic incongruence has been repeatedly observed in gymnosperms. We collated a comprehensive dataset from available genomes of 17 gymnosperms across the five major groups and added our own high-quality assembly of a species from Podocarpaceae (the second largest conifer family) to increase sampling width. We used these data to infer reconciled nuclear species phylogenies using two separate methods to ensure the robustness of our conclusions. We also reconstructed organelle phylogenomic trees from 42 mitochondrial and 82 plastid genes from 38 and 289 gymnosperm species across the five major groups, respectively. Our nuclear phylogeny consistently recovers the Ginkgo-cycads clade as the first lineage split from other gymnosperm clades and the Pinaceae as sister to gnetophytes (the Gnepines hypothesis). In contrast, the mitochondrial tree places cycads as the earliest lineage in gymnosperms and gnetophytes as sister to cupressophytes (the Gnecup hypothesis) while the plastomic tree supports the Ginkgo-cycads clade and gnetophytes as the sister to cupressophytes. We also examined the effect of mitochondrial RNA editing sites on the gymnosperm phylogeny by manipulating the nucleotide and amino acid sequences at these sites. Only complete removal of editing sites has an effect on phylogenetic inference, leading to a closer congruence between mitogenomic and nuclear phylogenies. This suggests that RNA editing sites carry a phylogenetic signal with distinct evolutionary traits.

Genome-wide association study uncovers key genomic regions governing agro-morphological and quality traits in Indian mustard [Brassica juncea (L.) Czern. and Coss.]

In Indian mustard, improving agro-morphological and quality traits through conventional methods are both cumbersome and resource-intensive. Marker-aided breeding presents a promising solution to these challenges. Hence, the present research aimed to identify genomic regions governing agro-morphological and quality traits using genome-wide association studies (GWAS). The GWAS panel comprised 142 diverse genotypes of Indian mustard were evaluated for 20 different agro-morphological and quality traits, revealing significant difference among genotypes. Subsequently, the GWAS panel genotyped using the Brassica 90K SNP array (Illumina). Structure and diversity analysis grouped the GWAS panel into 3 sub-populations or groups, and LD decay of 1.05 Mb was confirmed through genotypic analysis. GWAS using the BLINK model revealed a total of 49 marker-trait associations (MTAs), in which 28 and 21 MTAs were observed during rabi 2020-21 and rabi 2021-22 for various agro-morphological and quality traits, respectively. Amongst them, twelve MTAs demonstrated stable associations with the studied traits, including days to 50% flowering (DF), days to 100% flower termination (DFT), days to maturity (DM), plant height (PH), main shoot length (MSL), siliqua length (SL), seeds per siliqua (SPS), oil content (OC), and glucosinolates content (Glu) in both years. Moreover, in silico analysis of nearby regions of these stable SNPs revealed their association with 31 candidate genes known to be involved in various molecular, physiological, and biochemical pathways relevant to the studied traits. These genes can be further characterized and deciphered for more precise utilization in breeding programs in the future.

A chromosome-level genome assembly of the polyembryonic endoparasitoid Macrocentrus cingulum

Polyembryony is a unique reproductive pattern, where multiple genetically identical offspring develop from a single egg. Macrocentrus cingulum, a polyembryonic endoparasitoid, is one of the dominant parasitoids of the globally agricultural pest, Ostrinia furnacalis, and it can serve as a valuable model for investigating the mechanisms of polyembryony. However, the previously published genome of M. cingulum remains highly fragmented, limiting comprehensive studies of its biological characteristics. Here, we present a chromosome-level genome assembly of M. cingulum using PacBio HiFi and Hi-C sequencing. This genome assembly is approximately 151.93 Mb with a contig N50 of 8.42 Mb and a scaffold N50 of 16.93 Mb, organized into 9 chromosomes. The repeat sequences constitute 25.52% of the genome assembly. A total of 14,471 protein-coding genes with 98.2% BUSCO completeness were predicted, of which 12,500 genes have been annotated in public biological databases. In conclusion, this reported genome should be a valuable genomic resource for exploring macroevolutionary mechanisms underlying polyembryony, and the molecular mechanisms of polyembryony in M. cingulum.

Spiral phyllotaxis predicts left-right asymmetric growth and style deflection in mirror-image flowers of Cyanella alba

Many animals and plants show left-right (LR) asymmetry. The LR asymmetry of mirror-image flowers has clear functional significance, with the reciprocal placement of male and female organs in left- versus right-handed flowers promoting cross-pollination. Here, we study how handedness of mirror-image flowers is determined and elaborated during development in the South African geophyte Cyanella alba. Inflorescences of C. alba produce flowers with a largely consistent handedness. However, this handedness has no simple genetic basis and individual plants can switch their predominant handedness between years. Rather, it is the direction of the phyllotactic spiral that predicts floral handedness. Style deflection is driven by increased cell expansion in the adaxial carpel facing the next oldest flower compared to the other adaxial carpel. The more expanding carpel shows transcriptional signatures of increased auxin signaling and auxin application can reverse the orientation of style deflection. We propose that a recently described inherent LR auxin asymmetry in the initiating organs of spiral phyllotaxis determines handedness in C. alba, creating a stable yet non-genetic floral polymorphism. This mechanism links chirality across different levels of plant development and exploits a developmental constraint in a core patterning process to produce morphological variation of ecological relevance.

Decoding Specificity of Cyanobacterial MysDs in Mycosporine-like Amino Acid Biosynthesis through Heterologous Expression in Saccharomyces cerevisiae

Mycosporine-like amino acids (MAAs) are potent natural UV-protectants, but their industrial production is hindered by efficiency and sustainability issues in large-scale extraction of native hosts. Heterologous biosynthesis in chassis organisms provides a promising alternative route, although the substrate promiscuity of the ATP-grasp ligase MysD complicates the biosynthesis of specific MAAs. In this study, we developed a Saccharomyces cerevisiae strain with enhanced capacity of producing mycosporine-glycine (MG). This strain serves as an efficient MysD expression platform that converts MG into shinorine and porphyra-334. Through structural modeling, site-directed mutagenesis, and mutant characterization, we identified two residues on the omega-loop of MysD involved in determining product specificity. We further characterized the product specificity of 20 MysDs from diverse cyanobacterial lineages and confirmed the residue pattern-product specificity correlation. Our findings provide guidance for screening, selecting, and designing novel MysDs for industrial-scale MAA production through heterologous expression.

CMP-sialic acid synthetase in Drosophila requires N-glycosylation of a noncanonical site

Sialylation plays important roles in animals, affecting numerous molecular and cell interactions. In Drosophila, sialylation regulates neural transmission and mediates communication between neurons and glia. Drosophila CMP-sialic acid synthetase (CSAS), a key enzyme of the sialylation pathway, is localized to the Golgi and modified by N-glycosylation, suggesting that this modification can affect CSAS function. Here, we tested this hypothesis using in vitro and in vivo approaches. We found that CSAS proteins from divergent Drosophila species have two conserved N-glycosylation sites, including the rarely glycosylated noncanonical N-X-C sequon. We investigated CSAS glycosylation by generating CSAS "glycomutants" lacking glycosylation sites and analyzing them in vivo in transgenic rescue assays. The removal of noncanonical glycosylation significantly decreased CSAS activity, while the canonical site mutation did not affect CSAS function. Although all glycomutants were similarly localized to the Golgi, the non-canonical glycosylation, unlike the canonical one, affected CSAS stability in vivo and in vitro. Our results suggested that CSAS functions as a dimer, which was also supported by protein structure predictions that produced a dimer recapitulating the crystal structures of mammalian and bacterial counterparts, highlighting the evolutionary conservation of the CSAS structure-function relationship. This conclusion was supported by the rescue of CSAS mutants using the human ortholog. The noncanonical CSAS glycosylation was discussed in terms of a potential mechanism of temperature-dependent regulation of sialylation in poikilotherms that modulates neural activity in heat shock conditions. Taken together, we uncovered an important regulation of sialylation in Drosophila, highlighting a novel interplay between glycosylation pathways in neural regulation.

An S-acylated N-terminus and a conserved loop regulate the activity of the ABHD17 deacylase

The dynamic addition and removal of long-chain fatty acids modulate protein function and localization. The α/β hydrolase domain-containing (ABHD) 17 enzymes remove acyl chains from membrane-localized proteins such as the oncoprotein NRas, but how the ABHD17 proteins are regulated is unknown. Here, we used cell-based studies and molecular dynamics simulations to show that ABHD17 activity is controlled by two mobile elements-an S-acylated N-terminal helix and a loop-that flank the putative substrate-binding pocket. Multiple S-acylation events anchor the N-terminal helix in the membrane, enabling hydrophobic residues in the loop to engage with the bilayer. This stabilizes the conformation of both helix and loop, alters the conformation of the binding pocket, and optimally positions the enzyme for substrate engagement. S-acylation may be a general feature of acyl-protein thioesterases. By providing a mechanistic understanding of how the lipid modification of a lipid-removing enzyme promotes its enzymatic activity, this work contributes to our understanding of cellular S-acylation cycles.

Chromosome-scale and haplotype-resolved genome assembly of Populus trichocarpa

Populus trichocarpa, a pivotal model organism for woody transgenic research, not only garners substantial scientific interest but plays an integral role in forestry economics. Previous genomic assemblies of P. trichocarpa predominantly treated its heterozygous genome as homozygous, thereby neglecting crucial haplotypic diversity. Leveraging the high-fidelity (HiFi) sequencing capabilities of PacBio sequencing and the chromosome conformation capture insights provided by Illumina's Hi-C technique, this study is the first to achieve a near telomere-to-telomere assembly of both paternal and maternal haplotypes in P. trichocarpa. Comparative genomic analysis between these haplotypes has uncovered several allelic variants and pathways critical for trait determination through allele-specific expression. Furthermore, utilizing RNA-seq data from multiple tissues, this investigation has detailed the tissue-specific expression patterns of the leucine-rich repeat gene family, which are essential in mediating plant signal transduction and developmental regulation. Our results not only illuminate the functional genomics landscape of P. trichocarpa but also provide invaluable theoretical underpinnings for the genetic improvement of woody plants and a robust framework for exploring genetic variability and allelic expression disparities in arboreal species.

Structural insights into spliceosome fidelity: DHX35-GPATCH1- mediated rejection of aberrant splicing substrates

The spliceosome, a highly dynamic macromolecular assembly, catalyzes the precise removal of introns from pre-mRNAs. Recent studies have provided comprehensive structural insights into the step-wise assembly, catalytic splicing and final disassembly of the spliceosome. However, the molecular details of how the spliceosome recognizes and rejects suboptimal splicing substrates remained unclear. Here, we show cryo-electron microscopy structures of spliceosomal quality control complexes from a thermophilic eukaryote, Chaetomium thermophilum. The spliceosomes, henceforth termed B*Q, are stalled at a catalytically activated state but prior to the first splicing reaction due to an aberrant 5' splice site conformation. This state is recognized by G-patch protein GPATCH1, which is docked onto PRP8-EN and -RH domains and has recruited the cognate DHX35 helicase to its U2 snRNA substrate. In B*Q, DHX35 has dissociated the U2/branch site helix, while the disassembly helicase DHX15 is docked close to its U6 RNA 3'-end substrate. Our work thus provides mechanistic insights into the concerted action of two spliceosomal helicases in maintaining splicing fidelity by priming spliceosomes that are bound to aberrant splice substrates for disassembly.

Chromosome-level genome assembly of an endoparasitoid Cotesia ruficrus

Cotesia ruficrus is a gregarious larval endoparasitoid of many important agricultural pest moths, such as the fall armyworm, northern armyworm and rice leaffolder. While eight of Cotesia genomes had been reported, high-quality genomic data of C. ruficrus has yet to be performed. This study presents a chromosome-scale genome assembly of C. ruficrus, utilizing both PacBio HiFi-reads and Hi-C reads. The genome size is about 185.3 Mb harboring 10 chromosomes, and 209 unanchored scaffolds with a scaffold N50 length of 15.81 Mb. The anchored rate of sequences to chromosomes is 95.55%. A total of 14,001 genes have been annotated. This work delivers a first high-quality, chromosome-scale genome for C. ruficrus, enriching the datasets for comparative evolutionary studies within the Cotesia genus. Additionally, it offers a valuable genomic resource for future biological and genetic research, as well as for enhancing pest management strategies, such as biological control.

Evolutionary history of the DNA repair protein, Ku, in eukaryotes and prokaryotes

Ku is essential in non-homologous end-joining (NHEJ) across prokaryotes and eukaryotes, primarily in double-stranded breaks (DSBs) repair. It often presents as a multi-domain protein in eukaryotes, unlike their prokaryotic single-domain homologs. We systematically searched for Ku proteins across different domains of life. To elucidate the evolutionary history of the Ku protein, we constructed a maximum likelihood phylogenetic tree using Ku protein sequences from 100 representative eukaryotic, prokaryotic, and viral species. The resulting tree revealed a common node for eukaryotic Ku proteins, while viral and prokaryotic species clustered into a distinct clade. Our phylogenetic analysis reveals that the common ancestry of Ku70 and Ku80 likely resulted from a gene duplication event in the ancestral eukaryote. This inference is supported by BLASTp results, which indicate a close resemblance between archaeal Ku and eukaryotic Ku, particularly Ku70. The presence of both Ku protein paralogs in the Discoba group further supports the hypothesis that the gene duplication occurred early in eukaryotic evolution. It is plausible that archaea, which may have acted as intermediaries for Ku transfer, subsequently lost the Ku protein. Nonetheless, the extensive horizontal transfer of Ku among prokaryotes and its relatively higher prevalence in bacteria complicates our understanding of how Ku protein was inherited by early-branching eukaryotes.

Mutation of the LRG1 Rho-GAP gene is responsible for the hyper branching C-variant phenotype in the quorn mycoprotein fungus Fusarium venenatum A3/5

BACKGROUND

Quorn mycoprotein, a protein-rich meat alternative, is produced through large-scale fermentation of the fungus Fusarium venenatum. However, a major challenge during F. venenatum fermentation is the consistent appearance of mutants called colonial variants (C-variants). These C-variants have a highly branched morphology, which ultimately lead to a less desirable final product and early termination of the fermentation process. This study aimed to identify the genetic mutations responsible for C-variant morphology.

RESULTS

We first isolated both C-variant and wild-type strains from commercial fermentation samples and characterised radial growth rates on solid media. Whole genome sequencing facilitated the identification of mutations in a gene called jg4843 in 11 out of 12 C-variant isolates, which were not observed in the wild-type isolates. The jg4843 gene was identified as the ortholog of LRG1, a Rho-GTPase activating protein that regulates the Rho1 signalling pathway affecting fungal growth. Notably, the mutations in jg4843 were primarily located in the RhoGAP domain responsible for LRG1 activity. To confirm the role of these mutations, we used CRISPR/Cas9-mediated homology-directed recombination to introduce the C-variant mutations into the wild-type isolate, which successfully recapitulated the characteristic C-variant morphology.

CONCLUSIONS

This study identified mutations in the LRG1 ortholog jg4843 as the genetic cause of C-variant morphology in commercial fermentation F. venenatum isolates. Understanding this genetic basis paves the way for developing strategies to prevent C-variants arising, potentially leading to more efficient and sustainable production of Quorn mycoprotein.

Virus targeting as a dominant driver of interfacial evolution in the structurally resolved human-virus protein-protein interaction network

Regions on a host protein that interact with virus proteins (exogenous interfaces) frequently overlap with those that interact with other host proteins (endogenous interfaces), resulting in competition between hosts and viruses for these shared interfaces (mimic-targeted interfaces). Yet, the evolutionary consequences of this competitive relationship on the host are not well understood. Here, we integrate experimentally determined structures and homology-based templates of protein complexes with protein-protein interaction networks to construct a high-resolution human-virus structural interaction network. We perform site-specific evolutionary rate analyses on this structural interaction network and find that exogenous-specific interfaces evolve faster than endogenous-specific interfaces. Mimic-targeted interfaces evolve as fast as exogenous-specific interfaces, despite being targeted by both human and virus proteins. Our findings suggest that virus targeting plays a dominant role in host interfacial evolution within the context of domain-domain interactions and that mimic-targeted interfaces on human proteins are the key battleground for a mammalian-specific host-virus evolutionary arms race.

A chromosome-level genome assembly of eriophyoid mite Setoptus koraiensis

Eriophyoidea represents a highly diverse superfamily of herbivorous mites in the Acariformes, including over 5,000 named species that are distributed worldwide. However, the lack of chromosome-level genome prevents our understanding of the evolution in this group. Here, we report the first chromosome-level genome assembly of Setoptus koraiensis using Illumina, PacBio, and Hi-C sequencing technologies. The assembled genome has a size of 47 Mb with an N50 of 24.53 Mb, anchored into two chromosomes. The chromosome-level genome assembly had a BUSCO completeness of 89%. We identified 5,954 protein-coding genes, with 4,770 genes that could be functionally annotated. This genome provides resources to further understand the genetic and evolution of eriophyoid mites.

Characterization of HOL-30: a novel tandem-repeat galectin from the marine sponge Halichondria okadai

We here report the novel primary structure of a new member in the galectin family, the β-galactoside-binding lectin HOL-30, from the marine sponge Halichondria okadai, whose full-length sequence was determined thanks to the combination between Edman degradation and transcriptome analysis. The HOL-30 polypeptide is a tandem-repeat dimeric galectin, consisting of 281 amino acids, which includes two carbohydrate recognition domains (CRDs). Unlike most other galectins described in Porifera, HOL-30 did not have a signal peptide sequence for secretion. In solution, HOL-30 exhibited a molecular weight of 60 kDa, indicating a dimeric organization consisting of two 30 kDa tandem-repeat subunits stabilized by non-covalent interactions. Although the two CRDs had a similar predicted 3D structure, they displayed low pairwise sequence identity, approximately 20 %. HOL-30 exhibited glycan-binding affinities for type-1 (Galβ1-3GlcNAc) and type-2 (Galβ1-4GlcNAc) LacNAc. Furthermore, it also recognized blood type B-oligosaccharides on type-1 and type-2 LacNAc (Galα1-3Gal[Fucα1-2]β1-3/4GlcNAc), and blood type H-oligosaccharide on type-3 (Gal[Fucα1-2]β1-3GalNAcα). The glycan-binding properties of HOL-30 were compared with those of the hRTL galectin, previously identified in Chondrilla australiensis, consisting of tetrameric 15 kDa prototype subunits. The two sponge galectins displayed similar, but not identical, carbohydrate-binding properties, as evidenced by the fact that despite effectively binding to vertebrate cultured cells, HOL-30 had minimal impact on cell growth. Antiserum analysis revealed a mosaic distribution of HOL-30 in the parenchymal cells of sponge tissues within dense cell clusters surrounding the spicules.

The sodium leak channel NALCN is regulated by neuronal SNARE complex proteins

NALCN (sodium leak channel, nonselective) is vital for regulating electrical activity in neurons and other excitable cells, and mutations in the channel or its auxiliary proteins lead to severe neurodevelopmental disorders. Here, we show that the neuronal SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) complex proteins syntaxin and SNAP25 (synaptosome-associated protein 25), which enable synaptic transmission in the nervous system, inhibit the activity of the NALCN channel complex in both heterologous systems and primary neurons. The existence of this interaction suggests that the neurotransmitter release machinery can regulate electrical signaling directly and therefore modulate the threshold for its own activity. We further find that reduction of NALCN currents is sufficient to promote cell survival in syntaxin-depleted cells. This suggests that disinhibited NALCN may cause the puzzling phenomenon of rapid neuronal cell death in the absence of syntaxin. This interaction could offer opportunities for future drug development against genetic diseases linked to both NALCN- and SNARE protein-containing complexes.

Single specimen genome assembly of Culicoides stellifer shows evidence of a non-retroviral endogenous viral element

BACKGROUND

Advancing our knowledge of vector species genomes is a key step in our battle against the spread of diseases. Biting midges of the genus Culicoides are vectors of arboviruses that significantly affect livestock worldwide. Culicoides stellifer is a suspected vector with a wide range distribution in North America, for which cryptic diversity has been described.

RESULTS

With just one specimen of C. stellifer, we assembled and annotated the nuclear and mitochondrial genome using the ultra-low input DNA PacBio protocol. The genome assembly is 119 Mb in length with a contig N50 value of 479.3 kb, contains 11% repeat sequences and 18,895 annotated protein-coding genes. To further elucidate the role of this species as a vector, we provide genomic evidence of a non-retroviral endogenous viral element integrated into the genome that corresponds to rhabdovirus nucleocapsid proteins, the same family as the vesicular stomatitis virus.

CONCLUSIONS

This genomic information will pave the way for future investigations into this species's putative vector role. We also demonstrate the practicability of completing genomic studies in small dipterans using single specimens preserved in ethanol as well as introduce a workflow for data analysis that considers the challenges of insect genome assembly.

Chromosome-level genome of the brown lacewing Micromus angulatus (Stephens, 1836) (Neuroptera: Hemerobiidae)

The brown lacewing Micromus angulatus (Stephens), a member within the subfamily Microminae of the family Hemerobiidae, is a globally distributed species and an important predatory natural enemy of various agricultural and forestry crop pests. Despite its global distribution and agricultural significance, genomic resources for the development of novel pest management strategies for M. angulatus and the family Hemerobiidae remain limited. Here, we present the first high-quality chromosome-level reference genome for M. angulatus, achieved through PacBio HiFi and Hi-C technologies. The assembled genome spans 1.29 Gb, with scaffold N50 of 63.78 Mb, and complete BUSCO score of 97.00%, containing eight pseudochromosomes, demonstrates a high degree of continuity. Functional annotation identified 13,250 protein-coding genes, and repetitive sequences, which account for 83.65% of the genome, were also characterized. This comprehensive assembly offers a robust reference for in-depth research on the genetic basis of predation and adaptability in M. angulatus, supports further research into the understanding the genetic diversity within Hemerobiidae, contributing to the broader knowledge within the group and its related species.

Revisiting the Plasmodium falciparum druggable genome using predicted structures and data mining

Identification of novel drug targets is a key component of modern drug discovery. While antimalarial targets are often identified through the mechanism of action studies on phenotypically derived inhibitors, this method tends to be time- and resource-consuming. The discoverable target space is also constrained by existing compound libraries and phenotypic assay conditions. Leveraging recent advances in protein structure prediction, we systematically assessed the Plasmodium falciparum genome and identified 867 candidate protein targets with evidence of small-molecule binding and blood-stage essentiality. Of these, 540 proteins showed strong essentiality evidence and lack inhibitors that have progressed to clinical trials. Expert review and rubric-based scoring of this subset based on additional criteria such as selectivity, structural information, and assay developability yielded 27 high-priority antimalarial target candidates. This study also provides a genome-wide data resource for P. falciparum and implements a generalizable framework for systematically evaluating and prioritizing novel pathogenic disease targets.

Evolution of Transcription Factor-containing Superfamilies in Eukaryotes

Regulation of gene expression helps determine various phenotypes in most cellular life forms. It is orchestrated at different levels and at the point of transcription initiation by transcription factors (TFs). TFs bind to DNA through domains that are evolutionarily related, by shared membership of the same superfamilies (TF-SFs), to those found in other nucleic acid binding and protein-binding functions (nTFs for non-TFs). Here we ask how TF DNA binding sequence families in eukaryotes have evolved in relation to their nTF relatives. TF numbers scale by power law with the total number of protein-coding genes differently in different clades, with fungi usually showing sub-linear powers whereas chordates show super-linear scaling. The LECA probably encoded a complex regulatory machinery with both TFs and nTFs, but with an excess of nTFs when compared to the relative distribution of TFs and nTFs in extant organisms. Losses drive the evolution of TFs and nTFs, with the possible exception of TFs in animals for some tree topologies. TFs are highly dynamic in evolution, showing higher gain and loss rates than nTFs in some TF-SFs though both are conserved to similar extents. Gains of TFs and nTFs are driven by the appearance of a large number of new sequence clusters in a small number of nodes, which determine the presence of as many as a third of extant TFs and nTFs as well as the relative presence of TFs and nTFs. Whereas nodes showing explosion of TF numbers belong to multicellular clades, those for nTFs lie among the fungi and the protists.

Damsels in Disguise: Development of Ultraviolet Sensitivity and Colour Patterns in Damselfishes (Pomacentridae)

Damselfishes (Pomacentridae) are widespread and highly abundant on tropical coral reefs. They exhibit diverse body colouration within and between the ~250 species and across ontogenetic stages. In addition to human-visible colours (i.e., 400-700 nm), most adult damselfishes reflect ultraviolet (UV, 300-400 nm) colour patches. UV sensitivity and UV colour signals are essential for feeding and form the basis for a secret communication channel invisible to the many UV-blind predatory fish on the reef; however, how these traits develop across ontogenetic stages and their distribution across the damselfish family is poorly characterised. Here, we used UV photography, phylogenetic reconstructions of opsin genes, and differential gene expression analysis (DGE) of retinal samples to investigate the development of UV vision and colour patterns in three ontogenetic stages (pre-settlement larval, juvenile, and adult) of 11 damselfish species. Using DGE, we found similar gene expression between juveniles and adults, which strongly differed from larvae. All species and all stages expressed at least one UV-sensitive sws1 opsin gene. However, UV body colour patterns only started to appear at the juvenile stage. Moreover, Pomacentrus species displayed highly complex UV body patterns that were correlated with the expression of two sws1 copies. This could mean that some damselfishes can discriminate colours that change only in their UV component. We demonstrate dramatic shifts in both UV sensitivity and UV colouration across the development stages of damselfish while highlighting the importance of considering ontogeny when studying the coevolution of visual systems and colour signals.

A Conserved Somatic Sex Determination Cascade Instructs Trait-Specific Sexual Dimorphism in Horned Dung Beetles

Sex-specific trait expression represents a striking dimension of morphological variation within and across species. The mechanisms instructing sex-specific organ development have been well studied in a small number of insect model systems, suggesting striking conservation in some parts of the somatic sex determination pathway while hinting at possible evolutionary lability in others. However, further resolution of this phenomenon necessitates additional taxon sampling, particularly in groups in which sexual dimorphisms have undergone significant elaboration and diversification. Here, we functionally investigate the somatic sex determination pathway in the gazelle dung beetle Digitonthophagus gazella, an emerging model system in the study of the development and evolution of sexual dimorphisms. We find that RNA interference (RNAi) targeting transformer (tra) caused chromosomal females to develop morphological traits largely indistinguishable from those normally only observed in males, and that traRNAi is sufficient to induce splicing of the normally male-specific isoform of doublesex in chromosomal females, while leaving males unaffected. Further, intersexRNAi was found to phenocopy previously described RNAi phenotypes of doublesex in female but not male beetles. These findings match predictions derived from models of the sex determination cascade as developed largely through studies in Drosophila melanogaster. In contrast, efforts to target transformer2 via RNAi resulted in high juvenile mortality but did not appear to affect doublesex splicing, whereas RNAi targeting Sex-lethal and two putative orthologs of hermaphrodite yielded no obvious phenotypic modifications in either males or females, raising the possibility that the function of a subset of sex determination genes may be derived in select Diptera and thus nonrepresentative of their roles in other holometabolous orders. Our results help illuminate how the differential evolutionary lability of the somatic sex determination pathway has contributed to the extraordinary morphological diversification of sex-specific trait expression found in nature.

A pangenome reveals LTR repeat dynamics as a major driver of genome evolution in Chenopodium

The genus Chenopodium L. is characterized by its wide geographic distribution and ecological adaptability. Species such as quinoa (Chenopodium quinoa Willd.) have served as domesticated staple crops for centuries. Wild Chenopodium species exhibit diverse niche adaptations and are important genetic reservoirs for beneficial agronomic traits, including disease resistance and climate hardiness. To harness the potential of the wild taxa for crop improvement, we developed a Chenopodium pangenome through the assembly and comparative analyses of 12 Chenopodium species that encompass the eight known genome types (A-H). Six of the species are new chromosome-scale assemblies, and many are polyploids; thus, a total of 20 genomes were included in the pangenome analyses. We show that the genomes vary dramatically in size with the D genome being the smallest (∼370 Mb) and the B genome being the largest (∼700 Mb) and that genome size was correlated with independent expansions of the Copia and Gypsy LTR retrotransposon families, suggesting that transposable elements have played a critical role in the evolution of the Chenopodium genomes. We annotated a total of 33,457 pan-Chenopodium gene families, of which ∼65% were classified as shell (2% private). Phylogenetic analysis clarified the evolutionary relationships among the genome lineages, notably resolving the taxonomic placement of the F genome while highlighting the uniqueness of the A genome in the Western Hemisphere. These genomic resources are particularly important for understanding the secondary and tertiary gene pools available for the improvement of the domesticated chenopods while furthering our understanding of the evolution and complexity within the genus.

A conserved somatic sex determination cascade instructs trait-specific sexual dimorphism in horned dung beetles

Sex-specific trait expression represents a striking dimension of morphological variation within and across species. The mechanisms instructing sex-specific organ development have been well studied in a small number of insect model systems, suggesting striking conservation in some parts of the somatic sex determination pathway while hinting at possible evolutionary lability in others. However, further resolution of this phenomenon necessitates additional taxon sampling, particularly in groups in which sexual dimorphisms have undergone significant elaboration and diversification. Here, we functionally investigate the somatic sex determination pathway in the gazelle dung beetle Digitonthophagus gazella, an emerging model system in the study of the development and evolution of sexual dimorphisms. We find that RNA interference (RNAi) targeting transformer (tra) caused chromosomal females to develop morphological traits largely indistinguishable from those normally only observed in males, and that tra RNAi is sufficient to induce splicing of the normally male-specific isoform of doublesex in chromosomal females, while leaving males unaffected. Further, intersex RNAi was found to phenocopy previously described RNAi phenotypes of doublesex in female but not male beetles. These findings match predictions derived from models of the sex determination cascade as developed largely through studies in Drosophila melanogaster. In contrast, transformer2 NAi resulted in larval mortality and was not sufficient to affect doublesex splicing, whereas RNAi targeting Sex-lethal and two putative orthologs of hermaphrodite yielded no obvious phenotypic modifications in either males or females, raising the possibility that the function of a subset of sex determination genes may be derived in select Diptera and thus non-representative of their roles in other holometabolous orders. Our results help illuminate how the differential evolutionary lability of the somatic sex determination pathway has contributed to the extraordinary morphological diversification of sex-specific trait expression found in nature.

Molecular underpinnings of plasticity and supergene-mediated polymorphism in fire ant queens

Characterizing molecular underpinnings of plastic traits and balanced polymorphisms represent 2 important goals of evolutionary biology. Fire ant gynes (pre-reproductive queens) provide an ideal system to study potential links between these phenomena because they exhibit both supergene-mediated polymorphism and nutritional plasticity in weight and colony-founding behaviour. Gynes with the inversion supergene haplotype are lightweight and depend on existing workers to initiate reproduction. Gynes with only the ancestral, non-inverted gene arrangement accumulate more nutrient reserves as adults and, in a distinct colony-founding behaviour, initiate reproduction without help from workers. However, when such gynes overwinter in the natal nest they develop an environmentally induced lightweight phenotype and colony-founding behaviour, similar to gynes with the inversion haplotype that have not overwintered. To evaluate the extent of shared mechanisms between plasticity and balanced polymorphism in fire ant gyne traits, we assessed whether genes with expression variation linked to overwintering plasticity may be affected by the evolutionary divergence between supergene haplotypes. To do so, we first compared transcriptional profiles of brains and ovaries from overwintered and non-overwintered gynes to identify plasticity-associated genes. These genes were enriched for metabolic and behavioural functions. Next, we compared plasticity-associated genes to those differentially expressed by supergene genotype, revealing a significant overlap of the 2 sets in ovarian tissues. We also identified sequence substitutions between supergene variants of multiple plasticity-associated genes, consistent with a scenario in which an ancestrally plastic phenotype responsive to an environmental condition became increasingly genetically regulated.

Draft genome of the endangered visayan spotted deer (Rusa alfredi), a Philippine endemic species

The Visayan Spotted Deer (VSD), or Rusa alfredi, is an endangered and endemic species in the Philippines. Despite its status, genomic information on R. alfredi, and the genus Rusa in general, is missing. This study presents the first draft genome assembly of the VSD using the Illumina short-read sequencing technology. The resulting RusAlf_1.1 assembly has a 2.52 Gb total length, with a contig N50 of 46 Kb and scaffold N50 size of 75 Mb. The assembly has a BUSCO complete score of 95.5%, demonstrating the genome's completeness, and includes the annotation of 24,531 genes. Our phylogenetic analysis based on single-copy orthologs revealed a close evolutionary relationship between R. alfredi and the genus Cervus. RusAlf_1.1 represents a significant advancement in our understanding of the VSD. It opens opportunities for further research in population genetics and evolutionary biology, potentially contributing to more effective conservation and management strategies for this endangered species.

Chromosome-level genome assembly for the ecologically and economically important alga Saccharina japonica

Saccharina japonica is a major kelp species of brown algae with the highest production among aquaculture seaweeds and holds important ecological and economic value. Despite advancements in domestication, a high-quality chromosome-level genome assembly is needed to assist its genetic improvement. Previous genome assemblies of S. japonica were either on a draft-level or highly fragmented. Here, we generated a high-quality chromosome-level genome for the female sporophyte using PacBio sequencing and Hi-C. The genome is 516.11 Mb, with contig N50 length of 491.30 Kb and scaffold N50 length of 16.24 Mb, anchored into 32 pseudo-chromosomes. Repetitive elements constituted 45.07% of the genome, and 17,739 protein-coding genes were predicted, of which 82% were functionally annotated. This genome provides a crucial resource for biotechnological advances in S. japonica breeding and offers insights into the ecology and evolution of brown algae.

Exploring the genetics of social behaviour in C. calcarata

Studies investigating social evolution often focus on species that are obligately eusocial, where presumably all of the adaptive genetic changes associated with sociality have already been completed. To fully understand eusociality, we must study species with facultative social behaviour. The small carpenter bee Ceratina calcarata is an ideal model for studying the genetics and molecular biology of eusocial evolution as it can exhibit both subsocial behaviour with parental care and social behaviour facilitated by the altruistic dwarf eldest daughter. Here, we sequenced the genomes of subsocial and social C. calcarata to identify mutations and genes associated with social behaviour and used these data to test several hypotheses related to the evolution of eusociality. Many single nucleotide polymorphisms that had high levels of genetic differentiation (Fst) between social and subsocial C. calcarata were in or near genes or regions important for regulating gene expression. These results are consistent with the Genetic Toolkit Hypothesis of eusocial evolution. Our findings suggest that the low behavioural complexity observed in C. calcarata may involve modulation of existing regulatory genes and gene networks to generate phenotypes associated with social behaviour.

Antagonistic kinesin-14s within a single chromosomal drive haplotype

In maize, there are two meiotic drive systems that operate on large tandem repeat arrays called knobs that are found on chromosome arms. One meiotic drive haplotype, Abnormal chromosome 10 (Ab10), encodes two kinesin proteins that interact with two distinct tandem repeat arrays in a sequence-specific manner to confer meiotic drive. The kinesin KINDR associates with knob180 repeats while the kinesin TRKIN associates with TR-1 repeats. Prior data show that meiotic drive is conferred primarily by the KINDR/knob180 system, with the TRKIN/TR-1 system having little or no role. The second meiotic drive haplotype, K10L2, shows low levels of meiotic drive and only encodes the TRKIN/TR-1 system. Here we used long-read sequencing to assemble the K10L2 haplotype and showed that it has strong homology to an internal portion of the Ab10 haplotype. We also carried out CRISPR mutagenesis of Trkin to test the role of Trkin on Ab10 and K10L2. The data indicate that the Trkin gene on Ab10 does not improve drive or fitness but instead has a weak deleterious effect when paired with a normal chromosome 10. The deleterious effect is more severe when Ab10 is paired with K10L2: in this context functional Trkin on either chromosome nearly abolishes Ab10 drive. We modeled the effect of Trkin on Ab10 and found it should not persist in the population. We conclude that Trkin either confers an advantage to Ab10 in untested circumstances or that it is in the process of being purged from the Ab10 population.

Genome assembly, gene content, and plastic gene expression responses to salinity changes in the Brackishwater Clam (Corbicula japonica) from a dynamic estuarine environment

Estuaries are dynamic transition zones between marine and freshwater environments, where salinity varies greatly on spatial and temporal scales. The temporal salinity fluctuations of these habitats require organisms to rapidly regulate ionic concentrations and osmotic pressure to survive in these dynamic conditions. Understanding the extent of plasticity of euryhaline animals is vital for predicting their responses and resilience to salinity change. We generated the first high-resolution genome and transcriptome sequences of C. japonica. In comparison with 11 other molluscan genomes, the C. japonica genome displayed striking expansions of putative neuron-related genes and gene families. The involvement of these genes in the glutamate/GABA-glutamine and glycine cycle suggests a possible contribution to the excitation of neuronal networks, particularly under high salinity conditions. This study contributes to our understanding of mechanisms underlying the rapid responses of estuarine species to changing conditions and raises many intriguing hypotheses and questions for future investigation.

Tyrosine transfer RNA levels and modifications during blood-feeding and vitellogenesis in the mosquito, Aedes aegypti

Mosquitoes such as Aedes aegypti must consume a blood meal for the nutrients necessary for egg production. Several transcriptome and proteome changes occur post-blood meal that likely corresponds with codon usage alterations. Transfer RNA (tRNA) is the adapter molecule that reads messenger RNA codons to add the appropriate amino acid during protein synthesis. Chemical modifications to tRNA enhance codon decoding, improving the accuracy and efficiency of protein synthesis. Here, we examined tRNA modifications and transcripts associated with the blood meal and subsequent periods of vitellogenesis in A. aegypti. More specifically, we assessed tRNA transcript abundance and modification levels in the fat body at critical times post blood-feeding. Based on a combination of alternative codon usage and identification of particular modifications, we discovered that increased transcription of tyrosine tRNAs is likely critical during the synthesis of egg yolk proteins in the fat body following a blood meal. Altogether, changes in both the abundance and modification of tRNA are essential factors in the process of vitellogenin production after blood-feeding in mosquitoes.

Convergent Gene Duplication in Arctic and Antarctic Teleost Fishes

Teleost fishes have independently colonized polar regions multiple times, facing many physiological and biochemical challenges due to frigid temperatures. Although increased gene copy numbers can contribute to adaptive evolution in extreme environments, it remains unclear which categories of genes exhibit increased copy numbers associated with polar colonization. Using 104 species of ray-finned fishes, we systematically identified genes with a significant correlation between copy number and polar colonization after phylogenetic correction. Several genes encoding extracellular glycoproteins, including zona pellucida (ZP) proteins, which increase their copy number in Antarctic notothenioid fishes, exhibited elevated copy numbers across multiple polar fish lineages. Additionally, some genes reported to be highly expressed under cold stress, such as cold-inducible RNA-binding protein (CIRBP), had significantly increased copy numbers in polar fishes. Further analysis will provide a fundamental basis for understanding the role of gene duplication in polar adaptations.

A Mouse Model of Sporadic Alzheimer's Disease with Elements of Major Depression

After olfactory bulbectomy, animals are often used as a model of major depression or sporadic Alzheimer's disease and, hence, the status of this model is still disputable. To elucidate the nature of alterations in the expression of the genome after the operation, we analyzed transcriptomes of the cortex, hippocampus, and cerebellum of the olfactory bulbectomized (OBX) mice. Analysis of the functional significance of genes in the brain of OBX mice indicates that the balance of the GABA/glutamatergic systems is disturbed with hyperactivation of the latter in the hippocampus, leading to the development of excitotoxicity and induction of apoptosis in the background of severe mitochondrial dysfunction and astrogliosis. On top of this, the synthesis of neurotrophic factors decreases leading to the disruption of the cytoskeleton of neurons, an increase in the level of intracellular calcium, and the activation of tau protein hyperphosphorylation. Moreover, the acetylcholinergic system is deficient in the background of the hyperactivation of acetylcholinesterase. Importantly, the activity of the dopaminergic, endorphin, and opiate systems in OBX mice decreases, leading to hormonal dysfunction. On the other hand, genes responsible for the regulation of circadian rhythms, cell migration, and innate immunity are activated in OBX animals. All this takes place in the background of a drastic downregulation of ribosomal protein genes in the brain. The obtained results indicate that OBX mice represent a model of Alzheimer's disease with elements of major depression.