Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Single-nucleotide variation”

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.

At least 55 records · Page 3Linked to original sources

KDM3B Regulates Postradiation Fibrotic Responses in Prostate Stroma via N6-methyladenosine Modification of LOX.

PURPOSE: Genome-wide association studies have uncovered single-nucleotide polymorphisms (SNPs) linked to radiation therapy (RT)-induced toxicities in patients with prostate cancer. SNP rs17599026, located in intron 21 of the KDM3B gene, has been associated with late-onset urinary toxicity, with an increased frequency of urination observed 2 years post-RT compared with pretreatment conditions. This study aimed to explore the underlying mechanisms driving this association. METHODS AND MATERIALS: A clustered regularly interspaced short palindromic repeats-dead Cas9 prime editing system was used to mimic KDM3B genetic variants in prostate stromal cell lines. Murine models with wild-type and heterozygous Kdm3b genotypes were used to assess fibrosis following radiation. RNA immunoprecipitation, transcript stability assays, and protein analysis elucidated the role of N6-methyladenosine (m6A) modification in regulating lysyl oxidase (LOX) expression. α-ketoglutarate (α-KG) supplementation was tested for its effects on KDM3B protein stability, LOX expression, and fibrosis mitigation. RESULTS: The rs17599026 SNP reduced KDM3B protein expression via circular RNA and microRNA-mediated mechanisms, leading to decreased m6A modification and increased stability of LOX messenger RNA. Elevated LOX expression promoted collagen cross-linking and fibrosis in prostate stroma. α-KG supplementation restored KDM3B protein levels, reduced LOX expression, and mitigated fibrosis in vitro and in vivo. CONCLUSIONS: KDM3B genetic variations influence radiation-induced fibrosis through posttranscriptional regulation of LOX. Dietary α-KG supplementation may serve as a mechanism-based strategy to alleviate radiation toxicity in patients with prostate cancer, offering a potential therapeutic pathway to improve treatment outcomes.

Male↗

Identification and masking of artifactual and misleading within-host variants in deep-sequencing SARS-CoV-2 data.

Deep-sequencing data are increasingly used to study within-host viral diversity and to inform evolutionary inference. For SARS-CoV-2, analyses based on intra-host single-nucleotide variants (iSNVs) have been widely applied to quantify within-host diversity and infer transmission dynamics. However, these applications critically depend on the reliable identification of low-frequency variants, which remain vulnerable to systematic and technical artifacts. In this study, we show that recurrent artifactual iSNVs are common in large-scale SARS-CoV-2 sequencing data and can persist even under conservative minor allele frequency thresholds. Using data from the UK's Office for National Statistics COVID-19 Infection Survey, we demonstrate that such artifacts are predominantly sequencing center-specific rather than primer-specific. Each center exhibits a modest, distinct set of recurrent artifactual variants showing little overlap with sites routinely masked at the consensus level. To address this, we developed a systematic, dataset-aware framework that uses recurrence within sequencing datasets to identify small, noise-adapted sets of artifactual iSNVs to mask. Applying this framework reduces spurious sharing of low-frequency variants between samples and qualitatively alters downstream inferences, including estimates of within-host diversity and transmission bottleneck sizes. Although this study focused on SARS-CoV-2, it is likely that recurrent artifactual iSNVs will be problematic for other viruses as mass-sequencing becomes increasingly routine. Together, these findings highlight the importance of explicit, dataset-aware artifact control for robust inference from within-host variation, particularly as genomic studies increasingly seek to exploit sub-consensus diversity in rapidly evolving pathogens.

Humans↗

De novo rare EMX2 variants lead to idiopathic hypogonadotropic hypogonadism.

PURPOSE: The genetic etiology of infertility remains unknown. To identify genes for human infertility, we applied a de novo variant analysis in 142 parent-proband trios with idiopathic hypogonadotropic hypogonadism (IHH), an infertility disorder caused by gonadotropin-releasing hormone (GnRH) deficiency. METHODS: Rare de novo copy-number and single-nucleotide variants (CNVs and SNVs) were called from exome sequencing data of the IHH trios. An association study of common EMX2 variants and disease outcomes was performed in the Massachusetts General Brigham Biobank (N = 65,253). GnRH neuronal development and migration was studied in organotypic explants with knocked down of Emx2 and in a mouse model lacking Emx2. RESULTS: We identified that the gene EMX2 harbored both rare de novo CNVs and SNVs. Rare de novo EMX2 variants led to IHH, developmental delay, and hearing loss. Common EMX2 variants were linked to infertility, Parkinson disease, and hearing loss. Knockdown of Emx2 in nasal explants resulted in attenuated GnRH cell migration and GnRH cells were confined to nasal regions of Emx2 knockout (KO) mice, consistent with IHH pathogenesis. CONCLUSION: By utilizing a de novo variant analysis and cellular assays, EMX2 was uncovered as a gene for human infertility.

Humans↗

Association of FOXC1 Duplications With Juvenile Open-Angle Glaucoma.

IMPORTANCE: While FOXC1 single-nucleotide variants and deletions are well-established causes of Axenfeld-Rieger syndrome, few FOXC1 duplications have been reported. This study investigated families with duplications encompassing the FOXC1 gene to refine the associated phenotypic spectrum and contribution to glaucoma. OBJECTIVE: To investigate the prevalence and phenotype of FOXC1 duplications in 2 large glaucoma registries. DESIGN, SETTING, AND PARTICIPANTS: This retrospective observational genetic cohort study included participants recruited from the Australian & New Zealand Registry of Advanced Glaucoma (ANZRAG) and the Massachusetts Eye and Ear (MEE) cohort from 2008 through 2025. Participants with glaucoma, and available relatives, underwent genomic testing to identify duplications encompassing FOXC1 using exome sequencing and genotyping arrays (ANZRAG) or whole-genome sequencing (MEE). Data analyses were conducted from 2022 through 2025. MAIN OUTCOMES AND MEASURES: Prevalence of FOXC1 duplications, age at glaucoma onset, and phenotype, including ocular and systemic features. RESULTS: Twenty individuals from 10 families (50% female and 50% male; 70% self-described as broadly European [Australian/British, British, English/German, English/Polish, European, or Scottish], 25% as Asian [Chinese or Filipino], and 5% as Latin American [Salvadoran]) were identified with FOXC1 duplications. All genetically tested individuals were diagnosed with glaucoma, demonstrating high penetrance. Seventeen individuals were referred with juvenile open-angle glaucoma (JOAG), 1 with primary open-angle glaucoma, 1 with primary congenital glaucoma, and 1 with anterior segment dysgenesis. The diagnosis of 4 individuals from 1 family with ectropion uveae was revised to anterior segment dysgenesis. Systemic features were reported for 2 participants (10.5%), including subtle dental findings and mild facial dysmorphism. Duplications encompassing FOXC1 were among the most common monogenic contributors to JOAG. In the ANZRAG group, they accounted for 13.5% (95% CI, 6.7%-25.3%) of JOAG probands with a genetic diagnosis, second to MYOC (53.8%; 95% CI, 40.5%-66.7%). In the MEE group, FOXC1 duplications accounted for 9.5% (95% CI, 2.7%-28.9%) of JOAG probands with a genetic diagnosis. CONCLUSIONS AND RELEVANCE: These findings suggest FOXC1 duplications are an underrecognized, highly penetrant, but variably expressive, genetic variation associated with JOAG. Findings for the relatively modest number of individuals in the retrospective study were associated with wide confidence intervals. This limitation is often inherent to studies of JOAG, a rare condition for which individual genetic variants account for only a subset of cases. Despite this, the findings highlight the genetic heterogeneity of JOAG and support the potential importance of considering routine genetic copy-number variant analysis for individuals with JOAG.

Humans↗

Expression of ZNF804A in human brain and alterations in schizophrenia, bipolar disorder, and major depressive disorder: a novel transcript fetally regulated by the psychosis risk variant rs1344706.

IMPORTANCE: The single-nucleotide polymorphism rs1344706 in the zinc finger protein 804A gene (ZNF804A) shows genome-wide association with schizophrenia and bipolar disorder. Little is known regarding the expression of ZNF804A and the functionality of rs1344706. OBJECTIVES: To characterize ZNF804A expression in human brain and to investigate how it changes across the life span and how it is affected by rs1344706, schizophrenia, bipolar disorder, and major depressive disorder. DESIGN, SETTING, AND PARTICIPANTS: Molecular and immunochemical methods were used to study ZNF804A messenger RNA (mRNA) and ZNF804A protein, respectively. ZNF804A transcripts were investigated using next-generation sequencing and polymerase chain reaction-based methods, and ZNF804A protein was investigated using Western blots and immunohistochemistry. Samples of dorsolateral prefrontal cortex and inferior parietal lobe tissue were interrogated from 697 participants between 14 weeks' gestational age and age 85 years, including patients with schizophrenia, bipolar disorder, or major depressive disorder. MAIN OUTCOMES AND MEASURES: Quantitative measurements of ZNF804A mRNA and immunoreactivity, and the effect of diagnosis and rs1344706 genotype. RESULTS: ZNF804A was expressed across the life span, with highest expression prenatally. An abundant and developmentally regulated truncated ZNF804A transcript was identified, missing exons 1 and 2 (ZNF804AE3E4) and predicted to encode a protein lacking the zinc finger domain. rs1344706 influenced expression of ZNF804AE3E4 mRNA in fetal brain (P&#x2009;=&#x2009;.02). In contrast, full-length ZNF804A showed no association with genotype (P&#x2009;>&#x2009;.05). ZNF804AE3E4 mRNA expression was decreased in patients with schizophrenia (P&#x2009;=&#x2009;.006) and increased in those with major depressive disorder (P&#x2009;<&#x2009;.001), and there was a genotype-by-diagnosis interaction in bipolar disorder (P&#x2009;=&#x2009;.002). ZNF804A immunoreactivity was detected in fetal and adult human cerebral cortex. It was localized primarily to pyramidal neurons, with cytoplasmic as well as dendritic and nuclear staining. No differences in ZNF804A-immunoreactive neurons were seen in schizophrenia or related to rs1344706 (P&#x2009;>&#x2009;.05). CONCLUSIONS AND RELEVANCE: rs1344706 influences the expression of ZNF804AE3E4, a novel splice variant. The effect is limited to fetal brain and to this isoform. It may be part of the mechanism by which allelic variation in ZNF804A affects risk of psychosis. ZNF804A is translated in human brain, where its functions may extend beyond its predicted role as a transcription factor.

Adolescent↗

Variants in autophagy-related genes and clinical characteristics in melanoma: a population-based study.

Autophagy has been linked with melanoma risk and survival, but no polymorphisms in autophagy-related (ATG) genes have been investigated in relation to melanoma progression. We examined five single-nucleotide polymorphisms (SNPs) in three ATG genes (ATG5; ATG10; and ATG16L) with known or suspected impact on autophagic flux in an international population-based case-control study of melanoma. DNA from 911 melanoma patients was genotyped. An association was identified between (GG) (rs2241880) and earlier stage at diagnosis (OR 0.47; 95% Confidence Intervals (CI)&#xa0;=&#xa0;0.27-0.81, P&#xa0;=&#xa0;0.02) and a decrease in Breslow thickness (P&#xa0;=&#xa0;0.03). The ATG16L heterozygous genotype (AG) (rs2241880) was associated with younger age at diagnosis (P&#xa0;=&#xa0;0.02). Two SNPs in ATG5 were found to be associated with increased stage (rs2245214 CG, OR 1.47; 95% CI&#xa0;=&#xa0;1.11-1.94, P&#xa0;=&#xa0;0.03; rs510432 CC, OR 1.84; 95% CI&#xa0;=&#xa0;1.12-3.02, P&#xa0;=&#xa0;0.05). Finally, we identified inverse associations between ATG5 (GG rs2245214) and melanomas on the scalp or neck (OR 0.20, 95% CI = 0.05-0.86, P&#xa0;=&#xa0;0.03); ATG10 (CC) (rs1864182) and brisk tumor infiltrating lymphocytes (TILs) (OR 0.42; 95% CI&#xa0;=&#xa0;0.21-0.88, P&#xa0;=&#xa0;0.02), and ATG5 (CC) (rs510432) with nonbrisk TILs (OR 0.55; 95% CI&#xa0;=&#xa0;0.34-0.87, P&#xa0;=&#xa0;0.01). Our data suggest that ATG SNPs might be differentially associated with specific host and tumor characteristics including age at diagnosis, TILs, and stage. These associations may be critical to understanding the role of autophagy in cancer, and further investigation will help characterize the contribution of these variants to melanoma progression.

Adult↗

Copy number variants in BRCA1 and BRCA2 genes in Polish patients with breast and ovarian cancer.

PURPOSE: BRCA1 and BRCA2 are key susceptibility genes in hereditary breast and ovarian cancer (HBOC), with mutational status guiding PARP inhibitor therapy. While single-nucleotide variants (SNVs) predominate, the prevalence of copy number variants (CNVs) varies significantly across different populations. This study aims to determine the incidence of BRCA1/2 CNVs in the Polish population, where data remain scarce due to non-mandatory CNV testing. METHODS: We retrospectively analysed the results of genetic tests assessing the presence of BRCA1/2 CNVs performed in 2720 individuals tested at the Lower Silesian Oncology Centre (2021-2024), including 2702 breast/ovarian cancer patients and 18 unaffected relatives. The mean age was 54.7&#x2009;&#xb1;&#x2009;15.15&#xa0;years. Genetic testing involved DNA extraction, NGS, and MLPA for CNV confirmation. Variants were classified according to ACMG-AMP guidelines and verified through independent testing. RESULTS: In this study, no BRCA2 CNVs were identified, consistent with previous Central European findings. Pathogenic BRCA1 CNVs were detected in 0.85% of the analyzed cohort and in 0.52% of the cancer patient subgroup, affecting 23 individuals from 13 families. Eight distinct BRCA1 CNVs were detected, the most common being exon 21 deletion. Affected families exhibited a high incidence of HBOC-related cancers, with early-onset breast cancer and a notable proportion of triple-negative breast cancer cases. CONCLUSIONS: This study highlights the clinical significance of BRCA1 CNVs in Polish patients with HBOC-spectrum cancers and their families. Although rare, these variants were associated with aggressive cancer phenotypes and early onset. Given their diagnostic and therapeutic implications, BRCA1 CNVs should be routinely analysed in high-risk families to ensure accurate detection and personalised treatment planning.

Humans↗

Panmixia in a Widespread Butterfly: High Dispersal and Ecological Generalism Buffer Against Landscape Fragmentation.

Habitat fragmentation is widely expected to reduce population connectivity and increase genetic differentiation, although the strength of these effects depends on species-specific traits such as dispersal ability. Here, we investigated the population genetic structure of the cosmopolitan butterfly, Pieris rapae L. (Lepidoptera: Pieridae), across western Germany using genome-wide single-nucleotide polymorphism (SNP) data. To analyze the effects of landscape structure on genetic connectivity, we applied a paired study design comprising four landscape pairs, each consisting of a highly intensified, modern agricultural landscape and a more heterogeneous, traditional landscape. Our results revealed no evidence of genetic differentiation. Pairwise FST values were close to zero; we detected no isolation by distance, and clustering analyses supported a single genetic population. No meaningful associations between genetic variation and environmental variables were detected, with landscape effects explaining less than 0.4% of genomic variation. Consequently, we found no evidence for stronger genetic structuring in modern compared to more connected traditional landscapes. Our results suggest that extensive habitat fragmentation does not necessarily translate into reduced genetic connectivity in highly mobile, generalist species. In P. rapae , high dispersal ability and ecological generalism appear to buffer against the genetic consequences of landscape modification, resulting in panmictic population structure even across strongly contrasting agricultural landscapes.

Pieris rapae↗

TET2 promotes monocyte inflammatory activation in asthma via ALKBH5-m6A regulation and PI3K signaling: evidence from m6A-SNP and single-cell analyses.

Asthma is a complex inflammatory airway disease with strong genetic determinants, yet the functional relevance of most asthma-associated non-coding variants remains unclear. Emerging evidence suggests that N6-methyladenosine (m6A) modification may serve as a critical epitranscriptomic link between genetic variation and immune regulation. In this study, we aimed to systematically identify functionally relevant m6A-regulated genes in asthma by integrating large-scale GWAS data, m6A-SNP annotations, and single-cell transcriptomic analyses, and to investigate their roles in monocyte-driven airway inflammation. We identified TET2 as a key m6A-regulated gene associated with both asthma and lung function, which was selectively upregulated in monocytes during asthma and accompanied by activation of inflammatory and PI3K signaling pathways. Mechanistic experiments further demonstrated that inflammatory stimulation induced ALKBH5 expression, reduced m6A modification of TET2 mRNA, and increased TET2 protein levels, thereby promoting PI3K/AKT signaling and pro-inflammatory cytokine production, whereas inhibition of TET2 or ALKBH5 attenuated these effects. Collectively, these findings demonstrate that ALKBH5-mediated m6A regulation of TET2 enhances PI3K/AKT signaling in monocytes, thereby promoting inflammatory responses in asthma. Our study establishes TET2 as a key m6A-regulated gene linking genetic susceptibility to monocyte-driven inflammation, and highlights the ALKBH5-m6A-TET2 axis as a potential therapeutic target for modulating aberrant immune responses in asthma.

Humans↗

Genome-wide detection of human 5' UTR variants that impact protein translation.

The 5' untranslated region (5' UTR) of messenger RNAs (mRNAs) plays a central role in regulating protein synthesis initiation, particularly through the Kozak sequence and upstream open reading frames (uORFs). Genetic variants within these regulatory elements could affect translation, altering gene expression and contributing to clinical phenotypes in humans. We developed a computational method called 5ULTRA (5' Untranslated Region Annotation) for analysis of whole-exome sequencing and whole-genome sequencing data to detect, annotate, and prioritize 5' UTR variants with potential translation impact. 5ULTRA identifies single-nucleotide variants, indels, and splicing variants that affect uORFs by creating or disrupting start/stop codons and that alter Kozak sequence strength of either the uORFs or the main coding sequence. 5ULTRA incorporates recent uORF databases and provides comprehensive annotations. 5ULTRA implements a machine-learning score to prioritize candidate variants with predicted effects on translation and also provides specific mechanistic predictions. The score correlates strongly with experimentally measured protein-level effects of 5' UTR variants. We applied 5ULTRA to multiple genetics datasets across diverse disease contexts, identifying candidate variants including potential cancer-driving somatic mutations predicted to decrease ABI1 level or increase NRAS abundance; common variants associated with traits such as multiple sclerosis, lung function, and cardiovascular function, by altering protein levels of TAGAP, VRTN, and SPAAR, respectively; and rare germline variants in our cohort, including a splicing variant of RPSA leading to 5' UTR sequence alteration that causes congenital asplenia and a variant of TNF that could predispose to tuberculosis.

Humans↗

Combined somatic mutation and transcriptome analysis reveals region-specific differences in clonal architecture in human cortex.

The human cerebral cortex is specialized into regions, but little is known about how human cellular lineages shape cortical regional variation and neuronal cell-type distribution during development. Here, we map single-cell lineages of human cortical regions and neuronal subtypes using >1,000 somatic single-nucleotide variants (sSNVs) identified from deep bulk whole-genome sequencing and analyzed over 25 regions and >72,000 single cells. In the fronto-parietal cortex, sSNVs are rarely restricted, marking neuron-generating clones that disperse into neighboring regions. In contrast, the primary visual cortex harbors 30%-70% more sSNVs than the neighboring secondary visual cortex. Clones at this border exhibit more restricted dispersion, suggesting late developmental lineage segregation. Single-nucleus sSNV and whole-transcriptome analysis reveal glutamatergic neuron clones with modest regional restrictions that share low-mosaic sSNVs with some GABAergic neurons, suggesting a recent dorsal cortical progenitor. Our analysis reveals human-specific cortical lineage patterns, regional differences in clonal patterns, and late divergence of some glutamatergic/GABAergic lineages.

Humans↗

Cell-type-specific enrichment of somatic aneuploidy in the mammalian brain.

Somatic mutations alter the genomes of a subset of an individual's brain cells, impacting gene regulation and contributing to disease processes. Mosaic single-nucleotide variants have been characterized with single-cell resolution in the brain, but we have limited information about large-scale structural variation such as whole-chromosome duplication or loss. We used a dataset of over 415,000 single-cell DNA methylation and chromatin conformation profiles from the adult mouse brain to comprehensively identify and characterize aneuploid cells. Somatic trisomy events were strongly enriched on chromosome 16, which is syntenic with human chromosome 21. We also observed a specific enrichment of chromosome gain and loss events in specific cell types, including Pons neurons and oligodendrocyte precursor cells. Chromosome 16 trisomy occurred in multiple cell types and across brain regions, suggesting that nondisjunction is a recurrent feature of somatic structural variation in the brain.

Animals↗

Optimizing GRIDSS for clinical use: A targeted NGS filtering strategy for germline structural variant detection.

Detecting intermediate-sized structural variants (SVs) remains challenging in diagnostics, as tools for single-nucleotide and copy-number variants, particularly read-depth-based methods, are often insufficient. GRIDSS addresses this gap by integrating paired-end mapping, split-read analysis, and assembly-based approaches. However, its use in targeted sequencing and diagnostic workflows remains complex. NGS panel data from 9726 patients with suspected hereditary cancer were analyzed using GRIDSS. A filtering strategy was developed to prioritize clinically relevant germline SVs. Multiple parameter settings were tested to optimize performance. The initial dataset of 1,307,592 variants was reduced to 89 candidates after applying the selected filtering strategy. Of these, 24 had been previously detected by routine callers and were not further analyzed. Among the remaining 65, 13 were considered likely true positives after visual inspection using IGV. Experimental validation was performed by Sanger/Nanopore long-read sequencing for these variants, all of which were confirmed. Eight were classified as (likely) pathogenic, including two frameshift duplications in MSH6, one splicing variant in BARD1, and five mobile element insertions in APC, BRCA2, and PALB2. Altogether, GRIDSS implementation increased diagnostic yield while maintaining feasibility for diagnostic workflows. Comprehensive workflow scheme for germline structural variant detection and results in our diagnostic setting.

Humans↗

Determinants of functional burden pleiotropy and gene dosage responses across human traits.

Pleiotropic and monotonic effects of gene dosage are central to understanding comorbidities in developmental pediatric and psychiatric disorders, yet the underlying biological processes are not well characterized. Here we develop a functional burden analysis to investigate the association of all protein-coding copy-number variants, genome-wide, with 43 complex traits in approximately 500,000 UK Biobank participants. We test variant associations disrupting 172 tissue or cell-type gene sets, finding associations for all traits, which we replicate in the All of Us cohort. Functional burden pleiotropy, defined as the number of traits significantly associated with a gene set, correlates with genetic constraint and is higher for brain than non-brain functions, even after normalizing for genetic constraint. Levels of pleiotropy, measured by burden correlation, are similar in deletions and loss-of-function single-nucleotide variants, and higher than in common variants and duplications. Most gene dosage responses are non-monotonic, with deletions and duplications showing same-direction effects, and monotonic responses decrease with genetic constraint. We observe associations between functional gene sets and traits for either deletions or duplications, but rarely both, with negatively correlated effect sizes. Together, these results link genetic constraint and brain-specific mechanisms to the whole-body multimorbidity of neurodevelopmental and psychiatric conditions.

Humans↗

scSNViz: visualization and analysis of cell-specific expressed SNVs.

MOTIVATION: Accurately characterizing expressed genetic variation at the single-cell level is essential for understanding transcriptional heterogeneity, allelic regulation, and mutational dynamics within complex tissues. However, few tools enable comprehensive visualization and quantitative analysis of expressed variants across individual cells. RESULTS: scSNViz is an R package for the exploration, quantification, and visualization of expressed single-nucleotide variants (SNVs) from cell-barcoded single-cell RNA sequencing (scRNA-seq) data. The software supports estimation of variant allele fractions, clustering of SNV expression profiles, and 2D and 3D visualization of individual SNVs or user-defined SNV groups. Beyond visualization, scSNViz facilitates investigation of cell-, cluster-, or lineage-specific variant expression patterns, as well as allelic dynamics including imprinting, random allele inactivation, and transcriptional bursting. It interoperates seamlessly with established single-cell frameworks-Seurat for clustering, Slingshot for trajectory inference, scType for cell-type annotation, and CopyKat for copy-number profiling-enabling integrative multi-omic analyses of expressed variation. AVAILABILITY AND IMPLEMENTATION: scSNViz is implemented in R and freely available at https://github.com/HorvathLab/scSNViz (DOI: 10.5281/zenodo.17307516). The package includes comprehensive documentation and example workflows designed for users with limited bioinformatics experience.

Software↗

Genomic prediction and genome-wide association study for liver abscesses in crossbred beef cattle.

Liver abscesses are a concern in feedlot cattle, and little is known about the role of genetics in their development. This study aimed to estimate genetic parameters and to identify single-nucleotide polymorphisms (SNPs) associated with liver abscesses. Crossbred cattle representing 18 breeds in the U.S. Meat Animal Research Center Germplasm Evaluation Program were phenotyped for liver abscesses at slaughter (n&#x2005;=&#x2005;9,044). Seventeen percent of cattle had liver abscesses. These cattle had genotypes that were imputed to sequence variant genotypes. After filtering and quality control, 340,723 SNPs were used in the analysis. Liver abscess prevalence was modeled with a single-step genomic best linear unbiased prediction (ssGBLUP) threshold model using a Bayesian framework. The model included contemporary group (sex, treatment group, and slaughter date), additive genomic, and residual effects. Genomic heritability was 0.039 (95% highest posterior density&#x2005;=&#x2005;0.005, 0.081), which was very small. To assess prediction quality, a 5-fold random cross-validation structure was used. Method Linear Regression was used to assess accuracy, bias, and dispersion by comparing estimated breeding values (EBV) from full and reduced analyses. Cross-validation metrics showed EBV based on genotypes had 0.05 reliability (SD&#x2005;<&#x2005;0.01) with no bias relative to EBV based on genotypes and phenotypes. For the genome-wide association study, SNP effects were back calculated from the EBV solutions from ssGBLUP. No SNPs were associated with liver abscesses at a Benjamini-Hochberg adjusted 0.05 significance level. Although a large dataset was used, this result was because of the low genomic heritability and imprecise EBV used to calculate SNP effects. Based on these results, environmental factors contribute to most of the variation in liver abscesses. Genetic selection to reduce liver abscesses would be slow because of the low genomic heritability, measurement late in life, and inability to measure breeding animals. A faster approach would be finding additional environmental interventions that maintain animal performance.

Animals↗

Conservation genomics of a threatened subtropical Rhododendron species highlights the distinct conservation actions required in marginal and admixed populations.

With the impact of climate change and anthropogenic activities, the underlying threats facing populations with different evolutionary histories and distributions, and the associated conservation strategies necessary to ensure their survival, may vary within a species. This is particularly true for marginal populations and/or those showing admixture. Here, we re-sequence genomes of 102 individuals from 21 locations for Rhododendron vialii, a threatened species distributed in the subtropical forests of southwestern China that has suffered from habitat fragmentation due to deforestation. Population structure results revealed that R. vialii can be divided into five genetic lineages using neutral single-nucleotide polymorphisms (SNPs), whereas selected SNPs divide the species into six lineages. This is due to the Guigu (GG) population, which is identified as admixed using neutral SNPs, but is assigned to a distinct genetic cluster using non-neutral loci. R. vialii has experienced multiple genetic bottlenecks, and different demographic histories have been suggested among populations. Ecological niche modeling combined with genomic offset analysis suggests that the marginal population (Northeast, NE) harboring the highest genetic diversity is likely to have the highest risk of maladaptation in the future. The marginal population therefore needs urgent ex situ conservation in areas where the influence of future climate change is predicted to be well buffered. Alternatively, the GG population may have the potential for local adaptation, and will need in situ conservation. The Puer population, which carries the heaviest genetic load, needs genetic rescue. Our findings highlight how population genomics, genomic offset analysis, and ecological niche modeling can be integrated to inform targeted conservation.

Rhododendron↗

Comparative analysis of distinct genomic landscapes in young-onset gBRCA1/2 breast cancer.

Carriers of germline BRCA1/2 pathogenic variants (gBRCA1/2 PVs) have elevated young-onset breast cancer risk. To define the pretreatment genomic landscapes of young-onset gBRCA-associated breast cancer, we evaluated 136 treatment-naive tumors diagnosed before age 50 in the prospective POSH study and 66 noncarriers from The Cancer Genome Atlas. Using whole-exome sequencing, we analyzed somatic variation, allele-specific loss of heterozygosity (asLOH), homologous recombination deficiency (HRD), and single-base substitution (SBS) signatures. gBRCA1 and gBRCA2 breast cancers had high rates of asLOH but differed significantly in average HRD scores and median SBS composition of signatures SBS1 (aging-associated), SBS18 (ROS-associated), and SBS3 (HRD-associated). Compared with gBRCA2 tumors, gBRCA1 tumors with asLOH were significantly enriched for alterations in hallmark ROS, DNA repair, and epithelial-mesenchymal transition pathways. In ER-positive, HER2-negative tumors from gBRCA1/2 carriers compared with noncarriers, we found significant enrichment of RB1, TP53, FAT1, and MYC single-nucleotide variants, indels, and copy number variants associated with CDK4/6 inhibitor (CDK4/6i) resistance. Together, these findings demonstrate significant differences between gBRCA1- and gBRCA2-associated breast cancers, and preexisting CDK4/6i resistance mechanisms, supporting prospective trials comparing individualized therapy for gBRCA1 versus gBRCA2 carriers and comparing poly(ADP-ribose) polymerase inhibitors versus CDK4/6i for ER-positive gBRCA1/2-associated breast cancer.

Humans↗