GBSX: a toolkit for experimental design and demultiplexing genotyping by sequencing experiments

Background

Massive parallel sequencing is a powerful tool for variant discovery and genotyping. To reduce costs, sequencing of restriction enzyme based reduced representation libraries can be utilized. This technology is generally referred to as Genotyping By Sequencing (GBS). To deal with GBS experimental design and initial processing specific bioinformatic tools are needed.

Results

GBSX is a package that assists in selecting the appropriate enzyme and the design of compatible in-line barcodes. Post sequencing, it performs optimized demultiplexing using these barcodes to create fastq files per barcode which can easily be plugged into existing variant analysis pipelines. Here we demonstrate the usability of the GBSX toolkit and demonstrate improved in-line barcode demultiplexing and trimming performance compared to existing tools.

Conclusions

GBSX provides an easy to use suite of tools for designing and demultiplexing of GBS experiments.     

Phylogenetic visualization of the spread of H7 influenza A viruses

Viruses of influenza A subtype H7 can be highly pathogenic and periodically infect humans. For example, there have been numerous outbreaks of H7 in the Americas and Europe since 1996. More recently, a reassortant H7N9 has emerged among humans and birds during 2013–2014 in China, Taiwan and Hong Kong. This H7N9 genome consists of genetic segments that assort with H7 and H9 viruses previously circulating in chickens and wild birds in China and ducks in Korea. Epidemic risk modellers have used agricultural, climatic and demographic data to predict that the virus will spread to northern Vietnam via poultry. To shed light on the traffic of H7 viruses in general, we examine genetic segments of influenza that have assorted with many strains of H7 viruses dating back to 1902. We focus on use cases from the United States, Italy and China. We apply a novel metric, betweenness, an associated phylogenetic visualization technique, transmission networks, and compare these with another technique, route mapping. In contrast to traditional views, our results illustrate that segments that assort with H7 viruses are spread frequently between the Americas and Eurasia. In summary, genetic segments that historically assort with H7 influenza viruses have been spread from China to: Australia, Czech Republic, Denmark, Egypt, Germany, Hong Kong, Italy, Japan, Mongolia, the Netherlands, New Zealand, Pakistan, South Africa, South Korea, Spain, Sweden, the UK, the US, and Vietnam.     

Two soluble glycosyltransferases glycosylate less efficiently in vivo than their membrane bound counterparts

Many Golgi glycosyltransferases are type II membrane proteins which are cleaved to produce soluble forms that are released from cells. Cho and Cummings recently reported that a soluble form of alpha1, 3- galactosyltransferase was comparable to its membrane bound counterpart in its ability to galactosylate newly synthesized glycoproteins (Cho,S.K. and Cummings,R.D. (1997) J. Biol. Chem., 272, 13622-13628). To test the generality of their findings, we compared the activities of the full length and soluble forms of two such glycosyltransferases, ss1,4 N-Acetylgalactosaminyltransferase (GM2/GD2/ GA2 synthase; GalNAcT) and beta galactoside alpha2,6 sialyltransferase (alpha2,6-ST; ST6Gal I), for production of their glycoconjugate products in vivo . Unlike the full length form of GalNAcT which produced ganglioside GM2 in transfected cells, soluble GalNAcT did not produce detectable GM2 in vivo even though it possessed in vitro GalNAcT activity comparable to that of full length GalNAcT. When compared with cells expressing full length alpha2,6-ST, cells expressing a soluble form of alpha2,6-ST contained 3-fold higher alpha2,6-ST mRNA levels and secreted 7-fold greater alpha2,6-ST activity as measured in vitro , but in striking contrast contained 2- to 4-fold less of the alpha2,6-linked sialic acid moiety in cellular glycoproteins in vivo . In summary these results suggest that unlike alpha1,3-galactosyltransferase the soluble forms of these two glycosyltransferases are less efficient at glycosylation of membrane proteins and lipids in vivo than their membrane bound counterparts.      

Mutant p53R175H upregulates Twist1 expression and promotes epithelial–mesenchymal transition in immortalized prostate cells

A mutation within one allele of the p53 tumor suppressor gene can inactivate the remaining wild-type allele in a dominant-negative manner and in some cases can exert an additional oncogenic activity, known as mutant p53 ‘gain of function'' (GOF). To study the role of p53 mutations in prostate cancer and to discriminate between the dominant-negative effect and the GOF activity of mutant p53, we measured, using microarrays, the expression profiles of three immortalized prostate epithelial cultures expressing wild-type, inactivated p53 or mutated p53. Analysis of these gene expression profiles showed that both inactivated p53 and p53^R175H mutant expression resulted in the upregulation of cell cycle progression genes. A second group, which was upregulated exclusively by mutant p53^R175H, was predominantly enriched in developmental genes. This group of genes included the Twist1, a regulator of metastasis and epithelial–mesenchymal transition (EMT). Twist1 levels were also elevated in metastatic prostate cancer-derived cell line DU145, in immortalized lung fibroblasts and in a subset of lung cancer samples, all in a mutant p53-dependent manner. p53^R175H mutant bearing immortalized epithelial cells showed typical features of EMT, such as higher expression of mesenchymal markers, lower expression of epithelial markers and enhanced invasive properties in vitro. The mechanism by which p53^R175H mutant induces Twist1 expression involves alleviation of the epigenetic repression. Our data suggest that Twist1 expression might be upregulated following p53 mutation in cancer cells.     

Insertion of leader peptidase into the thylakoid membrane during synthesis in a chloroplast translation system

The mechanisms of targeting and insertion of chloroplast-encoded thylakoid membrane proteins are poorly understood. In this study, we have used a translation system isolated from chloroplasts to begin to investigate these mechanisms. The bacterial membrane protein leader peptidase (Lep) was used as a model protein because its targeting and insertion mechanisms are well understood for Escherichia coli and for the endoplasmic reticulum. Lep could thus provide insight into the functional homologies between the different membrane systems. Lep was efficiently expressed in the chloroplast translation system, and the protein could be inserted into thylakoid membranes with the same topology as in E. coli cytoplasmic membranes, following the positive-inside rule. Insertion of Lep into the thylakoid membrane was stimulated by the trans-thylakoid proton gradient and was strongly inhibited by azide, suggesting a requirement for SecA activity. Insertion most likely occurred in a cotranslational manner, because insertion could only be observed if thylakoid membranes were present during translation reactions but not when thylakoid membranes were added after translation reactions were terminated. To halt the elongation process at different stages, we translated truncated Lep mRNAs without a stop codon, resulting in the formation of stable ribosome nascent chain complexes. These complexes showed a strong, salt-resistant affinity for the thylakoid membrane, implying a functional interaction of the ribosome with the membrane and supporting a cotranslational insertion mechanism for Lep. Our study supports a functional homology for the insertion of Lep into the thylakoid membrane and the E. coli cytoplasmic membrane.     

Calibration and uncertainties in personal exposure measurements of radiofrequency electromagnetic fields

In the past 5 years radiofrequency personal exposure meters have been used to characterize the exposure during daily activities. We found from calibration tests for the 12 frequency bands of the EME Spy 121 exposimeter in a Gigahertz Transverse Electromagnetic cell and an Open Area Test Site, that these measurements tend to underestimate the actual exposure. Therefore, a maximum frequency‐dependent correction factor of 1.1–1.6 should be applied to the electric field. This correction factor consists of three multipliers correcting for calibration, elevation arrival angle, and influence of the body. The calibration correction factor should be determined per exposimeter, as the maximum range of response between exposimeters in a frequency band is 2.4 dB. Since the range of response for different elevation angles could reach 10.2 dB, a strict protocol for wearing the exposimeter during fieldwork should be followed to be able to compare and combine measurements made by different persons in the same microenvironments. Because the influence of the body depends on the azimuth angle of arrival, it may lead to an over‐ or underestimation. Thus, the body correction factor is an average over the angles and should only be applied in activities involving movement through the full 360° range of random angles of arrival. Bioelectromagnetics. Bioelectromagnetics 32:652–663, 2011. © 2011 Wiley Periodicals, Inc.     

MICU1 Controls Both the Threshold and Cooperative Activation of the Mitochondrial Ca2+ Uniporter

1. Download : Download high-res image (89KB)
2. Download : Download full-size image

     

Origin of the stability conferred upon collagen by fluorination

According to a prevailing theory, (2S,4R)-4-hydroxyproline (Hyp) residues stabilize the collagen triple helix via a stereoelectronic effect that preorganizes appropriate backbone torsion angles for triple-helix formation. This theory is consistent with the marked stability that results from replacing the hydroxyl group with the more electron-withdrawing fluoro group, as in (2S,4R)-4-fluoroproline (Flp). Nonetheless, the hyperstability of triple helices containing Flp has been attributed by others to the hydrophobic effect rather than a stereoelectronic effect. We tested this hypothesis by replacing Hyp with 4,4-difluoroproline (Dfp) in collagen-related peptides. Dfp retains the hydrophobicity of Flp, but lacks the ability of Flp to preorganize backbone torsion angles. Unlike Flp, Dfp does not endow triple helices with elevated stability, indicating that the hyperstability conferred by Flp is not due to the hydrophobic effect.     

Development of the ACTH and corticosterone response to acute hypoxia in the neonatal rat

Acute episodes of severe hypoxia are among the most common stressors in neonates. An understanding of the development of the physiological response to acute hypoxia will help improve clinical interventions. The present study measured ACTH and corticosterone responses to acute, severe hypoxia (8% inspired O(2) for 4 h) in neonatal rats at postnatal days (PD) 2, 5, and 8. Expression of specific hypothalamic, anterior pituitary, and adrenocortical mRNAs was assessed by real-time PCR, and expression of specific proteins in isolated adrenal mitochondria from adrenal zona fascisulata/reticularis was assessed by immunoblot analyses. Oxygen saturation, heart rate, and body temperature were also measured. Exposure to 8% O(2) for as little as 1 h elicited an increase in plasma corticosterone in all age groups studied, with PD2 pups showing the greatest response ( approximately 3 times greater than PD8 pups). Interestingly, the ACTH response to hypoxia was absent in PD2 pups, while plasma ACTH nearly tripled in PD8 pups. Analysis of adrenal mRNA expression revealed a hypoxia-induced increase in Ldlr mRNA at PD2, while both Ldlr and Star mRNA were increased at PD8. Acute hypoxia decreased arterial O(2) saturation (SPo(2)) to approximately 80% and also decreased body temperature by 5-6 degrees C. The hypoxic thermal response may contribute to the ACTH and corticosterone response to decreases in oxygen. The present data describe a developmentally regulated, differential corticosterone response to acute hypoxia, shifting from ACTH independence in early life (PD2) to ACTH dependence less than 1 wk later (PD8).