scholarly journals BraCeR: B-cell-receptor reconstruction and clonality inference from single-cell RNA-seq

2018 ◽  
Vol 15 (8) ◽  
pp. 563-565 ◽  
Author(s):  
Ida Lindeman ◽  
Guy Emerton ◽  
Lira Mamanova ◽  
Omri Snir ◽  
Krzysztof Polanski ◽  
...  
2017 ◽  
Author(s):  
Simone Rizzetto ◽  
David NP Koppstein ◽  
Jerome Samir ◽  
Mandeep Singh ◽  
Joanne H. Reed ◽  
...  

AbstractThe B-cell receptor (BCR) performs essential functions for the adaptive immune system including recognition of pathogen-derived antigens. Cell-to-cell variability of BCR sequences due to V(D)J recombination and somatic hypermutation (SHM) necessitates single-cell characterization of BCR sequences. Single-cell RNA sequencing (scRNA-seq) presents the opportunity for simultaneous capture of the BCR sequence and transcriptomic signature for a detailed understanding of the dynamics of an immune response.We developed VDJPuzzle 2.0, a bioinformatic tool that reconstructs productive, full-length B-cell receptor sequences of both heavy and light chains. VDJPuzzle successfully reconstructs BCRs from 98.3% (n=117) of human and 96.5% (n=200) from murine B cells. 92.0% of clonotypes and 90.3% of mutations were concordant with single-cell Sanger sequencing of the immunoglobulin chains. VDJPuzzle is available at https://bitbucket.org/kirbyvisp/vdjpuzzle2


2018 ◽  
Vol 34 (16) ◽  
pp. 2846-2847 ◽  
Author(s):  
Simone Rizzetto ◽  
David N P Koppstein ◽  
Jerome Samir ◽  
Mandeep Singh ◽  
Joanne H Reed ◽  
...  

2021 ◽  
Vol 12 ◽  
Author(s):  
Han Sun ◽  
Hu-Qin Yang ◽  
Kan Zhai ◽  
Zhao-Hui Tong

B cells play vital roles in host defense against Pneumocystis infection. However, the features of the B cell receptor (BCR) repertoire in disease progression remain unclear. Here, we integrated single-cell RNA sequencing and single-cell BCR sequencing of immune cells from mouse lungs in an uninfected state and 1–4 weeks post-infection in order to illustrate the dynamic nature of B cell responses during Pneumocystis infection. We identified continuously increased plasma cells and an elevated ratio of (IgA + IgG) to (IgD + IgM) after infection. Moreover, Pneumocystis infection was associated with an increasing naïve B subset characterized by elevated expression of the transcription factor ATF3. The proportion of clonal expanded cells progressively increased, while BCR diversity decreased. Plasma cells exhibited higher levels of somatic hypermutation than naïve B cells. Biased usage of V(D)J genes was observed, and the usage frequency of IGHV9-3 rose. Overall, these results present a detailed atlas of B cell transcriptional changes and BCR repertoire features in the context of Pneumocystis infection, which provides valuable information for finding diagnostic biomarkers and developing potential immunotherapeutic targets.


2017 ◽  
Vol 199 (2) ◽  
pp. 782-791 ◽  
Author(s):  
Bishnudeo Roy ◽  
Ralf S. Neumann ◽  
Omri Snir ◽  
Rasmus Iversen ◽  
Geir Kjetil Sandve ◽  
...  

Blood ◽  
2009 ◽  
Vol 114 (22) ◽  
pp. 1263-1263
Author(s):  
Erik Evensen ◽  
Adam Palazzo ◽  
Ying-Wen Huang ◽  
Alessandra Cesano ◽  
Laura Z. Rassenti ◽  
...  

Abstract Abstract 1263 Poster Board I-285 Background In conjunction with antigen-driven responses, ligand-independent signaling (termed tonic signaling) through both the pre-B cell receptor and B-cell receptor has an important role in B cell development, maturation and survival. In addition to the recognized role of CD79 alpha and CD79 beta BCR signaling, tyrosine phosphatases can impact tonic BCR signaling (Wienands et al. PNAS, 93 p.7865 (1996), Monroe Nat. Rev. Immunol. 6 p.283 (2006)). We previously subjected chronic lymphocytic leukemia (CLL) cells with modulators of BCR signaling and monitored their responses using flow cytometry-based Single Cell Network Profiling (SCNP). Of the many signaling modulators studied, hydrogen peroxide treatment (a general inhibitor of tyrosine phosphatase activity) augmented BCR signaling in a subset of CLL patient samples evaluated. In the remaining samples there was an apparent lack of response to hydrogen peroxide. These data suggested that differential phosphatase activity proximal to BCR signaling was driving the biology of these two patient groups. Objectives Studies were designed to evaluate whether there were any associations between tonic and/or ligand-dependent BCR signaling and in vitro sensitivity to fludarabine, as well as whether such response profiles showed a relationship to the hydrogen peroxide-dependent signaling we observed previously. Methods 23 CLL samples and 7 healthy PBMCs were treated with anti-m alone, hydrogen peroxide alone or the combination for 10 minutes. Separate aliquots of the same sample were exposed to F-ara-A for 48 hours. SCNP was carried out on gated B cells with quantitation of single cell measures of intracellular phosphorylated kinases and adaptor proteins downstream of the BCR. Additionally, the relative activation status of several protein markers of the apoptotic cascade (cytoplasmic cytochrome C, cleaved caspase 3, and cleaved PARP) was measured. Results As previously observed, CLL samples could be segregated into one of two groups exhibiting either responsive or refractory signaling after exposure to hydrogen peroxide alone. Moreover, responsive signaling in CLL cells was correlated in that all the measured components of the canonical B cell receptor network (p-Lyn, p-Syk, p-BLNK, p-PLC-gamma-2, p-Erk and p-Akt) showed the same phosphorylation response: either augmented in unison, or not activated at all. In vitro F-ara-A treatment (48 hours in the presence of 1mM F-ara-A) of parallel samples from these same CLL patients identified distinct populations of apoptosis responsive and refractory cells. Surprisingly, the capacity of patient samples to show augmented BCR signaling in response to hydrogen peroxide was associated prominently with the ability of cells in these patients to exhibit apoptotic proficiency to F-ara-A in vitro. This implies a link between mechanisms governing apoptosis in these CLL cells, survival pathways, and cell states that govern the role of phosphatase activity and BCR signaling potential. Conclusions This study reveals a link between tonic BCR signaling and regulation of apoptosis pathways. This suggests that the subgroup of CLL patients with active phosphatase activity (which suppresses BCR responses) have cell populations that are responsive to F-ara-A, a standard drug in CLL therapy. Conversely, the presence of CLL cells in a patient sample that remain unresponsive to hydrogen peroxide repression of phosphatase activity appear to identify patient samples which cannot undergo apoptosis in response to in vitro F-Ara-A exposure. The clinical implications of this work will be the focus of future translational studies. Disclosures Evensen: Nodality Inc.: Employment, Equity Ownership. Palazzo:Nodality Inc.: Employment, Equity Ownership. Huang:Nodality Inc.: Employment, Equity Ownership. Cesano:Nodality Inc.: Employment, Equity Ownership. Fantl:Nodality, Inc.: Employment, Equity Ownership.


2020 ◽  
Author(s):  
Anna M. Nia ◽  
Kamil Khanipov ◽  
Brooke L. Barnette ◽  
Robert L. Ullrich ◽  
George Golovko ◽  
...  

Abstract Background: One of the health risks posed to astronauts during deep space flights is exposure to high charge, high-energy (HZE) ions (Z>13), which can lead to induction of hepatocellular carcinoma (HCC). However, little is known on the molecular mechanisms of HZE irradiation induced HCC.Results: We performed comparative RNA-Seq transcriptomic analyses to assess the carcinogenic effects of 600 MeV/n 56 Fe (0.2 Gy), 1 GeV/n 16 O (0.2 Gy), and 350 MeV/n 28 Si (0.2 Gy) ions in a mouse model for irradiation-induced HCC. C3H/HeNCrl mice were subjected to total body irradiation to simulate space environment HZE-irradiation, and liver tissues were extracted at five different time points post-irradiation to investigate the time-dependent carcinogenic response at the transcriptomic level. Our data demonstrated a clear difference in the biological effects of these HZE ions, particularly immunological, such as Acute Phase Response Signaling, B Cell Receptor Signaling, IL-8 Signaling, and ROS Production in Macrophages. Also seen in this study were novel unannotated transcripts that were significantly affected by HZE. To investigate the biological functions of these novel transcripts, we used a machine learning technique known as self-organizing maps (SOMs) to characterize the transcriptome expression profiles of 60 samples (45 HZE-irradiated, 15 non-irradiated control) from liver tissues. A handful of localized modules in the maps emerged as groups of co-regulated and co-expressed transcripts. The functional context of these modules was discovered using overrepresentation analysis. We found that these spots typically contained enriched populations of transcripts related to specific immunological molecular processes (e.g., Acute Phase Response Signaling, B Cell Receptor Signaling, IL-3 Signaling), and RNA Transcription/Expression.Conclusions: A large number of transcripts were found differentially expressed post-HZE irradiation. These results provide valuable information for uncovering the differences in molecular mechanisms underlying HZE specific induced HCC carcinogenesis. Additionally, a handful of novel differentially expressed unannotated transcripts were discovered for each HZE ion. Taken together, these findings may provide a better understanding of biological mechanisms underlying risks for HCC after HZE irradiation, and may also have important implications for discovery of potential countermeasures against and identification of biomarkers for HZE-induced HCC.


Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 41-41
Author(s):  
Krystle Nomie ◽  
Liang Zhang ◽  
Yixin Yao ◽  
Yang Liu ◽  
Shaojun Zhang ◽  
...  

Abstract Introduction Mantle cell lymphoma (MCL) is an incurable B-cell lymphoma subtype and constitutive activation of the B-cell receptor pathway is a hallmark of B-cell lymphomas. Bruton's tyrosine kinase (BTK) is a critical component of the B-cell receptor pathway, and ibrutinib, a first-in-class, once-daily, and oral covalent inhibitor of BTK, was developed to reduce/silence B-cell receptor pathway activity, leading to clinically remarkable anti-tumor activity. In our prior multiple-center Phase II clinical trial, the overall response rate in relapsed/refractory MCL patients was 68% (Wang et al., NEJM, 2013), surpassing the effectiveness of other therapies. Although ibrutinib is extremely efficacious in patients with relapsed/refractory MCL, the one-year overall survival rate of ibrutinib-exposed patients who relapse is only 22%. Methods Patient primary cells were isolated from MCL patients treated with ibrutinib either prior to treatment or at treatment discontinuation. Whole exome sequencing (WES) was performed to determine the mutational landscape of ibrutinib resistance. RNA-seq was employed to compare the gene expression profiles between ibrutinib-sensitive and -resistant patient samples. Gene set enrichment analysis was utilized to identify dysregulated molecular pathways associated with the resistant phenotype. The RNA-seq data were then validated with reverse phase protein array (RPPA) analysis of ibrutinib-sensitive and -resistant MCL cell lines. Metabolic assays including the measurement of mitochondria respiration rates with the Seahorse analyzer and reactive oxygen species (ROS) levels, targeted metabolomics, and ATP analysis. Functional studies targeting this molecular pathway were conducted, including in vitro cell viability and apoptosis assays, as well as in vivo efficacy studies in an ibrutinib-resistant MCL patient-derived xenograft mouse model. Results WES data analysis identified frequent inactivating somatic alterations in ATM, KMT2D, and TP53 in both the ibrutinib-sensitive and -resistant tumors. CDKN2A (5/7, 71%) was frequently deleted, and the deletion was only observed in the ibrutinib-resistant tumors (p = 0.010). The RNA-seq analysis identified a total of 63 protein-coding genes as the most differentially expressed genes (DEGs) between the ibrutinib-resistant and -sensitive groups, with a fold change of ≥ 2 or ≤ -2 and the false discovery rate (FDR q-value) ≤ 0.01. Among the DEGs, 26 genes were upregulated in ibrutinib-resistant tumors. In addition, gene set enrichment analysis (GSEA) revealed the marked upregulation of oncogenic pathways including c-MYC, mTOR (mTORC1), Wnt, and NF-ĸb signaling, followed by cell cycle, apoptosis, BCR signaling and DNA repair in the ibrutinib-resistant tumors. Notably, in addition to these oncogenic pathways, the metabolic pathways, including oxidative phosphorylation (OXPHOS), were significantly enriched in the ibrutinib-resistant tumors (normalized enrichment score > 3 and FDR q-value < 1e-5). In further support of this finding, metabolomics analysis and the measurement of ATP production and mitochondrial respiration indicated that the OXPHOS pathway is the predominant metabolic pathway employed by ibrutinib-resistant MCL cells. To determine the effects of targeting these pathways, OXPHOS was inhibited with a novel electron transport complex I inhibitor (IACS-010759, developed by MD Anderson Cancer Center) in both MCL cell lines and ibrutinib-resistant MCL patient-derived xenograft (PDX) models. Single agent IACS-010759 treatment at 10 mg/kg oral gavage for 5 consecutive days/week completely prevented tumor growth compared with the vehicle control as shown by measuring tumor volume (n = 5, p < 0.0001) and human β2M levels (n = 5; p < 0.0001) throughout treatment. No apparent toxicities were observed in the IACS-010759-treated MCL PDX mice. Conclusion This current study warrants the exploitation of active cancer metabolic pathways, especially OXPHOS, to improve the clinical outcomes of MCL and additional lymphoma, which is actively being investigated in a Phase I lymphoma clinical trial (NCT03291938). Disclosures Wang: AstraZeneca: Consultancy, Research Funding; Celgene: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; Juno: Research Funding; Novartis: Research Funding; Janssen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding; MoreHealth: Consultancy; Acerta Pharma: Honoraria, Research Funding; Kite Pharma: Research Funding; Pharmacyclics: Honoraria, Research Funding; Dava Oncology: Honoraria.


Author(s):  
Wen Wen ◽  
Wenru Su ◽  
Hao Tang ◽  
Wenqing Le ◽  
Xiaopeng Zhang ◽  
...  

AbstractCOVID-19, caused by SARS-CoV-2, has recently affected over 300,000 people and killed more than 10,000. The manner in which the key immune cell subsets change and their states during the course of COVID-19 remain unclear. Here, we applied single-cell technology to comprehensively characterize transcriptional changes in peripheral blood mononuclear cells during the recovery stage of COVID-19. Compared with healthy controls, in patients in the early recovery stage (ERS) of COVID-19, T cells decreased remarkably, whereas monocytes increased. A detailed analysis of the monocytes revealed that there was an increased ratio of classical CD14++ monocytes with high inflammatory gene expression as well as a greater abundance of CD14++IL1B+ monocytes in the ERS. CD4+ and CD8+ T cells decreased significantly and expressed high levels of inflammatory genes in the ERS. Among the B cells, the plasma cells increased remarkably, whereas the naïve B cells decreased. Our study identified several novel B cell-receptor (BCR) changes, such as IGHV3-23 and IGHV3-7, and confirmed isotypes (IGHV3-15, IGHV3-30, and IGKV3-11) previously used for virus vaccine development. The strongest pairing frequencies, IGHV3-23-IGHJ4, indicated a monoclonal state associated with SARS-CoV-2 specificity. Furthermore, integrated analysis predicted that IL-1β and M-CSF may be novel candidate target genes for inflammatory storm and that TNFSF13, IL-18, IL-2 and IL-4 may be beneficial for the recovery of COVID-19 patients. Our study provides the first evidence of an inflammatory immune signature in the ERS, suggesting that COVID-19 patients are still vulnerable after hospital discharge. Our identification of novel BCR signaling may lead to the development of vaccines and antibodies for the treatment of COVID-19.Highlights-The immune response was sustained for more than 7 days in the early recovery stage of COVID-19, suggesting that COVID-19 patients are still vulnerable after hospital discharge.-Single-cell analysis revealed a predominant subset of CD14++ IL1β+ monocytes in patients in the ERS of COVID-19.-Newly identified virus-specific B cell-receptor changes, such as IGHV3-23, IGHV3-7, IGHV3-15, IGHV3-30, and IGKV3-11, could be helpful in the development of vaccines and antibodies against SARS-CoV-2.-IL-1β and M-CSF were discovered as novel mediators of inflammatory cytokine storm, and TNFSF13, IL-2, IL-4, and IL-18 may be beneficial for recovery.


Blood ◽  
2019 ◽  
Vol 134 (Supplement_1) ◽  
pp. 4275-4275
Author(s):  
Daniel Mertens ◽  
Christine Wolf ◽  
Carsten Maus ◽  
Michael Persicke ◽  
Katharina Filarsky ◽  
...  

B-cell receptor (BCR) signalling is central for the pathomechanism of chronic lymphocytic leukemia (CLL). Novel inhibitors of BCR signalling have recently substantially improved treatment of CLL, and a better characterization of the molecular circuitry of leukemic BCR signalling will allow a more refined targeting of this Achilles heel. In order to model malignant and non-malignant BCR signalling, we quantified after stimulation 5 components of BCR signaling (ZAP70/SYK, BTK, PLCy2, AKT, ERK1/2) in single cells from primary human leukemic and non-malignant tissue via phospho-specific flow cytometry over 6 time points. We stimulated cells from 11 patients and non-malignant CD19 negative enriched B-cells from 5 healthy donors by crosslinking the BCR with anti-IgM and/or anti-CD19 and synchronous inhibition of phosphatases with H2O2. As expected, we found more phosphorylation of all BCR signalling components after stimulation in malignant vs non-malignant cells and in IGHV non-mutated CLL cells compared to IGHV mutated CLL cells. Intriguingly, inhibition of phosphatases with H2O2 led to higher phosphorylation of BCR components in CLL cells with mutated IGHV genes compared to CLL cells with non-mutated IGHV genes, suggesting a stronger dampening of signalling activity in mutated IGHV CLL by phosphatases. In order to characterize the signalling circuitry, we modelled the connectivity of the cascade components by correlating signal intensities across single cells of the cell populations of single samples (Figure 1). Surprisingly, upon stimulation no substantial differences in network topology were observed between malignant and non-malignant cells. To additionally test for changes in network topology, we challenged the BCR signaling cascade with inhibitors for BTK (ibrutinib), PI3K (idelalisib). Ibrutinib and idelalisib acted complementary, but not synergistic, and were similarly effective in IGHV mutated and non-mutated CLL. Effects of idelalisib were the same on malignant and non-malignant cells, whereas ibrutinib was mostly active on CLL cells, not on non-malignant B-cells. Upon stimulation with combinations of IgM and CD19 crosslinking augmented with H2O2, phosphorylation of PLCy2 could not be significantly inhibited by idelalisib or ibrutinib on a timescale of 28mins. We therefore aimed to identify central activating nodes of the BCR signalling cascade using targeted inhibitors. In fact, we found that inhibition of LYN with dasatinib and inhibition of SYK with entospletinib could substantially reduce phosphorylation of PLCy2, BTK and ERK but not AKT after all combinations of BCR stimulation. This suggests additional signalling cascades modulating AKT and a strong impact of SYK/LYN activity on the regulation of PLCy2. In summary, our findings underline the importance of single cell analysis of the dynamic circuitry of B-cell receptor signalling to understand development of resistance mechanisms and potential vulnerabilities. Figure 1: Workflow scheme of the Bayesian network learning and averaging approach. After discretizing the continuous single cell data, an optimal network is derived from each of R bootstrap samples. The Bayesian network learning strategy uses the BDe scoring function and a greedy hill-climbing algorithm to find the network model that represents the resampled data best. An average arc strength for each connection between nodes is derived from the number of occurrences of the respective connection in the set of R best scoring networks. Further averaging among networks derived from different data sets was applied for identifying conditional, temporal, and group-specific differences. Figure 1 Disclosures Döhner: Novartis: Consultancy, Honoraria, Research Funding; Astex: Consultancy, Honoraria; Bristol Myers Swuibb: Research Funding; Amgen: Consultancy, Honoraria, Research Funding; Arog: Research Funding; Seattle Genetics: Consultancy, Honoraria; Astellas: Consultancy, Honoraria; Roche: Consultancy, Honoraria; Janssen: Consultancy, Honoraria; Agios: Consultancy, Honoraria; Pfizer: Research Funding; Celgene Corporation: Consultancy, Honoraria, Research Funding; Jazz: Consultancy, Honoraria, Research Funding; AbbVie: Consultancy, Honoraria. Stilgenbauer:GSK: Consultancy, Honoraria, Research Funding, Speakers Bureau; Janssen: Consultancy, Honoraria, Research Funding, Speakers Bureau; Pharmacyclics: Other: Travel support; Amgen: Consultancy, Honoraria, Research Funding, Speakers Bureau; Novartis: Consultancy, Honoraria, Research Funding, Speakers Bureau; Gilead: Consultancy, Honoraria, Research Funding, Speakers Bureau; AstraZeneca: Consultancy, Honoraria, Research Funding, Speakers Bureau; Celgene: Consultancy, Honoraria, Research Funding, Speakers Bureau; Hoffmann La-Roche: Consultancy, Honoraria, Research Funding, Speakers Bureau; AbbVie: Consultancy, Honoraria, Research Funding, Speakers Bureau.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. SCI-26-SCI-26
Author(s):  
Louis M. Staudt

Abstract Abstract SCI-26 We have developed loss-of-function, RNA interference-based screens to reveal genes essential for cancer cell proliferation and survival. In parallel, we are using high-throughput RNA resequencing (RNA-seq) to identify somatic mutations and other structural abnormalities in cancer. The intersection of these two data sets has helped us to discover novel pathogenetic pathways in lymphoma that are amenable to therapeutic attack. The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) has constitutive activation of the NF-κB pathway, which we traced to the signaling adapter CARD11. In some ABC DLBCL biopsies (∼10%), somatic mutations produce CARD11 isoforms that spontaneously activate NF-κB signaling. In ABC DLBCL tumors with wild-type CARD11, we defined a “chronic active” form of B-cell receptor (BCR) signaling that activates NF-κB. Such ABC DLBCLs are killed by knockdown of BCR signaling components, such as Bruton's tyrosine kinase (BTK), or components of the BCR itself. Over one-fifth of ABC DLBCLs have mutations in the CD79B or CD79A subunits of the BCR. In 18 percent of cases, mutations occur in a single tyrosine residue in the critical “ITAM” signaling motif, generating BCRs that avoid negative autoregulation by the LYN tyrosine kinase. Based on these findings, we are conducting clinical trials of ibrutinib in relapsed/refractory DLBCL. Ibrutinib is an irreversible and highly selective small-molecule inhibitor of BTK. Thus far, ibrutinib monotherapy has induced many complete and partial responses in patients with ABC DLBCL, including those with “primary refractory” tumors that had never responded to any prior therapy. One patient has been in a sustained complete response for over 19 months, taking ibrutinib daily with no discernible side effects. Of note, ABC DLBCL tumors with and without CD79B mutations have responded, suggesting that BCR pathway addiction may be a prevalent feature in this lymphoma subtype. More recently, we have uncovered a “tonic” form of BCR signaling in Burkitt lymphoma that engages the prosurvival PI(3) kinase pathway. Two-thirds of Burkitt lymphoma cell lines die upon knockdown of BCR subunits or the proximal kinase SYK, due to loss of PI(3) kinase signaling. Moreover, a gene expression signature of PI(3) kinase activity is more highly expressed in Burkitt lymphoma biopsies than in biopsies of other aggressive lymphomas. Tonic BCR signaling in Burkitt lymphoma is mechanistically distinct from chronic active BCR signaling in ABC DLBCL, since it does not engage BTK, CARD11, or NF-κB. RNA-seq revealed that 70 percent of sporadic Burkitt lymphoma cases harbor somatic mutations that potentiate the action of the transcription factor TCF3 by preventing its inhibitory heterodimerization with the DNA-binding inhibitor ID3. TCF3 promotes tonic BCR signaling and PI(3) kinase activity in Burkitt lymphoma by transactivating the immunoglobulin heavy- and light-chain genes, thereby increasing surface BCR expression, and by repressing the phosphatase SHP-1, a potent negative regulator of BCR signaling. Hence, inhibitors of proximal BCR signaling and the PI(3) kinase pathway should be evaluated in Burkitt lymphoma, especially in patients for whom high-dose chemotherapy is infeasible, such as older individuals and those with the endemic form of this lymphoma. Disclosures: Off Label Use: I will be discussing clinical trials of ibrutinib (PCI-32765) in lymphoma. Ibrutinib is an investigational drug that has not yet received FDA approval for any indication.


Sign in / Sign up

Export Citation Format

Share Document