Sequence diversity of hepatitis D virus in Mongolia.

Introduction: The Hepatitis Delta Virus (HDV) is a defective, single-stranded RNA virusoid encoding for a single protein, the Hepatitis Delta Antigen (HDAg), which requires the hepatitis B virus (HBV) envelope protein (HBsAg) for its transmission. Currently, hepatitis D is the most aggressive form of viral hepatitis and treatment options are limited. Worldwide 12 million people are chronically infected with HDV being at high risk for progression to cirrhosis and development of liver cancer. Objectives: Although it is well established that Mongolia is the country with the highest prevalence of HDV infections, the information on the molecular epidemiology and factors contributing to HDV sequence diversity are largely unclear. The aim of the study was to characterize the sequence diversity of HDV in rural areas from Mongolia and to determine the extent of HLA class I-associated selection pressure. Patients and methods: From the HepMongolia cohort from rural areas in Mongolia, 451 HBsAg-positive individuals were selected and anti-HDV, HDV-RNA and the sequence of the large HDAg was determined. For all individuals the HLA class I locus was genotyped. Residues under selection pressure in the presence of individual HLA class I types were identified with the recently published analysis tool HAMdetector. Results: Of 431 HBsAg positive patients, 281 were anti-HDV positive (65%), and HDV-RNA could be detected in 207 of 281 (74%) of patients. The complete large HDAg was successfully sequenced from 131 samples. Phylogenetic analysis revealed that all Mongolian HDV isolates belong to genotype 1, however, they separate into several different clusters without clear regional association. In turn, from phylogeny there is strong evidence for recent local transmission events. Importantly, we found multiple residues with strong support for HLA class I-associated selection pressure consistent with a functional CD8+ T cell response directed against HDV. Conclusion: HDV isolates from Mongolia are highly diverse. The molecular epidemiology suggests circulation of multiple subtypes and provides evidence for ongoing recent transmissions.

A novel cancer risk prediction score for the natural course of FA patients with biallelic BRCA2/FANCD1 mutations.

Biallelic germline mutations in BRCA2 occur in the Fanconi anemia (FA)-D1 subtype of the rare pediatric disorder, FA, characterized clinically by severe congenital abnormalities and a very high propensity to develop malignancies early in life. Clinical and genetic data from 96 FA-D1 patients with biallelic BRCA2 mutations were collected and used to develop a new cancer risk prediction score system based on the specific mutations in BRCA2. This score takes into account the location of frameshift/stop and missense mutations relative to exon 11 of BRCA2, which encodes the major sites for interaction with the RAD51 recombinase, and uses the MaxEnt and HBond splicing scores to analyze potential splice site perturbations. Among 75 FA-D1 patients with ascertained BRCA2 mutations, 66 patients developed 102 malignancies, ranging from one to three independent tumors per individual. The median age at the clinical presentation of peripheral embryonal tumors was 1.0, at the onset of hematologic malignancies 1.8 and at the manifestation of CNS tumors 2.7 years, respectively. Patients who received treatment lived longer than those without. Using our novel scoring system, we could distinguish three distinct cancer risk groups among FA-D1 patients: in the first, patients developed their initial malignancy at a median age of 1.3 years (n = 36, 95% CI = 0.9-1.8), in the second group at 2.3 years (n = 17, 95% CI = 1.4-4.4) and in the third group at 23.0 years (n = 22, 95% CI = 4.3-n/a). Therefore, this scoring system allows, for the first time, to predict the cancer manifestation of FA-D1 patients simply based on the type and position of the mutations in BRCA2.

Adjusted COVID-19 booster schedules balance age-dependent differences in antibody titers benefitting risk populations.

We provide follow-up data on the humoral immune response after COVID-19 vaccinations of two distinct cohorts aged below 60 and over 80 years to screen for age-related differences in the longevity and magnitude of the induction of the antibody responses post booster-vaccinations. While anti-SARS-CoV-2 spike-specific IgG and neutralization capacity waned rapidly after the initial vaccination schedule, additional boosters highly benefitted the humoral immune responses especially in the elderly cohort, including the neutralization of Omikron variants. Thus, adjusted COVID-19 booster vaccination schedules are an appropriate tool to overcome limitations in the success of vaccinations.

Modeling splicing outcome by combining 5'ss strength and splicing regulatory elements.

Correct pre-mRNA processing in higher eukaryotes vastly depends on splice site recognition. Beyond conserved 5'ss and 3'ss motifs, splicing regulatory elements (SREs) play a pivotal role in this recognition process. Here, we present in silico designed sequences with arbitrary a priori prescribed splicing regulatory HEXplorer properties that can be concatenated to arbitrary length without changing their regulatory properties. We experimentally validated in silico predictions in a massively parallel splicing reporter assay on more than 3000 sequences and exemplarily identified some SRE binding proteins. Aiming at a unified 'functional splice site strength' encompassing both U1 snRNA complementarity and impact from neighboring SREs, we developed a novel RNA-seq based 5'ss usage landscape, mapping the competition of pairs of high confidence 5'ss and neighboring exonic GT sites along HBond and HEXplorer score coordinate axes on human fibroblast and endothelium transcriptome datasets. These RNA-seq data served as basis for a logistic 5'ss usage prediction model, which greatly improved discrimination between strong but unused exonic GT sites and annotated highly used 5'ss. Our 5'ss usage landscape offers a unified view on 5'ss and SRE neighborhood impact on splice site recognition, and may contribute to improved mutation assessment in human genetics.

Selenoprotein T Protects Endothelial Cells against Lipopolysaccharide-Induced Activation and Apoptosis.

Sepsis is an exaggerated immune response upon infection with lipopolysaccharide (LPS) as the main causative agent. LPS-induced activation and apoptosis of endothelial cells (EC) can lead to organ dysfunction and finally organ failure. We previously demonstrated that the first twenty amino acids of the Apurinic/Apyrimidinic Endodeoxyribonuclease 1 (APEX1) are sufficient to inhibit EC apoptosis. To identify genes whose regulation by LPS is affected by this N-terminal APEX1 peptide, EC were transduced with an expression vector for the APEX1 peptide or an empty control vector and treated with LPS. Following RNA deep sequencing, genes upregulated in LPS-treated EC expressing the APEX1 peptide were identified bioinformatically. Selected candidates were validated by semi-quantitative real time PCR, a promising one was Selenoprotein T (SELENOT). For functional analyses, an expression vector for SELENOT was generated. To study the effect of SELENOT expression on LPS-induced EC activation and apoptosis, the SELENOT vector was transfected in EC. Immunostaining showed that SELENOT was expressed and localized in the ER. EC transfected with the SELENOT plasmid showed no activation and reduced apoptosis induced by LPS. SELENOT as well as APEX1(1-20) can protect EC against activation and apoptosis and could provide new therapeutic approaches in the treatment of sepsis.

Cilium induction triggers differentiation of glioma stem cells.

Glioblastoma multiforme (GBM) possesses glioma stem cells (GSCs) that promote self-renewal, tumor propagation, and relapse. Understanding the mechanisms of GSCs self-renewal can offer targeted therapeutic interventions. However, insufficient knowledge of GSCs' fundamental biology is a significant bottleneck hindering these efforts. Here, we show that patient-derived GSCs recruit elevated levels of proteins that ensure the temporal cilium disassembly, leading to suppressed ciliogenesis. Depleting the cilia disassembly complex components is sufficient to induce ciliogenesis in a subset of GSCs via relocating platelet-derived growth factor receptor-alpha (PDGFR-α) to a newly induced cilium. Importantly, restoring ciliogenesis enabled GSCs to switch from self-renewal to differentiation. Finally, using an organoid-based glioma invasion assay and brain xenografts in mice, we establish that ciliogenesis-induced differentiation can prevent the infiltration of GSCs into the brain. Our findings illustrate a role for cilium as a molecular switch in determining GSCs' fate and suggest cilium induction as an attractive strategy to intervene in GSCs proliferation.

Modifying splice site usage with ModCon: Maintaining the genetic code while changing the underlying mRNP code.

Codon degeneracy of amino acid sequences permits an additional "mRNP code" layer underlying the genetic code that is related to RNA processing. In pre-mRNA splicing, splice site usage is determined by both intrinsic strength and sequence context providing RNA binding sites for splicing regulatory proteins. In this study, we systematically examined modification of splicing regulatory properties in the neighborhood of a GT site, i.e. potential splice site, without altering the encoded amino acids. We quantified the splicing regulatory properties of the neighborhood around a potential splice site by its Splice Site HEXplorer Weight (SSHW) based on the HEXplorer score algorithm. To systematically modify GT site neighborhoods, either minimizing or maximizing their SSHW, we designed the novel stochastic optimization algorithm ModCon that applies a genetic algorithm with stochastic crossover, insertion and random mutation elements supplemented by a heuristic sliding window approach. To assess the achievable range in SSHW in human splice donors without altering the encoded amino acids, we applied ModCon to a set of 1000 randomly selected Ensembl annotated human splice donor sites, achieving substantial and accurate changes in SSHW. Using ModCon optimization, we successfully switched splice donor usage in a splice site competition reporter containing coding sequences from FANCA, FANCB or BRCA2, while retaining their amino acid coding information. The ModCon algorithm and its R package implementation can assist in reporter design by either introducing novel splice sites, silencing accidental, undesired splice sites, and by generally modifying the entire mRNP code while maintaining the genetic code.

Let It Go: HIV-1 cis-Acting Repressive Sequences.

After human immunodeficiency virus type 1 (HIV-1) was identified in the early 1980s, intensive work began to understand the molecular basis of HIV-1 gene expression. Subgenomic HIV-1 RNA regions, spread throughout the viral genome, were described to have a negative impact on the nuclear export of some viral transcripts. Those studies revealed an intrinsic RNA code as a new form of nuclear export regulation. Since such regulatory regions were later also identified in other viruses, as well as in cellular genes, it can be assumed that, during evolution, viruses took advantage of them to achieve more sophisticated replication mechanisms. Here, we review HIV-1 cis-acting repressive sequences that have been identified, and we discuss their possible underlying mechanisms and importance. Additionally, we show how current bioinformatic tools might allow more predictive approaches to identify and investigate them.

Age-dependent Immune Response to the Biontech/Pfizer BNT162b2 Coronavirus Disease 2019 Vaccination.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has led to the development of various vaccines. Real-life data on immune responses elicited in the most vulnerable group of vaccinees older than age 80 years old are still underrepresented despite the prioritization of the elderly in vaccination campaigns. METHODS: We conducted a cohort study with 2 age groups, young vaccinees below the age of 60 years and elderly vaccinees over the age of 80 years, to compare their antibody responses to the first and second dose of the BNT162b2 coronavirus disease 2019 vaccination. RESULTS: Although the majority of participants in both groups produced specific immunoglobulin G antibody titers against SARS-CoV-2 spike protein, titers were significantly lower in elderly participants. Although the increment of antibody levels after the second immunization was higher in elderly participants, the absolute mean titer of this group remained lower than the <60 years of age group. After the second vaccination, 31.3% of the elderly had no detectable neutralizing antibodies in contrast to the younger group, in which only 2.2% had no detectable neutralizing antibodies. CONCLUSIONS: Our data showed differences between the antibody responses raised after the first and second BNT162b2 vaccination, in particular lower frequencies of neutralizing antibodies in the elderly group. This suggests that this population needs to be closely monitored and may require earlier revaccination and/or an increased vaccine dose to ensure stronger long-lasting immunity and protection against infection.

Profiling of cis- and trans-acting factors supporting noncanonical splice site activation.

Recently, by combining transcriptomics with functional splicing reporter assays we were able to identify GT > GC > TT as the three highest ranked dinucleotides of human 5' splice sites (5'ss). Here, we have extended our investigations to the proteomic characterization of nuclear proteins that bind to canonical and noncanonical 5'ss. Surprisingly, we found that U1 snRNP binding to functional 5'ss sequences prevented components of the DNA damage response (DDR) from binding to the RNA, suggesting a close link between spliceosome arrangement and genome stability. We demonstrate that all tested noncanonical 5'ss sequences are bona-fide targets of the U2-type spliceosome and are bound by U1 snRNP, including U1-C, in the presence of splicing enhancers. The quantity of precipitated U1-C protein was similar for all noncanonical 5'ss dinucleotides, so that the highly different 5'ss usage was likely due to a later step after early U1 snRNP binding. In addition, we show that an internal GT at positions +5/+6 can be advantageous for splicing at position +1 of noncanonical splice sites. Likewise, and in agreement with previous observations, splicing inactive U1 snRNP binding sites could serve as splicing enhancers, which may also explain the higher abundance of U1 snRNPs compared to other U snRNPs. Finally, we observe that an arginine-serine (RS)-rich domain recruitment to stem loop I of the U1 snRNA is functionally sufficient to promote exon-definition and upstream 3'ss activation.

Altered HIV-1 mRNA Splicing Due to Drug-Resistance-Associated Mutations in Exon 2/2b.

The underlying molecular mechanism and their general effect on the replication capacity of HIV 1 drug-resistance-associated mutations is often poorly understood. To elucidate the effect of two such mutations located in a region with a high density of spicing regulatory elements on the HIV-1-splicing outcome, bioinformatic predictions were combined with transfection and infection experiments. Results show that the previously described R263K drug-resistance-associated integrase mutation has additionally a severe effect on the ESE2b splicing regulatory element (SRE) in exon 2b, which causes loss of SD2b recognition. This was confirmed by an R263R silent mutation with a similar predicted effect on the exon 2b SRE. In contrast, a V260I mutation and its silent counterpart with a lower effect on ESS2b did not exhibit any differences in the splicing pattern. Since HIV-1 highly relies on a balanced splicing reaction, changes in the splicing outcome can contribute to changes in viral replication and might add to the effect of escape mutations toward antiviral drugs. Thus, a classification of mutations purely addressing proteins is insufficient.

VarCon: An R Package for Retrieving Neighboring Nucleotides of an SNV.

Reporting of a single nucleotide variant (SNV) follows the Sequence Variant Nomenclature (http://varnomen.hgvs.org/), using an unambiguous numbering scheme specific for coding and noncoding DNA. However, the corresponding sequence neighborhood of a given SNV, which is required to assess its impact on splicing regulation, is not easily accessible from this nomenclature. Providing fast and easy access to this neighborhood just from a given SNV reference, the novel tool VarCon combines information of the Ensembl human reference genome and the corresponding transcript table for accurate retrieval. VarCon also displays splice site scores (HBond and MaxEnt scores) and HEXplorer profiles of an SNV neighborhood, reflecting position-dependent splice enhancing and silencing properties.

High Concentration of Low-Density Lipoprotein Results in Disturbances in Mitochondrial Transcription and Functionality in Endothelial Cells.

Concentrations of low-density lipoprotein (LDL) above 0.8 mg/ml have been associated with increased risk for cardiovascular diseases and impaired endothelial functionality. Here, we demonstrate that high concentrations of LDL (1 mg/ml) decreased NOS3 protein and RNA levels in primary human endothelial cells. In addition, RNA sequencing data, in particular splice site usage analysis, showed a shift in NOS3 exon-exon junction reads towards those specifically assigned to nonfunctional transcript isoforms further diminishing the functional NOS3 levels. The reduction in NOS3 was accompanied by decreased migratory capacity, which depends on intact mitochondria and ATP formation. In line with these findings, we also observed a reduced ATP content. While mitochondrial mass was unaffected by high LDL, we found an increase in mitochondrial DNA copy number and mitochondrial RNA transcripts but decreased expression of nuclear genes coding for respiratory chain proteins. Therefore, high LDL treatment most likely results in an imbalance between respiratory chain complex proteins encoded in the mitochondria and in the nucleus resulting in impaired respiratory chain function explaining the reduction in ATP content. In conclusion, high LDL treatment leads to a decrease in active NOS3 and dysregulation of mitochondrial transcription, which is entailed by reduced ATP content and migratory capacity and thus, impairment of endothelial cell functionality.

Context matters: Regulation of splice donor usage.

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Ranking noncanonical 5' splice site usage by genome-wide RNA-seq analysis and splicing reporter assays.

Most human pathogenic mutations in 5' splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5'ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5'ss are used. In this work, we systematically examined noncanonical 5' splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5'ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5'ss. In particular, we could confirm splicing at different positions (i.e., -1, +1, +5) of a splice site for all noncanonical dinucleotides "weaker" than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5'ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5'ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking.

Validation of Splicing Events in Transcriptome Sequencing Data.

Genomic alignments of sequenced cellular messenger RNA contain gapped alignments which are interpreted as consequence of intron removal. The resulting gap-sites, genomic locations of alignment gaps, are landmarks representing potential splice-sites. As alignment algorithms report gap-sites with a considerable false discovery rate, validations are required. We describe two quality scores, gap quality score (gqs) and weighted gap information score (wgis), developed for validation of putative splicing events: While gqs solely relies on alignment data wgis additionally considers information from the genomic sequence. FASTQ files obtained from 54 human dermal fibroblast samples were aligned against the human genome (GRCh38) using TopHat and STAR aligner. Statistical properties of gap-sites validated by gqs and wgis were evaluated by their sequence similarity to known exon-intron borders. Within the 54 samples, TopHat identifies 1,000,380 and STAR reports 6,487,577 gap-sites. Due to the lack of strand information, however, the percentage of identified GT-AG gap-sites is rather low. While gap-sites from TopHat contain ≈89% GT-AG, gap-sites from STAR only contain ≈42% GT-AG dinucleotide pairs in merged data from 54 fibroblast samples. Validation with gqs yields 156,251 gap-sites from TopHat alignments and 166,294 from STAR alignments. Validation with wgis yields 770,327 gap-sites from TopHat alignments and 1,065,596 from STAR alignments. Both alignment algorithms, TopHat and STAR, report gap-sites with considerable false discovery rate, which can drastically be reduced by validation with gqs and wgis.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality.

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5' end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bß-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3' exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3' exon end definition.