Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

Effect of oral fluid in urine samples on the analysis of selected erythropoietin receptor agonists and detection of saliva-specific markers for doping control purposes

Due to their performance-enhancing effect, erythropoiesis-stimulating agents (ESAs) are banned at all times by the World Anti-Doping Agency (WADA) in competitive sports. Doping control analyses for such compounds are routinely performed using gel electrophoretic and immunoblotting techniques, and degradation of the analytes can severely impair detection, results evaluation and interpretation. As oral fluid (OF) contains significant amounts of proteases, the question of whether its addition to a doping control urine sample may impede anti-doping analysis needs to be addressed. Intentional tampering attempts are likewise prohibited by WADA and require a detection strategy. It was observed that the addition of OF can indeed lead to impairments of ESA analyses, though the fraction of unidentifiable ESA signals varies depending on several factors, such as the individual composition of the OF, the sex of the OF donor, the time of sampling, the OF volume and the incubation conditions. Overall, 20 % of all generally valid analyses were classified as unidentifiable, 12 % as impaired, and 69 % as identifiable, highlighting the relevance for strategies that allow for the identification of OF in urine. While human salivary α-amylase was found insufficiently reliable as a marker, peptides of salivary proline rich proteins (saPRP) were shown to be both specific for OF and traceable with adequate sensitivity using a newly developed LC-HRMS/MS method. The approach was comprehensively characterized shown to be fit-for-purpose for routine doping controls where tampering attempts with OF are suspected.

Differential Expression of Hormone-Related Genes in the Heads of Adult and Nymphal Woodroaches (Cryptocercus)

Termites are eusocial cockroaches, but the crucial distinctions in gene expression during the evolution of eusociality remain unclear. One reason for the lack of this information is that comparative transcriptome analysis of termites with their sister group, the cockroach genus Cryptocercus, has not been conducted. We identified genes associated with three vital hormones (juvenile hormone [JH], 20-hydoroxyecdysone [20E], and insulin) from the genome sequence of Cryptocercus punctulatus and conducted RNA-seq analysis using the heads of female/male adults and nymphs to elucidate their expression levels. The comprehensive gene expression analysis revealed a multitude of genes exhibiting differences in expression between developmental stages rather than between sexes. Subsequently, we compared the differences in expression patterns of each hormone-related gene by combining the results of a previous RNA-seq study conducted on the heads of castes (reproductives, workers, and soldiers) in the termite Reticulitermes speratus. The results indicated that genes with expression differences among castes in R. speratus, particularly those related to JH and 20E, were significantly more abundant compared to genes with expression differences between adults and nymphs in C. punctulatus. While no significant difference was observed in the number of genes within the insulin signaling pathway, a trend of homologs highly expressed in adult woodroaches but not in adult termites was observed, and the expression patterns of positive and negative regulators in the pathway differed significantly between adults and nymphs. The differences in the expression patterns between Cryptocercus and termites are believed to reflect variations in hormone levels and signaling activities between adults and juveniles, the latter encompassing workers and soldiers in the case of termites.

Structure of a photosystem I supercomplex from Galdieria sulphuraria close to an ancestral red alga

Red algae exhibit unique photosynthetic adaptations, characterized by photosystem I (PSI) supercomplexes containing light-harvesting complexes (LHCs), forming PSI-LHCI supercomplexes. In this study, we solved the PSI-LHCI structure of Galdieria sulphuraria NIES-3638 at 2.19-angstrom resolution using cryo-electron microscopy, revealing a PSI monomer core associated with seven LHCI subunits. Structural analysis uncovered the absence of phylloquinones, the common secondary electron acceptor in PSI of photosynthetic organisms, suggesting adaptation to a benzoquinone-like molecule. Phylogenetic analysis suggests that G. sulphuraria retains traits characteristic of an ancestral red alga, including distinctive LHCI binding and interaction patterns. Variations in LHCI composition and interactions across red algae, particularly in red-lineage chlorophyll a/b-binding-like protein and red algal LHCs, highlight evolutionary divergence and specialization. These findings not only deepen our understanding of red algal PSI-LHCI diversification but also enable us to predict features of an ancestral red algal PSI-LHCI supercomplex, providing a framework to explore evolutionary adaptations from an ancestral red alga.

PIGOME: An Integrated and Comprehensive Multi-omics Database for Pig Functional Genomics Studies

In addition to being a major source of animal protein, pigs are an important model for studying development and diseases in humans. Over the past two decades, thousands of high-throughput sequencing studies in pigs have been performed using a variety of tissues from different breeds and developmental stages. However, multi-omics databases specifically designed for pig functional genomics research are still limited. Here, we present PIGOME, a user-friendly database of pig multi-omes. PIGOME currently contains seven types of pig omics datasets, including whole-genome sequencing (WGS), RNA sequencing (RNA-seq), microRNA sequencing (miRNA-seq), chromatin immunoprecipitation sequencing (ChIP-seq), assay for transposase-accessible chromatin sequencing (ATAC-seq), bisulfite sequencing (BS-seq), and methylated RNA immunoprecipitation sequencing (MeRIP-seq), from 6901 samples and 392 projects with manually curated metadata, integrated gene annotation, and quantitative trait locus information. Furthermore, various "Explore" and "Browse" functions have been established to provide user-friendly access to omics information. PIGOME implements several tools to visualize genomic variants, gene expression, and epigenetic signals of a given gene in the pig genome, enabling efficient exploration of spatiotemporal gene expression/epigenetic patterns, functions, regulatory mechanisms, and associated economic traits. Collectively, PIGOME provides valuable resources for pig breeding and is helpful for human biomedical research. PIGOME is available at https://pigome.com.

A rare stop-gain SNP mutation in BrGL2 causes aborted trichome development in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

KEY MESSAGE

A rare stop-gain SNP mutation in BrGL2 confers short hair phenotype of Chinese cabbage via bulked-segregant analysis sequencing, fine-mapping and gene silencing analysis. Trichomes negatively affect the quality of Chinese cabbage, a leafy vegetable crop in the family Brassicaceae. The short hair trait is caused by abnormal trichome development. In this study, the BraA07g025490.3C gene was identified as a candidate gene for the short hair trait in Chinese cabbage by BSA-seq and fine-mapping analyses. It was subsequently named BrGL2 because of its strong homology to AtGL2 (At1g79840). Sequence analysis indicated that a C to G single-nucleotide polymorphism (SNP) mutation in the sixth exon of BrGL2 produced a premature stop codon in the HCW (short hair) line, resulting in a loss-of-function mutation of BrGL2. This stop-gain SNP mutation was found exclusively in the HCW line, and not in 524 diverse B. rapa accessions. Further analysis by virus-induced gene silencing showed that the knock-down of BrGL2 in HN19-G lines (wild-type hair) reduced the size of leaf trichomes. BrGL2 affected trichome development probably by impacting the expression of downstream transcription factor genes and cell wall-related genes, as determined by comparative transcriptome analyses of wild type and short hair lines. On the basis of the identification and verification of the key stop-gain SNP mutation in BrGL2 resulting in aborted trichome development in Chinese cabbage, we propose a model for trichome development.

Large-scale transcriptome mining enables macrocyclic diversification and improved bioactivity of the stephanotic acid scaffold

Nearly 10,000 plant species are represented by RNA-seq datasets in the NCBI sequence read archive, which are difficult to search in unassembled format due to database size. Here, we optimize RNA-seq assembly to transform most of this public RNA-seq data to a searchable database for biosynthetic gene discovery. We test our transcriptome mining pipeline towards the diversification of moroidins, which are plant ribosomally-synthesized and posttranslationally-modified peptides (RiPPs) biosynthesized from copper-dependent peptide cyclases. Moroidins are bicyclic compounds with a conserved stephanotic acid scaffold, which becomes cytotoxic to non-small cell lung adenocarcinoma cells with an additional C-terminal macrocycle. We discover moroidin analogs with second ring structures diversified at the crosslink and the non-crosslinked residues including a moroidin analog from water chickweed, which exhibits higher cytotoxicity against lung adenocarcinoma cells than moroidin. Our study expands stephanotic acid-type peptides to grasses, Lowiaceae, mints, pinks, and spurges while demonstrating that large-scale transcriptome mining can broaden the medicinal chemistry toolbox for chemical and biological exploration of eukaryotic RiPP lead structures.

Genomic analysis of Proteus mirabilis: Unraveling global epidemiology and antimicrobial resistance dissemination - emerging challenges for public health and biosecurity

Given the escalating public health threat posed byProteus mirabilis(P. mirabilis) and its rapidly evolving drug resistance, it is imperative to elucidate its global epidemiology and resistance mechanisms through a comprehensive genomic lens. As of August 2024, 3,403 high-qualityP. mirabilisgenomes were retrieved from public databases (total 3,752), spanning 58 countries/regions, with the United States showing the highest report rate (52.51 %). Human-derived isolates, particularly from urine (34.47 %), were the primary source. A total of 239 antibiotic resistance genes (ARGs) were identified in P. mirabilis, with β-lactamase and carbapenemase genes being particularly widespread and isolates from China harboring the highest ARG counts. Globally,P. mirabilisisolates were categorized into 17 distinct clusters, with U.S. isolates showing the widest phylogenetic spread. Minimal SNP variations among isolates from different countries and hosts suggest transnational and cross-host clonal propagation. Frequent clonal transmission was also observed among diverse hosts and clinical sources.P. mirabiliscarries numerous integrative and conjugative elements (ICEs), some facilitating ARG dissemination (n = 215). Prophages, though ubiquitous, contributed minimally to ARG spread. Spearman's analysis revealed significant correlations between ARGs and insertion sequences (ISs), replicons, and ICEs. Ancestral state analysis indicated prophages were mainly acquired horizontally, while other mobile genetic elements (MGEs) were largely clonally transmitted. This study provides the first comprehensive genomic analysis ofP. mirabilis's global resistance landscape, highlighting the need to designate it as a novel antimicrobial resistance indicator and implement long-term surveillance.

Solanaceae pan-genomes reveal extensive fractionation and functional innovation of duplicated genes

The Solanaceae family contains many agriculturally important crops, including tomato, potato, pepper, and tobacco, as well as other species with potential for agricultural development, such as the orphan crops groundcherry, wolfberry, and pepino. Research progress varies greatly among these species, with model crops like tomato being far ahead. This disparity limits the broader agricultural application of other Solanaceae species. In this study, we constructed an interspecies pan-genome for the Solanaceae family and identified various gene retention patterns. Our findings reveal that the activity of specific transposable elements is closely associated with gene fractionation and transposition. The pan-genome was further resolved at the level of T subgenomes, which were generated by Solanaceae-specific paleo-hexaploidization (T event). We demonstrate substantial gene fractionation (loss) and divergence events following ancient duplications. For example, all class A and E flower model genes in Solanaceae originated from two tandemly duplicated genes, which expanded through the γ and T events before fractionating into 10 genes in tomato, each acquiring distinct functions critical for fruit development. Based on these results, we developed the Solanaceae Pan-Genome Database (SolPGD, http://www.bioinformaticslab.cn/SolPGD), which integrates datasets from both inter- and intra-species pan-genomes of Solanaceae. These findings and resources will facilitate future studies of solanaceous species, including orphan crops.

Phylogeny-metabolism dual-directed single-cell genomics for dissecting and mining ecosystem function by FISH-scRACS-seq

Microbiome-wide association studies (MWASs) have uncovered microbial markers linked to ecosystem traits, but the mechanisms underlying their functions can remain elusive. This is largely due to challenges in validating their in situ metabolic activities and tracing such activities to individual genomes. Here, we introduced a phylogeny-metabolism dual-directed single-cell genomics approach called fluorescence-in situ-hybridization-guided single-cell Raman-activated sorting and sequencing (FISH-scRACS-seq). It directly localizes individual cells from target taxon via an FISH probe for marker organism, profiles their in situ metabolic functions via single-cell Raman spectra, sorts cells of target taxonomy and target metabolism, and produces indexed, high-coverage, and precisely-one-cell genomes. From cyclohexane-contaminated seawater, cells representing the MWAS-derived marker taxon of γ-Proteobacteria and that are actively degrading cyclohexane in situ were directly identified via FISH and Raman, respectively, then sorted and sequenced for one-cell full genomes. In such a Pseudoalteromonas fuliginea cell, we discovered a three-component cytochrome P450 system that can convert cyclohexane to cyclohexanol in vitro, representing a previously unknown group of cyclohexane-degrading enzymes and organisms. Therefore, by unveiling enzymes, pathways, genomes, and their in situ cellular functions specifically for those organisms with ecological relevance at one-cell resolution, FISH-scRACS-seq is a rational and generally applicable approach to dissecting and mining microbiota functions.

Birmingham-group IncP-1α plasmids revisited: RP4, RP1 and RK2 are identical and their remnants can be detected in environmental isolates

RP4, RP1, RK2 and R68 were isolated from the multidrug-resistant bacterial wound isolates in 1969 in the Birmingham Accident Hospital, Birmingham, England, and collectively called Birmingham-group IncP-1α plasmids. These plasmids have been widely used as models to study different aspects of plasmid biology, develop genetic delivery systems and design plasmid vectors. Early studies showed that these plasmids conferred the same antibiotic resistance profile, had a similar size and were undistinguishable from each other using DNA heteroduplex electron microscopy and restriction endonuclease analyses. These observations have led to the widely held assumption that they are identical, although there has been no conclusive supporting evidence. In this work, we sequenced the plasmids RP1 and RP4 from our laboratory strain collection and compared these new sequences with the plasmids RP4 and RK2 assembled from a publicly available sequencing database, showing that the RP1, RP4 and RK2 plasmids are 60 095 bp in length and identical at the nucleotide resolution. Noteworthily, the plasmid sequence is highly conserved despite having been distributed to different labs over 50 years and propagated in different bacterial hosts, strengthening the previous observation that the bacterial host adapts to the RP4/RP1/RK2 plasmid rather than the opposite. In the updated RP4/RP1/RK2 sequence, we found a fusion gene, called pecM-orf2, that was formed putatively by a genetic deletion event. By searching for pecM-orf2 in the National Center for Biotechnology Information database, we detected remnants of the RP4/RP1/RK2 plasmid that carry features of laboratory-engineered vectors in bacterial environmental isolates, either in their chromosome or as a plasmid. This suggests a leak of these plasmids from the laboratory into the environment, which may subsequently impact bacterial evolution and raises concerns about the biocontainment of engineered plasmids when being handled in laboratory settings.

Blue-shifted ancyromonad channelrhodopsins for multiplex optogenetics

Light-gated ion channels from protists (channelrhodopsins or ChRs) are optogenetic tools widely used for controlling neurons and cardiomyocytes. Multiplex optogenetic applications require spectrally separated molecules that must be found in nature, as they are difficult to engineer without disrupting channel function. Scanning numerous sequence databases, we identified three robust naturally blue-shifted ChRs from ancyromonads. They form a separate branch on the phylogenetic tree and contain residue motifs characteristic of anion ChRs (ACRs). However, only two conduct chloride, whereas the close Nutomonas longa homolog (peak absorption at ~440 nm) generates inward cation currents in mammalian cells under physiological conditions, significantly exceeding those by previously known tools. Measurements of transient absorption changes and pH titration of purified Ancyromonas sigmoides ACR (AnsACR) combined with mutant analysis revealed the roles of the residues in the photoactive site. Both ancyromonad ACRs allowed optogenetic silencing of mouse cortical neurons in brain slices. AnsACR expression in the cholinergic neurons enabled photoinhibition of pharyngeal muscle contraction in live worms. AnsACR could be activated by near-infrared two-photon illumination, which is required to control specific neurons in thick tissue. Our results improved the mechanistic understanding of light-gated channel function and expanded the optogenetic toolkit.

Impact statement

Ancyromonad channelrhodopsins advance our understanding of ionic selectivity and wavelength regulation in light-gated ion channels and also expand the toolkit for all-optical electrophysiology.

Increasing rates of erm(B) and erm(N) in human Campylobacter coli and Campylobacter jejuni erythromycin-resistant isolates between 2018 and 2023 in France

Macrolides are the first-line compounds used for the treatment of campylobacteriosis. Macrolide resistance remains low in France, with mutations in 23S rDNA being the main associated resistance mechanism. However, two erythromycin methyltransferases have also been identified: erm(B), which is mainly described in animal reservoirs, and erm(N), which is strictly described in humans. In France, between 2018 and 2023, erythromycin-resistant Campylobacter species strains were systematically sequenced and analyzed via an in-house bioinformatics pipeline, leading to the identification of the resistomes, MLST and cgMLST, as well as the characterization of the source of contamination. In this study, the genomes of 280 erythromycin-resistant strains were sequenced over a 6-year period. The identification of erythromycin-associated resistance markers revealed a predominance of 23S rDNA mutations, in 90% of cases, but also erm-type methyltransferases in 10% of cases: 75% for erm(N) and 25% for erm(B). Over this period, an important increase in the rate of erm-positive isolates was observed: 2% in 2018 compared with 13% in 2023, with 10% for erm(N) and 3% for erm(B). erm(N) has been found exclusively within a CRISPR-Cas9 operon, whereas erm(B) has been found within diverse types of resistance genomic islands. Each erm(N)- or erm(B)-positive isolate had at least two other resistance markers (mostly ciprofloxacin, tetracycline, or ampicillin) and often carried aminoglycoside-associated resistance genes. The majority of the erm-positive isolates were obtained from chicken. The increasing rates of erm-positive and multiresistant isolates make the monitoring of erythromycin-resistant Campylobacter strains, specifically within the chicken meat production, a topic of serious importance.

High-Accuracy Long-Read Sequencing of Mycobacterium tuberculosis PSNK363 Isolated From the Democratic People's Republic of Korea

Long-read sequencing is a valuable technique for high-precision genome analysis. Despite the widespread use of the Mycobacterium tuberculosis H37Rv genome sequence as a reference for genetic variation analysis, its suitability for comparing clinical strains is limited. Therefore, we constructed the first known whole genome of a clinical M. tuberculosis strain, PSNK363, isolated from the Democratic People's Republic of Korea, using high-quality high-fidelity (HiFi) read sequencing and compared its genetic variations to those of H37Rv. PSNK363 was cultured to obtain genomic DNA, which was subjected to de novo whole-genome assembly using PacBio Sequel II with long-read HiFi sequencing. The sequences were compared to the reference genome H37Rv. HiFi long-read sequencing of M. tuberculosis PSNK363, with an accuracy of 99.99%, revealed a single circular chromosome of 4,422,110 bp, which is 10,578 bp longer than the H37Rv chromosome. The assembly had an average G + C content of 65.6%, 4079 protein-coding sequences, 53 tRNA genes, and 3 rRNA genes. Most genes (72.7%) were assigned as putative functions, whereas the remaining 27.3% were annotated as hypothetical. Comparison with H37Rv revealed a large inversion in the PSNK363 genome, which contains most of the deletion and insertion variants. M. tuberculosis PSNK363 had a longer genome sequence, more protein-coding genes, and a larger inversion region than H37Rv. High-accuracy whole-genome sequencing of PSNK363 holds the potential for enriching virulence databases and identifying informative loci for drug resistance analysis in M. tuberculosis isolates in the Democratic People's Republic of Korea.

Detoxome Capacity of the Adult Rumen Fluke Calicophoron daubneyi Extends into Its Secreted Extracellular Vesicles

Helminth parasites have long adapted to survive hostile host environments and can likely adapt against the chemical anthelmintic challenge. One proposed adaptation route is via Phase I and II xenobiotic metabolizing enzymes (XMEs). For successful Helminth pharmacotherapy discovery programs, a working understanding of Helminth-derived chemical detoxification, the Helminth detoxome, is a must. At present, the detoxome of a newly emerging Helminth parasite, the rumen fluke Calicophoron daubneyi, remains unexplored. Thus, a combined bioinformatics, sub-, and global-proteomic approach has been employed to examine the detoxome of adult C. daubneyi. Transcriptome analysis revealed a complement of Phase I (cytochrome P450s and monoamine oxygenases) and Phase II (glutathione transferases [GSTs] and sulfotransferases) XMEs. Affinity-led subproteomic exploration of the GSTs revealed six GST isoforms in adult rumen fluke (CdGST-Mu1-2, S1, and S3-5), with global approaches identifying additional GSTs (CdGST-O1-2, Z1, and S2) and a unique egg-specific variant (CdGST-S6). Examination of C. daubneyi extracellular vesicles revealed a GST profile replicating that of the adult with the absence of two isoforms (CdGST-S2 and S4), with an additional identification of a sulfotransferase. These data represent the first exploration into the complete rumen fluke detoxification capacity and will provide direction for future anthelmintic discovery programs.

sORFdb - a database for sORFs, small proteins, and small protein families in bacteria

Small proteins with fewer than 100, particularly fewer than 50, amino acids are still largely unexplored. Nonetheless, they represent an essential part of bacteria's often neglected genetic repertoire. In recent years, the development of ribosome profiling protocols has led to the detection of an increasing number of previously unknown small proteins. Despite this, they are overlooked in many cases by automated genome annotation pipelines, and often, no functional descriptions can be assigned due to a lack of known homologs. To understand and overcome these limitations, the current abundance of small proteins in existing databases was evaluated, and a new dedicated database for small proteins and their potential functions, called 'sORFdb', was created. To this end, small proteins were extracted from annotated bacterial genomes in the GenBank database. Subsequently, they were quality-filtered, compared, and complemented with proteins from Swiss-Prot, UniProt, and SmProt to ensure reliable identification and characterization of small proteins. Families of similar small proteins were created using bidirectional best BLAST hits followed by Markov clustering. Analysis of small proteins in public databases revealed that their number is still limited due to historical and technical constraints. Additionally, functional descriptions were often missing despite the presence of potential homologs. As expected, a taxonomic bias was evident in over-represented clinically relevant bacteria. This new and comprehensive database is accessible via a feature-rich website providing specialized search features for sORFs and small proteins of high quality. Additionally, small protein families with Hidden Markov Models and information on taxonomic distribution and other physicochemical properties are available. In conclusion, the novel small protein database sORFdb is a specialized, taxonomy-independent database that improves the findability and classification of sORFs, small proteins, and their functions in bacteria, thereby supporting their future detection and consistent annotation. All sORFdb data is freely accessible via https://sorfdb.computational.bio .

Genomic Resources for the Scuttle Fly Megaselia abdita: A Model Organism for Comparative Developmental Studies in Flies

The order Diptera (true flies) holds promise as a model taxon in evolutionary developmental biology due to the inclusion of the model organism, Drosophila melanogaster, and the ability to cost-effectively rear many species in laboratories. One of them, the scuttle fly Megaselia abdita (Phoridae) has been used in evolutionary developmental biology for 30 years and is an excellent phylogenetic intermediate between fruit flies and mosquitoes but remains underdeveloped in genomic resources. Here, we present a de novo chromosome-level assembly and annotation of M. abdita and transcriptomes of 9 embryonic and 4 postembryonic stages. We also compare 9 stage-matched embryonic transcriptomes between M. abdita and D. melanogaster. Our analysis of these resources reveals extensive chromosomal synteny with D. melanogaster, 28 orphan genes with embryo-specific expression including a novel F-box LRR gene in M. abdita, and conserved and diverged features of gene expression dynamics between M. abdita and D. melanogaster. Collectively, our results provide a new reference for studying the diversification of developmental processes in flies.

The HTIRDB: A resource containing a transcriptional atlas for 105 different tissues from each of seven species of domestic herbivore

Here, we describe the Herbivore Transcriptome Integrated Resource Database (HTIRDB, https://yanglab.hzau.edu.cn/HTIRDB#/). The HTIRDB comprises the self-generated transcriptomic data from 100 to 105 tissues from two female domestic herbivores from six species (cattle, donkey, goat, horse, rabbit, and sika deer) and two breeds of sheep, and an extra 28,710 related published datasets. The HTIRDB user-friendly interface provides tools and functionalities that facilitate the exploration of gene expression between tissues and species. The tools for comparative transcriptomics can be used to identify housekeeping genes, tissue-specific genes, species-specific genes, and species-conserved genes. To date, the HTIRDB is the most extensive transcriptome data resource for domestic herbivores that is freely available.

Conservation of symmetry breaking at the level of chromatin accessibility between fly species with unrelated anterior determinants

Establishing the anterior-posterior body axis is a fundamental process during embryogenesis, and the fruit fly, Drosophila melanogaster, provides one of the best-known case studies of this process. In Drosophila, localized mRNA of bicoid serves as anterior determinant (AD). Bicoid engages in a concentration-dependent competition with nucleosomes and initiates symmetry-breaking along the AP axis by promoting chromatin accessibility at the loci of transcription factor (TF) genes that are expressed in the anterior of the embryo. However, ADs of other fly species are unrelated and structurally distinct, and little is known about how they function. We addressed this question in the moth fly, Clogmia albipunctata, in which a maternally expressed transcript isoform of the pair-rule segmentation gene odd-paired is localized in the anterior egg and has been co-opted as AD. We provide a de novo assembly and annotation of the Clogmia genome and describe how knockdown of zelda and maternal odd-paired transcript affect chromatin accessibility and expression of TF-encoding loci. The results of these experiments suggest direct roles of Cal-Zld in opening and closing chromatin during nuclear cleavage cycles and show that Clogmia's maternal odd-paired activity promotes chromatin accessibility and anterior expression during the early phase of zygotic genome activation at Clogmia's homeobrain and sloppy-paired loci. We conclude that unrelated dipteran ADs initiate anterior-posterior axis-specification at the level of enhancer accessibility and that homeobrain and sloppy-paired homologs may serve a more widely conserved role in the initiation of anterior pattern formation given their early anterior expression and function in head development in other insects.

The prevalence of motility-related genes within the human oral microbiota

The human oral and nasal microbiota contains approximately 770 cultivable bacterial species. More than 2,000 genome sequences of these bacteria can be found in the expanded Human Oral Microbiome Database (eHOMD). We developed HOMDscrape, a freely available Python software tool to programmatically retrieve and process amino acid sequences and sequence identifiers from BLAST results acquired from the eHOMD website. Using the data obtained through HOMDscrape, the phylogeny of proteins involved in bacterial type 9 secretion system (T9SS)-driven gliding motility, flagellar motility, and type IV pilus-driven twitching motility was constructed. A comprehensive phylogenetic analysis was conducted for all components of the rotary T9SS, a machinery responsible for secreting various enzymes, virulence factors, and enabling bacterial gliding motility. Results revealed that the T9SS outer membrane β-barrel protein SprA of human oral bacteria underwent horizontal evolution. Overall, we catalog motile bacteria that inhabit the human oral microbiota and document their evolutionary connections. These results will serve as a guide for further studies exploring the impact of motility on the shaping of the human oral microbiota.IMPORTANCEThe human oral microbiota has been extensively studied, and many of the isolated bacteria have genome sequences stored on the human oral microbiome database (eHOMD). Spatial distribution and polymicrobial biofilms are observed in the oral microbiota, but little is understood on how they are influenced by motility. To bridge this gap, we developed a software tool to identify motile bacteria from eHOMD. The results enabled the cataloging of motile bacteria present in the oral microbiota but also provided insight into their evolutionary relationships. This information can guide future research to better understand how bacterial motility shapes the human oral microbiota.

20 years of the Bio-Analytic Resource for Plant Biology

The Bio-Analytic Resource for Plant Biology ('the BAR', at https://bar.utoronto.ca) is celebrating its 20th year in operation in 2025. The BAR encompasses and provides visualization tools for large 'omics data sets from plants. The BAR covers data from Arabidopsis, tomato, wheat, barley and 29 other plant species (with data for 2 others to be released soon). These data include nucleotide and protein sequence data, gene expression data, protein-protein and protein-DNA interactions, protein structures, subcellular localizations, and polymorphisms. The data are stored in more than 200 relational databases holding 186 GB of data and are presented to the researchers via web apps. These web apps provide data analysis and visualization tools. Some of the most popular tools are eFP ('electronic fluorescent pictograph') Browsers, ePlants and ThaleMine (an Arabidopsis-specific instance of InterMine). The BAR was designated a Global Core Biodata Resource in 2023. Like other GCBRs, the BAR has excellent operational stability, provides access without login requirement, and provides an API for researchers to be able to access BAR data programmatically. We present in this update a new overarching search tool called Gaia that permits easy access to all BAR data, powered by machine learning and artificial intelligence.

AMIR: a multi-omics data platform for Asteraceae plants genetics and breeding research

As the largest family of dicotyledon, the Asteraceae family comprises a variety of economically important crops, ornamental plants and numerous medicinal herbs. Advancements in genomics and transcriptomic have revolutionized research in Asteraceae species, generating extensive omics data that necessitate an efficient platform for data integration and analysis. However, existing databases face challenges in mining genes with specific functions and supporting cross-species studies. To address these gaps, we introduce the Asteraceae Multi-omics Information Resource (AMIR; https://yanglab.hzau.edu.cn/AMIR/), a multi-omics hub for the Asteraceae plant community. AMIR integrates diverse omics data from 74 species, encompassing 132 genomes, 4 408 432 genes annotated across seven different perspectives, 3897 transcriptome sequencing samples spanning 131 organs, tissues and stimuli, 42 765 290 unique variants and 15 662 metabolites genes. Leveraging these data, AMIR establishes the first pan-genome, comparative genomics and transcriptome system for the Asteraceae family. Furthermore, AMIR offers user-friendly tools designed to facilitate extensive customized bioinformatics analyses. Two case studies demonstrate AMIR's capability to provide rapid, reproducible and reliable analysis results. In summary, by integrating multi-omics data of Asteraceae species and developing powerful analytical tools, AMIR significantly advances functional genomics research and contributes to breeding practices of Asteraceae.

Identification of Six Novel Proteins Containing a ZP Module from Nemertean Species

During fertilization, a series of reactions between the eggs and spermatozoa proceed predominantly in a species-specific manner. The molecules mediating these species-specific reactions remain unknown except in a few organisms. In this study, we focused on two species belonging to the phylum Nemertea, Kulikovia alborostrata and K. fulva, and explored molecules involved in species-specific interactions between gametes. Orthologs of molecules known to be involved in species-specific reactions were not expressed in the ovaries of these two species. In contrast, we identified six novel proteins, named NeZPL1-NeZPL6, containing a ZP module. Among these, we found that NeZPL6 is located on the surface of an unfertilized egg and is suggested to be involved in its interaction with spermatozoa. Furthermore, we found an indel of three amino acids in the EGF-like domain of NeZPL6, which possibly confers species specificity to this interaction. Our results suggested the existence of a novel system for species recognition in animal gametes.

Gene loss in Antarctic icefish: evolutionary adaptations mimicking Fanconi Anemia?

BACKGROUND

The white-blooded Antarctic icefishes is a representative organism that survive under the stenothermal conditions of the Southern Ocean without the hemoglobin genes. To compensate for inefficient oxygen transport, distinct features such as increased heart size, greater blood volume, and reduced hematocrit density enhance the amount of dissolved oxygen and the velocity of blood flow.

RESULTS

Here, we investigated these unique characteristics by comparing high-quality genomic data between white-blooded and red-blooded fishes and identified the loss of FAAP20, which is implicated in anemia. Although the gene region containing FAAP20 is conserved in notothenioids as shown through collinear analysis, only remnants of FAAP20 persist in several icefish species. Additionally, we observed the loss of SOAT1, which plays a pivotal role in cholesterol metabolism, providing a clue for further investigations into the unique mitochondrial form of the icefish.

CONCLUSIONS

The loss of FAAP20, which is known to reduce erythrocyte counts under stress conditions in mice and humans, may provide a clue to understanding the genomic characteristics related to oxygen supply, such as low hematocrit, in Antarctic icefishes.

LukProt: A Database of Eukaryotic Predicted Proteins Designed for Investigations of Animal Origins

The origins and early evolution of animals are subjects with many outstanding questions. One problem faced by researchers trying to answer them is the absence of a comprehensive database with sequences from nonbilaterians. Publicly available data are plentiful but scattered and often not associated with proper metadata. A new database presented in this paper, LukProt, is an attempt at solving this issue. The database contains protein sequences obtained mostly from genomic, transcriptomic, and metagenomic studies and is an extension of EukProt (Richter DJ, Berney C, Strassert JFH, Poh Y-P, Herman EK, Muñoz-Gómez SA, Wideman JG, Burki F, de Vargas C. EukProt: a database of genome-scale predicted proteins across the diversity of eukaryotes. Peer Community J. 2022:2:e56. https://doi.org/10.24072/pcjournal.173). LukProt adopts the EukProt naming conventions and includes data from 216 additional animals. The database is associated with a taxonomic grouping (taxogroup) scheme suitable for studying early animal evolution. Minor updates to the database will contain species additions or metadata corrections, whereas major updates will synchronize LukProt to each new version of EukProt, and releases are permanently stored on Zenodo (https://doi.org/10.5281/zenodo.7089120). A BLAST server to search the database is available at: https://lukprot.hirszfeld.pl/. Users are invited to participate in maintaining and correcting LukProt. As it can be searched without downloading locally, the database aims to be a convenient resource not only for evolutionary biologists, but for the broader scientific community as well.

Comparison of genotypic features between two groups of antibiotic resistant Klebsiella pneumoniae clinical isolates obtained before and after the COVID-19 pandemic from Egypt

Klebsiella pneumoniae is a common pathogen capable of causing a wide range of infections. Antibiotic resistance complicates treatment of these infections significantly. We are comparing resistance levels and genotypes among two collections of K. pneumoniae clinical isolates from Alexandria Main University Hospital (AMUH). We used disc diffusion and Minimum Inhibitory Concentration (MIC) by microbroth dilution to assess resistance levels and performed whole genome sequencing (WGS) to describe multilocus sequence types (MLST) and resistance gene presence. Among a collection of 56 K. pneumoniae clinical isolates (19 from 2019 to 37 from 2021), multidrug resistance (MDR) was 33% and 10%, extended drug resistance (XDR) was 24% and 46% and pan-drug resistance (PDR) was 43% and 43%, respectively. We identified 15 MLST STs including two novel types (ST-6118 and ST-6119 ). ST-101 and ST-383 were common between the two collections; ST-101 was the most common genotype in 2019 (28.6%) and ST-147 was most common in 2021 (25%). Ampicillin/sulbactam, amikacin, cefepime, ceftriaxone and ertapenem MICs were significantly higher in 2021. Prevalence of aph(3') - Ia, aph(3')-VI, mphA was significantly higher in 2021. The increasing resistance levels and the persistence of some MDR/XDR genotypes is concerning. Understanding mechanisms of resistance will inform infection control and antimicrobial stewardship plans to prevent evolution and spread of XDR and PDR strains.

SesamumGDB: a comprehensive platform for Sesamum genetics and genomics analysis

Sesame (Sesamum indicum L., 2n = 26) is a crucial oilseed crop cultivated worldwide. The ancient evolutionary position of the Sesamum genus highlights its value for genomics and molecular genetics research among the angiosperms of other genera. However, Sesamum is considered a small orphan genus with only a few genomic databases for cultivated sesame to date. The urgent need to construct comprehensive, curated genome databases that include genus-specific gene resources for both cultivated and wild Sesamum species is being recognized. In response, we developed Sesamum Genomics Database (SesamumGDB), a user-friendly genomic database that integrates extensive genomic resources from two cultivated sesame varieties (S. indicum) and seven wild Sesamum species, covering all three chromosome groups (2n = 26, 32, and 64). This database showcases a total of 352 471 genes, including 6026 related to lipid metabolism and 17 625 transcription factors within Sesamum. Equipped with an array of bioinformatics tools such as BLAST (basic local alignment search tool) and JBrowse (the Javascript browser), SesamumGDB facilitates data downloading, screening, visualization, and analysis. As the first centralized Sesamum genome database, SesamumGDB offers extensive insights into the genomics and genetics of sesame, potentially enhancing the molecular breeding of sesame and other oilseed crops in the future. Database URL: http://www.sgbdb.com/sgdb/.

Gliding motility of the diatom Craspedostauros australis coincides with the intracellular movement of raphid-specific myosins

Raphid diatoms are one of the few eukaryotes capable of gliding motility, which is remarkably fast and allows for quasi-instantaneous directional reversals. Besides other mechanistic models, it has been suggested that an actomyosin system provides the force for diatom gliding. However, in vivo data on the dynamics of actin and myosin in diatoms are lacking. In this study, we demonstrate that the raphe-associated actin bundles required for diatom movement do not exhibit a directional turnover of subunits and thus their dynamics do not contribute directly to force generation. By phylogenomic analysis, we identified four raphid diatom-specific myosins in Craspedostauros australis (CaMyo51A-D) and investigated their in vivo localization and dynamics through GFP-tagging. Only CaMyo51B-D but not CaMyo51A exhibited coordinated movement during gliding, consistent with a role in force generation. The characterization of raphid diatom-specific myosins lays the foundation for unraveling the molecular mechanisms that underlie the gliding motility of diatoms.

Fundamental Patterns of Structural Evolution Revealed by Chromosome-Length Genomes of Cactophilic Drosophila

A thorough understanding of adaptation and speciation requires model organisms with both a history of ecological and phenotypic study as well as a complete set of genomic resources. In particular, high-quality genome assemblies of ecological model organisms are needed to assess the evolution of genome structure and its role in adaptation and speciation. Here, we generate new genomes of cactophilic Drosophila, a crucial model clade for understanding speciation and ecological adaptation in xeric environments. We generated chromosome-level genome assemblies and complete annotations for seven populations across Drosophila mojavensis, Drosophila arizonae, and Drosophila navojoa. We use these data first to establish the most robust phylogeny for this clade to date, and to assess patterns of molecular evolution across the phylogeny, showing concordance with a priori hypotheses regarding adaptive genes in this system. We then show that structural evolution occurs at constant rate across the phylogeny, varies by chromosome, and is correlated with molecular evolution. These results advance the understanding of the D. mojavensis clade by demonstrating core evolutionary genetic patterns and integrating those patterns to generate new gene-level hypotheses regarding adaptation. Our data are presented in a new public database (cactusflybase.arizona.edu), providing one of the most in-depth resources for the analysis of inter- and intraspecific evolutionary genomic data. Furthermore, we anticipate that the patterns of structural evolution identified here will serve as a baseline for future comparative studies to identify the factors that influence the evolution of genome structure across taxa.

The vacuolar sulfate transporter PsSULTR4 is a key determinant of seed yield and protein composition in pea

Pea is a grain legume crop with a high potential to accelerate the food transition due to its high seed protein content and relatively well-balanced amino acid composition. The critical role of external sulfur (S) supply in determining seed yield and seed quality in pea makes it essential to understand the impact of whole plant S management on the trade-off between these two traits. Here, we investigated the physiological relevance of vacuolar sulfate remobilization by targeting PsSULTR4, the only pea sulfate transporter showing substantial similarity to the vacuolar sulfate exporter AtSULTR4;1. Five mutations in PsSULTR4 were identified by TILLING (Targeting Induced Local Lesions IN Genomes), two of which, a loss of function (W78*) and a missense (E568K), significantly decreased seed yield under S deprivation. We demonstrate that PsSULTR4 triggers S distribution from source tissues, especially lower leaves, to reproductive organs to maintain seed yield under S deficiency. Under sufficient S supply, sultr4 seeds display lower levels of the S-rich storage protein PA1 at maturity. They also overaccumulate sulfate in the endosperm at the onset of seed filling. These findings uncover a role of PsSULTR4 in the remobilization of vacuolar sulfate during embryo development, allowing the efficient synthesis of S-rich proteins. Our study uncovers that PsSULTR4 functions (i) in source tissues to remobilize stored vacuolar sulfate for seed production under low S availability and (ii) in developing seeds well supplied with S to fine-tune sulfate remobilization from the endosperm as a critical control point for storage activities in the embryo.

An introduction to comparative genomics, EukProt, and the reciprocal best hit (RBH) method for bench biologists: Ancestral phosphorylation of Tom22 in eukaryotes as a case study

Comparative genomics is a useful approach for hypothesis generation for future functional investigations at the bench. However, most bench biologists shy away from computational methods. Here we reintroduce the simple but extremely effective Reciprocal Best Hit method for inferring protein orthologues. Because taxon set delimitation is perhaps the most important step in comparative genomics, we introduce The Comparative Set, a taxonomically representative subset of EukProt, a comprehensive eukaryotic predicted proteome database. After introducing the basic methods, we provide a step-by-step guide, including screen shots, for a case study on collecting Tom22 sequences from diverse eukaryotes. As an example of possible downstream analyses, we show that Tom22 proteins from diverse eukaryotes are likely regulated by conserved kinases at several sites. Though the sites evolve quickly, the processes and functions involved are likely ancestral and conserved across many eukaryotes.

SapBase: A central portal for functional and comparative genomics of Sapindaceae species

The Sapindaceae family, encompassing a wide range of plant forms such as herbs, vines, shrubs, and trees, is widely distributed across tropical and subtropical regions. This family includes economically important crops like litchi, longan, rambutan, and ackee. With the wide application of genomic technologies in recent years, several Sapindaceae plant genomes have been decoded, leading to an accumulation of substantial omics data in this field. This surge in data highlights the pressing need for a unified genomic data center capable of storing, sharing, and analyzing these data. Here, we introduced SapBase, that is, the Sapindaceae Genome Database. SapBase houses seven published plant genomes alongside their corresponding gene structure and functional annotations, small RNA annotations, gene expression profiles, gene pathways, and synteny block information. It offers user-friendly features for gene information mining, co-expression analysis, and inter-species comparative genomic analysis. Furthermore, we showcased SapBase's extensive capacities through a detailed bioinformatic analysis of a MYB gene in litchi. Thus, SapBase could serve as an integrative genomic resource and analysis platform for the scientific exploration of Sapinaceae species and their comparative studies with other plants.

OysterDB: A Genome Database for Ostreidae

The molluscan family Ostreidae, commonly known as oysters, is an important molluscan group due to its economic and ecological importance. In recent years, an abundance of genomic data of Ostreidae species has been generated and available in public domain. However, there is still a lack of a high-efficiency database platform to store and distribute these data with comprehensive tools. In this study, we developed an oyster genome database (OysterDB) to consolidate oyster genomic data. This database includes eight oyster genomes and 208,923 protein-coding gene annotations. Bioinformatic tools, such as BLAST and JBrowse, are integrated into the database to provide a user-friendly platform for homologous sequence searching, visualization of genomes, and screen for candidate gene information. Moreover, OysterDB will be continuously updated with ever-growing oyster genomic resources and facilitate future studies for comparative and functional genomic analysis of oysters ( http://oysterdb.com.cn/ ).

Annotation Vocabulary (Might Be) All You Need

Protein Language Models (pLMs) have revolutionized the computational modeling of protein systems, building numerical embeddings that are centered around structural features. To enhance the breadth of biochemically relevant properties available in protein embeddings, we engineered the Annotation Vocabulary, a transformer readable language of protein properties defined by structured ontologies. We trained Annotation Transformers (AT) from the ground up to recover masked protein property inputs without reference to amino acid sequences, building a new numerical feature space on protein descriptions alone. We leverage AT representations in various model architectures, for both protein representation and generation. To showcase the merit of Annotation Vocabulary integration, we performed 515 diverse downstream experiments. Using a novel loss function and only $3 in commercial compute, our premier representation model CAMP produces state-of-the-art embeddings for five out of 15 common datasets with competitive performance on the rest; highlighting the computational efficiency of latent space curation with Annotation Vocabulary. To standardize the comparison of de novo generated protein sequences, we suggest a new sequence alignment-based score that is more flexible and biologically relevant than traditional language modeling metrics. Our generative model, GSM, produces high alignment scores from annotation-only prompts with a BERT-like generation scheme. Of particular note, many GSM hallucinations return statistically significant BLAST hits, where enrichment analysis shows properties matching the annotation prompt - even when the ground truth has low sequence identity to the entire training set. Overall, the Annotation Vocabulary toolbox presents a promising pathway to replace traditional tokens with members of ontologies and knowledge graphs, enhancing transformer models in specific domains. The concise, accurate, and efficient descriptions of proteins by the Annotation Vocabulary offers a novel way to build numerical representations of proteins for protein annotation and design.

Genetic mapping of stripe rust resistance in a geographically diverse barley collection and selected biparental populations

Barley stripe or yellow rust (BYR) caused by Puccinia striiformis f. sp. hordei (Psh) is a significant constraint to barley production. The disease is best controlled by genetic resistance, which is considered the most economical and sustainable component of integrated disease management. In this study, we assessed the diversity of resistance to Psh in a panel of international barley genotypes (n = 266) under multiple disease environments (Ecuador, India, and Mexico) using genome-wide association studies (GWASs). Four quantitative trait loci (QTLs) (three on chromosome 1H and one on 7H) associated with resistance to Psh were identified. The QTLs were validated by mapping resistance to Psh in five biparental populations, which detected key genomic regions on chromosomes 1H (populations Pompadour/Zhoungdamei, Pompadour/Zug161, and CI9214/Baudin), 3H (Ricardo/Gus), and 7H (Fumai8/Baronesse). The QTL RpshQ.GWA.1H.1 detected by GWAS and RpshQ.Bau.1H detected using biparental mapping populations co-located were the most consistent and stable across environments and are likely the same resistance region. RpshQ.Bau.1H was saturated using population CI9214/Baudin by enriching the target region, which placed the resistance locus between 7.9 and 8.1 Mbp (flanked by markers sun_B1H_03, 0.7 cM proximal to Rpsh_1H and sun_B1H_KASP_02, 3.2 cM distal on 1HS) in the Morex reference genome v.2. A Kompetitive Allele Specific PCR (KASP) marker sun_B1H_KASP_01 that co-segregated for RpshQ.Bau.1H was developed. The marker was validated on 50 Australian barley cultivars, showing well-defined allelic discrimination and presence in six genotypes (Baudin, Fathom, Flagship, Grout, Sakurastar, and Shepherd). This marker can be used for reliable marker-assisted selection and pyramiding of resistance to Psh and in diversifying the genetic base of resistance to stripe rust.

A Quantitative Trait Locus with a Major Effect on Root-Lesion Nematode Resistance in Barley

Although the root-lesion nematode Pratylenchus thornei is known to affect barley (Hordeum vulgare L.), there have been no reports on the genetic control of P. thornei resistance in barley. In this research, P. thornei resistance was assessed for a panel of 46 barley mapping parents and for two mapping populations (Arapiles/Franklin and Denar/Baudin). With both populations, a highly significant quantitative trait locus (QTL) was mapped at the same position on the long arm of chromosome 7H. Single-nucleotide polymorphisms (SNPs) in this region were anchored to an RGT Planet pan-genome assembly and assayed on the mapping parents and other barley varieties. The results indicate that Arapiles, Denar, RGT Planet and several other varieties likely have the same resistance gene on chromosome 7H. Marker assays reported here could be used to select for P. thornei resistance in barley breeding. Analysis of existing barley pan-genomic and pan-transcriptomic data provided a list of candidate genes along with information on the expression and differential expression of some of those genes in barley root tissue. Further research is required to identify a specific barley gene that affects root-lesion nematode resistance.

HortGenome Search Engine, a universal genomic search engine for horticultural crops

Horticultural crops comprising fruit, vegetable, ornamental, beverage, medicinal and aromatic plants play essential roles in food security and human health, as well as landscaping. With the advances of sequencing technologies, genomes for hundreds of horticultural crops have been deciphered in recent years, providing a basis for understanding gene functions and regulatory networks and for the improvement of horticultural crops. However, these valuable genomic data are scattered in warehouses with various complex searching and displaying strategies, which increases learning and usage costs and makes comparative and functional genomic analyses across different horticultural crops very challenging. To this end, we have developed a lightweight universal search engine, HortGenome Search Engine (HSE; http://hort.moilab.net), which allows for the querying of genes, functional annotations, protein domains, homologs, and other gene-related functional information of more than 500 horticultural crops. In addition, four commonly used tools, including 'BLAST', 'Batch Query', 'Enrichment analysis', and 'Synteny Viewer' have been developed for efficient mining and analysis of these genomic data.

A mitochondrial carrier transports glycolytic intermediates to link cytosolic and mitochondrial glycolysis in the human gut parasite Blastocystis

Stramenopiles form a clade of diverse eukaryotic organisms, including multicellular algae, the fish and plant pathogenic oomycetes, such as the potato blight Phytophthora, and the human intestinal protozoan Blastocystis. In most eukaryotes, glycolysis is a strictly cytosolic metabolic pathway that converts glucose to pyruvate, resulting in the production of NADH and ATP (Adenosine triphosphate). In contrast, stramenopiles have a branched glycolysis in which the enzymes of the pay-off phase are located in both the cytosol and the mitochondrial matrix. Here, we identify a mitochondrial carrier in Blastocystis that can transport glycolytic intermediates, such as dihydroxyacetone phosphate and glyceraldehyde-3-phosphate, across the mitochondrial inner membrane, linking the cytosolic and mitochondrial branches of glycolysis. Comparative analyses with the phylogenetically related human mitochondrial oxoglutarate carrier (SLC25A11) and dicarboxylate carrier (SLC25A10) show that the glycolytic intermediate carrier has lost its ability to transport the canonical substrates malate and oxoglutarate. Blastocystis lacks several key components of oxidative phosphorylation required for the generation of mitochondrial ATP, such as complexes III and IV, ATP synthase, and ADP/ATP carriers. The presence of the glycolytic pay-off phase in the mitochondrial matrix generates ATP, which powers energy-requiring processes, such as macromolecular synthesis, as well as NADH, used by mitochondrial complex I to generate a proton motive force to drive the import of proteins and molecules. Given its unique substrate specificity and central role in carbon and energy metabolism, the carrier for glycolytic intermediates identified here represents a specific drug and pesticide target against stramenopile pathogens, which are of great economic importance.

Updates to the Alliance of Genome Resources central infrastructure

The Alliance of Genome Resources (Alliance) is an extensible coalition of knowledgebases focused on the genetics and genomics of intensively studied model organisms. The Alliance is organized as individual knowledge centers with strong connections to their research communities and a centralized software infrastructure, discussed here. Model organisms currently represented in the Alliance are budding yeast, Caenorhabditis elegans, Drosophila, zebrafish, frog, laboratory mouse, laboratory rat, and the Gene Ontology Consortium. The project is in a rapid development phase to harmonize knowledge, store it, analyze it, and present it to the community through a web portal, direct downloads, and application programming interfaces (APIs). Here, we focus on developments over the last 2 years. Specifically, we added and enhanced tools for browsing the genome (JBrowse), downloading sequences, mining complex data (AllianceMine), visualizing pathways, full-text searching of the literature (Textpresso), and sequence similarity searching (SequenceServer). We enhanced existing interactive data tables and added an interactive table of paralogs to complement our representation of orthology. To support individual model organism communities, we implemented species-specific "landing pages" and will add disease-specific portals soon; in addition, we support a common community forum implemented in Discourse software. We describe our progress toward a central persistent database to support curation, the data modeling that underpins harmonization, and progress toward a state-of-the-art literature curation system with integrated artificial intelligence and machine learning (AI/ML).

DELLA proteins recruit the Mediator complex subunit MED15 to coactivate transcription in land plants

DELLA proteins are negative regulators of the gibberellin response pathway in angiosperms, acting as central hubs that interact with hundreds of transcription factors (TFs) and regulators to modulate their activities. While the mechanism of TF sequestration by DELLAs to prevent DNA binding to downstream targets has been extensively documented, the mechanism that allows them to act as coactivators remains to be understood. Here, we demonstrate that DELLAs directly recruit the Mediator complex to specific loci in Arabidopsis, facilitating transcription. This recruitment involves DELLA amino-terminal domain and the conserved MED15 KIX domain. Accordingly, partial loss of MED15 function mainly disrupted processes known to rely on DELLA coactivation capacity, including cytokinin-dependent regulation of meristem function and skotomorphogenic response, gibberellin metabolism feedback, and flavonol production. We have also found that the single DELLA protein in the liverwort Marchantia polymorpha is capable of recruiting MpMED15 subunits, contributing to transcriptional coactivation. The conservation of Mediator-dependent transcriptional coactivation by DELLA between Arabidopsis and Marchantia implies that this mechanism is intrinsic to the emergence of DELLA in the last common ancestor of land plants.

Zygosity-based sex determination in a butterfly drives hypervariability of Masculinizer

Nature has devised many ways of producing males and females. Here, we report on a previously undescribed mechanism for Lepidoptera that functions without a female-specific gene. The number of alleles or allele heterozygosity in a single Z-linked gene (BaMasc) is the primary sex-determining switch in Bicyclus anynana butterflies. Embryos carrying a single BaMasc allele develop into WZ (or Z0) females, those carrying two distinct alleles develop into ZZ males, while (ZZ) homozygotes initiate female development, have mismatched dosage compensation, and die as embryos. Consequently, selection against homozygotes has favored the evolution of spectacular allelic diversity: 205 different coding sequences of BaMasc were detected in a sample of 246 females. The structural similarity of a hypervariable region (HVR) in BaMasc to the HVR in Apis mellifera csd suggests molecular convergence between deeply diverged insect lineages. Our discovery of this primary switch highlights the fascinating diversity of sex-determining mechanisms and underlying evolutionary drivers.

Analyzing the safety of the parasiticide fungus Mucor circinelloides: first insights on its virulence profile and interactions with the avian gut microbial community

Parasiticide fungi are considered an accurate, sustainable, and safe solution for the biocontrol of animal gastrointestinal (GI) parasites. This research provides an initial characterization of the virulence of the native parasiticide fungus Mucor circinelloides (FMV-FR1) and an assessment of its impact on birds' gut microbes. The genome of this fungus was sequenced to identify the genes coding for virulence factors. Also, this fungus was checked for the phenotypic expression of proteinase, lecithinase, DNase, gelatinase, hemolysin, and biofilm production. Finally, an in vivo trial was developed based on feeding M. circinelloides spores to laying hens and peacocks three times a week. Bird feces were collected for 3 months, with total genomic DNA being extracted and subjected to long-read 16S and 25S-28S sequencing. Genes coding for an iron permease (FTR1), iron receptors (FOB1 and FOB2), ADP-ribosylation factors (ARFs) (ARF2 and ARF6), and a GTPase (CDC42) were identified in this M. circinelloides genome. Also, this fungus was positive only for lecithinase activity. The field trial revealed a fecal microbiome dominated by Firmicutes and Proteobacteria in laying hens, and Firmicutes and Bacteroidetes in peacocks, whereas the fecal mycobiome of both bird species was mainly composed of Ascomycetes and Basidiomycetes fungi. Bacterial and fungal alpha-diversities did not differ between sampling time points after M. circinelloides administrations (P = 0.62 and P = 0.15, respectively). Although findings from this research suggest the lack of virulence of this M. circinelloides parasiticide isolate, more complementary in vitro and in vivo research is needed to conclude about the safety of its administration to birds, aiming at controlling their GI parasites.IMPORTANCEA previous study revealed that the native Mucor circinelloides isolate (FMV-FR1) can develop parasiticide activity toward coccidia oocysts, one of the most pathogenic GI parasites in birds. However, ensuring its safety for birds is of utmost importance, namely by studying its virulence profile and potential effect on commensal gut microbes. This initial study revealed that although this M. circinelloides isolate had genes coding for four types of virulence factors-iron permease, iron receptors, ADP-ribosylation factors, and GTPase-and only expressed phenotypically the enzyme lecithinase, the administration of its spores to laying hens and peacocks did not interfere with the abundances and diversities of their gut commensal bacteria and fungi. Although overall results suggest the lack of virulence of this M. circinelloides isolate, more complementary research is needed to conclude about the safety of its administration to birds in the scope of parasite biocontrol programs.

Phylogenomics and genetic analysis of solvent-producing Clostridium species

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Unique spatially and temporary-regulated/sex-specific expression of a long ncRNA, Nb-1, suggesting its pleiotropic functions associated with honey bee lifecycle

Honey bees are social insects, and each colony member has unique morphological and physiological traits associated with their social tasks. Previously, we identified a long non-coding RNA from honey bees, termed Nb-1, whose expression in the brain decreases associated with the age-polyethism of workers and is detected in some neurosecretory cells and octopaminergic neurons, suggesting its role in the regulation of worker labor transition. Herein, we investigated its spatially and temporary-regulated/sex-specific expression. Nb-1 was expressed as an abundant maternal RNA during oogenesis and embryogenesis in both sexes. In addition, Nb-1 was expressed preferentially in the proliferating neuroblasts of the mushroom bodies (a higher-order center of the insect brain) in the pupal brains, suggesting its role in embryogenesis and mushroom body development. On the contrary, Nb-1 was expressed in a drone-specific manner in the pupal and adult retina, suggesting its role in the drone visual development and/or sense. Subcellular localization of Nb-1 in the brain during development differed depending on the cell type. Considering that Nb-1 is conserved only in Apidae, our findings suggest that Nb-1 potentially has pleiotropic functions in the expression of multiple developmental, behavioral, and physiological traits, which are closely associated with the honey bee lifecycle.

Chromosomal DNA sequences of the Pacific saury genome: versatile resources for fishery science and comparative biology

Pacific saury (Cololabis saira) is a commercially important small pelagic fish species in Asia. In this study, we conducted the first-ever whole genome sequencing of this species, with single molecule, real-time (SMRT) sequencing technology. The obtained high-fidelity (HiFi) long-read sequence data, which amount to ~30-folds of its haploid genome size that was measured with quantitative PCR (1.17 Gb), were assembled into contigs. Scaffolding with Hi-C reads yielded a whole genome assembly containing 24 chromosome-scale sequences, with a scaffold N50 length of 47.7 Mb. Screening of repetitive elements including telomeric repeats was performed to characterize possible factors that need to be resolved towards 'telomere-to-telomere' sequencing. The larger genome size than in medaka, a close relative in Beloniformes, is at least partly explained by larger repetitive element quantity, which is reflected in more abundant tRNAs, in the Pacific saury genome. Protein-coding regions were predicted using transcriptome data, which resulted in 22,274 components. Retrieval of Pacific saury homologs of aquaporin (AQP) genes known from other teleost fishes validated high completeness and continuity of the genome assembly. These resources are available at https://treethinkers.nig.ac.jp/saira/ and will assist various molecular-level studies in fishery science and comparative biology.

The LaCLE35 peptide modifies rootlet density and length in cluster roots of white lupin

White lupin (lupinus albus L.) forms special bottlebrush-like root structures called cluster roots (CR) when phosphorus is low, to remobilise sparingly soluble phosphates in the soil. The molecular mechanisms that control the CR formation remain unknown. Root development in other plants is regulated by CLE  (CLAVATA3/ EMBRYO SURROUNDING REGION (ESR)-RELATED) peptides, which provide more precise control mechanisms than common phytohormones. This makes these peptides interesting candidates to be involved in CR formation, where fine tuning to environmental factors is required. In this study we present an analysis of CLE peptides in white lupin. The peptides LaCLE35 (RGVHy PSGANPLHN) and LaCLE55 (RRVHy PSCHy PDPLHN) reduced root growth and altered CR in hydroponically cultured white lupins. We demonstrate that rootlet density and rootlet length were locally, but not systemically, impaired by exogenously applied CLE35. The peptide was identified in the xylem sap. The inhibitory effect of CLE35 on root growth was attributed to arrested cell elongation in root tips. Taken together, CLE peptides affect both rootlet density and rootlet length, which are two critical factors for CR formation, and may be involved in fine tuning this peculiar root structure that is present in a few crops and many Proteaceae species, under low phosphorus availability.

Genetic structure, UV-vision, wing coloration and size coincide with colour polymorphism in Fabriciana adippe butterflies

Colour polymorphisms have long served as model systems in evolutionary studies and continue to inform about processes involved in the origin and dynamics of biodiversity. Modern sequencing tools allow for evaluating whether phenotypic differences between morphs reflect genetic differentiation rather than developmental plasticity, and for investigating whether polymorphisms represent intermediate stages of diversification towards speciation. We investigated phenotypic and genetic differentiation between two colour morphs of the butterfly Fabriciana adippe using a combination of ddRAD-sequencing and comparisons of body size, colour patterns and optical properties of bright wing spots. The silvery-spotted adippe form had larger and darker wings and reflected UV light, while the yellow cleodoxa form displayed more green scales and reflected very little UV, showcasing that they constitute distinct and alternative integrated phenotypes. Genomic analyses revealed genetic structuring according to source population, and to colour morph, suggesting that the phenotypic differentiation reflects evolutionary modifications. We report 17 outlier loci associated with colour morph, including ultraviolet-sensitive visual pigment (UVRh1), which is associated with intraspecific communication and mate choice in butterflies. Together with the demonstration that the wings of the adippe (but essentially not the cleodoxa) morph reflect UV light, that UV reflectance is higher in females than males and that morphs differ in wing size, this suggests that these colour morphs might represent genetically integrated phenotypes, possibly adapted to different microhabitats. We propose that non-random mating might contribute to the differentiation and maintenance of the polymorphism.

SoyMD: a platform combining multi-omics data with various tools for soybean research and breeding

Advanced multi-omics technologies offer much information that can uncover the regulatory mechanisms from genotype to phenotype. In soybean, numerous multi-omics databases have been published. Although they cover multiple omics, there are still limitations when it comes to the types and scales of omics datasets and analysis methods utilized. This study aims to address these limitations by collecting and integrating a comprehensive set of multi-omics datasets. This includes 38 genomes, transcriptomes from 435 tissue samples, 125 phenotypes from 6686 accessions, epigenome data involving histone modification, transcription factor binding, chromosomal accessibility and chromosomal interaction, as well as genetic variation data from 24 501 soybean accessions. Then, common analysis pipelines and statistical methods were applied to mine information from these multi-omics datasets, resulting in the successful establishment of a user-friendly multi-omics database called SoyMD (https://yanglab.hzau.edu.cn/SoyMD/#/). SoyMD provides researchers with efficient query options and analysis tools, allowing them to swiftly access relevant omics information and conduct comprehensive multi-omics data analyses. Another notable feature of SoyMD is its capability to facilitate the analysis of candidate genes, as demonstrated in the case study on seed oil content. This highlights the immense potential of SoyMD in soybean genetic breeding and functional genomics research.

Foxn1 is not essential for T-cell development in teleosts

In mammals, T-cell development depends on the activity of the Foxn1 transcription factor in the thymic epithelium; mutations in the vertebrate-specific Foxn1 gene are associated with profound T-cell lymphopenia and fatal immunodeficiency. Here, we examined the extent of T-cell development in teleosts lacking a functional foxn1 gene. In zebrafish carrying a deleterious internal deletion of foxn1, reduced but robust lymphopoietic activity is maintained in the mutant thymus. Moreover, pseudogenization or loss of foxn1 in the genomes of deep-sea anglerfishes is independent of the presence or absence of the canonical signatures of the T-cell lineage. Thus, in contrast to the situation in mammals, the teleost thymus can support foxn1-independent lymphopoiesis, most likely through the activity of the Foxn4, an ancient metazoan paralog of Foxn1. Our results imply that during the early stages of vertebrate evolution, genetic control of thymopoiesis was functionally redundant and thus robust; in mammals, the genetic network was reorganized to become uniquely dependent on the FOXN1 transcription factor.

Lymphocyte pathway analysis using naturally lymphocyte-deficient fish

Comparative phylogenetic analyses are of potential value to establish the essential components of genetic networks underlying physiological traits. For species that naturally lack particular lymphocyte lineages, we show here that this strategy readily distinguishes trait-specific actors from pleiotropic components of the genetic network governing lymphocyte differentiation. Previously, three of the four members of the DNA polymerase X family have been implicated in the junctional diversification process during the somatic assembly of antigen receptors. Our phylogenetic analysis indicates that the presence of terminal deoxynucleotidyl transferase is strictly associated with the facility of V(D)J recombination, whereas PolL and PolM genes are retained even in species lacking Rag-mediated somatic diversification of antigen receptor genes.