A single cell atlas of frozen shoulder capsule identifies features associated with inflammatory fibrosis resolution.

Frozen shoulder is a spontaneously self-resolving chronic inflammatory fibrotic human disease, which distinguishes the condition from most fibrotic diseases that are progressive and irreversible. Using single-cell analysis, we identify pro-inflammatory MERTKlowCD48+ macrophages and MERTK + LYVE1 + MRC1+ macrophages enriched for negative regulators of inflammation which co-exist in frozen shoulder capsule tissues. Micro-cultures of patient-derived cells identify integrin-mediated cell-matrix interactions between MERTK+ macrophages and pro-resolving DKK3+ and POSTN+ fibroblasts, suggesting that matrix remodelling plays a role in frozen shoulder resolution. Cross-tissue analysis reveals a shared gene expression cassette between shoulder capsule MERTK+ macrophages and a respective population enriched in synovial tissues of rheumatoid arthritis patients in disease remission, supporting the concept that MERTK+ macrophages mediate resolution of inflammation and fibrosis. Single-cell transcriptomic profiling and spatial analysis of human foetal shoulder tissues identify MERTK + LYVE1 + MRC1+ macrophages and DKK3+ and POSTN+ fibroblast populations analogous to those in frozen shoulder, suggesting that the template to resolve fibrosis is established during shoulder development. Crosstalk between MerTK+ macrophages and pro-resolving DKK3+ and POSTN+ fibroblasts could facilitate resolution of frozen shoulder, providing a basis for potential therapeutic resolution of persistent fibrotic diseases.

Single-cell insights into immune dysregulation in rheumatoid arthritis flare versus drug-free remission.

Immune-mediated inflammatory diseases (IMIDs) are typically characterised by relapsing and remitting flares of inflammation. However, the unpredictability of disease flares impedes their study. Addressing this critical knowledge gap, we use the experimental medicine approach of immunomodulatory drug withdrawal in rheumatoid arthritis (RA) remission to synchronise flare processes allowing detailed characterisation. Exploratory mass cytometry analyses reveal three circulating cellular subsets heralding the onset of arthritis flare - CD45RO+PD1hi CD4+ and CD8+ T cells, and CD27+CD86+CD21- B cells - further characterised by single-cell sequencing. Distinct lymphocyte subsets including cytotoxic and exhausted CD4+ memory T cells, memory CD8+CXCR5+ T cells, and IGHA1+ plasma cells are primed for activation in flare patients. Regulatory memory CD4+ T cells (Treg cells) increase at flare onset, but with dysfunctional regulatory marker expression compared to drug-free remission. Significant clonal expansion is observed in T cells, but not B cells, after drug cessation; this is widespread throughout memory CD8+ T cell subsets but limited to the granzyme-expressing cytotoxic subset within CD4+ memory T cells. Based on our observations, we suggest a model of immune dysregulation for understanding RA flare, with potential for further translational research towards novel avenues for its treatment and prevention.

Pro Nerve Growth Factor and Its Receptor p75NTR Activate Inflammatory Responses in Synovial Fibroblasts: A Novel Targetable Mechanism in Arthritis.

We have recently provided new evidence for a role of p75NTR receptor and its preferential ligand proNGF in amplifying inflammatory responses in synovial mononuclear cells of chronic arthritis patients. In the present study, to better investigate how activation of the p75NTR/proNGF axis impacts synovial inflammation, we have studied the effects of proNGF on fibroblast-like synoviocytes (FLS), which play a central role in modulating local immune responses and in activating pro-inflammatory pathways. Using single cell RNA sequencing in synovial tissues from active and treatment-naïve rheumatoid arthritis (RA) patients, we demonstrated that p75NTR and sortilin, which form a high affinity receptor complex for proNGF, are highly expressed in PRG4pos lining and THY1posCOL1A1pos sublining fibroblast clusters in RA synovia but decreased in RA patients in sustained clinical remission. In ex vivo experiments we found that FLS from rheumatoid arthritis patients (RA-FLS) retained in vitro a markedly higher expression of p75NTR and sortilin than FLS from osteoarthritis patients (OA-FLS). Inflammatory stimuli further up-regulated p75NTR expression and induced endogenous production of proNGF in RA-FLS, leading to an autocrine activation of the proNGF/p75NTR pathway that results in an increased release of pro-inflammatory cytokines. Our data on the inhibition of p75NTR receptor, which reduced the release of IL-1ß, IL-6 and TNF-α, further confirmed the key role of p75NTR activation in regulating inflammatory cytokine production. In a set of ex vivo experiments, we used RA-FLS and cultured them in the presence of synovial fluids obtained from arthritis patients that, as we demonstrated, are characterized by a high concentration of proNGF. Our data show that the high levels of proNGF present in inflamed synovial fluids induced pro-inflammatory cytokine production by RA-FLS. The blocking of NGF binding to p75NTR using specific inhibitors led instead to the disruption of this pro-inflammatory loop, reducing activation of the p38 and JNK intracellular pathways and decreasing inflammatory cytokine production. Overall, our data demonstrate that an active proNGF/p75NTR axis promotes pro-inflammatory responses in synovial fibroblasts, thereby contributing to chronic synovial inflammation, and point to the possible use of p75NTR inhibitors as a novel therapeutic approach in chronic arthritis.

Role of synovial fibroblast subsets across synovial pathotypes in rheumatoid arthritis: a deconvolution analysis.

OBJECTIVES: To integrate published single-cell RNA sequencing (scRNA-seq) data and assess the contribution of synovial fibroblast (SF) subsets to synovial pathotypes and respective clinical characteristics in treatment-naïve early arthritis. METHODS: In this in silico study, we integrated scRNA-seq data from published studies with additional unpublished in-house data. Standard Seurat, Harmony and Liger workflow was performed for integration and differential gene expression analysis. We estimated single cell type proportions in bulk RNA-seq data (deconvolution) from synovial tissue from 87 treatment-naïve early arthritis patients in the Pathobiology of Early Arthritis Cohort using MuSiC. SF proportions across synovial pathotypes (fibroid, lymphoid and myeloid) and relationship of disease activity measurements across different synovial pathotypes were assessed. RESULTS: We identified four SF clusters with respective marker genes: PRG4+ SF (CD55, MMP3, PRG4, THY1neg ); CXCL12+ SF (CXCL12, CCL2, ADAMTS1, THY1low ); POSTN+ SF (POSTN, collagen genes, THY1); CXCL14+ SF (CXCL14, C3, CD34, ASPN, THY1) that correspond to lining (PRG4+ SF) and sublining (CXCL12+ SF, POSTN+ + and CXCL14+ SF) SF subsets. CXCL12+ SF and POSTN+ + were most prominent in the fibroid while PRG4+ SF appeared highest in the myeloid pathotype. Corresponding, lining assessed by histology (assessed by Krenn-Score) was thicker in the myeloid, but also in the lymphoid pathotype + the fibroid pathotype. PRG4+ SF correlated positively with disease severity parameters in the fibroid, POSTN+ SF in the lymphoid pathotype whereas CXCL14+ SF showed negative association with disease severity in all pathotypes. CONCLUSION: This study shows a so far unexplored association between distinct synovial pathologies and SF subtypes defined by scRNA-seq. The knowledge of the diverse interplay of SF with immune cells will advance opportunities for tailored targeted treatments.

miR-155-overexpressing monocytes resemble HLAhighISG15+ synovial tissue macrophages from patients with rheumatoid arthritis and induce polyfunctional CD4+ T-cell activation.

MicroRNAs (miRs) are known to regulate pro-inflammatory effector functions of myeloid cells, and miR dysregulation is implicated in rheumatoid arthritis (RA), a condition characterized by inflammation and destruction of the joints. We showed previously that miR-155 is increased in myeloid cells in RA and induces pro-inflammatory activation of monocytes and macrophages; however, its role at the interface between innate and adaptive immunity was not defined. Here, RNA-sequencing revealed that overexpression of miR-155 in healthy donor monocytes conferred a specific gene profile which bears similarities to that of RA synovial fluid-derived CD14+ cells and HLAhighISG15+ synovial tissue macrophages, both of which are characterized by antigen-presenting pathways. In line with this, monocytes in which miR-155 was overexpressed, displayed increased expression of HLA-DR and both co-stimulatory and co-inhibitory molecules, and induced activation of polyfunctional T cells. Together, these data underpin the notion that miR-155-driven myeloid cell activation in the synovium contributes not only to inflammation but may also influence the adaptive immune response.

Inflammation causes remodeling of mitochondrial cytochrome c oxidase mediated by the bifunctional gene C15orf48.

Dysregulated mitochondrial function is a hallmark of immune-mediated inflammatory diseases. Cytochrome c oxidase (CcO), which mediates the rate-limiting step in mitochondrial respiration, is remodeled during development and in response to changes of oxygen availability, but there has been little study of CcO remodeling during inflammation. Here, we describe an elegant molecular switch mediated by the bifunctional transcript C15orf48, which orchestrates the substitution of the CcO subunit NDUFA4 by its paralog C15ORF48 in primary macrophages. Expression of C15orf48 is a conserved response to inflammatory signals and occurs in many immune-related pathologies. In rheumatoid arthritis, C15orf48 mRNA is elevated in peripheral monocytes and proinflammatory synovial tissue macrophages, and its expression positively correlates with disease severity and declines in remission. C15orf48 is also expressed by pathogenic macrophages in severe coronavirus disease 2019 (COVID-19). Study of a rare metabolic disease syndrome provides evidence that loss of the NDUFA4 subunit supports proinflammatory macrophage functions.

Translational targeting of inflammation and fibrosis in frozen shoulder: Molecular dissection of the T cell/IL-17A axis.

Frozen shoulder is a common fibroproliferative disease characterized by the insidious onset of pain and restricted range of shoulder movement with a significant socioeconomic impact. The pathophysiological mechanisms responsible for chronic inflammation and matrix remodeling in this prevalent fibrotic disorder remain unclear; however, increasing evidence implicates dysregulated immunobiology. IL-17A is a key cytokine associated with inflammation and tissue remodeling in numerous musculoskeletal diseases, and thus, we sought to determine the role of IL-17A in the immunopathogenesis of frozen shoulder. We demonstrate an immune cell landscape that switches from a predominantly macrophage population in nondiseased tissue to a T cell-rich environment in disease. Furthermore, we observed a subpopulation of IL-17A-producing T cells capable of inducing profibrotic and inflammatory responses in diseased fibroblasts through enhanced expression of the signaling receptor IL-17RA, rendering diseased cells more sensitive to IL-17A. We further established that the effects of IL-17A on diseased fibroblasts was TRAF-6/NF-κB dependent and could be inhibited by treatment with an IKKß inhibitor or anti-IL-17A antibody. Accordingly, targeting of the IL-17A pathway may provide future therapeutic approaches to the management of this common, debilitating disease.

COVID-19 and RA share an SPP1 myeloid pathway that drives PD-L1+ neutrophils and CD14+ monocytes.

We explored the potential link between chronic inflammatory arthritis and COVID-19 pathogenic and resolving macrophage pathways and their role in COVID-19 pathogenesis. We found that bronchoalveolar lavage fluid (BALF) macrophage clusters FCN1+ and FCN1+SPP1+ predominant in severe COVID-19 were transcriptionally related to synovial tissue macrophage (STM) clusters CD48hiS100A12+ and CD48+SPP1+ that drive rheumatoid arthritis (RA) synovitis. BALF macrophage cluster FABP4+ predominant in healthy lung was transcriptionally related to STM cluster TREM2+ that governs resolution of synovitis in RA remission. Plasma concentrations of SPP1 and S100A12 (key products of macrophage clusters shared with active RA) were high in severe COVID-19 and predicted the need for Intensive Care Unit transfer, and they remained high in the post-COVID-19 stage. High plasma levels of SPP1 were unique to severe COVID-19 when compared with other causes of severe pneumonia, and IHC localized SPP1+ macrophages in the alveoli of COVID-19 lung. Investigation into SPP1 mechanisms of action revealed that it drives proinflammatory activation of CD14+ monocytes and development of PD-L1+ neutrophils, both hallmarks of severe COVID-19. In summary, COVID-19 pneumonitis appears driven by similar pathogenic myeloid cell pathways as those in RA, and their mediators such as SPP1 might be an upstream activator of the aberrant innate response in severe COVID-19 and predictive of disease trajectory including post-COVID-19 pathology.

Mitochondria as Key Players in the Pathogenesis and Treatment of Rheumatoid Arthritis.

Mitochondria are major energy-producing organelles that have central roles in cellular metabolism. They also act as important signalling hubs, and their dynamic regulation in response to stress signals helps to dictate the stress response of the cell. Rheumatoid arthritis is an inflammatory and autoimmune disease with high prevalence and complex aetiology. Mitochondrial activity affects differentiation, activation and survival of immune and non-immune cells that contribute to the pathogenesis of this disease. This review outlines what is known about the role of mitochondria in rheumatoid arthritis pathogenesis, and how current and future therapeutic strategies can function through modulation of mitochondrial activity. We also highlight areas of this topic that warrant further study. As producers of energy and of metabolites such as succinate and citrate, mitochondria help to shape the inflammatory phenotype of leukocytes during disease. Mitochondrial components can directly stimulate immune receptors by acting as damage-associated molecular patterns, which could represent an initiating factor for the development of sterile inflammation. Mitochondria are also an important source of intracellular reactive oxygen species, and facilitate the activation of the NLRP3 inflammasome, which produces cytokines linked to disease symptoms in rheumatoid arthritis. The fact that mitochondria contain their own genetic material renders them susceptible to mutation, which can propagate their dysfunction and immunostimulatory potential. Several drugs currently used for the treatment of rheumatoid arthritis regulate mitochondrial function either directly or indirectly. These actions contribute to their immunomodulatory functions, but can also lead to adverse effects. Metabolic and mitochondrial pathways are attractive targets for future anti-rheumatic drugs, however many questions still remain about the precise role of mitochondrial activity in different cell types in rheumatoid arthritis.

Novel self-amplificatory loop between T cells and tenocytes as a driver of chronicity in tendon disease.

OBJECTIVES: Increasing evidence suggests that inflammatory mechanisms play a key role in chronic tendon disease. After observing T cell signatures in human tendinopathy, we explored the interaction between T cells and tendon stromal cells or tenocytes to define their functional contribution to tissue remodelling and inflammation amplification and hence disease perpetuation. METHODS: T cells were quantified and characterised in healthy and tendinopathic tissues by flow cytometry (FACS), imaging mass cytometry (IMC) and single cell RNA-seq. Tenocyte activation induced by conditioned media from primary damaged tendon or interleukin-1ß was evaluated by qPCR. The role of tenocytes in regulating T cell migration was interrogated in a standard transwell membrane system. T cell activation (cell surface markers by FACS and cytokine release by ELISA) and changes in gene expression in tenocytes (qPCR) were assessed in cocultures of T cells and explanted tenocytes. RESULTS: Significant quantitative differences were observed in healthy compared with tendinopathic tissues. IMC showed T cells in close proximity to tenocytes, suggesting tenocyte-T cell interactions. On activation, tenocytes upregulated inflammatory cytokines, chemokines and adhesion molecules implicated in T cell recruitment and activation. Conditioned media from activated tenocytes induced T cell migration and coculture of tenocytes with T cells resulted in reciprocal activation of T cells. In turn, these activated T cells upregulated production of inflammatory mediators in tenocytes, while increasing the pathogenic collagen 3/collagen 1 ratio. CONCLUSIONS: Interaction between T cells and tenocytes induces the expression of inflammatory cytokines/chemokines in tenocytes, alters collagen composition favouring collagen 3 and self-amplifies T cell activation via an auto-regulatory feedback loop. Selectively targeting this adaptive/stromal interface may provide novel translational strategies in the management of human tendon disorders.

Distinct synovial tissue macrophage subsets regulate inflammation and remission in rheumatoid arthritis.

Immune-regulatory mechanisms of drug-free remission in rheumatoid arthritis (RA) are unknown. We hypothesized that synovial tissue macrophages (STM), which persist in remission, contribute to joint homeostasis. We used single-cell transcriptomics to profile 32,000 STMs and identified phenotypic changes in patients with early/active RA, treatment-refractory/active RA and RA in sustained remission. Each clinical state was characterized by different frequencies of nine discrete phenotypic clusters within four distinct STM subpopulations with diverse homeostatic, regulatory and inflammatory functions. This cellular atlas, combined with deep-phenotypic, spatial and functional analyses of synovial biopsy fluorescent activated cell sorted STMs, revealed two STM subpopulations (MerTKposTREM2high and MerTKposLYVE1pos) with unique remission transcriptomic signatures enriched in negative regulators of inflammation. These STMs were potent producers of inflammation-resolving lipid mediators and induced the repair response of synovial fibroblasts in vitro. A low proportion of MerTKpos STMs in remission was associated with increased risk of disease flare after treatment cessation. Therapeutic modulation of MerTKpos STM subpopulations could therefore be a potential treatment strategy for RA.

Genomic selection for non-key traits in radiata pine when the documented pedigree is corrected using DNA marker information.

BACKGROUND: Non-key traits (NKTs) in radiata pine (Pinus radiata D. Don) refer to traits other than growth, wood density and stiffness, but still of interest to breeders. Branch-cluster frequency, stem straightness, external resin bleeding and internal checking are examples of such traits and are targeted for improvement in radiata pine research programmes. Genomic selection can be conducted before the performance of selection candidates is available so that generation intervals can be reduced. Radiata pine is a species with a long generation interval, which if reduced could significantly increase genetic gain per unit of time. The aim of this study was to evaluate the accuracy and predictive ability of genomic selection and its efficiency over traditional forward selection in radiata pine for the following NKTs: branch-cluster frequency, stem straightness, internal checking, and external resin bleeding. RESULTS: Nine hundred and eighty-eight individuals were genotyped using exome capture genotyping by sequencing (GBS) and 67,168 single nucleotide polymorphisms (SNPs) used to develop genomic estimated breeding values (GEBVs) with genomic best linear unbiased prediction (GBLUP). The documented pedigree was corrected using a subset of 704 SNPs. The percentage of trio parentage confirmed was about 49% and about 50% of parents were re-assigned. The accuracy of GEBVs was 0.55-0.75 when using the documented pedigree and 0.61-0.80 when using the SNP-corrected pedigree. A higher percentage of additive genetic variance was explained and a higher predictive ability was observed when using the SNP-corrected pedigree than using the documented pedigree. With the documented pedigree, genomic selection was similar to traditional forward selection when assuming a generation interval of 17 years, but worse than traditional forward selection when assuming a generation interval of 14 years. After the pedigree was corrected, genomic selection led to 37-115% and 13-77% additional genetic gain over traditional forward selection when generation intervals of 17 years and 14 years were assumed, respectively. CONCLUSION: It was concluded that genomic selection with a pedigree corrected by SNP information was an efficient way of improving non-key traits in radiata pine breeding.

A high-density exome capture genotype-by-sequencing panel for forestry breeding in Pinus radiata.

Development of genome-wide resources for application in genomic selection or genome-wide association studies, in the absence of full reference genomes, present a challenge to the forestry industry, where longer breeding cycles could benefit from the accelerated selection possible through marker-based breeding value predictions. In particular, large conifer megagenomes require a strategy to reduce complexity, whilst ensuring genome-wide coverage is achieved. Using a transcriptome-based reference template, we have successfully developed a high density exome capture genotype-by-sequencing panel for radiata pine (Pinus radiata D.Don), capable of capturing in excess of 80,000 single nucleotide polymorphism (SNP) markers with a minor allele frequency above 0.03 in the population tested. This represents approximately 29,000 gene models from a core set of 48,914 probes. A set of 704 SNP markers capable of pedigree reconstruction and differentiating individual genotypes were tested within two full-sib mapping populations. While as few as 70 markers could reconstruct parentage in almost all cases, the impact of missing genotypes was noticeable in several offspring. Therefore, 60 sets of 110 randomly selected SNP markers were compared for both parentage reconstruction and clone differentiation. The performance in parentage reconstruction showed little variation over 60 iterations. However, there was notable variation in discriminatory power between closely related individuals, indicating a higher density SNP marker panel may be required to elucidate hidden relationships in complex pedigrees.

Approaches to variant discovery for conifer transcriptome sequencing.

There is a wide diversity of bioinformatic tools available for the assembly of next generation sequence and subsequence variant calling to identify genetic markers at scale. Integration of genomics tools such as genomic selection, association studies, pedigree analysis and analysis of genetic diversity, into operational breeding is a goal for New Zealand's most widely planted exotic tree species, Pinus radiata. In the absence of full reference genomes for large megagenomes such as in conifers, RNA sequencing in a range of genotypes and tissue types, offers a rich source of genetic markers for downstream application. We compared nine different assembler and variant calling software combinations in a single transcriptomic library and found that Single Nucleotide Polymorphism (SNPs) discovery could vary by as much as an order of magnitude (8,061 SNPs up to 86,815 SNPs). The assembler with the best realignment of the packages trialled, Trinity, in combination with several variant callers was then applied to a much larger multi-genotype, multi-tissue transcriptome and identified 683,135 in silico SNPs across a predicted 449,951 exons when mapped to the Pinus taeda ver 1.01e reference.

Versatile catechol dioxygenases in Sphingobium scionense WP01T.

The objective was to understand the roles of multiple catechol dioxygenases in the type strain Sphingobium scionense WP01T (Liang and Lloyd-Jones in Int J Syst Evol Microbiol 60:413-416, 2010a) that was isolated from severely contaminated sawmill soil. The dioxygenases were identified by sequencing, examined by determining the substrate specificities of the recombinant enzymes, and by quantifying gene expression following exposure to model priority pollutants. Catechol dioxygenase genes encoding an extradiol xylE and two intradiol dioxygenases catA and clcA that are highly similar to sequences described in other sphingomonads are described in S. scionense WP01T. The distinct substrate specificities determined for the recombinant enzymes confirm the annotated gene functions and suggest different catabolic roles for each enzyme. The role of the three enzymes was evaluated by analysis of enzyme activity in crude cell extracts from cells grown on meta-toluate, benzoate, biphenyl, naphthalene and phenanthrene which revealed the co-induction of each enzyme by different substrates. This was corroborated by quantifying gene expression when cells were induced by biphenyl, naphthalene and pentachlorophenol. It is concluded that the ClcA and XylE enzymes are recruited in pathways that are involved in the degradation of chlorinated aromatic compounds such as pentachlorophenol, the XylE and ClcA enzymes will also play a role in degradation pathways that produce alkylcatechols, while the three enzymes ClcA, XylE and CatA will be simultaneously involved in pathways that generate catechol as a degradation pathway intermediate.

Golgi enrichment and proteomic analysis of developing Pinus radiata xylem by free-flow electrophoresis.

Our understanding of the contribution of Golgi proteins to cell wall and wood formation in any woody plant species is limited. Currently, little Golgi proteomics data exists for wood-forming tissues. In this study, we attempted to address this issue by generating and analyzing Golgi-enriched membrane preparations from developing xylem of compression wood from the conifer Pinus radiata. Developing xylem samples from 3-year-old pine trees were harvested for this purpose at a time of active growth and subjected to a combination of density centrifugation followed by free flow electrophoresis, a surface charge separation technique used in the enrichment of Golgi membranes. This combination of techniques was successful in achieving an approximately 200-fold increase in the activity of the Golgi marker galactan synthase and represents a significant improvement for proteomic analyses of the Golgi from conifers. A total of thirty known Golgi proteins were identified by mass spectrometry including glycosyltransferases from gene families involved in glucomannan and glucuronoxylan biosynthesis. The free flow electrophoresis fractions of enriched Golgi were highly abundant in structural proteins (actin and tubulin) indicating a role for the cytoskeleton during compression wood formation. The mass spectrometry proteomics data associated with this study have been deposited to the ProteomeXchange with identifier PXD000557.

Draft genome sequence of Novosphingobium nitrogenifigens Y88(T).

Novosphingobium nitrogenifigens was originally isolated from pulp and paper mill wastewater, a low-nitrogen, high-carbon environment. N. nitrogenifigens is the first known nitrogen-fixing, polyhydroxyalkanoate-accumulating sphingomonad, and we report the annotated draft genome sequence of the type strain Y88(T) here.

Biodiversity of active and inactive bacteria in the gut flora of wood-feeding huhu beetle larvae (Prionoplus reticularis).

Huhu grubs (Prionoplus reticularis) are wood-feeding beetle larvae endemic to New Zealand and belonging to the family Cerambycidae. Compared to the wood-feeding lower termites, very little is known about the diversity and activity of microorganisms associated with xylophagous cerambycid larvae. To address this, we used pyrosequencing to evaluate the diversity of metabolically active and inactive bacteria in the huhu larval gut. Our estimate, that the gut harbors at least 1,800 phylotypes, is based on 33,420 sequences amplified from genomic DNA and reverse-transcribed RNA. Analysis of genomic DNA- and RNA-derived data sets revealed that 71% of all phylotypes (representing 95% of all sequences) were metabolically active. Rare phylotypes contributed considerably to the richness of the community and were also largely metabolically active, indicating their participation in digestive processes in the gut. The dominant families in the active community (RNA data set) included Acidobacteriaceae (24.3%), Xanthomonadaceae (16.7%), Acetobacteraceae (15.8%), Burkholderiaceae (8.7%), and Enterobacteriaceae (4.1%). The most abundant phylotype comprised 14% of the active community and affiliated with Dyella ginsengisoli (Gammaproteobacteria), suggesting that a Dyella-related organism is a likely symbiont. This study provides new information on the diversity and activity of gut-associated microorganisms that are essential for the digestion of the nutritionally poor diet consumed by wood-feeding larvae. Many huhu gut phylotypes affiliated with insect symbionts or with bacteria present in acidic environments or associated with fungi.