htSeqTools: high-throughput sequencing quality control, processing and visualization in R

We provide a Bioconductor package with quality assessment, processing and visualization tools for high-throughput sequencing data, with emphasis in ChIP-seq and RNA-seq studies. It includes detection of outliers and biases, inefficient immuno-precipitation and overamplification artifacts, de novo identification of read-rich genomic regions and visualization of the location and coverage of genomic region lists. AVAILABILITY: www.bioconductor.org.     

The Evolution of the Major Hepatitis C Genotypes Correlates with Clinical Response to Interferon Therapy

Background

Patients chronically infected with hepatitis C virus (HCV) require significantly different durations of therapy and achieve substantially different sustained virologic response rates to interferon-based therapies, depending on the HCV genotype with which they are infected. There currently exists no systematic framework that explains these genotype-specific response rates. Since humans are the only known natural hosts for HCV–a virus that is at least hundreds of years old–one possibility is that over the time frame of this relationship, HCV accumulated adaptive mutations that confer increasing resistance to the human immune system. Given that interferon therapy functions by triggering an immune response, we hypothesized that clinical response rates are a reflection of viral evolutionary adaptations to the immune system.

Methods and Findings

We have performed the first phylogenetic analysis to include all available full-length HCV genomic sequences (n = 345). This resulted in a new cladogram of HCV. This tree establishes for the first time the relative evolutionary ages of the major HCV genotypes. The outcome data from prospective clinical trials that studied interferon and ribavirin therapy was then mapped onto this new tree. This mapping revealed a correlation between genotype-specific responses to therapy and respective genotype age. This correlation allows us to predict that genotypes 5 and 6, for which there currently are no published prospective trials, will likely have intermediate response rates, similar to genotype 3. Ancestral protein sequence reconstruction was also performed, which identified the HCV proteins E2 and NS5A as potential determinants of genotype-specific clinical outcome. Biochemical studies have independently identified these same two proteins as having genotype-specific abilities to inhibit the innate immune factor double-stranded RNA-dependent protein kinase (PKR).

Conclusion

An evolutionary analysis of all available HCV genomes supports the hypothesis that immune selection was a significant driving force in the divergence of the major HCV genotypes and that viral factors that acquired the ability to inhibit the immune response may play a role in determining genotype-specific response rates to interferon therapy.     

Genetic analysis of the requirement for flp-2, tadV, and rcpB in Actinobacillus actinomycetemcomitans biofilm formation

The tad locus of Actinobacillus actinomycetemcomitans encodes a molecular transport system required for tenacious, nonspecific adherence to surfaces and formation of extremely strong biofilms. This locus is dedicated to the biogenesis of Flp pili, which are required for colonization and virulence. We have previously shown that 11 of the 14 tad locus genes are required for adherence and Flp pilus production. Here, we present genetic and phylogenetic analyses of flp-2, tadV, and rcpB genes in biofilm formation. We show that tadV, predicted to encode prepilin peptidase, is required for adherence. In contrast, targeted insertional inactivation of flp-2, a gene closely related to the prepillin gene flp-1, did not abrogate biofilm formation. Expression studies did not detect Flp2-T7 protein under standard laboratory conditions. We present phylogenetic data showing that there is no significant evidence for natural selection in the available flp-2 sequences from A. actinomycetemcomitans, suggesting that flp-2 does not play a significant role in the biology of this organism. Mutants with insertions at the 3' end of rcpB formed biofilms equivalent to wild-type A. actinomycetemcomitans. Surprisingly, 5' end chromosomal insertion mutants in rcpB were obtained only when a wild-type copy of the rcpB gene was provided in trans or when the Tad secretion system was inactivated. Together, our results strongly suggest that A. actinomycetemcomitans rcpB is essential in the context of a functional tad locus. These data show three different phenotypes for the three genes.     

Negative Epistasis between Sickle and Foetal Haemoglobin Suggests a Reduction in Protection against Malaria

BackgroundHaemoglobin variants, Sickle (HbS) and foetal (HbF) have been associated with malaria protection. This study explores epistatic interactions between HbS and HbF on malaria infection.MethodsThe study was conducted between March 2004 and December 2013 within the sickle cell disease (SCD) programme at Muhimbili National Hospital, Tanzania. SCD status was categorized into HbAA, HbAS and HbSS using hemoglobin electrophoresis and High Performance Liquid Chromatography (HPLC). HbF levels were determined by HPLC. Malaria was diagnosed using rapid diagnostic test and/or blood film. Logistic regression and generalized estimating equations models were used to evaluate associations between SCD status, HbF and malaria.Findings2,049 individuals with age range 0-70 years, HbAA 311(15.2%), HbAS 241(11.8%) and HbSS 1,497(73.1%) were analysed. At enrolment, malaria prevalence was significantly higher in HbAA 13.2% compared to HbAS 1.24% and HbSS 1.34% (p<0.001). Mean HbF was lower in those with malaria compared to those without malaria in HbAA (0.43% vs 0.82%) but was the reverse in HbSS (8.10% vs 5.59%). An increase in HbF was associated with a decrease in risk of malaria OR=0.50 (95%CI: 0.28, 0.90; p=0.021) in HbAA, whereas for HbSS the risk of malaria increased OR=2.94 (1.44, 5.98; p=0.003). A similar pattern was seen during multiple visits; HbAA OR=0.52 (0.34, 0.80; p=0.003) vs HbSS OR=2.01 (1.27, 3.23; p=0.003).ConclusionHigher prevalence of malaria in HbAA compared to HbAS and HbSS confirmed the protective effect of HbS. Lower prevalence of malaria in HbAA with high HbF supports a protective effect of HbF. However, in HbSS, the higher prevalence of malaria with high levels of HbF suggests loss of malaria protection. This is the first epidemiological study to suggest a negative epistasis between HbF and HbS on malaria.     

Identification of Novel Type 2 Diabetes Candidate Genes Involved in the Crosstalk between the Mitochondrial and the Insulin Signaling Systems

Type 2 Diabetes (T2D) is a highly prevalent chronic metabolic disease with strong co-morbidity with obesity and cardiovascular diseases. There is growing evidence supporting the notion that a crosstalk between mitochondria and the insulin signaling cascade could be involved in the etiology of T2D and insulin resistance. In this study we investigated the molecular basis of this crosstalk by using systems biology approaches. We combined, filtered, and interrogated different types of functional interaction data, such as direct protein–protein interactions, co-expression analyses, and metabolic and signaling dependencies. As a result, we constructed the mitochondria-insulin (MITIN) network, which highlights 286 genes as candidate functional linkers between these two systems. The results of internal gene expression analysis of three independent experimental models of mitochondria and insulin signaling perturbations further support the connecting roles of these genes. In addition, we further assessed whether these genes are involved in the etiology of T2D using the genome-wide association study meta-analysis from the DIAGRAM consortium, involving 8,130 T2D cases and 38,987 controls. We found modest enrichment of genes associated with T2D amongst our linker genes (p = 0.0549), including three already validated T2D SNPs and 15 additional SNPs, which, when combined, were collectively associated to increased fasting glucose levels according to MAGIC genome wide meta-analysis (p = 8.12×10⁻⁵). This study highlights the potential of combining systems biology, experimental, and genome-wide association data mining for identifying novel genes and related variants that increase vulnerability to complex diseases.     

Identification of NOG as a specific breast cancer bone metastasis-supporting gene

Metastasis requires numerous biological functions that jointly provide tumor cells from a primary site to seed and colonize a distant organ. Some of these activities are selected for in the primary site, whereas others are acquired at the metastatic niche. We provide molecular evidence showing that the BMP inhibitor, NOG, provides metastatic breast cancer cells with the ability to colonize the bone. NOG expression is acquired during the late events of metastasis, once cells have departed from the primary site, because it is not enriched in primary tumors with high risk of bone relapse. On the contrary, breast cancer bone metastatic lesions do select for high levels of NOG expression when compared with metastasis to the lung, liver, and brain. Pivotal to the bone colonization functions is the contribution of NOG to metastatic autonomous and nonautonomous cell functions. Using genetic approaches, we show that when NOG is expressed in human breast cancer cells, it facilitates bone colonization by fostering osteoclast differentiation and bone degradation and also contributes to metastatic lesions reinitiation. These findings reveal how aggressive cancer cell autonomous and nonautonomous functions can be mechanistically coupled to greater bone metastatic potential.     

Nonspecific adherence by Actinobacillus actinomycetemcomitans requires genes widespread in bacteria and archaea

The gram-negative coccobacillus, Actinobacillus actinomycetemcomitans, is the putative agent for localized juvenile periodontitis, a particularly destructive form of periodontal disease in adolescents. This bacterium has also been isolated from a variety of other infections, notably endocarditis. Fresh clinical isolates of A. actinomycetemcomitans form tenacious biofilms, a property likely to be critical for colonization of teeth and other surfaces. Here we report the identification of a locus of seven genes required for nonspecific adherence of A. actinomycetemcomitans to surfaces. The recently developed transposon IS903phikan was used to isolate mutants of the rough clinical isolate CU1000 that are defective in tight adherence to surfaces (Tad(-)). Unlike wild-type cells, Tad(-) mutant cells adhere poorly to surfaces, fail to form large autoaggregates, and lack long, bundled fibrils. Nucleotide sequencing and genetic complementation analysis revealed a 6.7-kb region of the genome with seven adjacent genes (tadABCDEFG) required for tight adherence. The predicted TadA polypeptide is similar to VirB11, an ATPase involved in macromolecular transport. The predicted amino acid sequences of the other Tad polypeptides indicate membrane localization but no obvious functions. We suggest that the tad genes are involved in secretion of factors required for tight adherence of A. actinomycetemcomitans. Remarkably, complete and highly conserved tad gene clusters are present in the genomes of the bubonic plague bacillus Yersinia pestis and the human and animal pathogen Pasteurella multocida. Partial tad loci also occur in strikingly diverse Bacteria and Archaea. Our results show that the tad genes are required for tight adherence of A. actinomycetemcomitans to surfaces and are therefore likely to be essential for colonization and pathogenesis. The occurrence of similar genes in a wide array of microorganisms indicates that they have important functions. We propose that tad-like genes have a significant role in microbial colonization.     

caos software for use in character-based DNA barcoding

The success of character-based DNA barcoding depends on the efficient identification of diagnostic character states from molecular sequences that have been organized hierarchically (e.g. according to phylogenetic methods). Similarly, the reliability of these identified diagnostic character states must be assessed according to their ability to diagnose new sequences. Here, a set of software tools is presented that implement the previously described Characteristic Attribute Organization System for both diagnostic identification and diagnostic-based classification. The software is publicly available from http://sarkarlab.mbl.edu/CAOS.