SEARCH · Search PubMed

Results for “InDel markers”

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.

SourceResearch connectionFrom yearThrough year

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

17 recordsLinked to original sources

PubMed2026

Development of Genome-Derived InDel Markers and Genetic Diversity Analysis of Caragana acanthophylla in Xinjiang, China.

Caragana acanthophylla Kom. is an ecologically important drought-tolerant shrub in Xinjiang, China, but species-specific molecular markers for germplasm characterization remain limited. We sampled 93 individuals from 11 localities representing the currently known distribution of C. acanthophylla in Xinjiang. Three individuals per locality (33 in total) were whole-genome resequenced, yielding 2,873,410 high-quality SNPs and 5,679,915 InDels. Genome-wide SNP-based PCA and genetic relationship analysis provided an independent high-resolution assessment of the 33 resequenced individuals. From 34 candidate primer pairs, eight polymorphic InDel markers with stable amplification and clear genotyping profiles were retained and applied to all 93 individuals. The SNP dataset revealed clear regional differentiation and finer locality-associated relationships. Analysis of the same 33 individuals with the eight InDel loci recovered part of this broad pattern, particularly the differentiation of the western YL materials, but showed lower fine-scale resolution. Across all 93 individuals, the InDel panel revealed moderate to low marker-level genetic diversity and detectable regional differentiation. AMOVA attributed 67.00% of the variation to differences among the 11 original sampling localities, while the five exploratory analytical groups showed a similar among-group component (68.37%). The Mantel correlation detected across all 93 individuals (r = 0.801, p < 0.001) disappeared after YL was excluded (r = -0.032, p = 0.724), indicating that the overall spatial signal was largely driven by the geographic separation of YL. These results support the eight-marker panel as a practical, low-cost tool for preliminary germplasm characterization and broader sample screening, while genome-wide SNP data provide substantially greater resolution for population-level inference.

Caragana acanthophylla↗

PubMed2026

Pilot study of allele-specific multi-InDel markers for the detection of extremely unbalanced DNA mixtures.

Mixtures are common in forensic casework, and they represent one of the most challenging types of biological evidence. Traditional short tandem repeat analyses are often associated with limitations when dealing with extremely unbalanced mixtures because alleles from minor contributors can easily be masked by those of major contributors. Consequently, researchers have developed new technologies and methods for improving the analysis of mixtures, spanning upstream DNA extraction and downstream software analysis. Among these, strategies combining allele-specific amplification with compound markers have drawn particular interest because of their ability to selectively detect minor contributors in complex mixtures. In this study, we screened multi-InDels across the entire genome, designed allele-specific primers compatible with the capillary electrophoresis platform, and further explored their potential in unbalanced DNA mixtures and cell-free fetal DNA (cffDNA). Ultimately, a set comprising 10 multi-InDels was developed, and this included two groups of primers that separately amplified the long alleles (L primer set) and short alleles (S primer set). The results demonstrated that each primer pair could detect the minor component at a 1:1000 mixture ratio, whereas the L and S primer sets successfully detected the minor contributors at mixture ratios of 1:200 and 1:500, respectively. Furthermore, in the cffDNA analysis, 60 of 78 informative markers were successfully detected, with the complete detection of all informative markers achieved in 18 mother-child reference pairs. Overall, allele-specific amplification-based multi-InDel markers enabled the sensitive detection of minor contributors, providing a potential strategy for the analysis of unbalanced two-person mixtures.

Allelic-specific amplification↗

PubMed2026

Genome-Wide Identification of SSR and InDel Markers and Experimental Validation of SSR Markers for Distinguishing Cold-Tolerant and Cold-Sensitive Lily Cultivars.

In this study, whole-genome resequencing was performed on the cold-tolerant variety ND-6 and the cold-sensitive variety 'Sorbonne'. After evaluation, the Lilium davidii var. unicolor reference genome was selected to analyze SSR distribution characteristics. Whole-genome InDel identification and comparative analysis were conducted for the two varieties, yielding 34,812,909 and 24,497,857 InDels, respectively. Short InDels were predominant, with deletions slightly outnumbering insertions, mostly located in intergenic regions. Twenty pairs of SSR primers were screened and synthesized. Among them, 10 pairs amplified clearly, with a polymorphism rate of 82.6%, effectively distinguishing the two cultivars examined in this study. This study provides systematic data and a reliable marker resource for the analysis of lily genomic variation, laying a foundation for the identification of cold-tolerant germplasm; validation across additional cultivars and individuals will be required to extend their utility to broader germplasm.

cold resistant lilies↗

PubMed2026

Molecular dissection and functional characterization of the liguleless1 gene for manipulation of leaf angle in maize.

Recessive liguleless1 (lg1) gene significantly reduces leaf angle in maize and has become the choice in breeding for high plant density. Here, we sequenced the entire lg1 gene (5560 bp) among seven wild-type (Lg1) and one mutant (lg1) inbreds. The analysis revealed a total of 229 SNPs and 155 InDels within the Lg1 gene. The study also revealed the existence of three exons, with lg1-mutant having two exons. The lg1-mutant harboured an insertion of 130 bp Tourist MITE transposable element in exon-2 at 1663rd base, which deleted 157 amino acids of C-terminal region of the mutant LG1 protein. The mutant LG1 protein was 247 amino acids in length, in contrast to 399-404 amino acids in wild-type protein. The analysis with 26 paralogues and 66 orthologues of Lg1 revealed conservation of the squamosa promoter-binding (SBP) domain. A PCR-based co-dominant InDel marker (MGU-lg1-Tourist) specific to insertion of 130 bp was developed that differentiated the mutant allele (lg1) from the wild-type allele (Lg1). The marker was validated in two F2 populations, which showed a 1:2:1 ratio. F2 plants showed a 3 (wide angle: 41.58°) :1 (narrow angle: 5.88°) segregation for leaf angle. A set of 11 gene-based InDel markers (MGU-InDel1 to MGU-InDel11) specific to Lg1 was also developed, and along with MGU-lg1-Tourist, they classified a diverse set of 48 inbreds into 35 distinct haplotypes (hap1 to hap35) with lg1-based inbreds possessing hap1. This is the first report of the development and validation of a co-dominant gene-based marker specific to lg1, and information generated assumes great significance in maize breeding aimed to tailor the plant architecture suitable for high plant density..

Zea mays↗

PubMed2026

Development of molecular markers associated with saline-alkali tolerance in rapeseed (Brassica napus L.).

A total of 947 saline-alkaline tolerance-related molecular markers and a 5K cGPS genotyping chipwere developed, providing practical tools for marker-assisted selection and molecular design breeding of saline-alkaline-tolerant rapeseed. Rapeseed (Brassica napus L.) has relatively strong tolerance to saline-alkaline stress and shows great potential for the sustainable utilization and improvement of saline-alkaline soils. However, the breeding of highly tolerant cultivars still mainly depends on conventional hybridization combined with phenotype-based selection, which constrains breeding efficiency. In this study, previously reported saline-alkaline tolerance-related genes from rapeseed, rice, maize, wheat, sorghum, and Arabidopsis were collected. Candidate gene-based association analysis enabled the development of molecular markers and a genotyping chip. A total of 483 significantly associated genes were identified, among which 355 genes contained favorable haplotypes. Molecular markers were successfully developed for 275 genes, including 746 KASP and 201 InDel marker pairs, and four marker pairs were randomly selected for validation. In addition, a 5K cGPS liquid-phase chip (HZSW-cGPS-BRNAP-04), was developed and showed a high call rate and excellent reproducibility in genotyping. These markers and the chip are expected to improve the breeding efficiency of saline-alkaline-tolerant rapeseed cultivars. Overall, this study provides useful tools for early-generation evaluation and marker-assisted selection (MAS), and provides a foundation for molecular design breeding of saline-alkali-tolerant rapeseed.

Brassica napus↗

PubMed2025

Combining QTL mapping and RNA-Seq reveals candidate genes controlling flag leaf width in foxtail millet.

BACKGROUND: The flag leaf, a crucial component of plant architecture, significantly influences final grain yield in crops, including foxtail millet (Setaria italica L.). Optimizing flag leaf size is considered an effective strategy for enhancing grain yield potential under higher planting densities. However, the genetic mechanism underlying flag leaf size, particularly flag leaf width (FLW), remains largely unknown under varying planting densities in foxtail millet. RESULTS: An FLW phenotype variation analysis was conducted across multiple planting densities using a recombinant inbred line (RIL) population derived from Heizhigu (narrow leaf) and Changnong 35 (wide leaf). Based on a high-density genetic map with 3795 Bin markers, 11 flag leaf width (FLW) QTLs were identified on chromosomes 3, 5, and 6, explaining 2.35%-36.06%. Among these, qFLW5-2 was a major QTL, detected consistently across 3 environments and explaining a large proportion of FLW variation. The QTL was further validated with 9 InDel markers with its candidate region across different planting densities. Moreover, RNA-seq revealed 2,293 and 2,338 differentially expressed genes (DEGs) between biparents at heading stage and grain filling stage, respectively. There were 11 and 9 DEGs within the location range of qFLW5-2 among 2 comparison groups (HZG-H_vs_CN35-H and HZG-G_vs_CN35-G). Combining QTL mapping and RNA-seq, we speculated that Seita.5g134600 (encoding an auxin responsive protein Aux/IAA) and Seita.5G123900 (encoding a cytochrome P450 family protein) as key candidate genes for qFLW5-2. Furthermore, variation analysis confirmed that the lines or germplasm with Seita.5G1346005UTR277+ allele, both within the RIL population and natural populations, exhibited significantly wider leaves than those with Seita.5G1346005UTR277- allele. These findings advance our understanding of the genetic and molecular regulatory mechanisms governing flag leaf growth. CONCLUSIONS: This study elucidates genetic and molecular mechanism regulating flag leaf growth and development in foxtail millet. The results provide a theoretical foundation for improving plant architecture and facilitating molecular marker-assisted breeding in this crop.

Quantitative Trait Loci↗

PubMed2026

Validation of the TransplantTrace cfDNA Kidney assay for measurement of donor-derived cell-free DNA in transplant recipients.

INTRODUCTION: Donor-derived cell-free DNA (dd-cfDNA) has emerged as a promising non-invasive marker for assessing allograft status and guiding clinical management in transplant recipients. Its utility in kidney transplantation has repeatedly been demonstrated in large studies showing a strong association between elevated dd-cfDNA levels and allograft injury or rejection. This study evaluated the performance of a centralized next-generation sequencing (NGS)-based assay for measurement of dd-cfDNA in patients post-kidney transplantation. METHODS: The TransplantTrace cfDNA Kidney assay utilizes 50 insertion-deletion (indel) markers to discriminate dd-cfDNA. Evaluation of analytical performance included determination of input requirements, analytical sensitivity and specificity, as well as accuracy and precision parameters. Diagnostic performance was evaluated in a retrospective cohort of 104 post-transplantation samples by comparing dd-cfDNA results with biopsy-confirmed rejection. RESULTS: The assay required low DNA input (2 ng) and demonstrated high analytical sensitivity, with a verified limit of detection of 0.2% and limit of quantification of 0.3% dd-cfDNA. Analytical accuracy was excellent (R 2 = 1.00), with high repeatability and reproducibility across the reportable range of 0.2-30% dd-cfDNA. In the clinical validation, the assay showed high concordance with biopsy-confirmed rejection and excellent discriminatory performance for differentiating active from non-active rejection (AUC of 0.980). At a 1% cut-off, the assay exhibited a positive predictive value of 100%, supporting confident identification of patients likely to have treatable graft injury, while a high negative predictive value (95.6% at 15% prevalence) supports its reliability in ruling out active rejection. CONCLUSION: The TransplantTrace cfDNA Kidney assay demonstrated robust analytical performance and strong clinical concordance with biopsy-confirmed rejection status. Its high diagnostic accuracy supports reliable identification and exclusion of active rejection, with the potential to reduce reliance on invasive biopsy procedures in patients with elevated serum creatinine but low dd-cfDNA levels. In summary, the findings of this study support the implementation and use of this centralized assay for measurement of dd-cfDNA in patients post-kidney transplantation.

centralized↗

PubMed2021

Alien introgression and morpho-agronomic characterization of diploid progenies of Solanum lycopersicoides monosomic alien addition lines (MAALs) toward pre-breeding applications in tomato (S. lycopersicum).

Alien introgressions that were captured in the genome of diploid plants segregating from progenies of monosomic alien addition lines of S. lycopersicoides confer novel phenotypes with commercial and agronomic value in tomato breeding. Solanum lycopersicoides is a wild relative of tomato with a natural adaptation to a wide array of biotic and abiotic challenges. In this study, we identified and characterized diploid plants segregating from the progenies of monosomic alien addition lines (MAALs) of S. lycopersicoides to establish their potential as donors in breeding for target trait improvement in tomato. Molecular genotyping identified 28 of 38 MAAL progenies having the complete chromosome complement of the cultivated tomato parent and limited chromosome introgressions from the wild S. lycopersicoides parent. Analysis of SSR and indel marker profiles identified 34 unique alien introgressions in the 28 MAAL-derived introgression lines (MDILs) in the genetic background of tomato. Conserved patterns of alien introgressions were detected among sibs of MDILs 2, 3, 4 and 8. Across MDILs, a degree of preferential transmission of specific chromosome segments was also observed. Morphologically, the MDILs closely resembled the cultivated tomato more than S. lycopersicoides. The appearance of novel phenotypes in the MDILs that are lacking in the cultivated parent or the source MAALs indicates the capture of novel genetic variation by the diploid introgression lines that can add commercial and agronomic value to tomato. In particular, screening of representative MDILs for drought tolerance at the vegetative stage identified MDIL 2 and MDIL 11III as drought tolerant based on visual scoring. A regulated increase in stomatal conductance of MDIL 2 under drought stress indicates better water use efficiency that allowed it to survive for 7 days under 0% moisture level.

Chromosomes, Plant↗

PubMed2025

Cr3a, a candidate gene conferring fruit cracking resistance, was fine-mapped in an introgression line of Solanum lycopersicum L.

In the cultivation and production of tomato (Solanum lycopersicum L.), fruit cracking is a prevalent and detrimental issue that significantly impacts the esthetic quality and commercial value of the fruit. The complexity of the trait has resulted in a slow advancement in research aimed at identifying genes that influence tomato fruit cracking and the underlying regulatory mechanisms. In this study, a sub-introgression population for tomato crack-resistant fruit has been constructed from the cross between S. lycopersicum 1052 and Solanum pennellii LA0716, followed by 11 generations of selfing. Utilizing specifically designed InDel markers, the tomato crack-resistant gene, Cr3a, was fine-mapped, cloned, and its functionality was confirmed through transgenic and gene-knockout approaches. The precise localization of Cr3a was delineated to a 30 kb genomic region on chromosome 3, corresponding to the gene Sopen03g034650 in S. pennellii and Solyc03g115660.3 in the Heinz1706 variety. An integrated transcriptomic and metabolomic analysis of fruits with and without the Cr3a gene was finally conducted to elucidate the intricate regulatory mechanisms associated with Cr3a. The findings revealed a molecular regulatory network for tomato fruit crack resistance, characterized by 7 key metabolites, 13 pivotal genes, and 4 critical pathways: the phenylpropanoid biosynthesis pathway, the phenylalanine, tyrosine, and tryptophan biosynthesis pathway, the linolenic acid metabolism pathway, and the cysteine and methionine metabolism pathway. In summary, this research provides novel insights into the molecular underpinnings of tomato fruit crack resistance and holds substantial promise for accelerating the molecular breeding of tomatoes with enhanced fruit crack resistance.

Solanum lycopersicum↗

PubMed2025

Indel mutation in transcription factor PabHLH2 regulates amygdalin accumulation and kernel bitterness in apricot.

Amygdalin, the phytochemical responsible for the characteristic bitterness of apricot (Prunus armeniaca L.) kernels, also exhibits significant bioactive properties and therapeutic potential. Genetic regulation of amygdalin content is therefore a key objective in apricot breeding programs aimed at quality improvement. In this study, we conducted quantitative trait loci (QTL) mapping to uncover the genetic basis of sweet-bitter differentiation in apricot kernels. We identified a 15-bp insertion/deletion (indel) polymorphism strongly related to kernel bitterness, with marker validation achieving 100% concordance across 601 apricot germplasm accessions. Notably, this polymorphic site is located within the helix-loop-helix (HLH) domain of the basic HLH (bHLH) transcription factor PabHLH2. Protein interaction analyses revealed that the 15-bp deletion variant impaired dimerization capacity, reducing transcriptional activation of downstream targets. Using yeast one-hybrid screening and dual-luciferase reporter assays, we identified PaCYP71AN24 and PaCYP79D16 as direct transcriptional targets of PabHLH2. Functional characterization further indicated that the PabHLH2a variant (harboring the 15-bp insertion) significantly enhanced the promoter activity of these cytochrome P450 genes compared with the deletion variant. Transient overexpression and silencing experiments in apricot kernels further confirmed that the 15-bp insertion positively regulates both PaCYP71AN24/PaCYP79D16 expression and prunasin accumulation, the immediate biosynthetic precursor of amygdalin. Overall, these findings provide mechanistic insights into the allelic variation underlying kernel bitterness and delineate the molecular cascade of amygdalin biosynthesis. The identified molecular markers and functional characterization establish a basis for marker-assisted breeding of low-amygdalin apricot cultivars, supporting the dual-purpose utilization of kernels in food and pharmaceutical industries.

Amygdalin↗

PubMed2026

Integrating GWAS and Transcriptome Analysis Identifies Candidate Genes for Kernel Starch Quality Traits in Maize.

Maize (Zea mays L.) starch quality is a complex trait with significant implications for grain processing and industrial applications. However, the genetic basis underlying starch quality, particularly for gelatinization and thermodynamic properties, remains poorly understood. In this study, we evaluated 12 starch quality traits, including seven gelatinization characteristics, four thermodynamic traits, and kernel starch content (KSC) in a diverse panel of 335 maize inbred lines. Considerable phenotypic variation was observed for all traits. A total of 228 quantitative trait loci (QTLs) were significantly associated with 12 starch quality traits through genome-wide association studies (GWAS). By integrating a dynamic transcriptome analysis of two maize inbred lines with contrasting starch quality, we identified 60 candidate genes. One gene, waxy1, encoding a starch synthase, was found to be associated with enthalpy of gelatinization (ΔHgel) and pasting temperature (Ptemp). Six variants in waxy1 contributed to natural variation in ΔHgel and Ptemp, and a cost-effective InDel and two PARMS-based molecular markers were developed and validated in 144 maize inbred lines, enabling efficient marker-assisted selection. Our findings provide key genes and molecular markers for high-quality maize breeding with improved starch properties.

Zea mays↗

PubMed2026

Genome-wide association studies of plant traits and functional analysis of leaf development-related genes in citrus.

Labor-saving and high-light-efficiency tree architecture is a key breeding objective for woody fruit trees like citrus. However, population genetics information on these traits remains limited. In this study, tree architecture, thorn, and leaf traits were evaluated in 353 F2 progeny derived from a cross between Clementine mandarin and precocious trifoliate orange-an early-flowering variety. A random subset of 300 offspring was sequenced for a genome-wide association study (GWAS), which detected 10 216 significantly associated SNPs and defined several major quantitative trait loci (QTLs) for the target traits. Subsequent bulked segregant analysis (BSA) and GWAS on individuals with extreme compound leaf phenotypes mapped the causal gene(s) to a 0.8 Mb region (22.15-22.95 Mb) on chromosome 4. Genetic analysis across multiple hybrid combinations confirmed that the compound leaf trait in trifoliate orange is dominantly inherited and follows Mendelian segregation. Transcriptome profiling of parental leaves at different developmental stages identified a KNOX gene, CiKNAT6, as a candidate. Further validation using CAPS markers and Hi-Tom sequencing demonstrated tight linkage between an InDel polymorphism in CiKNAT6 and leaf shape across diverse citrus species and the F2 population, with co-segregation observed for the compound leaf trait. Due to alternative splicing producing seven splice variants, the CiKNAT6 DNA sequence was selected for genetic transformation experiments. Functional analysis revealed that the Clementine mandarin allele of CiKNAT6 is non-functional owing to an InDel, whereas ectopic expression of the trifoliate orange allele in tobacco and lemon induced leaf curling and reduced leaf size. CRISPR-Cas9 knockout of CiKNAT6 in trifoliate orange resulted in increased leaf area. These findings provide valuable genetic resources and insights for future studies on tree architecture and leaf morphology.

Plant Leaves↗

PubMed2025

Development of recombinant inbred lines and QTL analysis of plant height and fruit shape-related traits in Cucurbita pepo L.

UNLABELLED: Zucchini (Cucurbita pepo subsp. pepo) stands as an economically vital crop in China. In zucchini breeding, plant architectural patterns and fruit morphological characteristics serve as pivotal traits. In this study, we employed quantitative trait locus (QTL) analysis using recombinant inbred lines (RILs) derived from two distinct inbred lines, JinGL (subsp. ovifera) and HM-S2 (subsp. pepo), in conjunction with a high-density genetic map. Our investigation focused on ten QTLs associated with six horticulturally significant traits, including hypocotyl length (HL), plant height (PH), and four fruit-related traits: fruit length (FL), fruit diameter (FD), fruit shape index (FSI), and fruit weight (FW). The QTLs governing HL and PH were mapped to Chr03/LG10 and named qhl3.1 and qph3.1, respectively. The candidate gene Cp4.1LG10g05910/CpDw for qph3.1 was successfully identified. Additionally, three novel QTLs related to fruit size and shape were discovered. Among them, qfsi8.1/qfl8.1, demarcated by Marker238258 and Marker240069 on Chromosome 08/Linkage group 17 (Chr08/LG17), is a new major QTL regulating the fruit shape of zucchini. Through genomic insertion-deletion (InDel) and qRT-PCR analyses, we predicted genes within the qfsi8.1/qfl8.1 candidate interval, uncovering Cp4.1LG17g02030/CpIAA12 and Cp4.1LG17g02010/CpCalB as potential candidate genes. We developed molecular markers tightly linked to qph3.1 and qfl8.1 and validated them in 171 and 224 Cucurbita pepo germplasms, achieving accuracy rates of 96% and 100%, respectively. This study deepens our understanding of the genetic basis of key traits and provides valuable references for molecular breeding in Cucurbita pepo. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11032-025-01592-y.

Cucurbita pepo↗

PubMed2026

Genetic diversity, phylogenetic relationships, and marker development between Hydrangea serrata and H. macrophylla based on plastome and 45S nrDNA.

Ornamental hydrangeas (genus Hydrangea) are cultivated worldwide for their diverse flower colors and attractive morphology. Here, we assembled the complete plastid genome (plastome) and 45S nuclear ribosomal DNA (45S nrDNA) sequences of 22 individuals representing H. serrata, H. macrophylla, and related species (H. arborescens, H. paniculata, H. petiolaris, and H. hydrangeoides). The plastomes contained up to 2,344 single-nucleotide polymorphisms (SNPs) and 367 insertions/deletions (InDels) within the genus, whereas the assembled 45S nrDNA sequences showed 119 SNPs and 10 InDels. Phylogenetic analyses based on plastome and 45S nrDNA sequences clearly separated H. serrata and H. macrophylla from the other Hydrangea species. In the plastome-based tree, H. petiolaris was placed in the same clade as H. arborescens, whereas in the 45S nrDNA-based tree it showed a close relationship to H. hydrangeoides. The H. serrata and H. macrophylla samples were not always separated according to their species boundaries, as observed in samples Hse8-Hse12. Notably, one H. serrata sample (Hse8), collected from a wild mountainous region of Japan, exhibited a closer genetic relationship to H. macrophylla samples, indicating that cultivated hydrangeas may have originated from a specific wild lineage of H. serrata adapted to mountainous habitats. Using plastome-derived molecular markers, 66 Hydrangea samples were further classified into five groups, with Group II comprising both cultivated H. macrophylla and a subset of wild H. serrata samples, suggesting a close genetic affinity between this group and the ancestral gene pool of cultivated H. macrophylla. Based on these genomic resources, eight plastome-derived molecular markers were developed to differentiate cultivated hydrangeas from wild genotypes and to assess genetic diversity within H. serrata and H. macrophylla, providing practical tools for germplasm identification, breeding, and genetic resource management of Hydrangea species.

hydrangea↗

PubMed2025

Genome-scale multi-organ analysis of mutagenic effects of ethanol and acetaldehyde in Sprague Dawley rats.

Alcohol consumption is a major cancer risk factor, particularly for head and neck cancers, including the oral cavity. Acetaldehyde, the primary genotoxic metabolite of ethanol, may play key roles in oral carcinogenesis, though the mechanisms remain unclear. While mutational signatures SBS16, DBS4, and ID11 have been tentatively linked to alcohol use, they are not exclusive to alcohol-related cancers. In this study, we examined the genome-wide in vivo mutagenic effects of ethanol and acetaldehyde by analyzing tumors from the cheek, Zymbal gland, larynx, forestomach, and liver of rats chronically exposed to these compounds. Signature analysis revealed exposure-specific, early-onset formation of SBS17 in ~28% of head and neck tumors, suggesting inflammation and/or oxidative damage as potential mediators of carcinogenesis. Cancer driver gene analysis identified a relative enrichment of exposed tumors with mutations in the Tp53 and Mtor genes. No notable exposure-specific changes were observed in doublet-base substitutions, indel signatures, or copy number variants. Notably, SBS16, DBS4, and ID11 were absent. Our findings suggest direct mutagenicity may not be the main driver of alcohol-related cancer. Other harmful cellular effects, undetectable by whole genome sequencing, may be involved. Our findings suggest that SBS17 could function as a potential exposure-specific molecular marker of alcohol-related cancers in humans.

Journal Article↗

PubMed2026

Bridging the gap between legacy polymerase chain reaction-based microsatellite data with high-throughput sequencing data for conservation genomics.

Microsatellites are powerful markers for tracking genetic variation in wildlife populations due to their high polymorphism and genome-wide abundance. While polymerase chain reaction (PCR)-based fragment size analysis has been the standard for genotyping microsatellites, high-throughput sequencing offers greater resolution and the opportunity to sync historical datasets with modern analyses. We evaluated how genotypes from whole-genome sequencing align with PCR data for 15 microsatellite loci in 11 North American brown bears (Ursus arctos). Brown bear populations in the 48 contiguous United States have declined from approximately 50,000 to fewer than 2,000 over the past decades. Their endangered status has prompted extensive research and genetic monitoring, yielding large, multiyear microsatellite datasets upon which future conservation efforts can build. We achieved an overall microsatellite genotype concordance rate of 94.5% comparing high-throughput sequencing results to PCR based-fragment size results. All discrepancies occurred at complex loci containing multiple insertions and/or deletions (indels). Physically linked indels or single nucleotide polymorphisms (SNPs) occurring within the loci were misinterpreted as independent insertions, underscoring the need for genotyping tools that incorporate phasing when genotyping. To evaluate coverage effects, we downsampled high-throughput sequence data from 30x to 2x. Concordance remained high at 20 to 30x but dropped sharply at 10x, with 5x and 2x having discordant genotypes or insufficient coverage for genotyping. Accurate genotyping required both sufficient depth and number of reads spanning the entire repeat regions. Our results show that short-read whole-genome sequencing can recover microsatellite genotypes with high accuracy when paired with careful variant interpretation. By aligning historical PCR datasets with modern sequencing data, we can preserve decades of genetic insight and strengthen long-term monitoring of at-risk populations.

Animals↗

PubMed2026

Genome-wide variation analysis of two Salvia hispanica L. genotypes and implication for associations with metabolic and adaptive traits.

BACKGROUND: Advances in next-generation sequencing have accelerated genome-wide exploration of genetic diversity in underutilized oilseed crops. Salvia hispanica L. (chia), a high-nutrient pseudocereal rich in omega-3 fatty acids, is increasingly valued for its health benefits and commercial potential, yet it remains poorly characterized at the genomic level. Understanding the scale and nature of genomic variation is essential for improving complex traits such as oil yield, stress tolerance, and seed quality. METHODS: Two contrasting chia genotypes, Black-chia (CACH-B) and White- chia (CACH-W), were resequenced using the Bio-Resequencing Toolkit (BRT) pipeline. High-coverage sequencing, with a mapping rate exceeding 99% and an average depth of approximately 28×, facilitated the detection and annotation of single-nucleotide polymorphisms (SNPs), insertions and deletions (InDels), copy-number variations (CNVs), and structural variants (SVs). The functional classification of variant impacts enabled the identification of genes potentially linked to metabolic and adaptive traits. RESULTS: A total of 1.97 million SNPs, 401,493 InDels, 836 CNVs, and 15,288 SVs were identified across the chia genome. Notably, approximately 53% of exonic SNPs were non-synonymous (dN/dS ≈ 1.28), predominantly affecting lipid metabolism, transcriptional regulation, and stress response pathways, potentially altering key agronomic traits. In addition, CNV hotspots were concentrated in chromosomes 3 and 6, overlapping MYB, WRKY, and bZIP transcription factor loci, may potentially be involved in stress tolerance and yield. Furthermore, structural rearrangements, including inversions and duplications within the FAD2, FAD3, and CYP450 gene clusters, were potentially associated with seed pigmentation and omega-3 biosynthesis, pointing to their potential breeding relevance. Observed heterozygosity (Hₒ ≈ 0.71) and nucleotide diversity (π ≈ 7 × 10-3) indicated moderate to high allelic richness. In addition, the low FST value (0.038) indicates substantial genomic similarity between the two genotypes. CONCLUSION: This study presents the first comprehensive map integrating SNPs, CNVs, and SVs in S. hispanica L. The results reveal a structurally dynamic genome characterized by substantial sequence and structural variation, providing valuable insights into genomic diversity and potential adaptive mechanisms in chia. The coexistence of high SNP diversity and abundant structural variation underpins chia's nutritional specialization and environmental resilience. These results deliver a foundational genomic resource for marker-assisted breeding, genome-wide association studies, and the development of climate-resilient chia cultivars.

Copy-number variation, structural variation↗