BRCA1 and RNAi factors promote repair mediated by small RNAs and PALB2-RAD52.

Strong connections exist between R-loops (three-stranded structures harbouring an RNA:DNA hybrid and a displaced single-strand DNA), genome instability and human disease1-5. Indeed, R-loops are favoured in relevant genomic regions as regulators of certain physiological processes through which homeostasis is typically maintained. For example, transcription termination pause sites regulated by R-loops can induce the synthesis of antisense transcripts that enable the formation of local, RNA interference (RNAi)-driven heterochromation6. Pause sites are also protected against endogenous single-stranded DNA breaks by BRCA17. Hypotheses about how DNA repair is enacted at pause sites include a role for RNA, which is emerging as a normal, albeit unexplained, regulator of genome integrity8. Here we report that a species of single-stranded, DNA-damage-associated small RNA (sdRNA) is generated by a BRCA1-RNAi protein complex. sdRNAs promote DNA repair driven by the PALB2-RAD52 complex at transcriptional termination pause sites that form R-loops and are rich in single-stranded DNA breaks. sdRNA repair operates in both quiescent (G0) and proliferating cells. Thus, sdRNA repair can occur in intact tissue and/or stem cells, and may contribute to tumour suppression mediated by BRCA1.

Kinase-independent function of E-type cyclins in liver cancer.

E-type cyclins (cyclins E1 and E2) are components of the core cell cycle machinery and are overexpressed in many human tumor types. E cyclins are thought to drive tumor cell proliferation by activating the cyclin-dependent kinase 2 (CDK2). The cyclin E1 gene represents the site of recurrent integration of the hepatitis B virus in the pathogenesis of hepatocellular carcinoma, and this event is associated with strong up-regulation of cyclin E1 expression. Regardless of the underlying mechanism of tumorigenesis, the majority of liver cancers overexpress E-type cyclins. Here we used conditional cyclin E knockout mice and a liver cancer model to test the requirement for the function of E cyclins in liver tumorigenesis. We show that a ubiquitous, global shutdown of E cyclins did not visibly affect postnatal development or physiology of adult mice. However, an acute ablation of E cyclins halted liver cancer progression. We demonstrated that also human liver cancer cells critically depend on E cyclins for proliferation. In contrast, we found that the function of the cyclin E catalytic partner, CDK2, is dispensable in liver cancer cells. We observed that E cyclins drive proliferation of tumor cells in a CDK2- and kinase-independent mechanism. Our study suggests that compounds which degrade or inhibit cyclin E might represent a highly selective therapeutic strategy for patients with liver cancer, as these compounds would selectively cripple proliferation of tumor cells, while sparing normal tissues.

Loss of Tight Junction Protein Claudin 18 Promotes Progressive Neoplasia Development in Mouse Stomach.

BACKGROUND & AIMS: Loss of claudin 18 (CLDN18), a membrane-spanning tight junction protein, occurs during early stages of development of gastric cancer and associates with shorter survival times of patients. We investigated whether loss of CLDN18 occurs in mice that develop intraepithelial neoplasia with invasive glands due to infection with Helicobacter pylori, and whether loss is sufficient to promote the development of similar lesions in mice with or without H pylori infection. METHODS: We performed immunohistochemical analyses in levels of CLDN18 in archived tissues from B6:129 mice infected with H pylori for 6 to 15 months. We analyzed gastric tissues from B6:129S5-Cldn18tm1Lex/Mmucd mice, in which the CLDN18 gene was disrupted in gastric tissues (CLDN18-knockout mice), or from control mice with a full-length CLDN18 gene (CLDN18+/+; B6:129S5/SvEvBrd) or heterozygous disruption of CLDN18 (CLDN18+/-; B6:129S5/SvEvBrd) that were infected with H pylori SS1 or PMSS1 at 6 weeks of age and tissues collected for analysis at 20 and 30 weeks after infection. Tissues from CLDN18-knockout mice and control mice with full-length CLDN18 gene expression were also analyzed without infection at 7 weeks and 2 years after birth. Tissues from control and CLDN18-knockout mice were analyzed by electron microscopy, stained by conventional methods and analyzed for histopathology, prepared by laser capture microdissection and analyzed by RNAseq, and immunostained for lineage markers, proliferation markers, and stem cell markers and analyzed by super-resolution or conventional confocal microscopy. RESULTS: CLDN18 had a basolateral rather than apical tight junction localization in gastric epithelial cells. B6:129 mice infected with H pylori, which developed intraepithelial neoplasia with invasive glands, had increasing levels of CLDN18 loss over time compared with uninfected mice. In B6:129 mice infected with H pylori compared with uninfected mice, CLDN18 was first lost from most gastric glands followed by disrupted and reduced expression in the gastric neck and in surface cells. Gastric tissues from CLDN18-knockout mice had low levels of inflammation but increased cell proliferation, expressed markers of intestinalized proliferative spasmolytic polypeptide-expressing metaplasia, and had defects in signal transduction pathways including p53 and STAT signaling by 7 weeks after birth compared with full-length CLDN18 gene control mice. By 20 to 30 weeks after birth, gastric tissues from uninfected CLDN18-knockout mice developed intraepithelial neoplasia that invaded the submucosa; by 2 years, gastric tissues contained large and focally dysplastic polypoid tumors with invasive glands that invaded the serosa. CONCLUSIONS: H pylori infection of B6:129 mice reduced the expression of CLDN18 early in gastric cancer progression, similar to previous observations from human gastric tissues. CLDN18 regulates cell lineage differentiation and cellular signaling in mouse stomach; CLDN18-knockout mice develop intraepithelial neoplasia and then large and focally dysplastic polypoid tumors in the absence of H pylori infection.

Clinical and Molecular Analysis of Pathologic Fracture-associated Osteosarcoma: MicroRNA profile Is Different and Correlates with Prognosis.

BACKGROUND: MicroRNAs are small, noncoding RNAs that regulate the expression of posttranslational genes. The presence of some specific microRNAs has been associated with increased risk of both local recurrence and metastasis and worse survival in patients with osteosarcoma. Pathologic fractures in osteosarcoma are considered to be more the manifestation of a neoplasm with a more aggressive biological behavior than the cause itself of worse prognosis. However, this has not been proved at the biological or molecular level. Currently, there has not been a microRNA profiling study of patients who have osteosarcoma with and without pathologic fractures that has described differences in terms of microRNA profiling between these two groups and their correlation with biologic behavior. QUESTIONS/PURPOSES: (1) In patients with osteosarcoma of the extremities, how do the microRNA profiles of those with and without pathologic fractures compare? (2) What relationship do microRNAs have with local recurrence, risk of metastasis, disease-specific survival, and overall survival in osteosarcoma patients with pathologic fractures? METHODS: Between 1994 and 2013, 217 patients were diagnosed and treated at our institution for osteosarcoma of the extremities. Patients were excluded if (1) they underwent oncologic resection of the osteosarcoma at an outside institution (two patients) or (2) they were diagnosed with an extraskeletal osteosarcoma (29 patients) or (3) they had less than 1 year of clinical follow-up and no oncologic outcome (local recurrence, metastasis, or death) (four patients). A total of 182 patients were eligible. Of those, 143 were high-grade osteosarcomas. After evaluation of tumor samples before chemotherapy treatment, a total of 80 consecutive samples were selected for sequencing. Demographic and clinical comparison between the sequenced and non-sequenced patients did not demonstrate any differences, confirming that both groups were comparable. Diagnostic samples from the extremities of 80 patients with high-grade extremity osteosarcomas who had not yet received chemotherapy underwent microRNA sequencing for an ongoing large-scale osteosarcoma genome profiling project at our institution. Six samples were removed after a second look by a musculoskeletal pathologist who verified cellularity and quality of samples to be sequenced, leaving a total of 74 patients. Of these, two samples were removed as they were confirmed to be pelvic tumors in a second check after sequencing. The final study sample was 72 patients (11 patients with pathologic fractures and 61 without). Sequencing data were correlated with fractures and local recurrence, risk of metastasis, disease-specific survival, and overall survival through Kaplan-Meier analyses. RESULTS: Several microRNAs were expressed differently between the two groups. Among the markers with the highest differential expression (edgeR and DESeq algorithms), Hsa-mIR 656-3p, hsa-miR 493-5p, and hsa-miR 381-3p were upregulated in patients with pathologic fractures, whereas hsa-miR 363, hsa-miR 885-5p, and has-miR 20b-5p were downregulated. The highest differential expression fracture and nonfracture-associated microRNA markers also distinguished groups of patients with different metastasis risk, a well as different disease-specific and overall survival. Furthermore, the profile of pathologic fractures demonstrated a higher differential expression for microRNA markers that were previously associated with a higher risk of metastasis and lower survival rates in patients with osteosarcoma. CONCLUSIONS: In patients who have osteosarcoma, the microRNA profiles of those with pathologic fractures are different than of patients without pathologic fractures. The highest differential expression mircroRNA molecules in patients with pathologic fractures predict also higher risk of metastatic disease as well as worse disease-specific survival and overall survival. Furthermore, we found higher differential expression of microRNAs in the pathologic fracture group previously associated with poor prognosis. The higher risk of metastasis and poorer overall survival in patients with pathologic fractures is inherent to tumor aggressive biologic behavior. It is plausible that the fracture itself is not the direct cause of worse prognosis but another manifestation of tumor biologic aggressiveness. Identification of these molecules through liquid biopsies may help to determine which patients may benefit from surgery before fractures occur. The same technology can be applied to identify patterns of response to conventional chemotherapy, assisting in more specific and accurate systemic therapy. LEVEL OF EVIDENCE LEVEL: III, prognostic study.

Proteomic Landscape of Tissue-Specific Cyclin E Functions in Vivo.

E-type cyclins (cyclins E1 and E2) are components of the cell cycle machinery that has been conserved from yeast to humans. The major function of E-type cyclins is to drive cell division. It is unknown whether in addition to their 'core' cell cycle functions, E-type cyclins also perform unique tissue-specific roles. Here, we applied high-throughput mass spectrometric analyses of mouse organs to define the repertoire of cyclin E protein partners in vivo. We found that cyclin E interacts with distinct sets of proteins in different compartments. These cyclin E interactors are highly enriched for phosphorylation targets of cyclin E and its catalytic partner, the cyclin-dependent kinase 2 (Cdk2). Among cyclin E interactors we identified several novel tissue-specific substrates of cyclin E-Cdk2 kinase. In proliferating compartments, cyclin E-Cdk2 phosphorylates Lin proteins within the DREAM complex. In the testes, cyclin E-Cdk2 phosphorylates Mybl1 and Dmrtc2, two meiotic transcription factors that represent key regulators of spermatogenesis. In embryonic and adult brains cyclin E interacts with proteins involved in neurogenesis, while in adult brains also with proteins regulating microtubule-based processes and microtubule cytoskeleton. We also used quantitative proteomics to demonstrate re-wiring of the cyclin E interactome upon ablation of Cdk2. This approach can be used to study how protein interactome changes during development or in any pathological state such as aging or cancer.

A function for cyclin D1 in DNA repair uncovered by protein interactome analyses in human cancers.

Cyclin D1 is a component of the core cell cycle machinery. Abnormally high levels of cyclin D1 are detected in many human cancer types. To elucidate the molecular functions of cyclin D1 in human cancers, we performed a proteomic screen for cyclin D1 protein partners in several types of human tumours. Analyses of cyclin D1 interactors revealed a network of DNA repair proteins, including RAD51, a recombinase that drives the homologous recombination process. We found that cyclin D1 directly binds RAD51, and that cyclin D1-RAD51 interaction is induced by radiation. Like RAD51, cyclin D1 is recruited to DNA damage sites in a BRCA2-dependent fashion. Reduction of cyclin D1 levels in human cancer cells impaired recruitment of RAD51 to damaged DNA, impeded the homologous recombination-mediated DNA repair, and increased sensitivity of cells to radiation in vitro and in vivo. This effect was seen in cancer cells lacking the retinoblastoma protein, which do not require D-cyclins for proliferation. These findings reveal an unexpected function of a core cell cycle protein in DNA repair and suggest that targeting cyclin D1 may be beneficial also in retinoblastoma-negative cancers which are currently thought to be unaffected by cyclin D1 inhibition.

RNA polymerase II promoter-proximal pausing in mammalian long non-coding genes.

Mammalian genomes encode a large number of non-coding RNAs (ncRNAs) that greatly exceed mRNA genes. While the physiological and pathological roles of ncRNAs have been increasingly understood, the mechanisms of regulation of ncRNA expression are less clear. Here, our genomic study has shown that a significant number of long non-coding RNAs (lncRNAs, >1000 nucleotides) harbor RNA polymerase II (Pol II) engaged with the transcriptional start site. A pausing and transcriptional elongation factor for protein-coding genes, tripartite motif-containing 28 (TRIM28) regulates the transcription of a subset of lncRNAs in mammalian cells. In addition, the majority of lncRNAs in human and murine cells regulated by Pol II promoter-proximal pausing appear to function in stimulus-inducible biological pathways. Our findings suggest an important role of Pol II pausing for the transcription of mammalian lncRNA genes.

Analyses of clinicopathological, molecular, and prognostic associations of KRAS codon 61 and codon 146 mutations in colorectal cancer: cohort study and literature review.

BACKGROUND: KRAS mutations in codons 12 and 13 are established predictive biomarkers for anti-EGFR therapy in colorectal cancer. Previous studies suggest that KRAS codon 61 and 146 mutations may also predict resistance to anti-EGFR therapy in colorectal cancer. However, clinicopathological, molecular, and prognostic features of colorectal carcinoma with KRAS codon 61 or 146 mutation remain unclear. METHODS: We utilized a molecular pathological epidemiology database of 1267 colon and rectal cancers in the Nurse's Health Study and the Health Professionals Follow-up Study. We examined KRAS mutations in codons 12, 13, 61 and 146 (assessed by pyrosequencing), in relation to clinicopathological features, and tumor molecular markers, including BRAF and PIK3CA mutations, CpG island methylator phenotype (CIMP), LINE-1 methylation, and microsatellite instability (MSI). Survival analyses were performed in 1067 BRAF-wild-type cancers to avoid confounding by BRAF mutation. Cox proportional hazards models were used to compute mortality hazard ratio, adjusting for potential confounders, including disease stage, PIK3CA mutation, CIMP, LINE-1 hypomethylation, and MSI. RESULTS: KRAS codon 61 mutations were detected in 19 cases (1.5%), and codon 146 mutations in 40 cases (3.2%). Overall KRAS mutation prevalence in colorectal cancers was 40% (=505/1267). Of interest, compared to KRAS-wild-type, overall, KRAS-mutated cancers more frequently exhibited cecal location (24% vs. 12% in KRAS-wild-type; P < 0.0001), CIMP-low (49% vs. 32% in KRAS-wild-type; P < 0.0001), and PIK3CA mutations (24% vs. 11% in KRAS-wild-type; P < 0.0001). These trends were evident irrespective of mutated codon, though statistical power was limited for codon 61 mutants. Neither KRAS codon 61 nor codon 146 mutation was significantly associated with clinical outcome or prognosis in univariate or multivariate analysis [colorectal cancer-specific mortality hazard ratio (HR) = 0.81, 95% confidence interval (CI) = 0.29-2.26 for codon 61 mutation; colorectal cancer-specific mortality HR = 0.86, 95% CI = 0.42-1.78 for codon 146 mutation]. CONCLUSIONS: Tumors with KRAS mutations in codons 61 and 146 account for an appreciable proportion (approximately 5%) of colorectal cancers, and their clinicopathological and molecular features appear generally similar to KRAS codon 12 or 13 mutated cancers. To further assess clinical utility of KRAS codon 61 and 146 testing, large-scale trials are warranted.

CXCR4 pathway associated with family history of melanoma.

PURPOSE: Genetic predisposition plays a major role in the etiology of melanoma, but known genetic markers only account for a limited fraction of family-history-associated melanoma cases. Expression microarrays have offered the opportunity to identify further genomic profiles correlated with family history of melanoma. We aimed to distinguish mRNA expression signatures between melanoma cases with and without a family history of melanoma. METHODS: Based on the Nurses' Health Study, family history was defined as having one or more first-degree family members diagnosed with melanoma. Melanoma diagnosis was confirmed by reviewing pathology reports, and tumor blocks were collected by mail from across the USA. Genomic interrogation was accomplished through evaluating expression profiling of formalin-fixed paraffin-embedded tissues from 78 primary cutaneous invasive melanoma cases, on either a 6K or whole-genome (24K) Illumina gene chip. Gene set enrichment analysis was performed for each batch to determine the differentially enriched pathways and key contributing genes. RESULTS: The CXC chemokine receptor 4 (CXCR4) pathway was consistently up-regulated within cases of familial melanoma in both platforms. Leading edge analysis showed four genes from the CXCR4 pathway, including MAPK1, PLCG1, CRK, and PTK2, were among the core members that contributed to the enrichment of this pathway. There was no association between the enrichment of CXCR4 pathway and NRAS, BRAF mutation, or Breslow thickness of the primary melanoma cases. CONCLUSIONS: We found that the CXCR4 pathway might constitute a novel susceptibility pathway associated with family history of melanoma in first-degree relatives.

Next-generation sequencing and microarray-based interrogation of microRNAs from formalin-fixed, paraffin-embedded tissue: preliminary assessment of cross-platform concordance.

Next-generation sequencing is increasingly employed in biomedical investigations. Strong concordance between microarray and mRNA-seq levels has been reported in high quality specimens but information is lacking on formalin-fixed, paraffin-embedded (FFPE) tissues, and particularly for microRNA (miRNA) analysis. We conducted a preliminary examination of the concordance between miRNA-seq and cDNA-mediated annealing, selection, extension, and ligation (DASL) miRNA assays. Quantitative agreement between platforms is moderate (Spearman correlation 0.514-0.596) and there is discordance of detection calls on a subset of miRNAs. Quantitative PCR (q-RT-PCR) performed for several discordant miRNAs confirmed the presence of most sequences detected by miRNA-seq but not by DASL but also that miRNA-seq did not detect some sequences, which DASL confidently detected. Our results suggest that miRNA-seq is specific, with few false positive calls, but it may not detect certain abundant miRNAs in FFPE tissue. Further work is necessary to fully address these issues that are pertinent for translational research.

Compensatory mutations restore the replication defects caused by cytotoxic T lymphocyte escape mutations in hepatitis C virus polymerase.

While human leukocyte antigen B57 (HLA-B57) is associated with the spontaneous clearance of hepatitis C virus (HCV), the mechanisms behind this control remain unclear. Immunodominant CD8(+) T cell responses against the B57-restricted epitopes comprised of residues 2629 to 2637 of nonstructural protein 5B (NS5B(2629-2637)) (KSKKTPMGF) and E2(541-549) (NTRPPLGNW) were recently shown to be crucial in the control of HCV infection. Here, we investigated whether the selection of deleterious cytotoxic T lymphocyte (CTL) escape mutations in the NS5B KSKKTPMGF epitope might impair viral replication and contribute to the B57-mediated control of HCV. Common CTL escape mutations in this epitope were identified from a cohort of 374 HCV genotype 1a-infected subjects, and their impact on HCV replication assessed using a transient HCV replicon system. We demonstrate that while escape mutations at residue 2633 (position 5) of the epitope had little or no impact on HCV replication in vitro, mutations at residue 2629 (position 1) substantially impaired replication. Notably, the deleterious mutations at position 2629 were tightly linked in vivo to upstream mutations at residue 2626, which functioned to restore the replicative defects imparted by the deleterious escape mutations. These data suggest that the selection of costly escape mutations within the immunodominant NS5B KSKKTPMGF epitope may contribute in part to the control of HCV replication in B57-positive individuals and that persistence of HCV in B57-positive individuals may involve the development of specific secondary compensatory mutations. These findings are reminiscent of the selection of deleterious CTL escape and compensatory mutations by HLA-B57 in HIV-1 infection and, thus, may suggest a common mechanism by which alleles like HLA-B57 mediate protection against these highly variable pathogens.

Protective HLA class I alleles that restrict acute-phase CD8+ T-cell responses are associated with viral escape mutations located in highly conserved regions of human immunodeficiency virus type 1.

The control of human immunodeficiency virus type 1 (HIV-1) associated with particular HLA class I alleles suggests that some CD8(+) T-cell responses may be more effective than others at containing HIV-1. Unfortunately, substantial diversities in the breadth, magnitude, and function of these responses have impaired our ability to identify responses most critical to this control. It has been proposed that CD8 responses targeting conserved regions of the virus may be particularly effective, since the development of cytotoxic T-lymphocyte (CTL) escape mutations in these regions may significantly impair viral replication. To address this hypothesis at the population level, we derived near-full-length viral genomes from 98 chronically infected individuals and identified a total of 76 HLA class I-associated mutations across the genome, reflective of CD8 responses capable of selecting for sequence evolution. The majority of HLA-associated mutations were found in p24 Gag, Pol, and Nef. Reversion of HLA-associated mutations in the absence of the selecting HLA allele was also commonly observed, suggesting an impact of most CTL escape mutations on viral replication. Although no correlations were observed between the number or location of HLA-associated mutations and protective HLA alleles, limiting the analysis to mutations selected by acute-phase immunodominant responses revealed a strong positive correlation between mutations at conserved residues and protective HLA alleles. These data suggest that control of HIV-1 may be associated with acute-phase CD8 responses capable of selecting for viral escape mutations in highly conserved regions of the virus, supporting the inclusion of these regions in the design of an effective vaccine.

Transmission and long-term stability of compensated CD8 escape mutations.

Human immunodeficiency virus effectively evades CD8(+) T-cell responses through the development of CD8 escape mutations. Recent reports documenting reversion of transmitted mutations and the impact of specific escape mutations upon viral replication suggest that complex forces limit the accumulation of CD8 escape mutations at the population level. However, the presence of compensatory mutations capable of alleviating the impact of CD8 escape mutations on replication capacity may enable their persistence in an HLA-mismatched host. Herein, we illustrate the long-term stability of stereotypic escape mutations in the immunodominant HLA-B27-restricted epitope KK10 in p24/Gag following transmission when accompanied by a specific compensatory mutation.

Structural and functional constraints limit options for cytotoxic T-lymphocyte escape in the immunodominant HLA-B27-restricted epitope in human immunodeficiency virus type 1 capsid.

Control of human immunodeficiency virus type 1 (HIV-1) by HLA-B27-positive subjects has been linked to an immunodominant CD8(+) cytotoxic T-lymphocyte (CTL) response targeting the conserved KK10 epitope (KRWIILGLNK(263-272)) in p24/Gag. Viral escape in KK10 typically occurs through development of an R(264)K substitution in conjunction with the upstream compensatory mutation S(173)A, and the difficulty of the virus to escape from the immune response against the KK10 epitope until late in infection has been associated with slower clinical progression. Rare alternative escape mutations at R(264) have been observed, but factors dictating the preferential selection of R(264)K remain unclear. Here we illustrate that while all observed R(264) mutations (K, G, Q, and T) reduced peptide binding to HLA-B27 and impaired viral replication, the replicative defects of the alternative mutants were actually less pronounced than those for R(264)K. Importantly, however, none of these mutants replicated as well as an R(264)K variant containing the compensatory mutation S(173)A. In assessing the combined effects of viral replication and CTL escape using an in vitro coculture assay, we further observed that the compensated R(264)K mutant also displayed the highest replication capacity in the presence of KK10-specific CTLs. Comparisons of codon usage for the respective variants indicated that generation of the R(264)K mutation may also be favored due to a G-to-A bias in nucleotide substitutions during HIV-1 replication. Together, these data suggest that the preference for R(264)K is due primarily to the ability of the S(173)A-compensated virus to replicate better than alternative variants in the presence of CTLs, suggesting that viral fitness is a key contributor for the selection of immune escape variants.

Genetic characterization of human immunodeficiency virus type 1 in elite controllers: lack of gross genetic defects or common amino acid changes.

Despite reports of viral genetic defects in persons who control human immunodeficiency virus type 1 (HIV-1) in the absence of antiviral therapy, the extent to which such defects contribute to the long-term containment of viremia is not known. Most previous studies examining for such defects have involved small numbers of subjects, primarily focused on subjects expressing HLA-B57, or have examined single viral genes, and they have focused on cellular proviral DNA rather than plasma viral RNA sequences. Here, we attempted viral sequencing from 95 HIV-1 elite controllers (EC) who maintained plasma viral loads of <50 RNA copies/ml in the absence of therapy, the majority of whom did not express HLA-B57. HIV-1 gene fragments were obtained from 94% (89/95) of the EC, and plasma viral sequences were obtained from 78% (61/78), the latter indicating the presence of replicating virus in the majority of EC. Of 63 persons for whom nef was sequenced, only three cases of nef deletions were identified, and gross genetic defects were rarely observed in other HIV-1 coding genes. In a codon-by-codon comparison between EC and persons with progressive infection, correcting for HLA bias and coevolving secondary mutations, a significant difference was observed at only three codons in Gag, all three of which represented the historic population consensus amino acid at the time of infection. These results indicate that the spontaneous control of HIV replication is not attributable to shared viral genetic defects or shared viral polymorphisms.

Plasma Circulating Extracellular RNAs in Left Ventricular Remodeling Post-Myocardial Infarction.

Despite substantial declines in mortality following myocardial infarction (MI), subsequent left ventricular remodeling (LVRm) remains a significant long-term complication. Extracellular small non-coding RNAs (exRNAs) have been associated with cardiac inflammation and fibrosis and we hypothesized that they are associated with post-MI LVRm phenotypes. RNA sequencing of exRNAs was performed on plasma samples from patients with "beneficial" (decrease LVESVI ≥ 20%, n = 11) and "adverse" (increase LVESVI ≥ 15%, n = 11) LVRm. Selected differentially expressed exRNAs were validated by RT-qPCR (n = 331) and analyzed for their association with LVRm determined by cardiac MRI. Principal components of exRNAs were associated with LVRm phenotypes post-MI; specifically, LV mass, LV ejection fraction, LV end systolic volume index, and fibrosis. We then investigated the temporal regulation and cellular origin of exRNAs in murine and cell models and found that: 1) plasma and tissue miRNA expression was temporally regulated; 2) the majority of the miRNAs were increased acutely in tissue and at sub-acute or chronic time-points in plasma; 3) miRNA expression was cell-specific; and 4) cardiomyocytes release a subset of the identified miRNAs packaged in exosomes into culture media in response to hypoxia/reoxygenation. In conclusion, we find that plasma exRNAs are temporally regulated and are associated with measures of post-MI LVRm.

Comprehensive multi-center assessment of small RNA-seq methods for quantitative miRNA profiling.

RNA-seq is increasingly used for quantitative profiling of small RNAs (for example, microRNAs, piRNAs and snoRNAs) in diverse sample types, including isolated cells, tissues and cell-free biofluids. The accuracy and reproducibility of the currently used small RNA-seq library preparation methods have not been systematically tested. Here we report results obtained by a consortium of nine labs that independently sequenced reference, 'ground truth' samples of synthetic small RNAs and human plasma-derived RNA. We assessed three commercially available library preparation methods that use adapters of defined sequence and six methods using adapters with degenerate bases. Both protocol- and sequence-specific biases were identified, including biases that reduced the ability of small RNA-seq to accurately measure adenosine-to-inosine editing in microRNAs. We found that these biases were mitigated by library preparation methods that incorporate adapters with degenerate bases. MicroRNA relative quantification between samples using small RNA-seq was accurate and reproducible across laboratories and methods.

PFAAT version 2.0: a tool for editing, annotating, and analyzing multiple sequence alignments.

BACKGROUND: By virtue of their shared ancestry, homologous sequences are similar in their structure and function. Consequently, multiple sequence alignments are routinely used to identify trends that relate to function. This type of analysis is particularly productive when it is combined with structural and phylogenetic analysis. RESULTS: Here we describe the release of PFAAT version 2.0, a tool for editing, analyzing, and annotating multiple sequence alignments. Support for multiple annotations is a key component of this release as it provides a framework for most of the new functionalities. The sequence annotations are accessible from the alignment and tree, where they are typically used to label sequences or hyperlink them to related databases. Sequence annotations can be created manually or extracted automatically from UniProt entries. Once a multiple sequence alignment is populated with sequence annotations, sequences can be easily selected and sorted through a sophisticated search dialog. The selected sequences can be further analyzed using statistical methods that explicitly model relationships between the sequence annotations and residue properties. Residue annotations are accessible from the alignment viewer and are typically used to designate binding sites or properties for a particular residue. Residue annotations are also searchable, and allow one to quickly select alignment columns for further sequence analysis, e.g. computing percent identities. Other features include: novel algorithms to compute sequence conservation, mapping conservation scores to a 3D structure in Jmol, displaying secondary structure elements, and sorting sequences by residue composition. CONCLUSION: PFAAT provides a framework whereby end-users can specify knowledge for a protein family in the form of annotation. The annotations can be combined with sophisticated analysis to test hypothesis that relate to sequence, structure and function.

WebMeV: A Cloud Platform for Analyzing and Visualizing Cancer Genomic Data.

Although large, complex genomic datasets are increasingly easy to generate, and the number of publicly available datasets in cancer and other diseases is rapidly growing, the lack of intuitive, easy-to-use analysis tools has remained a barrier to the effective use of such data. WebMeV (http://mev.tm4.org) is an open-source, web-based tool that gives users access to sophisticated tools for analysis of RNA-Seq and other data in an interface designed to democratize data access. WebMeV combines cloud-based technologies with a simple user interface to allow users to access large public datasets, such as that from The Cancer Genome Atlas or to upload their own. The interface allows users to visualize data and to apply advanced data mining analysis methods to explore the data and draw biologically meaningful conclusions. We provide an overview of WebMeV and demonstrate two simple use cases that illustrate the value of putting data analysis in the hands of those looking to explore the underlying biology of the systems being studied. Cancer Res; 77(21); e11-14. ©2017 AACR.

High Throughput Sequencing of Extracellular RNA from Human Plasma.

The presence and relative stability of extracellular RNAs (exRNAs) in biofluids has led to an emerging recognition of their promise as 'liquid biopsies' for diseases. Most prior studies on discovery of exRNAs as disease-specific biomarkers have focused on microRNAs (miRNAs) using technologies such as qRT-PCR and microarrays. The recent application of next-generation sequencing to discovery of exRNA biomarkers has revealed the presence of potential novel miRNAs as well as other RNA species such as tRNAs, snoRNAs, piRNAs and lncRNAs in biofluids. At the same time, the use of RNA sequencing for biofluids poses unique challenges, including low amounts of input RNAs, the presence of exRNAs in different compartments with varying degrees of vulnerability to isolation techniques, and the high abundance of specific RNA species (thereby limiting the sensitivity of detection of less abundant species). Moreover, discovery in human diseases often relies on archival biospecimens of varying age and limiting amounts of samples. In this study, we have tested RNA isolation methods to optimize profiling exRNAs by RNA sequencing in individuals without any known diseases. Our findings are consistent with other recent studies that detect microRNAs and ribosomal RNAs as the major exRNA species in plasma. Similar to other recent studies, we found that the landscape of biofluid microRNA transcriptome is dominated by several abundant microRNAs that appear to comprise conserved extracellular miRNAs. There is reasonable correlation of sets of conserved miRNAs across biological replicates, and even across other data sets obtained at different investigative sites. Conversely, the detection of less abundant miRNAs is far more dependent on the exact methodology of RNA isolation and profiling. This study highlights the challenges in detecting and quantifying less abundant plasma miRNAs in health and disease using RNA sequencing platforms.