Composition-based statistics and translated nucleotide searches: improving the TBLASTN module of BLAST

Microprotein-encoding RNA regulation in cells treated with pro-inflammatory and pro-fibrotic stimuli

BACKGROUND

Recent analysis of the human proteome via proteogenomics and ribosome profiling of the transcriptome revealed the existence of thousands of previously unannotated microprotein-coding small open reading frames (smORFs). Most functional microproteins were chosen for characterization because of their evolutionary conservation. However, one example of a non-conserved immunomodulatory microprotein in mice suggests that strict sequence conservation misses some intriguing microproteins.

RESULTS

We examine the ability of gene regulation to identify human microproteins with potential roles in inflammation or fibrosis of the intestine. To do this, we collected ribosome profiling data of intestinal cell lines and peripheral blood mononuclear cells and used gene expression of microprotein-encoding transcripts to identify strongly regulated microproteins, including several examples of microproteins that are only conserved with primates.

CONCLUSION

This approach reveals a number of new microproteins worthy of additional functional characterization and provides a dataset that can be queried in different ways to find additional gut microproteins of interest.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

Two high quality chromosome-scale genome assemblies of female and male silver pomfret (Pampus argenteus)

Pampus argenteus is a highly commercial marine fish whose population is declining sharply. Here, we generated a female P. argenteus genome, spanning 536.33 Mb with contig N50 of 1.79 Mb; 24070 genes (99.50% of 24,182) were functionally annotated. To improve quality of it, we assembled a 553.79 Mb genome of male fish with contig N50 of 24.75 Mb through HiFi and ultra-long ONT sequence technologies; 550.82 Mb were anchored onto 24 gap-free chromosomes; 22,892 genes (98.1% of 23,346) were functionally annotated; the QV value was 51.55 with 98.9% of BUSCO and 99.39% coverage of Illumina reads. Finally, we compared this genome with previous published one, revealing 37,301 SVs. 52.82 Mb and 18.05 Mb SDs were characterized in our and published assemblies, respectively, and 48.96 Mb PURs were constructed. Thus, this genome assembly exhibits excellent completeness, continuity and accuracy comparing to the published one, which can be current preferred reference genome. Overall, these works help aquaculture and wild resources recovery of P. argenteus and provide a valuable genetic resource for study.

A unique cell division protein critical for the assembly of the bacterial divisome

Identification of unique essential bacterial genes is important for not only the understanding of their cell biology but also the development of new antimicrobials. Here, we report a previously unrecognized core component of the Acinetobacter baumannii divisome. Our results reveal that the protein, termed Aeg1 interacts with multiple cell division proteins, including FtsN, which is required for components of the divisome to localize to the midcell. We demonstrate that the FtsAE202K and FtsBE65A mutants effectively bypassed the need of Aeg1 by A. baumannii, as did the activation variants FtsWM254I and FtsWS274G. Our results suggest that Aeg1 is a cell division protein that arrives at the division site to initiate cell division by recruiting FtsN, which activates FtsQLB and FtsA to induce the septal peptidoglycan synthase FtsWI. The discovery of the new essential cell division protein has provided a new target for the development of antibacterial agents.

Chromosome-level genome provides insights into evolution and diving adaptability in the vulnerable common pochard (Aythya ferina)

The common pochard (Aythya ferina) is a freshwater diving duck found in the Palearctic region that has been classified as vulnerable by the IUCN due to continuous and rapid population declines across their distribution. To gain a better understanding of its genetic mechanism of adaptive evolution, we successfully sequenced and assembled the first high-quality chromosome-level genome of A. ferina using Illumina, Nanopore and Hi-C sequencing technologies. A total assembly length of 1,130.78 Mbp was obtained, with over 98.81% (1,117.37Mbp) of sequence anchored to 35 pseudo-chromosomes. We predicted 17,232 protein-coding genes, 95.9% of which were functionally annotated. We identified 339 expanded and 937 contracted gene families in the genome of A. ferina, and detected 95 genes that have been positively selected. The significantly enriched Gene Ontology and enriched pathways were related to energy metabolism, immune, nervous, and sensory systems, suggests that these factors likely played an important role in its evolution. Importantly, we recovered signatures of positive selection on genes related to vasoconstriction that may be associated with thermoregulatory adaptations of A. ferina for underwater diving. Overall, the high-quality genome assembly and annotation in this study provides valuable genomic resources for ecological and evolutionary studies, as well as toward the conservation of A. ferina.

Chromosome genome assembly and annotation of Adzuki Bean (Vigna angularis)

Adzuki bean (Vigna angularis) is a significant dietary legume crop that is prevalent in East Asia. It also holds traditional medicinal importance in China. In this study, we report a high-quality, chromosome-level genome assembly of adzuki bean obtained by employing Illumina short-read sequencing, PacBio long-read sequencing, and Hi-C technology. The assembly spans 447.8 Mb, encompassing 96.32% of the estimated genome, with contig and scaffold N50 values of 16.5 and 41.0 Mb, respectively. More than 98.2% of the 1,614 BUSCO genes were fully identified, and 25,939 genes were annotated, with 98.23% of them being functionally identifiable. Vigna angularis was estimated to diverge successively from Vigna unguiculata and Vigna radiata about 15.3 and 8.7 million years ago (Ma), respectively. This chromosome-level reference genome of Vigna angularis provides a robust foundation for exploring the functional genomics and genome evolution of adzuki bean, thereby facilitating advancements in molecular breeding of adzuki bean.

Identification of Key Genes Controlling Sugar and Organic Acid Accumulation in Wampee Fruit (Clausena lansium) via Genome Assembly and Genome-wide Association Analysis

Wampee (Clausena lansium) is an economically significant subtropical fruit tree widely cultivated in Southern China. To provide high-quality genomic resources for C. lansium, we report a chromosome-level genome sequence for the "JinFeng" cultivar. The 297.1 Mb C. lansium genome contained nine chromosomes with a scaffold N50 of 29.2 Mb and encoded 23,468 protein-coding genes. Selective sweep analysis between sweet and sour C. lansium varieties and genome-wide association analysis identified 14 candidate genes putatively involved in sugar and acid accumulation. ClERF061, encoding an ethylene response factor, and ClSWEET7, encoding a Sugars Will Eventually be Exported Transporters (SWEET) family protein, were proposed as key regulators of the sweet and sour tastes of the wampee fruit. ClERF061 and ClSWEET7 overexpression in tomatoes increased the total sugar and acid content in fruits. ClSWEET7 promoter activation by ClERF061 was confirmed via Nicotiana benthamiana transient expression. Our study provides valuable genomic resources for C. lansium genetics and breeding.

Chromosome-scale and haplotype-resolved genome assembly of the autotetraploid Misgurnus anguillicaudatus

In nature, diploids and tetraploids are two common types of polyploid evolution. Misgurnus anguillicaudatus (mud loach) is a remarkable fish species that exhibits both diploid and tetraploid forms. However, reconstructing the four haplotypes of its autotetraploid genome remains unresolved. Here, we generated the first haplotype-resolved, chromosome-level genome of autotetraploid M. anguillicaudatus with a size of 4.76 Gb, contig N50 of 6.78 Mb, and scaffold N50 of 44.11 Mb. We identified approximately 2.9 Gb (61.03% of genome) of repetitive sequences and predicted 91,485 protein-coding genes. Moreover, allelic gene expression levels indicated the absence of significant dominant haplotypes within the autotetraploid loach genome. This genome will provide a valuable biological model for unraveling the mechanisms of polyploid formation and evolution, adaptation to environmental changes, and benefit for aquaculture applications and biodiversity conservation.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes

Spartina alterniflora is an exo-recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high-quality, chromosome-level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high-throughput chromosome conformation capture (Hi-C) technology and Illumina-based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole-genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species-specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA-seq analysis identified several ion transporter genes including the high-affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up-regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ-ECS, GSH-S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte.

A genome assembly of decaploid Houttuynia cordata provides insights into the evolution of Houttuynia and the biosynthesis of alkaloids

Houttuynia cordata Thunb., commonly known as yuxingcao in China, is known for its characteristic fishy smell and is widely recognized as an important herb and vegetable in many parts of Asia. However, the lack of genomic information on H. cordata limits the understanding of its population structure, genetic diversity, and biosynthesis of medicinal compounds. Here we used single-molecule sequencing, Illumina paired-end sequencing, and chromosome conformation capture technology to construct the first chromosome-scale decaploid H. cordata reference genome. The genome assembly was 2.63 Gb in size, with 1348 contigs and a contig N50 of 21.94 Mb further clustered and ordered into 88 pseudochromosomes based on Hi-C analysis. The results of genome evolution analysis showed that H. cordata underwent a whole-genome duplication (WGD) event ~17 million years ago, and an additional WGD event occurred 3.3 million years ago, which may be the main factor leading to the high abundance of multiple copies of orthologous genes. Here, transcriptome sequencing across five different tissues revealed significant expansion and distinct expression patterns of key gene families, such as l-amino acid/l-tryptophan decarboxylase and strictosidine synthase, which are essential for the biosynthesis of isoquinoline and indole alkaloids, along with the identification of genes such as TTM3, which is critical for root development. This study constructed the first decaploid medicinal plant genome and revealed the genome evolution and polyploidization events of H. cordata .

Molecular Characterization of Odorant-Binding Protein Genes Associated with Host-Seeking Behavior in Oides leucomelaena

The identification of odorant-binding proteins (OBPs) involved in host location by Oides leucomelaena (O. leucomelaena Weise, 1922, Coleoptera, Galerucinae) is significant for its biological control. Tools in the NCBI database were used to compare and analyze the transcriptome sequences of O. leucomelaena with OBP and other chemosensory-related proteins of other Coleoptera insects. Subsequently, MEGA7 was utilized for OBP sequence alignment and the construction of a phylogenetic tree, combined with expression profiling to screen for candidate antennae-specific OBPs. In addition, fumigation experiments with star anise volatiles were conducted to assess the antennae specificity of the candidate OBPs. Finally, molecular docking was employed to speculate on the binding potential of antennae-specific OBPs with star anise volatiles. The study identified 42 candidate OBPs, 8 chemosensory proteins and 27 receptors. OleuOBP3, OleuOBP5, and OleuOBP6 were identified as classic OBP family members specific to the antennae, which was confirmed by volatile fumigation experiments. Molecular docking ultimately clarified that OleuOBP3, OleuOBP5, and OleuOBP6 all exhibit a high affinity for β-caryophyllene among the star anise volatiles. We successfully obtained three antennae-specific OBPs from O. leucomelaena and determined their high-affinity volatiles, providing a theoretical basis for the development of attractants in subsequent stages.

Telomere-to-telomere genome assembly of Eleocharis dulcis and expression profiles during corm development

Eleocharis dulcis (Burm. f.) Trin. ex Hensch., commonly known as Chinese water chestnut, is a traditional aquatic vegetable in China, and now is widely cultivated throughout the world because of its high nutritional value and unique tastes. Here, we report the assembly of a 493.24 Mb telomere-to-telomere (T2T) genome of E. dulcis accomplished by integrating ONT ultra-long reads, PacBio long reads and Hi-C data. The reference genome was anchored onto 111 gap-free chromosomes, containing 48.31% repeat elements and 33,493 predicted protein-coding genes. Whole genome duplication (WGD) and inter-genomic synteny analyses indicated that chromosome breakage and genome duplication in E. dulcis possibly occurred multiple times during genome evolution after its divergence from a common ancestor with Rhynchospora breviuscula at ca. 35.6 Mya. A comparative time-course transcriptome analysis of corm development revealed different patterns of gene expression between cultivated and wild accessions with the highest number of differentially expressed genes (DEGs, 15,870) at the middle swelling stage and some of the DEGs were significantly enriched for starch metabolic process.

Community Resource: Large-Scale Proteogenomics to Refine Wheat Genome Annotations

Triticum aestivum is an important crop whose reference genome (International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.1) offers a valuable resource for understanding wheat genetic structure, improving agronomic traits, and developing new cultivars. A key aspect of gene model annotation is protein-level evidence of gene expression obtained from proteomics studies, followed up by proteogenomics to physically map proteins to the genome. In this research, we have retrieved the largest recent wheat proteomics datasets publicly available and applied the Basic Local Alignment Search Tool (tBLASTn) algorithm to map the 861,759 identified unique peptides against IWGSC RefSeq v2.1. Of the 92,719 hits, 83,015 unique peptides aligned along 33,612 High Confidence (HC) genes, thus validating 31.4% of all wheat HC gene models. Furthermore, 6685 unique peptides were mapped against 3702 Low Confidence (LC) gene models, and we argue that these gene models should be considered for HC status. The remaining 2934 orphan peptides can be used for novel gene discovery, as exemplified here on chromosome 4D. We demonstrated that tBLASTn could not map peptides exhibiting mid-sequence frame shift. We supply all our proteogenomics results, Galaxy workflow and Python code, as well as Browser Extensible Data (BED) files as a resource for the wheat community via the Apollo Jbrowse, and GitHub repositories. Our workflow could be applied to other proteomics datasets to expand this resource with proteins and peptides from biotically and abiotically stressed samples. This would help tease out wheat gene expression under various environmental conditions, both spatially and temporally.

De Novo Emerged Gene Search in Eukaryotes with DENSE

The discovery of de novo emerged genes, originating from previously noncoding DNA regions, challenges traditional views of species evolution. Indeed, the hypothesis of neutrally evolving sequences giving rise to functional proteins is highly unlikely. This conundrum has sparked numerous studies to quantify and characterize these genes, aiming to understand their functional roles and contributions to genome evolution. Yet, no fully automated pipeline for their identification is available. Therefore, we introduce DENSE (DE Novo emerged gene SEarch), an automated Nextflow pipeline based on two distinct steps: detection of taxonomically restricted genes (TRGs) through phylostratigraphy, and filtering of TRGs for de novo emerged genes via genome comparisons and synteny search. DENSE is available as a user-friendly command-line tool, while the second step is accessible through a web server upon providing a list of TRGs. Highly flexible, DENSE provides various strategy and parameter combinations, enabling users to adapt to specific configurations or define their own strategy through a rational framework, facilitating protocol communication, and study interoperability. We apply DENSE to seven model organisms, exploring the impact of its strategies and parameters on de novo gene predictions. This thorough analysis across species with different evolutionary rates reveals useful metrics for users to define input datasets, identify favorable/unfavorable conditions for de novo gene detection, and control potential biases in genome annotations. Additionally, predictions made for the seven model organisms are compiled into a requestable database, which we hope will serve as a reference for de novo emerged gene lists generated with specific criteria combinations.

Geographic EBV variants confound disease-specific variant interpretation and predict variable immune therapy responses

ABSTRACT

Epstein-Barr virus (EBV) is a potent carcinogen linked to hematologic and solid malignancies and causes significant global morbidity and mortality. Therapy using allogeneic EBV-specific lymphocytes shows promise in certain populations, but the impact of EBV genome variation on these strategies remains unexplored. To address this, we sequenced 217 EBV genomes, including hematologic malignancies from Guatemala, Peru, Malawi, and Taiwan, and analyzed them alongside 1307 publicly available EBV genomes from cancer, nonmalignant diseases, and healthy individuals across Africa, Asia, Europe, North America, and South America. These included, to our knowledge, the first natural killer (NK)/T-cell lymphoma (NKTCL) EBV genomes reported outside of East Asia. Our findings indicate that previously proposed EBV genome variants specific to certain cancer types are more closely tied to geographic origin than to cancer histology. This included variants previously reported to be specific to NKTCL but were prevalent in EBV genomes from other cancer types and healthy individuals in East Asia. After controlling for geographic region, we did identify multiple NKTCL-specific variants associated with a 7.8-fold to 21.9-fold increased risk. We also observed frequent variations in EBV genomes that affected peptide sequences previously reported to bind common major histocompatibility complex alleles. Finally, we found several nonsynonymous variants spanning the coding sequences of current vaccine targets BALF4, BKRF2, BLLF1, BXLF2, BZLF1, and BZLF2. These results highlight the need to consider geographic variation in EBV genomes when devising strategies for exploiting adaptive immune responses against EBV-related cancers, ensuring greater global effectiveness and equity in prevention and treatment.

A chromosome-level genome assembly and evolutionary analysis of Coregonus ussuriensis Berg

Coregonus ussuriensis Berg, distributed widely in cold waters above 45° N latitude, is a savored freshwater whitefish that has been included in the list of endangered animals as a consequence of overfishing. Lack of genomic information seriously hampers evolutionary and genetic research on C. ussuriensis warranting the need to assemble a high-quality reference genome to promote its genetic breeding. We assembled and constructed a reference chromosome-level C. ussuriensis genome (sequence length, 2.51 Gb; contig N50 length, 4.27 Mb) using PacBio sequencing and Hi-C assembly technology, 3,109 contigs were assembled into scaffolds, resulting in a genome assembly with 40 chromosomes and a scaffold N50 length of 62.20 Mb. In addition, 43,320 protein-coding genes were annotated. The peak Ks position in the species comparison reflects the whole-genome replication event of C. ussuriensis. This chromosome-level genome provides reference data for further studies on the molecular breeding of C. ussuriensis.

The complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909)

Holocephali is a subclass of chondrichthyans with ample geographic distribution in marine ecosystems. Holocephalan species are organized into three families: Callorhinchidae, Chimaeridae, and Rhinochimaeridae. Despite the critical ecological and evolutionary importance, genomic information from holocephalans is still scarce, particularly from rhinochimaerids. The present study provides the first complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909). The whole mitogenome was sequenced from an R. atlantica specimen, collected on the Porcupine Bank (NE Atlantic), by Illumina high-throughput sequencing. The R. atlantica mitogenome has 17,852 nucleotides with 13 protein-coding genes, 22 transfer RNA, and two ribosomal RNA genes. Nine of these genes are in the complementary strand. This mitogenome has a GC content of 41.5% and an AT content of 58.5%. The phylogenetic reconstruction provided here, using all the available complete and partial Holocephali mitogenomes, places R. atlantica in the Rhinochimaeridae family, as expected. This genomic resource will be useful in the genomic characterization of this species.

A conserved phenylalanine motif among teleost fish provides insight for improving electromagnetic perception

Magnetoreceptive biology as a field remains relatively obscure; compared with the breadth of species believed to sense magnetic fields, it remains under-studied. Here, we present grounds for the expansion of magnetoreception studies among teleosts. We begin with the electromagnetic perceptive gene (EPG) from Kryptopterus vitreolus and expand to identify 72 teleosts with homologous proteins containing a conserved three-phenylalanine (3F) motif. Phylogenetic analysis provides insight as to how EPG may have evolved over time and indicates that certain clades may have experienced a loss of function driven by different fitness pressures. One potential factor is water type with freshwater fish significantly more likely to possess the functional motif version (FFF), and saltwater fish to have the non-functional variant (FXF). It was also revealed that when the 3F motif from the homologue of Brachyhypopomus gauderio (B.g.) is inserted into EPG-EPG(B.g.)-the response (as indicated by increased intracellular calcium) is faster. This indicates that EPG has the potential to be engineered to improve upon its response and increase its utility to be used as a controller for specific outcomes.

Restoration of the Functional nif Gene Cluster by Complex Recombination Events during Heterocyst Development in the Nitrogen-Fixing Cyanobacterium Calothrix sp. NIES-4101

In the genome of the heterocystous cyanobacterium Calothrix sp. NIES-4101 (NIES-4101), the four genes essential for nitrogen fixation (nifB, nifH, nifD and nifK) are highly fragmented into 13 parts in a 350-kb chromosomal region, and four of these parts are encoded in the reverse strand. Such a complex fragmentation feature makes it difficult to restore the intact nifBHDK genes by the excision mechanism found in the nifD gene of the Anabaena sp. PCC 7120 heterocyst. To examine the nitrogen-fixing ability of NIES-4101, we confirmed that NIES-4101 grew well on a combined nitrogen-free medium and showed high nitrogenase activity, which strongly suggested that the complete nifBHDK genes are restored by a complex recombination process in heterocysts. Next, we resequenced the genome prepared from cells grown under nitrogen-fixing conditions. Two contigs covering the complete nifHDK and nifB genes were found by de novo assembly of the sequencing reads. In addition, the DNA fragments covering the nifBHDK operon were successfully amplified by PCR. We propose that the process of nifBHDK restoration occurs as follows. First, the nifD-nifK genes are restored by four excision events. Then, the complete nifH and nifB genes are restored by two excision events followed by two successive inversion events between the inverted repeat sequences and one excision event, forming the functional nif gene cluster, nifB-fdxN-nifS-nifU-nifH-nifD-nifK. All genes coding recombinases responsible for these nine recombination events are located close to the terminal repeat sequences. The restoration of the nifBHDK genes in NIES-4101 is the most complex genome reorganization reported in heterocystous cyanobacteria.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species' adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.

Chromosome-Level Genome Assembly of Apoderus dimidiatus Voss (Coleoptera: Attelabidae): Insights into Evolution and Behavior

Attelabidae insects have attracted much attention due to their unique leaf rolling behavior before oviposition. However, the lack of genomic data makes it difficult to understand the molecular mechanism behind their behavior and their evolutionary relationship with other species. To address this gap, we utilized Illumina and Nanopore sequencing platforms along with Hi-C technology to establish a highly accurate whole genome of A. dimidiatus at the chromosome level. The resulting genome size was determined to be 619.26 Mb, with a contig N50 of 50.89 Mb and GC content of 33.89%. Moreover, a total of 12,572 genes were identified, with 82.59% being functionally annotated, and 64.78% designated as repeat sequences. Our subsequent phylogenetic tree analysis revealed that Attelabidae's divergence from Curculionidae occurred approximately 161.52 million years ago. Furthermore, the genome of A. dimidiatus contained 334 expanded gene families and 1718 contracted gene families. In addition, using Phylogenetic Analysis by Maximum Likelihood (PAML), we identified 106 rapidly evolved genes exhibiting significant signals and 540 positively selected genes. Our research endeavors to serve as an invaluable genomic data resource for the study of Attelabidae, offering fresh perspectives for the exploration of its leaf rolling behavior.

Identification and characterization of STAT family in silver pomfret (Pampus argenteus) involved in different exogenous stresses

Members of the Signal Transducer and Activator of Transcription (STAT) family function pivotally as transcriptional activators integral to the modulation of inflammatory responses. The aquaculture of silver pomfret is frequently compromised by the imposition of exogenous stressors, which include thermal fluctuations, notably low-temperatures, diminished oxygen levels, and the onslaught of bacterial pathogens. Notwithstanding the critical impact of these stressors, the scientific literature presents a notable gap in our understanding of the STAT pathway's role in the silver pomfret's adaptive response mechanisms. To address this lacuna, we identified stat genes in the silver pomfret-denominated as Pastat1, Pastat2, Pastat3, Pastat4, and Pastat5-through a thorough and systematic bioinformatics analysis. Further scrutiny of the gene configurations and constituent motifs has elucidated that STAT proteins possess analogous structural frameworks and exhibit significant evolutionary preservation. Subsequently, the expression patterns of five stat genes were verified by RT-qPCR in twelve different tissues and four growth periods in healthy fish, showing that the expression of Pastat genes was temporally and spatially specific, with most of the stat genes expressed at higher levels in the spleen, following muscle, gill, and liver. Transcriptomic analysis of exposure to exogenous stressors, specifically formaldehyde and low-temperature conditions, elucidated that Pastat1 and Pastat2 genes exhibited a heightened sensitivity to these environmental challenges. RT-qPCR assays demonstrated a marked alteration in the expression profiles of jak1 and Pastat gene suites in PaS upon prolonged bacterial infection subsequent to these exogenous insults. Moreover, the gene expression of the downstream effectors involved in innate immunity and apoptosis displayed marked deviations. This study additionally elucidated the Pastat gene family's role in modulating the innate immune response and apoptotic regulation within the silver pomfret during exogenous stressors and subsequent pathogenic incursions.

Metagenomic exploration of cold-active enzymes for detergent applications: Characterization of a novel, cold-active and alkali-stable GH8 endoglucanase from ikaite columns in SW Greenland

Microbial communities from extreme environments are largely understudied, but are essential as producers of metabolites, including enzymes, for industrial processes. As cultivation of most microorganisms remains a challenge, culture-independent approaches for enzyme discovery in the form of metagenomics to analyse the genetic potential of a community are rapidly becoming the way forward. This study focused on analysing a metagenome from the cold and alkaline ikaite columns in Greenland, identifying 282 open reading frames (ORFs) that encoded putative carbohydrate-modifying enzymes with potential applications in, for example detergents and other processes where activity at low temperature and high pH is desired. Seventeen selected ORFs, representing eight enzyme families were synthesized and expressed in two host organisms, Escherichia coli and Aliivibrio wodanis. Aliivibrio wodanis demonstrated expression of a more diverse range of enzyme classes compared to E. coli, emphasizing the importance of alternative expression systems for enzymes from extremophilic microorganisms. To demonstrate the validity of the screening strategy, we chose a recombinantly expressed cellulolytic enzyme from the metagenome for further characterization. The enzyme, Cel240, exhibited close to 40% of its relative activity at low temperatures (4°C) and demonstrated endoglucanase characteristics, with a preference for cellulose substrates. Despite low sequence similarity with known enzymes, computational analysis and structural modelling confirmed its cellulase-family affiliation. Cel240 displayed activity at low temperatures and good stability at 25°C, activity at alkaline pH and increased activity in the presence of CaCl2, making it a promising candidate for detergent and washing industry applications.

Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Metabolomic and transcriptomic analyses of yellow-flowered crocuses to infer alternative sources of saffron metabolites

BACKGROUND

The increasing demand for saffron metabolites in various commercial industries, including medicine, food, cosmetics, and dyeing, is driven by the discovery of their diverse applications. Saffron, derived from Crocus sativus stigmas, is the most expensive spice, and there is a need to explore additional sources to meet global consumption demands. In this study, we focused on yellow-flowering crocuses and examined their tepals to identify saffron-like compounds.

RESULTS

Through metabolomic and transcriptomic approaches, our investigation provides valuable insights into the biosynthesis of compounds in yellow-tepal crocuses that are similar to those found in saffron. The results of our study support the potential use of yellow-tepal crocuses as a source of various crocins (crocetin glycosylated derivatives) and flavonoids.

CONCLUSIONS

Our findings suggest that yellow-tepal crocuses have the potential to serve as a viable excessive source of some saffron metabolites. The identification of crocins and flavonoids in these crocuses highlights their suitability for meeting the demands of various industries that utilize saffron compounds. Further exploration and utilization of yellow-tepal crocuses could contribute to addressing the growing global demand for saffron-related products.

Development and evaluation of genomics informed real-time PCR assays for the detection and strain typing of Mycobacterium avium subsp. paratuberculosis

AIMS

This study aimed to identify specific genomic targets for the detection and strain typing of Map and analyse their sensitivity and specificity, and detect Map directly from faeces.

METHODS AND RESULTS

A comparative genomics approach was used to identify specific genomic targets for the detection and strain typing of Map. A Map specific qPCR using the primer pair 7132 that targets a DNA segregation ATPase protein was able to detect all strains of Map and is more sensitive than the current Johne's disease PCR assays with a sensitivity of 0.0002 fg µl-1. A strain specific qPCR using the Atsa primer pair that targets the arylsulfase gene was able to differentiate between Type S and Type C strains of Map and was more sensitive than the IS1311 PCR and REA with a sensitivity of 40 fg µl-1 and was specific for Type S Map. Both assays successfully detected Map directly from faeces.

CONCLUSION

This study developed and validated two genomics informed qPCR assays, 7132B Map and Atsa Type S and found both assays to be highly specific and sensitive for the detection of Map from culture and directly from faeces. This is the first time that a probe-based qPCR has been designed and developed for Map strain typing, which will greatly improve the response time during outbreak investigations.

The first chromosome-level genome of the stag beetle Dorcus hopei Saunders, 1854 (Coleoptera: Lucanidae)

Stag beetles (Coleoptera: Lucanidae) represent a significant saproxylic assemblage in forest ecosystems and are noted for their enlarged mandibles and male polymorphism. Despite their relevance as ideal models for the study of exaggerated mandibles that aid in attracting mates, the regulatory mechanisms associated with these traits remain understudied, and restricted by the lack of high-quality reference genomes for stag beetles. To address this limitation, we successfully assembled the first chromosome-level genome of a representative species Dorcus hopei. The genome was 496.58 Mb in length, with a scaffold N50 size of 54.61 Mb, BUSCO values of 99.8%, and 96.8% of scaffolds anchored to nine pairs of chromosomes. We identified 285.27 Mb (57.45%) of repeat sequences and annotated 11,231 protein-coding genes. This genome will be a valuable resource for further understanding the evolution and ecology of stag beetles, and provides a basis for studying the mechanisms of exaggerated mandibles through comparative analysis.

Genome assembly of Melilotus officinalis provides a new reference genome for functional genomics

BACKGROUND

Sweet yellow clover (Melilotus officinalis) is a diploid plant (2n = 16) that is native to Europe. It is an excellent legume forage. It can both fix nitrogen and serve as a medicine. A genome assembly of Melilotus officinalis that was collected from Best corporation in Beijing is available based on Nanopore sequencing. The genome of Melilotus officinalis was sequenced, assembled, and annotated.

RESULTS

The latest PacBio third generation HiFi assembly and sequencing strategies were used to produce a Melilotus officinalis genome assembly size of 1,066 Mbp, contig N50 = 5 Mbp, scaffold N50 = 130 Mbp, and complete benchmarking universal single-copy orthologs (BUSCOs) = 96.4%. This annotation produced 47,873 high-confidence gene models, which will substantially aid in our research on molecular breeding. A collinear analysis showed that Melilotus officinalis and Medicago truncatula shared conserved synteny. The expansion and contraction of gene families showed that Melilotus officinalis expanded by 565 gene families and shrank by 56 gene families. The contacted gene families were associated with response to stimulus, nucleotide binding, and small molecule binding. Thus, it is related to a family of genes associated with peptidase activity, which could lead to better stress tolerance in plants.

CONCLUSIONS

In this study, the latest PacBio technology was used to assemble and sequence the genome of the Melilotus officinalis and annotate its protein-coding genes. These results will expand the genomic resources available for Melilotus officinalis and should assist in subsequent research on sweet yellow clover plants.

A conserved phenylalanine motif among Teleost fish provides insight for improving electromagnetic perception

Magnetoreceptive biology as a field remains relatively obscure; compared to the breadth of species believed to sense magnetic fields, it remains under-studied. Here, we present grounds for the expansion of magnetoreception studies among Teleosts. We begin with the electromagnetic perceptive gene (EPG) from Kryptopterus vitreolus and expand to identify 72 Teleosts with homologous proteins containing a conserved three-phenylalanine (3F) motif. Phylogenetic analysis provides insight as to how EPG may have evolved over time, and indicates that certain clades may have experienced a loss of function driven by different fitness pressures. One potential factor is water type with freshwater fish significantly more likely to possess the functional motif version (FFF), and saltwater fish to have the non-functional variant (FXF). It was also revealed that when the 3F motif from the homolog of Brachyhypopomus gauderio (B.g.) is inserted into EPG - EPG(B.g.) - the response (as indicated by increased intracellular calcium) is faster. This indicates that EPG has the potential to be engineered to improve upon its response and increase its utility to be used as a controller for specific outcomes.

An inside out journey: biogenesis, ultrastructure and proteomic characterisation of the ectoparasitic flatworm Sparicotyle chrysophrii extracellular vesicles

BACKGROUND

Helminth extracellular vesicles (EVs) are known to have a three-way communication function among parasitic helminths, their host and the host-associated microbiota. They are considered biological containers that may carry virulence factors, being therefore appealing as therapeutic and prophylactic target candidates. This study aims to describe and characterise EVs secreted by Sparicotyle chrysophrii (Polyopisthocotyla: Microcotylidae), a blood-feeding gill parasite of gilthead seabream (Sparus aurata), causing significant economic losses in Mediterranean aquaculture.

METHODS

To identify proteins involved in extracellular vesicle biogenesis, genomic datasets from S. chrysophrii were mined in silico using known protein sequences from Clonorchis spp., Echinococcus spp., Fasciola spp., Fasciolopsis spp., Opisthorchis spp., Paragonimus spp. and Schistosoma spp. The location and ultrastructure of EVs were visualised by transmission electron microscopy after fixing adult S. chrysophrii specimens by high-pressure freezing and freeze substitution. EVs were isolated and purified from adult S. chrysophrii (n = 200) using a newly developed ultracentrifugation-size-exclusion chromatography protocol for Polyopisthocotyla, and EVs were characterised via nanoparticle tracking analysis and tandem mass spectrometry.

RESULTS

Fifty-nine proteins involved in EV biogenesis were identified in S. chrysophrii, and EVs compatible with ectosomes were observed in the syncytial layer of the haptoral region lining the clamps. The isolated and purified nanoparticles had a mean size of 251.8 nm and yielded 1.71 × 108 particles · mL-1. The protein composition analysis identified proteins related to peptide hydrolases, GTPases, EF-hand domain proteins, aerobic energy metabolism, anticoagulant/lipid-binding, haem detoxification, iron transport, EV biogenesis-related, vesicle-trafficking and other cytoskeletal-related proteins. Several identified proteins, such as leucyl and alanyl aminopeptidases, calpain, ferritin, dynein light chain, 14-3-3, heat shock protein 70, annexin, tubulin, glutathione S-transferase, superoxide dismutase, enolase and fructose-bisphosphate aldolase, have already been proposed as target candidates for therapeutic or prophylactic purposes.

CONCLUSIONS

We have unambiguously demonstrated for the first time to our knowledge the secretion of EVs by an ectoparasitic flatworm, inferring their biogenesis machinery at a genomic and transcriptomic level, and by identifying their location and protein composition. The identification of multiple therapeutic targets among EVs' protein repertoire provides opportunities for target-based drug discovery and vaccine development for the first time in Polyopisthocotyla (sensu Monogenea), and in a fish-ectoparasite model.

RrA, an enzyme from Rhodospirillum rubrum, is a prototype of a new family of short-chain L-asparaginases

L-Asparaginases (ASNases) catalyze the hydrolysis of L-Asn to L-Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300-400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum, reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N-terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L-Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short-chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild-type RrA, both with and without bound L-Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand-alone subfamily of ASNases that can catalyze the hydrolysis of L-Asn to L-Asp despite the lack of the C-terminal domain that is present in all ASNases described structurally to date.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Regional centromere configuration in the fungal pathogens of the Pneumocystis genus

Centromeres are constricted chromosomal regions that are essential for cell division. In eukaryotes, centromeres display a remarkable architectural and genetic diversity. The basis of centromere-accelerated evolution remains elusive. Here, we focused on Pneumocystis species, a group of mammalian-specific fungal pathogens that form a sister taxon with that of the Schizosaccharomyces pombe, an important genetic model for centromere biology research. Methods allowing reliable continuous culture of Pneumocystis species do not currently exist, precluding genetic manipulation. CENP-A, a variant of histone H3, is the epigenetic marker that defines centromeres in most eukaryotes. Using heterologous complementation, we show that the Pneumocystis CENP-A ortholog is functionally equivalent to CENP-ACnp1 of S. pombe. Using organisms from a short-term in vitro culture or infected animal models and chromatin immunoprecipitation (ChIP)-Seq, we identified CENP-A bound regions in two Pneumocystis species that diverged ~35 million years ago. Each species has a unique short regional centromere (<10 kb) flanked by heterochromatin in 16-17 monocentric chromosomes. They span active genes and lack conserved DNA sequence motifs and repeats. These features suggest an epigenetic specification of centromere function. Analysis of centromeric DNA across multiple Pneumocystis species suggests a vertical transmission at least 100 million years ago. The common ancestry of Pneumocystis and S. pombe centromeres is untraceable at the DNA level, but the overall architectural similarity could be the result of functional constraint for successful chromosomal segregation.IMPORTANCEPneumocystis species offer a suitable genetic system to study centromere evolution in pathogens because of their phylogenetic proximity with the non-pathogenic yeast S. pombe, a popular model for cell biology. We used this system to explore how centromeres have evolved after the divergence of the two clades ~ 460 million years ago. To address this question, we established a protocol combining short-term culture and ChIP-Seq to characterize centromeres in multiple Pneumocystis species. We show that Pneumocystis have short epigenetic centromeres that function differently from those in S. pombe.

Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum

BACKGROUND

Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China, the Philippines and Indonesia. Oncomelania hupensis (Gastropoda: Pomatiopsidae) is the unique intermediate host of S. japonicum. A complete genome sequence of O. hupensis will enable the fundamental understanding of snail biology as well as its co-evolution with the S. japonicum parasite. Assembling a high-quality reference genome of O. hupehensis will provide data for further research on the snail biology and controlling the spread of S. japonicum.

METHODS

The draft genome was de novo assembly using the long-read sequencing technology (PacBio Sequel II) and corrected with Illumina sequencing data. Then, using Hi-C sequencing data, the genome was assembled at the chromosomal level. CAFE was used to do analysis of contraction and expansion of the gene family and CodeML module in PAML was used for positive selection analysis in protein coding sequences.

RESULTS

A total length of 1.46 Gb high-quality O. hupensis genome with 17 unique full-length chromosomes (2n = 34) of the individual including a contig N50 of 1.35 Mb and a scaffold N50 of 75.08 Mb. Additionally, 95.03% of these contig sequences were anchored in 17 chromosomes. After scanning the assembled genome, a total of 30,604 protein-coding genes were predicted. Among them, 86.67% were functionally annotated. Further phylogenetic analysis revealed that O. hupensis was separated from a common ancestor of Pomacea canaliculata and Bellamya purificata approximately 170 million years ago. Comparing the genome of O. hupensis with its most recent common ancestor, it showed 266 significantly expanded and 58 significantly contracted gene families (P < 0.05). Functional enrichment of the expanded gene families indicated that they were mainly involved with intracellular, DNA-mediated transposition, DNA integration and transposase activity.

CONCLUSIONS

Integrated use of multiple sequencing technologies, we have successfully constructed the genome at the chromosomal-level of O. hupensis. These data will not only provide the compressive genomic information, but also benefit future work on population genetics of this snail as well as evolutional studies between S. japonicum and the snail host.

Hfe Permease and Haemophilus influenzae Manganese Homeostasis

Haemophilus influenzae is a commensal of the human upper respiratory tract that can infect diverse host niches due, at least in part, to its ability to withstand both endogenous and host-mediated oxidative stresses. Here, we show that hfeA, a gene previously linked to iron import, is essential for H. influenzae manganese recruitment via the HfeBCD transporter. Structural analyses show that metal binding in HfeA uses a unique mechanism that involves substantial rotation of the C-terminal lobe of the protein. Disruption of hfeA reduced H. influenzae manganese acquisition and was associated with decreased growth under aerobic conditions, impaired manganese-superoxide dismutase activity, reduced survival in macrophages, and changes in biofilm production in the presence of superoxide. Collectively, this work shows that HfeA contributes to H. influenzae manganese acquisition and virulence attributes. High conservation of the hfeABCD permease in Haemophilus species suggests that it may serve similar roles in other pathogenic Pasteurellaceae.

Chromosome-level genome assembly and annotation of eel goby (Odontamblyopus rebecca)

The eel gobies fascinate researchers with many important features, including its unique body structure, benthic lifestyle, and degenerated eyes. However, genome assembly and exploration of the unique genomic composition of the eel gobies are still in their infancy. This has severely limited research progress on gobies. In this study, multi-platform sequencing data were generated and used to assemble and annotate the genome of O. rebecca at the chromosome-level. The assembled genome size of O. rebecca is 918.57 Mbp, which is similar to the estimated genome size (903.03 Mbp) using 17-mer. The scaffold N50 is 41.67 Mbp, and 23 chromosomes were assembled using Hi-C technology with a mounting rate of 99.96%. Genome annotation indicates that 53.29% of the genome is repetitive sequences, and 22,999 protein-coding genes are predicted, of which 21,855 have functional annotations. The chromosome-level genome of O. rebecca will not only provide important genomic resources for comparative genomic studies of gobies, but also expand our knowledge of the genetic origin of their unique features fascinating researchers for decades.

Expression patterns of FGF and BMP pathway genes in the tardigrade Hypsibius exemplaris

A small number of conserved signaling pathways regulate development of most animals, yet we do not know where these pathways are deployed in most embryos. This includes tardigrades, a phylum with a unique body shape. We examined expression patterns of components of the BMP and FGF signaling pathways during embryonic segmentation and mesoderm development of the tardigrade Hypsibius exemplaris. Among the patterns examined, we found that an FGF ligand gene is expressed in ectodermal segment posteriors and an FGF receptor gene is expressed in underlying endomesodermal pouches, suggesting possible FGF signaling between these developing germ layers. We found that a BMP ligand gene is expressed in lateral ectoderm and dorsolateral bands along segment posteriors, while the BMP antagonist Sog gene is expressed in lateral ectoderm and also in a subset of endomesodermal cells, suggesting a possible role of BMP signaling in dorsal-ventral patterning of lateral ectoderm. In combination with known roles of these pathways during development of common model systems, we developed hypotheses for how the BMP and FGF pathways might regulate embryo segmentation and mesoderm formation of the tardigrade H. exemplaris. These results identify the expression patterns of genes from two conserved signaling pathways for the first time in the tardigrade phylum.

Whole genome sequencing of a novel sea anemone (Actinostola sp.) from a deep-sea hydrothermal vent

Deep-sea hydrothermal vents are usually considered as extreme environments with high pressure, high temperature, scarce food, and chemical toxicity, while many local inhabitants have evolved special adaptive mechanisms for residence in this representative ecosystem. In this study, we constructed a high-quality genome assembly for a novel deep-sea anemone species (Actinostola sp.) that was resident at a depth of 2,971 m in an Edmond vent along the central Indian Ocean ridge, with a total size of 424.3 Mb and a scaffold N50 of 383 kb. The assembled genome contained 265 Mb of repetitive sequences and 20,812 protein-coding genes. Taken together, our reference genome provides a valuable genetic resource for exploring the evolution and adaptive clues of this deep-sea anemone.

Establishing an itaconic acid production process with Ustilago species on the low-cost substrate starch

Ustilago maydis and Ustilago cynodontis are natural producers of a broad range of valuable molecules including itaconate, malate, glycolipids, and triacylglycerols. Both Ustilago species are insensitive toward medium impurities, and have previously been engineered for efficient itaconate production and stabilized yeast-like growth. Due to these features, these strains were already successfully used for the production of itaconate from different alternative feedstocks such as molasses, thick juice, and crude glycerol. Here, we analyzed the amylolytic capabilities of Ustilago species for metabolization of starch, a highly abundant and low-cost polymeric carbohydrate widely utilized as a substrate in several biotechnological processes. Ustilago cynodontis was found to utilize gelatinized potato starch for both growth and itaconate production, confirming the presence of extracellular amylolytic enzymes in Ustilago species. Starch was rapidly degraded by U. cynodontis, even though no α-amylase was detected. Further experiments indicate that starch hydrolysis is caused by the synergistic action of glucoamylase and α-glucosidase enzymes. The enzymes showed a maximum activity of around 0.5 U ml-1 at the fifth day after inoculation, and also released glucose from additional substrates, highlighting potential broader applications. In contrast to U. cynodontis, U. maydis showed no growth on starch accompanied with no detectable amylolytic activity.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

Host genetic regulation of human gut microbial structural variation

Although the impact of host genetics on gut microbial diversity and the abundance of specific taxa is well established1-6, little is known about how host genetics regulates the genetic diversity of gut microorganisms. Here we conducted a meta-analysis of associations between human genetic variation and gut microbial structural variation in 9,015 individuals from four Dutch cohorts. Strikingly, the presence rate of a structural variation segment in Faecalibacterium prausnitzii that harbours an N-acetylgalactosamine (GalNAc) utilization gene cluster is higher in individuals who secrete the type A oligosaccharide antigen terminating in GalNAc, a feature that is jointly determined by human ABO and FUT2 genotypes, and we could replicate this association in a Tanzanian cohort. In vitro experiments demonstrated that GalNAc can be used as the sole carbohydrate source for F. prausnitzii strains that carry the GalNAc-metabolizing pathway. Further in silico and in vitro studies demonstrated that other ABO-associated species can also utilize GalNAc, particularly Collinsella aerofaciens. The GalNAc utilization genes are also associated with the host's cardiometabolic health, particularly in individuals with mucosal A-antigen. Together, the findings of our study demonstrate that genetic associations across the human genome and bacterial metagenome can provide functional insights into the reciprocal host-microbiome relationship.

Acinetobacter baumannii coordinates central metabolism, plasmid dissemination, and virulence by sensing nutrient availability

IMPORTANCE

Plasmid conjugation is known to be an energy-expensive process, but our understanding of the molecular linkage between conjugation and metabolism is limited. Our finding reveals that Acinetobacter baumannii utilizes a two-component system to co-regulate metabolism, plasmid transfer, and virulence by sensing reaction intermediates of key metabolic pathways, which suggests that nutrient availability dictates not only bacterial proliferation but also horizontal gene transfer. The identification of Dot/Icm-like proteins as components of a conjugation system involved in the dissemination of antibiotic-resistance genes by A. baumannii has provided important targets for the development of agents capable of inhibiting virulence and the spread of anti-microbial-resistance genes in bacterial communities.

Whole genome assembly and annotation of the endangered Caribbean coral Acropora cervicornis

Coral species in the genus Acropora are key ecological components of coral reefs worldwide and represent the most diverse genus of scleractinian corals. While key species of Indo-Pacific Acropora have annotated genomes, no annotated genome has been published for either of the two species of Caribbean Acropora. Here we present the first fully annotated genome of the endangered Caribbean staghorn coral, Acropora cervicornis. We assembled and annotated this genome using high-fidelity nanopore long-read sequencing with gene annotations validated with mRNA sequencing. The assembled genome size is 318 Mb, with 28,059 validated genes. Comparative genomic analyses with other Acropora revealed unique features in A. cervicornis, including contractions in immune pathways and expansions in signaling pathways. Phylogenetic analysis confirms previous findings showing that A. cervicornis diverged from Indo-Pacific relatives around 41 million years ago, with the closure of the western Tethys Sea, prior to the primary radiation of Indo-Pacific Acropora. This new A. cervicornis genome enriches our understanding of the speciose Acropora and addresses evolutionary inquiries concerning speciation and hybridization in this diverse clade.

Bioprospecting Microalgae from Sewage Water: Assessment of Biochemicals for Biomass Utilization

Microalgal species from sewage treatment plant were identified by 18S rRNA sequencing and were explored for total lipids, carbohydrate, and protein contents, to serve as a potential candidate for biorefinery. Seven unicellular microalgae were identified as Chlorella sorokiniana, Dictyosphaerium sp., Graesiella emersonii belonging to Chlorellaceae and Scenedesmus sp., Desmodesmus sp., Tetranephris brasiliensis, and Coelastrella sp. belonging to Scenedesmaceae family. Biochemical assessment of all isolates revealed total lipid content from 17.49 ± 1.41 to 47.35 ± 0.61% w/w, total carbohydrate content from 12.82 ± 0.19 to 64.29 ± 0.63% w/w, and total protein content from 8.55 ± 0.19 to 16.65 ± 0.20% w/w. FAME analysis of extracted lipid was found to be rich in Hexadecane (C16:0), Tetradecane (C17:0), Octadecane (C18:0), Eicosane (C20:0), Tetracosane (C24:0), Pentacosane (C25:0) fatty acids, the presence of which makes excellent candidate for biodiesel. Being rich in lipid, microalgae Chlorella sorokiniana, Coelastrella sp., and Scenedesmus sp. have high potential for biofuels. Due to the presence of high protein content, Scenedesmus sp. and Chlorella sorokiniana can serve as food or feed supplement, whereas the high carbohydrate content of Dictyosphaerium sp., Coelastrella sp., and Scenedesmus sp. makes them an ideal candidate for fermentative production of alcohol and organic acids. Chlorella sp. and Scenedesmus sp., being dominant microalgae across all seasons, demonstrate remarkable resilience for their cultivation in sewage water and utilization of biomass in biorefineries.

Chromosome-level genome assembly of Przevalski's partridge (Alectoris magna)

Przevalski's partridge (Alectoris magna) is one of the birds in the genus Alectoris endemic to China. The distribution of A. magna was narrow, and it was only found in parts of the Qinghai, Gansu, and Ningxia provinces. A. magna was considered a monotypic species until it was distinguished into two subspecies. However, external morphological characteristics, rather than genetic differences or evolutionary relationships, are now commonly used as evidence of subspecies differentiation. In this study, a chromosome-level reference genome of A. magna has been constructed by combining Illumina, PacBio and Hi-C sequencing data. The 1135.01 Mb A. magna genome was ultimately assembled. The genome showed 96.9% completeness (BUSCO), with a contig N50 length of 23.34 Mb. The contigs were clustered and oriented on 20 chromosomes, covering approximately 99.96% of the genome assembly. Additionally, altogether 19,103 protein-coding genes were predicted, of which 95.10% were functionally annotated. This high-quality genome assembly could serve as a valuable genomic resource for future research on the functional genomics, genetic protection, and interspecific hybridization of A. magna.

De novo genome assembly of the invasive mosquito species Aedes japonicus and Aedes koreicus

BACKGROUND

Recently, two invasive Aedes mosquito species, Ae. japonicus and Ae. koreicus, are circulating in several European countries posing potential health risks to humans and animals. Vector control is the main option to prevent mosquito-borne diseases, and an accurate genome sequence of these mosquitoes is essential to better understand their biology and to develop effective control strategies.

METHODS

A de novo genome assembly of Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) has been produced based on a hybrid approach that combines Oxford Nanopore long-read and Illumina short-read data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed.

RESULTS

Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNA-seq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures, we also identified genes showing a differential temperature-dependent activation.

CONCLUSIONS

The assembly of Akor1 and Ajap1 genomes constitutes the first updated collective knowledge of the genomes of both mosquito species, providing the possibility of understanding key mechanisms of their biology such as the ability to adapt to harsh climates and to develop insecticide-resistance mechanisms.

Insights into RAG Evolution from the Identification of "Missing Link" Family A RAGL Transposons

A series of "molecular domestication" events are thought to have converted an invertebrate RAG-like (RAGL) transposase into the RAG1-RAG2 (RAG) recombinase, a critical enzyme for adaptive immunity in jawed vertebrates. The timing and order of these events are not well understood, in part because of a dearth of information regarding the invertebrate RAGL-A transposon family. In contrast to the abundant and divergent RAGL-B transposon family, RAGL-A most closely resembles RAG and is represented by a single orphan RAG1-like (RAG1L) gene in the genome of the hemichordate Ptychodera flava (PflRAG1L-A). Here, we provide evidence for the existence of complete RAGL-A transposons in the genomes of P. flava and several echinoderms. The predicted RAG1L-A and RAG2L-A proteins encoded by these transposons intermingle sequence features of jawed vertebrate RAG and RAGL-B transposases, leading to a prediction of DNA binding, catalytic, and transposition activities that are a hybrid of RAG and RAGL-B. Similarly, the terminal inverted repeats (TIRs) of the RAGL-A transposons combine features of both RAGL-B transposon TIRs and RAG recombination signal sequences. Unlike all previously described RAG2L proteins, RAG2L-A proteins contain an acidic hinge region, which we demonstrate is capable of efficiently inhibiting RAG-mediated transposition. Our findings provide evidence for a critical intermediate in RAG evolution and argue that certain adaptations thought to be specific to jawed vertebrates (e.g. the RAG2 acidic hinge) actually arose in invertebrates, thereby focusing attention on other adaptations as the pivotal steps in the completion of RAG domestication in jawed vertebrates.

Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization

Water caltrop (Trapa spp., Lythraceae) is a traditional but currently underutilized non-cereal crop. Here, we generated chromosome-level genome assemblies for the two diploid progenitors of allotetraploid Trapa. natans (4x, AABB), i.e., diploid T. natans (2x, AA) and Trapa incisa (2x, BB). In conjunction with four published (sub)genomes of Trapa, we used gene-based and graph-based pangenomic approaches and a pangenomic transposable element (TE) library to develop Trapa genomic resources. The pangenome displayed substantial gene-content variation with dispensable and private gene clusters occupying a large proportion (51.95%) of the total cluster sets in the six (sub)genomes. Genotyping of presence-absence variation (PAVs) identified 40 453 PAVs associated with 2570 genes specific to A- or B-lineages, of which 1428 were differentially expressed, and were enriched in organ development process, organic substance metabolic process and response to stimulus. Comparative genome analyses showed that the allotetraploid T. natans underwent asymmetric subgenome divergence, with the B-subgenome being more dominant than the A-subgenome. Multiple factors, including PAVs, asymmetrical amplification of TEs, homeologous exchanges (HEs), and homeolog expression divergence, together affected genome evolution after polyploidization. Overall, this study sheds lights on the genome architecture and evolution of Trapa, and facilitates its functional genomic studies and breeding program.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.