Search PubMedSearch

SEARCH · Search PubMed

Results for “Minerals-associated genes”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,556 records · Page 15Linked to original sources

The genomic origins and evolutionary path to a key innovation in the world's most venomous snakes.

Evolutionary innovation is a catalyst for the colonization of new environments and the adaptive radiations of major groups. Novel traits typically evolve through the modification of preexisting characters, but the genetic paths underlying their origin have been challenging to trace, and the general requirements for and relative order of different kinds of gene mutations have been difficult to assess. Here, we trace the genomic origins of four procoagulant venom toxins (factor X, factor V, group I phospholipase A2, and Kunitz-type toxins) that collectively underlie a novel, especially potent blood-clotting venom type in the recently evolved Australian brown snake and taipan clade. We find evidence for a previously unknown fifth toxin, coagulation factor VII, and show that the toxins evolved through two distinct genetic paths. The factor X and factor V toxins evolved through the sequential de novo co-option of ancestral clotting factor proteins that entailed their heterotopic expression in the venom gland, the fixation of segmental duplications containing each locus, and subsequent gain-of-function mutations that rendered factor X and factor V constitutively active. In contrast, the phospholipase A2 and Kunitz-type toxins evolved by modifying the functions of neurotoxins that were part of the venom arsenal. Our findings support models in which innovative mutations in single-copy genes precede gene duplication in the evolution of novel proteins and offer a rare view into the genesis of a complex trait that has played a central role in a major adaptive radiation.

Animals

Integrated single-cell and bulk transcriptomic analysis identifies a novel senescent fibroblast subtype associated with poor prognosis in acral melanoma.

BACKGROUND: Acral melanoma (AM) exhibits significant intratumoral heterogeneity, but its tumor microenvironment (TME) and immune regulation remain unclear. This study aims to dissect TME heterogeneity and establish a prognostic model based on key cell subpopulations. METHODS: We collected AM single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA). Unsupervised clustering, CellChat, and Scissor analysis were performed to characterize cellular heterogeneity, cell-cell communication, and prognosis-related cell subpopulations. Kaplan-Meier analysis was used to assess the prognostic value of key genes, which were further validated by multiplex immunohistochemistry (mIHC). RESULTS: In AM, Mel_C2, C7, and C9 with high SEMA6A and KIT expression were strongly linked to poor prognosis. We further identified a senescent fibroblast subpopulation (sCAF_CDKN2A) characterized by high fibroblast senescence signature (FSS) scores. Integrating Scissor analysis of fibroblast subtypes with bulk prognostic data, we identified COL3A1, VCAN, and KIT as prognosis-associated genes upregulated in poor-outcome-related fibroblast subsets. Cell-cell communication analysis revealed that sCAF_CDKN2A engages in an immunosuppressive network, interacting with regulatory T cells (Tregs) via MIF signaling and receiving signals from exhausted CD8+ T cells through PPIA-BSG interactions. Using transcription factor expression patterns from these fibroblast subtypes, we constructed a prognostic model that effectively stratified patients into distinct risk groups with significant differences in overall survival (OS). mIHC confirmed significantly higher protein levels of SEMA6A and COL3A1 in tumor tissues compared to matched normal tissues. CONCLUSIONS: We established a novel prognostic model for AM and identified sCAF_CDKN2A as an immunosuppressive senescent fibroblast subpopulation driving poor prognosis.

Acral melanoma

Gloeotrichia echinulata genomes from the United States are nontoxigenic and likely geosmin producers.

Six Gloeotrichia echinulata genomes derived from planktonic harmful algal blooms (HABs) with similar colonial morphology have been sequenced from lakes in the west and northeast regions of USA, four of them to completion. The c. 7 Mbp genomes exhibit a high level of conservation, with 98-99% pairwise genome-wide average nucleotide identity and high levels of synteny, representing a single species cluster. We observed strong conservation of gene clusters responsible for the synthesis of the secondary metabolites and bioactive peptides that are characteristic of HAB-forming cyanobacteria. All six G. echinulata genomes lack genes for the synthesis of classic cyanotoxins, including microcystin, but possess genes responsible for the synthesis of the taste and odor compound geosmin. Interestingly, the geoA geosmin synthase gene in three genomes is homologous to other cyanobacterial geoA genes, while the other three geoA genes are related to actinomyces geoA. Phylogenomic analysis places the G. echinulata genomes within a clade of benthic Nostocales, reflecting an ecological niche featuring extensive growth on the sediment surface before colonies disperse into the epilimnion for planktonic growth. We identify genes conserved in all six genomes that could represent physiological adaptations supporting active growth on sediments and pelagic recruitment independent of wind-driven mixing: phycoerythrin light harvesting complexes for optimal photosynthesis at depth; gliding motility to access patchy nutrient distributions; and gas vesicles with relatively small GvpC proteins that predict resistance to higher hydrostatic pressure. The strong genomic similarity across geographically distant populations suggests that G. echinulata in the United States is a tightly related non-toxigenic species group with predictable properties relevant to public health and drinking water management.

Cyanobacteria

CNOT1 is a potential YTHDF2 target that orchestrates maternal mRNA decay and zygotic genome activation during goat embryogenesis.

Timely and efficient degradation of maternal mRNA is essential for early embryonic development, which occurs from fertilization through the initiation of zygotic genome activation (ZGA). Yet, the regulatory mechanisms governing this process remain poorly characterized. In the present study, we investigated the function of CCR4-NOT transcription complex subunit 1 (CNOT1) during goat embryogenesis. We found that CNOT1 was upregulated during mammalian ZGA, and that its knockdown led to developmental arrest and a marked reduction in blastocyst formation. Moreover, CNOT1 knockdown impaired nascent RNA activity, resulting in 814 upregulated and 1014 downregulated genes, which were enriched for RNA splicing, regulation of chromosome organization, and RNA localization. RNA splicing analysis revealed differential splicing events in 2959 genes, of which 259 were downregulated following CNOT1 knockdown. Notably, CNOT1 was predicted to crosstalk with the m6A reader YTHDF2. Knockdown of YTHDF2 resulted in CNOT1 downregulation at the 8-cell stage in goats and increased transcription levels around polyadenylation sites during ZGA in mice. Together, these findings indicate that CNOT1 is a potential YTHDF2 target that orchestrates maternal mRNA decay and ZGA during goat embryogenesis. Our work provides new insight into the complex regulatory landscape underlying ZGA and may inform strategies to improve the efficiency of goat embryogenesis.

Animals

Comparative analysis of histological and transcriptomic characteristics in caudal muscles of nile crocodiles (Crocodylus niloticus), siamese crocodiles (Crocodylus siamensis), and their hybrids.

Crocodylus niloticus and Crocodylus siamensis are high-value aquaculture species. C. niloticus is large-bodied but less abundant, while C. siamensis grows fast but is small-sized. Their hybrids combine parental advantages, yet relevant research is scarce. This study compared the histological and transcriptomic characteristics of the caudal muscle across the three taxa. HE staining indicated that C. niloticus had significantly larger myofiber diameters (p&#xa0;<&#xa0;0.05); C. siamensis had the smallest, and the myofiber density of hybrids was much closer to that of C. siamensis. Masson's trichrome staining indicated that C. niloticus had the thickest collagen fibers (p&#xa0;<&#xa0;0.05), C. siamensis the thinnest, and hybrids exhibited highly similar histological traits to C. siamensis. C. niloticus had higher LDH and SDH activities in caudal muscles, whereas the hybrid crocodile indicated the highest CK activity. Transcriptomic analysis identified numerous differentially expressed genes (DEGs), which were enriched in growth, muscle metabolism, and energy allocation pathways via GO/KEGG annotations. PPI analysis screened 24 hub genes related to energy metabolism. This study systematically reveals caudal muscle differences, providing insights into growth-related molecular mechanisms and theoretical support for crocodile artificial breeding.

Animals

Functions of tandem-repeat galectins and domain coordination governs galectin-4 activity in grass carp (Ctenopharyngodon idella).

Galectins are &#x3b2;-galactoside-binding lectins that play essential roles in innate immunity. Among them, tandem-repeat galectins (TrGals), typically composed of two distinct carbohydrate-recognition domains (CRDs) connected by a linker peptide, are well established as key regulators of pathogen recognition and host defense in mammals. However, their structural diversity and immunological functions in teleost fish remain poorly understood. In this study, five TrGals (Gal-4, Gal-8a, Gal-8b, Gal-9, and Gal-9like) were identified in grass carp. Sequence and structural analysis revealed that Gal-8a/b, Gal-9, and Gal-9like possess the canonical two-CRD architecture, whereas Gal-4 uniquely contains four highly similar tandem-repeat domains. All five TrGals were broadly expressed across examined tissues, with predominant expression in the liver. Upon Aeromonas hydrophila infection, Gal-4, Gal-8a, Gal-8b, and Gal-9 were rapidly up-regulated at early time points (3-6&#x202f;h). To elucidate the functional significance of CRD number, recombinant full-length CiGal-4 (CiGal4-full) and three truncated variants containing one, two, or three CRDs (CiGal4-1CRD, CiGal4-2CRD, and CiGal4-3CRD) were generated and systematically characterized. All recombinant proteins contained the conserved &#x3b2;-sheet structure typical of galectin CRDs. Functional assays revealed that CiGal4-full displayed the strongest growth-inhibitory activity against all tested bacteria, whereas CiGal4-1CRD showed the weakest effect. Notably, CiGal4-2CRD exhibited the most potent bactericidal activity, surpassing the full-length protein, while CiGal4-3CRD showed no further enhancement. CiGal4-full and CiGal4-2CRD showed superior carbohydrate-binding activities compared with the other variants. Collectively, these results reveal that CRD copy number alone does not linearly determine galectin function. Instead, domain organization and conformational coordination are critical for optimizing antimicrobial activity. This study provides new insights into the structure-function relationships and evolutionary diversification of galectins in teleosts and highlights their potential as novel antimicrobial and immunomodulatory agents in aquaculture.

Animals

Transcription regulation of cell fate plasticity - from embryonic development to tissue regeneration.

Cell fate plasticity refers to the capacity of cells sharing the same genome to alter, reverse, or reconfigure their identity under physiological, pathological, or experimental conditions. This property underlies embryonic development, cellular reprogramming, and tissue regeneration, but becomes progressively restricted as lineage identity is stabilized. Embryonic development represents an intrinsic process of fate transitions, whereas reprogramming and regeneration reveal how differentiated cells can dedifferentiate or transdifferentiate under specific conditions. Across these contexts, plasticity is governed by multilayered regulatory networks involving transcription factors, epigenetic regulators, cofactors, and the core transcription machinery. Robust regulatory programs stabilize cell identity, whereas stochastic fluctuations in gene expression and chromatin state can prime cells for fate transitions, adding a tunable dimension to plasticity control. In this review, we synthesize recent advances in the regulation of cell fate plasticity across development, reprogramming, and regeneration, highlighting how transcription factors, epigenetic modifications, transcriptional cofactors, and core transcription machinery cooperate to control cell fate decisions and plasticity.

Animals

Maize ZmMYB59 inhibits post-germinative shoot and root elongation through ZmGA2ox3/10-mediated gibberellin catabolism.

Gibberellin (GA) promotes seed germination, but sustained or excessive GA signaling after germination can lead to aberrant root and shoot elongation. How GA homeostasis is transcriptionally restrained during post-germinative seedling development remains unclear. Using overexpression and gene-edited maize materials, we demonstrate that ZmMYB59 inhibits root and shoot elongation during post-germinative growth. Integrated RNA-Seq and CUT&Tag analyses identified the GA catabolism genes ZmGA2ox3 and ZmGA2ox10 as candidate direct targets of ZmMYB59. Hormone profiling analysis showed elevated bioactive GA1 and GA4 levels in the scutellum and aleurone layer cells of zmmyb59 mutants. Dual-luciferase assays, electrophoretic mobility shift assays, and ChIP-qPCR further confirmed that ZmMYB59 directly binds AC8 cis-elements in the ZmGA2ox3/10 promoters and activates their transcription. The zmga2ox3/10 double mutant, but neither single mutant, exhibited enhanced root and shoot elongation, accompanied by GA4 accumulation. This phenotype was suppressed by exogenous application of the GA biosynthesis inhibitor uniconazole. Transcriptomic and biochemical analyses further revealed enhanced starch degradation, reduced starch content, and increased soluble sugar accumulation in the double mutant. Taken together, these findings reveal that the ZmMYB59-ZmGA2ox3/10 module restrains GA accumulation and starch mobilization after germination, thereby coordinating reserve utilization with post-germinative root and shoot growth in maize.

Gibberellins

Fructophilic lactic acid bacteria as a window into multi-scale convergent evolution.

Fructophilic lactic acid bacteria (FLAB) are a group of lactic acid bacteria with unique growth characteristics, that is, poor growth on glucose. Their growth is enhanced in the presence of fructose or external electron acceptors. These organisms inhabit fructose-rich environments such as flowers, fruits, and pollinating insects, particularly honey bees. Apilactobacillus spp. and Fructobacillus spp. are representatives of FLAB, although they belong to phylogenetically distant clades. These organisms commonly possess markedly small genomes with a low number of coding DNA sequences. Furthermore, their genomes are characterized by a markedly reduced number of genes involved in carbohydrate transport and metabolism. Genome reduction in FLAB reflects convergent adaptation to fructose-rich environments rather than general genome streamlining. The two distinct FLAB genera, Fructobacillus and Apilactobacillus, independently lost more than 100 genes in statistically similar orders. In contrast, genes involved in carbohydrate and amino acid metabolism exhibited reversed orders of loss between the two genera. Furthermore, FLAB genomes lack an intact bifunctional alcohol/aldehyde dehydrogenase gene (adhE), which causes their poor growth on glucose. A comparative genomic study suggested the evolutionary process underlying adhE gene decay during adaptation to the fructose-rich environments, including pollinating insects. In conclusion, FLAB represent a unique example of habitat-driven convergent reductive evolution that can be investigated across multiple biological scales - from individual genes to whole genomes - in the diverse LAB group with a wide range of habitats, and partially share the fructophilic evolution with eukaryotic yeasts found in fructose-rich habitats.

Fructose

Identification and characterization of G protein-coupled receptors in the nocturnal halictid bee Megalopta genalis.

G protein-coupled receptors (GPCRs) are one of the largest families of membrane proteins in insects, regulating vision, neural signal transduction, and various physiological behaviors. Megalopta genalis exhibits a unique facultatively eusocial lifestyle and possesses adaptations for nocturnal activity; however, its GPCR family has not yet been systematically characterized. In this study, we performed genome-wide identification, phylogenetic analysis, and expression profiling of GPCRs in M. genalis by integrating genomic annotation and transcriptomic analysis. The results showed that a total of 99 GPCRs were identified in the genome of M. genalis, which were classified into four major families. Here, we show that M. genalis has undergone lineage-specific GPCR repertoire remodeling, marked by the expansion of novel orphan receptors and the systematic loss of multiple receptor subtypes, such as the neuropeptide receptors MIP-R and NPFR. Moreover, opsins have formed a diverse array of combinations and non-GPCR odorant receptors have undergone significant expansion via tandem duplication. Together, these features may represent part of the molecular repertoire associated with the adaptation of M. genalis to a nocturnal lifestyle. Furthermore, transcriptomic analysis revealed distinct spatiotemporal expression divergence within each of the Mth/Mthl and Fz GPCR families, suggesting functional specialization across development and adult tissues. This study provides the first systematic identification and initial functional characterization of GPCRs in M. genalis, revealing an evolutionary pattern characterized by the coexistence of contraction and expansion within the GPCR family. These findings lay a foundation for further studies aimed at elucidating the roles of these GPCRs in regulating M. genalis physiology and behavior.

Animals

Upscaling Genotyping by Amplicon Sequencing With GBAS-GUI.

Genotyping by amplicon sequencing (GBAS) is a relatively low-cost approach for generating genotypic data compared with established genomic methods, making it highly scalable and particularly suitable for large-scale genetic monitoring projects. However, most existing analytical pipelines are either marker-specific, insufficiently scalable, or lacking efficient data management systems for the long-term integration of genotypic information, limiting the full potential of GBAS. Here, we address this gap by introducing GBAS-GUI (https://github.com/sonnenbe-dot/GBAS-GUI), a pipeline capable of generating GBAS-based genotypic data for a wide variety of loci at scale. GBAS-GUI integrates a graphical user interface with multiple checkpoints to improve accessibility and robustness. It implements multiprocessing architecture and a relational database that links genotypic data with associated sample metadata to enhance scalability and data management. The pipeline further enables marker screening through automated calculation of polymorphism information content (PIC) and implements a strategy to recover homologous genotypic information from paralogous loci with non-overlapping amplicon length ranges. Using multiple empirical datasets, we demonstrate substantial improvements in processing speed, database management and handling artefacts related to co-amplification of unspecific regions and duplicates of the same genomic region. We further show that incorporating the full sequence information captured by an amplicon increases marker information content beyond what is achievable with length-based genotyping alone and expands the analytical versatility of GBAS. Overall, GBAS-GUI provides a robust, scalable and versatile framework that unlocks the potential of GBAS for large-scale population genetic and phylogeographic studies.

Genotyping Techniques

Molecular mechanisms of natural de novo shoot organogenesis and their applications.

Natural de novo shoot organogenesis (DNSO) is the spontaneous regeneration of shoots from wound sites outside the shoot apical region through endogenous developmental programs. This regenerative capacity enables plants to recover from severe tissue damage by re-establishing the shoot-root axis. Here, we review current knowledge about the molecular mechanisms of natural DNSO, focusing on transcriptomic and physiological studies in model plants. Accumulating evidence suggests that natural DNSO proceeds through three sequential phases: (i) early wound responses, characterized by the activation of the WIND1-ESR1 module and the establishment of apical-basal auxin asymmetry; (ii) cellular proliferation driven by metabolic and cell-cycle reprogramming; and (iii) cytokinin-mediated establishment of shoot apical meristem identity. We also discuss how these mechanistic insights have been harnessed for practical applications, including tissue culture-free transformation systems such as the cut-dip-budding (CDB) method, and developmental reprogramming strategies that employ ectopic expression of developmental regulator (DR) genes to induce DNSO in otherwise recalcitrant species. Together, these advances illustrate how understanding natural regeneration can guide the development of simplified, broadly applicable plant transformation technologies.

Plant Shoots

Transcriptomics reveals species-specific adaptive strategies to calorie restriction in two Argopecten scallops with distinct lifespans.

Calorie restriction (CR) is a well-established non-genetic intervention for lifespan extension in multiple model organisms. Seasonal food shortage in cold and temperate seas may mimic CR, inducing in bivalves a response similar to that in vertebrates and thereby prolonging life expectancy. However, the relationship and the mechanism underlying the food availability and lifespan in bivalves remain largely unexplored. Two closely related scallop species the short-lived warm-water Argopecten irradians (lifespan <2&#xa0;years) and the longer-lived cold-water Argopecten purpuratus (7-10&#xa0;years) provide an ideal comparative system to investigate species-specific adaptive strategies. In this study, we subjected both species to CR for 30 and 56&#xa0;days and performed comparative transcriptomic profiling, weighted gene co-expression network analysis (WGCNA), and physiological assays to elucidate their distinct molecular responses. Transcriptomic analysis revealed that A. purpuratus exhibited substantially more DEGs than A. irradians at both time points under CR, with both species showing downregulation of metabolic pathways but to different extents. A. irradians mounted an early nutrient-sensing response at 30&#xa0;days (IGF1R, PIK3R3, INSR suppression), indicating acute sensitivity to limitation; by contrast, A. purpuratus displayed delayed FoxO activation at 56&#xa0;days, along with its downstream effectors NFKBIA, CREB3L4, and SMAD4, suggesting a gradual adaptive program may link to its extended lifespan. WGCNA identified three negatively correlated modules in each species, with coral2 being the most prominent in A. irradians and darkolivegreen in A. purpuratus. The former was dominated by ciliary motility genes, whereas the latter featured coordinated repression of oxidative phosphorylation. Additionally, both species exhibited conserved suppression of mTOR/S6K growth signaling and activation of cellular maintenance programs. Collectively, these findings expand the understanding of CR-mediated longevity regulation in bivalves and provide candidate gene resources for future functional studies and breeding programs.

Pectinidae

Exercise-associated epigenetic remodeling and TCR repertoire dynamics in Lynch syndrome carriers.

Lynch syndrome (LS) carriers are at elevated cancer risk. Emerging evidence suggests that exercise may serve as a non-pharmacologic preventive strategy, yet the epigenetic and immunological mechanisms underlying its protective effects in this population remain unclear. Here, we perform integrative multi-omics profiling of DNA methylation, gene expression, and the T cell receptor (TCR) repertoire in LS carriers undergoing a 52-week aerobic cycling intervention. We identify compartment-specific DNA methylation changes, including innate immune activation in cfDNA and oncogenic pathway repression in tissue. Integrative transcriptomic analysis highlights ISL1 as a key exercise-repressed, epigenetically regulated gene, and identifies FLCN as a colorectal cancer (CRC)-associated methylation target. TCR analysis reveals an exercise-associated increase in systemic repertoire diversity and tissue-specific clonal convergence, thus suggesting antigen-driven recruitment. Collectively, these findings uncover epigenetic and immune remodeling as potential mechanisms of exercise-mediated protection in LS.

Lynch syndrome

Insights into the fate and dynamics of antibiotic resistance in multidrug-resistant Bacillus cereus during in vitro simulated gastrointestinal digestion.

Bacillus cereus, an important pathogen responsible for causing foodborne diseases worldwide, releases pore-forming enterotoxins, which target host epithelial cells, leading to osmotic lysis and ultimately manifesting as diarrheal syndrome. Moreover, some B. cereus strains carry antimicrobial resistance genes that confer multidrug resistance against a spectrum of antibiotics. Characterizing the survival traits of multidrug-resistant (MDR) B. cereus strains in the intestinal microenvironment is essential for developing targeted strategies to effectively manage diarrheal foodborne diseases caused by this pathogen. This study used whole-genome sequencing (WGS) to evaluate the pre- and post-digestion toxigenic potential, antimicrobial resistance profiles, and genetic diversity of MDR B. cereus strains isolated from food samples in Guangdong Province, China. The four B. cereus isolates investigated in this study exhibited a genetic diversity, as determined by multilocus sequence typing analysis of WGS data. All four isolates produced the diarrheal toxins Hbl, Nhe, and CytK to varying levels, indicative of their potential to cause outbreaks of foodborne diseases. Each of the four isolates exhibited resistance to more than three classes of antibiotics, fulfilling the criterion for multidrug resistance. At an initial concentration of 9 log colony-forming units (CFU)/mL, the intestinal concentration of these four isolates crossed the threshold required to induce widespread diarrhea in the general population. Under rice slurry protection, all tested isolates maintained intestinal concentration beyond the threshold when the initial concentration was increased to &#x2265;8 log CFU/mL. Moreover, the upregulations of genes associated with acid tolerance, bile tolerance and stress response were observed in the surviving MDR B. cereus isolates. Digestion markedly altered the antibiotic resistance profiles of the MDR B. cereus isolates. In the absence of a food matrix, the MDR isolates lost their resistance to imipenem, meropenem, amoxicillin-clavulanic acid, and trimethoprim-sulfamethoxazole post-digestion and was influenced by the initial concentration of the strains. In the presence of food matrix rice slurry, the effects of digestion on the antibiotic resistance of MDR B. cereus isolates can be mitigated, enabling them to maintain their antibiotic resistance to the greatest extent. Most remarkably, after digestion, the isolates Bce055 and Bce166 exhibited newly emergent resistance to cefotetan and trimethoprim-sulfamethoxazole, respectively. Our findings clarify the fate of MDR B. cereus isolates in the gastrointestinal tract and inform the development of prevention and control strategies for foodborne diseases caused by this pathogen.

Drug Resistance, Multiple, Bacterial

Comparison of paralog identification methods and their impact on species tree topologies in target capture phylogenomics within the Sindora clade (Detarioideae: Leguminosae).

Target capture is a common method of generating high throughput DNA sequencing data for phylogenetic reconstruction of species relationships, for which single copy genes are usually most informative. However, a pervasive problem with target capture is that putatively single copy genes may in fact be paralogs resulting from gene duplication, which are problematic for phylogenetic inference because their evolutionary history may differ from the divergence history of species. Here, we use as a case study a target enrichment dataset of 88 species of Detarioideae (Leguminosae) with a focus on the Sindora clade to examine approaches for handling paralogs, including the built-in paralog handling functions in HybPiper and CAPTUS, plus subsequent steps using Putative Paralog Detection and the tree-based Yang & Smith orthology inference approach. We compare the paralogs flagged using these methods and verify their performance with BLAST mapping against a reference genome sequence of Sindora glabra, and then subsequently compare the species tree topologies produced across these methods. Our comparisons of paralogs flagged across the Sindora clade show that the Putative Paralog Detection pipeline was the most accurate in identifying paralogs in terms of its similarity to the BLAST mapping, followed by the built-in paralog identification function of CAPTUS. However, the results we recovered for the Detarioideae subfamily suggest that the largest differences in species tree topology resulted from the use of paralog-filtered alignments (such as with the Putative Paralog Detection pipeline and the Yang & Smith orthology inference approaches) rather than just by removing the sequences of identified paralogous genes. This was the true for HybPiper-assembled datasets but was not seen in CAPTUS-assembled datasets. In all comparisons, the topological differences caused by different paralog handling methods tended to be confined to clades where processes such as hybridisation and introgression are prevalent. Our study provides a roadmap to establish the best approach to identify, eliminate or separate paralogs in the absence of a chromosomally contiguous reference genome for a study group, and highlights the importance of careful data inspection and processing in addition to understanding the extent of paralogy and paralog characteristics (e.g. sequence divergence between copies) for their study group.

Phylogeny

Genetic evidence for a causal relationship between melatonin metabolism and depression.

To investigate the causal relevance of melatonin metabolism, which provides the biological basis for circulating melatonin levels, to specific depression symptom subtypes, we performed a targeted systematic review of melatonin metabolism pathways in the human brain and liver. Using two-sample Mendelian randomization (MR), we assessed the causal effects of metabolism pathways and/or individual genes on major depressive disorder (MDD) and nine symptom subtypes derived from Patient Health Questionnaire-9 (PHQ-9). Instrumental variables (IVs) were expression quantitative trait loci (eQTL) for eight individual genes, one synthesis route, and three degradation routes. Results were assessed using Bayesian colocalization and phenome-wide association analyses. At the pathway-level, the genetically proxied synthesis-route signal was associated with PHQ-9 Assessment 5 (PHQ9A5, OR: 0.89, 95% CI: 0.85-0.93), but sensitivity analyses suggested this association was primarily driven by TPH1 and may reflect serotonin-related biology. In contrast, higher brain melatonin degradation raised the risk of both PHQ9A1 (OR: 1.03, 95% CI: 1.02-1.04) and PHQ9A7 (OR: 1.03, 95% CI: 1.02-1.03). Within degradation, up-regulation of the kynurenine sub-pathway increased the odds of PHQ9A3 (OR: 1.05, 95% CI: 1.02-1.07), PHQ9A4 (OR&#xa0;=&#xa0;1.04, 95% CI: 1.02-1.06) and PHQ9A7 (OR: 1.05, 95% CI: 1.02-1.07). Gene-level analyses were largely concordant, except for SULT1A1, whose higher expression was genetically protective for PHQ9A3 but risk-increased for PHQ9A1 and PHQ9A4. Overall, these results demonstrate that melatonin metabolism exerts symptom-specific and pathway-specific causal effects on depression. A stratified view of melatonin's role may help optimize the application of exogenous melatonin supplementation.

Melatonin

Hope is linked with more favorable tumor molecular signatures in serous ovarian cancer.

INTRODUCTION: Hope has been associated with improved quality of life and lower mortality in cancer, but the underlying biological mechanisms are poorly characterized. We previously reported that hope was associated with less inflammation and more normalized diurnal cortisol pre-treatment among women with ovarian cancer. We also reported associations of socio-environmental factors with pro-metastatic processes. Here, we used genome-wide transcriptional profiling to quantify associations between hope and tumor molecular signatures reflecting invasiveness, inflammation, and cellular immunity. METHOD: Participants were 74 women with serous ovarian cancer who provided demographic information and completed surveys pre-surgery. Hope was assessed using a face-valid item from the Center for Epidemiological Studies Depression Scale (CES-D). Depression was assessed using the full CES-D without the hope item. Illumina HT12 microarrays were used to assay tumor RNA, and associations between hope and tumor gene expression were quantified, adjusting for depression, age, BMI, grade, and stage. RESULTS: Adjusting for covariates, hope was associated with multiple favorable differences in RNA expression, including lower levels of mesenchymal differentiation (p&#x202f;=&#x202f;0.008) and pro-inflammatory gene regulation (NF-&#x3ba;B: p&#x202f;<&#x202f;0.001; IRF1: p&#x202f;=&#x202f;0.024; STAT: p&#x202f;=&#x202f;0.018), elevated epithelial differentiation (p&#x202f;=&#x202f;0.011), and elevated activity of the IRF7 transcription factor which promotes cellular immunity (p&#x202f;=&#x202f;0.016). CONCLUSIONS: These data suggest that hope is associated with an ovarian tumor gene expression profile characterized by reduced epithelial-mesenchymal transition (EMT) and inflammatory activity, and increased activity of a transcription factor promoting cellular immunity. These findings highlight potential biological implications of a resilience factor such as hope, but need replication with more robust assessments of hope.

Epithelial mesenchymal transition