Back to Authors

Chen, Xianfeng

https://orcid.org/0000-0002-7531-5862

Mayo Clinic

Mayo Clinic

Citations 21

Synergistic combination of cytotoxic chemotherapy and cyclin‐dependent kinase 4/6 inhibitors in biliary tract cancers
Arora, Mansi , Abdelrahman, Amro M. , Alva‐Ruiz, Roberto, Leiting, Jennifer L. , Chen, Xianfeng , Usón, Pedro Luiz Serrano , Dumbauld, Chelsae, Baker, Alexander T. , Gamb, Scott I., Egan, Jan B. , Zhou, Yumei , Nagalo, Bolni Marius , Meurice, Nathalie , Eskelinen, Eeva‐Liisa , Braggio, Esteban , Barrett, Michael T. , Buetow, Kenneth H. , Roberts, Lewis R. , Bekaii‐Saab, Tanios , Ahn, Daniel H. , Truty, Mark J. , Bogenberger, James M. , Yonkus, Jennifer, Alva-Ruiz, Roberto, Leiting, Jennifer L , Serrano Uson Junior, Pedro Luiz, Dumbauld, Chelsae R, Salomao, Marcela A, Roberts, Lewis R , Braggio, Esteban , Uson, Pedro L S , Bekaii-Saab, Tanios S , Nagalo, Bolni , Baker, Alexander , Abdelrahman, Amro M , Yonkus, Jennifer A. , Salomao, Marcela , Kosiorek, Heidi , Sonbol, Mohamad Bassam , Mansfield, Aaron S. , Borad, Mitesh J.

Hepatology | August 18, 2021 | Volume: 75 | Issue: 1 | Pages: 43-58

10.1002/hep.32102 Get It

Abstract Background and aims Biliary tract cancers (BTCs) are uncommon, but highly lethal, gastrointestinal malignancies. Gemcitabine/cisplatin is a standard‐of‐care systemic therapy, but has a modest impact on survival and harbors toxicities, including myelosuppression, nephropathy, neuropathy, and ototoxicity. Whereas BTCs are characterized by aberrations activating the cyclinD1/cyclin‐dependent kinase (CDK)4/6/CDK inhibitor 2a/retinoblastoma pathway, clinical use of CDK4/6 inhibitors as monotherapy is limited by lack of validated biomarkers, diffident preclinical efficacy, and development of acquired drug resistance. Emerging studies have explored therapeutic strategies to enhance the antitumor efficacy of CDK4/6 inhibitors by the combination with chemotherapy regimens, but their mechanism of action remains elusive. Approach and results Here, we report in vitro and in vivo synergy in BTC models, showing enhanced efficacy, reduced toxicity, and better survival with a combination comprising gemcitabine/cisplatin and CDK4/6 inhibitors. Furthermore, we demonstrated that abemaciclib monotherapy had only modest efficacy attributable to autophagy‐induced resistance. Notably, triplet therapy was able to potentiate efficacy through elimination of the autophagic flux. Correspondingly, abemaciclib potentiated ribonucleotide reductase catalytic subunit M1 reduction, resulting in sensitization to gemcitabine. Conclusions As such, these data provide robust preclinical mechanistic evidence of synergy between gemcitabine/cisplatin and CDK4/6 inhibitors and delineate a path forward for translation of these findings to preliminary clinical studies in advanced BTC patients.
Novel lincRNA Discovery and Tissue-Specific Gene Expression across 30 Normal Human Tissues
Chen, Xianfeng , Sun, Zhifu , Sun, Zhifu

Genes | April 21, 2021 | Volume: 12 | Issue: 5 | Pages: 614

10.3390/genes12050614 Get It

Long non-coding RNAs (lncRNAs) are a large class of gene transcripts that do not code proteins; however, their functions are largely unknown and many new lncRNAs are yet to be discovered. Taking advantage of our previously developed, super-fast, novel lncRNA discovery pipeline, UClncR, and rich resources of GTEx RNA-seq data, we performed systematic novel lincRNA discovery for over 8000 samples across 30 tissue types. We conducted novel detection for each major tissue type first and then consolidated the novel discoveries from all tissue types. These novel lincRNs were profiled and analyzed along with known genes to identify tissue-specific genes in 30 major human tissue types. Thirteen sub-brain regions were also analyzed in a similar manner. Our analysis revealed thousands to tens of thousands of novel lincRNAs for each tissue type. These lincRNAs could define each tissue type’s identity and demonstrated their reliability and tissue-specific expression. Tissue-specific genes were identified for each major tissue type and sub-brain region. The tissue-specific genes clearly defined each respective tissue’s unique function and could be used to expand the interpretation of non-coding SNPs from genome-wide association (GWAS) studies.
Placental growth factor promotes tumour desmoplasia and treatment resistance in intrahepatic cholangiocarcinoma
Aoki, Shuichi , Inoue, Koetsu , Klein, Sebastian , Halvorsen, Stefan , Matsui, Aya , Nikmaneshi, Mohammad R. , Hato, Tai , Chen, Xianfeng , Nia, Hadi T. , Chen, Ivy, Schanne, Daniel H. , Mamessier, Émilie , Shigeta, Kohei , Ramjiawan, Rakesh R. , Schmidt, Tyge CE , Iwasaki, M, Hong, Theodore S. , Plum, Patrick Sven , Dima, Simona , Popescu, Irinel , Bardeesy, Nabeel, Munn, Lance L. , Sassi, Slim , Jain, Rakesh K., Zhu, Andrew X. , Duda, Dan G. , Shigeta, Kohei , Yau, Thomas , Duda, Dan G , Chen, Jiang , Kitahara, Shuji , Kawakubo, Kazumichi , Kikuchi, Hiroto , Iwasaki, Masaaki, Quaas, Alexander , Borad, Mitesh J. , Jain, Rakesh K.

Gut | January 12, 2021 | Volume: 71 | Issue: 1 | Pages: 185-193

10.1136/gutjnl-2020-322493 Get It

Objective Intrahepatic cholangiocarcinoma (ICC)—a rare liver malignancy with limited therapeutic options—is characterised by aggressive progression, desmoplasia and vascular abnormalities. The aim of this study was to determine the role of placental growth factor (PlGF) in ICC progression. Design We evaluated the expression of PlGF in specimens from ICC patients and assessed the therapeutic effect of genetic or pharmacologic inhibition of PlGF in orthotopically grafted ICC mouse models. We evaluated the impact of PlGF stimulation or blockade in ICC cells and cancer-associated fibroblasts (CAFs) using in vitro 3-D coculture systems. Results PlGF levels were elevated in human ICC stromal cells and circulating blood plasma and were associated with disease progression. Single-cell RNA sequencing showed that the major impact of PlGF blockade in mice was enrichment of quiescent CAFs, characterised by high gene transcription levels related to the Akt pathway, glycolysis and hypoxia signalling. PlGF blockade suppressed Akt phosphorylation and myofibroblast activation in ICC-derived CAFs. PlGF blockade also reduced desmoplasia and tissue stiffness, which resulted in reopening of collapsed tumour vessels and improved blood perfusion, while reducing ICC cell invasion. Moreover, PlGF blockade enhanced the efficacy of standard chemotherapy in mice-bearing ICC. Conclusion PlGF blockade leads to a reduction in intratumorous hypoxia and metastatic dissemination, enhanced chemotherapy sensitivity and increased survival in mice-bearing aggressive ICC.
Over-Expression of Transferrin Receptor (TFRC/CD71) and Low Expression of Innate and Adaptive Immune Cell Subsets in HIV-Associated, GCB-DLBCL By Digital Gene Expression Profiling
Maguire, Alanna , Chen, Xianfeng , Wisner, Lee, Ramsower, Colleen , Glinsmann‐Gibson, Betty, Maguire, Alanna , Rimsza, Lisa M. , Maguire, Alanna

Blood | November 14, 2019 | Volume: 134 | Issue: Supplement_1 | Pages: 2783-2783

10.1182/blood-2019-129346 Get It

Introduction: HIV infected individuals are 17x more likely to receive a diagnosis of Diffuse Large B cell Lymphoma (DLBCL) compared to their uninfected counterparts. Moreover, DLBCL is more aggressive in HIV infected individuals. However, the molecular pathology driving the aggressive nature of HIV related DLBCL is poorly understood. Previously, we demonstrated that HIV-associated [HIV(+)] germinal center B-cell (GCB) DLBCL tumors are more proliferative, with enhanced genomic stability and increased expression of DNA repair genes compared to their HIV-not associated [HIV(-)] counterparts (Maguire et al, Int J Cancer,2019). Given the immunocompromised nature of these patients, herein we assessed whether specific immunological signaling pathways are also altered in HIV(+) GCB DLBCL tumors, that may work in conjunction with, or independently of, increased DNA repair, to drive the enhanced aggressive nature of HIV related DLBCL. Methods: We used our previously studied cohort of patient samples, which included a total of 39 molecularly defined GCB DLBCL cases. Of these, 18 were HIV(+) cases from the AIDS Cancer Specimen Resource Network (https://acsr.ucsf.edu/) and 21 were HIV(-) institutional cases. Samples underwent hematopathologist review for diagnosis confirmation and tumor content assessment. Samples with <60% tumor content were macrodissected, with a total of 4x 5µm FFPE sections per sample used for RNA/DNA extraction. Digital gene expression profiling was used to assess differential expression of 579 immune related genes using the NanoString PanCancer Immune Profiling panel. The data was normalized using nSolver, and differential expression analysis was performed using the R statistical software package NanoStringDiff. The resulting expression data were then analyzed using Gene Set Enrichment Analysis (GSEA), and GSEA Molecular Signatures Database (MSigDB) overlap analyses. Results: A total of 110 genes were differentially expressed between the HIV(+) and HIV(-) cohorts at p<0.05. As expected due to immunosuppression in HIV(+) individuals, 100 genes (91%) were significantly reduced and 10 genes (9%) were significantly increased in the HIV(+) cohort compared to the HIV(-) cohort. GSEA and MSigDB analyses revealed that the genes with reduced expression in the HIV(+) cohort were associated with losses in both adaptive and innate immune signaling including cytokine interactions, the IL1 pathway and B-cell receptor signaling amongst others. MSigDB overlap analyses of the 10 genes upregulated in the HIV(+) cohort identified 3 Gene Ontology molecular functions, all of which were receptor signaling related; including Receptor Binding (CD8A, CD160, CCL8, C1QBP, MIF), Receptor Activity (CD8A, CD160, CCRL1, TFRC, KLRC2), and Cargo Receptor Activity (CCRL1, TFRC). Gene overlaps between the subsets include genes associated with NK and cytolytic T cells (CD8A, CD160), as well as the transferrin receptor TRFC/CD71, which with a fold increase of 3.9, was the most significantly increased gene in the HIV(+) cohort (p=0.002, FDR=0.0025). Conclusions: As expected, the results demonstrate a loss in both adaptive and innate immune signaling in the HIV(+) cohort compared to the HIV negative cohort as well as alterations in receptor signaling. In line with our previous observation that HIV(+) GCB DLBCL tumors have higher expression of proliferation genes and the Ki67 antigen than their HIV(-) counterparts, the results of this study also reveal overexpression of TRFC/CD71 mRNA in HIV(+) GCB DLBCL tumors compared to the HIV(-) cases. TFRC/CD71 is well known as a mediator of iron uptake in erythroid cells. Iron uptake through TFRC/CD71 is also necessary for, and positively regulates T and B cell proliferation; and is commonly expressed on aggressive B cell lymphomas. We suggest that TFRC/CD71 may play a role in the enhanced proliferative capabilities of HIV(+) GCB DLBCL. These effects may be mediated through the HIV protein Nef, which has known effects on TFRC/CD71 protein recycling in lymphoid populations (Madrid et al, J Biol Chem 2005), and has been shown to be transferred from HIV infected macrophages to uninfected B-cells via contact-dependent cellular conduits (Xu et al, Nat Immunol 2009). Nef transfer between cells has also been described to occur via microvesicle transfer and trogocytosis (reviewed in Ellwanger et al, Infect Genet Evol 2017). Further analysis of these pathways is ongoing. Disclosures Rimsza: NanoSting: Patents & Royalties: Named inventor on a patent licensed to NanoSting [Institution].
Enhanced DNA Repair and Genomic Stability in HIV(+) Diffuse Large B Cell Lymphoma of Germinal Center Origin
Maguire, Alanna , Chen, Xianfeng , Wisner, Lee, Ramsower, Colleen , Glinsmann‐Gibson, Betty, McGrath, Michael S. , Maguire, Alanna , Rimsza, Lisa M. , Maguire, Alanna

Blood | October 14, 2019 | Volume: 132 | Issue: Supplement 1 | Pages: 1570-1570

10.1182/blood-2018-99-116571 Get It

Abstract Introduction: HIV infected individuals are 17x more likely to receive a diagnoses of Diffuse Large B cell Lymphoma (DLBCL) compared to their uninfected counterparts. Moreover, DLBCL is more aggressive in HIV infected individuals, with up to 70% of patients being primary refractory to chemotherapeutic regimens. However, the molecular pathology driving the aggressive nature of HIV related DLBCL is poorly understood. Here, we have assessed the genomic and transcriptional differences between HIV(+) and HIV(-) DLBCL in order to identify the mechanisms driving the enhanced aggressive and refractory nature of HIV related DLBCL. Methods: A total of 66 cases, including 27 HIV(+) from the AIDS Cancer Specimen Resource Network (https://acsr.ucsf.edu/) and 39 HIV(-) institutional cases were included in this study. Fresh H&Es were reviewed by a hematopathologist to validate diagnosis and determine tumor content. Samples with <60% tumor content were macro-dissected. A total of 4x 5µm FFPE sections per sample were used for DNA and RNA extraction. DNA damage at the gross aberration level was assessed using array comparative genomic hybridization (aCGH). NanoString digital gene expression profiling (GEP) was used to perform the Lymph2Cx DLBCL Cell of Origin molecular subtyping assay (Germinal Center B-cell/ GCB, Activated B-Cell/ ABC, Unclassifiable/UNC) and to assess the differential expression of genes with known roles in cancer using the PanCancer Pathways panel. Results: Both the HIV(+) and HIV(-) cohorts were found to be GCB enriched with GCB-ABC-UNC distributions of 70%(19/27)-11%(3/27)-19%(5/27) and 54%(21/39)-33%(13/39)-13%(5/39) respectively. For both statistical and biological reasons, analysis was restricted to the GCB-DLBCL subtype only. Array CGH revealed the HIV(+) tumors have less DNA aberrations than their HIV(-) counterparts, indicative of enhanced genomic stability. GEP GSEA analysis revealed significant differences in the HIV(+) cohort compared to the HIV(-) cohort; including increases in DNA replication (MCM2, MCM4, MCM7) and Cell Cycle progression (CDC25A, CCNA2, CCNB1, E2F1) related genes as well as decreases in cell cycle negative regulators (CDKN1A, CDKN1B, CDKN2B, RB1), pro-apoptotic BCL2 genes (BAX, BIM, BMF, PUMA, BNIP3) and pro-survival BCL2 genes (BCL2, BCLW). Analyses also revealed significant increases in DNA repair related genes, particularly those with roles in the Fanconi Anemia-Homologous Recombination-Translesion Synthesis axis (FANCA, FANCD1/BRCA2, FANCE, FANCG, FANCR/RAD51, FANCS/BRCA1, FANCT, FANCV/MAD2L2). Conclusions: Our results show, for the first time, that GCB DLBCL arising in HIV infected individuals possesses a distinct molecular pathology to that which arises in uninfected individuals. These results indicate that the aggressive nature of HIV(+) GCB DLBCL is mediated by increased proliferation in conjunction with reduced cell cycle inhibitory capabilities and potentiated by heightened DNA repair that promotes genomic stability. Moreover, it is probable that the bolstered DNA repair capabilities confer an innate ability to repair DNA damage resulting from the administration of genotoxic chemotherapeutic agents and may be the mechanism underlying the high primary refractory rate to chemotherapeutics in HIV related DLBCL. Disclosures Rimsza: NanoString: Other: Inventor on the patent for the Lymph2Cx assay.
Enhanced Expression of FGF Signaling in Primary Central Nervous System Lymphoma
Maguire, Alanna , Chen, Xianfeng , Maguire, Alanna , Hilal, Talal , Rosenthal, Allison , Rimsza, Lisa M. , Maguire, Alanna

Blood | October 14, 2019 | Volume: 132 | Issue: Supplement 1 | Pages: 2847-2847

10.1182/blood-2018-99-116416 Get It

Abstract Introduction: Primary central nervous system lymphoma (PCNSL) is a rare intracranial lymphoma that accounts for less than 1% of all non-Hodgkins lymphomas and 3% of all brain tumors. Histopathologically, approximately 90% of PCNSL cases are categorized as a diffuse large B-Cell lymphoma (DLBCL). DLBCL malignancies are subdivided by Cell of Origin (COO), with the vast majority of PCNSL categorized as non-germinal center B cell (non-GCB) by immunohistochemistry. Gene expression profiling (GEP), however, has shown that immunohistochemically defined non-GCB resolves into two distinct subtypes, namely activated B-cell (ABC) and unclassified (UNC) subtypes. Using the Lymph2Cx molecular COO subtyping assay, we have found that 91% of PCNSL are ABC (unpublished data). Unlike systemic-DLBCL, PCNSL is largely confined to and rarely metastasizes outside of the immune privileged central nervous system. Despite this, PCNSL is one of the most aggressive forms of DLBCL. Given the immune privileged milieu in which PCNSL arises, we hypothesized that this milieu elicits a transcriptional profile that contributes to the enhanced aggressive nature of PCNSL compared to systemic-DLBCL. To investigate this hypothesis, this study assessed the gene expression differences between ABC-PCNSL and ABC-systemic-DLBCL, in order to identify novel players in the pathogenesis of ABC-PCNSL. Methods: A total of 35 HIV negative samples, with proven ABC-subtype COO as per the GEP Lymph2Cx assay, were employed; including 10 ABC systemic-DLBCL and 25 primary ABC PCNSL cases with no concurrent or prior history of systemic DLBCL. Samples were reviewed by a hematopathologist to confirm diagnoses and determine tumor content. Samples with <60% tumor content were macro-dissected before nucleic acid extraction, which was performed using the Qiagen AllPrep DNA/RNA FFPE Kit. Extracted DNA and RNA were quantified using the Qubit HS-kit and NanoDrop respectively. Digital gene expression technology was used to perform the PanCancer Pathways panel (NanoString, Seattle, WA). Differential gene expression analysis was performed using the NanoString specific statistical method NanoStringDiff. Identified gene sets were analyzed using the online Gene Set Enrichment Analysis (GSEA) Molecular Signatures Database (MSigDB). Results: Of the 739 cancer related genes targeted by the PanCancer panel, 256 were found to be significantly differentially expressed in the ABC-PCNSL cohort compared to the ABC-DLBCL cohort (p<0.05). Fifty six genes were upregulated and 200 were downregulated. With a 4.9 fold change, the most significantly overexpressed gene was FGF1 (p=4.7E-11). FGF1 encodes a primary ligand for the fibroblast growth factor receptors (FGFR) -1, -2, -3 and -4; of which, FGFR2 (p=1.0E-7) and FGFR3 (p=0.003) were also significantly overexpressed. Moreover, MSigDB identified the FGF signaling pathway as enriched in the upregulated gene set (5 genes, p=7.4E-9, FDR=6.6E-7). FGFRs are a family of receptors that activate known mitogenic signaling pathways including MAPK signaling, which MSigDB identified as the most enriched pathway in the upregulated gene set (14 genes, p=1.65E-19, FDR=2.2E-16). MSigDB analysis of the 200 down regulated genes revealed that 5 of the top 20 enriched signaling pathways were immune related and included Signaling by interleukins (26 genes, p=2.9E-38, FDR=3.2E-36), Immune cytokine signaling (31 genes, p=1.1 E-34, FDR=1.1E-32), chemokine signaling (28 genes, p=1.5E-34, FDR=1.4E-32), T-cell receptor signaling (24 genes, p=2.4E-34, FDR=1.9E-32) and Toll-like receptor signaling (23 genes, p=4.1E-33, FDR=3.0E-31). Conclusions: We show, for the first time, that ABC-PCNSL and ABC-systemic-DLBCL possess significantly different transcriptional profiles despite identical, molecularly determined, COO status. A principle difference between these DLBCL malignancies is their anatomical location related immune privilege status which is reflected as reduced immune related signaling in the CNS-DLBCL cohort and may have important mitogenic signaling implications. Indeed, the results suggest that the enhanced aggressive nature of PCNSL compared to systemic-DLBCL is mediated, at least in part, by enhanced FGF signaling; a pathway with known roles in cell survival and proliferation. Disclosures Rimsza: NanoString: Other: Inventor on the patent for the Lymph2Cx assay.
Enhanced DNA repair and genomic stability identify a novel HIV‐related diffuse large B‐cell lymphoma signature
Maguire, Alanna , Chen, Xianfeng , Wisner, Lee, Malasi, Smriti, Ramsower, Colleen , Kendrick, Samantha , Barrett, Michael T. , Glinsmann‐Gibson, Betty, McGrath, Michael S. , Maguire, Alanna , Rimsza, Lisa M.

International Journal of Cancer | May 2, 2019 | Volume: 145 | Issue: 11 | Pages: 3078-3088

10.1002/ijc.32381 Get It

Diffuse large B‐cell lymphoma (DLBCL) is up to 17‐fold more likely to occur, follows a more aggressive clinical course and frequently presents at advanced stages in HIV infected (+) individuals compared to HIV negative (−) individuals. However, the molecular pathology underpinning the clinical features of DLBCL in HIV(+) patients relative to the general population is poorly understood. We performed a retrospective study examining the transcriptional, genomic and protein expression differences between HIV(+) and HIV(−) germinal center B‐cell (GCB) DLBCL cases using digital gene expression analysis, array comparative genomic hybridization (CGH) and immunohistochemistry (IHC). Genes associated with cell cycle progression ( CCNA2 , CCNB1 , CDC25A , E2F1 ), DNA replication ( MCM2 , MCM4 , MCM7 ) and DNA damage repair, including eight Fanconi anemia genes ( FANCA , FANCD1/BRCA2 , FANCE , FANCG , FANCR/RAD51 , FANCS/BRCA1 , FANCT/UBE2T , FANCV/MAD2L2 ), were significantly increased in HIV(+) GCB‐DLBCL tumors compared to HIV(−) tumors. In contrast, genes associated with cell cycle inhibition ( CDKN1A , CDKN1B ) as well as apoptosis regulating BCL2 family members ( BCL2 , BAX , BIM , BMF , PUMA ) were significantly decreased in the HIV(+) cohort. BCL2 IHC confirmed this expression. Array CGH data revealed that HIV(+) GCB‐DLBCL tumors have fewer copy number variations than their HIV(−) counterparts, indicating enhanced genomic stability. Together, the results show that HIV(+) GCB‐DLBCL is a distinct molecular malignancy from HIV(−) GCB‐DLBCL; with an increased proliferative capacity, confirmed by Ki67 IHC staining, and enhanced genomic stability, the latter of which is likely related to the enhanced expression of DNA repair genes.
An exhaustive algorithm for detecting copy number aberrations and large structural variants in whole-genome, mate-pair sequencing data.
Asmann, Yan W. , Necela, Brian M. , Chen, Xianfeng , Kocher, Jean‐Pierre , Feldman, Andew L., Novak, Anne J. , Perez, Edith A. , Thompson, E. Aubrey , Wang, Chen , Maurer, Matthew J. , Habermann, Thomas M. , Slager, Susan L. , Doǧan, Ahmet , Cerhan, James R.

Journal of Clinical Oncology | January 3, 2019 | Volume: 32 | Issue: 15_suppl | Pages: e22171-e22171

10.1200/jco.2014.32.15_suppl.e22171 Get It

Loss of SETD2 Induces a Metabolic Switch in Renal Cell Carcinoma Cell Lines toward Enhanced Oxidative Phosphorylation
Liu, Jingping , Hanavan, Paul D., Kras, Katon A., Ruiz, Yvette , Castle, Erik P., Lake, Douglas F. , O’Brien, Daniel R. , Luo, Huijun , Robertson, Keith D. , Gu, Haiwei , Ho, Thai H. , Chen, Xianfeng , Gu, Haiwei , Ho, Thai , Robertson, Keith D , Chen, Xianfeng

Journal of Proteome Research | November 8, 2018 | Volume: No volume available | Issue: No issue available | Pages: No page available

10.1021/acs.jproteome.8b00628 Get It

Hybrid-denovo: a de novo OTU-picking pipeline integrating single-end and paired-end 16S sequence tags
Chen, Xianfeng , Johnson, Stephen , Wang, Junwen , Kocher, Jean‐Pierre , Jeraldo, Patricio , Chia, Nicholas , Wang, Junwen , Chia, Nicholas , Chen, Jun

GigaScience | December 15, 2017 | Volume: 7 | Issue: 3 | Pages: No page available

10.1093/gigascience/gix129 Get It

Abstract Background Illumina paired-end sequencing has been increasingly popular for 16S rRNA gene-based microbiota profiling. It provides higher phylogenetic resolution than single-end reads due to a longer read length. However, the reverse read (R2) often has significant low base quality, and a large proportion of R2s will be discarded after quality control, resulting in a mixture of paired-end and single-end reads. A typical 16S analysis pipeline usually processes either paired-end or single-end reads but not a mixture. Thus, the quantification accuracy and statistical power will be reduced due to the loss of a large amount of reads. As a result, rare taxa may not be detectable with the paired-end approach, or low taxonomic resolution will result in a single-end approach. Results To have both the higher phylogenetic resolution provided by paired-end reads and the higher sequence coverage by single-end reads, we propose a novel OTU-picking pipeline, hybrid-denovo, that can process a hybrid of single-end and paired-end reads. Using high-quality paired-end reads as a gold standard, we show that hybrid-denovo achieved the highest correlation with the gold standard and performed better than the approaches based on paired-end or single-end reads in terms of quantifying the microbial diversity and taxonomic abundances. By applying our method to a rheumatoid arthritis (RA) data set, we demonstrated that hybrid-denovo captured more microbial diversity and identified more RA-associated taxa than a paired-end or single-end approach. Conclusions Hybrid-denovo utilizes both paired-end and single-end 16S sequencing reads and is recommended for 16S rRNA gene targeted paired-end sequencing data.