Back

Genomics

Elsevier BV

Preprints posted in the last 30 days, ranked by how well they match Genomics's content profile, based on 60 papers previously published here. The average preprint has a 0.08% match score for this journal, so anything above that is already an above-average fit.

1
Genome-wide analysis of Heavy metal ATPase (P1B-type ATPase) gene family in Mung bean and their expression analysis under heavy metal (Zn, Cd and Cu) stress

Panigrahi, J.; Panigrahy, D.; Rath, B.; Gupta, K.

2026-03-27 molecular biology 10.64898/2026.03.25.713876 medRxiv
Top 0.1%
6.8%
Show abstract

Heavy metal ATPases (HMAs) are important group of transmembrane proteins involved in homeostasis of metal ions in plant systems. In this study, a comprehensive analysis of genome assembly (VC1973A v7.1) resulted in the identification of nine HMA genes (VrHMA) and their corresponding proteins in Mungbean, an agronomically important legume crop known for its nutritional values. VrHMA proteins were also characterized based on their biomolecular features, conserved domains and motifs arrangement, transmembrane helices, pore-line helices, subcellular location and occurrence of signal peptides. Based on sequence homology, nine VrHMAs were clustered into two major substrate-specific groups: VrHMA1, VrHMA5 and VrHMA7 were categorized under the Zn/Co/Cd/Pb ATPase group, whereas the remaining six VrHMAs belong to the Cu/Ag subgroup. Gene structure analysis and promoter scanning revealed the structural divergence and presence of various stress-responsive cis-acting elements, respectively. The expression analysis of VrHMA genes in root and leaf tissues, in response to heavy metal (Zn, Cd and Cu) stress, indicates their role in the uptake, transport and sequestration of metal ions. Interestingly, VrHMA5 showed incremental upregulation in roots in response to all three heavy metal stresses, whereas its expression was only upregulated in the leaf tissues under Zn stress, which indicates its role in vascular transport in V. radiata. In addition, this study provides valuable insights into the functional roles of VrHMA genes and will lay a foundation for future genetic improvement in mung bean aimed at enhanced heavy metal stress tolerance and micronutrient homeostasis.

2
Genetic and heat-stress related environmental influences on pig whole-blood gene expression levels

Durante, A.; Feve, K.; Naylies, C.; Labrune, Y.; Gress, L.; Lippi, Y.; Legoueix, S.; Milan, D.; Gourdine, J.-L.; Gilbert, H.; Renaudeau, D.; Riquet, J.; Devailly, G.

2026-03-18 genomics 10.64898/2026.03.17.712411 medRxiv
Top 0.1%
5.0%
Show abstract

BackgroundGene expression levels are affected by genetics and environmental effects. However, quantification of the influence of genetics and environmental effects on gene expression remains limited, especially in farm animals. Here, the relative influence of genetic and heat-related environmental variations on gene expression levels was investigated in pigs, using a backcross herd of diverse heat adaptation levels. Backcross animals were raised in either a tropical or temperate environment. Animals raised in temperate environment were subjected to an experimental heat stress at the end of their growth. ResultsWe identified 1,967 differentially expressed genes (DEGs) between pigs raised in the tropical (n = 181) and temperate (n = 180) facilities, and 472 DEGs throughout a 3 weeks experimental heat stress. Transcriptome-wide association (TWAS) study identified 139 associations between gene expression levels and thermoregulation/production traits. We detected 6,014 expression quantitative trait loci (eQTLs) associated with the expression level of 3,297 genes. Genetic variance was estimated to explain 36.3% of gene expression variance on average, and was the main source of variance for 27.7% of transcripts. Most eQTLs found are located in proximal regions (cis-eQTLs) and few within distal regions (trans-eQTLs) to their assigned genes. A trans-eQTL hotspot highlighted a hematopoietic mechanism driven by GPATCH8. An integration of GWAS and TWAS pointed to TMCO1 and ZNF184 as candidate genes for backfat thickness. ConclusionsThis study provides a better understanding of the impact of climate, heat stress and genetic influences on the pig whole blood transcriptome.

3
Heat Stress Induces Locus-Specific DNA Hypomethylation Linked to Immune Regulation in Lactating Holstein Cows

Costa Monteiro Moreira, G.; Ruiz Gonzalez, A.; Joigner, M.; Costes, V.; Chaulot-Talmon, A.; Ali, F.; Bourgeois-Brunel, L.; Jammes, H.; Rico, D. E.

2026-03-26 genomics 10.64898/2026.03.23.713208 medRxiv
Top 0.4%
3.6%
Show abstract

Epigenetics may play a crucial role in livestock adaptation to environmental challenges like heat stress. In recent years, a growing number of studies have investigated the epigenetic mechanisms underlying dairy cow adaptation to heat stress. However, there is still limited knowledge about the effects of heat stress on immune cells and immune-related phenotypes. Herein we aim to identify heat-stress induced DNA methylation variations on blood methylome potentially affecting regulatory regions and associated phenotypes. Blood samples were collected and peripheral blood mononuclear cell (PBMC) isolated from four cows before (D0) and after (D14) a 14-d heat stress challenge (cyclical THI 72-82) and, from four cows kept in thermoneutral conditions (THI 61-64). Heat-stressed cows had ad libitum access to diets supplemented with adequate levels of vitamin D and Ca (12,000 IU/kg of vitamin D and 0.73% Ca, respectively). To eliminate confounding effects due to differences in nutrient intake, cows maintained under thermoneutral conditions were pair-fed (PF) to their heat-stressed counterparts and received adequate concentrations of vitamin D and Ca as well. Reduced representation bisulphite sequencing (RRBS) was used to profile PBMCs methylome. Differential methylation analysis was performed using methylKit and DSS softwares ({Delta}meth [&ge;] 25%, adjusted p-value < 0.01), retaining only commonly detected differentially methylated cytosines (DMCs). A total of 2,908 DMCs were identified when comparing pre- and post-heat stress samples. After excluding 649 DMCs that were also detected under thermoneutral conditions, as these changes were likely associated with feed restriction inherent to the pair-feeding design rather than with heat stress per se, 2,259 heat stress-specific DMCs remained, predominantly hypomethylated. About half of the DMCs are annotated in intronic and intergenic regions; known to harbor regulatory elements. By intersecting the DMRs with publicly available functional annotation data, we observed hypomethylation on regulatory regions putatively affecting cows immune system. As an example, we identified a loss of methylation within an enhancer region of the MSN gene, which is involved in lymphocyte homeostasis, and a loss of methylation in the promoter region of MECP2, a well-established epigenetic regulator with a central role in chromatin organization and gene expression. These findings highlight the impact of heat stress on dairy cow immunity and provide new insights into its epigenetic regulation under environmental stress. Interpretative summaryThis study examined DNA methylation changes induced by heat stress in dairy cows to elucidate epigenetic mechanisms of thermal adaptation. Using RRBS on PBMCs, 2,259 heat stress-specific differentially methylated cytosines were identified, predominantly hypomethylated and enriched in regulatory regions. Functional annotation highlighted immune-related pathways, including hypomethylated regulatory regions near genes (e.g., MSN, ZBTB33, SLC25A5, GNAS, FAM3A, and MECP2) associated with immune function. These findings indicate that heat stress induces targeted epigenetic modifications potentially affecting immune regulation in dairy cows.

4
Integrative Identification and Characterization of PCOS-Associated lncRNAs From the Interface of Genetic Association, Transcriptomics, and Gene Structure Evolution

He, Z.; Li, Y.; Shkurat, T. P.; Butenko, E. V.; Derevyanchuk, E. G.; Lomteva, S. V.; Chen, L.; Lipovich, L.

2026-04-02 genomics 10.64898/2026.03.31.715548 medRxiv
Top 0.4%
3.3%
Show abstract

BackgroundPolycystic ovary syndrome (PCOS) is a prevalent endocrine disorder and a leading cause of female infertility, with complex genetic, metabolic, and hormonal etiologies. Long non-coding RNAs (lncRNAs) have emerged as important regulators of diverse biological processes, yet their roles in PCOS remain underexplored. Here, we identified and characterized PCOS differentially expressed gene-associated lncRNAs (PDEGAL) with an integrative approach combining expression data, genetic association, and evolutionary analysis. MethodsThirty-three PCOS-associated protein-coding genes were obtained from our prior study, and all their nearby and overlapping lncRNAs were annotated. These candidates were analyzed using UCSC Genome Browser-mapped annotations and datasets, including NCBI RefSeq, GENCODE, GTEx, GWAS SNPs, and conservation, as well as the FANTOM5 cap analysis of gene expression (CAGE) promoter data, to assess their expression, regulatory potential, genetic variant overlaps, and evolutionary conservation. ResultsTwenty-three PDEGALs (18 antisense to, and 5 sharing bidirectional promoters with, known PCOS-associated protein-coding genes) were identified. 17 PDEGALs contained GWAS SNPs with statistically significant disease associations, 9 of which were associated with PCOS-related traits. 5 PDEGALs demonstrated expression in the KGN granulosa cell model of PCOS. Key gene structure element (KGSE) analysis revealed that most PDEGALs are primate-specific. Integrating four criteria--GTEx expression, GWAS SNPs, FANTOM promoterome, and KGSE conservation--highlighted HELLPAR as the only lncRNA fulfilling all four, while five others--PGR-AS1, MTOR-AS1, ENSG00000265179, ENSG00000256218, and LOC105377276--fulfilled three of the four criteria. ConclusionsWe have systematically identified candidate PCOS regulatory lncRNAs with convergent genetic, expression, and evolutionary evidence. These results provide a framework for functional validation and highlight lncRNAs as potential biomarkers and therapeutic targets in PCOS that function by regulating their nearby and overlapping protein-coding genes.

5
Computational Prediction of Plasmodium falciparum Antigen-T-cell Receptor Interactions via Molecular Docking: Implications for Malaria Vaccine Design

Kipkoech, G.; Kanda, W.; Irungu, B.; Nyangi, M.; Kimani, C.; Nyangacha, R.; Keter, L.; Atieno, D.; Gathirwa, J.; Kigondu, E.; Murungi, E.

2026-03-20 bioinformatics 10.64898/2026.03.18.712575 medRxiv
Top 0.4%
3.1%
Show abstract

Malaria is one of the deadliest diseases in sub-Saharan Africa and Southeast Asia. The majority of the fatalities occur mostly in children under 5 years and pregnant women and this is due to infection by Plasmodium spp, of which Plasmodium falciparum is the most virulent and is responsible for most of the morbidity and mortality. Despite various public health interventions such as use of insecticide-treated bed nets, spraying of homes with insecticides and use of WHO recommended artemisinin-based combination therapies (ACT), malaria prevention still faces major setback due to drug and insecticide resistance by P. falciparum and mosquitoes respectively. The study uses molecular docking and immunoinformatics to screen various Plasmodium spp antigens and evaluate their antigenicity and suitability as vaccine candidates. The P. falciparum antigens and T-cell receptor (TCR) structures were obtained from Protein Data Bank (PDB) based on a range of factors related to their role in the lifecycle of the parasite and their status as vaccine targets. Protein structures not available in the PDB were predicted using AlphaFold. The 3D structures of selected P. falciparum antigens and TCR structures were downloaded in PDB format then all water molecules, Hetatm, and bound ligands were deleted from the protein structures using BIOVIA Discovery Studio Visualizer. Subsequently, molecular docking was done using ClusPro v2.0 server and docked complexes were compared. The findings of this study gave valuable insights into the interaction of human immune response with P. falciparum antigens. The best three ranked antigen complexes are PfCyRPA, PfMSP10 and PfCSP and this confirm their use as potential candidates for vaccine development. This study highlights the usefulness of computational docking in identifying P. falciparum antigens of excellent immunogenic potential as vaccine candidates.

6
Comparative LUSZ Therapeutic Study (LUSZ_AVIST) of Antiviral, Antiretroviral, and Immunosuppressive Treatments in Hospitalized COVID-19 Patients with High-Risk Factors, Biomarkers, and Disease Progression.

Makdissy, N.; Makdessi, E. W.; Fenianos, F.; Nasreddine, N.; Daher, W.; El Hamoui, S.

2026-04-13 respiratory medicine 10.64898/2026.04.10.26350587 medRxiv
Top 0.5%
2.7%
Show abstract

COVID-19 has spread rapidly and caused a global pandemic making it one of the deadliest in history. Early identification of patients with coronavirus disease 2019 who may develop critical illness is of immense importance. Therefore, novel biomarkers were needed to identify patients who will suffer rapid disease progression to severe complications and death. Many treatments were adopted including the antiviral Remdesivir, the antiretroviral Lopinavir /Ritonavir and Tocilizumab. Our study aimed not only to specify high-risk factors and biomarkers of fatal outcome in hospitalized subjects with coronavirus but also to compare the efficacy of the three considered treatments to help clinicians better choose a therapeutic strategy and reduce mortality. We divided the population (n=711) into four main groups based according to the WHO ordinal severity scale. The percentage of mortality, in and out the hospital, the length of stay in the hospital, the pulmonary inflammatory lesion and its distribution, the SARS-CoV-2 IgM and IgG variations at admission, the inflammatory markers, the complete blood count, the coagulation factors and enzymes, proteins and electrolytes profile, glucose and lipid profile, and other relevant markers were measured. The significance of the observed variation was assessed by multivariate and ANOVA analyses. We succeeded to establish a novel predictive scoring model of disease progression based on a cohort of Lebanese hospitalized patients relying on the pulmonary inflammatory lesions, inflammation biomarkers such as LDH, D-Dimer, CRP, IL-6 and the lymphocyte count, the number of comorbidities and the age of the patient which all were significantly correlated with the illness severity showing best outcomes with immunomodulatory and anticoagulant treatments by the results. As top tier, Tocilizumab was more efficient than the two other treatments in non-severe cases but none of the used treatments was insanely effective alone to reduce mortality in severe cases.

7
lncOriL, a novel polyadenylated mitochondrial lncRNA common to zebrafish and human

Jorgensen, T. E.; Wardale, A.; Wolf Profant, S.; Amundsen, C.; Emblem, A.; Joakimsen, I. S.; Brekke, O.-L.; Karlsen, B. O.; Babiak, I.; Johansen, S. D.

2026-03-27 molecular biology 10.64898/2026.03.26.714394 medRxiv
Top 0.6%
2.1%
Show abstract

Even though teleost fish and mammals share the same mitochondrial gene content and organization, the teleost mitochondrial transcriptome is still poorly understood. We characterized the mitochondrial transcriptome during zebrafish (Danio rerio) early development by long-read direct RNA sequencing. All heavy-strand specific mRNAs were found to carry 3 poly-A tails of approximately 50-60 residues, and the transcriptome profile was distinctive but practically invariant between stages. Three unusual transcripts were however noted. These included two mRNAs (COI and ND5 mRNAs), with significant 3 untranslated regions corresponding to antisense gene sequences, and a previously not described noncoding RNA named here lncOriL. The ND5 mRNA was found to carry one third of all detected m6A methylation sites in the zebrafish mitochondrial transcriptome. The 313 nt-long lncOriL transcript had an abundance comparable to that of ND5 mRNA and it mapped to mitochondrial genome region covering the origin of light strand replication and four flanking antisense tRNAs. A mitochondrial tRNA-derived fragment (tiRNA5-Asn), with a 35 nt perfect pairing-potential to lncOriL, was present at all stages. Additional analyses including adult zebrafish, scissortail (Rasbora rasbora), and monkfish (Lophius piscatorius) strongly corroborate the results of COI mRNA, ND5 mRNA, and lncOriL transcript prevalence among teleost fish. Surprisingly, our findings in zebrafish were further supported by mitochondrial transcriptome analyses in domestic pig (Sus scrofa) and human (Homo sapiens), including tiRNA5-Asn commonly present in human tissues, suggesting that lncOriL is ubiquitously expressed and regulated in vertebrates. Author SummaryMitochondria contain their own genome and produce essential RNAs needed for energy production. Although fish and mammals share the same mitochondrial gene organization, less is known about how mitochondrial RNAs are processed and regulated in teleost. Using Nanopore direct RNA sequencing, we examined mitochondrial RNAs during early zebrafish development and discovered three unusual transcripts that include extended non-coding regions. Two of these molecules, COI and ND5 mRNAs, carry long 3' untranslated regions formed by antisense gene sequences, suggesting previously unrecognized regulatory potential. We also identified lncOriL, a highly structured long noncoding RNA that spans the origin of light-strand replication and is abundant during development. Strikingly, the same RNA feature, including lncOriL and a matching tRNA-derived small RNA (tiRNA5-Asn), was found not only in zebrafish but also in human mitochondrial transcriptomes. These findings support conservation of regulatory mitochondrial RNAs across main groups of vertebrate species. Our work reveals a new layer of mitochondrial RNA regulation and expands the current understanding of how mitochondrial gene expression is controlled.

8
Reference genomes of four miniature and non-miniature cypriniform fishes inhabiting acidic peat-swamp forest blackwaters of Southeast Asia

Sudasinghe, H.; Liu, Z.; Triginer-Llabres, L.; Hui Tan, H.; Britz, R.; Salzburger, W.; Peichel, C.; Rueber, L.

2026-03-24 genomics 10.64898/2026.03.21.713365 medRxiv
Top 0.8%
1.9%
Show abstract

The acidic blackwaters of Southeast Asias peat-swamp forests represent some of the most extreme freshwater environments on Earth. Despite their very low pH values, limited nutrients, and hypoxic conditions, these blackwater habitats harbor a remarkable diversity of freshwater fishes, including multiple lineages that have independently adapted to these extreme conditions and, in some cases, exhibiting extreme body miniaturization. These replicate evolutionary lineages therefore provide a powerful comparative framework to investigate adaptation to extreme environments and the genomic basis of miniaturization. Here, we present high-quality, annotated reference genomes for four cypriniform species endemic to these peat-swamp forest ecosystems: Paedocypris sp., Sundadanio atomus, Boraras brigittae, and Rasbora kalochroma. The first two are progenetic miniatures, including Paedocypris, comprising the smallest known fish, while B. brigittae represents a proportioned dwarf and R. kalochroma a non-miniature taxon. Genome sizes ranged from 401-1,290 Mb and heterozygosity from 0.34-1.7%. All genome assemblies achieved pseudo-chromosome-level contiguity, high k-mer completeness (>99%), and high BUSCO completeness (94.5-98.9%). Repeat analyses revealed lineage-specific differences in transposable element landscapes and abundances, while gene annotation identified notable intron length reduction in progenetic miniatures.

9
Comparative cross-species transcriptomic analysis identifies new candidates of Pooideae nitrate response

Gregoire, M.; Pateyron, S.; Brunaud, V.; Tamby, J. P.; Benghelima, L.; Martin, M.-L.; Girin, T.

2026-03-20 plant biology 10.64898/2026.03.18.712634 medRxiv
Top 0.8%
1.9%
Show abstract

AO_SCPLOWBSTRACTC_SCPLOWNitrogen fertilizers are essential for crop productivity but cause environmental harm, necessitating the development of cultivars that thrive under limited nitrogen. This study investigates the transcriptomic response to nitrate in Arabidopsis thaliana (a model dicot), Brachypodium distachyon (a model Pooideae), and Hordeum vulgare (barley, a domesticated Pooideae) to identify conserved and species-specific molecular mechanisms. Using RNA-seq after 1.5 and 3 hours of nitrate treatment, we found that core nitrate-responsive biological processes - such as nitrate transport, assimilation, carbon metabolism, and hormone signaling - are largely conserved across species. However, comparative analysis at gene level based on orthology revealed specificities between the species. For instance, rRNA processing was uniquely stimulated in Arabidopsis, while cysteine biosynthesis from serine and gibberellin biosynthesis were specifically regulated in Brachypodium and barley. Orthologs of key nitrate-responsive genes (e.g., NRT, NLP, TCP20) exhibited variable regulation, reflecting potential adaptations linked to domestication or nutrient acquisition strategies. These findings highlight the importance of integrating model and crop species to uncover targets for improving nitrogen use efficiency in cereals. The study provides a pipeline integrating gene ontology and orthology analyses to compare transcriptomic responses between species.

10
A comparative analysis of liver tissue and novel primary organoid cultures from ruminants reveals species-specific immune architecture and metabolic specialization

Garner, M. E.; Price, D. R. G.; McCarron, P.; Bartley, D. J.; Faber, M. N.; Quinn, B.; Robinson, M. W.; Smith, D.

2026-04-06 cell biology 10.64898/2026.04.01.715896 medRxiv
Top 1%
1.7%
Show abstract

The liver is widely considered to be one of the most conserved organs amongst vertebrates, with it being involved in blood detoxification, bile production and the metabolism of xenobiotic compounds. Liver organoids have previously been derived from several species and used as models of drug metabolism, toxicity, and fundamental tissue biology. To date, however, these models have not been developed from ruminant species, specifically cattle and sheep. Here we present the first report of the development and comprehensive characterisation of bovine and ovine liver organoids derived from primary liver tissue. When initially established, organoids from both species were comprised of KRT19- and KRT18-positive cholangiocytes. The capacity for organoids to differentiate into hepatocyte-enriched cultures was evaluated and it was noted that there was an increase in hepatocyte markers in bovine cultures. A comparative analysis of the liver tissue and organoids of both species revealed species-specific differences in gene expression, which were conserved within organoid cultures. Most notably, bovine liver tissue and organoids had enriched expression of genes associated with fatty acid uptake and storage whereas ovine samples had higher expression of genes associated with fatty acid conversion, highlighting fundamental differences between these two ruminant species. Differences in expression of cytochrome P450 family genes were identified alongside those associated with an inflammatory response specifically in bovine samples, whereas ovine samples had higher expression of genes associated with a protective immune response. Despite this, transcriptomic analysis of organoids from both species, cultured in both growth and differentiation media, revealed preserved expression of genes associated with key liver functions, including gluconeogenesis and xenobiotic metabolism. Transcripts associated with the flavin-containing monooxygenases (FMO) family were expressed in both organoid growth media and organoid development media (OGM and ODM respectively), and both species could metabolise triclabendazole into its primary metabolite triclabendazole sulfoxide, therefore validating the potential of the organoids to be applied as in vitro models of metabolism and/or toxicity. Overall, this study provides novel insights into differences in liver composition and function between ruminant species, as well as providing novel experimental models of the liver for both cattle and sheep.

11
The Effect of Vitamin-D Supplementation on HDAC2 Levels in Stable COPD Patients

Donastin, A.; Irawan, D.; Effendy, E.; Iryawan, R. D. A.; Nuari, N.; Oktaviana, B. M.; Yahya, D.; Muhammad, A. R.

2026-04-08 respiratory medicine 10.64898/2026.04.05.26348641 medRxiv
Top 1%
1.7%
Show abstract

Background: Chronic Obstructive Pulmonary Disease (COPD) is the third leading cause of global mortality, with persistent lung inflammation contributing to disease progression. This inflammation is partly associated with reduced levels of histone deacetylase 2 (HDAC2). Previous studies suggest that Vitamin D may modulate HDAC2 levels. This study aimed to evaluate the effect of Vitamin D supplementation on HDAC2 expression in stable COPD patients. This experimental study aimed to evaluate the effect of vitamin D supplementation on HDAC2 expression in stable COPD patients at Jemursari Islamic Hospital. Methods: Five COPD patients received a daily dose of 5000 IU of Vitamin D for three months. Serum levels of 25(OH)D3 and HDAC2 were measured before and after the intervention. Results: Vitamin D supplementation resulted in a significant increase in both 25(OH)D and HDAC2 levels. Pulmonary function parameters showed an increasing trend, however, no statistically significant differences were observed. Conclusion: Vitamin D supplementation was associated with increased HDAC2 levels, suggesting a potential anti-inflammatory effect. However, no significant improvement in pulmonary function was observed. Further studies are needed to determine its clinical impact.

12
The genome of the Delisea pulchra: a resource for the study of chemical host-microbe interactions in red algae

Dittami, S. M.; Hudson, J.; Brillet-Gueguen, L.; Ficko-Blean, E.; Tanguy, G.; Rousvoal, S.; Legeay, E.; Markov, G. V.; Delage, L.; Godfroy, O.; Corre, E.; Collen, J.; Leblanc, C.; Egan, S.

2026-04-02 genomics 10.64898/2026.03.31.715562 medRxiv
Top 1%
1.7%
Show abstract

BackgroundRed macroalgae (Rhodophyta) are ecologically and economically important marine primary producers, yet genomic resources for most species remain scarce. Delisea pulchra, a temperate red alga known for its halogenated furanone-based chemical defenses, serves as a model for studying algal-microbe interactions, antifouling mechanisms, and disease dynamics. ResultsHere we present a high-quality genome assembly of this species. The nuclear genome comprises 134 Mbp across 271 contigs with an N50 of 1.47 Mbp and encodes 13,387 predicted protein-coding genes. Comparative genomics with other red algae revealed expansions in gene families involved in DNA methylation, and oxidative stress responses, including glutathione S-transferases and superoxide dismutases. Analysis of glycosyltransferases, sulfatases, and sulfurylases implicated in galactan biosynthesis suggests D. pulchra possesses a complex and potentially novel extracellular matrix. We also identified several vanadium haloperoxidases (vHPOs), heme-dependent haloperoxidases (hHPOs), and two type III polyketide synthase (PKS) genes unique to the D. pulchra, which together represent promising candidate genes for bromofuranone production. ConclusionThe D. pulchra genome provides a foundation for molecular investigations into defense, signaling, and host-microbe interactions. It has been deposited at the European Nucleotide Archive under accession number PRJEB101077. All datasets, annotations, and interactive tools for exploring the genome are also available through the Rhodoexplorer portal at https://rhodoexplorer.sb-roscoff.fr.

13
Evolutionary history of ligand binding by the LRR domain of innate immunity receptors: the story of the TLR2 cavity

Namou, R.; Ichii, K.; Takkouche, A.; Jaroszewski, L.; Godzik, A.

2026-03-30 bioinformatics 10.64898/2026.03.26.714386 medRxiv
Top 1%
1.7%
Show abstract

Toll-like receptors (TLRs) are vital components of the innate immune system, recognizing both exogenous pathogens signals (PAMPs) and internal stress signals (DAMPs). TLR2 is unique among the human (Homo sapiens) TLR family members, as it contains a large cavity for binding hydrophobic ligands, such as lipoteichoic acid (LTA) and di/triacyl lipopeptides (Pam2/3CSK4). This study analyzed the structural phylogeny of cavity presence in the TLR2 lineage in vertebrates (vTLR) enabled by AI protein structure predictions and explored the potential convergent evolution of similar features in invertebrates (iTLRs). Analysis of AI models of TLR2s shows that this cavity is consistently present in TRL2 orthologs across jawed vertebrates (Gnathostomata). In jawless vertebrates (Cyclostomatha), these cavities were found in lamprey (Petromyzon marinus) TLR2 model, but only in some extant hagfish (Myxini), suggesting an ancestral origin in basal vertebrates followed by lineage-specific losses. TLR2 paralogs were found in several species, with a similar central cavity but potentially different ligand specificities. In silico ligand docking showed Pam2CSK4 binds to this cavity in all TLRs and paralogs consistently, demonstrating the conserved function of the ligand-binding pocket in gram-positive bacteria recognition across TLR2 branches. Changes in the TLR2 cavity size and shape in some vertebrate groups show the evolution of this DAMP recognition mechanism adapted to its respective pathogens. iTLRs form a separate phylogenetic branch with distinct structural features, but in literature some are considered to be TLR2 orthologs. Indeed, TLRs from some species of Helobdella and Ciona, contain a cavity with some similarity to that in the vTLR2 lineage. However, detailed structural comparisons of their location in the LRR domain and the structural details of the models suggest that their cavities have developed independently from that in TLR2s. Smaller cavities are present in other branches of the LRR family, but show different locations, shapes, and features, indicating that the binding of small ligands in the internal cavities within the LRR domains evolved multiple times in the LRR domain family history.

14
Distinctive Lacritin Cleavage-Potentiated Bactericidal Alteration of the P. aeruginosa Transcriptome

Sharifian Gh, M.; Norouzi, F.; Laurie, G. W.

2026-03-18 microbiology 10.64898/2026.03.17.712309 medRxiv
Top 1%
1.6%
Show abstract

Lacritin is a tear, saliva, plasma and cerebrospinal fluid glycoprotein with broad polypharmacology. Selective deficiency of its bioactive monomeric form appears to be deleterious for ocular surface health for which replacement therapy is beneficial. Its cleavage-potentiated C-terminus represented by the N-104 proteoform in tears is bactericidal and synergizes with the tear thrombin peptide GKY20. In the pathogenic and multidrug resistant PA14 strain of P. aeruginosa, we recently discovered that N-104 binds to the outer-membrane lipoprotein YaiW to gain access to the periplasm where it targets and inhibits the inner-membrane ferrous iron transporter FeoB (of FeoABC) as well as PotH, a subunit of the polyamine transporter PotFGHI. Further, PA14 gene expression shifts toward anaerobic respiratory pathways. Here we explore N-104-associated transcriptional changes over a broader range of functional categories pointing to a reduction in: (i) virulence by suppressed gene expression of virulence factors AprA and LasA; and Hcp1 and PsrA necessary for the respective assembly of type VI and III secretion systems, (ii) fitness (less AtsC, MgtA), (iii) metabolism (less AdhA, AtsC, GcvH2, GcvP2, FadE1, SsuD, SsuF, TauB, TauD, UspK, UspN), (iv) stress response (less UnG, PfpI, RmF), (v) proteostasis (less ClpB, GrpE, HtpG), (vii) quorum sensing (less CifR, GcvH2, GcvP2, PsrA, QuiP), and (viii) survival under anaerobic conditions (less AdhA, MhR, ModA, UspKLNO). Upregulated genes are directed towards enhancing PA14: (i) multidrug (more OprJ of MexCD-OprJ) and (ii) tellurite (more TerC) efflux, coupled with a seemingly PA14 survival attempt at (iii) anaerobic respiration (more NosR), (iv) translational fidelity (more QueE, RimP, TrmD) and (v) metabolism (CysT, MoaA1, Sbp, SsuA, SsuE). The overlap with aminoglycosides (4.3%), {beta}-lactams (0%), cyclic peptides (2.5%), fluoroquinilones (0%) and macrolide (1.9%) classes of antibiotics in P. aeruginosa was minimal. Thus, N-104 appears to widely perturb PA14 fundamental processes in a distinctive manner.

15
Investigating the dynamics of heat acclimation in pig through transcriptome analysis of blood samples

Huau, G.; Liaubet, L.; Labrune, Y.; Campos, P. H. R. F.; Gilbert, H.; Renaudeau, D.

2026-04-06 systems biology 10.64898/2026.04.01.715954 medRxiv
Top 1%
1.5%
Show abstract

This study aimed to investigate the dynamics of gene expression in pigs during heat stress (HS), focusing on both short-term (STHA) and long-term (LTHA) heat acclimation phases. A total of 12 castrated males were exposed to thermoneutral temperatures (24{degrees}C) for 14 days (TN) and then to a constant temperature of 30{degrees}C for 21 days. Rectal temperature measurements indicated a biphasic thermoregulatory response, with an initial peak followed by acclimation. Using whole blood transcriptome analysis at seven time points between day 5 before the initiation of HS challenge and day 13 post HS. A total of 525 genes were differentially expressed during the STHA (day 0-day 2) phase. A switch in the expression of most genes was observed around 20 hours after HS. Functional pathway enrichment analysis identified through shape-based clustering revealed the activation of the immune system, especially mediated through toll-like receptor signaling pathways. The LTHA phase (day 2-day 13) revealed 985 differentially expressed genes, with pathways associated with various metabolisms, including mitochondrial fatty acid beta-oxidation, and electron transport, ATP synthesis, and heat production by uncoupling proteins. Interestingly, oxidative phosphorylation was predicted to be activated during the LTHA, particularly in Complex V, whereas other complexes showed mixed regulation. Comparative pathway analysis indicated distinct metabolic adaptations between STHA and LTHA, with up-regulation of glucose and lipid metabolism in late STHA and down-regulation of lipid metabolism during LTHA. This study contributes to a better understanding of the time course of adaptation mechanisms in pigs to HS, underlying a coordinated regulation during STHA involving several stress-specific mechanisms (via the HSP) and metabolic variation to help pigs achieve homeothermy.

16
Identification of a microRNA with a mutation in the loop structure in the silkworm Bombyx mori

Harada, M.; Tabara, M.; Kuriyama, K.; Ito, K.; Bono, H.; Sakamoto, T.; Nakano, M.; Fukuhara, T.; Toyoda, A.; Fujiyama, A.; Tabunoki, H.

2026-03-27 molecular biology 10.64898/2026.03.24.714027 medRxiv
Top 1%
1.5%
Show abstract

MicroRNAs (miRNAs) play essential roles in the posttranscriptional regulation of gene expression in organisms. In the process of synthesizing mature miRNAs from miRNA precursors, the miRNA precursors are cleaved via Dicer at their loop structure, after which the miRNA precursors become mature and regulate transcription. However, the consequences of altering the loop sequence are not fully understood. The silkworm Bombyx mori is a lepidopteran insect with many genetic strains. We identified a mutant of the miRNA miR-3260 whose the part of the loop structure was lacking in a silkworm strain with translucent larval skin. Here, we aimed to analyze the role of wild-type miR-3260 and the influence of the mutation of the loop structure in B. mori. First, we identified the genomic region responsible for the translucent larval skin phenotype and determined that the mutated miR-3260 nucleotide sequences. Then, we predicted the binding partners of wild-type miR-3260 using the RNA hybrid tool and found two juvenile hormone (JH)-related genes as targets of wild-type miR-3260. Next, we assessed the relationships between miR-3260 and JH and found that miR-3260 was highly expressed in the Corpora allata and its expression responded to JH treatment. Meanwhile, miR-3260 mimic and inhibitor did not induce the typical phenotypes associated with JH in B. mori. Then, we compared the dicing products from wild-type and mutant miR-3260 precursors and observed that neither form underwent Dicer-mediated cleavage when the loop structure was altered. These results suggest that loop mutations in the miR-3260 precursor may not influence dicing activity, consistent with the lack of observable phenotypic effects.

17
Genome sequence of Tacca chantrieri reveals the genetic basis of floral pigmentation

de Oliveira, J. A. V. S.; Pucker, B.

2026-03-19 plant biology 10.64898/2026.03.17.712415 medRxiv
Top 1%
1.3%
Show abstract

Tacca chantrieri, black bat flower, has showy flowers often appearing almost black. Here, we present the genome sequence and corresponding annotation to identify the genetic basis of the pigmentation. Candidate genes associated with the anthocyanin biosynthesis were identified based on this genome sequence and investigated with respect to their properties. The best dihydroflavonol 4-reductase (DFR) candidate, which harbours all amino acid residues believed to be required for DFR activity, shows a threonine in the substrate preference determining position where most characterized DFRs display asparagine or aspartate. This amino acid residue appears to be frequent in the Dioscoreaceae family as a comprehensive investigation revealed.

18
Integrated Transcriptomic and Network-Based Identification of Prognostic Hub Genes in Oral Squamous Cell Carcinoma

Choudhary, S.; Guleria, V.

2026-04-06 bioinformatics 10.64898/2026.04.02.716250 medRxiv
Top 1%
1.3%
Show abstract

BackgroundThe most prevalent kind of oral cancer is oral squamous cell carcinoma (OSCC), which has a poor prognosis because of delayed detection and a lack of molecular indicators. MethodsTranscriptomic data from TCGA were analyzed to identify differentially expressed genes between OSCC and normal samples. Functional enrichment analysis was performed to determine biological pathways. A protein-protein interaction network was constructed using STRING and visualized in Cytoscape to identify hub genes. ResultsA total of 5732 differentially expressed genes were identified, including 2459 upregulated and 3273 downregulated genes. Network analysis revealed several highly connected hub genes such as CDK1, CCNB1, TOP2A, BUB1, and MMP9. Functional enrichment indicated significant involvement of cell cycle regulation and cancer-associated pathways. ConclusionThis integrative analysis identified key regulatory hub genes that may be involved in OSCC progression. These genes may serve as promising biomarkers and therapeutic targets for future studies.

19
GWAS Reveals Distinct Genetic Architecture of Schistosomiasis-Induced Hepatic Fibrosis with DGKG as a Key Mediator

Zhou, M.; Xue, C.; Zhang, L.; Hu, Y.; Ning, A.; Wang, L.; Shen, J.; Song, L.; Zhang, B.; Liu, J.; Liao, Y.; Chen, Z.; Khan, J.; Wu, Z.; Chen, C.; Sun, X.; Wu, X.; Li, M.

2026-03-24 genetic and genomic medicine 10.64898/2026.03.21.26348960 medRxiv
Top 2%
1.2%
Show abstract

Schistosomiasis is a major cause of hepatic fibrosis in endemic regions, yet the host genetic determinants of disease progression remain poorly defined. We aimed to identify genetic drivers and underlying mechanisms of schistosomiasis-induced hepatic fibrosis. We performed a genome-wide association study (GWAS) of 984 Schistosoma japonicum-infected individuals from hyperendemic areas in China followed by multi-omics integration and experimental validation to identify causal genes and fibrogenic pathways. Schistosomiasis-associated fibrosis exhibited a genetic architecture distinct from metabolic and viral liver fibrosis, supporting pathogen-specific mechanisms. Eight novel susceptibility loci were identified, including a genome-wide significant signal at 16p13 (rs73575170, P = 3.9 x 10-8). Integrative mapping linked these loci to 262 genes enriched in liver sinusoidal endothelial cells (P = 5.84 x 10-5) and sphingolipid metabolism pathways (P = 4.19 x 10-5). Notably, Diacylglycerol kinase gamma (DGKG, rs6762330, P = 4.37 x 10-6) emerged as a key candidate, with its expression in peri-granuloma and periportal hepatocytes strongly correlating with fibrosis severity (r = 0.816). In vivo, Dgkg knockout attenuated hepatic fibrosis and immunopathology while restoring cholesterol homeostasis, whereas Dgkg overexpression exacerbated fibrogenesis and increased TNF-{beta} levels tenfold. This study identifies DGKG as a key mediator linking lipid metabolism and immune signaling in schistosomiasis-induced fibrosis, uncovering a pathogen-specific genetic mechanism and providing a potential therapeutic target for infection-associated liver fibrosis.

20
Clarified an rDNA Gene Unit Pattern with (CTTT)n and (CT)n Microsatellites Aggregation Ahead of and Behind the Gene in Human Genome

Shen, J.; Tang, S.; Xia, Y.; Qin, J.; Xu, H.; Tan, Z.

2026-03-24 genetics 10.64898/2026.03.22.713381 medRxiv
Top 2%
1.2%
Show abstract

BackgroundConventional models of human ribosomal DNA (rDNA) array organization have historically depended on transcription-centric boundaries, partitioning the unit into a [~]13 kb rDNA transcription region and a monolithic [~]31 kb intergenic spacer (IGS). While our previous identification of Duplication Segment Units (DSUs) mapped these arrays based on an intuitive analysis of the microsatellite density landscape of the complete reference human genome, our present deep mining of this landscape has revealed a more accurate rDNA Gene Unit Pattern. Methods & ResultsIn this study, we conducted a deep mining analysis of our previously established microsatellite density landscape of the T2T-CHM13 assembly, focusing specifically on nucleolar organizing regions (NORs). We suggest a more accurate rDNA Gene Unit Pattern containing a (CTTT)n microsatellite aggregation ahead of the rDNA gene and a (CT)n microsatellite aggregation behind the gene, rather than a pattern featuring an IGS region inserted between two rDNA genes. ConclusionsA correct rDNA gene pattern of the human genome probably includes a (CTTT)n microsatellite aggregation ahead of the gene and a (CT)n microsatellite aggregation behind it, which possibly constitute cis- and trans-regulating regions; the (CTTT)n and (CT)n microsatellite aggregations may provide two different local stable DNA structures for regulatory protein binding.