The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999

Microbiome-metabolomic insights into the systemic regulation in Fangxian Huangjiu fermentation

Metabolic forces drive microecological succession in Huangjiu fermentation. This study investigates the dynamic metabolic-microbial interplay during Fangxian Huangjiu fermentation. Temporal changes of metabolome and microbiome revealed a syntropic relationship that purified the microbial community with convergent metabolic patterns. With species turnover driving microbial community structure, early-stage microbiomes exhibited great functional diversity. Functions related to energy and molecular building blocks were enriched at the end of early stage, and contributed greatly to microbial adaptation, highlighting the importance of metabolic forces in shaping community structure. Proteobacteria were identified as key facilitators of diverse metabolic activities, and Enterobacter emerged as a fundamental microbial community particularly for materials transformation. Correlation analysis enriched amino acid metabolism pathways. Further, Pantoea ananatis and Wickerhamomyces anomalus were isolated to enhance sphingosine-1-phosphate, γ-aminobutyric acid, and creatine levels without altering physicochemical properties. The study offers insights into the regulation of Huangjiu fermentation, and suggested potential micobiome manipulation to optimize characteristics.

Construction of the super pan-genome for the genus Actinidia reveals structural variations linked to phenotypic diversity

Kiwifruits, belonging to the genus Actinidia, are acknowledged as one of the most successfully domesticated fruits in the twentieth century. Despite the rich wild resources and diverse phenotypes within this genus, insights into the genomic changes are still limited. Here, we conducted whole-genome sequencing on seven representative materials from highly diversified sections of Actinidia, leading to the assembly and annotation of 14 haplotype genomes with sizes spanning from 602.0 to 699.6 Mb. By compiling these haplotype genomes, we constructed a super pan-genome for the genus. We identified numerous structural variations (SVs, including variations in gene copy number) and highly diverged regions in these genomes. Notably, significant SV variability was observed within the intronic regions of the MED25 and TTG1 genes across different materials, suggesting their potential roles in influencing fruit size and trichome formation. Intriguingly, our findings indicated a high genetic divergence between two haplotype genomes, with one individual, tentatively named Actinidia × leiocacarpae, from sect. Leiocacarpae. This likely hybrid with a heterozygous genome exhibited notable genetic adaptations related to resistance against bacterial canker, particularly through the upregulation of the RPM1 gene, which contains a specific SV, after infection by Pseudomonas syringae pv. actinidiae. In addition, we also discussed the interlineage hybridizations and taxonomic treatments of the genus Actinidia. Overall, the comprehensive pan-genome constructed here, along with our findings, lays a foundation for examining genetic compositions and markers, particularly those related to SVs, to facilitate hybrid breeding aimed at developing desired phenotypes in kiwifruits.

A chromosome-level genome assembly of Ficus benjamina, a fig tree with great ecological and ornamental value

Ficus benjamina, the weeping fig, is one of the most widely distributed and cultivated figs, with important ecological functions and landscape value. However, the lack of a reference genome has hindered molecular and functional research on this well-known fig-tree. Here we present a chromosome-scale genome assembly and annotation for F. benjamina, based on a combination of Illumina short-reads, PacBio subreads, and Hi-C sequencing data. The genome consists of 13 pseudochromosomes that contain 362.73 Mb of assembled sequences, with a contig N50 length of 25.76 Mb and a complete BUSCO score of 98.10%. In total, 28,840 protein-coding genes were identified, of which 96.22% were functionally annotated. Our study provides the first chromosome-level genome of F. benjamina, providing an important resource for exploring the genetic basis of its ecological and horticultural characters.

Chromosome-level genome assembly of the parasitoid wasp Aenasius arizonensis

Aenasius arizonensis is an important solitary endoparasitoid successfully used for biocontrol of cotton mealybug. However, lacking genomic resources has limited molecular-level investigations. Our exploration produced a superior genomic assembly of A. arizonensis from the chromosome level by combining MGISEQ short reads, Hi-C scaffolding, and PacBio Revio sequencing techniques. The genome measured 398.69 Mb, including a contig N50 of 4.73 Mb, a BUSCO completeness level of 97.07%, and a scaffold N50 of 35.96 Mb. Hi-C data were further utilized cluster and anchor 98.66% of the genome sequences into 11 chromosomes. Approximately, 165.90 Mb, representing about 41.61% of the genome, was identified as repeat elements. Non-coding sequence annotation identified 171 rRNAs, 117 small RNAs, 331 regulatory RNAs, and 872 tRNAs. Genome annotation reveals 11,727 protein-coding genes, with 10,842 (92.45%) genes functionally annotated. In summary, our chromosome-level genome assembly serves as a significant resource for advancing research on Encyrtidae parasitoids.

Comparative Genomics Reveals the Molecular Mechanisms of a Newly Isolated Pediococcus cellicola zy165 Strain and Its Adaptation in Corn Silage

Understanding how lactic acid bacteria (LAB) adapt to the silage environment is crucial for optimizing fermentation processes and developing efficient inoculants. In this study, Pediococcus cellicola zy165, isolated from fermented whole-crop corn, was subjected to whole-genome sequencing and comparative genomic analysis with two reference strains from NCBI (P. cellicola DSM 17757, and P. cellicola NBRC 106103, isolated from distilled-spirit-fermenting cellars), to elucidate its adaptation mechanisms in silage. The genome of P. cellicola zy165, which includes a circular plasmid and a CRISPR element, revealed enrichment in genes linked to carbohydrate metabolism, transport, and regulatory functions. Key adaptations for silage fermentation were evidenced by the presence of diverse phosphotransferase system (PTS) components, facilitating efficient sugar uptake and metabolism, alongside enzymes like phosphoglycerate mutase and L-lactate dehydrogenase, which are pivotal for glycolysis and lactic acid production, respectively. Additionally, the strain's genome encodes for acetate kinase, suggesting a strategic approach to pH management and energy conservation. Unique to P. cellicola zy165, genes encoding alpha-galactosidase and fructoselysine 6-phosphate deglycase were identified, indicating specialized capabilities for carbohydrate degradation in the silage niche. Structural variations and mutation analyses further highlighted adaptive genetic changes, including those in DNA metabolic processes, which could enhance survival under silage conditions. These genomic insights highlight the potential of P. cellicola zy165 as an effective silage inoculant, showcasing its evolutionary adaptations to the anaerobic, nutrient-rich corn silage environment.

DeepMFFGO: A Protein Function Prediction Method for Large-Scale Multifeature Fusion

Protein functional studies are crucial in the fields of drug target discovery and drug design. However, the existing methods have significant bottlenecks in utilizing multisource data fusion and Gene Ontology (GO) hierarchy. To this end, this study innovatively proposes the DeepMFFGO model designed for protein function prediction under large-scale multifeature fusion. A fine-tuning strategy using intermediate-level feature selection is proposed to reduce redundancy in protein sequences and mitigate distortion of the top-level features. A hierarchical progressive fusion structure is designed to explore feature connections, optimize complementarity through dynamic weight allocation, and reduce redundant interference. On the CAFA3 data set, the Fmax values of the DeepMFFGO model on the MF, BP, and CC ontologies reach 0.702, 0.599, and 0.704, respectively, which are improved by 4.2%, 2.4%, and 0.07%, respectively, compared with state-of-the-art multisource methods.

Telomere-to-telomere gapless genome assembly of Triplophysa yaopeizhii

The genus Triplophysa exhibits remarkable adaptability to the unique environment found at the Qinghai-Tibet Plateau (QTP). Higher quality genomes are helpful to the study of the adaptability to the extreme environment in the plateau. This study utilized PacBio HiFi, Ultra-long ONT, and Hi-C sequencing of Triplophysa yaopeizhii to construct the first telomere-to-telomere (T2T) gapless genome assembly of the genus Triplophysa. The genome size is 671.58 Mb, with a contig N50 length of 26.04 Mb. The sequences were anchored onto 25 chromosomes with all centromeres and telomeres. Furthermore, 293.98 Mb (43.77%) of repetitive sequences and 26,487 protein-coding genes were identified. Comparative analyses with the genomes of closely related species demonstrated high completeness, continuity, and accuracy of the genome. The genomic quality was further substantiated by the QV of 31.82 with 96.60% of BUSCO. This study provides a valuable genetic resource of the genus Triplophysa and serves as an essential reference for elucidating the adaptive genetic mechanisms of plateau fish to the high altitude.

A high-quality chromosome-level genome assembly of Antiaris toxicaria

OBJECTIVES

Antiaris toxicaria is a tall tree belonging to the Moraceae family, known for its medicinal value. Its latex contains various cardiac glycosides, which hold significant research and potential application value. However, the lack of genomic resources for A. toxicaria currently hinders molecular genetic studies on its medicinal components. For its effective conservation and elucidation of the distinctive genetic traits of and medical components, we present its chromosome-level genome assembly.

DATA DESCRIPTION

Here, we assembled two haplotypes of A. toxicaria, including a 671.73-Mb HapA subgenome containing 27,213 genes and a 666.41-Mb HapB subgenome containing 28,840 genes. Their contig N50 sizes were 90.18 and 90.29 Mb, respectively. The transposable elements represented 61.15% and 64.13% of the total assembled genome in HapA and HapB subgenome, respectively. A total of 27,213 and 28,840 genes were predicted in the two haplotypes. Hopefully, this chromosome-level genome of A. toxicaria will provide a valuable resource to enhance understanding of the biosynthesis of medicinal compounds.

Chromosome-level genome assembly of Jaguar guapote (Parachromis manguensis) by massive parallel sequencing

Parachromis managuensis is a native cichlid fish from Central America and is the most commonly traded species within its genus. This study presents the first chromosome-scale genome assembly of P. managuensis using PacBio HiFi long reads and Hi-C sequencing data. The size of the P. managuensis genome is approximately 896.66 Mb, with a scaffold N50 of 38.19 Mb. The assembled genome demonstrates high quality in terms of completeness and accuracy, with a BUSCO score of 98.85% and a quality value (QV) of 50.95. A total of 888.60 Mb (99.10%) sequences were anchored to 24 pseudochromosomes. Additionally, 21,145 protein-coding genes and 325.58 Mb (~36.31%) repetitive sequences were identified. This chromosome-level genome assembly provides a crucial reference for studying the evolution and ecological adaptability of P. managuensis.

Nitrogen metabolism of the highly ureolytic bacterium Proteus penneri S99 isolated from the rumen

BACKGROUND

The model rumen-dominant ureolytic bacterium P. penneri S99 exhibits high urease activity. It was cultivated using ammonia, urea, amino acids, or their combination as nitrogen sources. To identify differences in gene expression, the transcript abundances of various genes involved in nitrogen metabolism were analyzed by harvesting mRNA from cells during the exponential growth phases on different nitrogen sources.

RESULTS

P. penneri S99 can utilize ammonia, urea, or amino acids as the sole nitrogen sources for growth and shows a preference for utilizing urea. It exhibits similar growth rates and maximum biomass on ammonia and urea, but showed higher growth rates and maximum biomass on amino acids. Transcriptome sequencing analysis revealed different transcription patterns in response to different nitrogen sources. The urease gene expression was detected in all three different nitrogen sources, and complete hydrolysis of urea was also observed when other nitrogen sources were added to the medium containing urea. The regulation of urease in P. penneri S99 was characterized by constitutive expression, not by urea. The growth of P. penneri S99 on ammonia, ammonium acid, and urea was similar, with the only observed difference being an increase in urease transcript abundance.

CONCLUSIONS

The transcription patterns of nitrogen metabolism genes offer insights into how nitrogen is utilized in the rumen.

Characterization of agr-like Loci in Lactiplantibacillus plantarum and L. paraplantarum and Their Role in Quorum Sensing and Virulence Inhibition of Staphylococcus aureus

The pathogenicity of Staphylococcus aureus is largely regulated by the agr quorum sensing (QS) system encoded by agrBDCA, which coordinates virulence factor production through secretion and sensing of auto-inducing peptides (AIPs). agr-like systems are also present in coagulase-negative staphylococci, and several of these encode AIPs that inhibit S. aureus QS. In lactic acid bacteria, a similar locus was previously identified in Lactiplantibacillus plantarum WCSF1 termed lamBDCA. Here, we characterized the lamBDCA locus in L. plantarum LMG 13556 and L. paraplantarum CIRM-BIA 1870, and explored the effects on S. aureus QS. Notably, we found that co-cultivation with L. paraplantarum significantly inhibits S. aureus QS and hemolysin production, while less so for L. plantarum. The inhibition by L. paraplantarum was lost upon disruption of its lamBDCA locus, suggesting that the L. paraplantarum AIP mediates cross-species interference with S. aureus agr activation. Transcriptomic analysis revealed that lamBDCA in L. paraplantarum controls the expression of genes belonging to various functional categories, including stress response and metabolism. The latter includes genes encoding riboflavin (B2 vitamin) biosynthesis, which enabled the growth of the L. paraplantarum lamB mutant in the presence of roseoflavin, a toxic riboflavin analogue. Collectively, our results show that L. paraplantarum CIRM-BIA 1870 interferes with S. aureus virulence gene expression through QS suppression, and they implicate QS in the probiotic properties of L. paraplantarum.

PPI-CoAttNet: A Web Server for Protein-Protein Interaction Tasks Using a Coattention Model

Predicting protein-protein interactions (PPIs) is crucial for advancing drug discovery. Despite the proposal of numerous advanced computational methods, these approaches often suffer from poor usability for biologists and lack generalization. In this study, we designed a deep learning model based on a coattention mechanism that was capable of both PPI and site prediction and used this model as the foundation for PPI-CoAttNet, a user-friendly, multifunctional web server for PPI prediction. This platform provides comprehensive services for online PPI model training, PPI and site prediction, and prediction of interactions with proteins associated with highly prevalent cancers. In our Homo sapiens test set for PPI prediction, PPI-CoAttNet achieved an AUC of 0.9841 and an F1 score of 0.9440, outperforming most state-of-the-art models. Additionally, these results are generated in real time, delivering outcomes within minutes. We also evaluated PPI-CoAttNet for downstream tasks, including novel E3 ligase scoring, demonstrating outstanding accuracy. We believe that this tool will empower researchers, especially those without computational expertise, to leverage AI for accelerating drug development.

Chromosome-level genome assembly, annotation, and population genomic resource of argali (Ovis ammon)

Argali stands as the largest species among wild sheep in Central and East Asia, with a concerning rate of decline estimated at 30%. The intraspecific taxonomy of argali remains contentious due to limited genomic data and unclear geographic separation. In this study, we constructed a chromosome-level genome assembly and annotation for the Tibetan argali (O. a. hodgsoni), together with population genomic resequencing of 32 individuals representing four subspecies. The contig-level genome was 2.64 Gb in size, with a contig N50 length of 71.69 Mb and an estimated genomic completeness of 96.01%. Using Hi-C sequencing data scaffolding, 99.90% of initially assembled sequences were mapped and oriented onto 28 pseudo-chromosomes except the Y chromosome. Annotation uncovered 21,564 protein-coding genes and 46.38% repeat sequences. The average coverage of the population resequencing data was 23.74 with mean mapping ratio up to of 97.19%. The high-quality genome assembly and annotation of the Tibetan argali, coupled with the high-depth population genomic data, will serve as a valuable genetic resource for studies on the taxonomy and conservation of argali.

A near telomere-to-telomere genome assembly of the Jinhua pig: enabling more accurate genetic research

BACKGROUND

Pigs are crucial sources of meat and protein, valuable animal models, and potential donors for xenotransplantation. However, the existing reference genome for pigs is incomplete, with thousands of segments and centromeres and telomeres missing, which limits our understanding of the important traits in these genomic regions.

FINDINGS

We present a near-complete genome assembly for the Jinhua pig (JH-T2T) and provide a set of diploid Jinhua reference genomes, constructed using PacBio HiFi, ONT long reads, and Hi-C reads. This assembly includes all 18 autosomes and the X and Y sex chromosomes, with only 6 gaps. It features annotations of 46.90% repetitive sequences, 33 telomeres, 17 centromeres, and 23,924 high-confident genes. Compared to the Sscrofa11.1, JH-T2T closes nearly all gaps, extends sequences by 177 Mb, predicts more intact telomeres and centromeres, and gains 799 more genes and loses 114 genes. Moreover, it enhances the mapping rate for both Western and Chinese local pigs, outperforming Sscrofa11.1 as a reference genome. Additionally, this comprehensive genome assembly will facilitate large-scale variant detection.

CONCLUSIONS

This study produced a near-gapless assembly of the pig genome and provides a set of haploid Jinhua reference genomes. Our findings represent a significant advance in pig genomics, providing a robust resource that enhances genetic research, breeding programs, and biomedical applications.

Full-length transcriptome sequencing of the longissimus dorsi muscle of yak and cattle-yak using nanopore technology

Short-read RNA sequencing has been used to sequence the transcriptome of the skeletal muscle of yak and cattle-yak; however, full-length transcripts cannot be obtained and alternative splicing (AS) events cannot be inferred using this sequencing approach. Here, we used Oxford Nanopore Technologies (ONT) full-length sequencing to sequence the transcriptome of the longissimus dorsi of yak and cattle-yak. A total of 20,323 novel genes and 172,870 novel transcripts were identified, and 159,700 novel transcripts were successfully annotated. A total of 157,812 AS events, 58,073 simple sequence repeats, 57,468 complete open reading frames, 2296 transcription factors, and 20,404 lncRNAs were detected. Differentially expressed transcripts (DETs) in the longissimus dorsi muscle of yak and cattle-yak were involved in the MAPK and JAK-STAT signaling pathways related to muscle development and growth. Protein-protein interaction analysis of DETs suggested that TNNI2 might make a major contribution to differences in muscle growth and meat quality traits between yak and cattle-yak. The results have enriched the transcriptome data of dorsal muscles, providing new ideas for the study of transcriptional regulation processes, and also providing useful information for the production of higher yields of yak meat.

Chromosome-level genome assembly of the ratmouth barbel, Ptychidio jordani

The ratmouth barbel (Ptychidio jordani) is a critically endangered freshwater fish from the Cyprinidae family, primarily due to overfishing and habitat disruption. To address the challenges of its shrinking wild populations and the difficulties in artificial reproduction, we sequenced, assembled, and annotated a high-quality chromosome-level genome of P. jordani using next-generation short-read sequencing, third-generation long-read sequencing, and Hi-C sequencing. The final genome assembly was 1.14 Gb, consisting of 25 chromosomes with a contig N50 of 25.14 Mb and a scaffold N50 of 42.91 Mb. We identified 25,183 protein-coding genes, 751.75 Mb of repeats, and 19,373 ncRNAs. Methylation loci on most chromosomes ranged from 1,000 to 3,000 per 100 kb window. Gene expression levels across various tissues were analyzed, revealing 12,135 (caudal fin), 11,465 (liver), 14,438 (gill), 12,413 (heart), 8,301 (spleen), and 3,578 (kidney) differentially expressed genes compared to muscle. The comprehensive genomic and transcriptomic resources generated here will aid in understanding the ecology, adaptation, and environmental responses of P. jordani, supporting future research and conservation efforts.

Improving Antifreeze Proteins Prediction With Protein Language Models and Hybrid Feature Extraction Networks

Accurate identification of antifreeze proteins (AFPs) is crucial in developing biomimetic synthetic anti-icing materials and low-temperature organ preservation materials. Although numerous machine learning-based methods have been proposed for AFPs prediction, the complex and diverse nature of AFPs limits the prediction performance of existing methods. In this study, we propose AFP-Deep, a new deep learning method to predict antifreeze proteins by integrating embedding from protein sequences with pre-trained protein language models and evolutionary contexts with hybrid feature extraction networks. The experimental results demonstrated that the main advantage of AFP-Deep is its utilization of pre-trained protein language models, which can extract discriminative global contextual features from protein sequences. Additionally, the hybrid deep neural networks designed for protein language models and evolutionary context feature extraction enhance the correlation between embeddings and antifreeze pattern. The performance evaluation results show that AFP-Deep achieves superior performance compared to state-of-the-art models on benchmark datasets, achieving an AUPRC of 0.724 and 0.924, respectively.

The pineapple reference genome: Telomere-to-telomere assembly, manually curated annotation, and comparative analysis

Pineapple is the third most crucial tropical fruit worldwide and available in five varieties. Genomes of different pineapple varieties have been released to date; however, none of them are complete, with all exhibiting substantial gaps and representing only two of the five pineapple varieties. This significantly hinders the advancement of pineapple breeding efforts. In this study, we sequenced the genomes of three varieties: a wild pineapple variety, a fiber pineapple variety, and a globally cultivated edible pineapple variety. We constructed the first gap-free reference genome (Ref) for pineapple. By consolidating multiple sources of evidence and manually revising each gene structure annotation, we identified 26,656 protein-coding genes. The BUSCO evaluation indicated a completeness of 99.2%, demonstrating the high quality of the gene structure annotations in this genome. Utilizing these resources, we identified 7,209 structural variations across the three varieties. Approximately 30.8% of pineapple genes were located within ±5 kb of structural variations, including 30 genes associated with anthocyanin synthesis. Further analysis and functional experiments demonstrated that the high expression of AcMYB528 aligns with the accumulation of anthocyanins in the leaves, both of which may be affected by a 1.9-kb insertion fragment. In addition, we developed the Ananas Genome Database, which offers data browsing, retrieval, analysis, and download functions. The construction of this database addresses the lack of pineapple genome resource databases. In summary, we acquired a seamless pineapple reference genome with high-quality gene structure annotations, providing a solid foundation for pineapple genomics and a valuable reference for pineapple breeding.

Chromosome-level genome assembly and annotation of the Spinibarbus caldwelli

Spinibarbus caldwelli is an important freshwater economic fish in China. Owing to uncontrolled fishing, wild resources of S. caldwelli have decreased rapidly and may be on the verge of extinction. In this study, utilizing single-molecule real-time (SMRT) sequencing technology and chromatin interaction mapping (Hi-C) technologies, we assembled the first chromosome-scale genome for S. caldwelli about 1.77 Gb in size, with a contig N50 length of 11.83 Mb and scaffold N50 length of 33.91 Mb. In total 1.72 Gb (97.01%) of the contig sequences were anchored onto fifty chromosomes with the longest scaffold being 56.20 Mb. Furthermore, proximately 49.41% of the genome was composed of repetitive elements. In total, 49,377 protein-coding genes were predicted, of which 47,724 (96.65%) genes have been functionally annotated. The high-quality chromosome-level reference genome and annotation are vital for supporting basic genetic studies and will be contribute to genetic structure, functional elucidation, evolutionary inquiry, and germplasm conservation for S. caldwelli.

The Probiotic Enterococcus Lactis SF68 as a Potential Food Fermentation Microorganism for Safe Food Production

Due to the reports describing virulent and multidrug resistant enterococci, their use has become a topic of controversy despite most of them being safe and commonly used in traditionally fermented foods worldwide. We have characterized Enterococcus lactis SF68, a probiotic strain approved by the European Food Safety Authority (EFSA) for use in food and feed, and find that it has a remarkable potential in food fermentations. Genome analysis revealed the potential of SF68 to metabolize a multitude of carbohydrates, including lactose and sucrose, which was substantiated experimentally. Bacteriocin biosynthesis clusters were identified and SF68 was found to display a strong inhibitory effect against Listeria monocytogenes. Fermentation-wise, E. lactis SF68 was remarkably like Lactococcus lactis and displayed a clear mixed-acid shift on slowly fermented sugars. SF68 could produce the butter aroma compounds, acetoin and diacetyl, the production of which was enhanced under aerated conditions in a strain deficient in lactate dehydrogenase activity. Overall, E. lactis SF68 was found to be versatile, with a broad carbohydrate utilization capacity, a capacity for producing bacteriocins, and an ability to grow at elevated temperatures. This is key to eliminating pathogenic and spoilage microorganisms that are frequently associated with fermented foods.

Chromosome-level genome assembly and annotation of the Rhabdophis nuchalis (Hubei keelback)

Rhabdophis nuchalis, a snake widely distributed in China, possesses a unique trait: glands beneath the skin on its neck and back, known as nucho-dorsal glands. These features make it a valuable subject for studying genetic diversity and the evolution of complex traits. In this study, we obtained a high-quality chromosome-level reference genome of R. nuchalis using MGI short-read sequencing, PacBio Revio long-read sequencing, and Hi-C sequencing techniques. The final assembly comprised 1.92 Gb of the R. nuchalis genome, anchored to 20 chromosomes (including 9 macrochromosomes and 11 microchromosomes), with a contig N50 of 104.79 Mb, a scaffold N50 of 204.96 Mb, and a BUSCO completeness of 97.50%. Additionally, we annotated a total of 1.09 Gb of repetitive sequences (which constitute 56.51% of the entire genome) and identified 22,057 protein-coding genes. This high-quality reference genome of R. nuchalis furnishes essential genomic data for comprehending the genetic diversity and evolutionary history of the species, as well as for facilitating species conservation efforts and comparative genomics studies.

Telomere-to-telomere genome assembly of sorghum

"Cuohu Bazi" (CHBZ) is an ancient sorghum variety collected from the fields of China, known for its agronomic traits like dwarf stature, early maturation. In this study, we present the first telomere-to-telomere (T2T) and gap-free genome assembly of CHBZ using PacBio HiFi reads, Oxford Nanopore Technologies, and Hi-C data. The assembled genome comprises 724.85 Mb, effectively resolving all 3,913 gaps that were present in the previous sorghum BTx623 reference genome. Notably, the T2T assembly captures 10 centromeres and all 20 telomeres, providing strong support for their integrity. This assembly is of high quality in terms of contiguity (contig N50: 71.1 Mb), completeness (BUSCO score: 99.01%, k-mer completeness: 98.88%), and correctness (QV: 61.60). Repetitive sequences accounted for 70.41% of the genome and a total of 32,855 protein-coding genes have been annotated. Furthermore, 161 CHBZ-specific presence/absence variants genes have been identified when comparing to BTx623 genome. This study provides valuable insights for future research on sorghum genetics, genomics, and evolutionary history.

Whole genome analysis of Pseudomonas mandelii SW-3 and the insights into low-temperature adaptation

Pseudomonas mandelii SW-3, isolated from the Napahai plateau wetland, can survive in cold environments. The mechanisms underlying the survival of bacteria in low temperatures and high altitudes are not yet fully understood. In this study, the whole genome of SW-3 was sequenced to identify the genomic features that may contribute to survival in cold environments. The results showed that the genome size of strain SW-3 was 6,538,059 bp with a GC content of 59%. A total of 67 tRNAs, a 34,110 bp prophage sequence, and a large number of metabolic genes were found. Based on 16S rRNA gene phylogeny and average nucleotide identity analysis among P. mandelii, SW-3 was identified as a strain belonging to P. mandelii. In addition, we clarified the mechanisms by which SW-3 survived in a cold environment, providing a basis for further investigation of host-phage interaction. P. mandelii SW-3 showed stress resistance mechanisms, including glycogen and trehalose metabolic pathways, and antisense transcriptional silencing. Furthermore, cold shock proteins and glucose 6-phosphate dehydrogenase may play pivotal roles in facilitating adaptation to cold environmental conditions. The genome-wide analysis provided us with a deeper understanding of the cold-adapted bacterium.

Full-length agave transcriptome reveals candidate glycosyltransferase genes involved in hemicellulose biosynthesis

Agave species are typical crassulacean acid metabolism (CAM) plants commonly cultivated to produce beverages, fibers, and medicines. To date, few studies have examined hemicellulose biosynthesis in Agave H11648, which is the primary cultivar used for fiber production. We conducted PacBio sequencing to obtain full-length transcriptome of five agave tissues: leaves, shoots, roots, flowers, and fruits. A total of 41,807 genes were generated, with a mean length of 2394 bp and an annotation rate of 97.12 % using public databases. We identified 42 glycosyltransferase genes related to hemicellulose biosynthesis, including mixed-linkage glucan (1), glucomannan (5), xyloglucan (16), and xylan (20). Their expression patterns were examined during leaf development and fungal infection, together with hemicellulose content. The results revealed four candidate glycosyltransferase genes involved in xyloglucan and xylan biosynthesis, including glucan synthase (CSLC), xylosyl transferase (XXT), xylan glucuronyltransferase (GUX), and xylan α-1,3-arabinosyltransferase (XAT). These genes can be potential targets for manipulating xyloglucan and xylan traits in agaves, and can also be used as candidate enzymatic tools for enzyme engineering. We have provided the first full-length transcriptome of agave, which will be a useful resource for gene identification and characterization in agave species. We also elucidated the hemicellulose biosynthesis machinery, which will benefit future studies on hemicellulose traits in agave.

A near complete genome assembly of the East Friesian sheep genome

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

An immuno-informatics approach for annotation of hypothetical proteins and multi-epitope vaccine designed against the Mpox virus

A worrying new outbreak of Monkeypox (Mpox) in humans is caused by the Mpox virus (MpoxV). The pathogen has roughly 28 hypothetical proteins of unknown structure, function, and pathogenicity. Using reliable bioinformatics tools, we attempted to analyze the MpoxV genome, identify the role of hypothetical proteins (HPs), and design a potential candidate vaccine. Out of 28, we identified seven hypothetical proteins using multi-server validation with high confidence for the occurrence of conserved domains. Their physical, chemical, and functional characterizations, including molecular weight, theoretical isoelectric point, 3D structures, GRAVY value, subcellular localization, functional motifs, antigenicity, and virulence factors, were performed. We predicted possible cytotoxic T cell (CTL), helper T cell (HTL) and linear and conformational B cell epitopes, which were combined in a 219 amino acid multiepitope vaccine with human β defensin as a linker. This multi-epitopic vaccine was structurally modelled and docked with toll-like receptor-3 (TLR-3). The dynamical stability of the vaccine-TLR-3 docked complexes exhibited stable interactions based on RMSD and RMSF tests. Additionally, the modelled vaccine was cloned in-silico in an E. coli host to check the appropriate expression of the final vaccine built. Our results might conform to an immunogenic and safe vaccine, which would require further experimental validation.Communicated by Ramaswamy H. Sarma.

Haplotype-resolved assembly of a pig genome using single-sperm sequencing

Single gamete cell sequencing together with long-read sequencing can reliably produce chromosome-level phased genomes. In this study, we employed PacBio HiFi and Hi-C sequencing on a male Landrace pig, coupled with single-sperm sequencing of its 102 sperm cells. A haplotype assembly method was developed based on long-read sequencing and sperm-phased markers. The chromosome-level phased assembly showed higher phasing accuracy than methods that rely only on HiFi reads. The use of single-sperm sequencing data enabled the construction of a genetic map, successfully mapping the sperm motility trait to a specific region on chromosome 1 (105.40-110.70 Mb). Furthermore, with the assistance of Y chromosome-bearing sperm data, 26.16 Mb Y chromosome sequences were assembled. We report a reliable approach for assembling chromosome-level phased genomes and reveal the potential of sperm population in basic biology research and sperm phenotype research.

Evaluation of surfactant-aided polycyclic aromatic hydrocarbon biodegradation by molecular docking and molecular dynamic simulation in the marine environment

Marine oil spills directly cause polycyclic aromatic hydrocarbons (PAHs) pollution and affect marine organisms due to their toxic property. Chemical and bio-based dispersants composed of surfactants and solvents are considered effective oil spill-treating agents. Dispersants enhance oil biodegradation in the marine environment by rapidly increasing their solubility in the water column. However, the effect of dispersants, especially surfactants, on PAHs degradation by enzymes produced by microorganisms has not been studied at the molecular level. The role of the cytochrome P450 (CYP) enzyme in converting contaminants into reactive metabolites during the biodegradation process has been evidenced, but the activity in the presence of surfactants is still ambiguous. Thus, this study focused on the evaluation of the impact of chemical and bio-surfactants (i.e., Tween 80 (TWE) and Surfactin (SUC)) on the biodegradation of naphthalene (NAP), chrysene (CHR), and pyrene (PYR), the representative components of PAHs, with CYP enzyme from microalgae Parachlorella kessleri using molecular docking and molecular dynamics (MD) simulation. The molecular docking analysis revealed that PAHs bound to residues at the CYP active site through hydrophobic interactions for biodegradation. The MD simulation showed that the surfactant addition changed the enzyme conformation in the CYP-PAH complexes to provide more interactions between the enzyme and PAHs. This led to an increase in the enzyme's capability to degrade PAHs. Binding free energy (ΔG‌‌Bind) calculations confirmed that surfactant treatment could enhance PAHs degradation by the enzyme. The SUC gave a better result on NAP and PYR biodegradation based on ΔG‌‌Bind, while TWE facilitated the biodegradation of CHR. The research outputs could greatly facilitate evaluating the behaviors of oil spill-treating agents and oil spill response operations in the marine environment.

A high-quality chromosome-level genome assembly of Ficus hirta

Ficus species (Moraceae) play pivotal roles in tropical and subtropical ecosystems. Thriving across diverse habitats, from rainforests to deserts, they harbor a multitude of mutualistic and antagonistic interactions with insects, nematodes, and pathogens. Despite their ecological significance, knowledge about the genomic background of Ficus remains limited. In this study, we report a chromosome-level reference genome of F. hirta, with a total size of 297.27 Mb, containing 28,625 protein-coding genes and 44.67% repeat sequences. These findings illuminate the genetic basis of Ficus responses to environmental challenges, offering valuable genomic resources for understanding genome size, adaptive evolution, and co-evolution with natural enemies and mutualists within the genus.

Comprehensive Research on Druggable Proteins: From PSSM to Pre-Trained Language Models

Identification of druggable proteins can greatly reduce the cost of discovering new potential drugs. Traditional experimental approaches to exploring these proteins are often costly, slow, and labor-intensive, making them impractical for large-scale research. In response, recent decades have seen a rise in computational methods. These alternatives support drug discovery by creating advanced predictive models. In this study, we proposed a fast and precise classifier for the identification of druggable proteins using a protein language model (PLM) with fine-tuned evolutionary scale modeling 2 (ESM-2) embeddings, achieving 95.11% accuracy on the benchmark dataset. Furthermore, we made a careful comparison to examine the predictive abilities of ESM-2 embeddings and position-specific scoring matrix (PSSM) features by using the same classifiers. The results suggest that ESM-2 embeddings outperformed PSSM features in terms of accuracy and efficiency. Recognizing the potential of language models, we also developed an end-to-end model based on the generative pre-trained transformers 2 (GPT-2) with modifications. To our knowledge, this is the first time a large language model (LLM) GPT-2 has been deployed for the recognition of druggable proteins. Additionally, a more up-to-date dataset, known as Pharos, was adopted to further validate the performance of the proposed model.

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

A reference transcriptome for walnut anthracnose pathogen, Ophiognomonia leptostyla, guides the discovery of candidate virulence genes

Despite the economic losses due to the walnut anthracnose, Ophiognomonia leptostyla is an orphan fungus with respect to genomic resources. In the present study, the transcriptome of O. leptostyla was assembled for the first time. RNA sequencing was conducted for the fungal mycelia grown in a liquid media, and the inoculated leaf samples of walnut with the fungal conidia sampled at 48, 96 and 144 h post inoculation (hpi). The completeness, correctness, and contiguity of the de novo transcriptome assemblies generated with Trinity, Oases, SOAPdenovo-Trans and Bridger were compared to identify a single superior reference assembly. In most of the assessment criteria including N50, Transrate score, number of ORFs with known description in gene bank, the percentage of reads mapped back to the transcript (RMBT), BUSCO score, Swiss-Prot coverage bin and RESM-EVAL score, the Bridger assembly was the superior and thus used as a reference for profiling the O. leptostyla transcriptome in liquid media vs. during walnut infection. The k-means clustering of transcripts resulted in four distinct transcription patterns across the three sampling time points. Most of the detected CAZy transcripts had elevated transcription at 96 hpi that is hypothetically concurrent with the start of intracellular growth. The in-silico analysis revealed 103 candidate effectors of which six were members of Necrosis and Ethylene Inducing Like Protein (NLP) gene family belonging to three distinct k-means clusters. This study provided a complex and temporal pattern of the CAZys and candidate effectors transcription during six days post O. leptostyla inoculation on walnut leaves, introducing a list of candidate virulence genes for validation in future studies.

Multi-omics reveal mechanisms underlying chronic kidney disease of unknown etiology (CKDu) pathogenesis using zebrafish

Chronic kidney disease of unknown etiology (CKDu) is an endemic disease in the dry zone of farming communities, Sri Lanka. The drinking water in a CKDu prevalent area contains a high concentration of F-, hardness and other environmental pollutants, including heavy metals and microcystin, which are considered possible etiology of CKDu in these areas. Here, multi-omics data with host transcriptome, metabolome and gut microbiomes were obtained using simulated local drinking water of Sri Lanka after their exposure to adult zebrafish. Based on an integrated multi-omics analysis in the context of host physiology in the kidney injury samples with different pathologic grades, two common pathways necroptosis and purine metabolism were identified as potentially important pathways that affect kidney injury. The key metabolite acetyl adenylate in the purine metabolism pathway was significantly positively correlated with Comamonas (rho = 0.72) and significantly negatively correlated with Plesiomonas (rho = -0.58). This crucial metabolite and two key gut bacteria genera may not only be potential markers but also potential therapeutic targets in the uric acid metabolic pathway, which is an important factor in the pathogenesis of acute kidney injury (AKI) in general, as well as of chronic kidney disease (CKD). Based on this, we revealed the urea metabolism pathway of kidney injury in zebrafish and provided a new avenue for the treatment of CKDu in Sri Lanka.

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

Chromosome-Level Genome Assembly of Protosalanx chinensis and Response to Air Exposure Stress

Protosalanx chinensis is a suitable particular species for genetic studies on nearly scaleless skin, transparency and high sensitivity to hypoxia stress. Here, we generated a high-quality chromosome-level de novo assembly of P. chinensis. The final de novo assembly yielded 379.47 Mb with 28 pseudo-chromosomes and a scaffold N50 length of 14.52 Mb. In total, 21,074 protein-coding genes were predicted. P. chinensis, Esox lucius and Hypomesus transpacificus had formed a clade, which diverged about 115.5 million years ago. In the air exposure stress experiment, we found that some genes play an essential role during P. chinensis hypoxia, such as bhlh, Cry1, Clock, Arntl and Rorb in the circadian rhythm pathway. These genomic data offer a crucial foundation for P. chinensis ecology and adaptation studies, as well as a deeper understanding of the response to air exposure stress.

Molecular mechanism of thiamine in mitigating drought stress in Chinese wingnut (Pterocarya stenoptera): Insights from transcriptomics

Urban garden plants are frequently affected by drought, which can hinder their growth, development, and greening effect. Previous studies have indicated that Chinese wingnut (Pterocarya stenoptera) responds to drought stress by increasing the expression of thiamine synthesis genes. In this study, it was found that exogenous thiamine can effectively alleviate the negative effects of drought stress on plants. Forward transcriptome sequencing and physiological tests were further conducted to reveal the molecular mechanism of thiamine in alleviating drought stress. Results showed that exogenous thiamine activated the expression of eight chlorophyll synthesis genes in Chinese wingnut under drought stress. Moreover, physiological indicators proved that chlorophyll content increased in leaves of Chinese wingnut with thiamine treatment under drought stress. Photosynthesis genes were also activated in Chinese wingnut treated with exogenous thiamine under drought stress, as supported by photosynthetic indicators PIabs and PItotal. Additionally, exogenous thiamine stimulated the expression of genes in the auxin-activated signaling pathway, thus attenuating the effects of drought stress. This study demonstrates the molecular mechanism of thiamine in mitigating the effects of drought stress on non-model woody plants lacking transgenic systems. This study also provides an effective method to mitigate the negative impacts of drought stress on plants.

Cell-type specific defects in PTEN-mutant cortical organoids converge on abnormal circuit activity

De novo heterozygous loss-of-function mutations in phosphatase and tensin homolog (PTEN) are strongly associated with autism spectrum disorders; however, it is unclear how heterozygous mutations in this gene affect different cell types during human brain development and how these effects vary across individuals. Here, we used human cortical organoids from different donors to identify cell-type specific developmental events that are affected by heterozygous mutations in PTEN. We profiled individual organoids by single-cell RNA-seq, proteomics and spatial transcriptomics and revealed abnormalities in developmental timing in human outer radial glia progenitors and deep-layer cortical projection neurons, which varied with the donor genetic background. Calcium imaging in intact organoids showed that both accelerated and delayed neuronal development phenotypes resulted in similar abnormal activity of local circuits, irrespective of genetic background. The work reveals donor-dependent, cell-type specific developmental phenotypes of PTEN heterozygosity that later converge on disrupted neuronal activity.

Whole genome analysis and cold adaptation strategies of Pseudomonas sivasensis W-6 isolated from the Napahai plateau wetland

Microbial communities of wetlands play key roles in the earth's ecology and stability. To elucidate the cold adaptation mechanisms of bacteria in plateau wetlands, we conducted comparative genomic analyses of Pseudomonas sivasensis and closely related lineages. The genome of P. sivasensis W-6, a cold-adapted bacterium isolated from the Napahai plateau wetland, was sequenced and analyzed. The genome length was 6,109,123 bp with a G+C content of 59.5%. Gene prediction yielded 5360 protein-coding sequences, 70 tRNAs, 24 gene islands, and 2 CRISPR sequences. The isolate contained evidence of horizontal gene transfer events during its evolution. Two prophages were predicted and indicated that W-6 was a lysogen. The cold adaptation of the W-6 strain showed psychrophilic rather than psychrotrophic characteristics. Cold-adapted bacterium W-6 can utilize glycogen and trehalose as resources, associated with carbohydrate-active enzymes, and survive in a low-temperature environment. In addition, the cold-adapted mechanisms of the W-6 included membrane fluidity by changing the unsaturated fatty acid profile, the two-component regulatory systems, anti-sense transcription, the role played by rpsU genes in the translation process, etc. The genome-wide analysis of W-6 provided a deeper understanding of cold-adapted strategies of bacteria in environments. We elucidated the adaptive mechanism of the psychrophilic W-6 strain for survival in a cold environment, which provided a basis for further study on host-phage coevolution.

Functional viromic screens uncover regulatory RNA elements

The number of sequenced viral genomes has surged recently, presenting an opportunity to understand viral diversity and uncover unknown regulatory mechanisms. Here, we conducted a screening of 30,367 viral segments from 143 species representing 96 genera and 37 families. Using a library of viral segments in 3' UTR, we identified hundreds of elements impacting RNA abundance, translation, and nucleocytoplasmic distribution. To illustrate the power of this approach, we investigated K5, an element conserved in kobuviruses, and found its potent ability to enhance mRNA stability and translation in various contexts, including adeno-associated viral vectors and synthetic mRNAs. Moreover, we identified a previously uncharacterized protein, ZCCHC2, as a critical host factor for K5. ZCCHC2 recruits the terminal nucleotidyl transferase TENT4 to elongate poly(A) tails with mixed sequences, delaying deadenylation. This study provides a unique resource for virus and RNA research and highlights the potential of the virosphere for biological discoveries.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

Transcriptome Sequencing of Agave amaniensis Reveals Shoot-Related Expression Patterns of Expansin A Genes in Agave

Agave species are widely planted for fiber production. However, the molecular basis of agave fiber development has not been well understood. In this study, we performed a transcriptomic analysis in A. amaniensi, a well-known variety with high-quality fiber production. Approximately 43.87 million clean reads were obtained using Illumina sequencing. The de novo assembly produced 66,746 unigrams, 54% of which were annotated in a public database. In the Nr database, 21,490 unigenes of A. amaniensis were shown to be most closely related to Asparagus officinalis. Nine expansin A orthologs with full coding regions were obtained, which were named EXP1a, EXP1b, EXP2, EXP3, EXP4a, EXP4b, EXP11, EXP12, and EXP13. The maximum likelihood phylogenetic tree revealed the species-specific expansion of expansin genes in Arabidopsis, rice and agave. The expression analysis suggested the negative correlation between the expression of expansin genes and the leaf growth rate, except AhEXP11. Moreover, expansin genes were differentially affected by abiotic and biotic stresses. Notably, AhEXP2 expression level was highly upgraded after the infection of Phytophthora nicotiana. Nutrient deficiency also influent expansin genes expression. Together, our research will benefit future studies related to fiber development, disease resistance and nutrient usage in agave.

Deciphering the biodesulfurization pathway employing marine mangrove Bacillus aryabhattai strain NM1-A2 according to whole genome sequencing and transcriptome analyses

In the biogeochemical cycle, sulfur oxidation plays a vital role and is typically referred to as the elemental sulfur or reductive sulfide oxidation process. This study aimed to characterize a subtropical mangrove-isolated bacterial strain using biochemical, whole-genome, and transcriptome sequencing analyses to enhance our understanding of sulfur metabolism and biodegradation from a molecular genetic perspective. Strain NM1-A2 was characterized as Gram-positive and found to have a close molecular phylogenetic relationship with Bacillus aryabhattai. NM1-A2 efficiently converted dibenzothiophene (DBT) into 2-hydroxybiphenyl (2-HBP) via a 4S pathway with 95% efficiency, using enzymes encoded by the dsz operon (dszA, dszB, and dszC), which determine monooxygenases (DszA & DszC) and desulfinase (DszB). The whole-genome sequence of NM1-A2 had a length of approximately 5,257,678 bp and included 16 sulfur metabolism-related genes, featuring the ABC transport system, small subunit (ssu) and cysteine (cys) gene families, and adenosine 5'-phosphosulfate (APS) and 3'-phosphoadenosine-5'-phosphosulfate (PAPS) biosynthesis-related genes. Transcriptomic analysis of NM1-A2 using three sulfur groups-magnesium sulfate (MS), sulfur powder (SP), and sodium thiosulfate (ST) resulted in a significant number of differentially expressed genes (1200, 2304, and 2001, respectively). This analysis revealed that intracellular cysteine concentration directly regulated the expression of cys and ssu genes. Sulfate did not directly affect cys gene expression but repressed ssu gene expression. The cys gene expression levels decreased during the conversion of sulfate to sulfide and cysteine. The transcriptomic data was validated by analyzing the expression patterns of NM1-A2 using real-time quantitative PCR validation analysis. The expression levels of cysl, mccB, and nrnA were significantly upregulated, while cysH, metB, and sat were downregulated in the SP, ST, and MS groups, respectively. This research contributes to our understanding of marine mangrove microorganisms' bacterial efficiency through characterization, whole-genome, and transcriptome sequencing-based molecular degradation of organic compounds in the mangrove ecosystem, which may enhance nutrient availability.

Understanding the Relationship of the Human Bacteriome with COVID-19 Severity and Recovery

The Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) first emerged in 2019 in China and has resulted in millions of human morbidities and mortalities across the globe. Evidence has been provided that this novel virus originated in animals, mutated, and made the cross-species jump to humans. At the time of this communication, the Coronavirus disease (COVID-19) may be on its way to an endemic form; however, the threat of the virus is more for susceptible (older and immunocompromised) people. The human body has millions of bacterial cells that influence health and disease. As a consequence, the bacteriomes in the human body substantially influence human health and disease. The bacteriomes in the body and the immune system seem to be in constant association during bacterial and viral infections. In this review, we identify various bacterial spp. In major bacteriomes (oral, nasal, lung, and gut) of the body in healthy humans and compare them with dysbiotic bacteriomes of COVID-19 patients. We try to identify key bacterial spp. That have a positive effect on the functionality of the immune system and human health. These select bacterial spp. Could be used as potential probiotics to counter or prevent COVID-19 infections. In addition, we try to identify key metabolites produced by probiotic bacterial spp. That could have potential anti-viral effects against SARS-CoV-2. These metabolites could be subject to future therapeutic trials to determine their anti-viral efficacies.

A complete gap-free diploid genome in Saccharum complex and the genomic footprints of evolution in the highly polyploid Saccharum genus

A diploid genome in the Saccharum complex facilitates our understanding of evolution in the highly polyploid Saccharum genus. Here we have generated a complete, gap-free genome assembly of Erianthus rufipilus, a diploid species within the Saccharum complex. The complete assembly revealed that centromere satellite homogenization was accompanied by the insertions of Gypsy retrotransposons, which drove centromere diversification. An overall low rate of gene transcription was observed in the palaeo-duplicated chromosome EruChr05 similar to other grasses, which might be regulated by methylation patterns mediated by homologous 24 nt small RNAs, and potentially mediating the functions of many nucleotide-binding site genes. Sequencing data for 211 accessions in the Saccharum complex indicated that Saccharum probably originated in the trans-Himalayan region from a diploid ancestor (x = 10) around 1.9-2.5 million years ago. Our study provides new insights into the origin and evolution of Saccharum and accelerates translational research in cereal genetics and genomics.

Complet+: a computationally scalable method to improve completeness of large-scale protein sequence clustering

A major challenge for clustering algorithms is to balance the trade-off between homogeneity, i.e., the degree to which an individual cluster includes only related sequences, and completeness, the degree to which related sequences are broken up into multiple clusters. Most algorithms are conservative in grouping sequences with other sequences. Remote homologs may fail to be clustered together and instead form unnecessarily distinct clusters. The resulting clusters have high homogeneity but completeness that is too low. We propose Complet+, a computationally scalable post-processing method to increase the completeness of clusters without an undue cost in homogeneity. Complet+ proves to effectively merge closely-related clusters of protein that have verified structural relationships in the SCOPe classification scheme, improving the completeness of clustering results at little cost to homogeneity. Applying Complet+ to clusters obtained using MMseqs2's clusterupdate achieves an increased V-measure of 0.09 and 0.05 at the SCOPe superfamily and family levels, respectively. Complet+ also creates more biologically representative clusters, as shown by a substantial increase in Adjusted Mutual Information (AMI) and Adjusted Rand Index (ARI) metrics when comparing predicted clusters to biological classifications. Complet+ similarly improves clustering metrics when applied to other methods, such as CD-HIT and linclust. Finally, we show that Complet+ runtime scales linearly with respect to the number of clusters being post-processed on a COG dataset of over 3 million sequences. Code and supplementary information is available on Github: https://github.com/EESI/Complet-Plus.

Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. 'Donghong' (DH) and Actinidia latifolia 'Kuoye' (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Chromosome-level genome and population genomics reveal evolutionary characteristics and conservation status of Chinese indigenous geese

Geese are herbivorous birds that play an essential role in the agricultural economy. We construct the chromosome-level genome of a Chinese indigenous goose (the Xingguo gray goose, XGG; Anser cygnoides) and analyze the adaptation of fat storage capacity in the goose liver during the evolution of Anatidae. Genomic resequencing of 994 geese is used to investigate the genetic relationships of geese, which supports the dual origin of geese (Anser cygnoides and Anser anser). Chinese indigenous geese show higher genetic diversity than European geese, and a scientific conservation program can be established to preserve genetic variation for each breed. We also find that a 14-bp insertion in endothelin receptor B subtype 2 (EDNRB2) that determines the white plumage of Chinese domestic geese is a natural mutation, and the linkaged alleles rapidly increase in frequency as a result of genetic hitchhiking, leading to the formation of completely different haplotypes of white geese under strong artificial selection. These genomic resources and our findings will facilitate marker-assisted breeding of geese and provide a foundation for further research on geese genetics and evolution.