INTRODUCTION
The Orchidaceae is a cosmopolitan plant family that has its highest diversity in the tropics (Elórtegui & Novoa, 2008, Heywood et al. 2007). Tropical species are usually epiphytes (Harley & Smith 1983, Jones 1993, 2006), whereas temperate species are usually terrestrial or rupicolous (Cribb et al. 2003, Dressler 1993, Valadares et al. 2012). In Chile, only terrestrial orchid species are found, distributed in eight genera (Novoa et al. 2015). Most Chilean species are endemic, and they can be found from the northern-most region to the Magallanes Region in the south ranging from ∼18°S to 55°S with a gap in the Atacama Desert (Novoa et al. 2006, 2015). The genus Codonorchis includes two species; C. canisioi Mansf. found in Brazil, and the native C. lessonii (Brongn.) Lindl. (Freuler 2006, Novoa et al. 2015). C. lessonii is distributed only in Chile and Argentina, and is usually found in the understory of temperate forests in the central-southern part of South America (Romero 2012, Novoa et al. 2015). The phytogeographic origin of C. lessonii is still unclear and currently there are no proven theories to explain its presence in Chile where no close relatives are found among other Chilean orchid species (Novoa et al. 2015). This species is related to species of the Orchideae tribe, which include plants that grow mostly in Africa and the Old World (Elórtegui & Novoa 2008, Novoa et al. 2015, Szlachetko & Tukallo 2008). For many years it was associated to genera belonging to the Chloraeinae sub-tribe like Chloraea and Bipinnula (Elórtegui & Novoa, 2008, Freuler 2006, Szlachetko & Tukallo 2008), however it has been recently placed in the tribe Codonorchideae (Givnish et al. 2015).
Orchids reproduce sexually through seeds that once released depend on mycorrhizal fungi to germinate and sustain early embryo growth (Murguía & Lee 2007, Pereira et al. 2005a, Smith & Read 2008). This dependency is due to the lack of developed endosperm in orchid seeds. The fungi supply relevant nutrients to the embryo such as carbon, nitrogen, and phosphorus (Arditti & Ghani 2000, Smith & Read 2008). Many orchid species establish a close association with these fungi, sometimes considered parasitic, but the energetic cost of the interaction may be relatively low for the fungal partner (McKendrick et al. 2000). At adult stages, many species maintain symbiotic associations with fungi that colonize the cortical cells of their roots (Pereira et al. 2014, Atala et al. 2015, Herrera et al. 2017). However, the spectrum of fungal partners may change from germination to adult stages of an orchid (Bidartondo & Read 2008). Most of the fungal partners found associated with autotrophic orchids belong to or show Rhizoctonia traits (Atala et al. 2015, Otero & Bayman 2009, Mosquera-Espinosa et al. 2010, Pereira et al. 2014, Herrera et al. 2017). However, other Basidiomycetes belonging to Russulaceae and Thelephoraceae, and some Ascomycetes like Tuber spp. have been also found associated with heterotrophic orchids (Dearnaley et al. 2012, Karol et al. 2015, Selosse et al. 2004). Members of the Rhizoctonia complex do not form asexual spores and share certain distinctive vegetative characters (González et al. 2006, Nogueira et al. 2014). The Rhizoctonia complex includes four fungal teleomorphs: Ceratobasidium, Tanatephorus, Tulasnella and Sebacina (Agrios 2002, Goh et al. 1992, Mosquera-Espinosa et al. 2010). Moore (1987) proposed the subdivision of the complex into Ceratorhiza, Epulorhiza and Moniliopsis, based on the number of nuclei per hyphal segment and the ultrastructure of the hyphal septum. Other authors, however, have proposed its division into only two asexual genera, Ceratorhiza (including Ceratobasidium and Thanatephorus) and Epulorhiza (including Tulasnella and Sebacina) (Pereira et al. 2005b, González et al. 2006, Dearnaley et al. 2012). In the last decade, the number of studies on mycorrhizal fungi associated with tropical and terrestrial temperate orchids has increased significantly (Otero & Bayman, 2009, Otero et al. 2004, 2005, Pereira et al. 2005a, 2005b, Ruibal et al. 2017, Schiebold et al. 2017). However, this is not the case for South American terrestrial orchids, particularly temperate species, where studies are relatively scarce and have focused on a few genera like Bipinnula and Chloraea (Pereira et al. 2014, 2015, 2017, Atala et al. 2015). Recently, Fracchia et al. (2014) isolated a Tullasnella species from Codonorchis lessonii in Argentina and tested it in the symbiotic germination of Gavilea australis. Understanding the fungal association in orchids is of great importance, as fungal availability in the environment may play a key role in orchid distribution and diversity (Brundrett 2007, McCormick et al. 2012, Otero et al. 2002). On the other hand, mycorrhizal fungi can help with the efficient propagation of terrestrial orchids and their re-introduction into natural habitats (Brundrett 2007).
These orchid-fungi associations seem to be specific to some degree, with specificity being the phylogenetic breadth of fungal partners that a plant taxon may be associated with (Otero et al. 2004, Otero & Bayman 2009, Taylor et al. 2002). One orchid species can associate preferentially with a small range of phylogenetically close fungal partners, while in other cases, with just one particular fungal species or strain. Moreover, specificity can result in increased germination or plant growth when associating with a preferred fungus (Otero & Bayman 2009). Except for mycoheterotrophs, high specificity in orchid-fungi associations is not common in temperate or boreal areas (Zettler 1997, Zettler & Hofer 1998, Herrera et al. 2017). In Australia, however, temperate orchid-fungi associations are highly specific, often with one fungal strain associating with only one orchid species (Brundrett 2007, Dixon 1991). There are, nevertheless, other Australian species that associate with a wider diversity of fungi (Brundrett 2007). The degree of specificity may be related to their trophic strategies such as autotrophy or mycoheterotrophy, as well as to abiotic conditions (Otero & Bayman 2009).
To our knowledge, there is no information about the degrees of specificity in the plant-fungal associations of Chilean orchid species (but see Herrera et al. 2017), particularly for C. lessonii, which could have a significant impact on their conservation (Ramsay & Dixon 2003). Reduced distribution or ecological niches of highly specific orchids could be due to a limited distribution or availability of their fungal partners (Leake & Cameron 2012). Moreover, it has been suggested that a dependence on narrow specific interactions with fungi and pollinators may predispose many orchids to become rare (Dearnaley 2007, Swarts et al. 2010). For in situ conservation and restoration initiatives it is likely that specific fungi must also be reintroduced in order to have successful establishment of adult plants (Ramsay & Dixon 2003).
The identification and isolation of fungi associated with native and endemic Chilean orchids could help in the conservation and ecological restoration of orchid species, many of which have known conservation problems (Novoa et al. 2015, Atala et al. 2015, 2017,). It may also contribute to the improvement and development of Chilean orchids as ornamental species. The present study aims to isolate, characterize, and identify the fungal species symbiotically associated with adult individuals of Codonorchis lessonii growing in three populations in Central-South Chile.
MATERIALS AND METHODS
FIELD SAMPLING AND MORPHOLOGICAL IDENTIFICATION OF THE FUNGAL ISOLATES
Individuals of Codonorchis lessonii (Brongn.) Lindl. Were collected from three populations in Central-South Chile; 1) Los Guindos sector, road to Antuco Volcano (37º 21’ 31.26” S - 71º 51’ 33.28” W), 2) Las Chilcas sector, Parque Nacional Laguna del Laja (37º 23’ 20.37’’ S - 71º 24’ 22.10’’ W), and 3) Parque Nacional Villarrica (37º 23’ 20.37’’ S - 71º 24’ 22.10’’ W). Populations were located over 35 km apart from each other. The identification of the species in the field was based on its floral morphology (Freuler 2006, Mourgues 2015, Novoa et al. 2015), very distinctive among Chilean orchids (Novoa et al. 2015, Romero 2012). Three plants per population at least 1 m away from each other were sampled (a total of 9 individuals). Sampled plants ranged from 15 to 20 cm tall, considering the flower. For every sampled individual, roots were collected, stored in sealed plastic bags and transported to the laboratory. Roots were rinsed with tap water mixed with Tween 20. Transversal root segments of approximately 2 mm were cut and sterilized by serial ethanol baths (Otero & Bayman 2009, Pereira et al. 2014). Roots segments were examined for the presence of pelotons with a dissecting microscope at a 4X magnification, measuring the % of the total root section occupied by fungi following previously published methods (Bertolini et al. 2014). Root sections that showed presence of hyphae coils (pelotons) in their cells were placed in Petri dishes with autoclaved potato dextrose agar (PDA) medium at pH 5.8. The resulting plates were then placed in an oven and incubated in the dark at 24 ± 1 ºC for 10 days. Growing colonies were observed with a light microscope (Olympus CX30) at 40X and 100X magnification. Hyphal traits such as septum type, branching angle, constrictions of hyphae, and presence of monilioid cells were recorded and photographed using a digital camera attached to a microscope (Moticam 2000). Teleomorph identification was based on the number of nuclei per cell (Agrios 2002, Sneh et al. 1991, Valadares et al. 2012), by measuring 20 young fungal cells per plate. Radial growth was also determined in active fungal colonies isolated from the three orchid populations, averaging the obtained values from four replicates per isolate. Surface-sterilized root segments and mycelial isolates from randomly selected colonies sharing similar characteristics, were chosen from each plant and were individually stored in 300 μl of 2% CTAB at -20 ºC until DNA analysis.
MOLECULAR IDENTIFICATION OF THE FUNGAL ISOLATES
Genomic DNA from individual root sections and from 0.25 cm2 fragments of mycelia collected from the active growing edge of the colonies in PDA medium was extracted following a CTAB protocol (Gardes & Bruns 1993) using a GeneClean kit (Qbiogene) for DNA purification (Bidartondo et al. 2004). The ITS region of the nuclear ribosomal DNA was amplified using three primer pairs, the eukaryotic one ITS1-ITS4, the fungal specific ITS1F-ITS4, and ITS1-ITS4-tul, specific for the Tulasnellaceae family (Gardes & Bruns 1993, Taylor & McCormik 2008, White et al. 1990). PCR products were then purified using ExoSAP- it (USB) and bidirectionally sequenced using BigDye 3.1 in an ABI 3730 sequencer (Applied Biosystems).
Sequences were edited using Sequencher v.4.2 (GeneCodes, Ann Arbor, MI) and deposited in GenBank under the accession numbers KT003599 to KT003606 (Table 1). Edited ITS sequences were compared with available sequences in GenBank through BLAST searches (Altschul et al. 1997) to determine their closest known relatives. The top 50 sequences producing the closest hits with each sequence obtained in this study were included in the dataset for phylogenetic analysis. Sequences were aligned and trimmed using the RNA structure-based algorithm Q-INS-i implemented in MAFFT v7.221 (Katoh 2014). Phylogenetic analysis was performed under the maximum likelihood criterion using RAxML (Stamatakis 2014) implemented in raxmlGUI v1.3.1 (Silvestro & Michalak 2012). The GTRGAMMA model of evolution was used and branch support was assessed using 1000 nonparametric bootstrap replicates. Trees were midpoint rooted.
DATA ANALYSIS
The software Statistica v.6.0 (StatSoft Inc, OK, USA) was used for statistical analyses. Differences in fungal growth rate (Montgomery 1991) were tested using ANOVA and an a posteriori Tukey test.
RESULTS
Floral traits confirmed that the sampled individuals corresponded to Codonorchis lessonii in all studied populations (Fig. 1a). Root colonization through root hairs by septate hyphae (Fig. 1b, 1c) and pelotons inside cortical cells in all sampled individuals were observed. Over 85% of the cell area in the roots was occupied by hyphal coils (Fig. 1d, 1e).
Microscopical observations of the fungal cultures revealed that individuals of the three sampled populations of C. lessonii associated with fungi showing Rhizoctonia- like traits. The hyphae showed strait angle ramification, simple septa, constriction at the base of a branching hypha, and presence of monilioid cells (Fig. 1g, 1h).
Molecular analyses allowed the identification of the fungal pelotons in the root sections and confirmed the preliminary morphological identifications of the fungal isolates. Based on the GenBank searches, three different fungi belonging to Ceratobasidiaceae (from Antuco and P. N. Villarrica) and one to Tulasnellaceae (P. N. Laguna del Laja) were found (Table 1).
Phylogenetic analyses showed the placement of the sequences with closest matches to Ceratobasidiaceae sequences in three separated, well-supported clades (Fig. 2) and confirmed the identification of the isolates KT003603 and KT003604 as Tulasnellaceae sp. (Fig. 3).
Morphological traits and growth of isolated fungi associated with the roots of C. lessonii differed between the orchid populations from which the fungal strains were isolated (Table 2). Isolated strains differed in peloton size within plant cells, hyphae cell size, and monilioid cells size. All isolates showed monilioid cells in pure culture (Fig. 1f, 1h), which also varied in shape and frequency. We observed two nuclei per cell in all isolates (Fig. 1i,). Additionally, isolates identified as Ceratobasidiaceae showed a higher growth rate compared to isolates identified as Tulasnellaceae (Table 2).
DISCUSSION
We found mycorrhizal fungi associated with the roots of all studied individuals of C. lessonii. Intact pelotons were evident in the root cells of all plant populations, occupying over 85% of the total root section area. This, together with their poor root development (Fig. 1), suggest a close dependence of C. lessonii on these mycorrhizal fungi (Novoa et al. 2015). Intact pelotons are a source of inoculum that allows for the re-colonization by hyphae of the outer root surface (Pereira et al. 2005a, 2005b, Sathiyadash et al. 2012) and could be a key strategy for infection of new orchid seedlings in the field. Both, the beneficial mycorrhizal association in terms of plant fitness, and its role as a source of inoculum for new orchid seedlings, could explain the small high-density patches formed by this orchid in the field (personal observation).
Orchids frequently associate with Basidiomycetes belonging to the form-genus Rhizoctonia (Rasmussen 1995, Brundrett 2007, Otero et al. 2013, Pereira et al. 2014, Herrera et al. 2017). These fungi include saprotrophs, ectomycorrhizal fungi, and some parasites and plant pathogens (Dearnaley et al. 2012, Karol et al. 2015). The fungal isolates from all populations of C. lessonii showed binucleated young cells and, according to their morphological characteristics, they were preliminarily classified as belonging to Ceratobasidium (Ceratorhiza) and Tulasnella (Epulorhiza). Morphological features are important in classic taxonomy to identify genera and species of fungi. Orchid mycorrhizae have been studied for more than a century (Rasmussen 2002), yet the identification of the fungal species was severely limited until DNA sequencing techniques became more commonly used (McCormick & Jacquemyn 2014). Molecular techniques are required for greater precision, especially for species within the Rhizoctonia complex (Mosquera-Espinosa et al. 2010, Pereira et al. 2011, Valadares et al. 2012). Both methodological approaches have been used in this study, supporting our conclusions.
Phylogenetic analyses showed that three different fungal strains belonging to Ceratobasidiaceae and one to Tulasnellaceae associate with C. lessonii.Fracchia et al. (2014) identified Tulasnella calospora associated with the roots of an Argentinian population of C. lessonii. Similar associations with different Ceratobasidium and Tulasnella taxa had also been reported for other endemic Chilean orchids such as Gavilea araucana (Duran et al. 2007), Bipinnula fimbriata (Steinfort et al. 2010), Chloraea collicensis and C. gavilu (Pereira et al. 2014, Herrera et al. 2017), and C. cuneata (Atala et al. 2015). Fungi isolated from C. lessonii showed lower growth rates than isolates obtained from roots of other endemic orchids in Chile (Pereira et al. 2014, Steinfort et al. 2010) and the growth rate differed among populations. Further studies of the efficiency of different fungal strains in seed germination and seedling recruitment of this orchid species are needed.

FIGURE 1 a. Codonorchis lessonii individual. b. Roots of C. lessonii, the arrow indicates areas used for the fungal isolation. c. Hyphae entering the root through a root hair. d. Transversal section of a root showing pelotons inside parenchyma cells. e. Pelotons inside parenchyma cells. f. Isolated fungal colonies. g. Simple septum and 90° angle hyphal branching of Ceratobasidium. h. Monilioid cells of a Ceratobasidium isolate. i. Binucleated cells of members of Tulasnellaceae. c, e, g, h and i taken at 100X magnification. d taken at 40X magnification. / a. Individuo de Codonorchis lessonii. b. Raíces de C. lessonii, la flecha indica las áreas utilizadas para el aislamiento fúngico. c. Hifa entrando en la raíz a través de un pelo radical. d. Sección transversal de una raíz mostrando pelotones dentro de las células del parénquima. e. Pelotones dentro de las células del parénquima. f. Colonias fúngicas aisladas. g. El tabique simple y la ramificación de la hifa en ángulo de 90 ° de Ceratobasidium. h. Células monilioides de un aislamiento de Ceratobasidium. i. Células binucleadas de miembros de Tulasnellaceae. c, e, g, h e i tomadas a un aumento de 100X. d tomada a un aumento de 40X.

FIGURE 2 Maximum likelihood phylogram showing the phylogenetic placement of the obtained sequences belonging to Ceratobasidiaceae among the most similar sequences retrieved from GenBank as inferred from BLAST searches of the full ITS region. Tree was midpoint rooted. Only bootstrap supports above 70 % are shown. Blue arrows show the placement of the different isolates. / Filograma con máxima verosimilitud que muestra la posición filogenética de las secuencias obtenidas pertenecientes a las Ceratobasidiaceae entre las secuencias más similares obtenidas desde GenBank inferidas de las búsquedas BLAST de la región completa de ITS. El árbol se enraizó en el punto medio. Sólo se muestran soportes de bootstrap por encima del 70 %. Las flechas azules muestran la ubicación de los diferentes aislamientos.

FIGURE 3 Maximum likelihood phylogram showing the phylogenetic placement of the obtained Tulasnellaceae sequences among the most similar sequences retrieved from GenBank as inferred from BLAST searches of the full ITS region. Tree was midpoint rooted. Only bootstrap supports above 70 % are shown. Blue arrows show the placement of the obtained sequences. / Filograma con máxima verosimilitud mostrando la posición filogenética de secuencias de Tulasnellaceae obtenidas entre las secuencias más similares obtenidas desde GenBank inferidas de las búsquedas BLAST de la región completa de ITS. El árbol se enraizó en el punto medio. Sólo se muestran soportes de bootstrap por encima del 70 %. La flecha azul muestra la ubicación de las secuencias obtenidas de los aislamientos.
TABLE 1 Identity of mycorrhizal fungi associated with roots of Codonorchis lessonii individuals from different Chilean populations determined by similarity of their ITS region with previous GenBank accessions. / Identidad de los hongos micorrícicos asociados a las raíces de individuos de Codonorchis lessonii de diferentes poblaciones chilenas, determinada por la similitud de la región ITS con accesiones previas de GenBank.


TABLE 2 Dasometric characteristics of the fungal isolates from three Chilean populations of Codonorchis lessonii. Different letters correspond to statistical differences (Tukey test, p < 0.05). / Características dasométricas de hongos aislados de tres poblaciones chilenas de Codonorchis lessonii. Letras diferentes corresponden a diferencias estadísticas (test de Tukey, p < 0,05).

(+): Low frequency of monilioid cells. (++): Moderate frequency of moniliod cells. (+++): High frequency of moniliod cells.
The survival and distribution of an obligate mycorrhizal plant might depend on the distribution of its fungal partners in its habitat (Pandey et al. 2013). The specificity of the fungal-orchid association can vary in different conditions, and laboratory-based experiments do not always reflect the mycorrhizal associations that take place in the wild (Porras-Alfaro & Hayman 2007, Steinfort et al. 2010). Moreover, different fungi might be required in different stages of the orchid life cycle (Rasmussen 2002). Our results showed that the fungi associated with the roots of C. lessonii individuals from different populations are related to other strains found in different hosts and distant geographical areas (Table 1, Fig. 2, 3) more than to each other. Low fungal specificity of C. lessonni could explain their abundance in certain areas, since they are not restricted to only one fungal strain to survive and reproduce like some rare orchids that may be limited by specific fungal species availability (Graham & Dearnaley 2012, McCormick et al. 2012, Swarts et al. 2010). However, it is unknown if specific fungi are required for the natural germination of this species. Recent studies suggest that the distribution of orchids could be limited by the joint effect of suitable habitat availability and pollination limitation (McCormick & Jacquemyn 2014). Yet another possible factor limiting orchid distributions and acting as a selecting agent could be the orchid-mycorrhizal relationships (Brundrett 2007, Karol et al. 2015, Otero et al. 2004). Further studies with a higher number of individuals per population at different life cycle stages and reciprocal germination and inoculation experiments using different fungal strains are needed to broaden the knowledge about the specificity of these associations. Knowledge of this kind, could help refine laboratory protocols and propagation methods aimed at re-introduction of this orchid species in the field, in future restoration and conservation programs.
Orchid mycorrhizae are excellent models to reveal the general properties of mycorrhizal systems as well as specificity aspects such as ecological networks and evolution of the mycorrhizal status (Dearnaley et al. 2012). In the case of C. lessonii, two possibilities arise; 1) if at the adult stage this orchid show low fungal specificity, at least within a certain phylogenetic range, other factors might limit its distribution and abundance; while, 2) if fungal specificity exists at some stage, specific fungi will be required for the laboratory reproduction of individuals from different populations, and the species may be locally restrained by the availability of a suitable fungal partner. This specificity may change during ontogeny but, to our knowledge, no studies have addresses this issue. Further studied are needed to confirm these hypotheses, but our results suggests that C. lessonii associate with phylogenetically distant fungi, at least at the adult stage.
The conservation status of C. lessonii is currently unevaluated. This species usually inhabits the understory of forests dominated by Nothofagus species (Freuler 2006, Vidal & Reif 2011, Vidal et al. 2012) which in their natural state are disappearing due to land use change and anthropogenic activities. This orchid is intolerant to high light exposure and after thinning and clearing of forests, their presence can be severely reduced. Thus, the identification, isolation, and culture of the associated mycorrhizal fungi could be an essential element in their conservation in the long-term and in the understanding of their ecology and evolution.