The UCSC Genome Browser database: 2024 update

Sexually dimorphic renal expression of mouse Klotho is directed by a kidney-specific distal enhancer responsive to HNF1b

Transcription enhancers are genomic sequences regulating common and tissue-specific genes and their disruption can contribute to human disease development and progression. Klotho, a sexually dimorphic gene specifically expressed in kidney, is well-linked to kidney dysfunction and its deletion from the mouse genome leads to premature aging and death. However, the sexually dimorphic regulation of Klotho is not understood. Here, we characterize two candidate Klotho enhancers using H3K27ac epigenetic marks and transcription factor binding and investigate their functions, individually and combined, through CRISPR-Cas9 genome engineering. We discovered that only the distal (E1), but not the proximal (E2) candidate region constitutes a functional enhancer, with the double deletion not causing Klotho expression to further decrease. E1 activity is dependent on HNF1b transcription factor binding site within the enhancer. Further, E1 controls the sexual dimorphism of Klotho as evidenced by qPCR and RNA-seq. Despite the sharp reduction of Klotho mRNA, unlike germline Klotho knockouts, mutant mice present normal phenotype, including weight, lifespan, and serum biochemistry. Lastly, only males lacking E1 display more prominent acute, but not chronic kidney injury responses, indicating a remarkable range of potential adaptation to isolated Klotho loss, especially in female E1 knockouts, retaining renoprotection despite over 80% Klotho reduction.

Enhlink infers distal and context-specific enhancer-promoter linkages

Enhlink is a computational tool for scATAC-seq data analysis, facilitating precise interrogation of enhancer function at the single-cell level. It employs an ensemble approach incorporating technical and biological covariates to infer condition-specific regulatory DNA linkages. Enhlink can integrate multi-omic data for enhanced specificity, when available. Evaluation with simulated and real data, including multi-omic datasets from the mouse striatum and novel promoter capture Hi-C data, demonstrate that Enhlink outperfoms alternative methods. Coupled with eQTL analysis, it identified a putative super-enhancer in striatal neurons. Overall, Enhlink offers accuracy, power, and potential for revealing novel biological insights in gene regulation.

Unexpected nuclear hormone receptor and chromatin dynamics regulate estrous cycle dependent gene expression

Chromatin changes in response to estrogen and progesterone are well established in cultured cells, but how they control gene expression under physiological conditions is largely unknown. To address this question, we examined in vivo estrous cycle dynamics of mouse uterus hormone receptor occupancy, chromatin accessibility and chromatin structure by combining RNA-seq, ATAC-seq, HiC-seq and ChIP-seq. Two estrous cycle stages were chosen for these analyses, diestrus (highest estrogen) and estrus (highest progesterone). Unexpectedly, rather than alternating with each other, estrogen receptor alpha (ERα) and progesterone receptor (PGR) were co-bound during diestrus and lost during estrus. Motif analysis of open chromatin followed by hypoxia inducible factor 2A (HIF2A) ChIP-seq and conditional uterine deletion of this transcription factor revealed a novel role for HIF2A in regulating diestrus gene expression patterns that were independent of either ERα or PGR binding. Proteins in complex with ERα included PGR and cohesin, only during diestrus. Combined with HiC-seq analyses, we demonstrate that complex chromatin architecture changes including enhancer switching are coordinated with ERα and PGR co-binding during diestrus and non-hormone receptor transcription factors such as HIF2A during estrus to regulate most differential gene expression across the estrous cycle.

Haplotype analysis identifies functional elements in monoclonal gammopathy of unknown significance

Genome-wide association studies (GWASs) based on common single nucleotide polymorphisms (SNPs) have identified several loci associated with the risk of monoclonal gammopathy of unknown significance (MGUS), a precursor condition for multiple myeloma (MM). We hypothesized that analyzing haplotypes might be more useful than analyzing individual SNPs, as it could identify functional chromosomal units that collectively contribute to MGUS risk. To test this hypothesis, we used data from our previous GWAS on 992 MGUS cases and 2910 controls from three European populations. We identified 23 haplotypes that were associated with the risk of MGUS at the genome-wide significance level (p < 5 × 10-8) and showed consistent results among all three populations. In 10 genomic regions, strong promoter, enhancer and regulatory element-related histone marks and their connections to target genes as well as genome segmentation data supported the importance of these regions in MGUS susceptibility. Several associated haplotypes affected pathways important for MM cell survival such as ubiquitin-proteasome system (RNF186, OTUD3), PI3K/AKT/mTOR (HINT3), innate immunity (SEC14L1, ZBP1), cell death regulation (BID) and NOTCH signaling (RBPJ). These pathways are important current therapeutic targets for MM, which may highlight the advantage of the haplotype approach homing to functional units.

MethNet: a robust approach to identify regulatory hubs and their distal targets from cancer data

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover cis regulatory elements (CREs) in a 1 Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPR interference based single-cell RNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

IS-PRM-Based Peptide Targeting Informed by Long-Read Sequencing for Alternative Proteome Detection

Alternative splicing is a major contributor of transcriptomic complexity, but the extent to which transcript isoforms are translated into stable, functional protein isoforms is unclear. Furthermore, detection of relatively scarce isoform-specific peptides is challenging, with many protein isoforms remaining uncharted due to technical limitations. Recently, a family of advanced targeted MS strategies, termed internal standard parallel reaction monitoring (IS-PRM), have demonstrated multiplexed, sensitive detection of predefined peptides of interest. Such approaches have not yet been used to confirm existence of novel peptides. Here, we present a targeted proteogenomic approach that leverages sample-matched long-read RNA sequencing (lrRNA-seq) data to predict potential protein isoforms with prior transcript evidence. Predicted tryptic isoform-specific peptides, which are specific to individual gene product isoforms, serve as "triggers" and "targets" in the IS-PRM method, Tomahto. Using the model human stem cell line WTC11, LR RNaseq data were generated and used to inform the generation of synthetic standards for 192 isoform-specific peptides (114 isoforms from 55 genes). These synthetic "trigger" peptides were labeled with super heavy tandem mass tags (TMT) and spiked into TMT-labeled WTC11 tryptic digest, predicted to contain corresponding endogenous "target" peptides. Compared to DDA mode, Tomahto increased detectability of isoforms by 3.6-fold, resulting in the identification of five previously unannotated isoforms. Our method detected protein isoform expression for 43 out of 55 genes corresponding to 54 resolved isoforms. This lrRNA-seq-informed Tomahto targeted approach is a new modality for generating protein-level evidence of alternative isoforms─a critical first step in designing functional studies and eventually clinical assays.

Polymerization of ZBTB transcription factors regulates chromatin occupancy

BCL6, an oncogenic transcription factor (TF), forms polymers in the presence of a small-molecule molecular glue that stabilizes a complementary interface between homodimers of BCL6's broad-complex, tramtrack, and bric-à-brac (BTB) domain. The BTB domains of other proteins, including a large class of TFs, have similar architectures and symmetries, raising the possibility that additional BTB proteins self-assemble into higher-order structures. Here, we surveyed 189 human BTB proteins with a cellular fluorescent reporter assay and identified 18 ZBTB TFs that show evidence of polymerization. Through biochemical and cryoelectron microscopy (cryo-EM) studies, we demonstrate that these ZBTB TFs polymerize into filaments. We found that BTB-domain-mediated polymerization of ZBTB TFs enhances chromatin occupancy within regions containing homotypic clusters of TF binding sites, leading to repression of target genes. Our results reveal a role of higher-order structures in regulating ZBTB TFs and suggest an underappreciated role for TF polymerization in modulating gene expression.

ChIP-Atlas 3.0: a data-mining suite to explore chromosome architecture together with large-scale regulome data

ChIP-Atlas (https://chip-atlas.org/) presents a suite of data-mining tools for analyzing epigenomic landscapes, powered by the comprehensive integration of over 376 000 public ChIP-seq, ATAC-seq, DNase-seq and Bisulfite-seq experiments from six representative model organisms. To unravel the intricacies of chromatin architecture that mediates the regulome-initiated generation of transcriptional and phenotypic diversity within cells, we report ChIP-Atlas 3.0 that enhances clarity by incorporating additional tracks for genomic and epigenomic features within a newly consolidated 'annotation track' section. The tracks include chromosomal conformation (Hi-C and eQTL datasets), transcriptional regulatory elements (ChromHMM and FANTOM5 enhancers), and genomic variants associated with diseases and phenotypes (GWAS SNPs and ClinVar variants). These annotation tracks are easily accessible alongside other experimental tracks, facilitating better elucidation of chromatin architecture underlying the diversification of transcriptional and phenotypic traits. Furthermore, 'Diff Analysis,' a new online tool, compares the query epigenome data to identify differentially bound, accessible, and methylated regions using ChIP-seq, ATAC-seq and DNase-seq, and Bisulfite-seq datasets, respectively. The integration of annotation tracks and the Diff Analysis tool, coupled with continuous data expansion, renders ChIP-Atlas 3.0 a robust resource for mining the landscape of transcriptional regulatory mechanisms, thereby offering valuable perspectives, particularly for genetic disease research and drug discovery.

SpadaHC: a database to improve the classification of variants in hereditary cancer genes in the Spanish population

Accurate classification of genetic variants is crucial for clinical decision-making in hereditary cancer. In Spain, genetic diagnostic laboratories have traditionally approached this task independently due to the lack of a dedicated resource. Here we present SpadaHC, a web-based database for sharing variants in hereditary cancer genes in the Spanish population. SpadaHC is implemented using a three-tier architecture consisting of a relational database, a web tool and a bioinformatics pipeline. Contributing laboratories can share variant classifications and variants from individuals in Variant Calling Format (VCF) format. The platform supports open and restricted access, flexible dataset submissions, automatic pseudo-anonymization, VCF quality control, variant normalization and liftover between genome builds. Users can flexibly explore and search data, receive automatic discrepancy notifications and access SpadaHC population frequencies based on many criteria. In February 2024, SpadaHC included 18 laboratory members, storing 1.17 million variants from 4306 patients and 16 343 laboratory classifications. In the first analysis of the shared data, we identified 84 genetic variants with clinically relevant discrepancies in their classifications and addressed them through a three-phase resolution strategy. This work highlights the importance of data sharing to promote consistency in variant classifications among laboratories, so patients and family members can benefit from more accurate clinical management. Database URL: https://spadahc.ciberisciii.es/.

Variants in Vitamin D-related Genes and Prostate Cancer Risk in Black Men

BACKGROUND

The relationship between vitamin D and prostate cancer has primarily been characterized among White men. However, Black men have higher prostate cancer incidence and mortality rates, chronically low circulating vitamin D levels, and ancestry-specific genetic variants in vitamin D-related genes. Here, we examine six critical genes in the vitamin D pathway and prostate cancer risk in Black men.

METHODS

We assessed a total of 69 candidate variants in six genes ( GC, CYP27A1, CYP27B1, CYP24A1, VDR , and RXRA ) including functional variants previously associated with prostate cancer and circulating 25(OHD) in White men. Associations with prostate cancer risk were examined using genome-wide association study data for approximately 10,000 prostate cancer cases and 10,000 controls among Black men and over 85,000 cases and 91,000 controls among White men. A statistical significance threshold of 0.000724 was used to account for the 69 variants tested.

RESULTS

None of the variants examined were significantly associated with prostate cancer risk among Black men after multiple comparison adjustment. Four variants tested P<0.05 in Black men, including two in RXRA (rs41400444 OR=1.09, 95% CI: 1.01-1.17, P = 0.024 and rs10881574 OR = 0.93, 0.87-1.00, P = 0.046) and two in VDR (rs2853563 OR = 1.07, 1.01-1.13, P = 0.017 and rs1156882 OR = 1.06, 1.00-1.12, P = 0.045). Two variants in VDR were also positively associated with risk in White men (rs11568820 OR = 1.04, 1.02-1.06, P = 0.00024 and rs4516035 OR = 1.03, 1.01-1.04, P = 0.00055).

CONCLUSION

We observed suggestive non-significant associations between genetic variants in RXRA and VDR and prostate cancer risk in Black men. Future research exploring the relationship of vitamin D with cancer risk in Black men will need larger sample sizes to identify ancestry-specific variants relevant to risk in this population.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

Mutations causing premature termination codons discriminate and generate cellular and clinical variability in HHT

ABSTRACT

For monogenic diseases caused by pathogenic loss-of-function DNA variants, attention focuses on dysregulated gene-specific pathways, usually considering molecular subtypes together within causal genes. To better understand phenotypic variability in hereditary hemorrhagic telangiectasia (HHT), we subcategorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense-mediated decay. In 3 patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from patients with ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC genotypes, PTC-containing RNA transcripts persisted at low levels (8%-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni P < .05 in HHT+/PTC BOECs clustered significantly only to generic protein terms (isopeptide-bond/ubiquitin-like conjugation) and pulse-chase experiments detected subtle protein maturation differences but no evidence for PTC-truncated protein. BOECs displaying highest PTC persistence were discriminated in unsupervised hierarchical clustering of near-invariant housekeeper genes, with patterns compatible with higher cellular stress in BOECs with >11% PTC persistence. To test directionality, we used a HeLa reporter system to detect induction of activating transcription factor 4 (ATF4), which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modeled relevant ENG domains, with AlphaMissense suggesting that readthrough substitutions would be benign for ENG R93X and other less rare ENG nonsense variants but more damaging for Q436X. We conclude that PTCs should be distinguished from other loss-of-function variants, PTC transcript levels increase in stressed cells, and readthrough proteins and mechanisms provide promising research avenues.

Inhibition of miR-199b-5p reduces pathological alterations in osteoarthritis by potentially targeting Fzd6 and Gcnt2

Osteoarthritis (OA) is a degenerative disease with a high prevalence in the elderly population, but our understanding of its mechanisms remains incomplete. Analysis of serum exosomal small RNA sequencing data from clinical patients and gene expression data from OA patient serum and cartilage obtained from the GEO database revealed a common dysregulated miRNA, miR-199b-5p. In vitro cell experiments demonstrated that miR-199b-5p inhibits chondrocyte vitality and promotes extracellular matrix degradation. Conversely, inhibition of miR-199b-5p under inflammatory conditions exhibited protective effects against damage. Local viral injection of miR-199b-5p into mice induced a decrease in pain threshold and OA-like changes. In an OA model, inhibition of miR-199b-5p alleviated the pathological progression of OA. Furthermore, bioinformatics analysis and experimental validation identified Gcnt2 and Fzd6 as potential target genes of MiR-199b-5p. Thus, these results indicated that MiR-199b-5p/Gcnt2 and Fzd6 axis might be a novel therapeutic target for the treatment of OA.

Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Pathogenicity Prediction of Gene Fusion in Structural Variations: A Knowledge Graph-Infused Explainable Artificial Intelligence (XAI) Framework

When analyzing cancer sample genomes in clinical practice, many structural variants (SVs), other than single nucleotide variants (SNVs), have been identified. To identify driver variants, the leading candidates must be narrowed down. When fusion genes are involved, selection is particularly difficult, and highly accurate predictions from AI is important. Furthermore, we also wanted to determine how the prediction can make more reliable diagnoses. Here, we developed an explainable AI (XAI) suitable for SVs with gene fusions, based on the XAI technology we previously developed for the prediction of SNV pathogenicity. To cope with gene fusion variants, we added new data to the previous knowledge graph for SVs and we improved the algorithm. Its prediction accuracy was as high as that of existing tools. Moreover, our XAI could explain the reasons for these predictions. We used some variant examples to demonstrate that the reasons are plausible in terms of pathogenic basic mechanisms. These results can be seen as a hopeful step toward the future of genomic medicine, where efficient and correct decisions can be made with the support of AI.

Germline polymorphisms of the NOD2 pathway may predict the effectiveness of radioiodine in differentiated thyroid cancer treatment

PURPOSE

Differentiated thyroid cancer (DTC) presents a complex clinical challenge, especially in patients with distant metastases and resistance to standard treatments. This study aimed to investigate the influence of specific genes and their germline single nucleotide polymorphisms (SNPs) linked to both inflammatory processes and other neoplasms on the clinical and pathological characteristics of DTC, particularly their potential impact on radioiodine (RAI) treatment efficacy.

METHODS

This retrospective analysis involved a cohort of 646 patients diagnosed with DTC after thyroidectomy. Study covering 1998-2014, updated in 2023, included 567 women and 79 men (median age: 49; range: 7-83). SNP selection targeted functional significance, while mutational status was assessed by pyrosequencing for comprehensive characterization. Patient genetic profiles were assessed for associations with disease characteristics, RAI response, and cancer pathology.

RESULTS

Significant correlations emerged between certain SNPs and DTC features. Notably, the NOD2 c.802 T > C variant (rs2066842) was identified as a marker distinguishing between papillary thyroid cancer (PTC) and follicular thyroid cancer (FTC). Moreover, the c.802 T allele was associated with an enhanced response to RAI treatment, indicating a more substantial decrease in posttreatment stimulated thyroglobulin (sTg) concentrations. The NFKB1A allele c.126A (rs696) exhibited connections with lower FTC stages and a reduced probability of multifocality.

CONCLUSION

This study explored the molecular mechanisms of particular SNPs, highlighting the role of NOD2 in innate immunity and the stress response, and its potential impact on RAI efficacy. This research underscores the clinical promise of SNP analysis and contributes to personalized treatment strategies for DTC, emphasizing the relevance of genetic factors in cancer progression and treatment outcomes.

IS-PRM-based peptide targeting informed by long-read sequencing for alternative proteome detection

Alternative splicing is a major contributor of transcriptomic complexity, but the extent to which transcript isoforms are translated into stable, functional protein isoforms is unclear. Furthermore, detection of relatively scarce isoform-specific peptides is challenging, with many protein isoforms remaining uncharted due to technical limitations. Recently, a family of advanced targeted MS strategies, termed internal standard parallel reaction monitoring (IS-PRM), have demonstrated multiplexed, sensitive detection of pre-defined peptides of interest. Such approaches have not yet been used to confirm existence of novel peptides. Here, we present a targeted proteogenomic approach that leverages sample-matched long-read RNA sequencing (LR RNAseq) data to predict potential protein isoforms with prior transcript evidence. Predicted tryptic isoform-specific peptides, which are specific to individual gene product isoforms, serve as "triggers" and "targets" in the IS-PRM method, Tomahto. Using the model human stem cell line WTC11, LR RNAseq data were generated and used to inform the generation of synthetic standards for 192 isoform-specific peptides (114 isoforms from 55 genes). These synthetic "trigger" peptides were labeled with super heavy tandem mass tags (TMT) and spiked into TMT-labeled WTC11 tryptic digest, predicted to contain corresponding endogenous "target" peptides. Compared to DDA mode, Tomahto increased detectability of isoforms by 3.6-fold, resulting in the identification of five previously unannotated isoforms. Our method detected protein isoform expression for 43 out of 55 genes corresponding to 54 resolved isoforms. This LR RNA seq-informed Tomahto targeted approach, called LRP-IS-PRM, is a new modality for generating protein-level evidence of alternative isoforms - a critical first step in designing functional studies and eventually clinical assays.

AnnoView enables large-scale analysis, comparison, and visualization of microbial gene neighborhoods

The analysis and comparison of gene neighborhoods is a powerful approach for exploring microbial genome structure, function, and evolution. Although numerous tools exist for genome visualization and comparison, genome exploration across large genomic databases or user-generated datasets remains a challenge. Here, we introduce AnnoView, a web server designed for interactive exploration of gene neighborhoods across the bacterial and archaeal tree of life. Our server offers users the ability to identify, compare, and visualize gene neighborhoods of interest from 30 238 bacterial genomes and 1672 archaeal genomes, through integration with the comprehensive Genome Taxonomy Database and AnnoTree databases. Identified gene neighborhoods can be visualized using pre-computed functional annotations from different sources such as KEGG, Pfam and TIGRFAM, or clustered based on similarity. Alternatively, users can upload and explore their own custom genomic datasets in GBK, GFF or CSV format, or use AnnoView as a genome browser for relatively small genomes (e.g. viruses and plasmids). Ultimately, we anticipate that AnnoView will catalyze biological discovery by enabling user-friendly search, comparison, and visualization of genomic data. AnnoView is available at http://annoview.uwaterloo.ca.

Intron-lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envision and test a hypothesis for their formation using Saccharomyces cerevisiae, documenting full length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron-lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full length and processed circles. Post-splicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

The 2024 Nucleic Acids Research database issue and the online molecular biology database collection

The 2024 Nucleic Acids Research database issue contains 180 papers from across biology and neighbouring disciplines. There are 90 papers reporting on new databases and 83 updates from resources previously published in the Issue. Updates from databases most recently published elsewhere account for a further seven. Nucleic acid databases include the new NAKB for structural information and updates from Genbank, ENA, GEO, Tarbase and JASPAR. The Issue's Breakthrough Article concerns NMPFamsDB for novel prokaryotic protein families and the AlphaFold Protein Structure Database has an important update. Metabolism is covered by updates from Reactome, Wikipathways and Metabolights. Microbes are covered by RefSeq, UNITE, SPIRE and P10K; viruses by ViralZone and PhageScope. Medically-oriented databases include the familiar COSMIC, Drugbank and TTD. Genomics-related resources include Ensembl, UCSC Genome Browser and Monarch. New arrivals cover plant imaging (OPIA and PlantPAD) and crop plants (SoyMD, TCOD and CropGS-Hub). The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). Over the last year the NAR online Molecular Biology Database Collection has been updated, reviewing 1060 entries, adding 97 new resources and eliminating 388 discontinued URLs bringing the current total to 1959 databases. It is available at http://www.oxfordjournals.org/nar/database/c/.