Characterisation of a stable laboratory co-culture of acidophilic nanoorganisms

This study describes the laboratory cultivation of ARMAN (Archaeal Richmond Mine Acidophilic Nanoorganisms). After 2.5 years of successive transfers in an anoxic medium containing ferric sulfate as an electron acceptor, a consortium was attained that is comprised of two members of the order Thermoplasmatales, a member of a proposed ARMAN group, as well as a fungus. The 16S rRNA identity of one archaeon is only 91.6% compared to the most closely related isolate Thermogymnomonas acidicola. Hence, this organism is the first member of a new genus. The enrichment culture is dominated by this microorganism and the ARMAN. The third archaeon in the community seems to be present in minor quantities and has a 100% 16S rRNA identity to the recently isolated Cuniculiplasma divulgatum. The enriched ARMAN species is most probably incapable of sugar metabolism because the key genes for sugar catabolism and anabolism could not be identified in the metagenome. Metatranscriptomic analysis suggests that the TCA cycle funneled with amino acids is the main metabolic pathway used by the archaea of the community. Microscopic analysis revealed that growth of the ARMAN is supported by the formation of cell aggregates. These might enable feeding of the ARMAN by or on other community members.


Results and Discussion
Enrichment process. In 2013, Ziegler et al. described stalactite like biofilms that grow on the ceiling of an abundant pyrite mine in the Harz mountains in Germany 25 . Organisms belonging to the ARMAN as well as members of the Thermoplasmatales were detected within the anoxic core of the biofilms. The goal of this study was to produce a pure or at least highly enriched ARMAN culture to understand potential ecological dependencies of the organisms. Hence, the enrichment process used an anoxic medium containing casein and yeast extract as surrogates for biomass formed by chemolithoautotrophic primary producers of the biofilms. Furthermore, ferric sulfate was added as potential electron acceptor because iron and sulfate were prominent members of the water dripping from the biofilms 24 . Finally, the headspace of the hungate tubes was flushed with a hydrogen/ carbon dioxide mixture (80%/20%) because there was also evidence for potential methanogenic activity in the biofilm, which could be supported by the release of hydrogen from fermentative microorganisms. A mixture of four antibiotics was added to select for archaeal growth. After 1.5 years, bacteria could no longer be observed using CARD-FISH or PCR-analysis. Hence, the addition of antibiotics was omitted from this point onward.
After 2.5 years, the cultures showed a rather stable behavior and could be transferred every 10 weeks. Ferric iron reduction could be observed using the ferrozine assay within the individual growth intervals of the cultures. Moreover, addition of ferric sulfate was necessary for archaeal growth. Therefore, dissimilatory iron reduction seems to be an important trait of at least one member of the archaeal community. At this point, CARD-FISH analysis revealed the presence of archaea, ARMAN and a fungus (Supplementary Figure S1). Of note, the archaea probe does not hybridize to ARMAN species, which indicates the presence of at least one other archaeal species in the enrichment.
Metagenome and phylogenetic analysis. We could reconstruct the genomes of three different archaea.
Hence the CARD-FISH analysis using archaea probes displayed the presence of two different organisms. Supplementary Figure S2 shows the taxonomic profiles for metagenomes determined with Taxator-tk 26 . These profiles corroborate the CARD-FISH results and reveal that the majority of the available sequence (∼80%) can be assigned to the archaea, while the remaining part cannot be assigned to a taxonomic group. Unfortunately, we could not find clear evidence for DNA or RNA that could be assigned to the fungus in the enrichment cultures, possibly because the fungus was underrepresented at the point of sampling and due to the robust cell structure, which might have hampered DNA-isolation. The sequence data suggests a community that consists predominantly of a novel member of the Thermoplasmatales which is most closely related to Thermogymnomonas acidicola (91.6% identity, 1206 bp) and Cuniculiplasma divulgatum (91.7% identity, 1206 bp) (Supplementary Table S1 shows an overview of most closely related sequences). Moreover, an organism belonging to the ARMAN group was identified. Surprisingly, the ARMAN 16S rRNA gene sequence had a 99% identity to a sequence belonging to the ARMAN-1 group that was obtained from samples collected at the Iron Mountains in California (USA) by Baker et al. in ref. 19. A third and minor represented genome was also detected (Fig. 1). The 16S rRNA gene sequence within this part of the metagenome was identical to the published sequence of the recently isolated C. divulgatum 27 and Gplasma, one member of the alphabet-plasmas 28 . Hence, it seems that three archaea and a fungus together form the stable consortium that was isolated over 2.5 years of successive transfers.
In the following sections, we will refer to the ARMAN-1 organism as A_DKE (ARMAN_Drei_Kronen_ und_Ehrt), to the novel member of the Thermoplasmatales as B_DKE and to the organism with the perfect match to the 16S rRNA gene sequence of C. divulgatum as C_DKE (Cuniculiplasma_Drei_Kronen_und_Ehrt). We estimated the completeness for genome bins recovered from the metagenome samples and obtained values of 82% for A_DKE, 99% for B_DKE and 86% for C_DKE (Supplementary Table S2). Of note, all available ARMAN-related genomes known to date have a rather low predicted genome completeness by means of CheckM. This could indicate that ARMANs lack certain genes previously assumed to be present in all archaea. Hence, the completeness of the assembled ARMAN genome is probably higher than 82%. The expected genome sizes of the three organisms are 1, 1.9 and 1.6 Mbp, respectively.
Using the genomic sequences, we sought to determine differences in the occurrences of gene families between A/B/C_DKE and also in comparison to available ARMAN sequences and genomes of six related members of the Thermoplasmatales (Thermoplasma volcanium, T. acidophilum, Thermogymnomonas acidicola, Picrophilus torridus, Ferroplasma acidarmanus and Cuniculiplasma divulgatum S5) (Fig. 2). It is probably not surprising that there were significant differences in the overall number of gene families (Fig. 2a). Moreover, significant differences were detected within a number of functional gene families annotated as translation, amino acid metabolism and metabolism of cofactor and vitamins among others (Fig. 2b, Supplementary Table S3) (based on a two-sided Fisher's exact test, P = 1.9E-8, 3.6E-8 and 5.0E-6, respectively).
A_DKE and comparison to other ARMAN genomes. ARMAN genomes contain on average 1009 + −69 gene families (Fig. 2a). Within the ARMAN genomes, significant differences were detected in several functional gene families annotated as translation, amino acid metabolism and metabolism of cofactor and vitamins among others (Fig. 3b). Within the four ARMAN genomes, we identified 228 gene families that together build the core genome. The A_DKE can be differentiated from the other sequenced ARMAN species by 227 unique gene families compared to 9 unique gene families shared by ARMAN-2, -4 and -5 (Fig. 3c). We found that 15 of these 227 unique gene families could be successfully annotated (Supplementary Table S5). Compared to the other three ARMAN species, the genome of A_DKE shows more genes related to benzoate degradation, the phosphotransferase system, the metabolism of aromatic compounds as well as amino acid metabolism and degradation (Fig. 3b). Based on the 16S rRNA genes, A_DKE seems to be most closely related to ARMAN-2. This was also corroborated by further comparative genome analysis (Fig. 3). The GC-content of ARMAN-2 is rather similar to A_DKE, while ARMAN-4 and -5 share lower GC values. Furthermore, the number of similar genes -defined by expected E-values lower 1e-100 -also seems to cluster A_DKE with ARMAN-2 while ARMAN-4 and ARMAN-5 build via that analysis a second group (Fig. 3a). Using an automated assignment of genes to central metabolic pathways using the KEGG database we could highlight the presence and similarity within and between the different ARMAN genomes. It seems that all organisms share a similar set of genes involved in the citric acid cycle and in respiration. In contrary, genes involved in glycolysis/gluconeogenesis, the pentose phosphate pathway, the Enter-Doudoroff pathway or the ß-oxidation are not as similar and wide spread within the genomes of the organisms.
Overall, the genome of A_DKE differs from previously sequenced ARMAN species in terms of sugar metabolism. Baker et al. 18 described a pentose phosphate pathway for ARMAN-4 and -5 and an incomplete glycolysis for ARMAN-2 18 . In contrast, we could not identify genes encoding the necessary enzymes for the conversion of glucose to pyruvate. Moreover, we could identify only two potential genes that could be involved in the pentose phosphate   2) ORF density (annotated and hypothetical) within a sliding window over the genomes visualized in black. ORF density of hypothetical ORFs visualized in red based on a sliding window approach. (3)(4)(5)(6)(7)(8) Matches to KEGG terms in this order: beta-Oxidation (yellow), Oxidative phosphorylation (green), Entner-Doudoroff pathway (pink), pentose phosphate pathway (purple), citric acid cycle (orange), glycolysis/gluconeogenesis (blue). Lines within the inner segment indicate significant similarity between ORFs of the individual genomes and are colored according to KEGG terms if they fully overlap these. Contig sizes are visualized in Mb. (b) Bar plot showing the over-abundance of functions based on de novo inferred gene families based on orthoMCL analysis of three ARMAN genomes and the A_DKE genome. The length of the bars corresponds to the mean difference in the fraction of gene families annotated with the indicated functional group. Lines show the difference to the ARMAN genomes with the minimal and the maximal fraction. Only functional groups are shown that differ in more than 0.2% or less than −0.2%. (c) Venn diagram of gene families shared between indicated genomes. Gene families related to pathogenicity factors were not included in this overview. pathway. Moreover, the genes encoding for key reactions of the ED pathway were not detectable (Supplementary  Table S6). Hence, there seems to be a necessity for the uptake of hexose and pentose sugars from the environment. Furthermore, a ß-oxidation of fatty acids seems to be not possible, since the genome contains only genes for a 3-hydroxyacyl-CoA dehydrogenase and an acetyl-CoA thiolase. In contrary, we detected genes for a complete tricarboxylic acid cycle (TCA), with the exception that bioinformatic analysis revealed so far only one subunit of the succinyl-CoA synthetase in A_DKE. Archaeal electron transport complexes can have unusual compositions 29 . Nevertheless, bioinformatic analysis predicts so far only the existence of a NADH-dehydrogenase, a respiratory succinate-dehydrogenase and an ATP-synthase in the ARMAN genome. Hence, it is not clear how the energy metabolism of the organism might function. While the citric acid cycle would be a way to produce ATP/GTP, it also leads to the production of reducing-equivalents. How these reducing equivalents are transported to an electron acceptor and what this electron acceptor might be is not clear so far. Our data suggest an anoxic life style. Evidence for fermentative pathways was not found. Hence, we must assume that the organism can use either ferric iron and/ or sulfate as electron acceptor (the only electron acceptors present in the medium) or that it can conduct an interspecies electron transfer to one of the other community members. Therefore, we must also assume that ARMAN species might have evolved novel components for respiratory electron transfer that are currently not detectable by bioinformatics analyses and will necessitate biochemical analysis for characterisation.

B_ and C_DKE.
A dendrogram based upon the hierarchical clustering of gene families shared between genomes of five related members of the Thermoplasmatales (Fig. 2c) reveals that B_DKE and C_DKE cluster together with T. acidicola. The core genomes of the selected Thermoplasmatales seem to consist of 624 gene families while the overall average number of detected gene families was 1,935 + −557. A search for unique gene families (Fig. 4c, Supplementary Table S7) in B_DKE compared to C_DKE revealed 914 gene families. Most of these were annotated as ABC transporters, two-component systems or as involved in the oxidative phosphorylation (more abundant in B_DKE) (Fig. 4b). Comparative genome analysis corroborates the identity of the 16S rRNA genes of C. divulgatum and C_DKE. This similarity can also be observed on the genome level since most of the genes of C. dilvigatum have a similar counterpart in C_DKE (Fig. 4a). In contrary, the genome of B_DKE contains several areas with a rather low number of similar genes compared to C_DKE, C. divulgatum and T. acidicola. As expected, the genomes of the 4 members of the Thermoplasmatales have overall a low percentage of genes with hypothetical function compared to the ARMAN genomes (Fig. 3a).
Both B_ and C_DKE contain all necessary enzymes for pentose formation via the non-oxidative pentose phosphate pathway similar to other members of Thermoplasmatales 30 (Supplementary Table S6). Furthermore, the genome of B_DKE contains nearly all enzymes for the branched ED pathway including the key enzyme bifunctional KDG-/KDPG-aldolase. Nevertheless, we could not find a corresponding gene to a KDG-kinase. The genome of C_DKE shows all enzymes of the non-phosphorylated variant of the ED except its key enzyme, the KDG-aldolase. Nevertheless, a complete non-phosphorylative ED was detected in the genome of Gplasma 28 and other members of the Thermoplasmatales like Thermoplasma acidophilum 31 or Picrophilus torridus 32 . Hence, this potential lack of pathways for the catabolism of sugar compounds is not a general characteristic of similar organisms.
Sugar degradation catalyzed by the two Thermoplasmatales members via glycolysis can be excluded since both genomes lack evidence for an encoded phosphofructokinase, which corresponds well with the genomic data of C. divulgatum 33 . Still, B_DKE and C_DKE show all enzymes for gluconeogenesis, with the exception of the fructose 1,6-bisphosphate (FBP) aldolase. Nevertheless, a corresponding gene for this enzyme is missing in most archaeal genomes including Thermoplasma acidophilum 34 and Picrophilus torridus 35 . Results by Say and Fuchs (2010) revealed the presence of a bifunctional FBP aldolase-phosphatase with high FBP aldolase and FBP phosphatase activity, which guarantees an unidirectional gluconeogenesis pathway 36 . This exists in nearly all archaeal groups and may be the putative ancestral gluconeogenic enzyme. A corresponding gene was found in the genomes of B_DKE and C_DKE, which likely allows both organisms to produce sugar compounds. The key enzyme for glycolysis, the phosphofructokinase, is missing also in the genomic data of Gplasma, whereas a bifunctional FBP aldolase-phosphatase is detectable 28 .
B_ and C_DKE contain the genes for a complete tricarboxylic acid cycle (TCA). In contrast, genomic analyses of C. divulgatum PM4 and S5 revealed that 2-oxoglutarate dehydrogenase, fumarate reductase and fumarase might be missing 33 . Nevertheless, Golyshina et al. point out that these enzyme functions could be functionally replaced by enzymes catalyzing closely related reactions that were predicted by bioinformatic analysis. The presence of succinate dehydrogenase/fumarate reductase and 2-oxoglutarate/2-oxoacid ferredoxin oxidoreductase could also be evidence for the presence of a reverse TCA cycle in the three enriched archaea. Still, the ATP citrate lyase-the third indicator gene for the reverse TCA-is missing in all three genomes. Similar to the DKE_A genome, C_DKE and B_DKE do not seem to possess a complete ß-oxidation pathway. While a ß-hydroxyacyl-CoA dehydrogenase was not detectable in C_DKE, B_DKE does not seem to have ß-hydroxyacyl-CoA dehydratase and dehydrogenase genes.
Bioinformatic analysis predicts also in B_ and C_DKE the existence of a NADH-dehydrogenase, a respiratory succinate-dehydrogenase and an ATP-synthase. Furthermore, B_DKE and C_DKE contain corresponding genes for a cytochrome c-oxidase although genetic evidence for a cytochrome c protein is lacking. Both contain a Rieske FeS protein that is a characteristic component of complex III and it forms a gene cluster in B_DKE with a cytochrome b6-like protein that could have complex III function 29 . In addition, at least B_DKE has a gene which most probably encodes a copper protein belonging to the plastocyanin/azurin family, which could functionally replace cytochrome c 29 . We currently do not know how this electron transport core chain might be connected to ferric iron reduction, which is catalyzed at least by B_DKE within the community (see below). Furthermore, we detected in C_DKE only a potential part of complex III and no cytochrome c gene. Hence, it is unknown whether the latter might be a remnant of genomic reduction and selection towards a shorter electron transport chain involving only complex I and a cytochrome ubiquinol oxidase. Metatranscriptome. The ratio of sequence reads that can be referred to B_DKE, A_DKE and C_DKE is 1:0.5:0.27. Hence, all three archaea considerably contribute to the overall metatranscriptome of the study. The high content of RNA sequences that can be referred to A_DKE is surprising because all previous studies on ARMAN revealed a remarkably small number of ribosomes and because -as will be shown later -the number of ARMAN cells within the community is lower than the sum of B_DKE und C_DKE. We suggest that the considerably high content of A_DKE RNA within the metatranscriptome is evidence for active growth of the majority of ARMAN cells at sampling.
In A_DKE, 14 of the 20 most highly expressed genes encode for hypothetical proteins, which hampers the prediction of potential metabolic reactions that sustain growth of the organism. Nevertheless, the data suggests that at least glutamate could be a carbon and energy source for the organism. The glutamate dehydrogenase gene (gdh) is relatively high expressed (Fig. 5). The corresponding enzyme catalyzes the conversion of glutamate to 2-oxoglutarate, an intermediate of the citric acid cycle. Furthermore, transcripts could be identified for all detected genes of the citric acid cycle. Furthermore, complex I and V of the respiratory chain are expressed.
Similarly to A_DKE, the gdh gene is also highly expressed in B_DKE and C_DKE. In agreement with this, it can be observed in both transcriptomes that genes for TCA cycle enzymes are more highly expressed relative to the ED pathway (Fig. 5). The question about the variant of the ED remains unclear, because B_DKE expresses the 2-dehydro-3-deoxy-(phospho)gluconate (KD(P)G) aldolase as well as all enzymes for the non-phosphorylated variant of the ED and the glyceraldehyde-3-phosphate dehydrogenase (GAPN) that is involved in the semi-phosphorylated variant of the pathway and could not be detected in other members of the Thermoplasmatales 30 . The C_DKE and B_DKE also express the complete set of genes necessary for gluconeogenesis and the pentose phosphate pathway. Furthermore, all detected genes encoding for electron transfer processes are highly expressed.
An analysis of the 20 most highly expressed genes (Supplementary Table S8) showed that both Thermoplasmatales species have high expression rates for oxidative stress proteins as well as proteins linked to electron transport process. Glutamate dehydrogenase and 2-oxoacid ferredoxin oxidoreductase are also highly expressed in C_DKE. The detected high expression rates for oxidative stress proteins could be due to the presence of trace amounts of oxygen during the cultivation, which would also indicate that oxygen could be a potential electron acceptor for the organisms. Hence, growth experiments were conducted within a glove box containing a 95% N 2 /5% H 2 atmosphere. If it all, we saw slower growth of the fungus within the enrichment culture. The archaea did not seem to be affected by the strictly anoxic conditions (data not shown).

Enrichment of B_DKE and isolation of the fungus.
Interestingly, omitting the addition of antibiotics lead in some cultures to an increased growth rate of B_DKE. This resulted in cultures containing only B_DKE and the fungus. The fungus was isolated from these cultures using a solidified medium described by Baker et al. 37 . Fungal colonies on these plates were used for phylogenetic analysis using a primer set for ITS amplification. The phylogenetic classification of the 5.8S rRNA gene and relating ITS-DNA sequences of the fungus revealed a 99% identity to Acidothrix acidophila, which is an acidophilic fungus isolated from acidic soils in the Czech Republic 38 . Of note, growth of the fungus in the medium used for the enrichment process did not lead to a reduction of the ferric sulfate. In contrary, growth of a culture containing only the fungus and B_DKE was connected to ferric sulfate reduction. Hence, B_DKE is capable of ferric iron reduction.
Timeline experiments. The dynamics and interactions within the enrichment cultures were analyzed via timeline experiments. Hence, triplicate cultures were analyzed over eleven weeks of growth. The cultures were fixed for CARD-FISH analyses twice weekly, and DNA was isolated for qPCR-based quantitative community analyses once a week. Furthermore, samples were taken for measurements of ferric iron reduction as well as control of the pH values every week. Of note, as was emphasized before, we do not have genetic information regarding the fungus so far. Therefore, we could not conclude with certainty from gene to cell quantities. Hence, the qPCR analysis was conducted only with the archaeal members of the consortium.
Analysis of the qPCR and CARD-FISH data revealed a rather similar growth behavior of A_DKE as well as B_DKE and C_DKE (Fig. 6). In the lag phase, the cell numbers increased continuously with a slightly higher amount of A_DKE cells compared to that of B_DKE and C_DKE. By week six, all organisms seemed to enter the logarithmic growth phase. From this point on, B_DKE and C_DKE began to dominate the cultures. The highest cell concentrations were reached in week eight. Thereafter, a rapid decrease in cell numbers was detected. The CARD-FISH pictures also highlight that A_DKE cells are mostly part of B_DKE or C_DKE cell agglomerates. Assuming that agglomerates are formed by extracellular polymeric material, which usually consists to a large fraction out of sugar compounds, we speculated that the biofilm growth might ensure the access to sugar for A_DKE. We emphasized in the previous section that some cultures showed faster growth during the conducted transfers and that these cultures were composed only of B_DKE and the fungus. These cultures do not show growth in the form of flocks but rather are uniform planktonic cells. Hence, the change in the growth phenotype might be the reason for decreasing A_DKE cell numbers. Supplementary Figure S3 displays the difference between the two growth phenotypes of B_DKE.
Interestingly, the ferrous iron concentrations as well as the measured pH values correlate with the growth curve data (Fig. 6). Until week six, only a slightly increased content of ferrous iron could be measured, and the pH values remained constant. From week seven on, there was a rapid ferric iron reduction along with increased growth rates. Even the pH values showed a slight decrease from this point on representing the acidification associated with the ferric iron reduction.

Conclusion
We established a highly enriched consortium consisting of two acidophilic Thermoplasmatales species, an ARMAN, and a fungus. Of note, this is the only stable enrichment culture described so far containing ARMAN organisms. During the course of writing this manuscript, it was possible to obtain a culture that is composed only of A_DKE, B_DKE and the fungus (data not shown). Hence, we cannot say whether a specific archaeon is necessary to sustain growth of the ARMAN, but at least the interaction with B_DKE seems to be sufficient.
The cultivated ARMAN-1 has a complete or nearly complete TCA cycle and lacks nearly all enzymes necessary for glycolysis/glyconeogenesis and the pentose phosphate pathway. Moreover, we could not find evidence for gluconeogenesis and the pentosephosphate pathway, suggesting a dependency on synthesis of hexose and pentose sugars by other community members. This dependency seems to be corroborated by the identification of ARMAN cells within biofilm flocks that are likely built primarily by B_DKE. The colocalization of the cells might simplify access to suitable carbohydrates for the ARMAN. At least B_DKE seem to thrive using a respiratory metabolism with ferric iron as electron acceptor. Very recently, the sequence of another ARMAN-1 genome was published 39 . It will be interesting to study differences in the genomes obtained from the US and German field sites in future studies.
We tried to exclude oxygen from the growth medium and measurements of the oxygen content of the cultures by an optode suggest an anaerobic lifestyle. Moreover, similar growth of the archaea was also observed in an anoxic glove box. The reatively high expression levels of genes corresponding to enzymes involved in the elimination of oxidative stress might be due to the relatively short exposure to oxygen at the point of RNA sampling. In future experiments, we will combine heat inactivated B_DKE/fungus extracts with filtration based selection for pure ARMAN cultures. Nevertheless, B_DKE cells are also very small and pleomorphic, and the doubling time of the organisms is rather low. Hence, although future experiments seem straight forward they will require months or years of incubation time for several necessary culture transfers.

Experimental Procedures
Culturing conditions. Enrichment cultures of Thermoplasma spp. and ARMAN-1 originate from acidophilic stalactite-like biofilms of the former pyrite-mine "Drei Kronen und Ehrt" in the Harz Mountains, Germany (Ziegler et al. 24,25 ). Biofilms were inoculated under anoxic conditions in Hungate tubes and in a medium originally designed for the isolation of Picrophilus species 40 . The pH of the medium was adjusted to 2.5 using 0.5 M H 2 SO 4 . Afterwards, the medium was autoclaved and supplemented with sterile filtered solutions of yeast extract, casein and ferric sulfate to final concentrations of 0.1% each and 20 mM, respectively, as critical additives. The growth of the bacterial species was repressed with 150 µg ml −1 streptomycin, 50 µg ml −1 kanamycin, 30 µg ml −1 chloramphenicol and 2 µg ml −1 vancomycin. The headspace was flushed with an 80% H 2 /CO 2 20% gas mixture. Cultures were inoculated with 10% of a growing co-culture and incubated at 22 °C. Escherichia coli cells were routinely cultured in LB medium supplemented with 40 µg ml −1 kanamycin or 100 µg ml −1 ampicillin, if necessary. Isolation of a fungus was conducted with a medium based on the study by Baker et al. 37 . The pH of the medium was adjusted to 2 using 0.5 M H 2 SO 4 and it was supplemented with SL10 trace elements (DSMZ medium 320). The media plates were prepared with the addition of 2% (w/v) agar.
Ferrous iron quantification. Ferric iron reduction was determined by quantification of the ferrous iron content spectrophotometrically using the ferrozine assay as described previously 43 .

Isolation of fungal DNA and amplification.
A protocol based on Cenis (1992) 44 was used to isolate DNA from the fungus. For taxonomic classification, the internal transcribed spacer (ITS) region primers ITS1 and ITS4 45 were used to amplify the DNA with the iProof ™ High-Fidelity PCR Kit (Bio-Rad, Munich, Germany).

Quantification of cells using quantitative PCR.
Quantitative PCR was used to quantify growth of the archaeal species in the enrichment culture. Standard curves were developed by integrating plasmids containing the target sequences of the qPCR primer sets in the genome of E. coli. Primer sets were developed for the 23S rRNA gene of ARMAN-1 and the 23S rRNA of the two Thermoplasma spp. using NCBI/Primer-BLAST 46 as described in the supplemental information. The primer sequences are listed in Supplementary Table S9.
The organisms could contain more than one copy of the 23S rRNA gene, and this would hamper quantitative analysis. Thus, the normalized per-sample coverage of all genes was compared to the 23S rRNA gene coverage. The latter was always below the overall coverage (Supplementary Table S10). Hence, the assumption was made that the 23S rRNA was represented by one copy in the genome. DNA/RNA isolation for metagenomics/metatranscriptomics. Isolation of genomic DNA for metagenomic analysis was conducted according to Lo et al. 47 . The total RNA isolation was conducted according to the protocol of the TRIzol Max Bacterial RNA Isolation Kit (Thermo Fisher, Waltham, USA) with 8 ml of an enrichment culture as the starting material. The remaining DNA was hydrolyzed with the DNA-free ™ DNA Removal Kit (Thermo Fisher, Waltham, USA) following the manufacturer´s instructions.
Metagenomic/metatranscriptomic sequencing and analyses. Metagenomic DNA from two analogous enrichment cultures was used for Illumina-based sequencing using a 2 × 51 bp protocol. The metatranscriptomic analysis was based on an RNA isolation from 8 mL of the enrichment culture. This sample was pooled from analogous enrichment cultures with similar growth characteristics. The RNA was sequenced using the 2 × 100 bp protocol. Supplementary Table S11 shows the overall statistics of the sequencing experiments. Data availability. All DNA and RNA sequences that were retrieved within this study are publically available through NCBI BioProject: PRJNA358824.
16S rDNA gene sequence of B_DKE and 18S rDNA/ITS gene sequence of the fungus are available through NCBI GenBank accession numbers KY825129 and KY848528.