Characterization of Antibody Repertoires in Patients with HPV-Related HNSCC Undergoing Definitive Radiation with Immunotherapy

PURPOSE/OBJECTIVE(S)

Integrating immune checkpoint inhibitor (ICI) immunotherapy agents into the definitive treatment of head and neck squamous cell carcinoma (HNSCC) has been challenging. It is apparent that understanding how these therapies affect specific components of the immune response is critical to improve outcomes. The role of B cells has been increasingly recognized, especially in HPV+ HNSCC. As antibodies are one of the chief downstream products of B cells, we sought to evaluate the antibody repertoires of patients with HPV+ HNSCC undergoing definitive radiation therapy (RT) and ICI (RT-ICI).

MATERIALS/METHODS

Serum samples from patients with p16+ HNSCC undergoing RT-ICI were collected for the KEYCHAIN clinical trial (NCT03383094). We analyzed 8 samples from 4 patients collected pre-treatment and 3 months following treatment using the HuProt Human Proteome Microarray, which samples >20,000 human proteins. Genes encoding surface proteins (GESPs) were obtained from The Cancer Surfaceome Atlas. Mutational data were obtained from AACR Project GENIE. For murine studies, C3H mice were injected with 5x105 cells of AT-84-E7-OVA syngeneic SCC into the right flank. Anti-PD-L1 ICI (aPD-L1) was given via intraperitoneal injection every 3 days for 3 doses. RT was 12 Gy in 1 fraction. Cell surface markers in lymph nodes (LNs) were analyzed by flow cytometry.

RESULTS

In total, we detected antibodies against 10959 unique antigens in the pre-treatment serum, of which 14% and 11% were shared by 3 and 4 patients, respectively. Following completion of RT-ICI, we detected antibodies against 11019 unique antigens, of which 20% and 18% were shared by 3 and 4 patients, respectively. Of these, 5824 (53%) antigens were not detected in the pre-treatment serum, and therefore represent antibodies against "newly detected" antigens. We found 777 antigens that were newly detected and shared by at least 2 patients. Analysis of these antigens revealed enrichment in pathways such as, "response to oxygen levels." We next found that 114 (14.7%) of these represented GESPs. We integrated mutational analysis of the most frequently mutated genes in >1,800 HNSCC samples and found 2/114 GESPs were shared. Using a murine model of HPV+ HNSCC, we found that treatment with RT and aPD-L1 led to the greatest frequency of germinal center (GC) B cells in tumor-draining and non-tumor-draining LNs.

CONCLUSION

In patients with p16-positive HNSCC, proteomic analysis of antibody repertoires revealed many antigens against which antibodies were formed during RT-ICI that are shared between patients. Intriguingly, GC formation, which is the nidus for B cell responses, in locoregional LNs was greatest with RT-ICI. These findings support the further investigation of B-cell mediated responses in HPV+ HNSCC.

LysM-mediated signaling in Marchantia polymorpha highlights the conservation of pattern-triggered immunity in land plants

Pattern-recognition receptor (PRR)-triggered immunity (PTI) wards off a wide range of pathogenic microbes, playing a pivotal role in angiosperms. The model liverwort Marchantia polymorpha triggers defense-related gene expression upon sensing components of bacterial and fungal extracts, suggesting the existence of PTI in this plant model. However, the molecular components of the putative PTI in M. polymorpha and the significance of PTI in bryophytes have not yet been described. We here show that M. polymorpha has four lysin motif (LysM)-domain-containing receptor homologs, two of which, LysM-receptor-like kinase (LYK) MpLYK1 and LYK-related (LYR) MpLYR, are responsible for sensing chitin and peptidoglycan fragments, triggering a series of characteristic immune responses. Comprehensive phosphoproteomic analysis of M. polymorpha in response to chitin treatment identified regulatory proteins that potentially shape LysM-mediated PTI. The identified proteins included homologs of well-described PTI components in angiosperms as well as proteins whose roles in PTI are not yet determined, including the blue-light receptor phototropin MpPHOT. We revealed that MpPHOT is required for negative feedback of defense-related gene expression during PTI. Taken together, this study outlines the basic framework of LysM-mediated PTI in M. polymorpha and highlights conserved elements and new aspects of pattern-triggered immunity in land plants.

Cofunctioning of bacterial exometabolites drives root microbiota establishment

Soil-dwelling microbes are the principal inoculum for the root microbiota, but our understanding of microbe-microbe interactions in microbiota establishment remains fragmentary. We tested 39,204 binary interbacterial interactions for inhibitory activities in vitro, allowing us to identify taxonomic signatures in bacterial inhibition profiles. Using genetic and metabolomic approaches, we identified the antimicrobial 2,4-diacetylphloroglucinol (DAPG) and the iron chelator pyoverdine as exometabolites whose combined functions explain most of the inhibitory activity of the strongly antagonistic Pseudomonas brassicacearum R401. Microbiota reconstitution with a core of Arabidopsis thaliana root commensals in the presence of wild-type or mutant strains revealed a root niche-specific cofunction of these exometabolites as root competence determinants and drivers of predictable changes in the root-associated community. In natural environments, both the corresponding biosynthetic operons are enriched in roots, a pattern likely linked to their role as iron sinks, indicating that these cofunctioning exometabolites are adaptive traits contributing to pseudomonad pervasiveness throughout the root microbiota.

Functional genomics gives new insights into the ectomycorrhizal degradation of chitin

Ectomycorrhizal (EcM) fungi play a crucial role in the mineral nitrogen (N) nutrition of their host trees. While it has been proposed that several EcM species also mobilize organic N, studies reporting the EcM ability to degrade N-containing polymers, such as chitin, remain scarce. Here, we assessed the capacity of a representative collection of 16 EcM species to acquire ¹⁵ N from ¹⁵ N-chitin. In addition, we combined genomics and transcriptomics to identify pathways involved in exogenous chitin degradation between these fungal strains. Boletus edulis, Imleria badia, Suillus luteus, and Hebeloma cylindrosporum efficiently mobilized N from exogenous chitin. EcM genomes primarily contained genes encoding for the direct hydrolysis of chitin. Further, we found a significant relationship between the capacity of EcM fungi to assimilate organic N from chitin and their genomic and transcriptomic potentials for chitin degradation. These findings demonstrate that certain EcM fungal species depolymerize chitin using hydrolytic mechanisms and that endochitinases, but not exochitinases, represent the enzymatic bottleneck of chitin degradation. Finally, this study shows that the degradation of exogenous chitin by EcM fungi might be a key functional trait of nutrient cycling in forests dominated by EcM fungi.

Speciation Underpinned by Unexpected Molecular Diversity in the Mycorrhizal Fungal Genus Pisolithus

The mutualistic ectomycorrhizal (ECM) fungal genus Pisolithus comprises 19 species defined to date which colonize the roots of >50 hosts worldwide suggesting that substantial genomic and functional evolution occurred during speciation. To better understand this intra-genus variation, we undertook a comparative multi-omic study of nine Pisolithus species sampled from North America, South America, Asia, and Australasia. We found that there was a small core set of genes common to all species (13%), and that these genes were more likely to be significantly regulated during symbiosis with a host than accessory or species-specific genes. Thus, the genetic "toolbox" foundational to the symbiotic lifestyle in this genus is small. Transposable elements were located significantly closer to gene classes including effector-like small secreted proteins (SSPs). Poorly conserved SSPs were more likely to be induced by symbiosis, suggesting that they may be a class of protein that tune host specificity. The Pisolithus gene repertoire is characterized by divergent CAZyme profiles when compared with other fungi, both symbiotic and saprotrophic. This was driven by differences in enzymes associated with symbiotic sugar processing, although metabolomic analysis suggest that neither copy number nor expression of these genes is sufficient to predict sugar capture from a host plant or its metabolism in fungal hyphae. Our results demonstrate that intra-genus genomic and functional diversity within ECM fungi is greater than previously thought, underlining the importance of continued comparative studies within the fungal tree of life to refine our focus on pathways and evolutionary processes foundational to this symbiotic lifestyle.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

Phylogenomics and Comparative Genomics Highlight Specific Genetic Features in Ganoderma Species

The Ganoderma species in Polyporales are ecologically and economically relevant wood decayers used in traditional medicine, but their genomic traits are still poorly documented. In the present study, we carried out a phylogenomic and comparative genomic analyses to better understand the genetic blueprint of this fungal lineage. We investigated seven Ganoderma genomes, including three new genomes, G. australe, G. leucocontextum, and G. lingzhi. The size of the newly sequenced genomes ranged from 60.34 to 84.27 Mb and they encoded 15,007 to 20,460 genes. A total of 58 species, including 40 white-rot fungi, 11 brown-rot fungi, four ectomycorrhizal fungi, one endophyte fungus, and two pathogens in Basidiomycota, were used for phylogenomic analyses based on 143 single-copy genes. It confirmed that Ganoderma species belong to the core polyporoid clade. Comparing to the other selected species, the genomes of the Ganoderma species encoded a larger set of genes involved in terpene metabolism and coding for secreted proteins (CAZymes, lipases, proteases and SSPs). Of note, G. australe has the largest genome size with no obvious genome wide duplication, but showed transposable elements (TEs) expansion and the largest set of terpene gene clusters, suggesting a high ability to produce terpenoids for medicinal treatment. G. australe also encoded the largest set of proteins containing domains for cytochrome P450s, heterokaryon incompatibility and major facilitator families. Besides, the size of G. australe secretome is the largest, including CAZymes (AA9, GH18, A01A), proteases G01, and lipases GGGX, which may enhance the catabolism of cell wall carbohydrates, proteins, and fats during hosts colonization. The current genomic resource will be used to develop further biotechnology and medicinal applications, together with ecological studies of the Ganoderma species.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

The fungal root endophyte Serendipita vermifera displays inter-kingdom synergistic beneficial effects with the microbiota in Arabidopsis thaliana and barley

Plant root-associated bacteria can confer protection against pathogen infection. By contrast, the beneficial effects of root endophytic fungi and their synergistic interactions with bacteria remain poorly defined. We demonstrate that the combined action of a fungal root endophyte from a widespread taxon with core bacterial microbiota members provides synergistic protection against an aggressive soil-borne pathogen in Arabidopsis thaliana and barley. We additionally reveal early inter-kingdom growth promotion benefits which are host and microbiota composition dependent. Using RNA-sequencing, we show that these beneficial activities are not associated with extensive host transcriptional reprogramming but rather with the modulation of expression of microbial effectors and carbohydrate-active enzymes.

Evolutionary innovations through gain and loss of genes in the ectomycorrhizal Boletales

We aimed to identify genomic traits of transitions to ectomycorrhizal ecology within the Boletales by comparing the genomes of 21 symbiotrophic species with their saprotrophic brown-rot relatives. Gene duplication rate is constant along the backbone of Boletales phylogeny with large loss events in several lineages, while gene family expansion sharply increased in the late Miocene, mostly in the Boletaceae. Ectomycorrhizal Boletales have a reduced set of plant cell-wall-degrading enzymes (PCWDEs) compared with their brown-rot relatives. However, the various lineages retain distinct sets of PCWDEs, suggesting that, over their evolutionary history, symbiotic Boletales have become functionally diverse. A smaller PCWDE repertoire was found in Sclerodermatineae. The gene repertoire of several lignocellulose oxidoreductases (e.g. laccases) is similar in brown-rot and ectomycorrhizal species, suggesting that symbiotic Boletales are capable of mild lignocellulose decomposition. Transposable element (TE) proliferation contributed to the higher evolutionary rate of genes encoding effector-like small secreted proteins, proteases, and lipases. On the other hand, we showed that the loss of secreted CAZymes was not related to TE activity but to DNA decay. This study provides novel insights on our understanding of the mechanisms influencing the evolutionary diversification of symbiotic boletes.

Genetic determinants of endophytism in the Arabidopsis root mycobiome

The roots of Arabidopsis thaliana host diverse fungal communities that affect plant health and disease states. Here, we sequence the genomes of 41 fungal isolates representative of the A. thaliana root mycobiota for comparative analysis with other 79 plant-associated fungi. Our analyses indicate that root mycobiota members evolved from ancestors with diverse lifestyles and retain large repertoires of plant cell wall-degrading enzymes (PCWDEs) and effector-like small secreted proteins. We identify a set of 84 gene families associated with endophytism, including genes encoding PCWDEs acting on xylan (family GH10) and cellulose (family AA9). Transcripts encoding these enzymes are also part of a conserved transcriptional program activated by phylogenetically-distant mycobiota members upon host contact. Recolonization experiments with individual fungi indicate that strains with detrimental effects in mono-association with the host colonize roots more aggressively than those with beneficial activities, and dominate in natural root samples. Furthermore, we show that the pectin-degrading enzyme family PL1_7 links aggressiveness of endophytic colonization to plant health.

Evolution of the Mode of Nutrition in Symbiotic and Saprotrophic Fungi in Forest Ecosystems

In this review, we highlight the main insights that have been gathered from recent developments using large-scale genomics of fungal saprotrophs and symbiotrophs (including ectomycorrhizal and orchid and ericoid mycorrhizal fungi) inhabiting forest ecosystems. After assessing the goals and motivations underlying our approach, we explore our current understanding of the limits and future potential of using genomics to understand the ecological roles of these forest fungi. Comparative genomics unraveled the molecular machineries involved in lignocellulose decomposition in wood decayers, soil and litter saprotrophs, and mycorrhizal symbionts. They also showed that transitions from saprotrophy to mutualism entailed widespread losses of lignocellulose-degrading enzymes; diversification of novel, lineage-specific symbiosis-induced genes; and convergent evolution of genetic innovations that facilitate the accommodationof mutualistic symbionts within their plant hosts. We also identify the major questions that remain unanswered and propose new avenues of genome-based research to understand the role of soil fungi in sustainable forest ecosystems.

Genomic Analysis Enlightens Agaricales Lifestyle Evolution and Increasing Peroxidase Diversity

As actors of global carbon cycle, Agaricomycetes (Basidiomycota) have developed complex enzymatic machineries that allow them to decompose all plant polymers, including lignin. Among them, saprotrophic Agaricales are characterized by an unparalleled diversity of habitats and lifestyles. Comparative analysis of 52 Agaricomycetes genomes (14 of them sequenced de novo) reveals that Agaricales possess a large diversity of hydrolytic and oxidative enzymes for lignocellulose decay. Based on the gene families with the predicted highest evolutionary rates—namely cellulose-binding CBM1, glycoside hydrolase GH43, lytic polysaccharide monooxygenase AA9, class-II peroxidases, glucose–methanol–choline oxidase/dehydrogenases, laccases, and unspecific peroxygenases—we reconstructed the lifestyles of the ancestors that led to the extant lignocellulose-decomposing Agaricomycetes. The changes in the enzymatic toolkit of ancestral Agaricales are correlated with the evolution of their ability to grow not only on wood but also on leaf litter and decayed wood, with grass-litter decomposers as the most recent eco-physiological group. In this context, the above families were analyzed in detail in connection with lifestyle diversity. Peroxidases appear as a central component of the enzymatic toolkit of saprotrophic Agaricomycetes, consistent with their essential role in lignin degradation and high evolutionary rates. This includes not only expansions/losses in peroxidase genes common to other basidiomycetes but also the widespread presence in Agaricales (and Russulales) of new peroxidases types not found in wood-rotting Polyporales, and other Agaricomycetes orders. Therefore, we analyzed the peroxidase evolution in Agaricomycetes by ancestral-sequence reconstruction revealing several major evolutionary pathways and mapped the appearance of the different enzyme types in a time-calibrated species tree.

Conserved white-rot enzymatic mechanism for wood decay in the Basidiomycota genus Pycnoporus

White-rot (WR) fungi are pivotal decomposers of dead organic matter in forest ecosystems and typically use a large array of hydrolytic and oxidative enzymes to deconstruct lignocellulose. However, the extent of lignin and cellulose degradation may vary between species and wood type. Here, we combined comparative genomics, transcriptomics and secretome proteomics to identify conserved enzymatic signatures at the onset of wood-decaying activity within the Basidiomycota genus Pycnoporus. We observed a strong conservation in the genome structures and the repertoires of protein-coding genes across the four Pycnoporus species described to date, despite the species having distinct geographic distributions. We further analysed the early response of P. cinnabarinus, P. coccineus and P. sanguineus to diverse (ligno)-cellulosic substrates. We identified a conserved set of enzymes mobilized by the three species for breaking down cellulose, hemicellulose and pectin. The co-occurrence in the exo-proteomes of H₂O₂-producing enzymes with H₂O₂-consuming enzymes was a common feature of the three species, although each enzymatic partner displayed independent transcriptional regulation. Finally, cellobiose dehydrogenase-coding genes were systematically co-regulated with at least one AA9 lytic polysaccharide monooxygenase gene, indicative of enzymatic synergy in vivo. This study highlights a conserved core white-rot fungal enzymatic mechanism behind the wood-decaying process.

Desert truffle genomes reveal their reproductive modes and new insights into plant-fungal interaction and ectendomycorrhizal lifestyle

Desert truffles are edible hypogeous fungi forming ectendomycorrhizal symbiosis with plants of Cistaceae family. Knowledge about the reproductive modes of these fungi and the molecular mechanisms driving the ectendomycorrhizal interaction is lacking. Genomes of the highly appreciated edible desert truffles Terfezia claveryi Chatin and Tirmania nivea Trappe have been sequenced and compared with other Pezizomycetes. Transcriptomes of T. claveryi × Helianthemum almeriense mycorrhiza from well-watered and drought-stressed plants, when intracellular colonizations is promoted, were investigated. We have identified the fungal genes related to sexual reproduction in desert truffles and desert-truffles-specific genomic and secretomic features with respect to other Pezizomycetes, such as the expansion of a large set of gene families with unknown Pfam domains and a number of species or desert-truffle-specific small secreted proteins differentially regulated in symbiosis. A core set of plant genes, including carbohydrate, lipid-metabolism, and defence-related genes, differentially expressed in mycorrhiza under both conditions was found. Our results highlight the singularities of desert truffles with respect to other mycorrhizal fungi while providing a first glimpse on plant and fungal determinants involved in ecto to endo symbiotic switch that occurs in desert truffle under dry conditions.

Gene family expansions and transcriptome signatures uncover fungal adaptations to wood decay

Because they comprise some of the most efficient wood‐decayers, Polyporales fungi impact carbon cycling in forest environment. Despite continuous discoveries on the enzymatic machinery involved in wood decomposition, the vision on their evolutionary adaptation to wood decay and genome diversity remains incomplete. We combined the genome sequence information from 50 Polyporales species, including 26 newly sequenced genomes and sought for genomic and functional adaptations to wood decay through the analysis of genome composition and transcriptome responses to different carbon sources. The genomes of Polyporales from different phylogenetic clades showed poor conservation in macrosynteny, indicative of genome rearrangements. We observed different gene family expansion/contraction histories for plant cell wall degrading enzymes in core polyporoids and phlebioids and captured expansions for genes involved in signalling and regulation in the lineages of white rotters. Furthermore, we identified conserved cupredoxins, thaumatin‐like proteins and lytic polysaccharide monooxygenases with a yet uncharacterized appended module as new candidate players in wood decomposition. Given the current need for enzymatic toolkits dedicated to the transformation of renewable carbon sources, the observed genomic diversity among Polyporales strengthens the relevance of mining Polyporales biodiversity to understand the molecular mechanisms of wood decay.

AB0268 HUMAN T-CELL LEUKAEMIA VIRUS TYPE 1 MAY INVALIDATE T-SPOT.TB RESULTS AMONG RHEUMATOID ARTHRITIS PATIENTS: A RETROSPECTIVE OBSERVATIONAL STUDY

Background:

In clinical rheumatology, interferon-γ release assays (IGRAs) have been reported as a useful diagnostic test for latent tuberculosis infection (LTBI) before beginning the administration of biologics such as anti-TNF therapies (1). CD4-positive T cells are the main target in Human T-cell leukaemia virus type 1 (HTLV-1) infection. Several reports suggest that the reaction of tuberculin skin test (TST) is attenuated in HTLV-1-positive individuals compared with that in HTLV-1-negative individuals (2). However, it remains unclear whether IGRAs are reliable for detecting TB infection among HTLV-1-positive RA patients.

Objectives:

The present study aimed to investigate the usefulness of the T-SPOT.TBassay in HTLV-1-positive RA patients. In addition, the association between the existence of IFN-γ producing T cells and HTLV-1 proviral loads (PVLs) in HTLV-1-positive RA patients was analysed on the basis of the T-SPOT.TBassay results.

Methods:

We reviewed the medical records of 75 HTLV-1-negative and 29 HTLV-1-positive RA patients were suspected cases of LTBI and evaluated using the T-SPOT.TBassay as a clinical practice from April 2012 to July 2019. The results of T-SPOT.TBwere collected from medical records, retrospectively. Peripheral blood samples were obtained from HTLV-1-positive RA patients for the analysis of HTLV-1 PVLs values. The study protocol was approved by the research ethics committees of our hospitals.

Results:

Approximately 55% of the HTLV-1-positive RA patients showed invalid results for the T-SPOT.TBassay (p < 0.0001); the cause of invalid results was a spot-forming count of >10 spots in the negative controls of the T-SPOT.TBassay among HTLV-1-positive RA patients. Among HTLV-1-positive RA patients, HTLV-1 PVL values were significantly higher in 16 patients who showed invalid results than in 13 patients who did not (p = 0.003). There were no between-group differences in female patient ratio, age, RA disease activity and therapeutic regimens. IFN-γ producing cells were detected in the peripheral blood of HTLV-1-positive RA patients without stimulation with TB-specific antigens.

Conclusion:

The incidence of invalid results for the T-SPOT.TBassay has been reported to be as low as 0.6% (3). The results of this assay for screening of LTBI in HTLV-1-positive RA patients should be interpreted with caution. Furthermore, our results show that an increase in IFN-γ producing T cell numbers due to HTLV-1 infection in RA patients may affect the pathogenesis of RA.

References:

[1]Iannone, F., et al.J. Rheumatol. Suppl.91, 41-46 (2014).

[2]Tachibana, N., et al.Int. J. Cancer42, 829-831 (1988).

[3]Rego, K., et al.Tuberculosis (Edinb.)108, 178-185 (2018).

Acknowledgments:

We would like to thank Dr Yuki Hashikura and Ms Yuki Kaseda of the University of Miyazaki for their technical support in this work. We would also like to acknowledge Ms Yumiko Kai at the Institute of Rheumatology, Zenjinkai Shimin-no-Mori Hospital, for her help in data management.

A part this work was supported by a grant from the Practical Research Project for Rare/Intractable Diseases of the Japan Agency for Medical Research and Development (Grant No. JP19ek0109356), a Health and Labor Sciences Research Grant on Rare and Intractable Diseases from the Ministry of Health, Labor and Welfare of Japan (Grant No. 19FC1007), and a Grant-in-Aid for Clinical Research from Miyazaki University Hospital.

Disclosure of Interests:

Kunihiko Umekita Paid instructor for: Astellas Pharma Inc. Chugai Pharma Inc. Tanabe-Mitsubishi Pharma Inc., Speakers bureau: Bristol-Myers Squibb, Yayoi Hashiba: None declared, Risa Kudou: None declared, Shunichi Miyauchi: None declared, Masatoshi Kimura: None declared, Motohiro Matsuda: None declared, Chihiro Iwao: None declared, Yumi Kariya: None declared, Takeshi Kawaguchi: None declared, Katoko Takajo: None declared, Koushou Iwao: None declared, Yuuki Rikitake: None declared, Ichiro Takajo: None declared, Toshihiko Hidaka Paid instructor for: Astellas Pharma Inc. Chugai Pharma Inc. Tanabe-Mitsubishi Pharma Inc., Speakers bureau: Astellas Pharma Inc. Chugai Pharma Inc. Tanabe-Mitsubishi Pharma Inc., Akihiko Okayama: None declared

An ectomycorrhizal fungus alters sensitivity to jasmonate, salicylate, gibberellin, and ethylene in host roots

The phytohormones jasmonate, gibberellin, salicylate, and ethylene regulate an interconnected reprogramming network integrating root development with plant responses against microbes. The establishment of mutualistic ectomycorrhizal symbiosis requires the suppression of plant defense responses against fungi as well as the modification of root architecture and cortical cell wall properties. Here, we investigated the contribution of phytohormones and their crosstalk to the ontogenesis of ectomycorrhizae (ECM) between grey poplar (Populus tremula x alba) roots and the fungus Laccaria bicolor. To obtain the hormonal blueprint of developing ECM, we quantified the concentrations of jasmonates, gibberellins, and salicylate via liquid chromatography-tandem mass spectrometry. Subsequently, we assessed root architecture, mycorrhizal morphology, and gene expression levels (RNA sequencing) in phytohormone-treated poplar lateral roots in the presence or absence of L. bicolor. Salicylic acid accumulated in mid-stage ECM. Exogenous phytohormone treatment affected the fungal colonization rate and/or frequency of Hartig net formation. Colonized lateral roots displayed diminished responsiveness to jasmonate but regulated some genes, implicated in defense and cell wall remodelling, that were specifically differentially expressed after jasmonate treatment. Responses to salicylate, gibberellin, and ethylene were enhanced in ECM. The dynamics of phytohormone accumulation and response suggest that jasmonate, gibberellin, salicylate, and ethylene signalling play multifaceted roles in poplar L. bicolor ectomycorrhizal development.

P692Chemoradiation therapy to patients with lung cancer exacerbates thoracic aortic calcification

Introduction

Development of chemoradiation therapy (CRT) has improved mortality in patients with cancer. Whereas, it is emerging problem that cancer-survivors suffer from cardiovascular diseases, and the association between modern CRT and the increase in future cardiovascular events is suggested. Meanwhile, previous studies showed that thoracic aortic calcification (TAC) detected by computed tomography (CT), a marker of atherosclerosis, was associated with all-cause mortality and cardiovascular events. However, the influence of CRT on TAC progression remains unclear.

Purpose

The purpose of this study was to evaluate whether CRT would exacerbate TAC.

Methods

A total of 68 patients who treated lung cancer at our hospital between 2011 and 2015 were retrospectively analyzed (mean 62 year-old, male 78%): 35 patients underwent surgical treatment after induction CRT (CRT group) and 33 patients underwent surgical treatment alone (control group), extracted by propensity score matching by age, sex, smoking status, and diseased side. The volume of TAC between 2nd and 12th thoracic vertebrae was quantitatively measured with CT imaging, at baseline and at 1 year follow-up. The annual percent change in TAC was compared between the CRT and the control group. Moreover, the independent relationship between implementation of CRT and the progression of TAC was assessed by multivariate logistic regression analysis, adjusting for age, gender, conventional atherosclerotic risk factors and baseline aortic calcification volume.

Results

Patients in the CRT group received radiation (mean 47.3±4.0 Gy) and chemotherapy: 2 courses of cisplatin with docetaxel (34 cases) or vinorelbine (1 case). The prevalence of dyslipidemia, taking statins and diabetes drugs were significantly higher in the control groups (17% vs. 39%; p=0.041, 11% vs. 33%; p=0.029, 3% vs. 18%; p=0.044, respectively). Baseline C-reactive protein level was significantly higher in the CRT group (0.255 vs. 0.115; p=0.034). In univariate analysis, the annual percent change in TAC volume was significantly increased in the CRT group compared with the control group (37.6% vs. 23.3%; p=0.006). Multivariate logistic regression analysis demonstrated that CRT was an independent factor associated with the progression of TAC volume, even after adjustment for baseline calcification volume and coronary risk factors (OR, 3.90; 95% CI, 1.32–11.47; p=0.014).

Conclusion

CRT to patients with lung cancer exacerbates thoracic aortic calcification, which may result in future cardiovascular events.

Insights into an unusual Auxiliary Activity 9 family member lacking the histidine brace motif of lytic polysaccharide monooxygenases

Lytic polysaccharide monooxygenases (LPMOs) are redox-enzymes involved in biomass degradation. All characterized LPMOs possess an active site of two highly conserved histidine residues coordinating a copper ion (the histidine brace), which are essential for LPMO activity. However, some protein sequences that belong to the AA9 LPMO family display a natural N-terminal His to Arg substitution (Arg-AA9). These are found almost entirely in the phylogenetic fungal class Agaricomycetes, associated with wood decay, but no function has been demonstrated for any Arg-AA9. Through bioinformatics, transcriptomic, and proteomic analyses we present data, which suggest that Arg-AA9 proteins could have a hitherto unidentified role in fungal degradation of lignocellulosic biomass in conjunction with other secreted fungal enzymes. We present the first structure of an Arg-AA9, LsAA9B, a naturally occurring protein from Lentinus similis The LsAA9B structure reveals gross changes in the region equivalent to the canonical LPMO copper-binding site, whereas features implicated in carbohydrate binding in AA9 LPMOs have been maintained. We obtained a structure of LsAA9B with xylotetraose bound on the surface of the protein although with a considerably different binding mode compared with other AA9 complex structures. In addition, we have found indications of protein phosphorylation near the N-terminal Arg and the carbohydrate-binding site, for which the potential function is currently unknown. Our results are strong evidence that Arg-AA9s function markedly different from canonical AA9 LPMO, but nonetheless, may play a role in fungal conversion of lignocellulosic biomass.

Comparative genomics of Rhizophagus irregularis, R. cerebriforme, R. diaphanus and Gigaspora rosea highlights specific genetic features in Glomeromycotina

Glomeromycotina is a lineage of early diverging fungi that establish arbuscular mycorrhizal (AM) symbiosis with land plants. Despite their major ecological role, the genetic basis of their obligate mutualism remains largely unknown, hindering our understanding of their evolution and biology. We compared the genomes of Glomerales (Rhizophagus irregularis, Rhizophagus diaphanus, Rhizophagus cerebriforme) and Diversisporales (Gigaspora rosea) species, together with those of saprotrophic Mucoromycota, to identify gene families and processes associated with these lineages and to understand the molecular underpinning of their symbiotic lifestyle. Genomic features in Glomeromycotina appear to be very similar with a very high content in transposons and protein-coding genes, extensive duplications of protein kinase genes, and loss of genes coding for lignocellulose degradation, thiamin biosynthesis and cytosolic fatty acid synthase. Most symbiosis-related genes in R. irregularis and G. rosea are specific to Glomeromycotina. We also confirmed that the present species have a homokaryotic genome organisation. The high interspecific diversity of Glomeromycotina gene repertoires, affecting all known protein domains, as well as symbiosis-related orphan genes, may explain the known adaptation of Glomeromycotina to a wide range of environmental settings. Our findings contribute to an increasingly detailed portrait of genomic features defining the biology of AM fungi.

Dynamics of the Phanerochaete carnosa transcriptome during growth on aspen and spruce

BACKGROUND

The basidiomycete Phanerochaete carnosa is a white-rot species that has been mainly isolated from coniferous softwood. Given the particular recalcitrance of softwoods to bioconversion, we conducted a comparative transcriptomic analysis of P. carnosa following growth on wood powder from one softwood (spruce; Picea glauca) and one hardwood (aspen; Populus tremuloides). P. carnosa was grown on each substrate for over one month, and mycelia were harvested at five time points for total RNA sequencing. Residual wood powder was also analyzed for total sugar and lignin composition.

RESULTS

Following a slightly longer lag phase of growth on spruce, radial expansion of the P. carnosa colony was similar on spruce and aspen. Consistent with this observation, the pattern of gene expression by P. carnosa on each substrate converged following the initial adaptation. On both substrates, highest transcript abundances were attributed to genes predicted to encode manganese peroxidases (MnP), along with auxiliary activities from carbohydrate-active enzyme (CAZy) families AA3 and AA5. In addition, a lytic polysaccharide monooxygenase from family AA9 was steadily expressed throughout growth on both substrates. P450 sequences from clans CPY52 and CYP64 accounted for 50% or more of the most highly expressed P450s, which were also the P450 clans that were expanded in the P. carnosa genome relative to other white-rot fungi.

CONCLUSIONS

The inclusion of five growth points and two wood substrates was important to revealing differences in the expression profiles of specific sequences within large glycoside hydrolase families (e.g., GH5 and GH16), and permitted co-expression analyses that identified new targets for study, including non-catalytic proteins and proteins with unknown function.

Integrative visual omics of the white-rot fungus Polyporus brumalis exposes the biotechnological potential of its oxidative enzymes for delignifying raw plant biomass

BACKGROUND

Plant biomass conversion for green chemistry and bio-energy is a current challenge for a modern sustainable bioeconomy. The complex polyaromatic lignin polymers in raw biomass feedstocks (i.e., agriculture and forestry by-products) are major obstacles for biomass conversions. White-rot fungi are wood decayers able to degrade all polymers from lignocellulosic biomass including cellulose, hemicelluloses, and lignin. The white-rot fungus Polyporus brumalis efficiently breaks down lignin and is regarded as having a high potential for the initial treatment of plant biomass in its conversion to bio-energy. Here, we describe the extraordinary ability of P. brumalis for lignin degradation using its enzymatic arsenal to break down wheat straw, a lignocellulosic substrate that is considered as a biomass feedstock worldwide.

RESULTS

We performed integrative multi-omics analyses by combining data from the fungal genome, transcriptomes, and secretomes. We found that the fungus possessed an unexpectedly large set of genes coding for Class II peroxidases involved in lignin degradation (19 genes) and GMC oxidoreductases/dehydrogenases involved in generating the hydrogen peroxide required for lignin peroxidase activity and promoting redox cycling of the fungal enzymes involved in oxidative cleavage of lignocellulose polymers (36 genes). The examination of interrelated multi-omics patterns revealed that eleven Class II Peroxidases were secreted by the fungus during fermentation and eight of them where tightly co-regulated with redox cycling enzymatic partners.

CONCLUSION

As a peculiar feature of P. brumalis, we observed gene family extension, up-regulation and secretion of an abundant set of versatile peroxidases and manganese peroxidases, compared with other Polyporales species. The orchestrated secretion of an abundant set of these delignifying enzymes and redox cycling enzymatic partners could contribute to the delignification capabilities of the fungus. Our findings highlight the diversity of wood decay mechanisms present in Polyporales and the potentiality of further exploring this taxonomic order for enzymatic functions of biotechnological interest.

The integrative omics of white-rot fungus Pycnoporus coccineus reveals co-regulated CAZymes for orchestrated lignocellulose breakdown

Innovative green technologies are of importance for converting plant wastes into renewable sources for materials, chemicals and energy. However, recycling agricultural and forestry wastes is a challenge. A solution may be found in the forest. Saprotrophic white-rot fungi are able to convert dead plants into consumable carbon sources. Specialized fungal enzymes can be utilized for breaking down hard plant biopolymers. Thus, understanding the enzymatic machineries of such fungi gives us hints for the efficient decomposition of plant materials. Using the saprotrophic white-rot fungus Pycnoporus coccineus as a fungal model, we examined the dynamics of transcriptomic and secretomic responses to different types of lignocellulosic substrates at two time points. Our integrative omics pipeline (SHIN+GO) enabled us to compress layers of biological information into simple heatmaps, allowing for visual inspection of the data. We identified co-regulated genes with corresponding co-secreted enzymes, and the biological roles were extrapolated with the enriched Carbohydrate-Active Enzyme (CAZymes) and functional annotations. We observed the fungal early responses for the degradation of lignocellulosic substrates including; 1) simultaneous expression of CAZy genes and secretion of the enzymes acting on diverse glycosidic bonds in cellulose, hemicelluloses and their side chains or lignin (i.e. hydrolases, esterases and oxido-reductases); 2) the key role of lytic polysaccharide monooxygenases (LPMO); 3) the early transcriptional regulation of lignin active peroxidases; 4) the induction of detoxification processes dealing with biomass-derived compounds; and 5) the frequent attachments of the carbohydrate binding module 1 (CBM1) to enzymes from the lignocellulose-responsive genes. Our omics combining methods and related biological findings may contribute to the knowledge of fungal systems biology and facilitate the optimization of fungal enzyme cocktails for various industrial applications.

Sphingomyelin Phosphodiesterase 3 Enhances Cytodifferentiation of Periodontal Ligament Cells

Sphingomyelin phosphodiesterase 3 ( Smpd3), which encodes neutral sphingomyelinase 2 (nSMase2), is a key molecule for skeletal development as well as for the cytodifferentiation of odontoblasts and alveolar bone. However, the effects of nSMase2 on the cytodifferentiation of periodontal ligament (PDL) cells are still unclear. In this study, the authors analyzed the effects of Smpd3 on the cytodifferentiation of human PDL (HPDL) cells. The authors found that Smpd3 increases the mRNA expression of calcification-related genes, such as alkaline phosphatase (ALPase), type I collagen, osteopontin, Osterix (Osx), and runt-related transcription factor (Runx)-2 in HPDL cells. In contrast, GW4869, an inhibitor of nSMase2, clearly decreased the mRNA expression of ALPase, type I collagen, and osteocalcin in HPDL cells, suggesting that Smpd3 enhances HPDL cytodifferentiation. Next, the authors used exome sequencing to evaluate the genetic variants of Smpd3 in a Japanese population with aggressive periodontitis (AgP). Among 44 unrelated subjects, the authors identified a single nucleotide polymorphism (SNP), rs145616324, in Smpd3 as a putative genetic variant for AgP among Japanese people. Moreover, Smpd3 harboring this SNP did not increase the sphingomyelinase activity or mRNA expression of ALPase, type I collagen, osteopontin, Osx, or Runx2, suggesting that this SNP inhibits Smpd3 such that it has no effect on the cytodifferentiation of HPDL cells. These data suggest that Smpd3 plays a crucial role in maintaining the homeostasis of PDL tissue.

Visual Comparative Omics of Fungi for Plant Biomass Deconstruction

Wood-decay fungi contain the cellular mechanisms to decompose such plant cell wall components as cellulose, hemicellulose, and lignin. A multi-omics approach to the comparative analysis of wood-decay fungi gives not only new insights into their strategies for decomposing recalcitrant plant biomass, but also an understanding of how to exploit these mechanisms for biotechnological applications. We have developed an analytical workflow, Applied Biomass Conversion Design for Efficient Fungal Green Technology (ABCDEFGT), to simplify the analysis and interpretation of transcriptomic and secretomic data. ABCDEFGT utilizes self-organizing maps for grouping genes with similar transcription patterns, and an overlay with secreted proteins. The key feature of ABCDEFGT is simple graphic outputs of genome-wide transcriptomic and secretomic topographies, which enables visual inspection without a priori of the omics data and facilitates discoveries of co-regulated genes and proteins. Genome-wide omics landscapes were built with the newly sequenced fungal species Pycnoporus coccineus, Pycnoporus sanguineus, and Pycnoporus cinnabarinus grown on various carbon sources. Integration of the post-genomic data revealed a global overlap, confirming the pertinence of the genome-wide approach. ABCDEFGT was evaluated by comparison with the latest clustering method for ease of output interpretation, and ABCDEFGT gave a better biological representation of fungal behaviors. The genome-wide multi-omics strategy allowed us to determine the potential synergy of particular enzymes decomposing cellulose, hemicellulose, and lignin such as Lytic Polysaccharide Monooxygenases, modular enzymes associated with a cellulose binding module1, and Class II Peroxidase isoforms co-regulated with oxido-reductases. Overall, ABCDEFGT was capable of visualizing genome-wide transcriptional and secretomic profiles for intuitive interpretations and is suitable for exploration of newly-sequenced organisms.

Human Papilloma Virus Identification in Breast Cancer Patients with Previous Cervical Neoplasia

Purpose

Women with human papilloma virus (HPV)-associated cervical neoplasia have a higher risk of developing breast cancer than the general female population. The purpose of this study was to (i) identify high-risk HPVs in cervical neoplasia and subsequent HPV positive breast cancers which developed in the same patients and (ii) determine if these HPVs were biologically active.

Methods

A range of polymerase chain reaction and immunohistochemical techniques were used to conduct a retrospective cohort study of cervical precancers and subsequent breast cancers in the same patients.

Results

The same high-risk HPV types were identified in both the cervical and breast specimens in 13 (46%) of 28 patients. HPV type 18 was the most prevalent. HPVs appeared to be biologically active as demonstrated by the expression of HPV E7 proteins and the presence of HPV-associated koilocytes. The average age of these patients diagnosed with breast cancer following prior cervical precancer was 51 years, as compared to 60 years for all women with breast cancer (p for difference = 0.001).

Conclusion

These findings indicate that high-risk HPVs can be associated with cervical neoplasia and subsequent young age breast cancer. However, these associations are unusual and are a very small proportion of breast cancers. These outcomes confirm and extend the observations of two similar previous studies and offer one explanation for the increased prevalence of serious invasive breast cancer among young women.

Human Papilloma Viruses and Breast Cancer

PURPOSE

Human papillomaviruses (HPV) may have a role in some breast cancers. The purpose of this study is to fill important gaps in the evidence. These gaps are: (i) confirmation of the presence of high risk for cancer HPVs in breast cancers, (ii) evidence of HPV infections in benign breast tissues prior to the development of HPV-positive breast cancer in the same patients, (iii) evidence that HPVs are biologically active and not harmless passengers in breast cancer.

METHODS

RNA-seq data from The Cancer Genome Atlas (TCGA) was used to identify HPV RNA sequences in breast cancers. We also conducted a retrospective cohort study based on polymerase chain reaction (PCR) analyses to identify HPVs in archival specimens from Australian women with benign breast biopsies who later developed breast cancer. To assess whether HPVs in breast cancer were biologically active, the expression of the oncogenic protein HPV E7 was assessed by immunohistochemistry (IHC).

RESULTS

Thirty (3.5%) low-risk and 20 (2.3%) high-risk HPV types were identified in 855 breast cancers from the TCGA database. The high risk types were HPV 18 (48%), HPV 113 (24%), HPV 16 (10%), HPV 52 (10%). Data from the PCR cohort study indicated that HPV type 18 was the most common type identified in breast cancer specimens (55% of 40 breast cancer specimens) followed by HPV 16 (13%). The same HPV type was identified in both the benign and subsequent breast cancer in 15 patients. HPV E7 proteins were identified in 72% of benign breast specimens and 59% of invasive breast cancer specimens.

CONCLUSION

There were four observations of particular interest: (i) confirmation by both NGS and PCR of the presence of high-risk HPV gene sequences in breast cancers, (ii) a correlation between high-risk HPV in benign breast specimens and subsequent HPV-positive breast cancer in the same patient, (iii) HPVs in breast cancer are likely to be biologically active (as shown by transcription of HPV DNA to RNA plus the expression of HPV E7 proteins), (iv) HPV oncogenic influences may occur early in the development of breast cancer.