The UCSC Genome Browser database: 2023 update

kmerDB: A database encompassing the set of genomic and proteomic sequence information for each species

The decrease in sequencing expenses has facilitated the creation of reference genomes and proteomes for an expanding array of organisms. Nevertheless, no established repository that details organism-specific genomic and proteomic sequences of specific lengths, referred to as kmers, exists to our knowledge. In this article, we present kmerDB, a database accessible through an interactive web interface that provides kmer-based information from genomic and proteomic sequences in a systematic way. kmerDB currently contains 202,340,859,107 base pairs and 19,304,903,356 amino acids, spanning 54,039 and 21,865 reference genomes and proteomes, respectively, as well as 6,905,362 and 149,305,183 genomic and proteomic species-specific sequences, termed quasi-primes. Additionally, we provide access to 5,186,757 nucleic and 214,904,089 peptide sequences absent from every genome and proteome, termed primes. kmerDB features a user-friendly interface offering various search options and filters for easy parsing and searching. The service is available at: www.kmerdb.com.

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Comprehensive pan-cancer analysis of KLRB1-CLEC2D pair and identification of small molecule inhibitors to disrupt their interaction

The interplay between immune checkpoints KLRB1 and CLEC2D is crucial for tumor progression and immune evasion, yet the interaction dynamics are not fully understood. This study aims to elucidate the interaction across various cancers and identify small molecule inhibitors that can disrupt it. We perform a comprehensive pan-cancer analysis of the KLRB1-CLEC2D pair, including mRNA expression patterns, pathological stages, survival outcomes, and single-cell omics, immune infiltration, copy number variations, and DNA methylation profiles. Our findings reveal a consistently higher CLEC2D/KLRB1 ratio in most cancer types compared to normal tissues, and this ratio also increased with advancing pathological stages. Lower KLRB1 expression correlated with higher mortality in most cancers, opposite to CLEC2D. Expression variations were attributed to differential lymphocyte infiltration, CNV, and DNA methylation. Structure-based virtual screening analysis identified compounds including forsythiaside A and RGD peptides as effective inhibitors of the KLRB1-CLEC2D interaction, validated through microscale thermophoresis. This research advances understanding of the KLRB1-CLEC2D interaction within the tumor microenvironment and introduces novel therapeutic strategies to modulate this interaction.

Minimal phenotypes in transgenic mice with the human LOXL1/LOXL1-AS1 locus associated with exfoliation glaucoma

Exfoliation syndrome is a leading cause of secondary glaucoma worldwide. Among the risk-factors for exfoliation syndrome and exfoliation glaucoma that have been investigated, a genetic association with 15q24.1 is among the most striking. The leading candidates for the causal gene at this locus are LOXL1 and/or LOXL1-AS1, but studies have not yet coalesced in establishing, or ruling out, either candidate. Here, we contribute to studies of the 15q24.1 locus by making a partially humanized mouse model in which 166 kb of human genomic DNA from the 15q24.1 locus was introduced into the mouse genome via BAC transgenesis (B6-Tg(RP11-71M11)Andm). Transgenic expression of human genes in the BAC was only detectable for LOXL1-AS1. One cohort of 34 mice (21 experimental hemizygotes and 13 non-carrier control littermates) was assessed by slit-lamp exams and SD-OCT imaging at early (1-2 months) and mid (4-5 months) time points; fundus exams were performed at 5 months of age. A second smaller cohort (3 hemizygotes) were aged extensively (>12 months) to screen for overt abnormalities. Across all genotypes and ages, 136 slit-lamp exams, 128 SD-OCT exams, and 42 fundus exams detected no overt indices of exfoliation syndrome. Quantitatively, small, but statistically significant, age-related declines in ganglion cell complex thickness and total retinal thickness were detected in the hemizygotes at 4 months of age. Overall, this study demonstrates complexity in gene regulation from the 15q24.1 locus and suggests that LOXL1-AS1 is unlikely to be a monogenic cause of exfoliation syndrome but may contribute to glaucomatous retinal damage.

Control of epigenomic landscape and development of fetal male germ cells through L-serine metabolism

Sex-specific metabolic characteristics emerge in the mouse germ line after reaching the genital ridges around embryonic day 10.5, coinciding with sexual differentiation. However, the impact of such metabolic characteristics on germ cell development remains unclear. In this study, we observed the specific upregulation in male fetal germ cells of D-3-phosphoglycerate dehydrogenase (PHGDH), the primary enzyme in the serine-glycine-one-carbon metabolism, along with an increase in a downstream metabolite, S-adenosylmethionine (SAM), crucial for protein and nucleic acid methylation. Inhibiting PHGDH in fetal testes resulted in reduced SAM levels in germ cells, accompanied by increases in the number of mouse vasa homolog (MVH/VASA)-positive germ cells and the promyelocytic leukemia zinc finger (PLZF)-positive undifferentiated spermatogonia ratio. Furthermore, PHGDH inhibition led to a decrease in the methylation of histone H3 and DNA, resulting in aberrations in gene expression profiles. In summary, our findings underscore the significant role of certain metabolic mechanisms in the development of male germ cells.

A deep intronic splice-altering AIRE variant causes APECED syndrome through antisense oligonucleotide-targetable pseudoexon inclusion

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED) is a life-threatening monogenic autoimmune disorder primarily caused by biallelic deleterious variants in the autoimmune regulator (AIRE) gene. We prospectively evaluated 104 patients with clinically diagnosed APECED syndrome and identified 17 patients (16%) from 14 kindreds lacking biallelic AIRE variants in exons or flanking intronic regions; 15 had Puerto Rican ancestry. Through whole-genome sequencing, we identified a deep intronic AIRE variant (c.1504-818 G>A) cosegregating with the disease in all 17 patients. We developed a culture system of AIRE-expressing primary patient monocyte-derived dendric cells and demonstrated that c.1504-818 G>A creates a cryptic splice site and activates inclusion of a 109-base pair frame-shifting pseudoexon. We also found low-level AIRE expression in patient-derived lymphoblastoid cell lines (LCLs) and confirmed pseudoexon inclusion in independent extrathymic AIRE-expressing cell lines. Through protein modeling and transcriptomic analyses of AIRE-transfected human embryonic kidney 293 and thymic epithelial cell 4D6 cells, we showed that this variant alters the carboxyl terminus of the AIRE protein, abrogating its function. Last, we developed an antisense oligonucleotide (ASO) that reversed pseudoexon inclusion and restored the normal AIRE transcript sequence in LCLs. Thus, our findings revealed c.1504-818 G>A as a founder APECED-causing AIRE variant in the Puerto Rican population and uncovered pseudoexon inclusion as an ASO-reversible genetic mechanism underlying APECED.

Nearby and non-nested genes in the human genome have more similar genotype tissue expression

Neighboring genes within a shared promoter arrangement (i.e. opposite direction with the neighboring ends as the transcriptional start sites) are expected to have a high similarity in genotype tissue expression due to the potential overlap in the promoter region. This raises the question of whether similarity in expression profiles depends on orientation of the neighboring genes and whether there exist thresholds of locality where the similarity diminishes. Thus, in this work, we compared genotype tissue expression profiles at different genomic orientations and localities. Interestingly, there exist gene pairs in the human genome with very high or low expression similarity. Shorter chromosomes tend to have more similarly expressed genes. Also, a cluster of 3 adjacent genes within the average range of 20 to 60 kilobase pairs can have very similar expression profiles regardless of their orientations. However, when genes are nested and in opposite orientations, a lower than expected similarity was observed. Lastly, in cases where genotype tissue expression data does not exist or have low read counts (e.g. non-coding RNA), our identified influencing range can be a first estimate of the genotype tissue expression.

Phenotype wide association study links bronchopulmonary dysplasia with eosinophilia in children

Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth. Despite this, genetic drivers of BPD are poorly understood. The objective of this study is to better understand the impact of single nucleotide polymorphisms (SNPs) previously associated with BPD by examining associations with other phenotypes. We drew pediatric subjects from the biorepository at the Center for Applied Genomics to identify associations between these SNPs and 2,146 imputed phenotypes. Methylation data, external cohorts, and in silico validation methods were used to corroborate significant associations. We identified 60 SNPs that were previously associated with BPD. We found a significant association between rs3771150 and rs3771171 and mean eosinophil percentage in a European cohort of 6,999 patients and replicated this in external cohorts. Both SNPs were also associated with asthma, COPD and FEV1/FVC ratio. These SNPs displayed associations with methylation probes and were functionally linked to ST2 (IL1RL1) levels in blood and lung tissue. Our findings support a genetic justification for the epidemiological link between BPD and asthma. Given the well-established link between ST2 and type 2 inflammation in asthma, these findings provide a rationale for future studies exploring the role of type 2 inflammation in the pathogenesis of BPD.

Stress increases sperm respiration and motility in mice and men

Semen quality and fertility has declined over the last 50 years, corresponding to ever-increasing environmental stressors. However, the cellular mechanisms involved and their impact on sperm functions remain unknown. In a repeated sampling human cohort study, we identify a significant effect of prior perceived stress to increase sperm motility 2-3 months following stress, timing that expands upon our previous studies revealing significant stress-associated changes in sperm RNA important for fertility. We mechanistically examine this post-stress timing in mice using an in vitro stress model in the epididymal epithelial cells responsible for sperm maturation and find 7282 differentially H3K27me3 bound DNA regions involving genes critical for mitochondrial and metabolic pathways. Further, prior stress exposure significantly changes the composition and size of epithelial cell-secreted extracellular vesicles that when incubated with mouse sperm, increase mitochondrial respiration and sperm motility, adding to our prior work showing impacts on embryo development. Together, these studies identify a time-dependent, translational signaling pathway that communicates stress experience to sperm, ultimately affecting reproductive functions.

Systemic lupus erythematosus genetics: insights into pathogenesis and implications for therapy

Systemic lupus erythematosus (SLE) is a prime example of how the interplay between genetic and environmental factors can trigger systemic autoimmunity, particularly in young women. Although clinical disease can take years to manifest, risk is established by the unique genetic makeup of an individual. Genome-wide association studies have identified almost 200 SLE-associated risk loci, yet unravelling the functional effect of these loci remains a challenge. New analytic tools have enabled researchers to delve deeper, leveraging DNA sequencing and cell-specific and immune pathway analysis to elucidate the immunopathogenic mechanisms. Both common genetic variants and rare non-synonymous mutations can interact to increase SLE risk. Notably, variants strongly associated with SLE are often located in genome super-enhancers that regulate MHC class II gene expression. Additionally, the 3D conformations of DNA and RNA contribute to genome regulation and innate immune system activation. Improved therapies for SLE are urgently needed and current and future knowledge from genetic and genomic research should provide new tools to facilitate patient diagnosis, enhance the identification of therapeutic targets and optimize testing of agents.

High-risk screening for late-onset Pompe disease in China: An expanded multicenter study

Late-onset Pompe disease (LOPD) is caused by a genetic deficiency of the lysosomal enzyme acid alpha-glucosidase (GAA), leading to progressive limb-girdle weakness and respiratory impairment. The insidious onset of non-specific early symptoms often prohibits timely diagnosis. This study aimed to validate the high-risk screening criteria for LOPD in the Chinese population. A total of 726 patients were included, including 96 patients under 14 years of age. Dried blood spots (DBS) and tandem mass spectrometry (MS/MS) were employed to evaluate serum GAA activity. Forty-four patients exhibited a decreased GAA activity, 16 (2.2%) of which were confirmed as LOPD by genetic testing. Three previously unreported GAA mutations were also identified. The median diagnostic delay was shortened to 3 years, which excelled the previous retrospective studies. At diagnosis, most patients exhibited impaired respiratory function and/or limb-girdle weakness. Elevated serum creatine kinase (CK) levels were more frequently observed in patients who manifested before age 16. Overall, high-risk screening is a feasible and efficient method to identify LOPD patients at an early stage. Patients over 1 year of age with either weakness in axial and/or proximal limb muscles, or unexplained respiratory distress shall be subject to GAA enzymatic test, while CK levels above 2 times the upper normal limit shall be an additional criterion for patients under 16. This modified high-risk screening criteria for LOPD requires further validation in larger Chinese cohorts.

Deciphering regulation of FOXP3 expression in human conventional T cells

FOXP3 is a lineage-defining transcription factor that controls differentiation and maintenance of suppressive function of regulatory T cells (Tregs). Foxp3 is exclusively expressed in Tregs in mice. However, in humans, FOXP3 is not only constitutively expressed in Tregs; it is also transiently expressed in stimulated CD4+CD25- conventional T cells (Tconvs)1-3. Mechanisms governing the expression of FOXP3 in human Tconvs are not understood. Here, we performed CRISPR interference (CRISPRi) screens using a 15K-member gRNA library tiling 39 kb downstream of the FOXP3 transcriptional start site (TSS) to 85 kb upstream of the TSS in Treg and Tconvs. The FOXP3 promoter and conserved non-coding sequences (CNS0, CNS1, CNS2 and CNS3), characterized as enhancer elements in murine Tregs, were required for maintenance of FOXP3 in human Tregs. In contrast, FOXP3 in human Tconvs depended on regulation at CNS0 and a novel Tconv-specific noncoding sequence (TcNS+) located upstream of CNS0. Arrayed validations of these sites identified an additional repressive cis-element overlapping with the PPP1R3F promoter (TcNS-). Pooled CRISPR knockouts revealed multiple transcription factors required for proper expression of FOXP3 in Tconvs, including GATA3, STAT5, IRF4, ETS1 and DNA methylation-associated regulators DNMT1 and MBD2. Analysis of ChIP-seq and ATAC-seq paired with knock-out (KO) of GATA3, STAT5, IRF4, and ETS1 revealed regulation of CNS0 and TcNS+ accessibility. Collectively, this work identified Treg-shared and Tconv-specific cis-elements and the trans-factors that interact with them, building a network of regulators controlling FOXP3 expression in human Tconvs.

An epigenome-wide study of selenium status and DNA methylation in the Strong Heart Study

BACKGROUND

Selenium (Se) is an essential nutrient linked to adverse health endpoints at low and high levels. The mechanisms behind these relationships remain unclear and there is a need to further understand the epigenetic impacts of Se and their relationship to disease. We investigated the association between urinary Se levels and DNA methylation (DNAm) in the Strong Heart Study (SHS), a prospective study of cardiovascular disease (CVD) among American Indians adults.

METHODS

Selenium concentrations were measured in urine (collected in 1989-1991) using inductively coupled plasma mass spectrometry among 1,357 participants free of CVD and diabetes. DNAm in whole blood was measured cross-sectionally using the Illumina MethylationEPIC BeadChip (850 K) Array. We used epigenome-wide robust linear regressions and elastic net to identify differentially methylated cytosine-guanine dinucleotide (CpG) sites associated with urinary Se levels.

RESULTS

The mean (standard deviation) urinary Se concentration was 51.8 (25.1) μg/g creatinine. Across 788,368 CpG sites, five differentially methylated positions (DMP) (hypermethylated: cg00163554, cg18212762, cg11270656, and hypomethylated: cg25194720, cg00886293) were significantly associated with Se in linear regressions after accounting for multiple comparisons (false discovery rate p-value: 0.10). The top hypermethylated DMP (cg00163554) was annotated to the Disco Interacting Protein 2 Homolog C (DIP2C) gene, which relates to transcription factor binding. Elastic net models selected 425 hypo- and hyper-methylated DMPs associated with urinary Se, including three sites (cg00163554 [DIP2C], cg18212762 [MAP4K2], cg11270656 [GPIHBP1]) identified in linear regressions.

CONCLUSIONS

Urinary Se was associated with minimal changes in DNAm in adults from American Indian communities across the Southwest and the Great Plains in the United States, suggesting that other mechanisms may be driving health impacts. Future analyses should explore other mechanistic biomarkers in human populations, determine these relationships prospectively, and investigate the potential role of differentially methylated sites with disease endpoints.

Genetic drivers of age-related changes in urinary magnesium excretion

Although age-dependent alterations in urinary magnesium (Mg2+) excretion have been described, the underlying mechanism remains elusive. As heritability significantly contributes to variations in urinary Mg2+ excretion, we measured urinary Mg2+ excretion at different ages in a cohort of genetically variable Diversity Outbred (DO) mice. Compared with animals aged 6 mo, an increase in Mg2+ excretion was observed at 12 and 18 mo. Quantitative trait locus (QTL) analysis revealed an association of a locus on chromosome 10 with Mg2+ excretion at 6 mo of age, with Oit3 (encoding oncoprotein-induced transcript 3; OIT3) as our primary candidate gene. To study the possible role of OIT3 in renal Mg2+ handling, we generated and characterized Oit3 knockout (Oit3-/-) mice. Although a slightly lower serum Mg2+ concentration was present in male Oit3-/- mice, this effect was not observed in female Oit3-/- mice. In addition, urinary Mg2+ excretion and the expression of renal magnesiotropic genes were unaltered in Oit3-/- mice. For animals aged 12 and 18 mo, QTL analysis revealed an association with a locus on chromosome 19, which contains the gene encoding TRPM6, a known Mg2+ channel involved in renal Mg2+ reabsorption. Comparison with RNA sequencing (RNA-Seq) data revealed that Trpm6 mRNA expression is inversely correlated with the QTL effect, implying that TRPM6 may be involved in age-dependent changes in urinary Mg2+ excretion in mice. In conclusion, we show here that variants in Oit3 and Trpm6 are associated with urinary Mg2+ excretion at distinct periods of life, although OIT3 is unlikely to affect renal Mg2+ handling.NEW & NOTEWORTHY Aging increased urinary magnesium (Mg2+) excretion in mice. We show here that variation in Oit3, a candidate gene for the locus associated with Mg2+ excretion in young mice, is unlikely to be involved as knockout of Oit3 did not affect Mg2+ excretion. Differences in the expression of the renal Mg2+ channel TRPM6 may contribute to the variation in urinary Mg2+ excretion in older mice.

Deletions in the CDKL5 5' untranslated region lead to CDKL5 deficiency disorder

Pathogenic variants in the cyclin-dependent kinase-like 5 (CDKL5) gene are associated with CDKL5 deficiency disorder (CDD), a severe X-linked developmental and epileptic encephalopathy. Deletions affecting the 5' untranslated region (UTR) of CDKL5, which involve the noncoding exon 1 and/or alternatively spliced first exons (exons 1a-e), are uncommonly reported. We describe genetic and phenotypic characteristics for 15 individuals with CDKL5 partial gene deletions affecting the 5' UTR. All individuals presented characteristic features of CDD, including medically refractory infantile-onset epilepsy, global developmental delay, and visual impairment. We performed RNA sequencing on fibroblast samples from three individuals with small deletions involving exons 1 and/or 1a/1b only. Results demonstrated reduced CDKL5 mRNA expression with no evidence of expression from alternatively spliced first exons. Our study broadens the genotypic spectrum for CDD by adding to existing evidence that deletions affecting the 5' UTR of the CDKL5 gene are associated with the disorder. We propose that smaller 5' UTR deletions may require additional molecular testing approaches such as RNA sequencing to determine pathogenicity.

SpliceVarDB: A comprehensive database of experimentally validated human splicing variants

Variants that alter gene splicing are estimated to comprise up to a third of all disease-causing variants, yet they are hard to predict from DNA sequencing data alone. To overcome this, many groups are incorporating RNA-based analyses, which are resource intensive, particularly for diagnostic laboratories. There are thousands of functionally validated variants that induce mis-splicing; however, this information is not consolidated, and they are under-represented in ClinVar, which presents a barrier to variant interpretation and can result in duplication of validation efforts. To address this issue, we developed SpliceVarDB, an online database consolidating over 50,000 variants assayed for their effects on splicing in over 8,000 human genes. We evaluated over 500 published data sources and established a spliceogenicity scale to standardize, harmonize, and consolidate variant validation data generated by a range of experimental protocols. According to the strength of their supporting evidence, variants were classified as "splice-altering" (∼25%), "not splice-altering" (∼25%), and "low-frequency splice-altering" (∼50%), which correspond to weak or indeterminate evidence of spliceogenicity. Importantly, 55% of the splice-altering variants in SpliceVarDB are outside the canonical splice sites (5.6% are deep intronic). These variants can support the variant curation diagnostic pathway and can be used to provide the high-quality data necessary to develop more accurate in silico splicing predictors. The variants are accessible through an online platform, SpliceVarDB, with additional features for visualization, variant information, in silico predictions, and validation metrics. SpliceVarDB is a very large collection of splice-altering variants and is available at https://splicevardb.org.

Integrated clinical risk prediction of type 2 diabetes with a multifactorial polygenic risk score

Combining information from multiple GWASs for a disease and its risk factors has proven a powerful approach for development of polygenic risk scores (PRSs). This may be particularly useful for type 2 diabetes (T2D), a highly polygenic and heterogeneous disease where the additional predictive value of a PRS is unclear. Here, we use a meta-scoring approach to develop a metaPRS for T2D that incorporated genome-wide associations from both European and non-European genetic ancestries and T2D risk factors. We evaluated the performance of this metaPRS and benchmarked it against existing genome-wide PRS in 620,059 participants and 50,572 T2D cases amongst six diverse genetic ancestries from UK Biobank, INTERVAL, the All of Us Research Program, and the Singapore Multi-Ethnic Cohort. We show that our metaPRS was the most powerful PRS for predicting T2D in European population-based cohorts and had comparable performance to the top ancestry-specific PRS, highlighting its transferability. In UK Biobank, we show the metaPRS had stronger predictive power for 10-year risk than all individual risk factors apart from BMI and biomarkers of dysglycemia. The metaPRS modestly improved T2D risk stratification of QDiabetes risk scores for 10-year risk prediction, particularly when prioritising individuals for blood tests of dysglycemia. Overall, we present a highly predictive and transferrable PRS for T2D and demonstrate that the potential for PRS to incrementally improve T2D risk prediction when incorporated into UK guideline-recommended screening and risk prediction with a clinical risk score.

Trained immunity is regulated by T cell-induced CD40-TRAF6 signaling

Trained immunity is characterized by histone modifications and metabolic changes in innate immune cells following exposure to inflammatory signals, leading to heightened responsiveness to secondary stimuli. Although our understanding of the molecular regulation of trained immunity has increased, the role of adaptive immune cells herein remains largely unknown. Here, we show that T cells modulate trained immunity via cluster of differentiation 40-tissue necrosis factor receptor-associated factor 6 (CD40-TRAF6) signaling. CD40-TRAF6 inhibition modulates functional, transcriptomic, and metabolic reprogramming and modifies histone 3 lysine 4 trimethylation associated with trained immunity. Besides in vitro studies, we reveal that single-nucleotide polymorphisms in the proximity of CD40 are linked to trained immunity responses in vivo and that combining CD40-TRAF6 inhibition with cytotoxic T lymphocyte antigen 4-immunoglobulin (CTLA4-Ig)-mediated co-stimulatory blockade induces long-term graft acceptance in a murine heart transplantation model. Combined, our results reveal that trained immunity is modulated by CD40-TRAF6 signaling between myeloid and adaptive immune cells and that this can be leveraged for therapeutic purposes.

Multiple DNA Viruses and HPV Integration in Inverted Papilloma and Associated Sinonasal Carcinoma

OBJECTIVES

Sinonasal inverted papilloma (IP) has a locally destructive growth pattern, can relapse, and can undergo malignant transformation (IP-associated sinonasal squamous cell carcinoma (IP-SNSCC)). Human papillomaviruses (HPV)-6 and -16 are frequently detected in IPs. To clarify the possible roles of other DNA viruses in IPs, we explored viruses not studied in this context before. With the setting of pre- and post-malignant transformation samples, we investigated HPV genomes in depth to assess the integration of HPV into the human genome and the presence of minor intratypic variants.

MATERIALS AND METHODS

We analyzed 35 IP samples representing 28 individuals, of which six had IP-SNSCC. For virus screening, we applied qPCR to detect 16 different DNA viruses in three virus families, comprising herpesviruses, parvoviruses, and polyomaviruses. In addition, targeted next generation sequencing (NGS) was used for detailed HPV analysis.

RESULTS

We detected herpes-, parvo-, and polyomaviruses in 13/28 (46%) patients, with codetections of multiple viruses in six (21%) patients. NGS revealed HPV16 DNA in 2/6 IP-SNSCC and in their respective earlier benign IP samples, as well as in a plasma sample from one of these patients. HPV6 was detected in two IP samples without subsequent malignant transformation. We identified sequence reads containing junctions of HPV6 and HPV16 and host genome suggestive of viral integration. HPV6 and HPV16 minor intratypic variants were present across pre- and post-malignant transformation, with mostly nonsynonymous mutations.

CONCLUSIONS

Multiple DNA viruses were present in IPs. HPV16 was detected only in IP-SNSCCs or in tumors that later underwent malignant transformation.

LEVEL OF EVIDENCE

3 Laryngoscope, 2024.

Decoding the Role of NEIL1 Gene in DNA Repair and Lifespan: A Literature Review with Bioinformatics Analysis

Longevity, the length of an organism's lifespan, is impacted by environmental factors, metabolic processes, and genetic determinants. The base excision repair (BER) pathway is crucial for maintaining genomic integrity by repairing oxidatively modified base lesions. Nei-like DNA Glycosylase 1 (NEIL1), part of the BER pathway, is vital in repairing oxidative bases in G-rich DNA regions, such as telomeres and promoters. Hence, in this comprehensive review, it have undertaken a meticulous investigation of the intricate association between NEIL1 and longevity. The analysis delves into the multifaceted aspects of the NEIL1 gene, its various RNA transcripts, and the diverse protein isoforms. In addition, a combination of bioinformatic analysis is conducted to identify NEIL1 mutations, transcription factors, and epigenetic modifications, as well as its lncRNA/pseudogene/circRNA-miRNA-mRNA regulatory network. The findings suggest that the normal function of NEIL1 is a significant factor in human health and longevity, with defects in NEIL1 potentially leading to various cancers and related syndromes, Alzheimer's disease, obesity, and diabetes.

Identifying a gene signature of metastatic potential by linking pre-metastatic state to ultimate metastatic fate

Identifying the key molecular pathways that enable metastasis by analyzing the eventual metastatic tumor is challenging because the state of the founder subclone likely changes following metastatic colonization. To address this challenge, we labeled primary mouse pancreatic ductal adenocarcinoma (PDAC) subclones with DNA barcodes to characterize their pre-metastatic state using ATAC-seq and RNA-seq and determine their relative in vivo metastatic potential prospectively. We identified a gene signature separating metastasis-high and metastasis-low subclones orthogonal to the normal-to-PDAC and classical-to-basal axes. The metastasis-high subclones feature activation of IL-1 pathway genes and high NF-κB and Zeb/Snail family activity and the metastasis-low subclones feature activation of neuroendocrine, motility, and Wnt pathway genes and high CDX2 and HOXA13 activity. In a functional screen, we validated novel mediators of PDAC metastasis in the IL-1 pathway, including the NF-κB targets Fos and Il23a, and beyond the IL-1 pathway including Myo1b and Tmem40. We scored human PDAC tumors for our signature of metastatic potential from mouse and found that metastases have higher scores than primary tumors. Moreover, primary tumors with higher scores are associated with worse prognosis. We also found that our metastatic potential signature is enriched in other human carcinomas, suggesting that it is conserved across epithelial malignancies. This work establishes a strategy for linking cancer cell state to future behavior, reveals novel functional regulators of PDAC metastasis, and establishes a method for scoring human carcinomas based on metastatic potential.

C2CDB: an advanced platform integrating comprehensive information and analysis tools of cancer-related circRNAs

Motivation

Circular RNAs (circRNAs) play important roles in gene expression and their involvement in tumorigenesis is emerging. circRNA-related database is a powerful tool for researchers to investigate circRNAs. However, existing databases lack advanced platform integrating comprehensive information and analysis tools of cancer-related circRNAs.

Results

We developed a comprehensive platform called CircRNA to Cancer Database (C2CDB), encompassing 318 158 cancer-related circRNAs expressed in tumors and adjacent tissues across 30 types of cancers. C2CDB provides basic details such as sequence and expression levels of circRNAs, as well as crucial insights into biological mechanisms, including miRNA binding, RNA-binding protein interaction, coding potential, base modification, mutation, and secondary structure. Moreover, C2CDB collects an extensive compilation of published literature on cancer circRNAs, extracting and presenting pivotal content encompassing biological functions, underlying mechanisms, and molecular tools in these studies. Additionally, C2CDB offers integrated tools to analyse three potential mechanisms: circRNA-miRNA ceRNA interaction, circRNA encoding, and circRNA biogenesis, facilitating investigators with convenient access to highly reliable information. To enhance clarity and organization, C2CDB has meticulously curated and integrated the previously chaotic nomenclature of circRNAs, addressing the prevailing confusion and ambiguity surrounding their designations.

Availability and implementation

C2CDB is freely available at http://pengyonglab.com/c2cdb.

Identification of CNTN2 as a genetic modifier of PIGA-CDG through pedigree analysis of a family with incomplete penetrance and functional testing in Drosophila

Loss of function mutations in the X-linked PIGA gene lead to PIGA-CDG, an ultra-rare congenital disorder of glycosylation (CDG), typically presenting with seizures, hypotonia, and neurodevelopmental delay. We identified two brothers (probands) with PIGA-CDG, presenting with epilepsy and mild developmental delay. Both probands carry PIGA S132C , an ultra-rare variant predicted to be damaging. Strikingly, the maternal grandfather and a great-uncle also carry PIGA S132C , but neither presents with symptoms associated with PIGA-CDG. We hypothesized genetic modifiers may contribute to this reduced penetrance. Using whole genome sequencing and pedigree analysis, we identified possible susceptibility variants found in the probands and not in carriers and possible protective variants found in the carriers and not in the probands. Candidate variants included heterozygous, damaging variants in three genes also involved directly in GPI-anchor biosynthesis and a few genes involved in other glycosylation pathways or encoding GPI-anchored proteins. We functionally tested the predicted modifiers using a Drosophila eye-based model of PIGA-CDG. We found that loss of CNTN2 , a predicted protective modifier, rescues loss of PIGA in Drosophila eye-based model, like what we predict in the family. Further testing found that loss of CNTN2 also rescues patient-relevant phenotypes, including seizures and climbing defects in Drosophila neurological models of PIGA-CDG. By using pedigree information, genome sequencing, and in vivo testing, we identified CNTN2 as a strong candidate modifier that could explain the incomplete penetrance in this family. Identifying and studying rare disease modifier genes in human pedigrees may lead to pathways and targets that may be developed into therapies.

TIMP1 promotes thyroid cancer cell progression through macrophage phenotypic polarization via the PI3K/AKT signaling pathway

Increasing evidence suggests that tissue inhibitor of metalloproteinase 1 (TIMP1) played a pivotal role in immune regulation. Our study focused on examining the expression and function of TIMP1 in humans, particularly in its regulation of tumor-associated macrophages (TAMs) in papillary thyroid carcinoma (PTC). We observed an upregulation of TIMP1 in 16 different types of malignancies, including thyroid cancer. TIMP1 shaped the inflammatory TME in PTC. Inhibiting the expression of TIMP1 has been demonstrated to reduce the malignant biological traits of PTC cells. Furthermore, reducing TIMP1 expression impeded M2 macrophage polarization as well as facilitated M1 macrophage polarization in PTC. ELISA results demonstrated that downregulated TIMP1 expression correlated with decreased levels of IL10 and TGF-β in cell supernatants. Furthermore, the supernatant from polarized macrophages in the TIMP1-silenced group inhibited the motility of wild-type PTC cells. Therefore, TIMP1 may enhance the progression of PTC by stimulating the PI3K/AKT pathway via the secretion of IL10 and TGF-β, consequently influencing M2-type polarization in TAMs.

Quantitative comparison of single-cell RNA sequencing versus single-molecule RNA imaging for quantifying transcriptional noise

Stochastic fluctuations (noise) in transcription generate substantial cell-to-cell variability. However, how best to quantify genome-wide noise, remains unclear. Here we utilize a small-molecule perturbation (IdU) to amplify noise and assess noise quantification from numerous scRNA-seq algorithms on human and mouse datasets, and then compare to noise quantification from single-molecule RNA FISH (smFISH) for a panel of representative genes. We find that various scRNA-seq analyses report amplified noise, without altered mean-expression levels, for ~90% of genes and that smFISH analysis verifies noise amplification for the vast majority of genes tested. Collectively, the analyses suggest that most scRNA-seq algorithms are appropriate for quantifying noise including a simple normalization approach, although all of these systematically underestimate noise compared to smFISH. From a practical standpoint, this analysis argues that IdU is a globally penetrant noise-enhancer molecule-amplifying noise without altering mean-expression levels-which could enable investigations of the physiological impacts of transcriptional noise.

Upstream alternative polyadenylation in SCN5A produces a short transcript isoform encoding a mitochondria-localized NaV1.5 N-terminal fragment that influences cardiomyocyte respiration

SCN5A encodes the cardiac voltage-gated Na+ channel, NaV1.5, that initiates action potentials. SCN5A gene variants cause arrhythmias and increased heart failure risk. Mechanisms controlling NaV1.5 expression and activity are not fully understood. We recently found a well-conserved alternative polyadenylation (APA) signal downstream of the first SCN5A coding exon. This yields a SCN5A-short transcript isoform expressed in several species (e.g. human, pig, and cat), though rodents lack this upstream APA. Reanalysis of transcriptome-wide cardiac APA-seq and mRNA-seq data shows reductions in both upstream APA usage and short/full-length SCN5A mRNA ratios in failing hearts. Knock-in of the human SCN5A APA sequence into mice is sufficient to enable expression of SCN5A -short transcript, while significantly decreasing expression of full-length SCN5A mRNA. Notably, SCN5A -short transcript encodes a novel protein (NaV1.5-NT), composed of an N-terminus identical to NaV1.5 and a unique C-terminus derived from intronic sequence. AAV9 constructs were able to achieve stable NaV1.5-NT expression in mouse hearts, and western blot of human heart tissues showed bands co-migrating with NaV1.5-NT transgene-derived bands. NaV1.5-NT is predicted to contain a mitochondrial targeting sequence and localizes to mitochondria in cultured cardiomyocytes and in mouse hearts. NaV1.5-NT expression in cardiomyocytes led to elevations in basal oxygen consumption rate, ATP production, and mitochondrial ROS, while depleting NADH supply. Native PAGE analyses of mitochondria lysates revealed that NaV1.5-NT expression resulted in increased levels of disassembled complex V subunits and accumulation of complex I-containing supercomplexes. Overall, we discovered that APA-mediated regulation of SCN5A produces a short transcript encoding NaV1.5-NT. Our data support that NaV1.5-NT plays a multifaceted role in influencing mitochondrial physiology: 1) by increasing basal respiration likely through promoting complex V conformations that enhance proton leak, and 2) by increasing overall respiratory efficiency and NADH consumption by enhancing formation and/or stability of complex I-containing respiratory supercomplexes, though the specific molecular mechanisms underlying each of these remain unresolved.

Genes Co-Expressed with ESR2 Influence Clinical Outcomes in Cancer Patients: TCGA Data Analysis

ERβ has been assigned a tumor suppressor role in many cancer types. However, as conflicting findings emerge, ERβ's tissue-specific expression and functional role have remained elusive. There remains a notable gap in compact and comprehensive analyses of ESR2 mRNA expression levels across diverse tumor types coupled with an exploration of its potential gene network. In this study, we aim to address these gaps by presenting a comprehensive analysis of ESR2 transcriptomic data. We distinguished cancer types with significant changes in ESR2 expression levels compared to corresponding healthy tissue and concluded that ESR2 influences patient survival. Gene Set Enrichment Analysis (GSEA) distinguished molecular pathways affected by ESR2, including oxidative phosphorylation and epithelial-mesenchymal transition. Finally, we investigated genes displaying similar expression patterns as ESR2 in tumor tissues, identifying potential co-expressed genes that may exert a synergistic effect on clinical outcomes, with significant results, including the expression of ACIN1, SYNE2, TNFRSF13C, and MDM4. Collectively, our results highlight the significant influence of ESR2 mRNA expression on the transcriptomic landscape and the overall metabolism of cancerous cells across various tumor types.

Pathogenic variants in KMT2C result in a neurodevelopmental disorder distinct from Kleefstra and Kabuki syndromes

Trithorax-related H3K4 methyltransferases, KMT2C and KMT2D, are critical epigenetic modifiers. Haploinsufficiency of KMT2C was only recently recognized as a cause of neurodevelopmental disorder (NDD), so the clinical and molecular spectrums of the KMT2C-related NDD (now designated as Kleefstra syndrome 2) are largely unknown. We ascertained 98 individuals with rare KMT2C variants, including 75 with protein-truncating variants (PTVs). Notably, ∼15% of KMT2C PTVs were inherited. Although the most highly expressed KMT2C transcript consists of only the last four exons, pathogenic PTVs were found in almost all the exons of this large gene. KMT2C variant interpretation can be challenging due to segmental duplications and clonal hematopoesis-induced artifacts. Using samples from 27 affected individuals, divided into discovery and validation cohorts, we generated a moderate strength disorder-specific KMT2C DNA methylation (DNAm) signature and demonstrate its utility in classifying non-truncating variants. Based on 81 individuals with pathogenic/likely pathogenic variants, we demonstrate that the KMT2C-related NDD is characterized by developmental delay, intellectual disability, behavioral and psychiatric problems, hypotonia, seizures, short stature, and other comorbidities. The facial module of PhenoScore, applied to photographs of 34 affected individuals, reveals that the KMT2C-related facial gestalt is significantly different from the general NDD population. Finally, using PhenoScore and DNAm signatures, we demonstrate that the KMT2C-related NDD is clinically and epigenetically distinct from Kleefstra and Kabuki syndromes. Overall, we define the clinical features, molecular spectrum, and DNAm signature of the KMT2C-related NDD and demonstrate they are distinct from Kleefstra and Kabuki syndromes highlighting the need to rename this condition.

Sleep and diurnal alternative polyadenylation sites associated with human APA-linked brain disorders

Disruption of sleep and circadian rhythms are a comorbid feature of many pathologies, and can negatively influence many health conditions, including neurodegenerative disease, metabolic illness, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. The interaction between sleep and/or circadian rhythms with the use of Alternative Polyadenylation (APA) has been largely undescribed, particularly in the context of other disorders. APA is a process that generates various transcript isoforms of the same gene affecting its mRNA translation, stability, localization, and subsequent function. Here we identified unique APAs expressed in rat brain over time-of-day, immediately following sleep deprivation, and the subsequent recovery period. From these data, we performed a secondary analysis of these sleep- or time-of-day associated PASs with recently described APA-linked human brain disorder susceptibility genes.

βA3/A1-crystallin is an epigenetic regulator of histone deacetylase 3 (HDAC3) in the retinal pigmented epithelial (RPE) cells

The retinal pigmented epithelial (RPE) cells maintain retinal homeostasis, and alterations in their function contribute to non-exudative age-related macular degeneration (AMD) 1,2 . Here, we explore the intricate relationship between RPE cells, epigenetic modifications, and the development of AMD. Importantly, the study reveals a substantial decrease in histone deacetylase 3 (HDAC3) activity and elevated histone acetylation in the RPE of human AMD donor eyes. To investigate epigenetic mechanisms in AMD development, we used a mouse model with RPE-specific Cryba1 knockout 3-5 , revealing that the loss of βA3/A1-crystallin selectively reduces HDAC3 activity, resulting in increased histone acetylation. βA3/A1-crystallin activates HDAC3 by facilitating its interaction with the casein kinase II (CK2) and phosphorylating HDAC3, as well as by regulating intracellular InsP6 (phytic acid) levels, required for activating HDAC3. These findings highlight a novel function of βA3/A1-crystallin as an epigenetic regulator of HDAC3 in the RPE cells and provide insights into potential therapeutic strategies in non-exudative AMD.

FOSL1-mediated LINC01566 negatively regulates CD4+ T-cell activation in myasthenia gravis

BACKGROUND

Myasthenia gravis (MG) is an autoimmune disease characterized by pathogenic antibodies that target structures of the neuromuscular junction. The evidence suggests that the regulation of long noncoding RNAs (lncRNAs) that is mediated by transcription factors (TFs) plays a key role in the pathophysiology of MG. Nevertheless, the detailed molecular mechanisms of lncRNAs in MG remain largely undetermined.

METHODS

Using microarray analysis, we analyzed the lncRNA levels in MG. By bioinformatics analysis, LINC01566 was found to potentially play an important role in MG. First, qRT‒PCR was performed to verify the LINC1566 expressions in MG patients. Then, fluorescence in situ hybridization was conducted to determine the localization of LINC01566 in CD4 + T cells. Finally, the impact of LINC01566 knockdown or overexpression on CD4 + T-cell function was also analyzed using flow cytometry and CCK-8 assay. A dual-luciferase reporter assay was used to validate the binding of the TF FOSL1 to the LINC01566 promoter.

RESULTS

Based on the lncRNA microarray and differential expression analyses, we identified 563 differentially expressed (DE) lncRNAs, 450 DE mRNAs and 19 DE TFs in MG. We then constructed a lncRNA-TF-mRNA network. Through network analysis, we found that LINC01566 may play a crucial role in MG by regulating T-cell-related pathways. Further experiments indicated that LINC01566 is expressed at low levels in MG patients. Functionally, LINC01566 is primarily distributed in the nucleus and can facilitate CD4 + T-cell apoptosis and inhibit cell proliferation. Mechanistically, we hypothesized that LINC01566 may negatively regulate the expressions of DUSP3, CCR2, FADD, SIRPB1, LGALS3 and SIRPB1, which are involved in the T-cell activation pathway, to further influence the cellular proliferation and apoptosis in MG. Moreover, we found that the effect of LINC01566 on CD4 + T cells in MG was mediated by the TF FOSL1, and in vitro experiments indicated that FOSL1 can bind to the promoter region of LINC01566.

CONCLUSIONS

In summary, our research revealed the protective roles of LINC01566 in clinical samples and cellular experiments, illustrating the potential roles and mechanism by which FOSL1/LINC01566 negatively regulates CD4 + T-cell activation in MG.

Optimized rAAV8 targeting acinar KLF4 ameliorates fibrosis in chronic pancreatitis via exosomes-enriched let-7s suppressing pancreatic stellate cells activation

Chronic pancreatitis (CP) is marked by progressive fibrosis and the activation of pancreatic stellate cells (PSCs), accompanied by the destruction of pancreatic parenchyma, leading to the loss of acinar cells (ACs). Few research studies have explored the mechanism by which damaged ACs (DACs) contribute to PSCs activation and pancreatic fibrosis. Currently, there are no effective drugs for curing CP or limiting the progression of pancreatic fibrosis. In this research, co-culture with intact acinar cells (IACs) suppressed PSC activation, while co-culture with DACs did the opposite. Krüppel-like factor 4 (KLF4) was significantly upregulated in DACs and was established as the key molecule that switches ACs from PSCs-suppressor to PSCs-activator. We revealed the exosomes of IACs contributed to the anti-activated function of IACs-CS on PSCs. MiRNome profiling showed that let-7 family is significantly enriched in IAC-derived exosomes (>30% miRNome), which partially mediates IACs' suppressive impacts on PSCs. Furthermore, it has been observed that the enrichment of let-7 in exosomes was influenced by the expression level of KLF4. Mechanistic studies demonstrated that KLF4 in ACs upregulated Lin28A, thereby decreasing let-7 levels in AC-derived exosomes, and thus promoting PSCs activation. We utilized an adeno-associated virus specifically targeting KLF4 in ACs (shKLF4-pAAV) to suppress PSCs activation in CP, resulting in reduced pancreatic fibrosis. IAC-derived exosomes hold potential as potent weapons against PSCs activation via let-7s, while activated KLF4/Lin28A signaling in DACs diminished such functions. ShKLF4-pAAV holds promise as a novel therapeutic approach for CP.

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Optical genome mapping identifies a homozygous deletion in the non-coding region of the SCN9A gene in individuals from the same family with congenital insensitivity to pain

We report an index patient with complete insensitivity to pain and a history of painless fractures, joint hypermobility, and behavioral problems. The index patient descends from a family with notable cases among his maternal relatives, including his aunt and his mother's first cousin, both of whom suffer from congenital insensitivity to pain. The patient had normal results for prior genetic testing: fragile-X syndrome testing, chromosomal microarray analysis, and exome sequencing. Optical genome mapping detected a homozygous deletion affecting the noncoding 5' untranslated region (UTR) and the first non-coding exon of the SCN9A gene in all affected family members, compatible with recessive disease transmission. Pathogenic homozygous loss-of-function variants in the SCN9A gene are associated with impaired pain sensation in humans. Optical genome mapping can thus detect pathogenic structural variants in patients without molecular etiology by standard diagnostic procedures and is a more accessible diagnostic tool than short-read or long-read whole-genome sequencing.

BertSNR: an interpretable deep learning framework for single-nucleotide resolution identification of transcription factor binding sites based on DNA language model

MOTIVATION

Transcription factors are pivotal in the regulation of gene expression, and accurate identification of transcription factor binding sites (TFBSs) at high resolution is crucial for understanding the mechanisms underlying gene regulation. The task of identifying TFBSs from DNA sequences is a significant challenge in the field of computational biology today. To address this challenge, a variety of computational approaches have been developed. However, these methods face limitations in their ability to achieve high-resolution identification and often lack interpretability.

RESULTS

We propose BertSNR, an interpretable deep learning framework for identifying TFBSs at single-nucleotide resolution. BertSNR integrates sequence-level and token-level information by multi-task learning based on pre-trained DNA language models. Benchmarking comparisons show that our BertSNR outperforms the existing state-of-the-art methods in TFBS predictions. Importantly, we enhanced the interpretability of the model through attentional weight visualization and motif analysis, and discovered the subtle relationship between attention weight and motif. Moreover, BertSNR effectively identifies TFBSs in promoter regions, facilitating the study of intricate gene regulation.

AVAILABILITY AND IMPLEMENTATION

The BertSNR source code can be found at https://github.com/lhy0322/BertSNR.

N6-methyladenosine modified circPAK2 promotes lymph node metastasis via targeting IGF2BPs/VEGFA signaling in gastric cancer

Circular RNAs (circRNAs) have emerged as key regulators of cancer occurrence and progression, as well as promising biomarkers for cancer diagnosis and prognosis. However, the potential mechanisms of circRNAs implicated in lymph node (LN) metastasis of gastric cancer remain unclear. Herein, we identify a novel N6-methyladenosine (m6A) modified circRNA, circPAK2, which is significantly upregulated in gastric cancer tissues and metastatic LN tissues. Functionally, circPAK2 enhances the migration, invasion, lymphangiogenesis, angiogenesis, epithelial-mesenchymal transition (EMT), and metastasis of gastric cancer in vitro and in vivo. Mechanistically, circPAK2 is exported by YTH domain-containing protein 1 (YTHDC1) from the nucleus to the cytoplasm in an m6A methylation-dependent manner. Moreover, increased cytoplasmic circPAK2 interacts with Insulin-Like Growth Factor 2 mRNA-Binding Proteins (IGF2BPs) and forms a circPAK2/IGF2BPs/VEGFA complex to stabilize VEGFA mRNA, which contributes to gastric cancer vasculature formation and aggressiveness. Clinically, high circPAK2 expression is positively associated with LN metastasis and poor prognosis in gastric cancer. This study highlights m6A-modified circPAK2 as a key regulator of LN metastasis of gastric cancer, thus supporting circPAK2 as a promising therapeutic target and prognostic biomarker for gastric cancer.

A map of canine sequence variation relative to a Greenland wolf outgroup

For over 15 years, canine genetics research relied on a reference assembly from a Boxer breed dog named Tasha (i.e., canFam3.1). Recent advances in long-read sequencing and genome assembly have led to the development of numerous high-quality assemblies from diverse canines. These assemblies represent notable improvements in completeness, contiguity, and the representation of gene promoters and gene models. Although genome graph and pan-genome approaches have promise, most genetic analyses in canines rely upon the mapping of Illumina sequencing reads to a single reference. The Dog10K consortium, and others, have generated deep catalogs of genetic variation through an alignment of Illumina sequencing reads to a reference genome obtained from a German Shepherd Dog named Mischka (i.e., canFam4, UU_Cfam_GSD_1.0). However, alignment to a breed-derived genome may introduce bias in genotype calling across samples. Since the use of an outgroup reference genome may remove this effect, we have reprocessed 1929 samples analyzed by the Dog10K consortium using a Greenland wolf (mCanLor1.2) as the reference. We efficiently performed remapping and variant calling using a GPU-implementation of common analysis tools. The resulting call set removes the variability in genetic differences seen across samples and breed relationships revealed by principal component analysis are not affected by the choice of reference genome. Using this sequence data, we inferred the history of population sizes and found that village dog populations experienced a 9-13 fold reduction in historic effective population size relative to wolves.

Complete sequencing of ape genomes

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.

High expression of YEATS2 as a predictive factor of poor prognosis in patients with hepatocellular carcinoma

YEATS domain containing 2 (YEATS2), it may function as a proto-oncogene. This study aims to investigate if YEATS2 correlates with prognosis in hepatocellular carcinoma. The prognostic landscape of YEATS2 and its relationship with expression in hepatocellular carcinoma were deciphered with public databases, RT-qPCR and western-blot in tissue samples. The expression profiling and prognostic value of YEATS2 were explored using UALCAN, TIMER, OncoLnc database. Transcription and survival analyses of YEATS2 in hepatocellular carcinoma were investigated with cBioPortal database. The STRING database was explored to identify molecular functions and signaling pathways downstream of YEATS2. YEATS2 expression was significantly higher in hepatocellular carcinoma compared with adjacent non-malignant tissues. Promoter methylation of YEATS2 exhibited different patterns in hepatocellular carcinoma. High expression of YEATS2 was associated with poorer survival. Mechanistically, YEATS2 was involved in mediating multiple biological processes including morphogenesis and migration.

Impact of genome build on RNA-seq interpretation and diagnostics

Transcriptomics is a powerful tool for unraveling the molecular effects of genetic variants and disease diagnosis. Prior studies have demonstrated that choice of genome build impacts variant interpretation and diagnostic yield for genomic analyses. To identify the extent genome build also impacts transcriptomics analyses, we studied the effect of the hg19, hg38, and CHM13 genome builds on expression quantification and outlier detection in 386 rare disease and familial control samples from both the Undiagnosed Diseases Network and Genomics Research to Elucidate the Genetics of Rare Disease Consortium. Across six routinely collected biospecimens, 61% of quantified genes were not influenced by genome build. However, we identified 1,492 genes with build-dependent quantification, 3,377 genes with build-exclusive expression, and 9,077 genes with annotation-specific expression across six routinely collected biospecimens, including 566 clinically relevant and 512 known OMIM genes. Further, we demonstrate that between builds for a given gene, a larger difference in quantification is well correlated with a larger change in expression outlier calling. Combined, we provide a database of genes impacted by build choice and recommend that transcriptomics-guided analyses and diagnoses are cross referenced with these data for robustness.

Epigenetic Alteration of Smooth Muscle Cells Regulates Endothelin-Dependent Blood Pressure and Hypertensive Arterial Remodeling

Long-standing hypertension (HTN) affects multiple organ systems and leads to pathologic arterial remodeling, which is driven largely by smooth muscle cell (SMC) plasticity. Although genome wide association studies (GWAS) have identified numerous variants associated with changes in blood pressure in humans, only a small percentage of these variants actually cause HTN. In order to identify relevant genes important in SMC function in HTN, we screened three separate human GWAS and Mendelian randomization studies to identify SNPs located within non-coding gene regions, focusing on genes encoding epigenetic enzymes, as these have been recently identified to control SMC fate in cardiovascular disease. We identified SNPs rs62059712 and rs74480102 in the promoter of the human JMJD3 gene and show that the minor C allele increases JMJD3 transcription in SMCs via increased SP1 binding to the JMJD3 promoter. Using our novel SMC-specific Jmjd3-deficient murine model ( Jmjd3 flox/flox Myh11 CreERT ), we show that loss of Jmjd3 in SMCs results in HTN, mechanistically, due to decreased EDNRB expression and a compensatory increase in EDNRA expression. As a translational corollary, through single cell RNA-sequencing (scRNA-seq) of human arteries, we found strong correlation between JMJD3 and EDNRB expression in SMCs. Further, we identified that JMJD3 is required for SMC-specific gene expression, and loss of JMJD3 in SMCs in the setting of HTN results in increased arterial remodeling by promoting the SMC synthetic phenotype. Our findings link a HTN-associated human DNA variant with regulation of SMC plasticity, revealing therapeutic targets that may be used in the screening and/or personalized treatment of HTN.

IsoVis - a webserver for visualization and annotation of alternative RNA isoforms

Genes commonly express multiple RNA products (RNA isoforms), which differ in exonic content and can have different functions. Making sense of the plethora of known and novel RNA isoforms being identified by transcriptomic approaches requires a user-friendly way to visualize gene isoforms and how they differ in exonic content, expression levels and potential functions. Here we introduce IsoVis, a freely available webserver that accepts user-supplied transcriptomic data and visualizes the expressed isoforms in a clear, intuitive manner. IsoVis contains numerous features, including the ability to visualize all RNA isoforms of a gene and their expression levels; the annotation of known isoforms from external databases; mapping of protein domains and features to exons, allowing changes to protein sequence and function between isoforms to be established; and extensive species compatibility. Datasets visualised on IsoVis remain private to the user, allowing analysis of sensitive data. IsoVis visualisations can be downloaded to create publication-ready figures. The IsoVis webserver enables researchers to perform isoform analyses without requiring programming skills, is free to use, and available at https://isomix.org/isovis/.

MotifbreakR v2: extended capability and database integration

MotifbreakR is a software tool that scans genetic variants against position weight matrices of transcription factors (TF) to determine the potential for the disruption of TF binding at the site of the variant. It leverages the Bioconductor suite of software packages and annotations to operate across a diverse array of genomes and motif databases. Initially developed to interrogate the effect of single nucleotide variants (common and rare SNVs) on potential TF binding sites, in motifbreakR v2, we have updated the functionality. New features include the ability to query other types of more complex genetic variants, such as short insertions and deletions (indels). This function allows modeling a more extensive array of variants that may have more significant effects on TF binding. Additionally, while TF binding is based partly on sequence preference, predictions of TF binding based on sequence preference alone can indicate many more potential binding events than observed. Adding information from DNA-binding sequencing datasets lends confidence to motif disruption prediction by demonstrating TF binding in cell lines and tissue types. Therefore, motifbreakR implements querying the ReMap2022 database for evidence that a TF matching the disrupted motif binds over the disrupting variant. Finally, in motifbreakR, in addition to the existing interface, we have implemented an R/Shiny graphical user interface to simplify and enhance access to researchers with different skill sets.

Benchmarking Algorithms for Gene Set Scoring of Single-cell ATAC-seq Data

Gene set scoring (GSS) has been routinely conducted for gene expression analysis of bulk or single-cell RNA sequencing (RNA-seq) data, which helps to decipher single-cell heterogeneity and cell type-specific variability by incorporating prior knowledge from functional gene sets. Single-cell assay for transposase accessible chromatin using sequencing (scATAC-seq) is a powerful technique for interrogating single-cell chromatin-based gene regulation, and genes or gene sets with dynamic regulatory potentials can be regarded as cell type-specific markers as if in single-cell RNA-seq (scRNA-seq). However, there are few GSS tools specifically designed for scATAC-seq, and the applicability and performance of RNA-seq GSS tools on scATAC-seq data remain to be investigated. Here, we systematically benchmarked ten GSS tools, including four bulk RNA-seq tools, five scRNA-seq tools, and one scATAC-seq method. First, using matched scATAC-seq and scRNA-seq datasets, we found that the performance of GSS tools on scATAC-seq data was comparable to that on scRNA-seq, suggesting their applicability to scATAC-seq. Then, the performance of different GSS tools was extensively evaluated using up to ten scATAC-seq datasets. Moreover, we evaluated the impact of gene activity conversion, dropout imputation, and gene set collections on the results of GSS. Results show that dropout imputation can significantly promote the performance of almost all GSS tools, while the impact of gene activity conversion methods or gene set collections on GSS performance is more dependent on GSS tools or datasets. Finally, we provided practical guidelines for choosing appropriate preprocessing methods and GSS tools in different application scenarios.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Identification of Polymorphisms in EAAT1 Glutamate Transporter Gene SLC1A3 Associated with Reduced Migraine Risk

Dysfunction in ion channels or processes involved in maintaining ionic homeostasis is thought to lower the threshold for cortical spreading depression (CSD), and plays a role in susceptibility to associated neurological disorders, including pathogenesis of a migraine. Rare pathogenic variants in specific ion channels have been implicated in monogenic migraine subtypes. In this study, we further examined the channelopathic nature of a migraine through the analysis of common genetic variants in three selected ion channel or transporter genes: SLC4A4, SLC1A3, and CHRNA4. Using the Agena MassARRAY platform, 28 single-nucleotide polymorphisms (SNPs) across the three candidate genes were genotyped in a case-control cohort comprised of 182 migraine cases and 179 matched controls. Initial results identified significant associations between migraine and rs3776578 (p = 0.04) and rs16903247 (p = 0.05) genotypes within the SLC1A3 gene, which encodes the EAAT1 glutamate transporter. These SNPs were subsequently genotyped in an independent cohort of 258 migraine cases and 290 controls using a high-resolution melt assay, and association testing supported the replication of initial findings-rs3776578 (p = 0.0041) and rs16903247 (p = 0.0127). The polymorphisms are in linkage disequilibrium and localise within a putative intronic enhancer region of SLC1A3. The minor alleles of both SNPs show a protective effect on migraine risk, which may be conferred via influencing the expression of SLC1A3.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

A Transcriptomic Dataset of Embryonic Murine Telencephalon

Sex bias is known in the prevalence/pathology of neurodevelopmental disorders. Sex-dependent differences of the certain brain areas are known to emerge perinatally through the exposure to sex hormones, while gene expression patterns in the rodent embryonic brain does not seem to be completely the same between male and female. To investigate potential sex differences in gene expression and cortical organization during the embryonic period in mice, we conducted a comprehensive analysis of gene expression for the telencephalon at embryonic day (E) 11.5 (a peak of neural stem cell expansion) and E14.5 (a peak of neurogenesis) using bulk RNA-seq data. As a result, our data showed the existence of notable sex differences in gene expression patterns not obviously at E11.5, but clearly at E14.5 when neurogenesis has become its peak. These data can be useful for exploring potential contribution of genes exhibiting sex differences to the divergence in brain development. Additionally, our data underscore the significance of studying the embryonic period to gain a deeper understanding of sex differences in brain development.

The Maintenance of Deleterious Variation in Wild Chinese Rhesus Macaques

Understanding how deleterious variation is shaped and maintained in natural populations is important in conservation and evolutionary biology, as decreased fitness caused by these deleterious mutations can potentially lead to an increase in extinction risk. It is known that demographic processes can influence these patterns. For example, population bottlenecks and inbreeding increase the probability of inheriting identical-by-descent haplotypes from a recent common ancestor, creating long tracts of homozygous genotypes called runs of homozygosity (ROH), which have been associated with an accumulation of mildly deleterious homozygotes. Counterintuitively, positive selection can also maintain deleterious variants in a population through genetic hitchhiking. Here, we analyze the whole genomes of 79 wild Chinese rhesus macaques across five subspecies and characterize patterns of deleterious variation with respect to ROH and signals of recent positive selection. We show that the fraction of homozygotes occurring in long ROH is significantly higher for deleterious homozygotes than tolerated ones, whereas this trend is not observed for short and medium ROH. This confirms that inbreeding, by generating these long tracts of homozygosity, is the main driver of the high burden of homozygous deleterious alleles in wild macaque populations. Furthermore, we show evidence that homozygous LOF variants are being purged. Next, we identify seven deleterious variants at high frequency in regions putatively under selection near genes involved with olfaction and other processes. Our results shed light on how evolutionary processes can shape the distribution of deleterious variation in wild nonhuman primates.

Severity of effect considerations regarding the use of mutation as a toxicological endpoint for risk assessment: A report from the 8th International Workshop on Genotoxicity Testing (IWGT)

Exposure levels without appreciable human health risk may be determined by dividing a point of departure on a dose-response curve (e.g., benchmark dose) by a composite adjustment factor (AF). An "effect severity" AF (ESAF) is employed in some regulatory contexts. An ESAF of 10 may be incorporated in the derivation of a health-based guidance value (HBGV) when a "severe" toxicological endpoint, such as teratogenicity, irreversible reproductive effects, neurotoxicity, or cancer was observed in the reference study. Although mutation data have been used historically for hazard identification, this endpoint is suitable for quantitative dose-response modeling and risk assessment. As part of the 8th International Workshops on Genotoxicity Testing, a sub-group of the Quantitative Analysis Work Group (WG) explored how the concept of effect severity could be applied to mutation. To approach this question, the WG reviewed the prevailing regulatory guidance on how an ESAF is incorporated into risk assessments, evaluated current knowledge of associations between germline or somatic mutation and severe disease risk, and mined available data on the fraction of human germline mutations expected to cause severe disease. Based on this review and given that mutations are irreversible and some cause severe human disease, in regulatory settings where an ESAF is used, a majority of the WG recommends applying an ESAF value between 2 and 10 when deriving a HBGV from mutation data. This recommendation may need to be revisited in the future if direct measurement of disease-causing mutations by error-corrected next generation sequencing clarifies selection of ESAF values.