TISIGNER.com: web services for improving recombinant protein production.

Experiments that are planned using accurate prediction algorithms will mitigate failures in recombinant protein production. We have developed TISIGNER (https://tisigner.com) with the aim of addressing technical challenges to recombinant protein production. We offer three web services, TIsigner (Translation Initiation coding region designer), SoDoPE (Soluble Domain for Protein Expression) and Razor, which are specialised in synonymous optimisation of recombinant protein expression, solubility and signal peptide analysis, respectively. Importantly, TIsigner, SoDoPE and Razor are linked, which allows users to switch between the tools when optimising genes of interest.

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species.

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron-exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

Analysis of 11,430 recombinant protein production experiments reveals that protein yield is tunable by synonymous codon changes of translation initiation sites.

Recombinant protein production is a key process in generating proteins of interest in the pharmaceutical industry and biomedical research. However, about 50% of recombinant proteins fail to be expressed in a variety of host cells. Here we show that the accessibility of translation initiation sites modelled using the mRNA base-unpairing across the Boltzmann's ensemble significantly outperforms alternative features. This approach accurately predicts the successes or failures of expression experiments, which utilised Escherichia coli cells to express 11,430 recombinant proteins from over 189 diverse species. On this basis, we develop TIsigner that uses simulated annealing to modify up to the first nine codons of mRNAs with synonymous substitutions. We show that accessibility captures the key propensity beyond the target region (initiation sites in this case), as a modest number of synonymous changes is sufficient to tune the recombinant protein expression levels. We build a stochastic simulation model and show that higher accessibility leads to higher protein production and slower cell growth, supporting the idea of protein cost, where cell growth is constrained by protein circuits during overexpression.

Solubility-Weighted Index: fast and accurate prediction of protein solubility.

MOTIVATION: Recombinant protein production is a widely used technique in the biotechnology and biomedical industries, yet only a quarter of target proteins are soluble and can therefore be purified. RESULTS: We have discovered that global structural flexibility, which can be modeled by normalized B-factors, accurately predicts the solubility of 12 216 recombinant proteins expressed in Escherichia coli. We have optimized these B-factors, and derived a new set of values for solubility scoring that further improves prediction accuracy. We call this new predictor the 'Solubility-Weighted Index' (SWI). Importantly, SWI outperforms many existing protein solubility prediction tools. Furthermore, we have developed 'SoDoPE' (Soluble Domain for Protein Expression), a web interface that allows users to choose a protein region of interest for predicting and maximizing both protein expression and solubility. AVAILABILITY AND IMPLEMENTATION: The SoDoPE web server and source code are freely available at https://tisigner.com/sodope and https://github.com/Gardner-BinfLab/TISIGNER-ReactJS, respectively. The code and data for reproducing our analysis can be found at https://github.com/Gardner-BinfLab/SoDoPE_paper_2020. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A crowdsourced set of curated structural variants for the human genome.

A high quality benchmark for small variants encompassing 88 to 90% of the reference genome has been developed for seven Genome in a Bottle (GIAB) reference samples. However a reliable benchmark for large indels and structural variants (SVs) is more challenging. In this study, we manually curated 1235 SVs, which can ultimately be used to evaluate SV callers or train machine learning models. We developed a crowdsourcing app-SVCurator-to help GIAB curators manually review large indels and SVs within the human genome, and report their genotype and size accuracy. SVCurator displays images from short, long, and linked read sequencing data from the GIAB Ashkenazi Jewish Trio son [NIST RM 8391/HG002]. We asked curators to assign labels describing SV type (deletion or insertion), size accuracy, and genotype for 1235 putative insertions and deletions sampled from different size bins between 20 and 892,149 bp. 'Expert' curators were 93% concordant with each other, and 37 of the 61 curators had at least 78% concordance with a set of 'expert' curators. The curators were least concordant for complex SVs and SVs that had inaccurate breakpoints or size predictions. After filtering events with low concordance among curators, we produced high confidence labels for 935 events. The SVCurator crowdsourced labels were 94.5% concordant with the heuristic-based draft benchmark SV callset from GIAB. We found that curators can successfully evaluate putative SVs when given evidence from multiple sequencing technologies.

Systematic comparison of plasma EBV DNA, anti-EBV antibodies and miRNA levels for early detection and prognosis of nasopharyngeal carcinoma.

Nasopharyngeal carcinoma (NPC) is originated from the epithelial cells of nasopharynx, Epstein-Barr virus (EBV)-associated and has the highest incidence and mortality rates in Southeast Asia. Late presentation is a common issue and early detection could be the key to reduce the disease burden. Sensitivity of plasma EBV DNA, an established NPC biomarker, for Stage I NPC is controversial. Most newly reported NPC biomarkers have neither been externally validated nor compared to the established ones. This causes difficulty in planning for cost-effective early detection strategies. Our study systematically evaluated six established and four new biomarkers in NPC cases, population controls and hospital controls. We showed that BamHI-W 76 bp remains the most sensitive plasma biomarker, with 96.7% (29/30), 96.7% (58/60) and 97.4% (226/232) sensitivity to detect Stage I, early stage and all NPC, respectively. Its specificity was 94.2% (113/120) against population controls and 90.4% (113/125) against hospital controls. Diagnostic accuracy of BamHI-W 121 bp and ebv-miR-BART7-3p were validated. Hsa-miR-29a-3p and hsa-miR-103a-3p were not, possibly due to lower number of advanced stage NPC cases included in this subset. Decision tree modeling suggested that combination of BamHI-W 76 bp and VCA IgA or EA IgG may increase the specificity or sensitivity to detect NPC. EBNA1 99 bp could identify NPC patients with poor prognosis in early and advanced stage NPC. Our findings provided evidence for improvement in NPC screening strategies, covering considerations of opportunistic screening, combining biomarkers to increase sensitivity or specificity and testing biomarkers from single sampled specimen to avoid logistic problems of resampling.

The exon-intron gene structure upstream of the initiation codon predicts translation efficiency.

Introns in mRNA leaders are common in complex eukaryotes, but often overlooked. These introns are spliced out before translation, leaving exon-exon junctions in the mRNA leaders (leader EEJs). Our multi-omic approach shows that the number of leader EEJs inversely correlates with the main protein translation, as does the number of upstream open reading frames (uORFs). Across the five species studied, the lowest levels of translation were observed for mRNAs with both leader EEJs and uORFs (29%). This class of mRNAs also have ribosome footprints on uORFs, with strong triplet periodicity indicating uORF translation. Furthermore, the positions of both leader EEJ and uORF are conserved between human and mouse. Thus, the uORF, in combination with leader EEJ predicts lower expression for nearly one-third of eukaryotic proteins.

Hepatitis B virus nuclear export elements: RNA stem-loop α and ß, key parts of the HBV post-transcriptional regulatory element.

Many viruses contain RNA elements that modulate splicing and/or promote nuclear export of their RNAs. The RNAs of the major human pathogen, hepatitis B virus (HBV) contain a large (~600 bases) composite cis-acting 'post-transcriptional regulatory element' (PRE). This element promotes expression from these naturally intronless transcripts. Indeed, the related woodchuck hepadnavirus PRE (WPRE) is used to enhance expression in gene therapy and other expression vectors. These PRE are likely to act through a combination of mechanisms, including promotion of RNA nuclear export. Functional components of both the HBV PRE and WPRE are 2 conserved RNA cis-acting stem-loop (SL) structures, SLα and SLß. They are within the coding regions of polymerase (P) gene, and both P and X genes, respectively. Based on previous studies using mutagenesis and/or nuclear magnetic resonance (NMR), here we propose 2 covariance models for SLα and SLß. The model for the 30-nucleotide SLα contains a G-bulge and a CNGG(U) apical loop of which the first and the fourth loop residues form a CG pair and the fifth loop residue is bulged out, as observed in the NMR structure. The model for the 23-nucleotide SLß contains a 7-base-pair stem and a 9-nucleotide loop. Comparison of the models with other RNA structural elements, as well as similarity searches of human transcriptome and viral genomes demonstrate that SLα and SLß are specific to HBV transcripts. However, they are well conserved among the hepadnaviruses of non-human primates, the woodchuck and ground squirrel.

Recombinant production of Epstein-Barr virus BZLF1 trans-activator and characterization of its DNA-binding specificity.

This paper describes the recombinant production of a biologically active Epstein-Barr virus BZLF1 trans-activator, i.e., Z-encoded broadly reactive activator (ZEBRA), that recognized specific DNA motifs. We used auto-induction for histidine-tagged BZLF1 expression in Escherichia coli and immobilized cobalt affinity membrane chromatography for protein purification under native conditions. We obtained the purified BZLF1 at a yield of 5.4mg per gram of wet weight cells at 75% purity, in which 27% of the recombinant BZLF1 remained biologically active. The recombinant BZLF1 bound to oligonucleotides containing ZEBRA response elements, either AP-1 or ZIIIB, but not a ZIIIB mutant. The recombinant BZLF1 showed a specific DNA-binding activity which could be useful for functional studies.

Quantitative analysis of the splice variants expressed by the major hepatitis B virus genotypes.

Hepatitis B virus (HBV) is a major human pathogen that causes liver diseases. The main HBV RNAs are unspliced transcripts that encode the key viral proteins. Recent studies have shown that some of the HBV spliced transcript isoforms are predictive of liver cancer, yet the roles of these spliced transcripts remain elusive. Furthermore, there are nine major HBV genotypes common in different regions of the world, these genotypes may express different spliced transcript isoforms. To systematically study the HBV splice variants, we transfected human hepatoma cells, Huh7, with four HBV genotypes (A2, B2, C2 and D3), followed by deep RNA-sequencing. We found that 13-28 % of HBV RNAs were splice variants, which were reproducibly detected across independent biological replicates. These comprised 6 novel and 10 previously identified splice variants. In particular, a novel, singly spliced transcript was detected in genotypes A2 and D3 at high levels. The biological relevance of these splice variants was supported by their identification in HBV-positive liver biopsy and serum samples, and in HBV-infected primary human hepatocytes. Interestingly the levels of HBV splice variants varied across the genotypes, but the spliced pregenomic RNA SP1 and SP9 were the two most abundant splice variants. Counterintuitively, these singly spliced SP1 and SP9 variants had a suboptimal 5' splice site, supporting the idea that splicing of HBV RNAs is tightly controlled by the viral post-transcriptional regulatory RNA element.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs.

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

A new class of ribozyme from hepatitis B virus.

In this issue, Chakraborty and Ghosh present initial in vitro data indicating a 53-base part of the hepatitis B virus (HBV) epsilon region is a ribozyme with unusual properties. Self-cleavage first occurred at the HBV non-canonical polyadenylation signal UAUAAA and the "centre bulge" of the epsilon region, releasing a ribozyme. The released ribozyme then acted as a new class of trans-cleaving ribozyme cleaving other RNAs.

Inclusion bodies of recombinant Epstein-Barr virus capsid antigen p18 as potential immobilized antigens in enzyme immunoassays for detection of nasopharyngeal carcinoma.

BACKGROUND: Development of indirect enzyme-linked immunosorbent assays (ELISAs) often utilizes synthetic peptides or recombinant proteins from Escherichia coli as immobilized antigens. Because inclusion bodies (IBs) formed during recombinant protein expression in E. coli are commonly thought as misfolded aggregates, only refolded proteins from IBs are used to develop new or in-house diagnostic assays. However, the promising utilities of IBs as nanomaterials and immobilized enzymes as shown in recent studies have led us to explore the potential use of IBs of recombinant Epstein-Barr virus viral capsid antigen p18 (VCA p18) as immobilized antigens in ELISAs for serologic detection of nasopharyngeal carcinoma (NPC). METHODS: Thioredoxin fusion VCA p18 (VCA-Trx) and IBs of VCA p18 without fusion tags (VCA-IBs) were purified from E. coli. The diagnostic performances of IgG/VCA-IBs, IgG/VCA-Denat-IBs (using VCA-IBs coated in 8mol/l urea), IgG/VCA-Trx, and IgG/VCA-Peptide assays were compared by screening 100 NPC case-control pairs. RESULTS: The IgG/VCA-Denat-IBs assay showed the best area under the receiver operating characteristic curve (AUC: 0.802; p<0.05), while the AUCs for the IgG/VCA-IBs, IgG/VCA-Trx, and IgG/VCA-Peptide assays were comparable (AUC: 0.740, 0.727, and 0.741, respectively). CONCLUSION: We improved the diagnostic performance of the ELISA significantly using IBs of recombinant VCA p18.

R54C Mutation of NOTCH3 Gene in the First Rungus Family with CADASIL.

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a rare hereditary stroke caused by mutations in NOTCH3 gene. We report the first case of CADASIL in an indigenous Rungus (Kadazan-Dusun) family in Kudat, Sabah, Malaysia confirmed by a R54C (c.160C>T, p.Arg54Cys) mutation in the NOTCH3. This mutation was previously reported in a Caucasian and two Korean cases of CADASIL. We recruited two generations of the affected Rungus family (n = 9) and found a missense mutation (c.160C>T) in exon 2 of NOTCH3 in three siblings. Two of the three siblings had severe white matter abnormalities in their brain MRI (Scheltens score 33 and 50 respectively), one of whom had a young stroke at the age of 38. The remaining sibling, however, did not show any clinical features of CADASIL and had only minimal changes in her brain MRI (Scheltens score 17). This further emphasized the phenotype variability among family members with the same mutation in CADASIL. This is the first reported family with CADASIL in Rungus subtribe of Kadazan-Dusun ethnicity with a known mutation at exon 2 of NOTCH3. The penetrance of this mutation was not complete during the course of this study.

On optimizing the blocking step of indirect enzyme-linked immunosorbent assay for Epstein-Barr virus serology.

BACKGROUND: Because blocking agent occupies most binding surface of a solid phase, its ability to prevent nonspecific binding determines the signal-to-noise ratio (SNR) and reliability of an enzyme-linked immunosorbent assay (ELISA). METHODS: We demonstrate a stepwise approach to seek a compatible blocking buffer for indirect ELISA, via a case-control study (n=176) of Epstein-Barr virus (EBV)-associated nasopharyngeal carcinoma (NPC). RESULTS: Regardless of case-control status, we found that synthetic polymer blocking agents, mainly Ficoll and poly(vinyl alcohol) (PVA) were able to provide homogeneous backgrounds among samples, as opposed to commonly used blocking agents, notably nonfat dry milk (NFDM). The SNRs for NPC samples that correspond to blocking using PVA were approximately 3-fold, on average, higher than those blocking using NFDM. Both intra- and inter-assay precisions of PVA-based assays were <14%. CONCLUSION: A blocking agent of choice should have tolerable sample backgrounds from both cases and controls to ensure the reliability of an immunoassay.