PubMed HealthSearch

SEARCH · PubMed Health

Results for “Genomic Structural Variation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original 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

Integrative Genotyping and Analysis of Canine Structural Variation Using Long-read and Short-read Data.

Structural variation makes an important contribution to canine evolution and phenotypic differences. Although recent advances in long-read sequencing have enabled the generation of multiple canine genome assemblies, most prior analyses of structural variation have relied on short-read sequencing. To offer a more complete assessment of structural variation in canines, we performed an integrative analysis of structural variants present in 12 canine samples with available long-read and short-read sequencing data along with genome assemblies. Use of long-reads permits the discovery of heterozygous variation that is absent in existing haploid assembly representations while offering a marked increase in the ability to identify insertion variants relative to short-read approaches. Examination of the size spectrum of structural variants shows that dimorphic LINE-1 and SINE variants account for over 45% of all deletions and identified 1,410 LINE-1s with intact open reading frames that show presence-absence dimorphism. Using a graph-based approach, we genotype newly discovered structural variants in an existing collection of 1,879 resequenced dogs and wolves, generating a variant catalog containing a 56.5% increase in the number of deletions and 705% increase in the number of insertions previously found in the analyzed samples. Examination of allele frequencies across admixture components present across breed clades identified 283 structural variants evolving with a signature of selection.

Animals

A chromosome-level, haplotype-resolved genome assembly for the barn owl, Tyto alba.

Recent advances in long-read sequencing have enabled near telomere-to-telomere (T2T) assemblies across diverse taxa. However, avian genomes remain challenging due to numerous microchromosomes, small, typically < 20Mb, DNA molecules that are gene-, GC-, and repeat-rich. As a consequence, microchromosomes are often missing from genome assemblies. Here, we present a chromosome-level, haplotype-resolved genome assembly for the Western barn owl (Tyto alba). Using a trio-binning strategy with Illumina parental reads combined with PacBio HiFi and Oxford Nanopore Technologies data, we generated two phased contig sets. These were scaffolded into 40 linkage groups using a linkage map. Comparative analyses identified unplaced HiFi scaffolds corresponding to microchromosomes, which we integrated into six additional microchromosomes using long reads information. The two assemblies present 46 chromosomes, matching the karyotype of the species. They exhibit strong synteny between parental haplotypes, except for a &#x223c;38 Mb complex region on chromosome 7 containing nested inversions. This high-quality reference provides a haplotype-resolved and chromosome-level genome for Strigiformes, enabling fine-scale studies of structural variation and avian genome evolution.

Tyto alba

Reproductive Isolation due to Divergent Ecological Selection Is Accompanied by Vast Genomic Instability in Experimentally Evolved Yeast Populations.

Populations evolving independently in divergent environments accumulate genetic differences and potentially evolve reproductive isolation as a by-product of divergence. The speed and mechanisms underlying this process are difficult to investigate because we rarely get the opportunity to witness them in natural settings, and histories of selection and gene flow between populations are often unknown. Here, we experimentally evolved yeast for 1000 generations of evolution in both divergent and parallel environments. At regular time points during experimental evolution, we made crosses between parallel- and divergent-evolving populations to measure postzygotic reproductive isolation (gamete viability). We used whole genome population sequencing to determine the mutational load, the number and types of structural variation, and other genomic features of the parent, F1 and F2 intraspecific hybrids. We found evidence for large-scale phenotypic and genome-wide differentiation in response to divergent laboratory selection. Divergent-selected populations produced hybrids with reduced gamete viability-a classic signature of postzygotic reproductive isolation in the form of hybrid breakdown. Parallel-selected populations, on the other hand, remained more reproductively compatible (with exceptions). We found that F2 hybrid genomes contained vast genomic instability, that is, new structural variants (especially insertions, deletions and interchromosomal translocations) that were not observed in parent and F1 genomes, which is likely a result of chromosome missegregation and recombination errors in hybrid meiosis. Our results provide phenotypic and genomic evidence that partial reproductive isolation evolved due to adaptation to divergent environments, consistent with predictions of ecological speciation theory.

Reproductive Isolation

The human growth hormone gene locus: structure, evolution, and allelic variations.

Genomic clones containing the closely related genes for human growth hormone (hGH) and chorionic somatomammotropin (hCS) were obtained from genomic bacteriophage lambda and cosmid libraries. The entire GH/CS chromosomal locus was reconstructed utilizing overlapping restriction fragments characterized from the isolated clones. The hGH/hCS locus contains two GH genes and three CS genes spanning 48 kb of DNA in the order: 5'-(hGH-1/hCS-5/hCS-1/hGH-2/hCS-2)-3', confirming analysis of cosmid clones obtained from a different human library (Barsh et al., 1983). To complete the characterization of the hCS genes, the nucleotide sequence of the hCS-5 gene was determined. Sequence analysis revealed a mutation of the 5' splice site at the exon II-intron B boundary, suggesting that the hCS-5 gene is a pseudogene. The nucleotide sequence of an allelic variant of the hCS-2 gene was determined and found to contain a single amino acid substitution and the deletion of a single codon. The hGH/hCS gene locus was further characterized by the localization of at least 27 Alu-type repetitive sequences and identification of three unique sequences in the vicinity of several hGH and hCS genes which define the probable breakpoints of the evolutionary duplication units. These data, combined with the nucleotide sequences of all five GH and CS genes, indicate that the hGH/hCS gene locus has evolved by duplication mechanisms. Evidence for the occurrence of at least one gene conversion event involving the hCS-1 gene precursor and the hCS-2 gene was found, indicating that the hGH/hCS gene locus has evolved by concerted mechanisms. The structure of the hCS genes is discussed in light of recent studies of CS genes from other mammalian species.

Alleles

Comparative Analysis of Mammalian Adaptive Immune Loci Revealed Spectacular Divergence and Common Genetic Patterns.

Adaptive immune responses are mediated by the production of adaptive immune receptors, antibodies, and T-cell receptors, which bind antigens, thus causing their neutralization. Unlike other proteins, adaptive immune receptors are not fully encoded in the germline genome and result from a complex of somatic processes collectively called V(D)J recombination affecting germline immunoglobulin (IG) and T-cell receptor (TR) loci consisting of template genes. While various existing studies report extreme diversity of antibodies and T-cell receptors, little is known about the diversity of germline IG and TR loci. To overcome this gap, the first comparative analysis of full-length sequences of IG/TR loci across 46 mammalian species from 13 taxonomic orders was performed. First, germline gene counts were shown to correlate in immunoglobulin heavy chain immunoglobulin heavy chain (IGH)/immunoglobulin lambda (IGL) loci and T-cell receptor alpha (TRA)/T-cell receptor beta (TRB) and anticorrelate in immunoglobulin kappa (IGK)/IGL, possibly indicating coevolution between corresponding chains. Second, structures of IG/TR loci were analyzed, and it was shown that IG/TR loci formed by long arrays of high multiplicity repeats are more common for species that have experienced population bottlenecks. Finally, haplotypes of IG/TR loci with little or no sequence similarity within a species were found, suggesting that they may have a limited potential for homologous recombination. These results demonstrate that IG/TR loci are rapidly evolving genomic regions whose structural variation is shaped by the population history of the species and open new perspectives for immunogenomics studies.

Animals

PARTAGE: Parallel analysis of replication timing and gene expression.

The human genome is partitioned into functional compartments that replicate at specific times during the S-phase. This temporal program, referred to as replication timing (RT), is co-regulated with the 3D genome organization, is cell type-specific, and changes during development in coordination with gene expression. Moreover, RT alterations are linked to abnormal gene expression, genome instability, and structural variation in multiple diseases, including cancer. However, mechanistic links between RT, large-scale 3D genome architecture, and transcriptional regulation remain poorly understood. A major limitation is that current approaches require the separate profiling of RT and transcriptomes from independent batches of samples, obscuring the complex co-regulation between the epigenome and transcriptome. Here, we developed PARTAGE, a multiomics approach that enables joint profiling of copy number variation (CNV), RT, and gene expression from the same sample, providing a more accurate integrative view of the complex relationships between RT and gene regulation.

Journal Article

Global variation in G+C content along vertebrate genome DNA. Possible correlation with chromosome band structures.

The global, rather than local, variation in G+C content along the nuclear DNA sequences of various organisms was studied using GenBank sequence data. When long DNA sequences of the genomes of Escherichia coli and Saccharomyces cerevisiae were examined, the levels of their G+C content (G+C%) were found to be within a narrow range around that of the whole genome. The G+C% levels for sequences of vertebrate genomes, however, were found to cover a wide range, showing that their genome is a mosaic of sequences with different G+C% levels, in each of which the sequence is fairly homogeneous in its G+C% for a very long distance. Through surveying a human genetic map and GenBank DNA sequences, the global variations in G+C% along the human genome DNA were found to be correlated with chromosome band structures.

Animals

Graph-based pan-genome reveals structural and functional diversity across oil palm domestication gradients.

BACKGROUND: Oil palm (Elaeis guineensis Jacq.), the world's most land-efficient oil crop, underpins global vegetable oil supply yet faces mounting constraints from limited expansion, climate stress, and disease pressure. These challenges highlight the urgent need for genomic resources that capture species-wide diversity to support sustainable improvement. While recent reference assemblies have advanced trait discovery, single linear genomes fail to represent the full spectrum of structural and gene-content variation, limiting resolution of agronomic alleles. RESULTS: Here, we constructed a graph-based pan-genome from 30 diverse oil palm assemblies representing wild, semi-domesticated, and commercial accessions. We characterized structural variants, gene presence-absence variation, and copy-number gains, with focusing on functional stratification and resistance gene dynamics. The graph-based pan-genome revealed extensive structural and gene-content variation, including a large conserved core, complemented by shell and unique fractions enriched or biased toward regulatory, stress-responsive, and defense-related functions. Structural variation and duplication-derived copy-number gains contributed substantially to gene-content diversity, with semi-domesticated accessions exhibiting the greatest variability. Resistance gene repertoires showed contrasting patterns: receptor-like kinases remained comparatively stable, whereas the CNL subclass of NLR genes contributed disproportionately to shell-genome variation and duplication-associated turnover. CONCLUSIONS: This graph-based pan-genome provides a curated multi-assembly reference and comparative framework for oil palm genomics. By capturing structural variants, gene-content variations, copy-number gains, and resistance gene dynamics across domestication gradients, it establishes a foundation for future pan-GWAS analysis, functional genomics, and molecular breeding strategies aimed at improving resilience and productivity in this globally important crop.

Arecaceae

Genome-wide SNP data reveal geographic structure and landscape-associated genomic differentiation in a widespread lizard in arid Eastern Central Asia.

Arid landscapes provide important systems for examining how geographic structure and environmental heterogeneity shape genomic differentiation. In topographically complex desert regions, however, it remains challenging to determine whether population structure primarily reflects landscape resistance, geographic distance, or contemporary environmental variation. Here, we use genome-wide SNP data to investigate population structure, phylogenetic relationships, historical gene flow, demographic history, and landscape correlates of genomic differentiation in the variegated racerunner (Eremias vermiculata), a widespread lacertid lizard across arid Eastern Central Asia. Analyses of 164 individuals recovered six geographically structured nuclear clusters associated with major desert basins and mountain-bounded regions. Nuclear phylogenies resolved two broad regional clades corresponding to northeastern and southwestern parts of the species' range, while PCA and ADMIXTURE analyses recovered six finer-scale genetic clusters. Mitochondrial phylogenies, based on combined NCBI-derived Cyt b and COI sequences from the same individuals, recovered four deeper maternal lineages. These patterns indicate overall phylogeographic agreement between nuclear and mitochondrial datasets, with genome-wide SNPs providing finer-scale resolution of population structure. Demographic reconstructions further uncovered regionally heterogeneous Late Pleistocene histories among clusters, including signals of expansion, stability, and decline. Landscape genomic analyses revealed that genomic differentiation is primarily associated with landscape resistance, particularly elevation and land cover, as well as geographic distance, whereas contemporary environmental variables explained comparatively little variation after controlling for spatial structure. Together, our results suggest that genomic differentiation in E. vermiculata reflects the interplay of persistent landscape configuration, historical connectivity, and region-specific demographic histories across arid Eastern Central Asia. More broadly, this study highlights the value of integrating phylogeographic and landscape genomic approaches for understanding population differentiation and evolutionary history in topographically heterogeneous desert ecosystems.

Arid Eastern Central Asia

Genetic Differentiation is Constrained to Chromosomal Inversions and Putative Centromeres in Locally Adapted Populations With Higher Gene Flow.

The impact of genome structure on adaptation is a growing focus in evolutionary biology, revealing an important role for structural variation and recombination landscapes in shaping genetic diversity across genomes and among populations. This is particularly relevant when local adaptation occurs despite gene flow, where clustering of differentiated loci can maintain locally adapted variants by reducing recombination between them. However, the limited genomic resources for nonmodel species, including reference genomes and recombination maps, have constrained our understanding of these patterns. In this study, we leverage the Atlantic silverside-a nonmodel fish with extensive local adaptation across a steep latitudinal gradient-as an ideal system to explore how genome structure influences adaptation under varying levels of gene flow, using a newly available reference genome and multiple recombination maps. Analyzing 168 genomes from four populations, we found a continuum of genome-wide differentiation increasing from south to north, reflecting higher connectivity among southern populations and reduced gene flow at northern latitudes. With increasing gene flow, the number and clustering of FST outlier loci also increased, with differentiated loci found exclusively within large haploblocks harboring inversions and smaller peaks overlapping putative centromeric regions. Notably, sequence divergence was only evident in inversions, supporting their role in adaptive divergence with gene flow, whereas centromeric regions appeared differentiated because of low recombination and diversity, with no indication of elevated divergence. Our results support the hypothesis that clustered genomic architectures evolve with high gene flow and enhance our understanding of how inversions and centromeres are linked to different evolutionary processes.

Gene Flow

Rapid emergence of novel antigenic and genetic variants of equine infectious anemia virus during persistent infection.

Previous results from our laboratory have demonstrated that equine infectious anemia virus displays structural variations in its surface glycoproteins and RNA genome during passage and chronic infections in experimentally infected Shetland ponies (Montelaro et al., J. Biol. Chem. 259:10539-10544, 1984; Payne et al., J. Gen. Virol. 65:1395-1399, 1984). The present study was undertaken to obtain an antigenic and biochemical characterization of equine infectious anemia virus isolates recovered from an experimentally infected pony during sequential disease episodes, each separated by intervals of only 4 to 8 weeks. The virus isolates could be distinguished antigenically by neutralization assays with serum from the infected pony and by Western blot analysis with a monoclonal antibody against the major surface glycoprotein gp90, thus demonstrating that novel antigenic variants of equine infectious anemia virus predominate during each clinical episode. The respective virion glycoproteins displayed different electrophoretic mobilities on sodium dodecyl sulfate-polyacrylamide gels, indicating structural variation. Tryptic peptide and glycopeptide maps of the viral proteins of each virus isolate revealed biochemical alterations involving amino acid sequence and glycosylation patterns in the virion surface glycoproteins gp90 and gp45. In contrast, no structural variation was observed in the internal viral proteins pp15, p26, and p9 from any of the four virus isolates. Oligonucleotide mapping experiments revealed similar but unique RNase T1-resistant oligonucleotide fingerprints of the RNA genomes of each of the virus isolates. Localization of altered oligonucleotides for one virus isolate placed two of three unique oligonucleotides within the predicted env gene region of the genome, perhaps correlating with the structural variation observed in the envelope glycoproteins. Thus these results support the concept that equine infectious anemia virus is indeed capable of relatively rapid genomic variations during replication, some of which result in altered glycoprotein structures and antigenic variants which are responsible for the unique periodic disease nature observed in persistently infected animals. The findings of envelope specific differences in isolates of visna virus and of human T-cell lymphotropic virus III (acquired immune deficiency syndrome-related virus) suggest that this variation may be a common characteristic of the subfamily Lentivirinae.

Animals

Globin gene structure and the nature of mutation.

In this chapter I have tried to relate the salient features of globin gene structure, genomic organization, normal variation, and mutations affecting gene expression. The lessons learned from the normal beta-globin gene cluster and mutations producing beta-thalassemia should be highly applicable to studies of other inherited diseases. As more and more gene probes become available, which have relevance to the study of human disease, the striking extent of genetic heterogeneity producing single gene disorders of man will be illuminated.

Gene Expression Regulation

Low-pass whole-genome sequencing reveals genomic diversity and ecotype-specific adaptation in indigenous Tigrayan chickens.

Indigenous chickens play a critical role in food security and climate resilience in smallholder systems, yet their genomic diversity and adaptive potential remain insufficiently characterised. This study employed low-pass whole-genome sequencing (LP-WGS; 0.2-1.99&#xd7;) to investigate genomic diversity, population structure, inbreeding and candidate environment-associated genomic variation in 33 chickens from highland, midland, and lowland agroecologies in the Tigray region of northern Ethiopia. After imputation and stringent filtering, 23.4 million high-confidence SNPs were retained, including&#x2009;~&#x2009;17% novel variants, indicating substantial uncharacterised genetic diversity in these populations. SNP density (13.8&#x2009;&#xb1;&#x2009;8.6 SNPs/kb) was comparable to values reported from high-coverage Ethiopian chicken datasets, demonstrating the suitability of LP-WGS for population genomics in resource-limited settings. Marked differences in genomic diversity were observed among ecotypes: midland chickens showed the highest nucleotide diversity (&#x3c0;&#x2009;=&#x2009;0.00267), followed by lowland (&#x3c0;&#x2009;=&#x2009;0.00233), whereas highland chickens showed the lowest diversity (&#x3c0;&#x2009;=&#x2009;0.00203) and elevated genomic inbreeding (FROH and FHOM &#x2248; 0.18). Population structure analyses revealed clear genetic separation among ecotypes. PCA (13.91% variation explained) distinguished lowland chickens along PC1 and separated highland from midland along PC2, while ADMIXTURE and FST patterns supported three major ancestral genomic backgrounds. Functional annotation of private missense variants uncovered distinct adaptive signatures reflecting the contrasting agroecological conditions. Highland chickens showed enrichment of candidate genes potentially involved in physiological processes relevant to high-altitude environments, including cold response, angiogenesis, cardiovascular regulation and metabolic homeostasis (eg., PARP1, ACOX2, ITGB3, EDNRB, SOX8, and SOX10). Midland chickens exhibited candidate signals of selection in genes with known roles in innate antiviral immunity, bacterial defence and inflammatory regulation (eg., BAK1, CLSTN1, CYSLTR1, CYSLTR2, CXCR7, GIPR, DSCAM, GDAP1, TLR3, TLR4, TLR7, IFIH1, ADORA1, EPHB1, and TMPRSS2). Lowland chickens displayed candidate variants associated with heat-stress response, DNA damage repair, oxidative balance and cardiovascular support under extreme temperatures (e.g., MLH1, BDKRB1, GPR19, FLT1, CCL18, TGM2, and RAMP3). Overall, the results indicate substantial genomic differentiation among ecotypes and suggest candidate environment-associated genetic divergence across Tigray's diverse agroecological zones. These populations may represent important reservoirs of adaptive genetic variation for climate-resilient poultry breeding, warranting further functional validation and conservation-oriented management.

Animals

Whole-Genome Sequencing Reveals Population Structure, Genetic Diversity, and Selection Signatures in Kazakh Dromedary and Bactrian Camels.

Understanding the genomic basis of environmental adaptation is essential for the conservation and genetic improvement of domestic camels. In this study, we investigated the population structure, genetic diversity, and genomic variation potentially associated with environmental adaptation of Kazakh dromedary and Bactrian camels using whole-genome sequencing. Whole-genome sequencing data were generated for Kazakh camels (15 dromedaries and 16 Bactrian camels) and integrated with 131 publicly available genomes representing camel populations from the Arabian Peninsula, Iran, Xinjiang, Inner Mongolia, and Mongolian wild camels. Population structure, genetic diversity, and genome-wide selection were evaluated using principal component analysis, ADMIXTURE, nucleotide diversity, linkage disequilibrium, runs of homozygosity, genomic inbreeding (FROH), and selection scans based on FST, &#x3b8;&#x3c0; ratio, and XP-EHH. Population genomic analyses revealed clear differentiation between dromedary and Bactrian camels, whereas Kazakh camel populations exhibited higher nucleotide diversity (&#x3b8;&#x3c0; = 1.307-1.551 &#xd7; 10-3), and lower genomic inbreeding (median FROH: 0.037-0.056) than Arabian populations. Genome-wide selection analyses identified MC4R as the prominent candidate gene in Kazakh dromedaries and RYR1 as a prominent candidate gene in Kazakh Bactrian camels. Functional enrichment analyses highlighted pathways related to energy metabolism, thermogenesis, calcium signaling, skeletal muscle function, mitochondrial activity, and oxidative stress response. These findings provide new insights into genomic variation potentially associated with environmental adaptation in Kazakh camels and offer valuable genomic resources for future conservation, breeding, and evolutionary studies.

MC4R

Hidden genomic structure and widespread structural polymorphism across environmental gradients in the spiny sea star Marthasterias glacialis.

Genomic regions of reduced recombination can preserve linkage among co-adapted alleles, facilitating local adaptation despite high connectivity. Such regions-often generated by chromosomal inversions-may be especially important in highly dispersive marine taxa yet remain poorly documented in echinoderms. Here, we combined a chromosome-level reference genome with genome-wide ddRAD-seq from 296 Marthasterias glacialis individuals across 19 Atlantic-Mediterranean locations to quantify population structure and scan for recombination-suppressed haploblocks. Genome-wide neutral markers showed significant population differentiation together with evidence of high connectivity, revealed by the presence of inter-ecoregion migrants. Additionally, we identified 16 polymorphic haploblocks with patterns consistent with putative chromosomal inversions spanning 18.6% of the genome. Haploblock haplotypes were strongly environmentally and geographically structured and contained genes with key functions in stress response, osmoregulation and thermal tolerance. Haplotype distributions also paralleled previously described mitochondrial lineages despite nuclear gene flow, consistent with a model of ancient divergence followed by secondary contact. Overall, our results suggest a role for widespread structural polymorphism in adaptive differentiation in Echinodermata, providing a framework for linking echinoderm genome rearrangements to ecological divergence. Marthasterias glacialis thus emerges as a promising system to explore how structural variation contributes to adaptation and genome evolution in highly dispersive organisms.

Animals

Long-read sequencing to interrogate strain-level variation among adherent-invasive Escherichia coli isolated from human intestinal tissue.

Adherent-invasive Escherichia coli (AIEC) is a pathovar linked to inflammatory bowel diseases (IBD), especially Crohn's disease, and colorectal cancer. AIEC are genetically diverse, and in the absence of a universal molecular signature, are defined by in vitro functional attributes. The relative ability of difference AIEC strains to colonize, persist, and induce inflammation in an IBD-susceptible host is unresolved. To evaluate strain-level variation among tissue-associated E. coli in the intestines, we develop a long-read sequencing approach to identify AIEC by strain that excludes host DNA. We use this approach to distinguish genetically similar strains and assess their fitness in colonizing the intestine. Here we have assembled complete genomes using long-read nanopore sequencing for a model AIEC strain, NC101, and seven strains isolated from the intestinal mucosa of Crohn's disease and non-Crohn's tissues. We show these strains can colonize the intestine of IBD susceptible mice and induce inflammatory cytokines from cultured macrophages. We demonstrate that these strains can be quantified and distinguished in the presence of 99.5% mammalian DNA and from within a fecal population. Analysis of global genomic structure and specific sequence variation within the ribosomal RNA operon provides a framework for efficiently tracking strain-level variation of closely-related E. coli and likely other commensal/pathogenic bacteria impacting intestinal inflammation in experimental settings and IBD patients.

Animals

Whole-genome sequencing implicates rare, low-frequency and structural non-coding variation at the SCN5A locus in Brugada syndrome.

Brugada syndrome (BrS) is an inherited cardiac condition characterized by a hallmark ECG pattern and an increased risk of sudden cardiac death. Central to the aetiology of BrS, the SCN5A region harbours both common non-coding risk variants and rare coding variants that are causative in approximately 20% of patients. However, rare non-coding genetic variation in this region remains largely unexplored. Here, we used whole-genome sequencing (WGS) of 752 European-ancestry BrS cases and 1,827 ancestry-matched controls to identify BrS-associated rare non-coding genetic variation at the SCN5A locus. Sliding-window and cis-regulatory element (CRE)-based rare-variant aggregate testing implicated three conserved CREs, including a dense aggregation of case singleton variants within a 178 bp enhancer in intron 17 of SCN5A which replicated in an independent BrS cohort. Prioritised BrS-associated rare and low-frequency non-coding variants within these elements were predicted to alter cardiac transcription factor motifs, and altered CRE activity in hiPSC-CM luciferase assays or were associated with BrS-relevant ECG endophenotypes in the UK Biobank. Single-variant analysis across the region identified a Bonferroni-significant five-fold case-enriched low-frequency variant within a known CRE in intron 1 of SCN5A, which replicated, was associated with slower cardiac conduction in the UK Biobank and accounted for part of the BrS GWAS signal at this locus. Structural variant analyses identified a 10.5 kb deletion upstream of SCN5A in a BrS case that encompassed a cardiac CRE and reduced sodium current density in a hiPSC-CM model, as well as a 6 kb BrS-enriched retrotransposon insertion in SCN5A that appeared to underlie part of the GWAS signal in this region. Together, these findings implicate rare and low-frequency non-coding variation at the SCN5A locus in BrS susceptibility and demonstrate the value of targeted WGS analysis of key disease loci.

Journal Article

In silico prediction of the impact of genomic variations in the small conductance calcium activated potassium channel SK3 structure and function.

The small-conductance calcium-activated potassium channel SK3, encoded by the KCNN3 gene, plays a critical role in regulating dopaminergic neuron (DN) firing patterns by modulating after hyperpolarization currents. SK3 dysfunction has been implicated in neuropsychiatric and neurodegenerative disorders. We analyzed structural and functional consequences of KCNN3 splicing and genetic variation. Alternative splicing variants of the KCNN3 gene were retrieved from the Ensembl database and aligned using T-Coffee, manually inspected and curated. Protein domains were identified with Pfam 35.0, SMART 9.0, and InterPro 98.0, and visualized. An AlphaFold2 model of SK3 full-length protein (UniProt: Q9UGI6) used as reference and structural models of its splicing variants were predicted with ColabFold. Functional domains (S1-S6 transmembrane helices, H5 pore loop, and calmodulin-binding) were defined and superimposed onto the AlphaFold2 reference. Domain integrity was assessed based on completeness of all expected residue indices within each functional region. SNPs and CNVs across all coding KCNN3 splicing variants were analyzed, classified, and filtered to isolate pathogenic variants prioritizing non-synonymous amino acid substitutions. Differential variant impacts across splicing isoforms were assessed by mapping variant positions to individual transcript protein sequences and used to predict functional consequences. Two long and two short splicing variants are known. Short variants lack the motif required for potassium channels. Pathogenic variants result from missense mutations resulting in amino acid substitutions. In all cases, the consequential effects depend on the specific location and role of the amino acid being changed.

SK3 channels