Wolfgang Kaisers

An intronic G-run within HIV-1 intron 2 is critical for splicing regulation of vif-mRNA

ABSTRACT Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3′ splice site (3′ss) A1 but lack splicing at 5′ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3′ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3′ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5′ss D2. Here we show that an intronic G run (GI2-1) represses the use of a second 5′ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of GI2-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

Identification of amino acid determinants in CYP4B1 for optimal catalytic processing of 4-ipomeanol.

Mammalian CYP4B1 enzymes are cytochrome P450 mono-oxygenases that are responsible for the bioactivation of several exogenous pro-toxins including 4-ipomeanol (4-IPO). In contrast with the orthologous rabbit enzyme, we show here that native human CYP4B1 with a serine residue at position 427 is unable to bioactivate 4-IPO and does not cause cytotoxicity in HepG2 cells and primary human T-cells that overexpress these enzymes. We also demonstrate that a proline residue in the meander region at position 427 in human CYP4B1 and 422 in rabbit CYP4B1 is important for protein stability and rescues the 4-IPO bioactivation of the human enzyme, but is not essential for the catalytic activity of the rabbit CYP4B1 protein. Systematic substitution of native and p.S427P human CYP4B1 with peptide regions from the highly active rabbit enzyme reveals that 18 amino acids in the wild-type rabbit CYP4B1 protein are key for conferring high 4-IPO metabolizing activity. Introduction of 12 of the 18 amino acids that are also present at corresponding positions in other human CYP4 family members into the p.S427P human CYP4B1 protein results in a mutant human enzyme (P+12) that is as stable and as active as the rabbit wild-type CYP4B1 protein. These 12 mutations cluster in the predicted B-C loop through F-helix regions and reveal new amino acid regions important to P450 enzyme stability. Finally, by minimally re-engineering the human CYP4B1 enzyme for efficient activation of 4-IPO, we have developed a novel human suicide gene system that is a candidate for adoptive cellular therapies in humans.

Succession of splicing regulatory elements determines cryptic 5′ss functionality

Abstract 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.

Age, gender and UV-exposition related effects on gene expression in in vivo aged short term cultivated human dermal fibroblasts

Ageing, the progressive functional decline of virtually all tissues, affects numerous living organisms. Main phenotypic alterations of human skin during the ageing process include reduced skin thickness and elasticity which are related to extracellular matrix proteins. Dermal fibroblasts, the main source of extracellular fibrillar proteins, exhibit complex alterations during in vivo ageing and any of these are likely to be accompanied or caused by changes in gene expression. We investigated gene expression of short term cultivated in vivo aged human dermal fibroblasts using RNA-seq. Therefore, fibroblast samples derived from unaffected skin were obtained from 30 human donors. The donors were grouped by gender and age (Young: 19 to 25 years, Middle: 36 to 45 years, Old: 60 to 66 years). Two samples were taken from each donor, one from a sun-exposed and one from a sun-unexposed site. In our data, no consistently changed gene expression associated with donor age can be asserted. Instead, highly correlated expression of a small number of genes associated with transforming growth factor beta signalling was observed. Also, known gene expression alterations of in vivo aged dermal fibroblasts seem to be non-detectable in cultured fibroblasts.

rbamtools: an R interface to samtools enabling fast accumulative tabulation of splicing events over multiple RNA-seq samples.

The open source environment R isf the most widely used software to statistically explore biological data sets including sequence alignments. BAM is the de facto standard file format for sequence alignment. With rbamtools, we provide now a full spectrum of accessibility to BAM for R users such as reading, writing, extraction of subsets and plotting of alignment depth where the script syntax closely follows the SAM/BAM format. Additionally, rbamtools enables fast accumulative tabulation of splicing events over multiple BAM files.

Differences in the Early Development of Human and Mouse Embryonic Stem Cells

We performed a systematic analysis of gene expression features in early (10–21 days) development of human vs mouse embryonic cells (hESCs vs mESCs). Many development features were found to be conserved, and a majority of differentially regulated genes have similar expression change in both organisms. The similarity is especially evident, when gene expression profiles are clustered together and properties of clustered groups of genes are compared. First 10 days of mESC development match the features of hESC development within 21 days, in accordance with the differences in population doubling time in human and mouse ESCs. At the same time, several important differences are seen. There is a clear difference in initial expression change of transcription factors and stimulus responsive genes, which may be caused by the difference in experimental procedures. However, we also found that some biological processes develop differently; this can clearly be shown, for example, for neuron and sensory organ development. Some groups of genes show peaks of the expression levels during the development and these peaks cannot be claimed to happen at the same time points in the two organisms, as well as for the same groups of (orthologous) genes. We also detected a larger number of upregulated genes during development of mESCs as compared to hESCs. The differences were quantified by comparing promoters of related genes. Most of gene groups behave similarly and have similar transcription factor (TF) binding sites on their promoters. A few groups of genes have similar promoters, but are expressed differently in two species. Interestingly, there are groups of genes expressed similarly, although they have different promoters, which can be shown by comparing their TF binding sites. Namely, a large group of similarly expressed cell cycle-related genes is found to have discrepant TF binding properties in mouse vs human.

Sample Size Estimation for Detection of Splicing Events in Transcriptome Sequencing Data

Merging data from multiple samples is required to detect low expressed transcripts or splicing events that might be present only in a subset of samples. However, the exact number of required replicates enabling the detection of such rare events often remains a mystery but can be approached through probability theory. Here, we describe a probabilistic model, relating the number of observed events in a batch of samples with observation probabilities. Therein, samples appear as a heterogeneous collection of events, which are observed with some probability. The model is evaluated in a batch of 54 transcriptomes of human dermal fibroblast samples. The majority of putative splice-sites (alignment gap-sites) are detected in (almost) all samples or only sporadically, resulting in an U-shaped pattern for observation probabilities. The probabilistic model systematically underestimates event numbers due to a bias resulting from finite sampling. However, using an additional assumption, the probabilistic model can predict observed event numbers within a <10% deviation from the median. Single samples contain a considerable amount of uniquely observed putative splicing events (mean 7122 in alignments from TopHat alignments and 86,215 in alignments from STAR). We conclude that the probabilistic model provides an adequate description for observation of gap-sites in transcriptome data. Thus, the calculation of required sample sizes can be done by application of a simple binomial model to sporadically observed random events. Due to the large number of uniquely observed putative splice-sites and the known stochastic noise in the splicing machinery, it appears advisable to include observation of rare splicing events into analysis objectives. Therefore, it is beneficial to take scores for the validation of gap-sites into account.

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.

Biomechanical assessment of healthy and keratoconic corneas (with/without crosslinking) using dynamic ultrahigh-speed Scheimpflug technology and the relevance of the parameter (A1L−A2L)

Aims To examine corneal biomechanics in healthy and keratoconic eyes, with or without crosslinking obtained by ultrahigh-speed Scheimpflug measurements (Corvis ST). Methods One hundred and seventeen eyes were studied in three groups: group 1 (n=39) contained keratoconic eyes without crosslinking. Group 2 (CXL; n=28) comprised keratoconic eyes after crosslinking. These were compared with a control group (n=50 matched healthy eyes). In addition, 10 keratoconus patients, before and after CXL treatment, respectively, were examined. Results The novel parameter A1L–A2L demonstrated highly significant differences between crosslinked corneas and untreated keratoconic or healthy corneas. Velocity during second applanation (A2V) and deformation amplitude (DA) were significantly increased in crosslinked keratoconic eyes both compared with untreated keratoconic eyes and with healthy controls. Radius at highest curvature also was significant among all groups. Inward applanation length (A1L) was significantly increased in controls, whereas outward applanation length (A2L) was significantly reduced in crosslinked keratoconic eyes compared with both other groups. The follow-up analysis revealed statistically significant changes in pachymetry and intraocular pressure and showed tendencies towards significance in applanation times 1 and 2 and in DA. Conclusions Both A2V and A2L are viable parameters to discriminate healthy from keratoconic but also crosslinked from non-crosslinked keratoconic corneas. The difference of A1L−A2L could reliably discriminate crosslinked from non-crosslinked and healthy corneas. Follow-up examination in a small cohort allows distinction between crosslinked and untreated keratoconus in follow-up examinations. The difference of A1L-A2L could reliably discriminate crosslinked from non-crosslinked and healthy corneas. Measurements of corneal deformation using dynamic ultrahigh-speed Scheimpflug technology are reproducible and provide useful information about keratoconus assessment and biomechanics. Therefore, the Corvis ST seems to provide useful technology to monitor therapeutic success of crosslinking treatment.

Immunohistochemical Detection of Lymph Node‐DTCs in Patients with Node‐Negative HNSCC

This study was performed to systematically assess the prevalence, topography and prognostic impact of disseminated tumor cells (DTCs) in lymph nodes (LN) of patients with primary, regional and distant metastasis‐free head and neck squamous cell carcinoma (HNSCC) who underwent resection with elective neck dissection. From the routinely processed resection specimen, we could prospectively analyze a total of 1.137 exactly mapped LNs of 50 pN0‐HNSCC patients, classified as tumor free by routine histopathology. Three immunohistochemistry (IHC) assays using antibodies directed against CK5/14, a broad spectrum of CKs (1–8, 10, 14–16 and 19), and CD44v6, respectively, were applied on 4.190 LN sections to detect DTCs. The IHC results were correlated with clinicopathologic parameters and clinical follow‐up data. We detected seven micrometastases (MM) in five patients and 31 DTCs in 12 patients. Overall, 15 (30%) patients were positive for DTCs or MMs. Strikingly, the anatomical distribution of LN affected with DTCs was not random, but was dependent on the lateralization of the primary tumor and clustered significantly most proximal to the primary tumor. None of the investigated patients developed loco‐regional lymphatic or distant metastasis during the mean follow‐up period of 71 months. Our results reveal clinically occult tumor cell dissemination as an early and frequent event in HNSCC. Considering that higher rates of recurrences in therapeutic LN dissection concepts have been reported than in elective neck dissection strategies, our DTC‐data support to perform elective neck dissections, since they appear to be effective in preventing loco‐regional lymphatic recurrence from LN DTCs or MMs.