LLGL2 rescues nutrient stress by promoting leucine uptake in ER+ breast cancer.

Drosophila Lgl and its mammalian homologues, LLGL1 and LLGL2, are scaffolding proteins that regulate the establishment of apical-basal polarity in epithelial cells1,2. Whereas Lgl functions as a tumour suppressor in Drosophila1, the roles of mammalian LLGL1 and LLGL2 in cancer are unclear. The majority (about 75%) of breast cancers express oestrogen receptors (ERs)3, and patients with these tumours receive endocrine treatment4. However, the development of resistance to endocrine therapy and metastatic progression are leading causes of death for patients with ER+ disease4. Here we report that, unlike LLGL1, LLGL2 is overexpressed in ER+ breast cancer and promotes cell proliferation under nutrient stress. LLGL2 regulates cell surface levels of a leucine transporter, SLC7A5, by forming a trimeric complex with SLC7A5 and a regulator of membrane fusion, YKT6, to promote leucine uptake and cell proliferation. The oestrogen receptor targets LLGL2 expression. Resistance to endocrine treatment in breast cancer cells was associated with SLC7A5- and LLGL2-dependent adaption to nutrient stress. SLC7A5 was necessary and sufficient to confer resistance to tamoxifen treatment, identifying SLC7A5 as a potential therapeutic target for overcoming resistance to endocrine treatments in breast cancer. Thus, LLGL2 functions as a promoter of tumour growth and not as a tumour suppressor in ER+ breast cancer. Beyond breast cancer, adaptation to nutrient stress is critically important5, and our findings identify an unexpected role for LLGL2 in this process.

Estrogen receptor signaling is reprogrammed during breast tumorigenesis.

Limited knowledge of the changes in estrogen receptor (ER) signaling during the transformation of the normal mammary gland to breast cancer hinders the development of effective prevention and treatment strategies. Differences in estrogen signaling between normal human primary breast epithelial cells and primary breast tumors obtained immediately following surgical excision were explored. Transcriptional profiling of normal ER+ mature luminal mammary epithelial cells and ER+ breast tumors revealed significant difference in the response to estrogen stimulation. Consistent with these differences in gene expression, the normal and tumor ER cistromes were distinct and sufficient to segregate normal breast tissues from breast tumors. The selective enrichment of the DNA binding motif GRHL2 in the breast cancer-specific ER cistrome suggests that it may play a role in the differential function of ER in breast cancer. Depletion of GRHL2 resulted in altered ER binding and differential transcriptional responses to estrogen stimulation. Furthermore, GRHL2 was demonstrated to be essential for estrogen-stimulated proliferation of ER+ breast cancer cells. DLC1 was also identified as an estrogen-induced tumor suppressor in the normal mammary gland with decreased expression in breast cancer. In clinical cohorts, loss of DLC1 and gain of GRHL2 expression are associated with ER+ breast cancer and are independently predictive for worse survival. This study suggests that normal ER signaling is lost and tumor-specific ER signaling is gained during breast tumorigenesis. Unraveling these changes in ER signaling during breast cancer progression should aid the development of more effective prevention strategies and targeted therapeutics.

ChiLin: a comprehensive ChIP-seq and DNase-seq quality control and analysis pipeline.

BACKGROUND: Transcription factor binding, histone modification, and chromatin accessibility studies are important approaches to understanding the biology of gene regulation. ChIP-seq and DNase-seq have become the standard techniques for studying protein-DNA interactions and chromatin accessibility respectively, and comprehensive quality control (QC) and analysis tools are critical to extracting the most value from these assay types. Although many analysis and QC tools have been reported, few combine ChIP-seq and DNase-seq data analysis and quality control in a unified framework with a comprehensive and unbiased reference of data quality metrics. RESULTS: ChiLin is a computational pipeline that automates the quality control and data analyses of ChIP-seq and DNase-seq data. It is developed using a flexible and modular software framework that can be easily extended and modified. ChiLin is ideal for batch processing of many datasets and is well suited for large collaborative projects involving ChIP-seq and DNase-seq from different designs. ChiLin generates comprehensive quality control reports that include comparisons with historical data derived from over 23,677 public ChIP-seq and DNase-seq samples (11,265 datasets) from eight literature-based classified categories. To the best of our knowledge, this atlas represents the most comprehensive ChIP-seq and DNase-seq related quality metric resource currently available. These historical metrics provide useful heuristic quality references for experiment across all commonly used assay types. Using representative datasets, we demonstrate the versatility of the pipeline by applying it to different assay types of ChIP-seq data. The pipeline software is available open source at https://github.com/cfce/chilin . CONCLUSION: ChiLin is a scalable and powerful tool to process large batches of ChIP-seq and DNase-seq datasets. The analysis output and quality metrics have been structured into user-friendly directories and reports. We have successfully compiled 23,677 profiles into a comprehensive quality atlas with fine classification for users.

In vivo emergence of a novel mutant L159F/L320F in the NS5B polymerase confers low-level resistance to the HCV polymerase inhibitors mericitabine and sofosbuvir.

BACKGROUND: Resistance to mericitabine (prodrug of HCV NS5B polymerase inhibitor PSI-6130) is rare and conferred by the NS5B S282T mutation. METHODS: Serum HCV RNA from patients who experienced viral breakthrough, partial response, or nonresponse in 2 clinical trials in which patients received mericitabine plus peginterferon alfa-2a (40KD)/ribavirin were analyzed by population and clonal sequence analysis as well as phenotypic assay for assessment of in vivo mericitabine resistance. RESULTS: Among 405 patients treated with mericitabine plus peginterferon alfa-2a/ribavirin in PROPEL and JUMP-C, virologic breakthrough or nonresponse were not observed; 12 patients experienced a partial response. The NS5B S282T resistance mutation was not observed in any patient. A number of treatment-associated NS5B changes were observed and characterized. A novel double mutant (L159F/L320F) with impaired replication capacity was detected in one HCV genotype 1b-infected patient. Introduction of double mutant L159F/L320F into genotype 1a (H77) and 1b (Con-1) replicons, respectively, increased the EC50 for mericitabine by 3.1- and 5.5-fold and the EC90 by 3.1- and 8.9-fold. The double mutant also decreased susceptibility to sofosbuvir (GS-7977) and GS-938 but not setrobuvir, relative to wild-type. CONCLUSIONS: A novel and replication-deficient double mutation (L159F/L320F) confers low-level resistance to mericitabine and cross-resistance to both sofosbuvir and GS-938. CLINICAL TRIALS REGISTRATION: NCT00869661, NCT01057667.

Identification of the NS5B S282T resistant variant and two novel amino acid substitutions that affect replication capacity in hepatitis C virus-infected patients treated with mericitabine and danoprevir.

Baseline and posttreatment samples from hepatitis C virus (HCV) genotype (GT) 1-infected patients who received a combination of danoprevir and mericitabine from a phase II clinical study (INFORM-SVR) were analyzed. In addition to resistance monitoring, sequencing and phenotypic assays were combined with statistical analysis to identify potential novel amino acid substitutions associated with treatment outcome. The NS5B S282T substitution associated with mericitabine resistance was identified in 2/30 viral breakthrough patients and was replaced by wild-type viruses after cessation of drug treatment (during follow-up). The NS3 R155K substitution associated with danoprevir resistance was also observed in these 2 patients. All 69 GT 1a-infected patients who experienced viral breakthrough on treatment or relapsed during follow-up (relapsers) developed NS3 R155K. Among GT 1b-infected patients, substitutions at the danoprevir resistance locus NS3 D168 were observed in 15/20 subjects, whereas substitutions at the danoprevir resistance locus NS3 R155 were observed in 5/20 subjects. Interestingly, the baseline polymorphism NS5B Q47H was more prevalent in GT 1a-infected patients who achieved a sustained virologic response at follow-up week 24 (SVR24) than in non-SVR24 patients (2/13 versus 0/72), and a postbaseline NS3 S122G substitution was more prevalent in GT 1a-infected patients with viral breakthrough than in relapsers (4/22 versus 0/47). Neither substitution conferred resistance to danoprevir or mericitabine, but the substitutions reduced (NS5B Q47H) or improved (NS3 S122G) replication capacity by 2- to 4-fold. The NS5B S282T mericitabine-resistant variant was rare and did not persist once drug was discontinued. NS5B Q47H and NS3 S122G are two newly identified substitutions that affected replication capacity and were enriched in distinct treatment response groups. (This study has been registered at ClinicalTrials.gov under registration no. NCT01278134.).

The genome sequence of the leaf-cutter ant Atta cephalotes reveals insights into its obligate symbiotic lifestyle.

Leaf-cutter ants are one of the most important herbivorous insects in the Neotropics, harvesting vast quantities of fresh leaf material. The ants use leaves to cultivate a fungus that serves as the colony's primary food source. This obligate ant-fungus mutualism is one of the few occurrences of farming by non-humans and likely facilitated the formation of their massive colonies. Mature leaf-cutter ant colonies contain millions of workers ranging in size from small garden tenders to large soldiers, resulting in one of the most complex polymorphic caste systems within ants. To begin uncovering the genomic underpinnings of this system, we sequenced the genome of Atta cephalotes using 454 pyrosequencing. One prediction from this ant's lifestyle is that it has undergone genetic modifications that reflect its obligate dependence on the fungus for nutrients. Analysis of this genome sequence is consistent with this hypothesis, as we find evidence for reductions in genes related to nutrient acquisition. These include extensive reductions in serine proteases (which are likely unnecessary because proteolysis is not a primary mechanism used to process nutrients obtained from the fungus), a loss of genes involved in arginine biosynthesis (suggesting that this amino acid is obtained from the fungus), and the absence of a hexamerin (which sequesters amino acids during larval development in other insects). Following recent reports of genome sequences from other insects that engage in symbioses with beneficial microbes, the A. cephalotes genome provides new insights into the symbiotic lifestyle of this ant and advances our understanding of host-microbe symbioses.

Complete mitochondrial genome phylogeographic analysis of killer whales (Orcinus orca) indicates multiple species.

Killer whales (Orcinus orca) currently comprise a single, cosmopolitan species with a diverse diet. However, studies over the last 30 yr have revealed populations of sympatric "ecotypes" with discrete prey preferences, morphology, and behaviors. Although these ecotypes avoid social interactions and are not known to interbreed, genetic studies to date have found extremely low levels of diversity in the mitochondrial control region, and few clear phylogeographic patterns worldwide. This low level of diversity is likely due to low mitochondrial mutation rates that are common to cetaceans. Using killer whales as a case study, we have developed a method to readily sequence, assemble, and analyze complete mitochondrial genomes from large numbers of samples to more accurately assess phylogeography and estimate divergence times. This represents an important tool for wildlife management, not only for killer whales but for many marine taxa. We used high-throughput sequencing to survey whole mitochondrial genome variation of 139 samples from the North Pacific, North Atlantic, and southern oceans. Phylogenetic analysis indicated that each of the known ecotypes represents a strongly supported clade with divergence times ranging from approximately 150,000 to 700,000 yr ago. We recommend that three named ecotypes be elevated to full species, and that the remaining types be recognized as subspecies pending additional data. Establishing appropriate taxonomic designations will greatly aid in understanding the ecological impacts and conservation needs of these important marine predators. We predict that phylogeographic mitogenomics will become an important tool for improved statistical phylogeography and more precise estimates of divergence times.

An insect herbivore microbiome with high plant biomass-degrading capacity.

Herbivores can gain indirect access to recalcitrant carbon present in plant cell walls through symbiotic associations with lignocellulolytic microbes. A paradigmatic example is the leaf-cutter ant (Tribe: Attini), which uses fresh leaves to cultivate a fungus for food in specialized gardens. Using a combination of sugar composition analyses, metagenomics, and whole-genome sequencing, we reveal that the fungus garden microbiome of leaf-cutter ants is composed of a diverse community of bacteria with high plant biomass-degrading capacity. Comparison of this microbiome's predicted carbohydrate-degrading enzyme profile with other metagenomes shows closest similarity to the bovine rumen, indicating evolutionary convergence of plant biomass degrading potential between two important herbivorous animals. Genomic and physiological characterization of two dominant bacteria in the fungus garden microbiome provides evidence of their capacity to degrade cellulose. Given the recent interest in cellulosic biofuels, understanding how large-scale and rapid plant biomass degradation occurs in a highly evolved insect herbivore is of particular relevance for bioenergy.

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes).

BACKGROUND: Two strains of the silver fox (Vulpes vulpes), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed. RESULTS: cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome. CONCLUSIONS: Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.

Prostate cancer reactivates developmental epigenomic programs during metastatic progression.

Epigenetic processes govern prostate cancer (PCa) biology, as evidenced by the dependency of PCa cells on the androgen receptor (AR), a prostate master transcription factor. We generated 268 epigenomic datasets spanning two state transitions-from normal prostate epithelium to localized PCa to metastases-in specimens derived from human tissue. We discovered that reprogrammed AR sites in metastatic PCa are not created de novo; rather, they are prepopulated by the transcription factors FOXA1 and HOXB13 in normal prostate epithelium. Reprogrammed regulatory elements commissioned in metastatic disease hijack latent developmental programs, accessing sites that are implicated in prostate organogenesis. Analysis of reactivated regulatory elements enabled the identification and functional validation of previously unknown metastasis-specific enhancers at HOXB13, FOXA1 and NKX3-1. Finally, we observed that prostate lineage-specific regulatory elements were strongly associated with PCa risk heritability and somatic mutation density. Examining prostate biology through an epigenomic lens is fundamental for understanding the mechanisms underlying tumor progression.

Super-Enhancer-Associated LncRNA UCA1 Interacts Directly with AMOT to Activate YAP Target Genes in Epithelial Ovarian Cancer.

Long noncoding RNAs (lncRNAs) have emerged as critical regulators of tumorigenesis, and yet their mechanistic roles remain challenging to characterize. Here, we integrate functional proteomics with lncRNA-interactome profiling to characterize Urothelial Cancer Associated 1 (UCA1), a candidate driver of ovarian cancer development. Reverse phase protein array (RPPA) analysis indicates that UCA1 activates transcription coactivator YAP and its target genes. In vivo RNA antisense purification (iRAP) of UCA1 interacting proteins identified angiomotin (AMOT), a known YAP regulator, as a direct binding partner. Loss-of-function experiments show that AMOT mediates YAP activation by UCA1, as UCA1 enhances the AMOT-YAP interaction to promote YAP dephosphorylation and nuclear translocation. Together, we characterize UCA1 as a lncRNA regulator of Hippo-YAP signaling and highlight the UCA1-AMOT-YAP signaling axis in ovarian cancer development.

ARv7 Represses Tumor-Suppressor Genes in Castration-Resistant Prostate Cancer.

Androgen deprivation therapy for prostate cancer (PCa) benefits patients with early disease, but becomes ineffective as PCa progresses to a castration-resistant state (CRPC). Initially CRPC remains dependent on androgen receptor (AR) signaling, often through increased expression of full-length AR (ARfl) or expression of dominantly active splice variants such as ARv7. We show in ARv7-dependent CRPC models that ARv7 binds together with ARfl to repress transcription of a set of growth-suppressive genes. Expression of the ARv7-repressed targets and ARv7 protein expression are negatively correlated and predicts for outcome in PCa patients. Our results provide insights into the role of ARv7 in CRPC and define a set of potential biomarkers for tumors dependent on ARv7.

High sensitivity to transmission of attenuation in the metal activating effect of platinum on sequential H-D exchange at the four C8 sites of inosines in tetrakis(inosine)platinum(II) chloride.

Hydrogen-deuterium exchange of the carbon-bound C(8)-H protons of the inosine residues in tetrakis(inosine)platinum(ii) chloride, S, with Pt binding at N(7), was studied in aqueous buffer solutions at 60 degrees C by (1)H NMR spectroscopy. The kinetics at all four C(8) sites as a function of pD of the D(2)O/OD(-) medium was measured through the disappearance of the C(8)-H signal, which yielded the pseudo first-order rate constant for exchange, k(obs). Plots of k(obs)versus [OD(-)] showed curvature reminiscent of saturation type kinetics and indicative of competitive deprotonation of N(1)-H sites. In contrast, the analogous N(1)-methylated cis-bis(1-methylinosine)diammineplatinum(ii) chloride leads to a linear k(obs)versus [OD(-)] plot. The potentiometrically determined macroscopic composite N(1)-H ionization constant was further dissected into the successive microscopic N(1)-H acidity constants of the four inosine residues of the complex S. The k(obs) values were also deconvoluted into individual rate constants k(ex) (i.e.k(0), k(1), k(2), k(3) for exchange of the successively deprotonated inosine moieties, S, S(1), S(2), S(3), it being assumed that S(4) where all four inosine ligands are deprotonated at N(1) is unreactive ("immunized") to exchange. The k(ex) values show a progressive attenuation in Pt activation of H-D exchange along the series, k(0), k(1), k(2), k(3). The k(ex) data thus generated, together with the deconvoluted individual pK(a) values allow the construction of the plot, log k(ex) [C(8)-H] vs. pK(a) [N(H)-1]. Remarkably, this plot exhibits good linearity (R(2) = 0.99), which accords this as a linear free energy relationship (LFER). The large negative slope value (-2.3) of this LFER reflects the high sensitivity of transmission of electron density from the ionized N(1) via Pt and/or through space to the remaining C(8)-H sites. This is to our knowledge the first instance in which a LFER is generated through modulation of a structure in a single molecule. One can anticipate that this approach may lead to: (1) predicting N-H acidity; (2) C-H H-D exchange susceptibility in a range of metal-biomolecule complexes; (3) their carbon acidity.

Mechanical unfolding intermediates observed by single-molecule force spectroscopy in a fibronectin type III module.

Domain 10 of type III fibronectin (10FNIII) is known to play a pivotal role in the mechanical interactions between cell surface integrins and the extracellular matrix. Recent molecular dynamics simulations have predicted that 10FNIII, when exposed to a stretching force, unfolds along two pathways, each with a distinct, mechanically stable intermediate. Here, we use single-molecule force spectroscopy combined with protein engineering to test these predictions by probing the mechanical unfolding pathway of 10FNIII. Stretching single polyproteins containing the 10FNIII module resulted in sawtooth patterns where 10FNIII was seen unfolding in two consecutive steps. The native state unfolded at 100(+/-20) pN, elongating (10)FNIII by 12(+/-2) nm and reaching a clearly marked intermediate that unfolded at 50(+/-20) pN. Unfolding of the intermediate completed the elongation of the molecule by extending another 19(+/-2) nm. Site-directed mutagenesis of residues in the A and B beta-strands (E9P and L19P) resulted in sawtooth patterns with all-or-none unfolding events that elongated the molecule by 19(+/-2) nm. In contrast, mutating residues in the G beta-strand gave results that were dependent on amino acid position. The mutation I88P in the middle of the G beta-strand resulted in native like unfolding sawtooth patterns showing an intact intermediate state. The mutation Y92P, which is near the end of G beta-strand, produced sawtooth patterns with all-or-none unfolding events that lengthened the molecule by 17(+/-2) nm. These results are consistent with the view that 10FNIII can unfold in two different ways. Along one pathway, the detachment of the A and B beta-strands from the body of the folded module constitute the first unfolding event, followed by the unfolding of the remaining beta-sandwich structure. Along the second pathway, the detachment of the G beta-strands is involved in the first unfolding event. These results are in excellent agreement with the sequence of events predicted by molecular dynamics simulations of the 10FNIII module.

Chromatin immunoprecipitation from fixed clinical tissues reveals tumor-specific enhancer profiles.

Extensive cross-linking introduced during routine tissue fixation of clinical pathology specimens severely hampers chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) analysis from archived tissue samples. This limits the ability to study the epigenomes of valuable, clinically annotated tissue resources. Here we describe fixed-tissue chromatin immunoprecipitation sequencing (FiT-seq), a method that enables reliable extraction of soluble chromatin from formalin-fixed paraffin-embedded (FFPE) tissue samples for accurate detection of histone marks. We demonstrate that FiT-seq data from FFPE specimens are concordant with ChIP-seq data from fresh-frozen samples of the same tumors. By using multiple histone marks, we generate chromatin-state maps and identify cis-regulatory elements in clinical samples from various tumor types that can readily allow us to distinguish between cancers by the tissue of origin. Tumor-specific enhancers and superenhancers that are elucidated by FiT-seq analysis correlate with known oncogenic drivers in different tissues and can assist in the understanding of how chromatin states affect gene regulation.

Oncogenic Deregulation of EZH2 as an Opportunity for Targeted Therapy in Lung Cancer.

UNLABELLED: As a master regulator of chromatin function, the lysine methyltransferase EZH2 orchestrates transcriptional silencing of developmental gene networks. Overexpression of EZH2 is commonly observed in human epithelial cancers, such as non-small cell lung carcinoma (NSCLC), yet definitive demonstration of malignant transformation by deregulated EZH2 remains elusive. Here, we demonstrate the causal role of EZH2 overexpression in NSCLC with new genetically engineered mouse models of lung adenocarcinoma. Deregulated EZH2 silences normal developmental pathways, leading to epigenetic transformation independent of canonical growth factor pathway activation. As such, tumors feature a transcriptional program distinct from KRAS- and EGFR-mutant mouse lung cancers, but shared with human lung adenocarcinomas exhibiting high EZH2 expression. To target EZH2-dependent cancers, we developed a potent open-source EZH2 inhibitor, JQEZ5, that promoted the regression of EZH2-driven tumors in vivo, confirming oncogenic addiction to EZH2 in established tumors and providing the rationale for epigenetic therapy in a subset of lung cancer. SIGNIFICANCE: EZH2 overexpression induces murine lung cancers that are similar to human NSCLC with high EZH2 expression and low levels of phosphorylated AKT and ERK, implicating biomarkers for EZH2 inhibitor sensitivity. Our EZH2 inhibitor, JQEZ5, promotes regression of these tumors, revealing a potential role for anti-EZH2 therapy in lung cancer. Cancer Discov; 6(9); 1006-21. ©2016 AACR.See related commentary by Frankel et al., p. 949This article is highlighted in the In This Issue feature, p. 932.

Somatic Cell Fusions Reveal Extensive Heterogeneity in Basal-like Breast Cancer.

Basal-like and luminal breast tumors have distinct clinical behavior and molecular profiles, yet the underlying mechanisms are poorly defined. To interrogate processes that determine these distinct phenotypes and their inheritance pattern, we generated somatic cell fusions and performed integrated genetic and epigenetic (DNA methylation and chromatin) profiling. We found that the basal-like trait is generally dominant and is largely defined by epigenetic repression of luminal transcription factors. Definition of super-enhancers highlighted a core program common in luminal cells but a high degree of heterogeneity in basal-like breast cancers that correlates with clinical outcome. We also found that protein extracts of basal-like cells are sufficient to induce a luminal-to-basal phenotypic switch, implying a trigger of basal-like autoregulatory circuits. We determined that KDM6A might be required for luminal-basal fusions, and we identified EN1, TBX18, and TCF4 as candidate transcriptional regulators of the luminal-to-basal switch. Our findings highlight the remarkable epigenetic plasticity of breast cancer cells.

Similar prevalence of low-abundance drug-resistant variants in treatment-naive patients with genotype 1a and 1b hepatitis C virus infections as determined by ultradeep pyrosequencing.

BACKGROUND AND OBJECTIVES: Hepatitis C virus (HCV) variants that confer resistance to direct-acting-antiviral agents (DAA) have been detected by standard sequencing technology in genotype (G) 1 viruses from DAA-naive patients. It has recently been shown that virological response rates are higher and breakthrough rates are lower in G1b infected patients than in G1a infected patients treated with certain classes of HCV DAAs. It is not known whether this corresponds to a difference in the composition of G1a and G1b HCV quasispecies in regards to the proportion of naturally occurring DAA-resistant variants before treatment. METHODS: We used ultradeep pyrosequencing to determine the prevalence of low-abundance (<25% of the sequence reads) DAA-resistant variants in 191 NS3 and 116 NS5B isolates from 208 DAA-naive G1-infected patients. RESULTS: A total of 3.5 million high-quality reads of ≥ 200 nucleotides were generated. The median coverage depth was 4150x and 4470x per NS3 and NS5B amplicon, respectively. Both G1a and G1b populations showed Shannon entropy distributions, with no difference between G1a and G1b in NS3 or NS5B region at the nucleotide level. A higher number of substitutions that confer resistance to protease inhibitors were observed in G1a isolates (mainly at amino acid 80 of the NS3 region). The prevalence of amino acid substitutions that confer resistance to NS5B non-nucleoside inhibitors was similar in G1a and G1b isolates. The NS5B S282T variant, which confers resistance to the polymerase inhibitors mericitabine and sofosbuvir, was not detected in any sample. CONCLUSION: The quasispecies genetic diversity and prevalence of DAA-resistant variants was similar in G1a and G1b isolates and in both NS3 and NS5B regions, suggesting that this is not a determinant for the higher level of DAA resistance observed across G1a HCV infected patients upon treatment.

Transcriptome profiling of the intoxication response of Tenebrio molitor larvae to Bacillus thuringiensis Cry3Aa protoxin.

Bacillus thuringiensis (Bt) crystal (Cry) proteins are effective against a select number of insect pests, but improvements are needed to increase efficacy and decrease time to mortality for coleopteran pests. To gain insight into the Bt intoxication process in Coleoptera, we performed RNA-Seq on cDNA generated from the guts of Tenebrio molitor larvae that consumed either a control diet or a diet containing Cry3Aa protoxin. Approximately 134,090 and 124,287 sequence reads from the control and Cry3Aa-treated groups were assembled into 1,318 and 1,140 contigs, respectively. Enrichment analyses indicated that functions associated with mitochondrial respiration, signalling, maintenance of cell structure, membrane integrity, protein recycling/synthesis, and glycosyl hydrolases were significantly increased in Cry3Aa-treated larvae, whereas functions associated with many metabolic processes were reduced, especially glycolysis, tricarboxylic acid cycle, and fatty acid synthesis. Microarray analysis was used to evaluate temporal changes in gene expression after 6, 12 or 24 h of Cry3Aa exposure. Overall, microarray analysis indicated that transcripts related to allergens, chitin-binding proteins, glycosyl hydrolases, and tubulins were induced, and those related to immunity and metabolism were repressed in Cry3Aa-intoxicated larvae. The 24 h microarray data validated most of the RNA-Seq data. Of the three intoxication intervals, larvae demonstrated more differential expression of transcripts after 12 h exposure to Cry3Aa. Gene expression examined by three different methods in control vs. Cry3Aa-treated larvae at the 24 h time point indicated that transcripts encoding proteins with chitin-binding domain 3 were the most differentially expressed in Cry3Aa-intoxicated larvae. Overall, the data suggest that T. molitor larvae mount a complex response to Cry3Aa during the initial 24 h of intoxication. Data from this study represent the largest genetic sequence dataset for T. molitor to date. Furthermore, the methods in this study are useful for comparative analyses in organisms lacking a sequenced genome.

Contour length and refolding rate of a small protein controlled by engineered disulfide bonds.

The introduction of disulfide bonds into proteins creates additional mechanical barriers and limits the unfolded contour length (i.e., the maximal extension) measured by single-molecule force spectroscopy. Here, we engineer single disulfide bonds into four different locations of the human cardiac titin module (I27) to control the contour length while keeping the distance to the transition state unchanged. This enables the study of several biologically important parameters. First, we are able to precisely determine the end-to-end length of the transition state before unfolding (53 Angstrom), which is longer than the end-to-end length of the protein obtained from NMR spectroscopy (43 Angstrom). Second, the measured contour length per amino acid from five different methods (4.0 +/- 0.2 Angstrom) is longer than the end-to-end length obtained from the crystal structure (3.6 Angstrom). Our measurement of the contour length takes into account all the internal degrees of freedom of the polypeptide chain, whereas crystallography measures the end-to-end length within the "frozen" protein structure. Furthermore, the control of contour length and therefore the number of amino acids unraveled before reaching the disulfide bond (n) facilitates the test of the chain length dependence on the folding time (tau(F)). We find that both a power law scaling tau(F) lambda n(lambda) with lambda = 4.4, and an exponential scaling with n(0.6) fit the data range, in support of different protein-folding scenarios.