Shared brain transcriptomic signature in TDP-43 type A FTLD patients with or without GRN mutations.

Frontotemporal lobar degeneration with TDP-43 inclusions (FTLD-TDP) is a complex heterogeneous neurodegenerative disorder for which mechanisms are poorly understood. To explore transcriptional changes underlying FTLD-TDP, we performed RNA-sequencing on 66 genetically unexplained FTLD-TDP patients, 24 FTLD-TDP patients with GRN mutations and 24 control participants. Using principal component analysis, hierarchical clustering, differential expression and coexpression network analyses, we showed that GRN mutation carriers and FTLD-TDP-A patients without a known mutation shared a common transcriptional signature that is independent of GRN loss-of-function. After combining both groups, differential expression as compared to the control group and coexpression analyses revealed alteration of processes related to immune response, synaptic transmission, RNA metabolism, angiogenesis and vesicle-mediated transport. Deconvolution of the data highlighted strong cellular alterations that were similar in FTLD-TDP-A and GRN mutation carriers with NSF as a potentially important player in both groups. We propose several potentially druggable pathways such as the GABAergic, GDNF and sphingolipid pathways. Our findings underline new disease mechanisms and strongly suggest that affected pathways in GRN mutation carriers extend beyond GRN and contribute to genetically unexplained forms of FTLD-TDP-A.

ABCA7 haplodeficiency disturbs microglial immune responses in the mouse brain.

Carrying premature termination codons in 1 allele of the ABCA7 gene is associated with an increased risk for Alzheimer's disease (AD). While the primary function of ABCA7 is to regulate the transport of phospholipids and cholesterol, ABCA7 is also involved in maintaining homeostasis of the immune system. Since inflammatory pathways causatively or consequently participate in AD pathogenesis, we studied the effects of Abca7 haplodeficiency in mice on brain immune responses under acute and chronic conditions. When acute inflammation was induced through peripheral lipopolysaccharide injection in control or heterozygous Abca7 knockout mice, partial ABCA7 deficiency diminished proinflammatory responses by impairing CD14 expression in the brain. On breeding to AppNL-G-F knockin mice, we observed increased amyloid-ß (Aß) accumulation and abnormal endosomal morphology in microglia. Taken together, our results demonstrate that ABCA7 loss of function may contribute to AD pathogenesis by altering proper microglial responses to acute inflammatory challenges and during the development of amyloid pathology, providing insight into disease mechanisms and possible treatment strategies.

Somatic selection distinguishes oncogenes and tumor suppressor genes.

MOTIVATION: Functions of cancer driver genes vary substantially across tissues and organs. Distinguishing passenger genes, oncogenes (OGs) and tumor-suppressor genes (TSGs) for each cancer type is critical for understanding tumor biology and identifying clinically actionable targets. Although many computational tools are available to predict putative cancer driver genes, resources for context-aware classifications of OGs and TSGs are limited. RESULTS: We show that the direction and magnitude of somatic selection of protein-coding mutations are significantly different for passenger genes, OGs and TSGs. Based on these patterns, we develop a new method (genes under selection in tumors) to discover OGs and TSGs in a cancer-type specific manner. Genes under selection in tumors shows a high accuracy (92%) when evaluated via strict cross-validations. Its application to 10 172 tumor exomes found known and novel cancer drivers with high tissue-specificities. In 11 out of 13 OGs shared among multiple cancer types, we found functional domains selectively engaged in different cancers, suggesting differences in disease mechanisms. AVAILABILITY AND IMPLEMENTATION: An R implementation of the GUST algorithm is available at https://github.com/liliulab/gust. A database with pre-computed results is available at https://liliulab.shinyapps.io/gust. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Apolipoprotein E regulates lipid metabolism and α-synuclein pathology in human iPSC-derived cerebral organoids.

APOE4 is a strong genetic risk factor for Alzheimer's disease and Dementia with Lewy bodies; however, how its expression impacts pathogenic pathways in a human-relevant system is not clear. Here using human iPSC-derived cerebral organoid models, we find that APOE deletion increases α-synuclein (αSyn) accumulation accompanied with synaptic loss, reduction of GBA levels, lipid droplet accumulation and dysregulation of intracellular organelles. These phenotypes are partially rescued by exogenous apoE2 and apoE3, but not apoE4. Lipidomics analysis detects the increased fatty acid utilization and cholesterol ester accumulation in apoE-deficient cerebral organoids. Furthermore, APOE4 cerebral organoids have increased αSyn accumulation compared to those with APOE3. Carrying APOE4 also increases apoE association with Lewy bodies in postmortem brains from patients with Lewy body disease. Our findings reveal the predominant role of apoE in lipid metabolism and αSyn pathology in iPSC-derived cerebral organoids, providing mechanistic insights into how APOE4 drives the risk for synucleinopathies.

ApoE (Apolipoprotein E) in Brain Pericytes Regulates Endothelial Function in an Isoform-Dependent Manner by Modulating Basement Membrane Components.

OBJECTIVE: The ε4 allele of the APOE gene (APOE4) is the strongest genetic risk factor for Alzheimer disease when compared with the common ε3 allele. Although there has been significant progress in understanding how apoE4 (apolipoprotein E4) drives amyloid pathology, its effects on amyloid-independent pathways, in particular cerebrovascular integrity and function, are less clear. Approach and Results: Here, we show that brain pericytes, the mural cells of the capillary walls, differentially modulate endothelial cell phenotype in an apoE isoform-dependent manner. Extracellular matrix protein induction, tube-like structure formation, and barrier formation were lower with endothelial cells cocultured with pericytes isolated from apoE4-targeted replacement (TR) mice compared with those from apoE3-TR mice. Importantly, aged apoE4-targeted replacement mice had decreased extracellular matrix protein expression and increased plasma protein leakages compared with apoE3-TR mice. CONCLUSIONS: ApoE4 impairs pericyte-mediated basement membrane formation, potentially contributing to the cerebrovascular effects of apoE4.

NucMap: a database of genome-wide nucleosome positioning map across species.

Dynamics of nucleosome positioning affects chromatin state, transcription and all other biological processes occurring on genomic DNA. While MNase-Seq has been used to depict nucleosome positioning map in eukaryote in the past years, nucleosome positioning data is increasing dramatically. To facilitate the usage of published data across studies, we developed a database named nucleosome positioning map (NucMap, http://bigd.big.ac.cn/nucmap). NucMap includes 798 experimental data from 477 samples across 15 species. With a series of functional modules, users can search profile of nucleosome positioning at the promoter region of each gene across all samples and make enrichment analysis on nucleosome positioning data in all genomic regions. Nucleosome browser was built to visualize the profiles of nucleosome positioning. Users can also visualize multiple sources of omics data with the nucleosome browser and make side-by-side comparisons. All processed data in the database are freely available. NucMap is the first comprehensive nucleosome positioning platform and it will serve as an important resource to facilitate the understanding of chromatin regulation.

Identification of missing variants by combining multiple analytic pipelines.

BACKGROUND: After decades of identifying risk factors using array-based genome-wide association studies (GWAS), genetic research of complex diseases has shifted to sequencing-based rare variants discovery. This requires large sample sizes for statistical power and has brought up questions about whether the current variant calling practices are adequate for large cohorts. It is well-known that there are discrepancies between variants called by different pipelines, and that using a single pipeline always misses true variants exclusively identifiable by other pipelines. Nonetheless, it is common practice today to call variants by one pipeline due to computational cost and assume that false negative calls are a small percent of total. RESULTS: We analyzed 10,000 exomes from the Alzheimer's Disease Sequencing Project (ADSP) using multiple analytic pipelines consisting of different read aligners and variant calling strategies. We compared variants identified by using two aligners in 50,100, 200, 500, 1000, and 1952 samples; and compared variants identified by adding single-sample genotyping to the default multi-sample joint genotyping in 50,100, 500, 2000, 5000 and 10,000 samples. We found that using a single pipeline missed increasing numbers of high-quality variants correlated with sample sizes. By combining two read aligners and two variant calling strategies, we rescued 30% of pass-QC variants at sample size of 2000, and 56% at 10,000 samples. The rescued variants had higher proportions of low frequency (minor allele frequency [MAF] 1-5%) and rare (MAF < 1%) variants, which are the very type of variants of interest. In 660 Alzheimer's disease cases with earlier onset ages of ≤65, 4 out of 13 (31%) previously-published rare pathogenic and protective mutations in APP, PSEN1, and PSEN2 genes were undetected by the default one-pipeline approach but recovered by the multi-pipeline approach. CONCLUSIONS: Identification of the complete variant set from sequencing data is the prerequisite of genetic association analyses. The current analytic practice of calling genetic variants from sequencing data using a single bioinformatics pipeline is no longer adequate with the increasingly large projects. The number and percentage of quality variants that passed quality filters but are missed by the one-pipeline approach rapidly increased with sample size.

Genetically predicted longer telomere length is associated with increased risk of B-cell lymphoma subtypes.

Evidence from a small number of studies suggests that longer telomere length measured in peripheral leukocytes is associated with an increased risk of non-Hodgkin lymphoma (NHL). However, these studies may be biased by reverse causation, confounded by unmeasured environmental exposures and might miss time points for which prospective telomere measurement would best reveal a relationship between telomere length and NHL risk. We performed an analysis of genetically inferred telomere length and NHL risk in a study of 10 102 NHL cases of the four most common B-cell histologic types and 9562 controls using a genetic risk score (GRS) comprising nine telomere length-associated single-nucleotide polymorphisms. This approach uses existing genotype data and estimates telomere length by weighing the number of telomere length-associated variant alleles an individual carries with the published change in kb of telomere length. The analysis of the telomere length GRS resulted in an association between longer telomere length and increased NHL risk [four B-cell histologic types combined; odds ratio (OR) = 1.49, 95% CI 1.22-1.82,P-value = 8.5 × 10(-5)]. Subtype-specific analyses indicated that chronic lymphocytic leukemia or small lymphocytic lymphoma (CLL/SLL) was the principal NHL subtype contributing to this association (OR = 2.60, 95% CI 1.93-3.51,P-value = 4.0 × 10(-10)). Significant interactions were observed across strata of sex for CLL/SLL and marginal zone lymphoma subtypes as well as age for the follicular lymphoma subtype. Our results indicate that a genetic background that favors longer telomere length may increase NHL risk, particularly risk of CLL/SLL, and are consistent with earlier studies relating longer telomere length with increased NHL risk.

Human Pegivirus infection and lymphoma risk and prognosis: a North American study.

We evaluated the association of Human Pegivirus (HPgV) viraemia with risk of developing lymphoma, overall and by major subtypes. Because this virus has also been associated with better prognosis in the setting of co-infection with human immunodeficiency virus, we further assessed the association of HPgV with prognosis. We used risk factor data and banked plasma samples from 2094 lymphoma cases newly diagnosed between 2002 and 2009 and 1572 frequency-matched controls. Plasma samples were tested for HPgV RNA by reverse transcription polymerase chain reaction (RT-PCR), and those with RNA concentrations <5000 genome equivalents/ml were confirmed using nested RT-PCR methods. To assess the role of HPgV in lymphoma prognosis, we used 2948 cases from a cohort study of newly diagnosed lymphoma patients (included all cases from the case-control study). There was a positive association of HPgV viraemia with risk of lymphoma overall (Odds ratio = 2·14; 95% confidence interval [CI] 1·63-2·80; P < 0·0001), and for all major subtypes except Hodgkin lymphoma and chronic lymphocytic leukaemia/small lymphocytic lymphoma, and this was not confounded by other lymphoma risk factors. In contrast, there was no association of HPgV viraemia with event-free survival (Hazard ratio [HR] = 1·00; 95% CI 0·85-1·18) or overall survival (HR = 0·97; 95% CI 0·79-1·20) for lymphoma overall, or any of the subtypes. These data support the hypothesis for a role of HPgV in the aetiology of multiple lymphoma subtypes.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Divergent brain gene expression patterns associate with distinct cell-specific tau neuropathology traits in progressive supranuclear palsy.

Progressive supranuclear palsy (PSP) is a neurodegenerative parkinsonian disorder characterized by tau pathology in neurons and glial cells. Transcriptional regulation has been implicated as a potential mechanism in conferring disease risk and neuropathology for some PSP genetic risk variants. However, the role of transcriptional changes as potential drivers of distinct cell-specific tau lesions has not been explored. In this study, we integrated brain gene expression measurements, quantitative neuropathology traits and genome-wide genotypes from 268 autopsy-confirmed PSP patients to identify transcriptional associations with unique cell-specific tau pathologies. We provide individual transcript and transcriptional network associations for quantitative oligodendroglial (coiled bodies = CB), neuronal (neurofibrillary tangles = NFT), astrocytic (tufted astrocytes = TA) tau pathology, and tau threads and genomic annotations of these findings. We identified divergent patterns of transcriptional associations for the distinct tau lesions, with the neuronal and astrocytic neuropathologies being the most different. We determined that NFT are positively associated with a brain co-expression network enriched for synaptic and PSP candidate risk genes, whereas TA are positively associated with a microglial gene-enriched immune network. In contrast, TA is negatively associated with synaptic and NFT with immune system transcripts. Our findings have implications for the diverse molecular mechanisms that underlie cell-specific vulnerability and disease risk in PSP.

Chromoanasynthesis is a common mechanism that leads to ERBB2 amplifications in a cohort of early stage HER2+ breast cancer samples.

BACKGROUND: HER2 positive (HER2+) breast cancers involve chromosomal structural alterations that act as oncogenic driver events. METHODS: We interrogated the genomic structure of 18 clinically-defined HER2+ breast tumors through integrated analysis of whole genome and transcriptome sequencing, coupled with clinical information. RESULTS: ERBB2 overexpression in 15 of these tumors was associated with ERBB2 amplification due to chromoanasynthesis with six of them containing single events and the other nine exhibiting multiple events. Two of the more complex cases had adverse clinical outcomes. Chromosomes 8 was commonly involved in the same chromoanasynthesis with 17. In ten cases where chromosome 8 was involved we observed NRG1 fusions (two cases), NRG1 amplification (one case), FGFR1 amplification and ADAM32 or ADAM5 fusions. ERBB3 over-expression was associated with NRG1 fusions and EGFR and ERBB3 expressions were anti-correlated. Of the remaining three cases, one had a small duplication fully encompassing ERBB2 and was accompanied with a pathogenic mutation. CONCLUSION: Chromoanasynthesis involving chromosome 17 can lead to ERBB2 amplifications in HER2+ breast cancer. However, additional large genomic alterations contribute to a high level of genomic complexity, generating the hypothesis that worse outcome could be associated with multiple chromoanasynthetic events.

Conserved brain myelination networks are altered in Alzheimer's and other neurodegenerative diseases.

INTRODUCTION: Comparative transcriptome analyses in Alzheimer's disease (AD) and other neurodegenerative proteinopathies can uncover both shared and distinct disease pathways. METHODS: We analyzed 940 brain transcriptomes including patients with AD, progressive supranuclear palsy (PSP; a primary tauopathy), and control subjects. RESULTS: We identified transcriptional coexpression networks implicated in myelination, which were lower in PSP temporal cortex (TCX) compared with AD. Some of these associations were retained even after adjustments for brain cell population changes. These TCX myelination network structures were preserved in cerebellum but they were not differentially expressed in cerebellum between AD and PSP. Myelination networks were downregulated in both AD and PSP, when compared with control TCX samples. DISCUSSION: Downregulation of myelination networks may underlie both PSP and AD pathophysiology, but may be more pronounced in PSP. These data also highlight conservation of transcriptional networks across brain regions and the influence of cell type changes on these networks.

Genome-wide association study of cardiotoxicity in the NCCTG N9831 (Alliance) adjuvant trastuzumab trial.

OBJECTIVES: The major clinical side effect of the ERBB2-targeted breast cancer therapy, trastuzumab, is a decline in the left ventricular ejection fraction (LVEF). Improved markers are needed to better identify patients susceptible to cardiotoxicity. METHODS: The NCCTG N9831 trial compared adjuvant doxorubicin and cyclophosphamide followed by either weekly paclitaxel (arm A); paclitaxel then trastuzumab (arm B); or concurrent paclitaxel and trastuzumab (arm C) in patients with HER2-positive breast cancer. A genome-wide association study was performed on all patients with available DNA (N=1446). We used linear regression to identify single nucleotide polymorphisms (SNPs) associated with decline in LVEF, adjusting for age, baseline LVEF, antihypertensive medications, and the first two principle components. RESULTS: In total, 618 863 SNPs passed quality control and DNA from 1191 patients passed genotyping quality control and were identified as Whites of non-Hispanic origin. SNPs at six loci were associated with a decline in LVEF (P=7.73×10 to 8.93×10), LDB2, BRINP1, chr6 intergenic, RAB22A, TRPC6, and LINC01060, in patients who received chemotherapy plus trastuzumab (arms BC, N=800). None of these loci were significant in patients who received chemotherapy only (arm A, N=391) and did not increase in significance in the combined analysis of all patients. We did not observe association, P<0.05, with SNPs previously associated with trastuzumab-induced cardiotoxicity at ERBB2, I655V, and P1170A. We replicated association, P<0.05, with SNPs previously associated with anthracycline-induced cardiotoxicity at CBR3 and ABCB1. CONCLUSION: Our study identified six putative novel cardiotoxicity loci in patients treated with combination chemotherapy and trastuzumab that require further investigation and confirmed known associations of anthracycline-induced cardiotoxicity.

In-depth clinico-pathological examination of RNA foci in a large cohort of C9ORF72 expansion carriers.

A growing body of evidence suggests that a loss of chromosome 9 open reading frame 72 (C9ORF72) expression, formation of dipeptide-repeat proteins, and generation of RNA foci contribute to disease pathogenesis in amyotrophic lateral sclerosis and frontotemporal dementia. Although the levels of C9ORF72 transcripts and dipeptide-repeat proteins have already been examined thoroughly, much remains unknown about the role of RNA foci in C9ORF72-linked diseases. As such, we performed a comprehensive RNA foci study in an extensive pathological cohort of C9ORF72 expansion carriers (n = 63). We evaluated two brain regions using a newly developed computer-automated pipeline allowing recognition of cell nuclei and RNA foci (sense and antisense) supplemented by manual counting. In the frontal cortex, the percentage of cells with sense or antisense RNA foci was 26 or 12%, respectively. In the cerebellum, 23% of granule cells contained sense RNA foci and 1% antisense RNA foci. Interestingly, the highest percentage of cells with RNA foci was observed in cerebellar Purkinje cells (~70%). In general, more cells contained sense RNA foci than antisense RNA foci; however, when antisense RNA foci were present, they were usually more abundant. We also observed that an increase in the percentage of cells with antisense RNA foci was associated with a delayed age at onset in the frontal cortex (r = 0.43, p = 0.003), whereas no other associations with clinico-pathological features were seen. Importantly, our large-scale study is the first to provide conclusive evidence that RNA foci are not the determining factor of the clinico-pathological variability observed in C9ORF72 expansion carriers and it emphasizes that the distribution of RNA foci does not follow the pattern of neurodegeneration, stressing the complex interplay between different aspects of C9ORF72-related diseases.

TP53 mutations, tetraploidy and homologous recombination repair defects in early stage high-grade serous ovarian cancer.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

A candidate regulatory variant at the TREM gene cluster associates with decreased Alzheimer's disease risk and increased TREML1 and TREM2 brain gene expression.

INTRODUCTION: We hypothesized that common Alzheimer's disease (AD)-associated variants within the triggering receptor expressed on myeloid (TREM) gene cluster influence disease through gene expression. METHODS: Expression microarrays on temporal cortex and cerebellum from ∼400 neuropathologically diagnosed subjects and two independent RNAseq replication cohorts were used for expression quantitative trait locus analysis. RESULTS: A variant within a DNase hypersensitive site 5' of TREM2, rs9357347-C, associates with reduced AD risk and increased TREML1 and TREM2 levels (uncorrected P = 6.3 × 10-3 and 4.6 × 10-2, respectively). Meta-analysis on expression quantitative trait locus results from three independent data sets (n = 1006) confirmed these associations (uncorrected P = 3.4 × 10-2 and 3.5 × 10-3, Bonferroni-corrected P = 6.7 × 10-2 and 7.1 × 10-3, respectively). DISCUSSION: Our findings point to rs9357347 as a functional regulatory variant that contributes to a protective effect observed at the TREM locus in the International Genomics of Alzheimer's Project genome-wide association study meta-analysis and suggest concomitant increase in TREML1 and TREM2 brain levels as a potential mechanism for protection from AD.

The relationship between quantitative human epidermal growth factor receptor 2 gene expression by the 21-gene reverse transcriptase polymerase chain reaction assay and adjuvant trastuzumab benefit in Alliance N9831.

INTRODUCTION: The N9831 trial demonstrated the efficacy of adjuvant trastuzumab for patients with human epidermal growth factor receptor 2 (HER2) locally positive tumors by protein or gene analysis. We used the 21-gene assay to examine the association of quantitative HER2 messenger RNA (mRNA) gene expression and benefit from trastuzumab. METHODS: N9831 tested the addition of trastuzumab to chemotherapy in stage I-III HER2-positive breast cancer. For two of the arms of the trial, doxorubicin and cyclophosphamide followed by paclitaxel (AC-T) and doxorubicin and cyclophosphamide followed by paclitaxel and trastuzumab concurrent chemotherapy-trastuzumab (AC-TH), recurrence score (RS) and HER2 mRNA expression were determined by the 21-gene assay (Oncotype DX®) (negative <10.7, equivocal 10.7 to <11.5, and positive ≥11.5 log2 expression units). Cox regression was used to assess the association of HER2 expression with trastuzumab benefit in preventing distant recurrence. RESULTS: Median follow-up was 7.4 years. Of 1,940 total patients, 901 had consent and sufficient tissue. HER2 by reverse transcriptase polymerase chain reaction (RT-PCR) was negative in 130 (14 %), equivocal in 85 (9 %), and positive in 686 (76 %) patients. Concordance between HER2 assessments was 95 % for RT-PCR versus central immunohistochemistry (IHC) (>10 % positive cells = positive), 91 % for RT-PCR versus central fluorescence in situ hybridization (FISH) (≥2.0 = positive) and 94 % for central IHC versus central FISH. In the primary analysis, the association of HER2 expression by 21-gene assay with trastuzumab benefit was marginally nonsignificant (nonlinear p = 0.057). In hormone receptor-positive patients (local IHC) the association was significant (p = 0.002). The association was nonlinear with the greatest estimated benefit at lower and higher HER2 expression levels. CONCLUSIONS: Concordance among HER2 assessments by central IHC, FISH, and RT-PCR were similar and high. Association of HER2 mRNA expression with trastuzumab benefit as measured by time to distant recurrence was nonsignificant. A consistent benefit of trastuzumab irrespective of mHER2 levels was observed in patients with either IHC-positive or FISH-positive tumors. Trend for benefit was observed also for the small groups of patients with negative results by any or all of the central assays. TRIAL REGISTRATION: Clinicaltrials.gov NCT00005970 . Registered 5 July 2000.

RVboost: RNA-seq variants prioritization using a boosting method.

MOTIVATION: RNA-seq has become the method of choice to quantify genes and exons, discover novel transcripts and detect fusion genes. However, reliable variant identification from RNA-seq data remains challenging because of the complexities of the transcriptome, the challenges of accurately mapping exon boundary spanning reads and the bias introduced during the sequencing library preparation. METHOD: We developed RVboost, a novel method specific for RNA variant prioritization. RVboost uses several attributes unique in the process of RNA library preparation, sequencing and RNA-seq data analyses. It uses a boosting method to train a model of 'good quality' variants using common variants from HapMap, and prioritizes and calls the RNA variants based on the trained model. We packaged RVboost in a comprehensive workflow, which integrates tools of variant calling, annotation and filtering. RESULTS: RVboost consistently outperforms the variant quality score recalibration from the Genome Analysis Tool Kit and the RNA-seq variant-calling pipeline SNPiR in 12 RNA-seq samples using ground-truth variants from paired exome sequencing data. Several RNA-seq-specific attributes were identified as critical to differentiate true and false variants, including the distance of the variant positions to exon boundaries, and the percent of the reads supporting the variant in the first six base pairs. The latter identifies false variants introduced by the random hexamer priming during the library construction. AVAILABILITY AND IMPLEMENTATION: The RVboost package is implemented to readily run in Mac or Linux environments. The software and user manual are available at http://bioinformaticstools.mayo.edu/research/rvboost/.