Polyadenylation Signal Research Articles

The full transcription map of cottontail rabbit papillomavirus in tumor tissues.

Cottontail rabbit papillomavirus (CRPV), the first papillomavirus associated with tumor development, has been used as a powerful model to study papillomavirus pathogenesis for more than 90 years. However, lack of a comprehensive analysis of the CRPV transcriptome has impeded the understanding of CRPV biology and molecular pathogenesis. Here, we report the construction of a complete CRPV transcription map from Hershey CRPV-induced skin tumor tissues. By using RNA-seq in combination with long-reads PacBio Iso-seq, 5' and 3' RACE, primer-walking RT-PCR, Northern blot, and RNA in situ hybridization, we demonstrated that the CRPV genome transcribes its early and late RNA transcripts unidirectionally from at least five distinct major promoters (P) and polyadenylates its transcripts at two major polyadenylation (pA) sites. The viral early transcripts are primarily transcribed from three "early" promoters, P90, P156, and P907 and polyadenylated at nt 4368 by using an early polyadenylation signal (PAS) at nt 4351. Like other low-risk human papillomaviruses and animal papillomaviruses, CRPV E6 and E7 transcripts are transcribed from three separate early promoters. Transcripts from two "late" promoters, P7525, and P1225, utilize either an early PAS for E1^E4 or a late PAS at 7399 for L2 and L1 RNA polyadenylation at nt 7415 to express capsid L2 and L1 proteins respectively. By using the mapped four 5' splice sites and three 3' splice sites, CRPV RNA transcripts undergo extensive alternative splicing to produce more than 33 viral RNA isoforms for production of at least 12 viral proteins, some of which without codon optimization are expressible in rabbit RK13 and human HEK293T cells. The constructed full CRPV transcription map in this study for the first time will enhance our understanding of the structures and expressions of CRPV genes and their contribution to molecular pathogenesis and tumorigenesis.

Dynamic landscape of m6A modifications and related post-transcriptional events in muscle-invasive bladder cancer

BackgroundMuscle-invasive bladder carcinoma (MIBC) is a serious and more advanced stage of bladder carcinoma. N6-Methyladenosine (m6A) is a dynamic and reversible modifications that primarily affects RNA stability and alternative splicing. The dysregulation of m6A in MIBC can be potential target for clinical interventions, but there have been limited studies on m6A modifications in MIBC and their associations with post-transcriptional regulatory processes.MethodsPaired tumor and adjacent-normal tissues were obtained from three patients with MIBC following radical cystectomy. The additional paired tissues for validation were obtained from patients underwent transurethral resection. Utilizing Nanopore direct-RNA sequencing, we characterized the m6A RNA methylation landscape in MIBC, with a focus on identifying post-transcriptional events potentially affected by changes in m6A sites. This included an examination of differential transcript usage, polyadenylation signal sites, and variations in poly(A) tail length, providing insights into the broader impact of m6A alterations on RNA processing in MIBC.ResultsThe prognostic-related m6A genes and m6A-risk model constructed by machine learning enables the stratification of high and low-risk patients with precision. A novel m6A modification site in the 3’ untranslated region (3’UTR) of IGLL5 gene were identified, characterized by a lower m6A methylation ratio, elongated poly(A) tails, and a notable bias in transcript usage. Furthermore, we discovered two particular transcripts, VWA1-203 and CEBPB-201. VWA1-203 displayed diminished m6A methylation levels, a truncated 3’UTR, and an elongated poly(A) tail, whereas CEBPB-201 showed opposite trends, highlighting the complex interplay between m6A modifications and RNA processing. Source code was provided on GitHub (https://github.com/lelelililele/Nanopore-m6A-analysis).ConclusionsThe state-of-the-art Nanopore direct-RNA sequencing and machine learning techniques enables comprehensive identification of m6A modification and provided insights into the potential post-transcriptional regulation mechanisms on the development and progression in MIBC.

The expressions of nr0b1/2 genes in the Chinese giant salamander and their response to estradiol benzoate/17α-methyltestosterone treatment

Previously, the nr0b1 and nr0b2 genes have been identified in various animals and found to be closely associated with sex determination. In order to investigate the regulatory roles of nr0b1 and nr0b2 genes in the sex determination of Chinese giant salamander (CGS, Andrias davidianus), we obtained the full-length cDNA sequences of nr0b1 and nr0b2 through cloning, determined the characteristics of their tissue expressions and the responses of nr0b1 and nr0b2 to estradiol benzoate (EB) and 17α-methyltestosterone (MT) by using quantitative real-time PCR (qRT-PCR), and evaluated gonad development through histological sections. The full-lengths of cDNA of nr0b1/2 in C GS are separately 1278bp and 789bp. They have complete open reading frames (ORF) of 864bp and 774bp and encode 287 and 257 amino acids respectively. The 3' untranslated regions (UTR) are 276bp and 15bp, and the 5' UTR of nr0b1 is 138bp. In the 3' UTR of nr0b1, a putative polyadenylation signal (AATAAA) is identified. Tissue distribution analysis indicates that nr0b1 is mainly expressed in the pituitary, kidney, and testes, with significant differences between females (ZW-F) and males (ZZ-M) expression in the testes. On the other hand, nr0b2 is mainly expressed in the liver, followed by gonadal tissues. After immersion and feeding with EB and MT, EB's promotion on nr0b1 expression in CGS ovaries and testes is shown, by contrast MT only has a positive effect on nr0b2 in CGS testes. Histological sections show that after EB treatment, both genetic females (ZW-F) and sex-reversed females (ZZ-F) of CGS contain oogonia, indicating that EB can cause sex reversal in CGS; however, after MT treatment, no cellular morphological changes are observed either in genetic females (ZW-F) or in genetic males (ZW-M), indicating that this concentration cannot induce a shift in the sex ratio of CGS. These findings lay the foundation for further investigating the molecular mechanisms by which nr0b1/2 influence the differentiation and maintenance of CGS gonads.

Evolutionary dynamics of polyadenylation signals and their recognition strategies in protists.

The poly(A) signal, together with auxiliary elements, directs cleavage of a pre-mRNA and thus determines the 3' end of the mature transcript. In many species, including humans, the poly(A) signal is an AAUAAA hexamer, but we recently found that the deeply branching eukaryote Giardia lamblia uses a distinct hexamer (AGURAA) and lacks any known auxiliary elements. Our discovery prompted us to explore the evolutionary dynamics of poly(A) signals and auxiliary elements in the eukaryotic kingdom. We use direct RNA sequencing to determine poly(A) signals for four protists within the Metamonada clade (which also contains G. lamblia) and two outgroup protists. These experiments reveal that the AAUAAA hexamer serves as the poly(A) signal in at least four different eukaryotic clades, indicating that it is likely the ancestral signal, whereas the unusual Giardia version is derived. We find that the use and relative strengths of auxiliary elements are also plastic; in fact, within Metamonada, species like G. lamblia make use of a previously unrecognized auxiliary element where nucleotides flanking the poly(A) signal itself specify genuine cleavage sites. Thus, despite the fundamental nature of pre-mRNA cleavage for the expression of all protein-coding genes, the motifs controlling this process are dynamic on evolutionary timescales, providing motivation for future biochemical and structural studies as well as new therapeutic angles to target eukaryotic pathogens.

MGCN-PolyA: An Integrated Computational Framework for Predicting Poly(A) Signals with Multiscale-gated Convolutional Networks

Background: The accurate recognition of the polyadenylation signal (PAS) from DNA sequences is essential for understanding gene transcriptional regulation. A variety of machine learning-based computational methods have been developed to predict PAS in recent years; however, their performance and their generalization ability are unsatisfactory. It is highly desirable to design more preferable computational approaches for PAS prediction. Methods: In this work, we developed an integrated framework MGCN-PolyA for PAS prediction across four species, including Homo sapiens, Bos taurus, Mus musculus, and Drosophila melanogaster. MGCN-Poly(A) benefits from the diversity of feature engineering and the effectiveness of the model architecture. We combined features from different perspectives, such as word embedding, One-hot encoding, K-mer frequency, and Enhanced Nucleic Acid Composition (ENAC), which complement each other and provide rich and comprehensive information for model learning. In model architecture, MGCN-Poly(A) leverages a two-channel multi-scale gated convolutional network to effectively learn high-level feature representations at different scales, and then combines the statistical features to predict PAS using random forest algorithm. These designs not only speed up network training, but also improves the generalization ability Results: The benchmarking experiments on the independent test datasets demonstrate that MGCNPolyA outperforms other state-of-the-art algorithms in identifying PAS. MGCN-PolyA has the highest accuracy on all test datasets, and its excellent performance on cross-species validation also demonstrates the robustness of our model. Conclusion: Extracting features from different perspectives is important for PAS recognition, and the integration of DNNs and shallow machine learning algorithms can improve the model performance.

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.

Comprehensive analysis of single-nucleotide variants and alternative polyadenylation between inbred and outbred pigs

Inbreeding can lead to the accumulation of homozygous single nucleotide polymorphisms (SNPs) in the genome, which can significantly affect gene expression and phenotype. In this study, we examined the impact of homozygous SNPs resulting from inbreeding on alternative polyadenylation (APA) site selection and the underlying genetic mechanisms using inbred Luchuan pigs. Genome resequencing revealed that inbreeding results in a high accumulation of homozygous SNPs within the pig genome. 3′ mRNA-seq on leg muscle, submandibular lymph node, and liver tissues was performed to identify differences in APA events between inbred and outbred Luchuan pigs. We revealed different tissue-specific APA usage caused by inbreeding, which were associated with different biological processes. Furthermore, we explored the role of polyadenylation signal (PAS) SNPs in APA regulation under inbreeding and identified key genes such as PUM1, SCARF1, RIPOR2, C1D, and LRRK2 that are involved in biological processes regulation. This study provides resources and sheds light on the impact of genomic homozygosity on APA regulation, offering insights into genetic characteristics and biological processes associated with inbreeding.

SIZ1-mediated SUMOylation of CPSF100 promotes plant thermomorphogenesis by controlling alternative polyadenylation

Under warm temperatures, plants adjust their morphologies for environmental adaption via precise gene expression regulation. However, the function and regulation of alternative polyadenylation (APA), an important fine-tuning of gene expression, remains unknown in plant thermomorphogenesis. In this study, we found that SUMOylation, a critical post-translational modification, is induced by a long-term treatment at warm temperatures via a SUMO ligase SIZ1 in Arabidopsis. Disruption of SIZ1 altered the global usage of polyadenylation signals and affected the APA dynamic of thermomorphogenesis-related genes. CPSF100, a key subunit of the CPSF complex for polyadenylation regulation, is SUMOylated by SIZ1. Importantly, we demonstrated that SUMOylation is essential for the function of CPSF100 in genome-wide polyadenylation site choice during thermomorphogenesis. Further analyses revealed that the SUMO conjugation on CPSF100 attenuates its interaction with two isoforms of its partner CPSF30, increasing the nuclear accumulation of CPSF100 for polyadenylation regulation. In summary, our study uncovers a regulatory mechanism of APA via SIZ1-mediated SUMOylation in plant thermomorphogenesis.

Ere, a Family of Short Interspersed Elements in the Genomes of Odd-Toed Ungulates (Perissodactyla).

Short Interspersed Elements (SINEs) are eukaryotic retrotransposons transcribed by RNA polymerase III (pol III). Many mammalian SINEs (T+ SINEs) contain a polyadenylation signal (AATAAA), a pol III transcription terminator, and an A-rich tail in their 3'-end. The RNAs of such SINEs have the capacity for AAUAAA-dependent polyadenylation, which is unique to pol III-generated transcripts. The structure, evolution, and polyadenylation of the Ere SINE of ungulates (horses, rhinos, and tapirs) were investigated in this study. A bioinformatics analysis revealed the presence of up to ~4 × 105 Ere copies in representatives of all three families. These copies can be classified into two large subfamilies, EreA and EreB, the former distinguished by an additional 60 bp sequence. The 3'-end of numerous EreA and all EreB copies exhibit a 50 bp sequence designated as a terminal domain (TD). The Ere family can be further subdivided into subfamilies EreA_0TD, EreA_1TD, EreB_1TD, and EreB_2TD, depending on the presence and number of terminal domains (TDs). Only EreA_0TD copies can be assigned to T+ SINEs as they contain the AATAAA signal and the TCTTT transcription terminator. The analysis of young Ere copies identified by comparison with related perissodactyl genomes revealed that EreA_0TD and, to a much lesser extent, EreB_2TD have retained retrotranspositional activity in the recent evolution of equids and rhinoceroses. The targeted mutagenesis and transfection of HeLa cells were used to identify sequences in equine EreA_0TD that are critical for the polyadenylation of its pol III transcripts. In addition to AATAAA and the transcription terminator, two sites in the 3' half of EreA, termed the β and τ signals, were found to be essential for this process. The evolution of Ere, with a particular focus on the emergence of T+ SINEs, as well as the polyadenylation signals are discussed in comparison with other T+ SINEs.

Delineating yeast cleavage and polyadenylation signals using deep learning.

3'-end cleavage and polyadenylation is an essential process for eukaryotic mRNA maturation. In yeast species, the polyadenylation signals that recruit the processing machinery are degenerate and remain poorly characterized compared with the well-defined regulatory elements in mammals. Here we address this issue by developing deep learning models to deconvolute degenerate cis-regulatory elements and quantify their positional importance in mediating yeast poly(A) site formation, cleavage heterogeneity, and strength. In S. cerevisiae, cleavage heterogeneity is promoted by the depletion of U-rich elements around poly(A) sites as well as multiple occurrences of upstream UA-rich elements. Sites with high cleavage heterogeneity show overall lower strength. The site strength and tandem site distances modulate alternative polyadenylation (APA) under the diauxic stress. Finally, we develop a deep learning model to reveal the distinct motif configuration of S. pombe poly(A) sites, which show more precise cleavage than S. cerevisiae Altogether, our deep learning models provide unprecedented insights into poly(A) site formation of yeast species, and our results highlight divergent poly(A) signals across distantly related species.

Determinants of the interaction between the 5′-leader of HIV-1 genome and human lysyl-tRNA synthetase in reverse transcription primer release process

Reverse transcription of human immunodeficiency virus type 1 (HIV-1) initiates from the 3′ end of human tRNALys3. The primer tRNALys3 is selectively packaged into the virus in the form of a complex with human lysyl-tRNA synthetase (LysRS). To facilitate reverse transcription initiation, part of the 5′ leader (5′L) of HIV-1 genomic RNA (gRNA) evolves a tRNA anticodon-like element (TLE), which binds LysRS and releases tRNALys3 for primer annealing and reverse transcription initiation. Although TLE has been identified as a key element in 5′L responsible for LysRS binding, how the conformations and various hairpin structures of 5′L regulate 5′L-LysRS interaction is not fully understood. Here, these factors have been individually investigated using direct and competitive fluorescence anisotropy binding experiments. Our data showed that the conformation of 5′L significantly influences its binding affinity with LysRS. The 5′L conformation favoring gRNA dimerization and packaging exhibits much weaker binding affinity with LysRS compared to the alternative 5′L conformation that is not selected for packaging. Additionally, dimerization of 5′L impairs LysRS-5′L interaction. Furthermore, among various regions of 5′L, both the primer binding site/TLE domain and the stem-loop 3 are important for LysRS interaction, whereas the dimerization initiation site and the splicing donor plays a minor role. In contrast, the presence of the transacting responsive and the polyadenylation signal hairpins slightly inhibit LysRS binding. These findings reveal that the conformation and various regions of the 5′L of HIV-1 genome regulate its interaction with human LysRS and the reverse transcription primer release process.

Differential expression of NAT1 pharmacogene in hormone receptor positive vs. negative female breast tumors may affect drug treatment.

Studies have reported overexpression of NAT1 gene for xenobiotic metabolizing arylamine N -acetyltransferase type 1 in estrogen receptor positive breast tumors, and this association has been linked to patient chemoresistance and response to tamoxifen. We probed the expression of NAT1 , using quantitative reverse transcription PCR to screen clinically characterized breast cancer tissue cDNA arrays. Primers detecting all NAT1 alternative transcripts were used, and the protocol and results are reported according to consensus guidelines. The clinical information about 166 tumor samples screened is provided, including tumor stage, estrogen and progesterone receptor status and HER2 expression. NAT1 was found to be significantly ( P < 0.001) upregulated in hormone receptor positive vs. negative tumors. No correlation was apparent between NAT1 and tumor stage or HER2 expression. Our findings demonstrate a strong correlation between the expression of NAT1 and steroid hormone receptors in breast tumors, supporting its possible utility as a pharmacogenetic biomarker or drug target. Of the two polymorphic NAT genes, NAT1 is the one primarily expressed in breast tissue, and is subjected to regulation by two differential promoters and more than one polyadenylation signal. Hormonal factors may enhance NAT1 gene expression at the transcriptional or epigenetic level, and tamoxifen has additionally been shown to inhibit NAT1 enzymatic activity. The outcome of tamoxifen treatment is also more favorable in patients with NAT1 overexpressing tumors. The study adds to the growing body of evidence implicating NAT1 in breast cancer and its pharmacological treatment.

Premarital Counseling on the Alpha Thalassemia Allele HBA2:c.*94A&gt;G

The mutation HBA2:c.*94A&gt;G (AATAAA&gt;AATAAG; rs63751269) is a 3′-UTR (3 prime untranslated region) single-nucleotide substitution in the polyadenylation (PA) signal of HBA2 (αPA:A→G). This pathogenic (CADD score, 14.92) variant is sporadic in the Arabian Peninsula. It results in inefficient mRNA processing, transcription termination, and possibly using an alternate cryptic downstream polyadenylation signal. As a result, the allele αT (or αT-Saudi) poses challenges in premarital counseling with respect to fetal risk of hemoglobin H disease. Homozygous HBA2:c.*94A&gt;G (αTα/αTα) results in moderate-to-severe microcytosis (mean red cell volume, MCV, 55 to 65 fL), reflecting markedly impaired hemoglobin synthesis (hemoglobin H disease). Homozygous rightward -α3.7 (a 3804-neocleotide deletion allele, NM_000517.4:c.[-2_-3delAC; −α3.7]), on the other hand, results in mild microcytosis (MCV, 70 to 75 fL, alpha-thalassemia trait). Thus, HBA2:c.*94A&gt;G is more damaging than -α3.7. Consistently, the value of MCV in compound heterozygosity, HBA2:c.*94A&gt;G and −α3.7, is 65 to 70 fL. We report here a healthy couple who presented for premarital counseling on their hemoglobinopathy. The man has homozygous HBA2:c.*94A&gt;G (αTα/αTα), and the woman has compound heterozygous (−α3.7/αTα, also annotated as: −3.7α/αTα). As a result, the genotype of their offspring would be that of the father (αTα/αTα) or the mother (−α3.7/αTα). The counseling was mainly based on the benign phenotypes of the parents. As both were asymptomatic and their anemia was clinically insignificant, they proceeded with the marriage.

Investigating the Genetic Diversity of Hepatitis Delta Virus in Hepatocellular Carcinoma (HCC): Impact on Viral Evolution and Oncogenesis in HCC.

Hepatitis delta virus (HDV), an RNA virus with two forms of the delta antigen (HDAg), relies on hepatitis B virus (HBV) for envelope proteins essential for hepatocyte entry. Hepatocellular carcinoma (HCC) ranks third in global cancer deaths, yet HDV's involvement remains uncertain. Among 300 HBV-associated HCC serum samples from Taiwan's National Health Research Institutes, 2.7% (8/300) tested anti-HDV positive, with 62.7% (5/8) of these also HDV RNA positive. Genotyping revealed HDV-2 in one sample, HDV-4 in two, and two samples showed mixed HDV-2/HDV-4 infection with RNA recombination. A mixed-genotype infection revealed novel mutations at the polyadenylation signal, coinciding with the ochre termination codon for the L-HDAg. To delve deeper into the possible oncogenic properties of HDV-2, the predominant genotype in Taiwan, which was previously thought to be less associated with severe disease outcomes, an HDV-2 cDNA clone was isolated from HCC for study. It demonstrated a replication level reaching up to 74% of that observed for a widely used HDV-1 strain in transfected cultured cells. Surprisingly, both forms of HDV-2 HDAg promoted cell migration and invasion, affecting the rearrangement of actin cytoskeleton and the expression of epithelial-mesenchymal transition markers. In summary, this study underscores the prevalence of HDV-2, HDV-4, and their mixed infections in HCC, highlighting the genetic diversity in HCC as well as the potential role of both forms of the HDAg in HCC oncogenesis.

Epitranscriptome insights into Riccia fluitans L. (Marchantiophyta) aquatic transition using nanopore direct RNA sequencing

BackgroundRiccia fluitans, an amphibious liverwort, exhibits a fascinating adaptation mechanism to transition between terrestrial and aquatic environments. Utilizing nanopore direct RNA sequencing, we try to capture the complex epitranscriptomic changes undergone in response to land-water transition.ResultsA significant finding is the identification of 45 differentially expressed genes (DEGs), with a split of 33 downregulated in terrestrial forms and 12 upregulated in aquatic forms, indicating a robust transcriptional response to environmental changes. Analysis of N6-methyladenosine (m6A) modifications revealed 173 m6A sites in aquatic and only 27 sites in the terrestrial forms, indicating a significant increase in methylation in the former, which could facilitate rapid adaptation to changing environments. The aquatic form showed a global elongation bias in poly(A) tails, which is associated with increased mRNA stability and efficient translation, enhancing the plant’s resilience to water stress. Significant differences in polyadenylation signals were observed between the two forms, with nine transcripts showing notable changes in tail length, suggesting an adaptive mechanism to modulate mRNA stability and translational efficiency in response to environmental conditions. This differential methylation and polyadenylation underline a sophisticated layer of post-transcriptional regulation, enabling Riccia fluitans to fine-tune gene expression in response to its living conditions.ConclusionsThese insights into transcriptome dynamics offer a deeper understanding of plant adaptation strategies at the molecular level, contributing to the broader knowledge of plant biology and evolution. These findings underscore the sophisticated post-transcriptional regulatory strategies Riccia fluitans employs to navigate the challenges of aquatic versus terrestrial living, highlighting the plant’s dynamic adaptation to environmental stresses and its utility as a model for studying adaptation mechanisms in amphibious plants.

Exchange of subtelomeric regions between chromosomes 4q and 10q reverts the FSHD genotype and phenotype.

The most common form of facioscapulohumeral dystrophy (FSHD1) is caused by a partial loss of the D4Z4 macrosatellite repeat array in the subtelomeric region of chromosome 4. Patients with FSHD1 typically carry 1 to 10 D4Z4 repeats, whereas nonaffected individuals have 11 to 150 repeats. The ~150-kilobyte subtelomeric region of the chromosome 10q exhibits a ~99% sequence identity to the 4q, including the D4Z4 array. Nevertheless, contractions of the chr10 array do not cause FSHD or any known disease, as in most people D4Z4 array on chr10 is flanked by the nonfunctional polyadenylation signal, not permitting the DUX4 expression. Here, we attempted to correct the FSHD genotype by a CRISPR-Cas9-induced exchange of the chr4 and chr10 subtelomeric regions. We demonstrated that the induced t(4;10) translocation can generate recombinant genotypes translated into improved FSHD phenotype. FSHD myoblasts with the t(4;10) exhibited reduced expression of the DUX4 targets, restored PAX7 target expression, reduced sensitivity to oxidative stress, and improved differentiation capacity.

Helicases DDX5 and DDX17 promote heterogeneity in HBV transcription termination in infected human hepatocytes

Transcription termination fine-tunes gene expression and contributes to the specification of RNA function in eukaryotic cells. Transcription termination of HBV is subject to the recognition of the canonical polyadenylation signal (cPAS) common to all viral transcripts. However, the regulation of this cPAS and its impact on viral gene expression and replication is currently unknown. To unravel the regulation of HBV transcript termination, we implemented a 3' RACE (rapid amplification of cDNA ends)-PCR assay coupled to single molecule sequencing both in invitro-infected hepatocytes and in chronically infected patients. The detection of a previously unidentified transcriptional readthrough indicated that the cPAS was not systematically recognized during HBV replication invitro and invivo. Gene expression downregulation experiments demonstrated a role for the RNA helicases DDX5 and DDX17 in promoting viral transcriptional readthrough, which was, in turn, associated with HBV RNA destabilization and decreased HBx protein expression. RNA and chromatin immunoprecipitation, together with mutation of the cPAS sequence, suggested a direct role of DDX5 and DDX17 in functionally linking cPAS recognition to transcriptional readthrough, HBV RNA stability and replication. Our findings identify DDX5 and DDX17 as crucial determinants of HBV transcriptional fidelity and as host restriction factors for HBV replication. HBV covalently closed circular (ccc)DNA degradation or functional inactivation remains the holy grail for the achievement of HBV cure. Transcriptional fidelity is a cornerstone in the regulation of gene expression. Here, we demonstrate that two helicases, DDX5 and DDX17, inhibit recognition of the HBV polyadenylation signal and thereby transcriptional termination, thus decreasing HBV RNA stability and acting as restriction factors for efficient cccDNA transcription and viral replication. The observation that DDX5 and DDX17 are downregulated in patients chronically infected with HBV suggests a role for these helicases in HBV persistence invivo. These results open new perspectives for researchers aiming at identifying new targets to neutralise cccDNA transcription.

An AluYa5 Insertion in the 3′UTR of COL4A1 and Cerebral Small Vessel Disease

Cerebral small vessel diseases (CSVDs) account for one-fifth of stroke cases. Numerous familial cases remain unresolved after routine screening of known CSVD genes. To identify novel genes and mechanisms associated with familial CSVD. This 2-stage study involved linkage analysis and a case-control study; linkage analysis and whole exome and genome sequencing were used to identify candidate gene variants in 2 large families with CSVD (9 patients with CSVD). Then, a case-control analysis was conducted on 246 unrelated probands, including probands from these 2 families and 244 additional probands. All probands (clinical onset <age 55 years and ≥1 first-degree relative with CSVD) were referred to the French cerebrovascular referral center between 2013 and 2023. The large-scale gnomAD structural variant database and 467 healthy individuals of French ancestry were used as a control group. A pathogenic AluYa5 insertion was identified within the COL4A1 3'UTR in the 2 large families with CSVD. Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), Western blot, and long-read RNA sequencing were used to investigate outcomes associated with the insertion using patient fibroblasts. Clinical and magnetic resonance imaging features of probands with variants and available relatives were assessed. Among 246 probands (141 females [57.3%]; median [IQR] age at referral, 56 [49-64] years), 7 patients of French ancestry carried the insertion. This insertion was absent in 467 healthy French individuals in a control group (odds ratio, ∞; 95% CI, 2.78 to ∞; P = 5 × 10-4) and 10 847 individuals from the gnomAD structural variant database (odds ratio, ∞; 95% CI, 64.77 to ∞; P = 2.42 × 10-12). In these 7 patients' families, 19 family members with CSVD carried the insertion. RT-qPCR and Western blot showed an upregulation of COL4A1 mRNA (10.6-fold increase; 95% CI, 1.4-fold to 17.1-fold increase) and protein levels (2.8-fold increase; 95% CI, 2.1-fold to 3.5-fold increase) in patient vs control group fibroblasts. Long-read RNA sequencing data showed that the insertion was associated with perturbation in the use of canonical COL4A1 polyadenylation signals (approximately 87% of isoforms transcribed from the wild type allele vs 5% of isoforms transcribed from the allele with the insertion used the 2 distal canonical polyadenylation signals). The main clinical feature of individuals with CSVD was the recurrence of pontine ischemic lesions starting at an early age (17 of 19 patients [89.5%]). This study found a novel mechanism associated with COL4A1 upregulation and a highly penetrant adult-onset CSVD. These findings suggest that quantitative alterations of the cerebrovascular matrisome are associated with CSVD pathogenesis, with diagnostic and therapeutic implications.

HPV oncogenes expressed from only one of multiple integrated HPV DNA copies drive clonal cell expansion in cervical cancer.

Persistent oncogenic HPV infections lead to viral DNA integration into the human genome and the development of cervical, anogenital, and oropharyngeal cancers. The expression of the viral E6 and E7 oncogenes plays a key role in cell transformation and tumorigenesis. However, how E6 and E7 could be expressed from the integrated viral DNA which often lacks a viral polyadenylation signal in the cancer cells remains unknown. By analyzing the integrated HPV DNA sites and expressed HPV RNAs in cervical cancer tissues and cell lines, we show that HPV oncogenes are expressed from only one of multiple chromosomal HPV DNA integrated copies. A host polyadenylation signal downstream of the integrated viral DNA is used for polyadenylation and stabilization of the virus-host chimeric RNAs, making the oncogenic transcripts targetable by siRNAs. This observation provides further understanding of the tumorigenic mechanism of HPV integration and suggests possible therapeutic strategies for the development of precision medicine for HPV cancers.

Splicing Switching of Alternative Last Exons Due to a Deletion Including Canonical Polyadenylation Site in COL6A2 Gene Causes Recessive UCMD.

Collagen VI-related myopathy spans a clinical continuum from severe Ullrich congenital muscular dystrophy to milder Bethlem myopathy, caused by genetic variants in COL6A1, COL6A2, and COL6A3 genes. Our objective was to report a newly identified patient with the pathogenic variants restricted to a polyadenylation signal in the 3'-untranslated region, which have not been reported in hereditary muscle disease. We performed clinicopathologic diagnosis and analysis using whole-genome and RNA sequencing. We report Ullrich congenital muscular dystrophy caused by a homozygous deletion, c.*198_*466del, which includes a polyadenylation signal in the canonical last exon of the COL6A2 gene. The parents were consanguineously married and had the heterozygous variant, but they were completely asymptomatic. In the patient's muscles, collagen VI was deficient in the sarcolemma, but present in the interstitium, showing the pattern of sarcolemma-specific collagen VI deficiency rather than a pattern of complete deficiency despite the lack of a polyadenylation signal. The RNA sequencing of the patient's muscle showed that alternative last exons were raised in COL6A2 transcript. Our case provides a valuable example of the mechanism of alternative splicing switches for polyadenylation selection.