InterProScan: protein domains identifier

pyRBDome: a comprehensive computational platform for enhancing RNA-binding proteome data

High-throughput proteomics approaches have revolutionised the identification of RNA-binding proteins (RBPome) and RNA-binding sequences (RBDome) across organisms. Yet, the extent of noise, including false positives, associated with these methodologies, is difficult to quantify as experimental approaches for validating the results are generally low throughput. To address this, we introduce pyRBDome, a pipeline for enhancing RNA-binding proteome data in silico. It aligns the experimental results with RNA-binding site (RBS) predictions from distinct machine-learning tools and integrates high-resolution structural data when available. Its statistical evaluation of RBDome data enables quick identification of likely genuine RNA-binders in experimental datasets. Furthermore, by leveraging the pyRBDome results, we have enhanced the sensitivity and specificity of RBS detection through training new ensemble machine-learning models. pyRBDome analysis of a human RBDome dataset, compared with known structural data, revealed that although UV-cross-linked amino acids were more likely to contain predicted RBSs, they infrequently bind RNA in high-resolution structures. This discrepancy underscores the limitations of structural data as benchmarks, positioning pyRBDome as a valuable alternative for increasing confidence in RBDome datasets.

Upstream regulator of genomic imprinting in rice endosperm is a small RNA-associated chromatin remodeler

Genomic imprinting is observed in endosperm, a placenta-like seed tissue, where transposable elements (TEs) and repeat-derived small RNAs (sRNAs) mediate epigenetic changes in plants. In imprinting, uniparental gene expression arises due to parent-specific epigenetic marks on one allele but not on the other. The importance of sRNAs and their regulation in endosperm development or in imprinting is poorly understood in crops. Here we show that a previously uncharacterized CLASSY (CLSY)-family chromatin remodeler named OsCLSY3 is essential for rice endosperm development and imprinting, acting as an upstream player in the sRNA pathway. Comparative transcriptome and genetic analysis indicated its endosperm-preferred expression and its likely paternal imprinted nature. These important features are modulated by RNA-directed DNA methylation (RdDM) of tandemly arranged TEs in its promoter. Upon perturbation of OsCLSY3 in transgenic lines, we observe defects in endosperm development and a loss of around 70% of all sRNAs. Interestingly, well-conserved endosperm-specific sRNAs (siren) that are vital for reproductive fitness in angiosperms are also dependent on OsCLSY3. We observed that many imprinted genes and seed development-associated genes are under the control of OsCLSY3. These results support an essential role of OsCLSY3 in rice endosperm development and imprinting, and propose similar regulatory strategies involving CLSY3 homologs among other cereals.

Chromosome-scale genome assembly and annotation of the tetraploid potato cultivar Diacol Capiro adapted to the Andean region

Potato (Solanum tuberosum) is an essential crop for food security and is ranked as the third most important crop worldwide for human consumption. The Diacol Capiro cultivar holds the dominant position in Colombian cultivation, primarily catering to the food processing industry. This highly heterozygous, autotetraploid cultivar belongs to the Andigenum group and it stands out for its adaptation to a wide variety of environments spanning altitudes from 1,800 to 3,200 meters above sea level. Here, a chromosome-scale assembly, referred to as DC, is presented for this cultivar. The assembly was generated by combining circular consensus sequencing with proximity ligation Hi-C for the scaffolding and represents 2.369 Gb with 48 pseudochromosomes covering 2,091 Gb and an anchor rate of 88.26%. The reference genome metrics, including an N50 of 50.5 Mb, a BUSCO (Benchmarking Universal Single-Copy Orthologue) score of 99.38%, and an Long Terminal Repeat Assembly Index score of 13.53, collectively signal the achieved high assembly quality. A comprehensive annotation yielded a total of 154,114 genes, and the associated BUSCO score of 95.78% for the annotated sequences attests to their completeness. The number of predicted NLR (Nucleotide-Binding and Leucine-Rich-Repeat genes) was 2107 with a large representation of NBARC (for nucleotide binding domain shared by Apaf-1, certain R gene products, and CED-4) containing domains (99.85%). Further comparative analysis of the proposed annotation-based assembly with high-quality known potato genomes, showed a similar genome metrics with differences in total gene numbers related to the ploidy status. The genome assembly and annotation of DC presented in this study represent a valuable asset for comprehending potato genetics. This resource aids in targeted breeding initiatives and contributes to the creation of enhanced, resilient, and more productive potato varieties, particularly beneficial for countries in Latin America.

HY5 and COP1 function antagonistically in the regulation of nicotine biosynthesis in Nicotiana tabacum

Nicotine constitutes approximately 90% of the total alkaloid content in leaves within the Nicotiana species, rendering it the most prevalent alkaloid. While the majority of genes responsible for nicotine biosynthesis express in root tissue, the influence of light on this process through shoot-to-root mobile ELONGATED HYPOCOTYL 5 (HY5) has been recognized. CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1), a key regulator of light-associated responses, known for its role in modulating HY5 accumulation, remains largely unexplored in its relationship to light-dependent nicotine accumulation. Here, we identified NtCOP1, a COP1 homolog in Nicotiana tabacum, and demonstrated its ability to complement the cop1-4 mutant in Arabidopsis thaliana at molecular, morphological, and biochemical levels. Through the development of NtCOP1 overexpression (NtCOP1OX) plants, we observed a significant reduction in nicotine and flavonol content, inversely correlated with the down-regulation of nicotine and phenylpropanoid pathway. Conversely, CRISPR/Cas9-based knockout mutant plants (NtCOP1CR) exhibited an increase in nicotine levels. Further investigations, including yeast-two hybrid assays, grafting experiments, and Western blot analyses, revealed that NtCOP1 modulates nicotine biosynthesis by targeting NtHY5, thereby impeding its transport from shoot-to-root. We conclude that the interplay between HY5 and COP1 functions antagonistically in the light-dependent regulation of nicotine biosynthesis in tobacco.

The study of iron- and copper-binding proteome of Fusarium oxysporum and its effector candidates

Fusarium oxysporum f.sp. lycopersici is a phytopathogen which causes vascular wilt disease in tomato plants. The survival tactics of both pathogens and hosts depend on intricate interactions between host plants and pathogenic microbes. Iron-binding proteins (IBPs) and copper-binding proteins (CBPs) play a crucial role in these interactions by participating in enzyme reactions, virulence, metabolism, and transport processes. We employed high-throughput computational tools at the sequence and structural levels to investigate the IBPs and CBPs of F. oxysporum. A total of 124 IBPs and 37 CBPs were identified in the proteome of Fusarium. The ranking of amino acids based on their affinity for binding with iron is Glu > His> Asp > Asn > Cys, and for copper is His > Asp > Cys respectively. The functional annotation, determination of subcellular localization, and Gene Ontology analysis of these putative IBPs and CBPs have unveiled their potential involvement in a diverse array of cellular and biological processes. Three iron-binding glycosyl hydrolase family proteins, along with four CBPs with carbohydrate-binding domains, have been identified as potential effector candidates. These proteins are distinct from the host Solanum lycopersicum proteome. Moreover, they are known to be located extracellularly and function as enzymes that degrade the host cell wall during pathogen-host interactions. The insights gained from this report on the role of metal ions in plant-pathogen interactions can help develop a better understanding of their fundamental biology and control vascular wilt disease in tomato plants.

Spirocerca lupi draft genome, vaccine and anthelmintic targets

Spirocerca lupi is a parasitic nematode affecting predominantly domestic dogs. It causes spirocercosis, a disease that is often fatal. The assembled draft genome of S. lupi consists of 13,627 predicted protein-coding genes and is approximately 150 Mb in length. Several known anthelmintic gene targets such as for β-Tubulin, glutamate, and GABA receptors as well as known vaccine gene targets such as cysteine protease inhibitor and cytokines were identified in S. lupi by comparing orthologs of C. elegans anthelmintic gene targets as well as orthologs to known vaccine candidates. New anthelmintic targets were predicted through an inclusion-exclusion strategy and new vaccine targets were predicted through an immunoinformatics approach. New anthelminthic targets include DNA-directed RNA polymerases, chitin synthase, polymerases, and other enzymes. New vaccine targets include cuticle collagens. These gene targets provide a starting platform for new drug identification and vaccine design.

Diverse RNA viruses discovered in multiple seagrass species

Seagrasses are marine angiosperms that form highly productive and diverse ecosystems. These ecosystems, however, are declining worldwide. Plant-associated microbes affect critical functions like nutrient uptake and pathogen resistance, which has led to an interest in the seagrass microbiome. However, despite their significant role in plant ecology, viruses have only recently garnered attention in seagrass species. In this study, we produced original data and mined publicly available transcriptomes to advance our understanding of RNA viral diversity in Zostera marina, Zostera muelleri, Zostera japonica, and Cymodocea nodosa. In Z. marina, we present evidence for additional Zostera marina amalgavirus 1 and 2 genotypes, and a complete genome for an alphaendornavirus previously evidenced by an RNA-dependent RNA polymerase gene fragment. In Z. muelleri, we present evidence for a second complete alphaendornavirus and near complete furovirus. Both are novel, and, to the best of our knowledge, this marks the first report of a furovirus infection naturally occurring outside of cereal grasses. In Z. japonica, we discovered genome fragments that belong to a novel strain of cucumber mosaic virus, a prolific pathogen that depends largely on aphid vectoring for host-to-host transmission. Lastly, in C. nodosa, we discovered two contigs that belong to a novel virus in the family Betaflexiviridae. These findings expand our knowledge of viral diversity in seagrasses and provide insight into seagrass viral ecology.

DeepMineLys: Deep mining of phage lysins from human microbiome

Vast shotgun metagenomics data remain an underutilized resource for novel enzymes. Artificial intelligence (AI) has increasingly been applied to protein mining, but its conventional performance evaluation is interpolative in nature, and these trained models often struggle to extrapolate effectively when challenged with unknown data. In this study, we present a framework (DeepMineLys [deep mining of phage lysins from human microbiome]) based on the convolutional neural network (CNN) to identify phage lysins from three human microbiome datasets. When validated with an independent dataset, our method achieved an F1-score of 84.00%, surpassing existing methods by 20.84%. We expressed 16 lysin candidates from the top 100 sequences in E. coli, confirming 11 as active. The best one displayed an activity 6.2-fold that of lysozyme derived from hen egg white, establishing it as the most potent lysin from the human microbiome. Our study also underscores several important issues when applying AI to biology questions. This framework should be applicable for mining other proteins.

Parallel loss of sexual reproduction in field populations of a brown alga sheds light on the mechanisms underlying the emergence of asexuality

Sexual reproduction is widespread, but asexual lineages have repeatedly arisen from sexual ancestors across a wide range of eukaryotic taxa. The molecular changes underpinning the switch to asexuality remain elusive, particularly in organisms with haploid sexual systems. Here we explore independent events of loss of sex in the brown alga Scytosiphon, examine the proximate and evolutionary mechanisms involved, and test the importance of sexual conflict on gene expression changes following loss of sex. We find that asexual females ('Amazons') lose ability to produce sex pheromone and, consequently, are incapable of attracting males, whereas they gain rapid parthenogenic development from large, unfertilized eggs. These phenotypic changes are accompanied by convergent changes in gene expression. Decay of female functions, rather than relaxation of sexual antagonism, may be a dominant force at play during the emergence of asexuality in haploid sexual systems. Moreover, we show that haploid purifying selection plays a key role in limiting the accumulation of deleterious alleles in Amazons, and we identify an autosomal locus associated with the Amazon phenotype. The sex chromosome, together with this autosomal locus, may underlie the switch to obligate asexuality in the Amazon populations.

Genome-Wide Identification and Molecular Evolutionary History of the Whirly Family Genes in Brassica napus

Whirly transcription factors are unique to plants, playing pivotal roles in managing leaf senescence and DNA repair. While present in various species, their identification in Brassica napus L. (B. napus) and their differences during hybridization and polyploidy has been elusive. Addressing this, our study delves into the functional and evolutionary aspects of the Whirly gene family during the emergence of B. napus, applying bioinformatics and comparative genomics. We identified six Whirly genes in B. napus. In Brassica rapa L. (B. rapa), three Whirly genes were identified, while four were found in Brassica oleracea L. (B. oleracea). The results show that the identified Whirly genes not only have homology but also share the same chromosomal positions. Phylogenetic analysis indicates that Whirly genes in monocots and dicots exhibit high conservation. In the evolutionary process, the Whirly gene family in B. napus experienced events of intron/exon loss. Collinearity insights point to intense purifying selection post-duplication. Promoter regions housed diverse cis-acting elements linked to photoresponse, anaerobic initiation, and methyl jasmonate responsiveness. Notably, elements tied to abscisic acid signaling and meristem expression were prominent in diploid ancestors but subdued in tetraploid B. napus. Tissue-specific expression unveiled analogous patterns within subfamily genes. Subsequent qRT-PCR analysis spotlighted BnAWHY1b's potential significance in abiotic stress response, particularly drought. These findings can be used as theoretical foundations to understand the functions and effects of the Whirly gene family in B. napus, providing references for the molecular mechanism of gene evolution between this species and its diploid ancestors.

Artificial intelligence-driven reverse vaccinology for Neisseria gonorrhoeae vaccine: Prioritizing epitope-based candidates

Neisseria gonorrhoeae is the causative agent of the sexually transmitted disease gonorrhea. The increasing prevalence of this disease worldwide, the rise of antibiotic-resistant strains, and the difficulties in treatment necessitate the development of a vaccine, highlighting the significance of preventative measures to control and eradicate the infection. Currently, there is no widely available vaccine, partly due to the bacterium's ability to evade natural immunity and the limited research investment in gonorrhea compared to other diseases. To identify distinct vaccine candidates, we chose to focus on the uncharacterized, hypothetical proteins (HPs) as our initial approach. Using the in silico method, we first carried out a comprehensive assessment of hypothetical proteins of Neisseria gonorrhoeae, encompassing assessments of physicochemical properties, cellular localization, secretary pathways, transmembrane regions, antigenicity, toxicity, and prediction of B-cell and T-cell epitopes, among other analyses. Detailed analysis of all HPs resulted in the functional annotation of twenty proteins with a great degree of confidence. Further, using the immuno-informatics approach, the prediction pipeline identified one CD8+ restricted T-cell epitope, seven linear B-cell epitopes, and seven conformational B-cell epitopes as putative epitope-based peptide vaccine candidates which certainly require further validation in laboratory settings. The study accentuates the promise of functional annotation and immuno-informatics in the systematic design of epitope-based peptide vaccines targeting Neisseria gonorrhoeae.

Draft genome sequence of Penicillium citrinum B9, a plant growth-promoting symbiont from barley rhizosphere

Penicillium citrinum strain B9 is a plant growth-promoting fungus isolated from Barley (Hordeum vulgare) rhizosphere. We report the first draft genome of P. citrinum B9 assembled using single-molecule real-time sequencing and Illumina reads. The assembled genome spans 31.3 Mb comprising nine contigs and 10,106 protein-encoding genes.

The Velvet transcription factor PnVeA regulates necrotrophic effectors and secondary metabolism in the wheat pathogen Parastagonospora nodorum

The fungus Parastagonospora nodorum causes septoria nodorum blotch on wheat. The role of the fungal Velvet-family transcription factor VeA in P. nodorum development and virulence was investigated here. Deletion of the P. nodorum VeA ortholog, PnVeA, resulted in growth abnormalities including pigmentation, abolished asexual sporulation and highly reduced virulence on wheat. Comparative RNA-Seq and RT-PCR analyses revealed that the deletion of PnVeA also decoupled the expression of major necrotrophic effector genes. In addition, the deletion of PnVeA resulted in an up-regulation of four predicted secondary metabolite (SM) gene clusters. Using liquid-chromatography mass-spectrometry, it was observed that one of the SM gene clusters led to an accumulation of the mycotoxin alternariol. PnVeA is essential for asexual sporulation, full virulence, secondary metabolism and necrotrophic effector regulation.

Uncovering Endolysins against Methicillin-Resistant Staphylococcus aureus Using a Microbial Single-Cell Genome Database

Endolysins, peptidoglycan hydrolases derived from bacteriophages (phages), are being developed as a promising alternative to conventional antibiotics. To obtain highly active endolysins, a diverse library of these endolysins is vital. We propose here microbial single-cell genome sequencing as an efficient tool to discover dozens of previously unknown endolysins, owing to its culture-independent sequencing method. As a proof of concept, we analyzed and recovered endolysin genes within prophage regions of Staphylococcus single-amplified genomes in human skin microbiome samples. We constructed a library of chimeric endolysins by shuffling domains of the natural endolysins and performed high-throughput screening against Staphylococcus aureus. One of the lead endolysins, bbst1027, exhibited desirable antimicrobial properties, such as rapid bactericidal activity, no detectable resistance development, and in vivo efficacy. We foresee that this endolysin discovery pipeline is in principle applicable to any bacterial target and boost the development of novel antimicrobial agents.

Identification and expression profiles of olfactory-related genes in the antennal transcriptome of Graphosoma rubrolineatum (Hemiptera: Pentatomidae)

Graphosoma rubrolineatum (Hemiptera: Pentatomidae) is an important pest of vegetables and herbs (e.g., Umbelliferae and Cruciferae) in China, Siberia, Korea, and Japan. Insects are highly dependent on their olfactory system to detect odorants. However, no molecular-mediated olfactory genes in G. rubrolineatum have yet been identified. In this study, we first established the antennal transcriptome of G. rubrolineatum and identified 189 candidate olfactory genes, including 31 odorant-binding proteins (OBPs), 15 chemosensory proteins (CSPs), four sensory neuron membrane proteins (SNMPs),94 odorant receptors (ORs), 23 ionotropic receptors (IRs), and 22 gustatory receptors (GRs). Additionally, phylogenetic trees were constructed for olfactory genes between G. rubrolineatum and other hemipteran insects. We also detected the expression profiles of ten OBPs, five CSPs, two SNMPs, five ORs, four IRs, and four GRs by real-time quantitative PCR. The results revealed that most genes (GrubOBP1/11/31, GrubCSP3/8, GrubSNMP1a/1b, GrubOrco/OR9/11/13, GrubGR1/4/22, GrubIR25/75h/76b/GluR1) were highly expressed in the antennae, GrubOBP13/31 and GrubCSP4/11/12 were highly expressed in the legs, while GrubOBP20 and GrubGR19 were highly expressed in the wings. Our results will enrich the gene inventory of G. rubrolineatum and provide further insight into the molecular chemosensory mechanisms of G. rubrolineatum.

Integrated proteomics, transcriptomics, and metabolomics offer novel insights into Cd resistance and accumulation in Poa pratensis

Kentucky bluegrass (Poa pratensis L., KB) demonstrates superior performance in both cadmium (Cd) accumulation and tolerance; however, the regulatory mechanisms and detoxification pathways in this species remain unclear. Therefore, phenotype, root ultrastructure, cell wall components, proteomics, transcriptomics, and metabolomics were analyzed under the hydroponic system to investigate the Cd tolerance and accumulation mechanisms in the Cd-tolerant KB variety 'Midnight (M)' and the Cd-sensitive variety 'Rugby II (R)' under Cd stress. The M variety exhibited higher levels of hydroxyl and carboxyl groups as revealed by Fourier transform infrared spectroscopy spectral analysis. Additionally, a reduced abundance of polysaccharide degradation proteins was observed in the M variety. The higher abundance of glutathione S-transferase and content of L-cysteine-glutathione disulfide and oxidized glutathione in the M variety may contribute to better performance of the M variety under Cd stress. Additionally, the R variety had an enhanced content of carboxylic acids and derivatives, increasing the Cd translocation capacity. Collectively, the down-regulation of cell wall polysaccharide degradation genes coupled with the up-regulation of glutathione metabolism genes enhances the tolerance to Cd stress in KB. Additionally, lignification of the endodermis and the increase in carboxylic acids and derivatives play crucial roles in the redistribution of Cd in KB.

Densovirus infection facilitates plant-virus transmission by an aphid

The interactions among plant viruses, insect vectors, and host plants have been well studied; however, the roles of insect viruses in this system have largely been neglected. We investigated the effects of MpnDV infection on aphid and PVY transmission using bioassays, RNA interference (RNAi), and GC-MS methods and green peach aphid (Myzus persicae (Sulzer)), potato virus Y (PVY), and densovirus (Myzus persicae nicotianae densovirus, MpnDV) as model systems. MpnDV increased the activities of its host, promoting population dispersal and leading to significant proliferation in tobacco plants by significantly enhancing the titer of the sesquiterpene (E)-β-farnesene (EβF) via up-regulation of expression levels of the MpFPPS1 gene. The proliferation and dispersal of MpnDV-positive individuals were faster than that of MpnDV-negative individuals in PVY-infected tobacco plants, which promoted the transmission of PVY. These results combined showed that an insect virus may facilitate the transmission of a plant virus by enhancing the locomotor activity and population proliferation of insect vectors. These findings provide novel opportunities for controlling insect vectors and plant viruses, which can be used in the development of novel management strategies.

Transcriptomic dataset of the development and maturation of the Rhipicephalus microplus ovary

To conduct differential gene expression analysis, ovaries from the cattle tick Rhipicephalus microplus were dissected at three distinct developmental stages (preingurgitated, immature ingurgitated, and mature ingurgitated). Additionally, undissected intact mature males and complete ingurgitated female ticks without ovaries (carcasses) were also collected to serve as reference samples for analysis. To perform total RNA purification, tissue from ten individuals representing each of the five previously described conditions was pooled. mRNA was isolated from the purified total RNA using the oligo (dT) method. Following fragmentation, double stranded cDNA was synthesized and ligated to sequencing adapters. Suitable-sized fragments were subsequently used for PCR amplification. Libraries were analyzed and quantified using an Agilent 2100 Bioanalyzer and an ABI StepOnePlus Real-Time PCR System. A total of 45.64 Gb bases were sequenced using the Illumina HiSeq sequencing platform. After assembling the samples and correcting for abundance, we obtained 82,877 unigenes. The total length, average length, N50, and GC content of the unigenes were 89,754,828 bp,1,082 bp,2,068 bp and 49.04 % respectively. For functional annotation, the unigenes were aligned with 7 functional databases. The number of unigenes identified in the functional databases were as follows: 32,518 (NR:39.24 %), 10,259 (NT:12.38 %), 23,624 (Swissprot:28.50 %), 22,203 (KOG:26.79 %), 25,072 (KEGG:30.25 %), 17,435(GO:21.04 %), and 23,220 (InterPro:28.02 %). Unigene candidate coding regions (CDS) among the unigenes were predicted using TransDecoder software and 42,143 CDS were detected. We also detected 10,522 simple sequence repeats (SSRs) distributed on 8,126 unigenes, and predicted 4,672 transcription factors (TF) coding unigenes. Our data can be used to identify genes that are important for male and female tick and arachnid reproduction and tick general physiology.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

Detection of Known and Novel Virus Sequences in the Black Solider Fly and Expression of Host Antiviral Pathways

The mass rearing of animals in close quarters can be highly conducive to microbe transmission, including pathogens. This has been shown multiple times in the case of important industrial insects such as crickets, silkworms, and honeybees. One industrial insect of increasing importance is the black soldier fly (Diptera: Hermetia illucens), as it can convert organic waste into high-quality protein and fatty acids. Along with this, they take up far less space than traditional protein sources, as millions of black soldier flies can be reared in a relatively small facility. Because of this, there is a growing interest in the pathogens that could impact black soldier fly-rearing efforts. So far, only three black soldier fly-associated viruses have been identified. We used metatranscriptomic sequencing to survey black soldier fly guts, frass, and diet for viruses. We detected sequences from two novel viruses. One, which we name Hermetia illucens sigma-like virus 1, is phylogenetically related to viruses of the genus Sigmavirus, which have been highly studied in Drosophila. The other novel virus, which we name Hermetia illucens inse-like virus 1, is the second double-stranded RNA virus of the order Ghabrivirales described in the black soldier fly, and groups within a new family of insect viruses called the Inseviridae. We also detected two black soldier fly-associated viruses previously identified by our group: BSF nairo-like virus 1 and BSF uncharacterized bunyavirus-like 1. Consistent with our previous study, these two viruses are found primarily in frass samples and occur together more often than expected at random. When analyzing host transcription, we found significant differences in gene expression for eight candidate antiviral genes in the black soldier fly when comparing samples with and without viral sequences. Our results suggest that black soldier fly-virus interactions are ongoing, and they could be of interest to black soldier fly producers.

Exploring oak processionary caterpillar induced lepidopterism (Part 1): unveiling molecular insights through transcriptomics and proteomics

Lepidopterism, a skin inflammation condition caused by direct or airborne exposure to irritating hairs (setae) from processionary caterpillars, is becoming a significant public health concern. Recent outbreaks of the oak processionary caterpillar (Thaumetopoea processionea) have caused noteworthy health and economic consequences, with a rising frequency expected in the future, exacerbated by global warming promoting the survival of the caterpillar. Current medical treatments focus on symptom relief due to the lack of an effective therapy. While the source is known, understanding the precise causes of symptoms remain incomplete understood. In this study, we employed an advanced method to extract venom from the setae and identify the venom components through high-quality de novo transcriptomics, venom proteomics, and bioinformatic analysis. A total of 171 venom components were identified, including allergens, odorant binding proteins, small peptides, enzymes, enzyme inhibitors, and chitin biosynthesis products, potentially responsible for inflammatory and allergic reactions. This work presents the first comprehensive proteotranscriptomic database of T. processionea, contributing to understanding the complexity of lepidopterism. Furthermore, these findings hold promise for advancing therapeutic approaches to mitigate the global health impact of T. processionea and related caterpillars.

Single-molecule imaging reveals the mechanism of bidirectional replication initiation in metazoa

Metazoan genomes are copied bidirectionally from thousands of replication origins. Replication initiation entails the assembly and activation of two CMG helicases (Cdc45⋅Mcm2-7⋅GINS) at each origin. This requires several replication firing factors (including TopBP1, RecQL4, and DONSON) whose exact roles are still under debate. How two helicases are correctly assembled and activated at each origin is a long-standing question. By visualizing the recruitment of GINS, Cdc45, TopBP1, RecQL4, and DONSON in real time, we uncovered that replication initiation is surprisingly dynamic. First, TopBP1 transiently binds to the origin and dissociates before the start of DNA synthesis. Second, two Cdc45 are recruited together, even though Cdc45 alone cannot dimerize. Next, two copies of DONSON and two GINS simultaneously arrive at the origin, completing the assembly of two CMG helicases. Finally, RecQL4 is recruited to the CMG⋅DONSON⋅DONSON⋅CMG complex and promotes DONSON dissociation and CMG activation via its ATPase activity.

The acetyltransferase SCO0988 controls positively specialized metabolism and morphological differentiation in the model strains Streptomyces coelicolor and Streptomyces lividans

Streptomycetes are well-known antibiotic producers possessing in their genomes numerous silent biosynthetic pathways that might direct the biosynthesis of novel bio-active specialized metabolites. It is thus of great interest to find ways to enhance the expression of these pathways to discover most needed novel antibiotics. In this study, we demonstrated that the over-expression of acetyltransferase SCO0988 up-regulated the production of specialized metabolites and accelerated sporulation of the weak antibiotic producer, Streptomyces lividans and that the deletion of this gene had opposite effects in the strong antibiotic producer, Streptomyces coelicolor. The comparative analysis of the acetylome of a S. lividans strain over-expressing sco0988 with that of the original strain revealed that SCO0988 acetylates a broad range of proteins of various pathways including BldKB/SCO5113, the extracellular solute-binding protein of an ABC-transporter involved in the up-take of a signal oligopeptide of the quorum sensing pathway. The up-take of this oligopeptide triggers the "bald cascade" that regulates positively specialized metabolism, aerial mycelium formation and sporulation in S. coelicolor. Interestingly, BldKB/SCO5113 was over-acetylated on four Lysine residues, including Lys425, upon SCO0988 over-expression. The bald phenotype of a bldKB mutant could be complemented by native bldKB but not by variant of bldKB in which the Lys425 was replaced by arginine, an amino acid that could not be acetylated or by glutamine, an amino acid that is expected to mimic acetylated lysine. Our study demonstrated that Lys425 was a critical residue for BldKB function but was inconclusive concerning the impact of acetylation of Lys425 on BldKB function.

Genetic factors acting prior to dormancy in sour cherry influence bloom time the following spring

Understanding the process of Prunus species floral development is crucial for developing strategies to manipulate bloom time and prevent crop loss due to climate change. Here, we present a detailed examination of flower development from initiation until bloom for early- and late-blooming sour cherries (Prunus cerasus) from a population segregating for a major bloom time QTL on chromosome 4. Using a new staging system, we show floral buds from early-blooming trees were persistently more advanced than those from late-blooming siblings. A genomic DNA coverage analysis revealed the late-blooming haplotype of this QTL, k, is located on a subgenome originating from the late-blooming P. fruticosa progenitor. Transcriptome analyses identified many genes within this QTL as differentially expressed between early- and late-blooming trees during the vegetative-to-floral transition. From these, we identified candidate genes for the late bloom phenotype, including multiple transcription factors homologous to Reproductive Meristem B3 domain-containing proteins. Additionally, we determined that the basis of k in sour cherry is likely separate from candidate genes found in sweet cherry-suggesting several major regulators of bloom time are located on Prunus chromosome 4.

Label-free-based proteomics analysis reveals differential proteins of sheep, goat and cow milk

Regarding the limited information on species protein differences between sheep, goat, and cow milk, the differentially expressed proteins in sheep, goat, and cow milk and their functional differences are analyzed using label-free proteomics technology to identify potential biomarkers. 770 proteins and 2914 peptide segments were identified. The statistical analysis showed significant differences in the relative abundances of the 74 proteins among the sheep, goat, and cow milk. CSN3 and LALBA can be used as potential biomarkers for goat milk, XDH can be used as potential biomarkers for cow milk, and CTSB and BPIFB1 can be used as potential biomarkers for sheep milk. The functional analysis using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes showed that these significantly different proteins were enriched by different pathways including thyroid hormone synthesis and glycerol phospholipid metabolism. The data revealed differences in the amounts and physiological functions of the milk proteins of different species, which may provide an important basis for research on the nutritional composition of dairy products and adulteration identification technology.

Saccharopolyspora mangrovi sp. nov., a novel mangrove soil actinobacterium with distinct metabolic potential revealed by comparative genomic analysis

A novel mangrove soil-derived actinomycete, strain S2-29T, was found to be most closely related to Saccharopolyspora karakumensis 5K548T based on 16 S rRNA sequence (99.24% similarity) and genomic phylogenetic analyses. However, significant divergence in digital DNA-DNA hybridization, average nucleotide identity, and unique biosynthetic gene cluster possession distinguished S2-29T as a distinct Saccharopolyspora species. Pan genome evaluation revealed exceptional genomic flexibility in genus Saccharopolyspora, with > 95% accessory genome content. Strain S2-29T harbored 718 unique genes, largely implicated in energetic metabolisms, indicating different metabolic capacities from its close relatives. Several uncharacterized biosynthetic gene clusters in strain S2-29T highlighted the strain's untapped capacity to produce novel functional compounds with potential biotechnological applications. Designation as novel species Saccharopolyspora mangrovi sp. nov. (type strain S2-29T = JCM 34,548T = CGMCC 4.7716T) was warranted, expanding the known Saccharopolyspora diversity and ecology. The discovery of this mangrove-adapted strain advances understanding of the genus while highlighting an untapped source of chemical diversity.

Tidings from the Tides-De novo transcriptome assembly of the endemic estuarine bivalve Villorita cyprinoides

The Indian black clam Villorita cyprinoides Gray, 1825, is an economically valuable estuarine bivalve that faces challenges from multiple stressors and anthropogenic pressures. However, limited genomic resources have hindered molecular investigations into the impact of these stressors on clam populations. Here, we have generated the first transcriptomic reference datasets for V. cyprinoides to address this knowledge gap. A total of 25,040,592 and 22,486,217 million Illumina paired-end reads generated from two individuals were assembled using Trinity and rnaSPAdes. From the 47,607 transcripts identified as Coding Domain Sequences, 37,487 returned positive BLAST hits against six different databases. Additionally, a total of 14,063 Single Sequence Repeats were identified using GMATA. This study significantly enhances the genetic understanding of V. cyprinoides, a potential candidate for aquaculture that supports the livelihoods of many people dependent on small-scale fisheries. The data generated provides insights into broader genealogical connections within the family Cyrenidae through comparative transcriptomics. Furthermore, this transcriptional profile serves as baseline data for future studies in toxicological and conservation genetics.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

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.

Comparative transcriptome database for Camellia sinensis reveals genes related to the cold sensitivity and albino mechanism of 'Anji Baicha'

Tea, a globally popular beverage, contains various beneficial secondary metabolites. Tea plants (Camellia sinensis) exhibit diverse genetic traits across cultivars, impacting yield, adaptability, morphology, and secondary metabolite composition. Many tea cultivars have been the subject of much research interest, which have led to the accumulation of publicly available RNA-seq data. As such, it has become possible to systematically summarize the characteristics of different cultivars at the transcriptomic level, identify functional genes, and infer gene functions through co-expression analysis. Here, the transcriptomes of 9 tea cultivars were assembled, and comparative analysis was conducted on the coding sequences of 13 cultivars. To give access to this data, we present TeaNekT (https://teanekt.sbs.ntu.edu.sg/), a web resource that facilitates the prediction of gene functions of various tea cultivars. We used TeaNekT to perform a cross-cultivar comparison of co-expressed gene clusters and tissue-specific gene expression. We observed that 'Anji Baicha' possesses the highest number of cultivar-specific genes and the second-highest number of expanded genes. These genes in 'Anji Baicha' tend to be enriched in functions associated with cold stress response, chloroplast thylakoid structure, and nitrogen metabolism. Notably, we identified three significantly expanded homologous genes in 'Anji Baicha' encoding the ICE1, SIZ1, and MAPKK2, which are closely associated with the cold sensitivity of 'Anji Baicha'. Additionally, one significantly expanded homologous gene in 'Anji Baicha' encoding regulatory factor RIQ may play a crucial role in the abnormal chloroplast structure and absence of thylakoid membranes in 'Anji Baicha'.

Genome-resolved metagenomics: a game changer for microbiome medicine

Recent substantial evidence implicating commensal bacteria in human diseases has given rise to a new domain in biomedical research: microbiome medicine. This emerging field aims to understand and leverage the human microbiota and derivative molecules for disease prevention and treatment. Despite the complex and hierarchical organization of this ecosystem, most research over the years has relied on 16S amplicon sequencing, a legacy of bacterial phylogeny and taxonomy. Although advanced sequencing technologies have enabled cost-effective analysis of entire microbiota, translating the relatively short nucleotide information into the functional and taxonomic organization of the microbiome has posed challenges until recently. In the last decade, genome-resolved metagenomics, which aims to reconstruct microbial genomes directly from whole-metagenome sequencing data, has made significant strides and continues to unveil the mysteries of various human-associated microbial communities. There has been a rapid increase in the volume of whole metagenome sequencing data and in the compilation of novel metagenome-assembled genomes and protein sequences in public depositories. This review provides an overview of the capabilities and methods of genome-resolved metagenomics for studying the human microbiome, with a focus on investigating the prokaryotic microbiota of the human gut. Just as decoding the human genome and its variations marked the beginning of the genomic medicine era, unraveling the genomes of commensal microbes and their sequence variations is ushering us into the era of microbiome medicine. Genome-resolved metagenomics stands as a pivotal tool in this transition and can accelerate our journey toward achieving these scientific and medical milestones.

Structure-aware deep learning model for peptide toxicity prediction

Antimicrobial resistance is a critical public health concern, necessitating the exploration of alternative treatments. While antimicrobial peptides (AMPs) show promise, assessing their toxicity using traditional wet lab methods is both time-consuming and costly. We introduce tAMPer, a novel multi-modal deep learning model designed to predict peptide toxicity by integrating the underlying amino acid sequence composition and the three-dimensional structure of peptides. tAMPer adopts a graph-based representation for peptides, encoding ColabFold-predicted structures, where nodes represent amino acids and edges represent spatial interactions. Structural features are extracted using graph neural networks, and recurrent neural networks capture sequential dependencies. tAMPer's performance was assessed on a publicly available protein toxicity benchmark and an AMP hemolysis data we generated. On the latter, tAMPer achieves an F1-score of 68.7%, outperforming the second-best method by 23.4%. On the protein benchmark, tAMPer exhibited an improvement of over 3.0% in the F1-score compared to current state-of-the-art methods. We anticipate tAMPer to accelerate AMP discovery and development by reducing the reliance on laborious toxicity screening experiments.

Clinical Genetic Testing for Hearing Loss: Implications for Genetic Counseling and Gene-Based Therapies

Genetic factors contribute significantly to congenital hearing loss, with non-syndromic cases being more prevalent and genetically heterogeneous. Currently, 150 genes have been associated with non-syndromic hearing loss, and their identification has improved our understanding of auditory physiology and potential therapeutic targets. Hearing loss gene panels offer comprehensive genetic testing for hereditary hearing loss, and advancements in sequencing technology have made genetic testing more accessible and affordable. Currently, genetic panel tests available at a relatively lower cost are offered to patients who face financial barriers. In this study, clinical and audiometric data were collected from six pediatric patients who underwent genetic panel testing. Known pathogenic variants in MYO15A, GJB2, and USH2A were most likely to be causal of hearing loss. Novel pathogenic variants in the MYO7A and TECTA genes were also identified. Variable hearing phenotypes and inheritance patterns were observed amongst individuals with different pathogenic variants. The identification of these variants contributes to the continually expanding knowledge base on genetic hearing loss and lays the groundwork for personalized treatment options in the future.

Temporal gene expression during asexual development of the apicomplexan Sarcocystis neurona

Asexual replication in the apicomplexan Sarcocystis neurona involves two main developmental stages: the motile extracellular merozoite and the sessile intracellular schizont. Merozoites invade host cells and transform into schizonts that undergo replication via endopolygeny to form multiple (64) daughter merozoites that are invasive to new host cells. Given that the capabilities of the merozoite vary significantly from the schizont, the patterns of transcript levels throughout the asexual lifecycle were determined and compared in this study. RNA-Seq data were generated from extracellular merozoites and four intracellular schizont development time points. Of the 6,938 genes annotated in the S. neurona genome, 6,784 were identified in the transcriptome. Of these, 4,111 genes exhibited significant differential expression between the merozoite and at least one schizont development time point. Transcript levels were significantly higher for 2,338 genes in the merozoite and 1,773 genes in the schizont stages. Included in this list were genes encoding the secretory pathogenesis determinants (SPDs), which encompass the surface antigen and SAG-related sequence (SAG/SRS) and the secretory organelle proteins of the invasive zoite stage (micronemes, rhoptries, and dense granules). As anticipated, many of the S. neurona SPD gene transcripts were abundant in merozoites. However, several SPD transcripts were elevated in intracellular schizonts, suggesting roles unrelated to host cell invasion and the initial establishment of the intracellular niche. The hypothetical genes that are potentially unique to the genus Sarcocystis are of particular interest. Their conserved expression patterns are instructive for future investigations into the possible functions of these putative Sarcocystis-unique genes.

IMPORTANCE

The genus Sarcocystis is an expansive clade within the Apicomplexa, with the species S. neurona being an important cause of neurological disease in horses. Research to decipher the biology of S. neurona and its host-pathogen interactions can be enhanced by gene expression data. This study has identified conserved apicomplexan orthologs in S. neurona, putative Sarcocystis-unique genes, and gene transcripts abundant in the merozoite and schizont stages. Importantly, we have identified distinct clusters of genes with transcript levels peaking during different intracellular schizont development time points, reflecting active gene expression changes across endopolygeny. Each cluster also has subsets of transcripts with unknown functions, and investigation of these seemingly Sarcocystis-unique transcripts will provide insights into the interesting biology of this parasite genus.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

Genome-Wide Identification of NDPK Family Genes and Expression Analysis under Abiotic Stress in Brassica napus

The NDPK gene family is an important group of genes in plants, playing a crucial role in regulating energy metabolism, growth, and differentiation, cell signal transduction, and response to abiotic stress. However, our understanding of the NDPK gene family in Brassica napus L. remains limited. This paper systematically analyzes the NDPK gene family in B. napus, particularly focusing on the evolutionary differences within the species. In this study, sixteen, nine, and eight NDPK genes were identified in B. napus and its diploid ancestors, respectively. These genes are not only homologous but also highly similar in their chromosomal locations. Phylogenetic analysis showed that the identified NDPK proteins were divided into four clades, each containing unique motif sequences, with most NDPKs experiencing a loss of introns/exons during evolution. Collinearity analysis revealed that the NDPK genes underwent whole-genome duplication (WGD) events, resulting in duplicate copies, and most of these duplicate genes were subjected to purifying selection. Cis-acting element analysis identified in the promoters of most NDPK genes elements related to a light response, methyl jasmonate response, and abscisic acid response, especially with an increased number of abscisic acid response elements in B. napus. RNA-Seq results indicated that NDPK genes in B. napus exhibited different expression patterns across various tissues. Further analysis through qRT-PCR revealed that BnNDPK genes responded significantly to stress conditions such as salt, drought, and methyl jasmonate. This study enhances our understanding of the NDPK gene family in B. napus, providing a preliminary theoretical basis for the functional study of NDPK genes and offering some references for further revealing the phenomenon of polyploidization in plants.

Comparative proteomics reveals that fatty acid metabolism is involved in myocardial adaptation to chronic hypoxic injury

Congenital heart disease (CHD) is the most serious form of heart disease, and chronic hypoxia is the basic physiological process underlying CHD. Some patients with CHD do not undergo surgery, and thus, they remain susceptible to chronic hypoxia, suggesting that some protective mechanism might exist in CHD patients. However, the mechanism underlying myocardial adaptation to chronic hypoxia remains unclear. Proteomics was used to identify the differentially expressed proteins in cardiomyocytes cultured under hypoxia for different durations. Western blotting assays were used to verify protein expression. A Real-Time Cell Analyzer (RTCA) was used to analyze cell growth. In this study, 3881 proteins were identified by proteomics. Subsequent bioinformatics analysis revealed that proteins were enriched in regulating oxidoreductase activity. Functional similarity cluster analyses showed that chronic hypoxia resulted in proteins enrichment in the mitochondrial metabolic pathway. Further KEGG analyses found that the proteins involved in fatty acid metabolism, the TCA cycle and oxidative phosphorylation were markedly upregulated. Moreover, knockdown of CPT1A or ECI1, which is critical for fatty acid degradation, suppressed the growth of cardiomyocytes under chronic hypoxia. The results of our study revealed that chronic hypoxia activates fatty acid metabolism to maintain the growth of cardiomyocytes.

Unveiling biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2

BACKGROUND

A growing number of studies have demonstrated that the polar regions have the potential to be a significant repository of microbial resources and a potential source of active ingredients. Genome mining strategy plays a key role in the discovery of bioactive secondary metabolites (SMs) from microorganisms. This work highlighted deciphering the biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2 by a combination of whole genome analysis and antiSMASH as well as feature-based molecular networking (MN) in the Global Natural Products Social Molecular Networking (GNPS).

RESULTS

In this study, a high-quality whole genome sequence of an Arctic marine strain MNP-2, with a size of 34.9 Mb was successfully obtained. Its total number of genes predicted by BRAKER software was 13,218, and that of non-coding RNAs (rRNA, sRNA, snRNA, and tRNA) predicted by using INFERNAL software was 204. AntiSMASH results indicated that strain MNP-2 harbors 56 biosynthetic gene clusters (BGCs), including 18 NRPS/NRPS-like gene clusters, 10 PKS/PKS-like gene clusters, 8 terpene synthse gene clusters, 5 indole synthase gene clusters, 10 hybrid gene clusters, and 5 fungal-RiPP gene clusters. Metabolic analyses of strain MNP-2 grown on various media using GNPS networking revealed its great potential for the biosynthesis of bioactive SMs containing a variety of heterocyclic and bridge-ring structures. For example, compound G-8 exhibited a potent anti-HIV effect with an IC50 value of 7.2 nM and an EC50 value of 0.9 nM. Compound G-6 had excellent in vitro cytotoxicities against the K562, MCF-7, Hela, DU145, U1975, SGC-7901, A549, MOLT-4, and HL60 cell lines, with IC50 values ranging from 0.10 to 3.3 µM, and showed significant anti-viral (H1N1 and H3N2) activities with IC50 values of 15.9 and 30.0 µM, respectively.

CONCLUSIONS

These findings definitely improve our knowledge about the molecular biology of genus A. sydowii and would effectively unveil the biosynthetic potential of strain MNP-2 using genomics and metabolomics techniques.

Identification of antennal alternative splicing by combining genome and full-length transcriptome analysis in Bactrocera dorsalis

Alternative splicing is an essential post-transcriptional regulatory mechanism that diversifies gene function by generating multiple protein isoforms from a single gene and act as a crucial role in insect environmental adaptation. Olfaction, a key sense for insect adaptation, relies heavily on the antennae, which are the primary olfactory organs expressing most of the olfactory genes. Despite the extensive annotation of olfactory genes within insect antennal tissues facilitated by high-throughput sequencing technology advancements, systematic analyses of alternative splicing are still relatively less. In this study, we focused on the oriental fruit fly (Bactrocera dorsalis), a significant pest of fruit crops. We performed a detailed analysis of alternative splicing in its antennae by utilizing the full-length transcriptome of its antennal tissue and the insect's genome. The results revealed 8600 non-redundant full-length transcripts identified in the oriental fruit fly antennal full-length transcriptome, spanning 4,145 gene loci. Over 40% of these loci exhibited multiple isoforms. Among these, 161 genes showed sex-biased isoform switching, involving seven different types of alternative splicing. Notably, events involving alternative transcription start sites (ATSS) and alternative transcription termination sites (ATTS) were the most common. Of all the genes undergoing ATSS and ATTS alternative splicing between male and female, 32 genes were alternatively spliced in protein coding regions, potentially affecting protein function. These genes were categorized based on the length of the sex-biased isoforms, with the highest difference in isoform fraction (dIF) associated with the ATSS type, including genes such as BdorABCA13, BdorCAT2, and BdorTSN3. Additionally, transcription factor binding sites for doublesex were identified upstream of both BdorABCA13 and BdorCAT2. Besides being expressed in the antennal tissues, BdorABCA13 and BdorCAT2 are also expressed in the mouthparts, legs, and genitalia of both female and male adults, suggesting their functional diversity. This study reveals alternative splicing events in the antennae of Bactrophora dorsalis from two aspects: odorant receptor genes and other types of genes expressed in the antennae. This study not only provides a research foundation for understanding the regulation of gene function by alternative splicing in the oriental fruit fly but also offers new insights for utilizing olfaction-based behavioral manipulation techniques to manage this pest.

Multi-omics Analysis of Volatile Flavor Components in Pacific Chub and Spanish Mackerel during Freezing using GC-MS-O

This study employed gas chromatography-mass spectrometry with olfactory (GC-MS-O) and multi-omics methods to investigate the changes in volatile flavor compounds during the freezing process of Pacific chub mackerel (Scomber japonicus) from Japan and China, and Spanish mackerel (Scomberomorus niphonius). A total of 18 volatile flavor compounds were identified, and significant differences in volatile flavor components were observed among samples frozen for 1 week, 1 year, and 2 years. The results of the Partial least squares regression (PLSR) indicated that the fishy odor was correlated with independent variables such as fatty acids (FA 22:4, FA 28:6, FA 24:4), differentially expressed genes (Gene.2425 (NDUFA5), Gene.38 (GPX1), and Gene.2844 (DAD1)). Classification and regression tree (CART) analysis revealed that the peak area values of fatty acids (FA 22:5, FA 20:4) and fatty acid esters of hydroxy fatty acids (FAHFA 18:0/22:3) were the main differentiating factors for fishy odor perception.

Chromosome-scale genome assembly of medicinal plant Tinospora sagittata (Oliv.) Gagnep. from the Menispermaceae family

Tinospora sagittata (Oliv.) Gagnep. is an important medicinal tetraploid plant in the Menispermaceae family. Its tuber, Radix Tinosporae, used in traditional Chinese medicine, is rich in diterpenoids and benzylisoquinoline alkaloids (BIAs). To enhance our understanding of medicinal compounds' biosynthesis and Menispermaceae's evolution, we herein report assembling a high-quality chromosome-scale genome with both PacBio HiFi and Illumina sequencing technologies. PacBio Sequel II generated 2.5 million circular consensus sequencing (CCS) reads, and a hybrid assembly strategy with Illumina sequencing resulted in 4483 contigs. The assembled genome size was 2.33 Gb, consisting of 4070 scaffolds (N50 = 42.06 Mb), of which 92.05% were assigned to 26 pseudochromosomes. T. sagittata's chromosomal-scale genome assembly, the first species in Menispermaceae, aids Menispermaceae evolution and T. sagittata's secondary metabolites biosynthesis understanding.

Isolation and Characterization of Phenylalanine Ammonia Lyase (PAL) Genes in Ferula pseudalliacea: Insights into the Phenylpropanoid Pathway

Phenylalanine ammonia lyase (PAL) is a key enzyme regulating the biosynthesis of the compounds of the phenylpropanoid pathway. This study aimed to isolate and characterize PAL genes from Ferula pseudalliacea Rech.f. (Apiales: Apiaceae) to better understand the regulation of metabolite production. Three PAL gene isoforms (FpPAL1-3) were identified and cloned using the 3'-RACE technique and confirmed by sequencing. Bioinformatics analysis revealed important structural features, such as phosphorylation sites, physicochemical properties, and evolutionary relationships. Expression analysis by qPCR demonstrated the differential transcription profiles of each FpPAL isoform across roots, stems, leaves, flowers, and seeds. FpPAL1 showed the highest expression in stems, FpPAL2 in roots and flowers, and FpPAL3 in flowers. The presence of three isoforms of PAL in F. pseudalliacea, along with the diversity of PAL genes and their tissue-specific expression profiles, suggests that complex modes of regulation exist for phenylpropanoid biosynthesis in this important medicinal plant. The predicted interaction network revealed associations with key metabolic pathways, emphasizing the multifaceted roles of these PAL genes. In silico biochemical analyses revealed the hydrophilicity of the FpPAL isozyme; however, further analysis of substrate specificity and enzyme kinetics can clarify the specific role of each FpPAL isozyme. These comprehensive results increase the understanding of PAL genes in F. pseudalliacea, helping to characterize their contributions to secondary metabolite biosynthesis.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Structural and Biochemical Characterization of a New Phage-Encoded Muramidase, KTN6 Gp46

Background

Endolysins are phage-encoded lytic enzymes that degrade bacterial peptidoglycan at the end of phage lytic cycles to release new phage particles. These enzymes are being explored as an alternative to small-molecule antibiotics.

Methods

The crystal structure of KTN6 Gp46 was determined and compared with a ColabFold model. Cleavage specificity was examined using a peptidoglycan digest and reversed-phase high-performance liquid chromatography coupled to mass spectrometry (HPLC/MS).

Results

The structure of KTN6 Gp46 could be determined at 1.4 Å resolution, and key differences in loops of the putative peptidoglycan binding domain were identified in comparison with its closest known homologue, the endolysin of phage SPN1S. Reversed-phase HPLC/MS analysis of the reaction products following peptidoglycan digestion confirmed the muramidase activity of Gp46, consistent with structural predictions.

Conclusion

These insights into the structure and function of endolysins further expand the toolbox for endolysin engineering and explore their potential in enzyme-based antibacterial design strategies.

A multi-tissue de novo transcriptome assembly and relative gene expression of the vulnerable freshwater salmonid Thymallus ligericus

Freshwater ecosystems are among the most endangered ecosystems worldwide. While numerous taxa are on the verge of extinction as a result of global changes and direct or indirect anthropogenic activity, genomic and transcriptomic resources represent a key tool for comprehending species' adaptability and serve as the foundation for conservation initiatives. The Loire grayling, Thymallus ligericus, is a freshwater European salmonid endemic to the upper Loire River basin. The species is comprised of fragmented populations that are dispersed over a small area and it has been identified as a vulnerable species. Here, we provide a multi-tissue de novo transcriptome assembly of T. ligericus. The completeness and integrity of the transcriptome were assessed before and after redundancy removal with lineage-specific libraries from Eukaryota, Metazoa, Vertebrata, and Actinopterygii. Relative gene expression was assessed for each of the analyzed tissues, using the de novo assembled transcriptome and a genome-based analysis using the available T. thymallus genome as a reference. The final assembly, with a contig N50 of 1221 and Benchmarking Universal Single-Copy Orthologs (BUSCO) scores above 94%, is made accessible along with structural and functional annotations and relative gene expression of the five tissues (NCBI SRA and FigShare databases). This is the first transcriptomic resource for this species, which provides a foundation for future research on this and other salmonid species that are increasingly exposed to environmental stressors.

Haplotype-resolved chromosomal-level genome assembly reveals regulatory variations in mulberry fruit anthocyanin content

Understanding the intricate regulatory mechanisms underlying the anthocyanin content (AC) in fruits and vegetables is crucial for advanced biotechnological customization. In this study, we generated high-quality haplotype-resolved genome assemblies for two mulberry cultivars: the high-AC 'Zhongsang5801' (ZS5801) and the low-AC 'Zhenzhubai' (ZZB). Additionally, we conducted a comprehensive analysis of genes associated with AC production. Through genome-wide association studies (GWAS) on 112 mulberry fruits, we identified MaVHAG3, which encodes a vacuolar-type H+-ATPase G3 subunit, as a key gene linked to purple pigmentation. To gain deeper insights into the genetic and molecular processes underlying high AC, we compared the genomes of ZS5801 and ZZB, along with fruit transcriptome data across five developmental stages, and quantified the accumulation of metabolic substances. Compared to ZZB, ZS5801 exhibited significantly more differentially expressed genes (DEGs) related to anthocyanin metabolism and higher levels of anthocyanins and flavonoids. Comparative analyses revealed expansions and contractions in the flavonol synthase (FLS) and dihydroflavonol 4-reductase (DFR) genes, resulting in altered carbon flow. Co-expression analysis demonstrated that ZS5801 displayed more significant alterations in genes involved in late-stage AC regulation compared to ZZB, particularly during the phase stage. In summary, our findings provide valuable insights into the regulation of mulberry fruit AC, offering genetic resources to enhance cultivars with higher AC traits.

Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Detection of serum major histocompatibility complex I (HLA-1) and β2-microglobulin (β2M) in pre-eclampsia using isobaric tags for relative and absolute quantitation (iTRAQ)

OBJECTIVE

The purpose of this preliminary investigation into the pathogenesis of pre-eclampsia was to screen the differential proteins in the serum of pregnant women with normal pregnancy and early-onset pre-eclampsia using isobaric tags for relative and absolute quantitation (iTRAQ), so as to identify serum biomarkers for the early diagnosis of pre-eclampsia.

METHODS

We examined the peripheral serum of 58 normal pregnant women and 42 pregnant women with early-onset pre-eclampsia using iTRAQ; the differentially expressed proteins were screened for bioinformatics analysis; and the expression of candidate proteins human leukocyte antigen-1 (HLA-1) and β2-microglobulin (β2M) in placental tissues was detected using western blot.

RESULTS

We identified a total of 63 differential proteins in the serum of patients from the normal control group and the pre-eclampsia group, and this included 24 up-regulated proteins and 39 down-regulated proteins. The western blot results of placental tissue showed reduced HLA-1 expression (1.12 ± 0.23) in the placenta in the pre-eclampsia group as compared with the normal control group (1.34 ± 0.22). Consistent with the results observed in the serum, β2M in the placenta in the pre-eclampsia group was significantly elevated (1.05 ± 0.47) in comparison with the normal group (0.75 ± 0.33) (P < 0.05).

CONCLUSION

In this study, we found that iTRAQ technology was useful for identifying differentially expressed proteins in the peripheral serum of pregnant women with pre-eclampsia, and that HLA-1 and β2M, which may be involved in the occurrence of pre-eclampsia, show promise as predictive markers of pre-eclampsia.

Decoding the Arsenal: Protist Effectors and Their Impact on Photosynthetic Hosts

Interactions between various microbial pathogens including viruses, bacteria, fungi, oomycetes, and their plant hosts have traditionally been the focus of phytopathology. In recent years, a significant and growing interest in the study of eukaryotic microorganisms not classified among fungi or oomycetes has emerged. Many of these protists establish complex interactions with photosynthetic hosts, and understanding these interactions is crucial in understanding the dynamics of these parasites within traditional and emerging types of farming, including marine aquaculture. Many phytopathogenic protists are biotrophs with complex polyphasic life cycles, which makes them difficult or impossible to culture, a fact reflected in a wide gap in the availability of comprehensive genomic data when compared to fungal and oomycete plant pathogens. Furthermore, our ability to use available genomic resources for these protists is limited by the broad taxonomic distance that these organisms span, which makes comparisons with other genomic datasets difficult. The current rapid progress in genomics and computational tools for the prediction of protein functions and interactions is revolutionizing the landscape in plant pathology. This is also opening novel possibilities, specifically for a deeper understanding of protist effectors. Tools like AlphaFold2 enable structure-based function prediction of effector candidates with divergent protein sequences. In turn, this allows us to ask better biological questions and, coupled with innovative experimental strategies, will lead into a new era of effector research, especially for protists, to expand our knowledge on these elusive pathogens and their interactions with photosynthetic hosts. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.