Search

Iris cangshanensis X. F Jiang, Y. Yang, J. F Liu & Y. P Xie, 2024, is a newly described perennial herb with a small population size. This study utilized Illumina sequencing technology to construct the complete chloroplast genome of Iris cangshanensis. The chloroplast genome is 154,008 base pairs (bp) in size with 37.32% overall GC content, consisting of an 83,559 bp (35.35%) large single-copy region (LSC), a 17,937 bp (30.56%) small single-copy region (SSC), and a pair of 26,256 bp (42.78%) inverted repeat (IR) sequences. The chloroplast genome has 130 genes, including 84 protein-coding genes, 38 tRNA genes, and eight rRNA genes. A maximum-likelihood (ML) phylogenetic analysis based on the plastomes of 15 Iridaceae species revealed that Iris cangshanensis is closely related to the other 12 species within the Iris genus. Owing to the limited chloroplast genome sampling of the genus Iris in the present study, determining the position of Iris cangshanensis in this genus is difficult. Overall, the chloroplast genome reported in this study will be useful for further taxonomic and evolutionary studies of Iris. This study constructs the complete chloroplast genome of Iris cangshanensis. The chloroplast genome is 154,008 base pairs (bp) in size with 37.32% overall GC content, consisting of an 83,559 bp (35.35%) large single-copy region (LSC), a 17,937 bp (30.56%) small single-copy region (SSC), and a pair of 26,256 bp (42.78%) inverted repeat (IR) sequences. The chloroplast genome has 130 genes, including 84 protein-coding genes, 38 tRNA genes, and eight rRNA genes. A maximum-likelihood (ML) phylogenetic analysis based on 15 Iridaceae species revealed that Iris cangshanensis is closely related to the species within the Iris genus.

The Leuciscidae (minnows, shiners and relatives) is a diverse family of freshwater fishes with many species endangered due to anthropogenic stressors. Notropis oxyrhynchus and Notropis buccula are two shiners found only in the upper Brazos River basin in Texas, USA and listed as endangered due to contracted habitat. The complete mitochondrial genome was sequenced for two vouchered specimens for each species; Notropis oxyrhynchus having a total mitogenome length of 16,711 bp and N. buccula having a total mitogenome length 16685–16686 bp, with both including 13 protein-coding genes, 22 transfer RNAs genes, and 2 ribosomal RNA genes. Phylogenetic analysis supports previous hypotheses regarding placement of these species.

In present study, the mitogenome of Rosa laevigata was assembled and characterized, with a total length of 281,693 bp and a GC content of 45.4%. The genome encompasses 53 genes, including 32 protein-coding genes, 3 rRNA genes and 18 tRNA genes. Furthermore, 19 MTPTs were identified, ranging from 48 to 1585 bp, covering 3.1% (8764 bp) of the mitogenome. Phylogenetic analysis of 34 Rosaceae species based on 21 common conserved protein-coding genes detected the monophyly of Rosa, with R. laevigata and R. chinensis forming a sister clade to R. rugosa. The mitogenome provides valuable genetic resources for R. laevigata utilization and for further phylogeny reconstruction of Rosa and Rosaceae.

Agata Rocka,1 Maria Suchcicka,2 Aleksandra M Jankowska,2 Magdalena M Woźniak,1 Monika Lejman3 1Pediatric Radiology, Medical University of Lublin, Medical University of Lublin, Prof. Antoni Gębali 6, Lublin, 20-093, Poland; 2Student Scientific Society of Independent Laboratory of Genetic Diagnostics, Medical University of Lublin, Prof. Antoni Gębali 6, Lublin, 20-093, Poland; 3Independent Laboratory of Genetic Diagnostics, Medical University of Lublin, Prof. Antoni Gębali 6, Lublin, 20-093, PolandCorrespondence: Agata Rocka, Pediatric Radiology, Medical University of Lublin, Medical University of Lublin, Prof. Antoni Gębali 6, Lublin, 20-093, Poland, Email agatarocka2@gmail.comAbstract: The aim of this study is to analyze available research on targeting signaling pathways for the development of new drugs in patients with T-cell acute lymphoblastic leukemia (T-ALL). This analysis focuses specifically on the role of LCK tyrosine kinase and mTOR signaling pathways in pediatric patients. Outcome: Current literature suggests that these pathways play a significant role in the regulation of T-cell cycles, making them potential therapeutic targets. However, despite promising findings, there remains a need for further research, particularly in pediatric populations, to fully understand the therapeutic implications and to optimize drug development. The conclusion drawn from this analysis highlights the significant influence of LCK and mTOR on T-cell cycle regulation, underscoring the importance of continued investigation in this area.Keywords: T-ALL, pediatric, signaling pathways, targeted therapy

Structural variations (SVs), a newly discovered genetic variation, have gained increasing recognition for their importance, yet much about them remains unknown. With the completion of whole-genome sequencing projects in oil crops, more SVs have been identified, revealing their types, genomic distribution, and characteristics. These findings have demonstrated the crucial roles of SVs in regulating gene expression, driving trait innovation, facilitating domestication, making this an opportune time for a systematic review. We summarized the progress of SV-related studies in oil crops, focusing on the types of SVs and their mechanisms of occurrence, the strategies and methods for SV detection, and the SVs identified in oil crops such as rapeseed, soybean, peanut, and sesame. The various types of SVs, such as presence-absence variations (PAVs), copy number variations (CNVs), and homeologous exchanges (HEs), have been shown. Along with their genomic characterization, their roles in crop domestication and breeding, and regulatory impact on gene expression and agronomic traits have also been demonstrated. This review will provide an overview of the SV research process in oil crops, enabling researchers to quickly understand key information and apply this knowledge in future studies and crop breeding.

In this study, we present the complete mitochondrial genome of Rhinogobius shennongensis, sequenced using the Illumina HiSeq platform. The result is a circular mitogenome of 16,500 base pairs (bp) with a pronounced A + T nucleotide bias. Phylogenetic analysis was conducted using both Bayesian inference (BI) and maximum-likelihood (ML) methods, utilizing the complete mitochondrial genomes of R. shennongensis and 19 additional Rhinogobius species. The analysis indicated a close mitochondrial relationship between R. maculagenys and R. shennongensis. This research not only enriches genetic resources but also supports phylogenetic studies within the genus Rhinogobius. To further clarify evolutionary dynamics, additional mitochondrial genome sequencing of other Rhinogobius species is essential.

Abstract Background Genotype imputation is a cost-effective method for obtaining sequence genotypes for downstream analyses such as genome-wide association studies (GWAS). However, low imputation accuracy can increase the risk of false positives, so it is important to pre-filter data or at least assess the potential limitations due to imputation accuracy. In this study, we benchmarked three different imputation programs (Beagle 5.2, Minimac4 and IMPUTE5) and compared the empirical accuracy of imputation with the software estimated accuracy of imputation (Rsqsoft). We also tested the accuracy of imputation in cattle for autosomal and X chromosomes, SNP and INDEL, when imputing from either low-density or high-density genotypes. Results The accuracy of imputing sequence variants from real high-density genotypes was higher than from low-density genotypes. In our software benchmark, all programs performed well with only minor differences in accuracy. While there was a close relationship between empirical imputation accuracy and the imputation Rsqsoft, this differed considerably for Minimac4 compared to Beagle 5.2 and IMPUTE5. We found that the Rsqsoft threshold for removing poorly imputed variants must be customised according to the software and this should be accounted for when merging data from multiple studies, such as in meta-GWAS studies. We also found that imposing an Rsqsoft filter has a positive impact on genomic regions with poor imputation accuracy due to large segmental duplications that are susceptible to error-prone alignment. Overall, our results showed that on average the imputation accuracy for INDEL was approximately 6% lower than SNP for all software programs. Importantly, the imputation accuracy for the non-PAR (non-Pseudo-Autosomal Region) of the X chromosome was comparable to autosomal imputation accuracy, while for the PAR it was substantially lower, particularly when starting from low-density genotypes. Conclusions This study provides an empirically derived approach to apply customised software-specific Rsqsoft thresholds for downstream analyses of imputed variants, such as needed for a meta-GWAS. The very poor empirical imputation accuracy for variants on the PAR when starting from low density genotypes demonstrates that this region should be imputed starting from a higher density of real genotypes.

Abstract The Human Phenotype Ontology (HPO) is widely used for annotating clinical text data, and sufficient annotation is crucial for the effective utilization of clinical texts. It was known that the use of LLMs can successfully extract symptoms and findings, but cannot annotate them with the HPO. We hypothesized that one of the potential issue for this is the lack of appropriate terms in the HPO. Therefore, during the Biomedical Linked Annotation Hackathon 8 (BLAH8), we attempted the following two tasks in order to grasp the overall picture of HPO. (1) Extract all HPO terms for each of the 23 HPO subclasses (defined as categories) directly under the HPO "Phenotypic abnormality" and then (2) search for major attributes in each of 23 categories. We employed LLM for these two tasks related to examining HPO and, at the same time, found that LLM didn't work well without ingenuity for tasks that lacked sentences and context. A manual search for terms within each category revealed that the HPO contains a mix of terms with four major attributes: (1) Disease Name, (2) Condition, (3) Test Data, and (4) Symptoms and Findings. Manual curation showed that the ratio of symptoms and findings varied from 0 to 93.1% across categories. For clinicians, who are end-users of medical terminology including HPO, it is difficult to understand ontologies. However, for good quality ontology is also important for good-quality data, and a clinician’s help is essential. It is also important to make the overall picture and limitations of ontologies easy to understand in order to bring out the explanatory power of LLMs and artificial intelligence.

Abstract Background Compatibility between plant parasites and their hosts is genetically determined {Citation}both interacting organisms. For example, plants may carry resistance (R) genes or deploy chemical defences. Aphid saliva contains many proteins that are secreted into host tissues. Subsets of these proteins are predicted to act as effectors, either subverting or triggering host immunity. However, associating particular effectors with virulence or avirulence outcomes presents challenges due to the combinatorial complexity. Here we use defined aphid and host genetics to test for co-segregation of expressed aphid transcripts and proteins with virulent or avirulent phenotypes. Results We compared virulent and avirulent pea aphid parental genotypes, and their bulk segregant F1 progeny on Medicago truncatula genotypes carrying or lacking the RAP1 (Resistance to Acyrthosiphon pisum 1) resistance quantitative trait locus. Differential gene expression analysis of whole body and head samples, in combination with proteomics of saliva and salivary glands, enabled us to pinpoint proteins associated with virulence or avirulence phenotypes. There was relatively little impact of host genotype, whereas large numbers of transcripts and proteins were differentially expressed between parental aphids, likely a reflection of their classification as divergent biotypes within the pea aphid species complex. Many fewer transcripts intersected with the equivalent differential expression patterns in the bulked F1 progeny, providing an effective filter for removing genomic background effects. Overall, there were more upregulated genes detected in the F1 avirulent dataset compared with the virulent one. Some genes were differentially expressed both in the transcriptome and in the proteome datasets, with aminopeptidase N proteins being the most frequent differentially expressed family. In addition, a substantial proportion (27%) of salivary proteins lack annotations, suggesting that many novel functions remain to be discovered. Conclusions Especially when combined with tightly controlled genetics of both insect and host plant, multi-omics approaches are powerful tools for revealing and filtering candidate lists down to plausible genes for further functional analysis as putative aphid effectors.

Abstract Background Nematodes of the genus Steinernema and their Xenorhabdus bacterial symbionts are lethal entomopathogens that are useful in the biocontrol of insect pests, as sources of diverse natural products, and as research models for mutualism and parasitism. Xenorhabdus play a central role in all aspects of the Steinernema lifecycle, and a deeper understanding of their genomes therefore has the potential to spur advances in each of these applications. Results Here, we report a comparative genomics analysis of Xenorhabdus griffiniae, including the symbiont of Steinernema hermaphroditum nematodes, for which genetic and genomic tools are being developed. We sequenced and assembled circularized genomes for three Xenorhabdus strains: HGB2511, ID10 and TH1. We then determined their relationships to other Xenorhabdus and delineated their species via phylogenomic analyses, concluding that HGB2511 and ID10 are Xenorhabdus griffiniae while TH1 is a novel species. These additions to the existing X. griffiniae landscape further allowed for the identification of two subspecies within the clade. Consistent with other Xenorhabdus, the analysed X. griffiniae genomes each encode a wide array of antimicrobials and virulence-related proteins. Comparative genomic analyses, including the creation of a pangenome, revealed that a large amount of the intraspecies variation in X. griffiniae is contained within the mobilome and attributable to prophage loci. In addition, CRISPR arrays, secondary metabolite potential and toxin genes all varied among strains within the X. griffiniae species. Conclusions Our findings suggest that phage-related genes drive the genomic diversity in closely related Xenorhabdus symbionts, and that these may underlie some of the traits most associated with the lifestyle and survival of entomopathogenic nematodes and their bacteria: virulence and competition. This study establishes a broad knowledge base for further exploration of not only the relationships between X. griffiniae species and their nematode hosts but also the molecular mechanisms that underlie their entomopathogenic lifestyle.

Abstract Background Metabarcoding targeting the internal transcribed spacer (ITS) region is commonly used to characterize fungal communities of various environments. Given their size and complexity, raw ITS sequences are necessarily processed and quality-filtered with bioinformatic pipelines. However, such pipelines are not yet standardized, especially for fungal communities, and those available may produce contrasting results. While some pipelines cluster sequences based on a specified percentage of base pair similarity into operational taxonomic units (OTUs), others utilize denoising techniques to infer amplicon sequencing variants (ASVs). While ASVs are now considered a more accurate representation of taxonomic diversity for prokaryote communities based on 16S rRNA amplicon sequencing, the applicability of this method for fungal ITS sequences is still debated. Results Here we compared the performance of two commonly used pipelines DADA2 (inferring ASVs) and mothur (clustering OTUs) on fungal metabarcoding sequences originating from two different environmental sample types (fresh bovine feces and pasture soil). At a 99% OTU similarity threshold, mothur consistently identified a higher fungal richness compared to DADA2. In addition, mothur generated homogenous relative abundances across multiple technical replicates (n = 18), while DADA2 results for the same replicates were highly heterogeneous. Conclusions Our study highlights a potential pipeline-associated bias in fungal metabarcoding data analysis of environmental samples. Based on the homogeneity of relative abundances across replicates and the capacity to detect OTUs/ASVs, we suggest using OTU clustering with a similarity of 97% as the most appropriate option for processing fungal metabarcoding data.

Abstract The over 4,100 species of bird lice are classified into 214 genera in the parvorders Amblycera and Ischnocera. Congeneric species of bird lice usually share much similarity in morphology and in mitochondrial (mt) genome organization. Two recent studies, however, reported substantial intra-genus variation in mt genome organization in bird lice. Both the ancestral single-chromosome mt genome and a fragmented mt genome with two or three minichromosomes were observed in the genera Austromenopon and Laemobothrion. To better understand intra-genus variation in mt genome organization, we sequenced the complete mt genome of the white spoonbill louse Ibidoecus plataleae and compared it with that of the glossy ibis feather louse Ibidoecus bisignatus reported previously. We found that I. plataleae had a fragmented mt genome with 12 minichromosomes; each minichromosome was 2,798 to 3,628 bp in size and had 2 to 6 genes. This is in stark contrast to the mt genome of I. bisignatus, which has all genes on a single chromosome, 14,909 bp in size. This is the most drastic intra-genus variation in mt genome organization observed to date in animals, indicating an unprecedented rapid process of mt genome fragmentation in the genus Ibidoecus. The divergence time between I. plataleae and I. bisignatus is currently unknown but is estimated to be less than 23 million years. Either many minichromosal split events occurred after I. plataleae diverged from I. bisignatus, or one minichromosome splits into multiple minichromosomes in a single event. Sequencing and comparing more Ibidoecusi species will help understand the unusual mt genome fragmentation in this genus.

Abstract Background High-quality goatskins are valuable byproducts usually produced by indigenous goat breeds with poorer production performance in Asia and Africa. However, the genetic and molecular mechanisms underpinning goatskin’s biomechanical properties (e.g., tensile strength) remain elusive. Mechanistic exploration of these traits could greatly aid the genetic improvement and genetic resource conservation of native breeds in these regions. To fulfill this purpose, we collected skin tissues from three goat breeds: Huai goat (HG), a Chinese native variety producing high-quality goatskins; Yudong meat goat (YDMG), a crossbreed of HG and Boer goat; Henan dairy goat (HNDG), a dairy goat breed. Results Scanning electronic microscopy analysis of skin tissues found that the collagen fiber diameters, collagen fibril diameters, and crimps significantly differed among the three goat breeds; however, collagen fibril diameters are similar in HG and HNDG. A sum of 230, 775, and 86 differentially expressed genes (DEGs) were identified from YDMG versus HNDG, HG versus HNDG, and YDMG versus HG, respectively. Functional enrichment analysis suggested that signaling pathways involved in fatty acid, retinol, steroid metabolisms, and GO items related to the physical properties of the skin (e.g., collagen-containing extracellular matrix) are significantly overrepresented in DEGs identified from meat versus dairy goats. Furthermore, 106 DEGs (e.g., COL1A1, COL1A2, and SPARC) showed specific expression patterns in HG and YDMG versus HNDG. Items about biophysical features of skin (e.g., extracellular matrix organization and ECM proteoglycans) are markedly enriched. Protein-protein interaction analysis suggested that two growth factors (IGF1 and PDGFD) are latent collagen and other ECM protein expression modulators. Conclusion Ultrastructural analysis of goat skin tissues suggested that collagen fibril diameter is not a major factor affecting goatskin quality. Transcriptomic profiles unveiled core genes and associated biological processes potentially involved in regulating goatskin quality. These discoveries shined new light on deeper understanding the mechanisms of hide-related traits in goat and other livestock.

Abstract Background The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) persists, giving rise to new variants characterized by mutations in the spike protein. However, public data regarding the virus's evolutionary trend is not widely available after the downgrade of coronavirus disease 2019(COVID-19). Therefore, this study aimed to investigate the applicability of an in-house Sanger-based method for identifying SARS-CoV-2 variants, particularly focusing on newly emerged Omicron variants, and updating the epidemiology of COVID-19 during the 8th wave in Hiroshima Prefecture. Results A total of 639 saliva samples of individuals who had tested positive for COVID-19, received from Hiroshima City Medical Association Clinical Laboratory Center between February 01, 2023, and March 12, 2024, were included in the study. SARS-CoV-2 variants were identified in 69.3% (443/639) with the mean viral titer 2 × 106 copies/mL, and high viral titer in Omicron variant XBC.1.6* (5 × 108 copies/mL) using RT-qPCR. By partial Spike gene-based sequencing using the Sanger Sequencing strategy, Omicron sub-lineages XXB.1, BA.5, and EG.1 were identified during different periods. A comprehensive phylogenetic analysis of 7383 SARS-CoV-2 strains retrieved from GISAID, collected in Hiroshima from the onset of the COVID-19 pandemic in early 2020 until July 2024, revealed the dynamic evolution of SARS-CoV-2 variants over time. The study found a similar pattern of variant distribution between the full genomes from GISAID, and the partial genomes obtained from our screening strategy during the same period. Conclusions Our study revealed that all SARS-CoV-2 viruses circulated in Hiroshima were Omicron variants and their sub-lineages during the 8th wave outbreak in Hiroshima. Persistent molecular surveillance of SARS-CoV-2 is needed for the decision-making and strategic planning of the public promptly. Our study added evidence for the usefulness of SARS-CoV-2 spike gene partial sequencing-based SARS-CoV-2 variant identification strategy for mass screening and molecular surveillance even though the evolution of newly emerged various SARS-CoV-2 Omicron variants.

Abstract Background Tea is an important cash crop that significantly contributes to rural development, poverty reduction and food security in many developing countries. It provides livelihoods for millions of smallholder producers and aids their economic stability. Anthocyanins in tea leaves provides excellent commercial quality and germplasm exploration potential. These compounds give tea leaves vibrant colors and increase health benefits. The current understanding of the synergistic regulation mechanisms responsible for color changes in purple tea, attributed to anthocyanin degradation, remains unclear. Results In this study, we have identified 30 gene families within the genome that are associated to with anthocyanin metabolism from tea. These gene families play distinct roles in the biosynthesis of anthocyanin including the formation of the core, structure, modification of the molecular framework, facilitation of transport process, regulation of gene expression, breakdown pathways, sugar transportation and iron ion respectively. Subsequently, we investigated the synergistic mechanisms of anthocyanin metabolism related gene families within tea leaves using structural equation modeling. The results showed that sugar transport positively affects anthocyanin transportation, and promotes anthocyanin degradation during leaf pigmentation, whereas, it inhibits anthocyanin degradation during the fading of leaf color. Further, Iron ions facilitate the degradation of anthocyanins during their deposition and conversely, impede this degradation process during digestion. These finding suggests that tea plants may regulate the synthesis and degradation of anthocyanins through sugar transport and iron ions ensure healthy levels and vibrant colors. Conclusions Our study contributes valuable information into the dynamic equilibrium anthocyanin mechanism and sheds light on complex regulatory mechanisms that govern the synthesis, transport and degradation of these pigments. These insights could be further used to develop strategies for enhancing anthocyanins content in unique tea germplasm to aid tea industry in producing new tea products with increased health benefits and aesthetic appeals.

Abstract Background Pepper (Capsicum pubescens), one of five domesticated pepper species, has unique characteristics, such as numerous hairs on the epidermis of its leaves and stems, black seeds, and vibrant purple flowers. To date, no studies have reported on the complete assembly of the mitochondrial genome (mitogenome) of C. pubescens. Understanding the mitogenome is crucial for further research on C. pubescens. Results In our study, we successfully assembled the first mitogenome of C. pubescens, which was assigned the GenBank accession number OP957066. This mitogenome has a length of 454,165 bp and exhibits the typical circular structure observed in most mitogenomes. We annotated a total of 70 genes, including 35 protein-coding genes (PCGs), 30 tRNA genes, 3 rRNA genes, and 2 pseudogenes. Compared to the other three pepper mitogenomes (KJ865409, KJ865410, and MN196478), C. pubescens OP957066 exhibited four unique PCGs (atp4, atp8, mttB, and rps1), while two PCGs (rpl10 and rps3) were absent. Notably, each of the three pepper mitogenomes from C. annuum (KJ865409, KJ865410, and MN196478) experienced the loss of four PCGs (atp4, atp8, mttB, and rps1). To further explore the evolutionary relationships, we reconstructed a phylogenetic tree using the mitogenomes of C. pubescens and fourteen other species. Structural comparison and synteny analysis of the above four pepper mitogenomes revealed that C. pubescens shares high sequence similarity with KJ865409 and that C. pubescens has rearranged with the other three pepper mitogenomes. Interestingly, we observed 72 similar sequences between the mitochondrial and chloroplast genomes, which accounted for 12.60% of the mitogenome, with a total length of 57,207 bp. These sequences encompassed 12 tRNA genes and the rRNA gene (rrn18). Remarkably, selective pressure analysis suggested that the nad5 gene underwent obvious positive selection. Furthermore, a single-base mutation in three genes (nad1, nad2, and nad4) resulted in an amino acid change. Conclusion This study provides a high-quality mitogenome of pepper, providing valuable molecular data for future investigations into the exchange of genetic information between pepper organelle genomes.

Abstract In the rapidly expanding domain of scientific research, tracking and synthesizing information from the rapidly increasing volume of publications pose significant challenges. To address this, we introduce a novel high-throughput pipeline that employs ChatGPT to systematically extract and analyze connectivity information from the full-texts and abstracts of 24,237 and 150,538 research publications concerning Caenorhabditis elegans and Drosophila melanogaster, respectively. This approach has effectively identified 200,219 and 1,194,587 interactions within the C. elegans and Drosophila biomaps, respectively. Utilizing Cytoscape Web, we have developed a searchable online biomaps that link relevant keywords to their corresponding PubMed IDs, thus providing seamless access to an extensive knowledge network encompassing C. elegans and Drosophila. Our work highlights the transformative potential of integrating artificial intelligence with bioinformatics to deepen our understanding of complex biological systems. By revealing the intricate web of relationships among key entities in C. elegans and Drosophila, we offer invaluable insights that promise to propel advancements in genetics, developmental biology, neuroscience, longevity, and beyond. We also provide details and discuss significant nodes within both biomaps, including the insulin/IGF-1 signaling (IIS) and the notch pathways. Our innovative methodology sets a robust foundation for future research aimed at unravelling complex biological networks across diverse organisms. The two databases are available at worm.bio-map.com and drosophila.bio-map.com.

Abstract Background Novel technologies are needed to combat anopheline vectors of malaria parasites as the reductions in worldwide disease incidence has stalled in recent years. Gene drive-based approaches utilizing Cas9/guide RNA (gRNA) systems are being developed to suppress anopheline populations or modify them by increasing their refractoriness to the parasites. These systems rely on the successful cleavage of a chromosomal DNA target site followed by homology-directed repair (HDR) in germline cells to bias inheritance of the drive system. An optimal drive system should be highly efficient for HDR-mediated gene conversion with minimal error rates. A gene-drive system, AgNosCd-1, with these attributes has been developed in the Anopheles gambiae G3 strain and serves as a framework for further development of population modification strains. To validate AgNosCd-1 as a versatile platform, it must perform well in a variety of genetic backgrounds. Results We introduced or introgressed AgNosCd-1 into different genetic backgrounds, three in geographically-diverse Anopheles gambiae strains, and one each in an An. coluzzii and An. arabiensis strain. The overall drive inheritance, determined by presence of a dominant marker gene in the F2 hybrids, far exceeded Mendelian inheritance ratios in all genetic backgrounds that produced viable progeny. Haldane’s rule was confirmed for AgNosCd-1 introgression into the An. arabiensis Dongola strain and sterility of the F1 hybrid males prevented production of F2 hybrid offspring. Back-crosses of F1 hybrid females were not performed to keep the experimental design consistent across all the genetic backgrounds and to avoid maternally-generated mutant alleles that might confound the drive dynamics. DNA sequencing of the target site in F1 and F2 mosquitoes with exceptional phenotypes revealed drive system-generated mutations resulting from non-homologous end joining events (NHEJ), which formed at rates similar to AgNosCd-1 in the G3 genetic background and were generated via the same maternal-effect mechanism. Conclusions These findings support the conclusion that the AgNosCd-1 drive system is robust and has high drive inheritance and gene conversion efficiency accompanied by low NHEJ mutation rates in diverse An. gambiae s.l. laboratory strains.

Abstract The Fasin rainbow fish, scientifically named Melanotaenia fasinensis, is highly prized by aquarium enthusiasts for its vibrant colors and adaptability to artificial aquatic environments. This species is endemic to the karst landscape of the Bird’s Head region in Papua, Indonesia, and belongs to the family Melanotaeniidae. Discovered relatively recently in 2010, this species was designated as endangered by the International Union for Conservation of Nature (IUCN) in 2021. However, there is currently insufficient data regarding its phylogenetic position. To address this gap, our study employed next-generation sequencing (NGS) to analyze the entire mitochondrial genome of M. fasinensis. The mitochondrial genome comprises 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes, with a total length of 16,731 base pairs. The base composition of the mitogenome revealed percentages of 27.76% adenine (A), 27.34% thymine (T), 16.15% guanine (G), and 28.75% cytosine (C) residues. Our phylogenetic analysis based on sequence data indicated that all species of the Melanotaeniidae family clustered together on the same branch. Furthermore, the intergeneric and interspecific taxonomic positions were explicit and clear. Phylogenetically, Melanotaeniidae were more closely related to the family Isonidae than to the family Atherinomorus. The phylogenetic position of M. fasinensis was relatively basal within the genus Melanotaenia. This study provides valuable molecular insights for further exploration of the phylogeography and evolutionary history of M. fasinensis and other members of the genus Melanotaenia.

Abstract Background Visium is a widely-used spatially-resolved transcriptomics assay available from 10x Genomics. Standard Visium capture areas (6.5mm by 6.5mm) limit the survey of larger tissue structures, but combining overlapping images and associated gene expression data allow for more complex study designs. Current software can handle nested or partial image overlaps, but is designed for merging up to two capture areas, and cannot account for some technical scenarios related to capture area alignment. Results We generated Visium data from a postmortem human tissue sample such that two capture areas were partially overlapping and a third one was adjacent. We developed the R/Bioconductor package visiumStitched, which facilitates stitching the images together with Fiji (ImageJ), and constructing SpatialExperiment R objects with the stitched images and gene expression data. visiumStitched constructs an artificial hexagonal array grid which allows seamless downstream analyses such as spatially-aware clustering without discarding data from overlapping spots. Data stitched with visiumStitched can then be interactively visualized with spatialLIBD. Conclusions visiumStitched provides a simple, but flexible framework to handle various multi-capture area study design scenarios. Specifically, it resolves a data processing step without disrupting analysis workflows and without discarding data from overlapping spots. visiumStitched relies on affine transformations by Fiji, which have limitations and are less accurate when aligning against an atlas or other situations. visiumStitched provides an easy-to-use solution which expands possibilities for designing multi-capture area study designs.

Abstract Background Heat stress has deleterious effects on physiological and performance traits in livestock. Within this context, using tropically adapted cattle breeds in pure herds or terminal crossbreeding schemes to explore heterosis is attractive for increasing animal production in warmer climate regions. This study aimed to identify biological processes, pathways, and potential biomarkers related to thermotolerance in Caracu, a tropically adapted beef cattle breed, by proteomic analysis of blood plasma. To achieve this goal, 61 bulls had their thermotolerance evaluated through a heat tolerance index. A subset of 14 extreme animals, including the seven most thermotolerant (HIGH group) and the seven least thermotolerant (LOW group), had their blood plasma samples used for proteomic analysis by liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). The differentially regulated proteins detected between HIGH and LOW groups were used to perform functional enrichment analysis and a protein-protein interaction network analysis. Results A total of 217 proteins were detected only in the HIGH thermotolerant group and 51 only in the LOW thermotolerant group. In addition, 81 and 87 proteins had significantly higher and lower abundancies in the HIGH group, respectively. Regarding proteins with the highest absolute log-fold change values, we highlighted those encoded by DUSP5, IGFALS, ROCK2, RTN4, IRAG1, and NNT genes based on their functions. The functional enrichment analysis detected several biological processes, molecular functions, and pathways related to cellular responses to stress, immune system, complement system, and hemostasis in both HIGH and LOW groups, in addition to terms and pathways related to lipids and calcium only in the HIGH group. Protein-protein interaction (PPI) network revealed as important nodes many proteins with roles in response to stress, hemostasis, immune system, inflammation, and homeostasis. Additionally, proteins with high absolute log-fold change values and proteins detected as essential nodes by PPI analysis highlighted herein are potential biomarkers for thermotolerance, such as ADRA1A, APOA1, APOB, APOC3, C4BPA, CAT, CFB, CFH, CLU, CXADR, DNAJB1, DNAJC13, DUSP5, FGA, FGB, FGG, HBA, HBB, HP, HSPD1, IGFALS, IRAG1, KNG1, NNT, OSGIN1, PROC, PROS1, ROCK2, RTN4, RYR1, TGFB2, VLDLR, VTN, and VWF. Conclusions Identifying potential biomarkers, molecular mechanisms and pathways that act in response to heat stress in tropically adapted beef cattle contributes to developing strategies to improve performance and welfare traits in livestock under tropical climates.

Abstract Background Bisphenol S (BPS) is the main substitute for bisphenol A (BPA), a well-known plasticiser and endocrine disruptor. BPS disrupts ovarian function in several species. Moreover, a few studies have reported that the effects of BPS might be modulated by the metabolic status, and none have characterised the granulosa cell (GC) proteome after chronic BPS exposure. Objectives This study aimed to decipher the mechanisms of action of chronic BPS exposure on the proteome of ewe GCs while considering the interaction between a deliberate contrasted metabolism and reproductive function. Methods Forty ewes were split into two groups with contrasted diets: restricted (R, n = 20) and well-fed (WF, n = 20). The R and WF ewes were subdivided according to the dose of BPS administered through the diet (0–50 µg/kg/day), forming four groups: R0, R50, WF0 and WF50. After 3-month BPS daily exposure, GCs were recovered during the pre-ovulatory stage and proteins were analysed by nano-liquid chromatography coupled with tandem mass spectrometry. Results Chronic exposure to BPS affected the GC proteome differently according to the ewe metabolic status. Fifty-nine out of 958 quantified proteins were differentially abundant between groups and are mainly involved in carbohydrate and lipid pathways. Unsupervised hierarchical clustering of differentially abundant proteins (DAPs) identified four clusters of 34, 6, 5 and 14 proteins according to the BPS exposure and diet interaction. Pairwise comparisons between groups also revealed a strong effect of BPS exposure and diet interaction. Functional analysis of DAPs highlighted that BPS upregulated β-glucuronidase (GUSB; p = 0.002), a protein especially able to deconjugate bisphenol glucuronides (BP-g). Moreover, among unexposed ewes, GUSB was detected only in well-fed ewes. Discussion Conjugation of glucuronides inhibits the oestrogenic activity of bisphenols. Upregulation of GUSB in ewes dosed with BPS would prolong the oestrogenic effects of BPS by deconjugating BPS-g into free BPS. In addition, literature has reported an up-regulation of GUSB in people suffering from obesity. Therefore, people suffering from obesity could be subjected to prolonged and aggravated exposure to BPS. These data highlighted the deleterious effects of BPS and its interaction with metabolic status.

Abstract Background In recent years, the use of endocrine-disrupting chemicals (EDCs) has become increasingly common, leading to severe environmental pollution and harm to aquatic organisms. 17α-Methyltestosterone (MT) is a synthetic androgen that can cause immunotoxicity in aquaculture, affecting fish health. To address this issue, this study aimed to investigate the effect of Vitamin C (VC) on MT-induced immunotoxicity and determine the optimal VC supplementation. Results Carassius auratus was exposed to 50 ng/L MT and treated with 25, 50, and 150 mg/kg VC for 7, 14, and 21 d. Morphological indicators, histological characteristics, hepatic antioxidant capacity, and immune-related gene expression were analyzed. Additionally, RNA-seq was performed on the liver tissues of the control, MT, and MT + 25 mg/kg VC groups after 21 d. Results showed that, MT treatment significantly increased liver malondialdehyde content and inhibited immune-related gene expression (TNF-α, IL-8, INF-γ, IL-10, Caspase-9, and IGF-I), causing oxidative stress and immunotoxicity, leading to hepatic steatosis. However, supplementation with 25–50 mg/kg VC effectively alleviated the MT-induced damage to the hepatic structure and immune system. RNA-seq revealed significant enrichment of differentially expressed genes in multiple signaling pathways, including the mTOR, MAPK, and Wnt pathways. Conclusions In summary, 25–50 mg/kg VC alleviated inhibitory effect of MT on immune-related genes in C. auratus liver, reducing MT-induced tissue damage. VC not only alleviated inflammation, oxidative stress, and immunotoxicity induced by MT through the regulation of the mTOR, MAPK, and Wnt signaling pathways, but also indirectly enhanced cellular antioxidant defense mechanisms by regulating the NRF2 pathway. This provides a theoretical basis for VC application in aquaculture, improving fish health and increasing efficiency.

Abstract Bergenia purpurascens is an important medicinal, edible and ornamental plant. The lack of omics information hinders the study of its metabolic pathways and related genes. In order to investigate candidate genes and pathways involved in secondary metabolism in B. purpurascens, roots, stems and leaves of B. purpurascens were subjected to metabolomic and transcriptomic analyses in this study. A total of 351 differentially accumulated secondary metabolites were identified. We identified 120 candidate genes involved in phenylpropanoid and flavonoid biosynthesis pathway, from which 29 key candidate genes were obtained by WGCNA. Five UDP-Glycosyltransferases and four O-methyltransferases were suggested to be the candidate enzymes involved in synthetic pathway from gallic acid to bergenin by correlation analysis between transcriptional and metabolic levels and phylogenetic analysis. This study provides data resources and new insights for further studies on the biosynthesis of major active components in B. purpurascens.

Abstract Background DNA methylation is a covalent bond modification that is observed mainly at cytosine bases in the context of CG pairs. DNA methylation patterns reflect the status of individual tissues, such as cell composition, age, and the local environment, in mammals. Genetic factors also impact DNA methylation, and the genetic diversity among various dog breeds provides a valuable platform for exploring this topic. Compared to those in the human genome, studies on the profiling of methylation in the dog genome have been less comprehensive. Results Our study provides extensive profiling of DNA methylation in the whole blood of three dog breeds using whole-genome bisulfite sequencing. The difference in DNA methylation between breeds was moderate after removing CpGs overlapping with potential genetic variation. However, variance in methylation between individuals was common and often occurred in promoters and CpG islands (CGIs). Moreover, we adopted contextual awareness methodology to characterize DNA primary sequences using natural language processing (NLP). This method could be used to effectively separate unmethylated CGIs from highly methylated CGIs in the sequences that are identified by the conventional criteria. Conclusions This study presents a comprehensive DNA methylation landscape in the dog blood. Our observations reveal the similar methylation patterns across dog breeds, while CGI regions showed high variations in DNA methylation level between individuals. Our study also highlights the potential of NLP approach for analyzing low-complexity DNA sequences, such as CGIs.

Abstract Background Outbreaks of the coral predator Crown-of-Thorns Starfish (CoTS) pose a severe threat to coral reefs in the Indo-Pacific Ocean. In 2018, the South China Sea (SCS) experienced significant CoTS outbreaks, leading to extensive coral mortality across the Xisha, Zhongsha, Dongsha, and Nansha Islands, severely impacting the coral reef ecosystem. Results To explore the origins of these outbreaks, we conducted a comprehensive genomic analysis using data from genomic single nucleotide polymorphism sites (SNPs) and mitochondrial haplotypes. Our analysis reveals that CoTS populations in the SCS, which exhibit moderate genetic diversity and may have undergone positive selection or population expansion. There was limited genetic differentiation among CoTS populations from XS, ZS, and NS groups. Especially between the XS and ZS groups, there was almost no genetic differentiation. The populations from XS, ZS, and NS groups have strong genetic connections with populations in Vietnam and the Philippines. There was high gene flow from Vietnam to the Xisha Islands and from the Philippines to the Nansha Islands, suggesting that the CoTS populations in these regions primarily originate from these neighboring countries. Conclusion The comprehensive analyses of SNP and mitochondrial genomes have provided valuable insights into the population genetics of CoTS. This research has generated significant genomic resources and facilitated important studies on the genetics of the CoTS species. By identifying potential source populations and understanding the genetic basis of their spread, managers can develop more effective conservation strategies to protect vulnerable coral reef ecosystems in the SCS.

Abstract Background Indigenous buffaloes, as the important livestock species contributing to economy of the country, are the lifeline of livelihood in Nepal. They are distributed across diverse geographical regions of the country and have adapted to various feeding, breeding, and management conditions. The larger group of these native buffalo breeds are present in narrow and stiff hilly terrains. Their dispersal indicates a possible environmental adaptation mechanism, which is crucial for the conservation of these breeds. Results We utilized whole-genome sequencing (WGS) to investigate the genetic diversity, population structure, and selection signatures of Nepalese indigenous buffaloes. We compared 66 whole-genome sequences with 118 publicly available sequences from six river and five swamp buffalo breeds. Genomic diversity parameters indicated genetic variability level in the Nepalese buffaloes comparable to those of Indian breeds, and population genetic structure revealed distinct geography-mediated genetic differentiation among these breeds. We used locus-specific branch length analysis (LSBL) for genome-wide scan, which revealed a list of potentially selected genes in Lime and Parkote breeds that inhabit the hilly region. A gene ontology (GO) analysis discovered that many GO terms were associated with cardiac function regulation. Furthermore, complementary analyses of local selection signatures, tissue expression profiles, and haplotype differences identified candidate genes, including KCNE1, CSF1R, and PDGFRB, related to the regulation of cardiac and pulmonary functions. Conclusions This study is a comprehensive WGS-based genetic analysis of the native Nepalese buffalo breeds. Our study suggested that the Nepalese “hilly” buffaloes, especially the Lime and Parkote breeds, have undergone some characteristic genetic changes and evolved increased cardiac and pulmonary function for their adaptation to the steep hilly terrains of the country.

Abstract Background As the inflorescence of wheat, spike architecture largely determines grain productivity. Dissecting the genetic basis for the spike morphology of wheat can contribute to the designation of ideal spike morphology to improve grain production. Results The present study characterizes a dense spike1 (ds1) mutant, derived from Nongda3753, induced by EMS treatment, which exhibits a dense spike and reduced plant height. Through bulked segregant analysis sequencing (BSA-Seq) of two segregating populations, ds1 was mapped to the short arm of chromosome 7B. Further genotypic and phenotypic analyses of the residual heterozygous lines from F3 to F6 of Yong3002×ds1 revealed that there was a 0-135 Mb deletion in chromosome 7B associated with the dense spike phenotype. The reads count analysis of the two bulks in BSA-Seq, along with the cytological analysis of ds1, ND3753, NIL-ds1 and NIL-Y3002, confirmed that the partial unidirectional translocation of 5AL (543–713 Mb) to 7BS (0-135 Mb) exists in ds1. This translocation led to an increase in both copy number and expression of the Q gene, which is one of the reasons for the dense spike phenotype observed in ds1. Conclusion Partial unidirectional translocation from 5AL to 7BS was identified in the EMS-induced mutant ds1, which exhibits dense spike phenotype. This research illustrates the effect of one chromosome structure variation on wheat spike morphology, and provides new materials with several chromosome structure variations for future wheat breeding.

Abstract Background Global warming-induced environmental stresses have diverse effects on gene expression and regulation in the life processes of various aquatic organisms. N6 adenylate methylation (m6A) modifications are known to influence mRNA transcription, localization, translation, stability, splicing, and nuclear export, which are pivotal in mediating stress responses. Apostichopus japonicus is a significant species in aquaculture and a representative of benthic organisms in ecosystems, thus there is a growing need for research on its heat stress mechanism. Results In this study, m6A-modified whole transcriptome profiles of the respiratory tree tissues of A. japonicus in the control (T18) and high-temperature stress (T32) groups were obtained using MeRIP-seq technology. The results showed that 7,211 common m6A peaks, and 9,459 genes containing common m6A were identified in three replicates T18 and T32 groups. The m6A peaks were found to be highly enriched in the 3’ untranslated region, and the common sequence of the m6A peak was also enriched, which was shown as RRACH (R = G or A; H = A, C, or U). A total of 1,200 peaks were identified as significantly differentially enriched in the T32 group compared with the T18 group. Among them, 245 peaks were upregulated and 955 were downregulated, which indicated that high temperature stress significantly altered the methylation pattern of m6A, and there were more demethylation sites in the T32 group. Conjoint analysis of the m6A methylation modification and the transcript expression level (the MeRIP-seq and RNA-seq data) showed co-differentiated 395 genes were identified, which were subsequently divided into four groups with a predominant pattern that more genes with decreased m6A modification and up-regulated expression, including HSP70IV, EIF2AK1, etc. GO enrichment and KEGG analyses of differential m6A peak related genes and co-differentiated genes showed the genes were significantly associated with transcription process and pathways such as protein processing in the endoplasmic reticulum, Wnt signaling pathway, and mTOR signaling pathway, etc. Conclusion The comparisons of m6A modification patterns and conjoint analyses of m6A modification and gene expression profiles suggest that m6A modification was involved in the regulation of heat stress-responsive genes and important functional pathways in A. japonicus in response to high-temperature stress. The study will contribute to elucidate the regulatory mechanism of m6A modification in the response of A. japonicus to environmental stress, as well as the conservation and utilization of sea cucumber resources in the context of environmental changes.

Abstract Background Microbes produce diverse bioactive natural products with applications in fields such as medicine and agriculture. In their genomes, these natural products are encoded by physically clustered genes known as biosynthetic gene clusters (BGCs). Genome and metagenome sequencing advances have enabled high-throughput identification of BGCs as a promising avenue for natural product discovery. BGC mining from (meta)genomes using in silico tools has allowed access to a vast diversity of potentially novel natural products. However, a fundamental limitation has been the ability to assemble complete BGCs, especially from complex metagenomes. With their fragmented assemblies, short-read technologies struggle to recover complete BGCs, such as the long and repetitive nonribosomal peptide synthetase (NRPS) and polyketide synthase (PKS). Recent advances in long-read sequencing, such as the High Fidelity (HiFi) technology from PacBio, have reduced this limitation and can help retrieve both accurate and complete BGCs from metagenomes, warranting improvement in the existing BGC identification approach for better utilization of HiFi data. Results Here, we present HiFiBGC, a command-line-based workflow to identify BGCs in PacBio HiFi metagenomes. HiFiBGC leverages an ensemble of assemblies from three HiFi-tailored metagenome assemblers and the reads not represented in these assemblies. Based on our analyses of four HiFi metagenomic datasets from four different environments, we show that HiFiBGC identifies, on average, 78% more BGCs than the top-performing single-assembler-based method. This increase is due to HiFiBGC’s ensemble assembly approach, which improves recovery by 25%, as well as from the inclusion of mostly fragmented BGCs identified in the unmapped reads. Conclusions HiFiBGC is a computational workflow for identifying BGCs in long-read HiFi metagenomes, implemented majorly using Python programming language and workflow manager Snakemake. HiFiBGC is available on GitHub at https://github.com/ay-amityadav/HiFiBGC under the MIT license. The code related to the figures and analyses presented in the manuscript is available at https://github.com/ay-amityadav/HiFiBGC_analyses .

Abstract Background Cephalanthera nanchuanica is a terrestrial orchid species and has been red listed as a second-grade protected plant due to its limited distributions in China. Initially classified within a monotypic genus Tangtsinia, this species was later reassigned to Cephalanthera based on morphological and molecular data. However, previous phylogenetic analyses of Cephalanthera using several segment sequences exhibited a low discriminatory power in delineating its relationships. Results In this study, we characterized and comparatively analyzed the complete chloroplast (cp) genome of C. nanchuanica with those of six previously reported Cephalanthera species. Our findings revealed that the cp genome of C. nanchuanica had the typical quadripartite structure, with a size of 161,365 bp and a GC content of 37.27%. A total of 113 unique genes were annotated, among which nearly half of protein-encoding genes (RSCU > 1) showed a preference in codon usage. No structural rearrangements were observed among the cp genomes of Cephalanthera species, except for C. humilis, which displayed structural alterations due to gene loss, relocation, and inverted repeat (IR) expansion/contraction. The cp genomes of Cephalanthera species were highly conserved, with only a small number of SSRs detected, most of which preferred A/T bases. Comparative analysis of cp genomes indicated that IR and coding regions were less divergent than single copy and non-coding regions and eight mutational hotspots were identified. Phylogenetic analysis suggested that the tribe Neottieae was a monophyletic group, divided into five clades. Palmorchis was the earliest-diverging lineage, followed by Cephalanthera. Diplandrorchis was deeply nested within Neottia, forming a clade. Aphyllorchis and Limodorum formed another clade, sister to Epipactis. Within the Cephalanthera clade, C. nanchuanica was sister to C. falcata with a strong support. Conclusions This study demonstrated that the cp genome characters of C. nanchuanica are highly similar to those of other Cephalanthera species, except for the mycoheterotrophic species C. humilis. Although the cp genomes of Cephalanthera species (excluding C. humilis) exhibited conservation in genome structure and sequence, SSR repeats and mutational hotspots were identified, which could potentially serve as as molecular markers for distinguishing Cephalanthera species. The phylogenetic analysis based on the protein-coding genes provided high-resolution support for the infrageneric classification. Therefore, cp genome data will be instrumental in resolving the phylogeny of the genus Cephalanthera.

Abstract Background Phloem protein 2 (PP2), a dimeric lectin, is known for its involvement in plant responses to biotic and abiotic stresses. However, research on PP2 proteins in Moso bamboo is lacking. Results In this study, comprehensive genome-wide analysis of the PP2-like gene family was conducted in Moso bamboo (Phyllostachys edulis), which has a significant economic and ecological value. Using HMMER3 search and InterPro domain analysis, 23 PP2-like genes (PhePP2-1 to PhePP2-23) were identified in the P. edulis genome. These genes were distributed across 12 chromosomal scaffolds, with proteins ranging from 216 to 556 amino acids in length. Phylogenetic analysis, including 163 PP2 proteins from eight plant species, revealed six distinct groups, with Group III and Group V being the largest. Gene structure and motif analyses indicated conserved domains across the PhePP2 proteins. In addition, Cis-element analysis of the promoter regions highlighted their potential regulatory roles in hormone, stress, and light responses. Expression pattern analysis using RNA-seq data showed differential expression of PhePP2 genes under drought, salt, salicylic acid, and abscisic acid treatments, indicating their involvement in stress response pathways. Furthermore, qPCR validation in different tissues and organs of Moso bamboo confirmed the expression profiles of the selected PhePP2 genes. Conclusions This study provides a comprehensive understanding of the functional roles of PP2-like genes in Moso bamboo and insights into their potential applications in enhancing stress tolerance and growth in plants.

Abstract The heat shock transcription factor (HSF) family is one of the most widely studied transcription factor families in plants; HSFs can participate in the response to various stressors, such as heat stress, high salt, and drought stress. Based on garlic transcriptome data, we screened and identified 22 garlic HSFs. The HSF proteins of garlic and Arabidopsis can be divided into three (A, B, C) subfamilies. The phylogenetic relationship, chromosome localization, sequence characteristics, conserved motifs, and promoter analysis of the HSF family were analyzed through bioinformatics methods. RT-qPCR analysis showed that the nine selected genes had different degrees of response to heat stress. In addition, we isolated and identified a class B HSF gene, AsHSFB1, from garlic variety ‘Xusuan No.6’. Subsequently, the AsHSFB1 gene was overexpressed in Arabidopsis thaliana. Under heat stress, the germination rate and growth of wild-type plants were better than that of transgenic plants. Moreover, after heat treatment, the contents of peroxidase, catalase, and chlorophyll a and b of transgenic plants were lower, but the contents of malondialdehyde (MDA) and leaf conductivity were higher. Nitroblue tetrazolium (NBT) staining showed that the stained area of transgenic plant leaves was larger than that of the wild type. Further studies showed that AsHSFB1 overexpression inhibited the expression of related reverse resistance genes. These results indicate that AsHSFB1 might play a negative regulatory role in garlic resistance under high stress. Altogether, these findings provide valuable data for revealing the function of HSF genes and lay a foundation for the subsequent selection of heat-resistant garlic varieties.

Abstract Background One of the most promising approaches for early and more precise disease prediction and diagnosis is through the inclusion of proteomics data augmented with clinical data. Clinical proteomics data is often characterized by its high dimensionality and extremely limited sample size, posing a significant challenge when employing machine learning techniques for extracting only the most relevant information. Although there is a wide array of statistical techniques and numerous analysis pipelines employed in proteomics data analysis, it is unclear which of these methods produce the most efficient, reproducible, and clinically meaningful results. Results In this study, we compared 9 unique analysis schemes comprised of different machine learning and dimensionality reduction methods for the analysis of simulated proteomics data consisting of 1317 proteins measured in 26 subjects (i.e., 13 controls and 13 cases). In scenarios where the sample size is extremely small (i.e., n

Abstract Background The locus for naked neck (Na) in chickens reduces feather coverage and leads to increased heat dissipation from the body surface resulting in better adaptability to hot conditions. However, the Na gene is linked to significantly lower hatchability due to an increased late embryonic mortality. It has been argued that the causative gene GDF7 may have a direct pleiotropic effect on hatchability via its effect on muscle development. Thus, the study presented here analyses the transcriptome of the hatching muscle (M. complexus) and shows how GDF7 impacts development leading to reduced hatching rates in Na chickens. Results Using 12 chicken embryos (6 x wildtype (Wt) and 6 x Na) RNA was extracted from the M. complexus of each embryo and sequenced. The resulting differential expression analyses led to the discovery of 461 differentially expressed (DE) genes in the M. complexus of the experimental group. Among those, 77 genes were of uncertain function (LOC symbols), with 31 were classified as uncharacterised. The regulation of a number of pathways involved in normal embryonic development, were found to be negatively influenced by the Na genotype. Further pathways involved in cell-cell adhesion, cell signalling pathways, and amino acid (AA) metabolism/transport were also observed. GDF7 (alias BMP12), whose localised overexpression in the neck skin causes the Na/Na phenotype, was significantly overexpressed in the M. complexus of Na/Na embryos, and shows a significant increase in the number of binding sites for the transcription factor PITX2 was also observed. Conclusion In Na chickens, GDF7 is under the control of a mutated cis-regulatory element, whose actions are known to suppress the development and distribution of feathers through the sensitizing action of retinoic acid. In this study, a number of DE genes with over 10 retinoic acid response elements (RAREs) in close proximity were observed, indicating changes to the retinol metabolism. With the understanding that the Na/Na mutation leads to increased retinoic acid activity, this indicates a high likelihood of GDF7 excerpting a direct pleiotropic effect, not just in the observed reduction in feather patterning, but also impacting the development of the M.complexus, and consequently leading to the reduced hatchability observed in birds with the Na/Na genotype. Furthermore, the enrichment of PITX2 binding sites in proximity to DE genes in the M. complexus, also indicates that muscle development is still ongoing in Na embryos. This suggests that the M. complexus is not yet fully developed, further increasing the potential for late embryonic mortality in Na chicks at hatching.

Abstract Background Circular RNAs (circRNAs) function as essential regulatory elements with pivotal roles in various biological processes. However, their expression profiles and functional regulation during the differentiation of goat myoblasts have not been thoroughly explored. This study conducts an analysis of circRNA expression profiles during the proliferation phase (cultured in growth medium, GM) and differentiation phase (cultured in differentiation medium, DM1/DM5) of skeletal muscle satellite cells (MuSCs) in goats. Results A total of 2,094 circRNAs were identified, among which 84 were differentially expressed as determined by pairwise comparisons across three distinct groups. Validation of the expression levels of six randomly selected circRNAs was performed using reverse transcription PCR (RT-PCR) and quantitative RT-PCR (qRT-PCR), with confirmation of their back-splicing junction sites. Enrichment analysis of the host genes associated with differentially expressed circRNAs (DEcircRNAs) indicated significant involvement in biological processes such as muscle contraction, muscle hypertrophy, and muscle tissue development. Additionally, these host genes were implicated in key signaling pathways, including Hippo, TGF-beta, and MAPK pathways. Subsequently, employing Cytoscape, we developed a circRNA-miRNA interaction network to elucidate the complex regulatory mechanisms underlying goat muscle development, encompassing 21 circRNAs and 47 miRNAs. Functional assays demonstrated that circTGFβ2 enhances myogenic differentiation in goats, potentially through a miRNA sponge mechanism. Conclusion In conclusion, we identified the genome-wide expression profiles of circRNAs in goat MuSCs during both proliferation and differentiation phases, and established that circTGFβ2 plays a role in the regulation of myogenesis. This study offers a significant resource for the advanced exploration of the biological functions and mechanisms of circRNAs in the myogenesis of goats.

Abstract Background Genomic structural variations (GSVs), notably copy number variations (CNVs), significantly shape genetic diversity and facilitate adaptation in cattle populations. Despite their importance, the genome-wide characterization of CNVs in indigenous Ethiopian cattle breeds—Abigar, Fellata, and Gojjam-Highland remains largely unexplored. In this study, we applied a read-depth approach to whole genome sequencing (WGS) data to conduct the first comprehensive analysis of CNVs in these populations. Results We identified 3,893 CNV regions (CNVRs) covering 19.15 Mb (0.71% of the cattle genome). These CNVRs ranged from 1.60 kb to 488.0 kb, with an average size of 4.92 kb. These CNVRs included deletions (1713), duplications (1929), and mixed events (251) showing notable differences in distribution among the breeds. Four out of five randomly selected CNVRs were successfully validated using real time polymerase chain reaction (qPCR). Further analyses identified candidate genes associated with high-altitude adaptation (GBE1 and SOD1), heat stress adaptation (HSPA13, DNAJC18, and DNAJC8) and resistance to tick infestations (BoLA and KRT33A). In addition, variance stabilizing transformation (V ST ) statistics highlighted population-specific CNVRs, emphasizing the unique genetic signatures of high-altitude adaptation in the Gojjam-Highland cattle breed. Among the detected CNVRs, 4.93% (192 out of 3,893) overlapped with 520 quantitative traits loci (QTLs) associated with six economically important trait categories suggesting that these CNVRs may significantly contribute to the genetic variation underlying these traits. Conclusions Our comprehensive analysis reveals significant CNVRs associated with key adaptive traits in Ethiopian cattle breeds highlighting their genetic diversity and resilience. These findings offer valuable insights into the genetic basis of adaptability and can inform sustainable breeding practices and conservation efforts. Future research should prioritize the functional validation of these CNVRs and their integration into breeding programs to enhance traits such as disease resistance and environmental adaptability.

Abstract Mangiferin, a C-glucosyl xanthone, is a biologically active glycoside naturally synthesized in mango. Glycosyltransferase can catalyze the biosynthesis of mangiferin. In this study, we identified 221 members of the UGT glycosyltransferase family in mango. The 221 MiUGT genes were grouped into 13 subfamilies through phylogenetic tree analysis with Arabidopsis, Chinese bayberry, and mango. All UGT family members in mango were unevenly distributed on 17 chromosomes and found that tandem duplication dominated the expansion of UGT family members in mango. Purification selection primarily influenced the evolution of the mango UGT family members. In addition, cis-element analysis of the mango UGT gene family revealed the presence of MYB binding sites, which are involved in flavonoid biosynthesis; which further supports the role of UGT family members in the synthesis of flavonoids. To verify these results, we analyzed the expression of UGT family members in mango leaves, stems, and different developmental stages of fruit peel. The RNA-seq and qRT-PCR results showed significant differences in the expression patterns of MiUGT genes in various tissues and developmental stages of mango. We identified MiUGT gene-specific expression at different stages of fruit development. These results lay a theoretical foundation for research on the relationship between members of the mango UGT family and the synthesis of flavonoids, mangiferin.

Abstract Alterations in splicing patterns of leukemic cells have a functional impact and influence most cellular processes since aberrantly spliced isoforms can provide a proliferative advantage, enable to evade apoptosis, induce metabolic reprogramming, change cell signaling and antitumor immune response, or develop drug resistance. In this Review, we first characterize the general mechanism of mRNA processing regulation with a focus on the role of splicing factors, which are commonly mutated in blood neoplasms. Next, we provide a comprehensive summary on the current understanding of alternative splicing events, which confer resistance to targeted treatment strategies and immunotherapy. We introduce the functional consequences of mis-spliced variants (CD19-∆ex2, CD22-∆ex2, CD22-∆ex5-6, CD33-∆ex2, PIK3CD-S, BCR-ABL 35INS, BIM-γ, FPGS-8PR, dCK-∆ex2-3, and SLC29A1-∆ex13) production in leukemic cells. Of therapeutic relevance, we summarize novel strategies focused on pharmacological correction of aberrant splicing, including small-molecule splicing modulators and splice-switching oligonucleotides. We also include the findings of recent preclinical investigation of the antisense strategies based on modified oligonucleotides. Finally, we discuss the potential of emerging combination therapies for the treatment of hematological disorders with disrupted splicing.

Abstract Background Bivalves represent a vital economic resource in aquaculture for their high productivity and extensive market demand. Growth is one of the most important and desired aquaculture traits for bivalves, regulated by multiple levels, notably intricate metabolic processes. However, the understanding of the metabolic profiles that influence bivalve growth is limited, particularly from a multi-tissue perspective. Results In this study, metabolic profiles of multiple tissues of Chlamys farreri with different growth performance were systematically investigated by ultraperformance liquid chromatography quadrupole time-of-flight mass spectrometry (UPLC-Q-TOF-MS). Through comparing the metabolic variation between fast-growing (FG) scallops and slow-growing (SG) scallops, 613, 509, 105, and 192 significantly different metabolites (SDMs) were identified in the mantle, gill, adductor muscle, and digestive gland, respectively. Growth-related metabolic pathways including sphingolipid metabolism, fatty acid biosynthesis, and ABC transporter pathway, along with 11 SDMs associated with growth traits were identified in all four tissues, implying they were involved in the growth of multiple tissues in scallops. Tissue-specific metabolic profiling indicated that sulfur-containing amino acid metabolism in the mantle potentially contributed to shell growth, while the gill synergistically participated with the mantle through various metabolic processes, such as tyrosine metabolism, glycine, serine, and threonine metabolism and melanogenesis; energy metabolism was crucial for adductor muscle growth; and nutrients digestion and absorption in the digestive gland were linked to scallop growth. Conclusions Our results represent the first comprehensive analysis of the crucial pathways and metabolites associated with the growth of C. farreri, offering valuable insights for future bivalve aquaculture production.

Abstract Background Single-cell RNA sequencing experiments commonly use 10x Genomics (10x) kits due to their high-throughput capacity and standardized protocols. Recently, Parse Biosciences (Parse) introduced an alternative technology that uses multiple in-situ barcoding rounds within standard 96-well plates. Parse enables the analysis of more cells from multiple samples in a single run without the need for additional reagents or specialized microfluidics equipment. To evaluate the performance of both platforms, we conducted a benchmark study using biological and technical replicates of mouse thymus as a complex immune tissue. Results We found that Parse detected nearly twice the number of genes compared to 10x, with each platform detecting a distinct set of genes. The comparison of multiplexed samples generated from 10x and Parse techniques showed 10x data to have lower technical variability and more precise annotation of biological states in the thymus compared to Parse. Conclusion Our results provide a comprehensive comparison of the suitability of both single-cell platforms for immunological studies.

Abstract Background Atrial natriuretic peptide (ANP) has been associated with ischemic stroke (IS), but the underlying molecular mechanisms are not clear. The objective of this study was to examine whether DNA methylation and hydroxymethylation in the coding gene of ANP (NPPA) were associated with IS, as well as its functional outcome. Methods and results DNA methylation and hydroxymethylation in NPPA promoter region were quantified by targeted bisulfite sequencing and APOBEC-coupled epigenetic sequencing, respectively, in 615 IS patients and 610 age- and sex-matched healthy controls. True methylation was calculated as the quantified methylation level minus the quantified hydroxymethylation level. The functional outcome of IS was evaluated by the modified Rankin Scale (mRS) at a 3-month follow-up visit after onset. Multiple testing was controlled by the false discovery rate approach. Of the nine CpG sites assayed, hypomethylation at eight CpGs was not only associated with an increased risk of having IS but also predicted a higher mRS score after onset (all q

Abstract Background Inhomogeneous patterns of chromatin-chromatin contacts within 10–100-kb-sized regions of the genome are a generic feature of chromatin spatial organization. These features, termed topologically associating domains (TADs), have led to the loop extrusion factor (LEF) model. Currently, our ability to model TADs relies on the observation that in vertebrates TAD boundaries are correlated with DNA sequences that bind CTCF, which therefore is inferred to block loop extrusion. However, although TADs feature prominently in their Hi-C maps, non-vertebrate eukaryotes either do not express CTCF or show few TAD boundaries that correlate with CTCF sites. In all of these organisms, the counterparts of CTCF remain unknown, frustrating comparisons between Hi-C data and simulations. Results To extend the LEF model across the tree of life, here, we propose the conserved-current loop extrusion (CCLE) model that interprets loop-extruding cohesin as a nearly conserved probability current. From cohesin ChIP-seq data alone, we derive a position-dependent loop extrusion rate, allowing for a modified paradigm for loop extrusion, that goes beyond solely localized barriers to also include loop extrusion rates that vary continuously. We show that CCLE accurately predicts the TAD-scale Hi-C maps of interphase Schizosaccharomyces pombe, as well as those of meiotic and mitotic Saccharomyces cerevisiae, demonstrating its utility in organisms lacking CTCF. Conclusions The success of CCLE in yeasts suggests that loop extrusion by cohesin is indeed the primary mechanism underlying TADs in these systems. CCLE allows us to obtain loop extrusion parameters such as the LEF density and processivity, which compare well to independent estimates.

Abstract Background The safety of CRISPR-based gene editing methods is of the utmost priority in clinical applications. Previous studies have reported that Cas9 cleavage induced frequent aneuploidy in primary human T cells, but whether cleavage-mediated editing of base editors would generate off-target structure variations remains unknown. Here, we investigate the potential off-target structural variations associated with CRISPR/Cas9, ABE, and CBE editing in mouse embryos and primary human T cells by whole-genome sequencing and single-cell RNA-seq analyses. Results The results show that both Cas9 and ABE generate off-target structural variations (SVs) in mouse embryos, while CBE induces rare SVs. In addition, off-target large deletions are detected in 32.74% of primary human T cells transfected with Cas9 and 9.17% of cells transfected with ABE. Moreover, Cas9-induced aneuploid cells activate the P53 and apoptosis pathways, whereas ABE-associated aneuploid cells significantly upregulate cell cycle-related genes and are arrested in the G0 phase. A percentage of 16.59% and 4.29% aneuploid cells are still observable at 3 weeks post transfection of Cas9 or ABE. These off-target phenomena in ABE are universal as observed in other cell types such as B cells and Huh7. Furthermore, the off-target SVs are significantly reduced in cells treated with high-fidelity ABE (ABE-V106W). Conclusions This study shows both CRISPR/Cas9 and ABE induce off-target SVs in mouse embryos and primary human T cells, raising an urgent need for the development of high-fidelity gene editing tools.

Abstract Spatial transcriptomics is a cutting-edge technique that combines gene expression with spatial information, allowing researchers to study molecular patterns within tissue architecture. Here, we present IAMSAM, a user-friendly web-based tool for analyzing spatial transcriptomics data focusing on morphological features. IAMSAM accurately segments tissue images using the Segment Anything Model, allowing for the semi-automatic selection of regions of interest based on morphological signatures. Furthermore, IAMSAM provides downstream analysis, such as identifying differentially expressed genes, enrichment analysis, and cell type prediction within the selected regions. With its simple interface, IAMSAM empowers researchers to explore and interpret heterogeneous tissues in a streamlined manner.

Abstract Unveiling functional relationships between various molecular cell phenotypes from data using machine learning models is a key promise of multiomics. Existing methods either use flexible but hard-to-interpret models or simpler, misspecified models. VI-VS (Variational Inference for Variable Selection) balances flexibility and interpretability to identify relevant feature relationships in multiomic data. It uses deep generative models to identify conditionally dependent features, with false discovery rate control. VI-VS is available as an open-source Python package, providing a robust solution to identify features more likely representing genuine causal relationships.

Abstract Spatial transcriptomics is revolutionizing the exploration of intratissue heterogeneity in cancer, yet capturing cellular niches and their spatial relationships remains challenging. We introduce SpottedPy, a Python package designed to identify tumor hotspots and map spatial interactions within the cancer ecosystem. Using SpottedPy, we examine epithelial-mesenchymal plasticity in breast cancer and highlight stable niches associated with angiogenic and hypoxic regions, shielded by CAFs and macrophages. Hybrid and mesenchymal hotspot distribution follows transformation gradients reflecting progressive immunosuppression. Our method offers flexibility to explore spatial relationships at different scales, from immediate neighbors to broader tissue modules, providing new insights into tumor microenvironment dynamics.

Abstract Background and objectives Behçet’s disease (BD) is an ongoing systemic vasculitis with mouth and genital ulceration and eye, skin, and systemic affection. It has considerable morbidity and mortality, and a high incidence and a higher male-to-female affection characterize Egypt. We aimed to evaluate liver x receptor messenger RNA beta (LXR β) and microRNA-146a (miR-146a) gene expression in Behçet’s disease cases in Egypt to relate them with clinicopathological features. Subject and method Eighty Egyptian individuals were split into two groups for the current study: Forty Behçet’s disease cases and forty healthy controls matched by age and gender from the Rheumatology Department at Helwan University Hospital, Egypt. The transformed score, also known as the Behçet’s Disease Current Activity Form (BDCAF), was utilized to measure disease activity. To assess the amounts of LXR β and miR-146a serum expression via real-time PCR, a blood sample was obtained. Results There was a downregulation of both LXR β and miR-146a levels that significantly differed between the BD group and the control group (p = 0.0001 and 0.0001, respectively). There is a noteworthy inverse relationship between the expression level of LXR β and BDCAF Patients Index (r = − 0.79 and p = 0.0001) was found. Regarding miR-146a, it had a reverse correlation with BDCAF Patients Index (r = − 0.89 and p = 0.0001). Conclusion LXR β and miR-146a were found to be significant non-invasive predictor biomarkers for Behçet’s disease and can indicate disease activity.

Abstract Background Stage 4 neuroblastoma (NBL), a solid tumor of childhood, has a poor prognosis. Despite intensive molecular genetic studies, no targetable gene abnormalities have been identified. Stage 4S NBL has a characteristic of spontaneous regression, and elucidation of the mechanistic differences between stages 4 and 4S may improve treatment. Conventional NBL studies have mainly focused on the detection of abnormalities in individual genes and have rarely examined abnormalities in gene networks. While the gene coexpression network is expected to contribute to the detection of network abnormalities, the fragility of the network due to data noise and the extraction of arbitrary topological structures for the high-dimensional network are issues. Results The present paper concerns the classification method of stages 4 and 4S NBL patients using highly accurate gene coexpression network analysis based on RNA-sequencing data of transcription factors (TFs). In particular, after applying a noise reduction method RECODE, generalized topological overlapping measure (GTOM), which weighs the connections of nodes in the network structure, succeeded in extracting a cluster of TFs that showed high classification performance for stages 4 and 4S. In addition, we investigated how these clusters correspond to clinical information and to TFs which control the normal adrenal tissue and NBL characters. Conclusions A clustering method is presented for finding intermediate-scale clusters of TFs that give considerable separation performance for distinguishing between stages 4 and 4S. It is suggested that this method is useful as a way to extract factors that contribute to the separation of groups from multiple pieces of information such as gene expression levels.

Abstract Background Innate immune deficiencies can impair both cellular and humoral immune responses. In contrast, other immune functions may appear normal, leading to increased susceptibility to specific pathogens, such as severe viral infections or Mendelian Susceptibility to Mycobacterial Disease (MSMD). Studying these deficiencies is essential for understanding the pathophysiology of these infectious diseases. Main body While primary immunodeficiencies (PIDs) generally cause vulnerability to multiple infections, innate immunodeficiencies increase susceptibility to specific pathogens, despite normal immune responses to others. Patients with these deficiencies show normal immunoglobulins and lymphocyte subpopulations, complicating diagnosis. This review highlights genetic susceptibility to mycobacteria, pneumococci, herpes simplex virus, and candidiasis, emphasizing recognizing this subset of PIDs. Conclusion This review highlights the diverse spectrum of genetic mutations contributing to defects in innate and intrinsic immunity, including Mendelian susceptibility to mycobacterial disease (MSMD), chronic mucocutaneous candidiasis, and predispositions to invasive bacterial and viral infections. Identifying key mutations in pathprovideh such as TLR3, IFN signaling, and IL-17A/F immunity provides valuable insights into the pathogenesis of these conditions. Our findings underscore the need for early genetic diagnosis and targeted interventions, particularly in regions with high undiagnosed cases, to reduce the morbidity and mortality associated with defects in innate and intrinsic immunity.