Search PubMedSearch

SEARCH · Search PubMed

Results for “ChIP-seq”

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.

56 recordsLinked to original sources

ChIP-seq profiling identifies diapause-regulated H3K27me3 targets in the fat body of Culex pipiens.

Culex pipiens, a principal vector of significant arboviruses, survives winter through diapause, a hormonally controlled inactive phase that enhances endurance under severe cold circumstances. Recent data suggests that epigenetic processes, namely histone post-translational modifications (hPTMs), play a crucial role in regulating seasonal dormancy. Prior studies from our laboratory indicated a decrease in the methylation of Histone 3 (H3K27me3) in diapausing fat body tissue, associated with elevated expression of the histone demethylase UTX. Nonetheless, the precise genomic areas impacted by these chromatin alterations remained unidentified. We used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) to delineate the genome-wide distribution of H3K27me3 across fat body chromatin in diapausing (D) and non-diapausing (ND) female Cx. pipiens. Notably, the higher signal at transcription start sites (TSSs) reflects localized redistribution rather than a global decrease, as diapausing fat bodies retain less H3K27me3 overall but concentrate it at promoters. To investigate the functional significance of these chromatin alterations, we confirmed a number of target loci via ChIP-qPCR and assessed gene expression with qRT-PCR. We identified many critical genes that were markedly increased in diapausing mosquitoes, exhibiting an inverse relation to H3K27me3 enrichment. Our data demonstrates different H3K27me3 chromatin landscapes between diapausing and non-diapausing Cx. pipiens, corroborating a hypothesis of selective, locus-specific repression in the non-diapause state and its targeted removal during diapause to permit activation of dormancy-associated genes. These results suggest that chromatin remodeling is a core driver of the diapause switch.

Animals

Clonal haematopoiesis of indeterminate potential and epigenetic age acceleration: Systematic review and meta-analysis.

Clonal haematopoiesis of indeterminate potential (CHIP) represents somatic mutations in haematopoietic stem cells that drive clonal expansion. Epigenetic age acceleration (EAA), estimated from DNA methylation (DNAm) clocks, may capture age-related changes in haematopoiesis. This systematic review and meta-analysis was conducted to synthesise evidence on associations between CHIP and EAA and explore shared biological mechanisms that may underlie this relationship. Six databases were searched from January 1, 2011, to June 6, 2025, adhering to PRISMA 2020. Random-effects meta-analyses were performed. Five studies comprising 7483 individuals (ages 55-79, 67.1% female) assessing associations between CHIP and DNAm clocks were included. Across studies, CHIP individuals had higher EAA than no-CHIP individuals, and larger clones were associated with higher EAA. Meta-analysis of three cross-sectional studies (n = 6946) showed that CHIP had higher EAA versus no-CHIP for Horvath1Age IEAA (mean difference, MD=2.84 years, 95% confidence interval, CI: 1.49-4.19), HannumAge EEAA (MD=2.31 years, 95% CI: 1.14-3.49), PhenoAge (MD=1.84 years, 95% CI: 0.96-2.71), and GrimAge (MD=1.20 years, 95% CI: 0.80-1.61). Both DNMT3A- and TET2-mutated CHIP were associated with higher EAA with TET2-mutated CHIP showing larger effect sizes and more consistent associations than DNMT3A-mutated CHIP across DNAm clocks tested. Higher EAA may also act as an effect modifier for morbidity and mortality in CHIP. Larger longitudinal studies are needed to verify a temporal relationship and determine whether EAA provides incremental prognostic value for morbidity and mortality in CHIP.

Humans

First insights into the role of evolutionary history in shaping venom composition of Vipera ammodytes.

Understanding intraspecific venom variation requires distinguishing the contributions of neutral population history from natural selection. This study aims to determine whether venom variation in Vipera ammodytes species complex is structured across eight phylogenetic groups. Despite a complex evolutionary history, venom composition did not differ among phylogenetic groups within the analytical framework used, suggesting that shared ancestry alone does not explain venom variation. Whether local adaptation to environmental conditions explains the observed variation remains an open question for future studies.

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

Clinical performance of monolithic and veneered zirconia three-unit posterior FDPs: A five-year multicenter randomized controlled trial.

AIM: This randomized controlled clinical study compared monolithic, partially veneered, and fully veneered zirconia FDPs over a 5-year period with respect to survival, technical and biological complications, and patient-reported outcome measures (PROMs). MATERIALS AND METHODS: Sixty-four patients requiring three-unit posterior FDPs were randomly allocated to monolithic (MONO-FDP), partially veneered (PV-FDP), or fully veneered (FV-FDP) groups. All FDPs were fabricated from 4 mol% Y₂O₃ partially stabilized zirconia (4Y-TZP). Follow-up examinations were conducted at baseline, 1, 3, and 5 years. Technical parameters were evaluated using modified USPHS criteria. Periodontal measurements (PPD, BOP, PI) and patient satisfaction were assessed at all time points. RESULTS: A total of 63 FDPs were evaluated at baseline, 57 at 3 years, and 54 at 5 years. Survival at 5 years was 94.7% for MONO-FDPs, 100% for FV-FDPs and 100% for PV-FDPs. Technical complications occurred exclusively in PV-FDP (33%) and FV-FDP (38%) groups and consisted of minor, polishable chipping; no chipping or fractures were recorded in MONO-FDPs (0%) with statistically significant difference between MONO-FDP and the other two groups (p ≤ 0.014). Biological parameters remained stable across all groups, with no significant differences in PPD, BOP, or PI. PV-FDPs and FV-FDPs tended to receive more favorable professional color ratings, whereas PROMs were similar among the groups. CONCLUSIONS: All three FDP designs-monolithic, partially veneered, and fully veneered-fabricated from 4Y-TZP zirconia demonstrated excellent 5-year clinical performance. Technical complications were limited to PV- and FV-FDPs and consisted of minor chipping. Biological outcomes and patient satisfaction were similar across groups. CLINICAL SIGNIFICANCE: The five-year outcomes suggest that veneered, partially veneered, and monolithic FDPs can be used with high clinical reliability; however, restorations incorporating veneering ceramic may present an elevated risk of ceramic chipping.

Humans

ScRNA-seq analysis reveals the effects of nitrite stress on the endocrine system of the eyestalk in Litopenaeus vannamei.

Nitrite is a harmful substance generated in Litopenaeus vannamei farming systems, largely originating from the inadequate breakdown of surplus feed and shrimp feces. Its accumulation in the water can affect the growth and physiological functions of shrimp, damage the immune system, and even cause mass mortality, thus becoming a key environmental factor restricting the green development of the industry. Under nitrite stress, the eyestalk, as an important neuroendocrine regulatory center in crustaceans, participates in the stress adaptation of the organism and exerts a protective effect by regulating energy metabolism and immune function. However, the molecular regulatory mechanism of the eyestalk in response to nitrite stress remains unclear. In this study, single-cell RNA sequencing (scRNA-seq) technology was used to analyze the heterogeneity of eyestalk cells in L. vannamei under nitrite stress. A total of 18, 394 high-quality cells were obtained, and six major cell subpopulations, including Neurosecretory cell, Motor neuron, Sensory neuron, Interneuron, Neurogliocyte, and Support cell, were identified. Differential expression analysis identified 839 differentially expressed genes, and different cell types showed distinct specific responses to nitrite stress. Functional enrichment analysis indicated that pathways such as glycolysis, oxidative phosphorylation, ribosome function, and endoplasmic reticulum protein processing were significantly activated, while signal transduction and DNA repair-related pathways were inhibited. Further analysis revealed that nitrite stress could induce mitochondrial function changes and trigger oxidative stress, thereby affecting the neuroendocrine system function of the eyestalk. This study provided insights into transcriptomic responses of the eyestalk to nitrite stress at the single-cell level, laying a theoretical foundation for the management of aquaculture environments.

Animals

A contextual activity score (CAS) for inferring ADAR-associated transcriptional activity across RNA-seq, single-cell, and spatial transcriptomics.

BACKGROUND AND OBJECTIVE: Adenosine-to-inosine RNA editing, catalyzed by Adenosine Deaminases Acting on RNA (ADARs), is a widespread modification involved in neural function, immune regulation, and cancer. The Alu Editing Index (AEI) is the standard metric to estimate ADAR activity but requires raw sequencing reads and is poorly suited for single-cell and spatial transcriptomic data. This study aimed to develop an alternative framework for inferring ADAR-associated transcriptional activity from gene expression data across diverse transcriptomic technologies. METHODS: We developed the Contextual Activity Score (CAS), a framework based on transcriptional signatures from ADAR perturbation experiments. Context-specific signatures were generated for human neurons, mouse neurons, and cancer models to infer ADAR1 and ADAR2 activity. CAS was computed from normalized gene expression matrices using regulon-based enrichment analysis. Performance was evaluated by comparing with the Alu Editing Index across bulk RNA sequencing datasets, simulated sequencing depths, and library preparation protocols. RESULTS: CAS showed strong concordance with the Alu Editing Index across multiple datasets, while remaining robust to reduced sequencing depth and different library protocols. Unlike the Alu Editing Index, CAS can be applied to single-cell and spatial transcriptomic data and enables the independent assessment of ADAR2 activity. In cancer and neuronal contexts, CAS captured biologically meaningful variations in ADAR-associated transcriptional activity at sample, cell-type, and spatial levels. CONCLUSION: CAS provides a scalable approach applicable across multiple RNA-seq protocols for estimating ADAR-associated transcriptional activity using gene expression data. This method, implemented in an open-source R package for broad adoption, expands the ability to study ADAR-associated transcriptional activity across transcriptomic modalities where direct editing quantification is challenging, such as single-cell and spatial transcriptomics.

Adenosine Deaminase

Functional role and regulatory network of miR-22-3p in chicken hepatic lipid metabolism.

Although microRNA-22-3p (miR-22-3p) is abundantly expressed in the avian liver, its epigenetic role in lipid homeostasis remains largely uncharacterized. To elucidate its in vivo function, 14-day-old female Qingyuan Partridge chickens were intravenously injected with lentiviral vectors to establish miR-22-3p overexpression and knockdown models. Phenotypic analysis demonstrated that miR-22-3p knockdown significantly elevated hepatic triglyceride (TG) levels (p&#xa0;<&#xa0;0.05) and drove marked steatosis, whereas its overexpression reduced TG content. Transcriptome sequencing (RNA-Seq) revealed profound metabolic remodeling, identifying 23 core lipid-associated genes (e.g., ELOVL6, FADS2, ACSBG2, and PTGIS) heavily enriched in steroid biosynthesis, fatty acid metabolism, and elongation pathways. In conclusion, miR-22-3p functions as a bidirectional epigenetic rheostat that negatively regulates hepatic lipid deposition by orchestrating a multilayered polygenic network, providing novel molecular targets for mitigating avian metabolic disorders and optimizing production traits in indigenous poultry breeds.

Animals

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

Transcriptomic analysis provides molecular insights into the innate immune defense of Mactra veneriformis against Vibrio alginolyticus infection.

Mactra veneriformis is an economically important bivalve mollusc in China, but its aquaculture is frequently threatened by Vibrio infections, particularly Vibrio alginolyticus. To investigate the molecular immune response of M. veneriformis to V. alginolyticus, we performed RNA-seq analysis of hepatopancreatic tissues collected at 48&#xa0;h post-infection, the peak mortality time point, with PBS-injected individuals used as controls. Infection with V. alginolyticus caused severe histopathological damage in the hepatopancreas and resulted in a cumulative mortality of 53.3% over 14 d, compared with 3.3% in the control group. Transcriptomic analysis identified 2623 differentially expressed genes (DEGs), including 1585 significantly up-regulated genes and 1038 down-regulated genes. KEGG enrichment analysis demonstrated that DEGs were significantly enriched in immune related and metabolism pathways, including the JAK-STAT signaling pathway, RIG-I-like receptor (RLR) signaling pathway, and cytochrome P450 (CYP450) signaling pathway. Collectively, these findings revealed candidate immune related genes (tlr3, tlr5, myd88, nfkb1, il-17d, and ifi44l), a putative TLR-MyD88-NF-&#x3ba;B signaling axis, and KEGG signaling pathways, including JAK-STAT, RLR and CYP450, that may be involved in the innate immune response of M. veneriformis to V. alginolyticus infection. These results provide a transcriptomic basis for understanding host-pathogen interactions in this species and highlight candidate genes and pathways for future functional validation and potential application in disease-resistance breeding.

Animals

Transcriptomic insights into thermal stress reveal physiological trade-off between thermal stress adaptation and reproductive investment in Spodoptera litura.

Spodoptera litura, a highly polyphagous lepidopteran pest, poses a major threat to agricultural productivity due to its remarkable adaptability to diverse environmental conditions. Although heat stress is known to trigger transcriptional reprogramming in insects, the molecular mechanisms underlying thermal stress responses in S. litura remain poorly understood. In the present study, fourth-instar larvae were exposed to acute heat stress (44&#xa0;&#xb0;C) and compared with control conditions (27&#xa0;&#xb1;&#xa0;1&#xa0;&#xb0;C) to investigate heat-induced transcriptional alterations affecting physiology and reproduction. High-quality RNA-Seq data achieved more than 80% mapping efficiency, with a total of 15,782 transcripts were identified. Transcriptome analysis of S. litura larvae showed 323 differentially expressed genes (DEGs), of which 262 genes were significantly upregulated and 61 were downregulated in heat-stressed larvae compared to the control group. The DEGs were associated with stress response, reproduction, signalling, proteostasis, detoxification, oxidative stress, metabolism, development, and chromatin regulation. Heat shock proteins genes, including HSP70, HSP90, and HSP27, together with co-chaperones such as TRET-1, STIP1, and Starvin, were strongly upregulated, indicating enhanced cellular protection against protein damage and oxidative stress under heat stress. Conversely, key reproductive and cell cycle-related genes, including BARR, CAPD2, FEO, CDK2 and MORULA, were significantly downregulated, suggesting reproductive impairment and developmental arrest. RT-qPCR validation corroborated the RNA-Seq findings, demonstrating a heat-induced physiological trade-off that prioritizes survival over reproduction. Consistent with these molecular responses, heat-stressed insects exhibited marked reproductive impairment, including significant reductions in gonadosomatic index, eupyrene sperm bundle count, mating frequency, mating success, female calling behaviour, copulation duration, fecundity, and egg fertility. Collectively, these findings provide comprehensive insights into the molecular basis of thermal adaptation in S. litura and demonstrate that acute heat stress compromises reproductive fitness while activating conserved stress-response pathways that promote short-term survival.

Animals

Single-cell RNA sequencing provides further insights into the immunostimulatory action of freeze-dried Lactiplantibacillus plantarum on Penaeus vannamei shrimp.

Immunostimulation through dietary interventions opened new avenues in developing disease control and prevention tools for shrimp aquaculture. We have previously shown that feeding with freeze-dried Lactiplantibacillus plantarum (LAB) increased disease resistance of Penaeus vannamei against both Vibrio parahaemolyticus and white spot syndrome virus (WSSV) based on bulk RNA sequencing of shrimp gills. This tissue participates in ion transport and serves as a first line of defense against environmental stressors and pathogenic infections. However, characterization of their cell composition and functions remains limited. Here, we implemented a single-cell RNA sequencing approach to further gather insights into how feeding with freeze-dried LAB modulates host immunity which may not be evident with bulk RNA sequencing approach. A total of five clusters with unique transcriptional signatures were identified, corresponding to pillar cells, septal cells, and sessile hemocytes. Pseudo-bulk analyses at global- and cluster-levels showed differential expression of genes related to host immunity and metabolism. We further revealed how overall transcriptomic changes are not exclusively caused by gene expression changes but may also be driven by cell population dynamics. This study highlighted how single-cell RNA sequencing approach may shed light on the mechanisms of action of immunostimulants which may be masked in bulk transcriptome analyses.

Animals

Integrative analysis of transcriptome and DNA methylome dynamics during caudal fin regeneration in silver pomfret (Pampus argenteus).

Caudal fin regeneration in teleost fish is a complex, multi-stage process involving coordinated molecular and cellular changes. While the role of epigenetic regulation particularly DNA methylation has been studied in model freshwater species such as zebrafish, its contribution to regeneration in marine teleosts remains largely unexplored. In this study, we integrated transcriptomic and DNA methylomic data to characterize the temporal dynamics of gene expression and methylation during caudal fin regeneration in the silver pomfret (Pampus argenteus). Using RNA-sequencing and reduced representation bisulfite sequencing (RRBS) at three biologically critical time points 1, 3, and 7&#xa0;days post-amputation (dpa), we characterized the spatiotemporal molecular landscape of caudal fin regeneration. These time points capture the key transitional phases of wound healing and inflammation (1 dpa), blastema formation and progenitor proliferation (3 dpa), and regenerative outgrowth with tissue remodeling (7 dpa), enabling robust detection of the major molecular programs underlying epimorphic regeneration. Concurrently, CG-methylome analysis identified thousands of dynamically changing differentially methylated regions (DMRs). A strong global inverse correlation was observed between promoter methylation and gene expression. Integrative analysis pinpointed key regeneration genes (fgf20a, msxb, sox9b) whose expression was associated with dynamic methylation changes in their promoters or gene bodies. We conclude that DNA methylation is a dynamic and key regulatory layer that acts in concert with transcriptional reprogramming to coordinate tissue regeneration, providing new insights into the epigenetic mechanisms underlying complex regenerative processes in teleosts.

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

Transcriptomic insights into temperature regulation of proliferation and differentiation in skeletal muscle cells of Nibea albiflora.

Myogenesis involves sequential stages of muscle satellite cell activation, myoblast proliferation, differentiation, and fusion into multinucleated myotubes. Teleost muscle exhibits indeterminate growth and is highly sensitive to environmental temperature, yet the underlying mechanisms by which temperature regulate proliferation and differentiation remain poorly understood. In this study, we established a primary skeletal muscle cell culture from the yellow drum (Nibea albiflora), an economically important marine fish, and integrated morphological observations with comparative transcriptomics analysis to characterize cellular and molecular responses at 28&#xa0;&#xb0;C and 20&#xa0;&#xb0;C during both proliferation and differentiation stages. Phenotypic analysis revealed that 28&#xa0;&#xb0;C significantly enhanced both myoblast proliferation and myogenic differentiation ability compared with 20&#xa0;&#xb0;C. Transcriptomic profiling revealed that at 28&#xa0;&#xb0;C, differentiation upregulated extracellular matrix(ECM) organization, calcium signaling, and sarcomere assembly, while proliferation enhanced focal adhesion, growth factor signaling, and lipid metabolism. At 20&#xa0;&#xb0;C, differentiation was characterized by glutathione metabolism, and ferroptosis, while proliferation involved cytokine-cytokine receptor interaction and negative regulation of signal transduction. Core myogenic regulatory factors (MRFs), particularly myogenin, were strongly upregulated at 28&#xa0;&#xb0;C during the differentiation stage, serving as an internal control. Based on these findings, we propose a testable model that elevated temperature coordinates Ca2+-dependent MRF activation with ECM-integrin signaling to drive sarcomere assembly and muscle growth. Key differentially expressed genes (DEGs) regulating myogenesis in N. albiflora were also identified. This study provides a mechanistic framework for temperature adaptation in teleost skeletal muscle and identifies candidate genes for functional validation and marker-assisted selection, as well as a rationale for temperature management strategies to improve aquaculture yield of N. albiflora.

Animals

In vivo genome-wide CRISPR screens identify FOXR1 as a suppressor of CD8+ T cell antitumor immunity.

T cell dysfunction critically limits the efficacy of T cell-based immunotherapies in solid tumors, yet the intrinsic regulators of T cell dysfunction remain incompletely understood. Through an in vivo genome-wide CRISPR screen in tumor-infiltrating CD8+ T cells, we identified Forkhead Box R1 (FOXR1) as a potent transcriptional suppressor of CD8+ T cell effector functions. Genetic ablation of FOXR1 significantly enhanced cytokine production and cytotoxic capacity in both murine and human CD8+ T cells, whereas its overexpression impaired T cell activation and effector molecule expression. Mechanistically, multiomics integration of RNA-seq, CUT&Tag-seq, and ATAC-seq revealed that FOXR1 binds directly to promoter regions of key effector genes, including IL2, GZMB, and PRF1, and represses their expression. Importantly, FOXR1 deletion in human anti-CD19 CAR T cells improved their efficacy against solid tumors, demonstrating that FOXR1 is a checkpoint of T cell effector function and targeting FOXR1 is a promising strategy to enhance CAR T cell efficacy against solid tumors.

Animals

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