Bosiite, NaFe3+3(Al4Mg2)(Si6O18)(BO3)3(OH)3O, is a new mineral species of the tourmaline supergroup from the Darasun gold deposit (Darasun mine), Vershino-Darasunskiy, Transbaikal Krai, Eastern-Siberian Region, Russia (52°20′24″N, 115°29′23″E). Bosiite formed as a hydrothermal phase in a gold-bearing quartz-vein spatially related to the Amudzhikan–Sretensky subvolcanic K-rich granodiorite-porphyry intrusion. Ores of this deposit are enriched in sulfides (up to 60%). Bosiite is intimately associated with other tourmalines. The first tourmaline generation is bosiite, which is followed by a second generation of oxy-dravite and a third generation of dravite. Bosiite also coexists with quartz and pyrite; further associated minerals in the vein are gangue minerals (quartz, calcite, and dolomite), sulfides (pyrite, arsenopyrite, chalcopyrite, pyrrhotite, tetrahedrite, sphalerite, and galena) and native gold. Crystals of bosiite are dark brown to black with a pale-brown streak. Bosiite is brittle and has a Mohs hardness of 7; it is non-fluorescent, has no observable parting and cleavage. It has a measured density of 3.23(3) g/cm3 (by pycnometry) and a calculated density of 3.26(1) g/cm3. In plane-polarized light, it is pleochroic, O = yellow-brown, E = red-brown. Bosiite is uniaxial negative, ω = 1.760(5), ϵ = 1.687(5). The mineral is trigonal, space group R3m, a = 16.101(3), c = 7.327(2) Å, V = 1645.0(6) Å3. The eight strongest X-ray diffraction lines in the (calculated) powder pattern [d in Å(I)hkl] are: 2.606(100)(50-1), 8.051(58)(100), 3.008(58)(3-1-2), 4.025(57)(4-20), 3.543(50)(10-2), 4.279(46) (3-11), 2.068(45)(6-1-2), 4.648(28)(300). Analysis by a combination of electron microprobe (EMPA), inductively coupled plasma mass spectrometry (ICP-MS), Mössbauer spectroscopic data and crystal-structure refinement results in the empirical structural formula:
According to the IMA-CNMNC guidelines, the dominant valence at the Y site is 3+ and the dominant cation is Fe3+. To accommodate the disorder and allocating cations to the Z and Y sites, the recommended procedure leads to the optimized empirical formula (based on 31 O): X(Na0.73Ca0.230.04) Y(Fe3+2.40Fe2+0.59Ti4+0.01)Z(Al3.36Mg1.69Fe3+0.95) T(Si5.92Al0.08O18) (BO3)3V(OH)3W[O0.85(OH)0.15]. Bosiite, ideally NaFe3+3(Al4Mg2)(Si6O18)(BO3)3(OH)3O, is related to end-member povondraite, ideally NaFe3+3(Fe3+4Mg2)(Si6O18)(BO3)3(OH)3O, by the substitution ZAl4ZFe3+4. Further, bosiite is related to oxy-dravite, ideally Na(Al2Mg)(Al5Mg)(Si6O18)(BO3)3(OH)3O, by the substitutions [6]Fe3+3[6]Al3.

Bosiite is named after Dr. Ferdinando Bosi, researcher at the University of Rome La Sapienza, Italy, and an expert on the crystallography and mineralogy of the tourmaline-supergroup minerals and the spinels.

The general chemical formula of the tourmaline-super-group minerals can be written XY3Z6[T6O18] (BO3)3V3W, as proposed by Henry et al. (2011). These authors and Hawthorne (1996, 2002) suggest occupancies by the following most common cations:
Crystal-chemical relations in the tourmaline supergroup and the crystal chemistry of tourmaline-supergroup minerals have been investigated by many authors in the last 40 years (e.g., Povondra & Čech, 1976; Foit & Rosenberg, 1979; Deer et al., 1986; Foit, 1989; Hawthorne et al., 1993; Hawthorne, 1996, 2002; Henry & Dutrow, 1996; Ertl et al., 1997, 2012a, Bosi & Lucchesi, 2004, 2007; Bosi et al., 2004, 2005, 2013, 2014; Hughes et al., 2004, 2011; Henry et al., 2011; Ertl & Tillmanns, 2012; Reznitskii et al., 2014; 2015, 2016). Substitutions at the Y site are more constrained by size than by valence and hence the variety of cations is larger than at the Z site, where substitutions are more constrained by valence than by size (Grice & Ercit, 1993).

Tourmaline, rich in Fe3+, has been recorded from many different localities: from the type locality of povondraite (Žáček et al., 1998), from the cap rock of a salt dome in the Gulf of Mexico (Henry et al., 1999), and from evaporite formations in Namibia (Henry et al., 2008). Tourmalines from the Larderello geothermal field, Italy (Cavaretta & Puxeddu, 1990), from northern Scandinavia, and from porphyry-copper deposits formed during hydrothermal processes (Baksheev et al., 2011). The origin of Al-poor tourmaline from the Western Carpathians is still not clear (Bačík et al., 2008). Tourmaline from the Darasun gold deposit is, except for the Mg content, similar to tourmaline formed in evaporites and similar to tourmaline from the Western Carpathians, which are also enriched in Ca, but have a different compositional evolution: from Fe3+-rich tourmaline (bosiite) to Fe-bearing dravite in Darasun, whereas in the Western Carpathians it evolved from schorl–dravite to povondraite. Samples from Darasun are very similar in composition to those from the Vetka porphyry-copper deposit, Chukotka, Russia, and also in Mg and Ca contents (Baksheev et al., 2010, 2011). The maximum Fetot content of tourmaline from Darasun (~5 apfu) is higher than in tourmalines from deposits in porphyry-copper deposits (~3 apfu) (Lynch & Ortega, 1997). In addition, the Fe3+/Fetot value in tourmalines from the Darasun gold deposit is significantly higher. Tourmaline from porphyry tin deposits also displays a homovalent AlFe−1 substitution, but is characterized by a lower Fe3+/Fetot value (Baksheev et al., 2010, 2011, 2012). The latter authors concluded that the homovalent Fe3+→Al substitution in tourmaline (also observed in tourmalines from Darasun gold deposit) may be used as prospecting guide for porphyry-style deposits and that tourmalines from Cu, Sn and Au porphyry-type deposits can be distinguished from each other on the basis of Li concentration and the ratio Fe3+/Fetot.

The Darasun gold deposit (Darasun mine), Vershino-Darasunskiy, Transbaikal Krai, Eastern-Siberian Region, Russia (52°20′24″N, 115°29′23″E), was the second mine in terms of gold reserves of the Transbaikal Krai. This deposit is hosted in basic and intermediate igneous rocks. The most ancient rock is metamorphosed Early Palaeozoic gabbro, which was intruded by Middle Palaeozoic and Middle Mesozoic granodiorite, diorite, granite, granosyenite and syenite, and also by granitic rocks of the Amanansky complex of Early Jurassic age. The formation of gold deposits in Transbaikalia seems to be controlled by small intermediate and felsic subvolcanic intrusions enriched in alkalis, which belong to the Middle to Late Jurassic Amudzhikan–Sretensky complex and comprise granite and granodiorite porphyry, fracture-related intrusions and stocks, quartz diorite porphyry, granophyre, felsite and pitchstone dykes (Baksheev et al., 2011). The Darasun deposit comprises extensive steeply dipping gold-bearing quartz veins and mineralized zones spatially related to the Amudzhikan–Sretensky subvolcanic K-rich granodiorite porphyry intrusion. Ores of this deposit are rich in sulfides (up to 60%). Gangue minerals are quartz, tourmaline, calcite and dolomite.

Bosiite was found in a sulfide- and gold-bearing quartz vein. Native gold, calcite, dolomite, arsenopyrite, chalcopyrite, pyrrhotite, tetrahedrite, sphalerite and galena occur in the same vein. Bosiite coexists with quartz and pyrite (Fig. 1 and Fig. S1 in Supplementary Material, linked to this article and freely available online at the GSW website of the journal: and is intimately associated with other tourmalines (Figs 2 and 3). The first tourmaline generation is bosiite, which is followed by a second generation (Fig. 4) of a new tourmaline (recently described as oxy-dravite by Bosi & Skogby, 2013), and a third generation, which is represented by dravite (~20 μm thick veins). Depletion of Fe in the second and third generations is caused by deposition of early sulfides.

Baksheev et al. (2011) have studied 73 primary fluid inclusions in quartz associated with tourmaline from the bosiite type locality. At room temperature, they found vapour-dominant and aqueous inclusions, indicating boiling in the fluid. Microthermometry indicated that vapour-dominant inclusions homogenize at ~360–380°C and contain a fluid with a salinity of about 1.7 wt% equiv. NaCl. Aqueous inclusions homogenize to liquid at ~350–390°C and contain fluid with a salinity ranging from ~7 to 10 wt% equiv. NaCl. The calculated pressure from vapour-dominated inclusions is 180–230 bars, which corresponds to a depth of 600–800 m. These authors concluded that the fluid inclusions in the associated quartz indicate that first-stage tourmaline (bosiite) precipitated from boiling fluid at low pressure. Boiling favoured the increasing oxygen activity (causing the loss of H2 into the vapour phase), and resulted in the increase of Fe3+ in that tourmaline. Based on these data, bosiite crystallized probably during a hydrothermal event at PT conditions of ~0.02 GPa/380°C.

The holotype specimen of bosiite (crystal used for the structural study and coarse-grained sample) is deposited in the collection of the Naturhistorisches Museum Wien, Vienna, Austria, under the catalogue number N 9793. Parts of the cotype material (coarse-grained sample) have also been deposited in the collection of the Fersman Mineralogical Museum, Moscow, Russia (catalogue no. 94999).

4.1. Crystal-structure refinement

A fragment of bosiite was separated from the same hand specimen from which the chemistry had been determined; it was subsequently studied on a Bruker AXS Kappa APEXII single-crystal diffractometer equipped with a charge-coupled device (CCD) area detector and an Incoatec Microfocus Source IμS (30 W, multilayer mirror, MoKα). Single-crystal X-ray diffraction data were collected at room temperature (up to 80° 2θ), integrated and corrected for Lorentz and polarization factors with an absorption correction by evaluation of partial multiscans. The structure was refined with SHELXL97 (Sheldrick, 1998) using scattering factors for neutral atoms and a tourmaline starting model from Ertl et al. (2010). The H atom bonded to the O3 atom was located from a difference-Fourier map and subsequently refined. Refinement was done with anisotropic displacement parameters for all non-H atoms. Table 1 provides crystal data and details of the structure refinement. Site occupancies were refined according to well-known characteristics of the tourmaline structure (Na and Ca were refined at the X site, Fe and Al were refined at the Y and Z sites; for other details, see Table 2). The refinement converged at a R1(F) value of ~3.8% (Table 1). Atom parameters and equivalent isotropic-displacement parameters are given in Table 2, selected interatomic distances in Table 3.

4.2. Chemical analysis

A tourmaline fragment was prepared as a section (polished on one side of the sample) for chemical analysis. Concentrations of all elements except B, Li, Be and H were determined with a JEOL JSM–6480LV electron microscope equipped with an Inca Energy–350 energy-dispersion system (EDS) (Laboratory of Analytical Techniques of High Spatial Resolution, Department of Petrology, Lomonosov Moscow State University, Moscow, Russia). Chemical analyses were done by a combination of EDS Inca Energy–350 and WD (wavelength-dispersion) Inca Wave–500 (accelerating potential 15 kV, 64 ± 0.1nA, ~5 μm beam diameter). The concentration of F (MgF2 was used as a standard) was measured with a WD TAP crystal and the content of the other elements with EDS. For the EDS analyses, line profiles of characteristic X-ray radiation were optimized and normalized using natural silicates as standards (Jarosewich et al., 1980). The error for elements with a concentration of ≥5 wt% does not exceed 1.5 relative %.

The Li content was determined by ICP–MS (Institute of Geology of Ore Deposits, Petrography, Mineralogy and Geochemistry (IGEM), Russian Academy of Sciences, Moscow, Russia). Initially, a sample was decomposed with H2SO4 and HF in a MILESTONE microwave device according to a four-step program. After the solution was evaporated, treatment by water and dissolution in HNO3 followed. The measurements were made by using a PLASMA QUAD mass spectrometer (VG Instruments Company).

The crystal-chemical formula was calculated on the basis of 31 (OH, O, F) atoms. The B2O3 content was calculated to fill the B site completely. Because of the small amount of available sample, H2O was not analysed, but was calculated for a charge-balanced formula assuming that there are no vacancies at the Y, Z, and T sites. Analytical data are given in Table 4. Back-scattered electron images were obtained at an accelerating potential of 15 kV and a current intensity of ~15 nA.

4.3. Mössbauer spectroscopy

Approximately 5 mg of tourmaline, consisting mainly of bosiite, were gently crushed under acetone, then mixed with a sugar–acetone solution designed to form sugar coatings around each grain and prevent preferred orientation. The 57Fe spectrum was acquired at 295 K using a SM–2201 Mössbauer spectrometer operating in constant acceleration mode with a 57Co (in Cr) source (Institute of Mineralogy, Ural Division, Russian Academy of Sciences, Miass, Russia). Calibration was done with a standard sample of sodium nitroprusside. Isomer shift refers to an α-Fe absorber at 293 K. Results of the measurements were processed by a least-squares procedure using the SPECTR program in fitting a thin absorber (using a Lorentzian line-shape).

Peak areas were not corrected for differential recoil-free fractions for Fe2+ and Fe3+ because the appropriate correction factors do not exist.

5.1. Optical properties and infrared spectroscopy

In plane-polarized transmitted light, bosiite is pleochroic, O = yellow-brown, E = red-brown. Bosiite is uniaxial negative, ω = 1.760(5), ϵ = 1.687(5). An infrared spectrum of a first-generation tourmaline single crystal (bosiite) was recorded by Baksheev et al. (2011, their Fig. 3B). These authors concluded that the spectrum indicates that a part of Fe3+ occurs at the Z site. The interactions of three OH groups at the V sites with one Y-site cation and two Z-site cations are indicated by a band in the range 3580–3480 cm−1, with a complex asymmetric profile. In the single-crystal absorption spectrum of bosiite, this band is located at a lower frequency with a maximum at 3552 cm−1. The interaction of the hydroxyl groups at the W site with three Y-site cations is represented by a single weak band at 3733 cm−1 (Baksheev et al., 2011).

5.2. Physical properties

Bosiite is brittle and has a Mohs hardness of 7; it is non-fluorescent, has no observable parting and cleavage, and has a vitreous lustre. The megascopic colour of bosiite is black and the streak is dark brown. It has a measured density of 3.23(3) g/cm3 (by pycnometry) and a calculated density of 3.26(1) g/cm3 from the empirical formula. The fracture is irregular to uneven and sub-conchoidal. It has a prismatic habit and only unspecified prismatic forms are developed. Twinning was not observed. The c:a ratio calculated from the refined unit-cell parameters is 0.455.

5.3. Mössbauer spectrum

The Mössbauer spectrum of bosiite (Fig. 5 in Baksheev et al., 2011) was fitted to three doublets (Table 5) that were assigned according to Fuchs et al. (1998); Béziat et al. (1999) and Dyar et al. (1998). The first two doublets ([Y1]Fe2+, [Y2]Fe2+) are subcomponents of a distribution corresponding to Fe2+ at the Y site and cover a total area of about 15%. This is also in agreement with the findings of Andreozzi et al. (2008), who modelled Fe2+ at the Y site in tourmalines of the elbaite–schorl–dravite series by three doublets (ΔEQ = 2.45, 2.19, and 1.72 mm/s). The third strongest doublet, with an isomer shift of 0.57 mm/s, is about 85%, implying that most of the iron in the tourmaline sample examined occurs as Fe3+ in octahedral coordination (Baksheev et al., 2011). However, the real amount of Fe3+ can be higher or lower because of the chemical heterogeneity of the studied sample, of the significant overlap of the peaks in the spectrum and of the additional uncertainty introduced in assuming the same recoil-free fraction for Fe2+ and Fe3+. As was shown in previous studies (e.g., Dyar et al., 1998; Andreozzi et al., 2008), the assignment of Fe3+ to Y or Z sites on the basis of Mössbauer data alone is problematic and sometimes not possible. Because of the limited resolution of our Mössbauer data (Table 5), we are not able to use it to establish whether Fe3+ occupies both the Y and Z sites. It is possible that the two expected Fe3+ doublets are hidden in a single doublet.

5.4. Crystal chemistry and structure analysis

The site populations at the Y and Z sites have been calculated assuming that Fe2+ and the small amount of Ti4+ occupy the Y site, using the method of Wright et al. (2000), although we cannot exclude that small amounts of Fe2+ and/or Ti4+ also occupy the Z site (e.g., Bosi, 2008; Ertl et al., 2012b). The amounts of Al, Fe3+ and Mg over the Y and Z site were optimized (Wright et al., 2000) according to mean atomic number and bond valences. Chemical characterization based on electron-microprobe analysis, ICP-MS, single-crystal structure refinement and Mössbauer spectroscopy resulted in the empirical structural formula:
Bosiite is trigonal, with space group R3m, a = 16.101 (3), c = 7.327(2) Å, V = 1645.0(6) Å3, Z = 3. The eight strongest X-ray diffraction lines in the (calculated) powder pattern [d in Å (I)hkl] are: 2.606(100)(50 - 1), 8.051(58)(100), 3.008(58)(3-1-2), 4.025(57)(4 - 20), 3.543(50)(10 - 2), 4.279(46)(3 - 11), 2.068(45)(6-1-2), 4.648(28)(300) (Table 6). The powder pattern had to be calculated because only a very small amount of bosiite was available, and it was intimately intergrown with oxy-dravite. The Flack parameter of the single-crystal structure shows that no racemic twinning is present (Table 1).

5.5. X-site occupancy

The X site in bosiite is occupied mainly by Na (0.73 apfu), but contains also significant amounts of Ca (0.23 apfu) (Table 4). The content of K is at the detection limit (Table 4). The refinement is in excellent agreement with the chemistry (Na0.75(2)Ca0.25(2); Table 2).

5.6. Y- and Z-site occupancy

The Y site is occupied by Fe3+ (~1.5 apfu), Mg (~0.8 apfu), Fe2+ (~0.6 apfu) and Al (~0.1 apfu). Small amounts of Ti4+ (0.01 apfu) might also occupy this site. Only traces of Li occupy the Y site (Li content is ~3 ppm; Table 4). The Y-site occupancy is consistent with the refinement results (Table 7) and is also in good agreement with the relatively large <Y–O> distance of 2.054 Å (Table 3).

The Z site is mainly occupied by Al (~3.2 apfu), but contains also significant amounts of Fe3+ (~1.9 apfu) and Mg (~0.9 apfu). The Z-site occupancy is consistent with the refinement results (Table 7) and is also in good agreement with the enlarged <Z–O> distance of 1.969 Å (Table 3). Hawthorne (1996) noted that the occurrence of Al3+ at the Y site and Mg2+ at the Z site is due to disorder rather than differences in chemical composition, and that the order–disorder reaction can be expressed as YMg + ZAl↔ YAl + ZMg [YMgZAl(YAlZMg)−1]. This disorder is driven by the short-range requirements of O2− at the W site (=O1 site) that is associated with 2Al + Mg configurations at the coordinating Y sites so as to satisfy its bond-valence requirements (Hawthorne, 1996, 2002). Taylor et al. (1995) showed that disorder of Al and Mg over the Y and Z sites in a calcic tourmaline occurs in tandem with significant O2− at the W site. Also in bosiite, the W site is mainly occupied by O2− (Table 4). All tourmalines from the type locality of bosiite have an Mg content of ~2 apfu, the W site is dominated by O2−, and chemical variations are in accord with the FeAl−1 or AlOFe−1OH−1 exchange vector (or both). In the suite, the composition of bosiite is substantially deficient in Al (<6 apfu), consistent with a predominance of the FeAl−1 exchange (see also Baksheev et al., 2011).

5.7. T-site occupancy

The T site in bosiite is mainly occupied by Si and contains relatively small amounts of Al (~0.1 apfu), consistent with the slightly enlarged <T–O> distance of 1.621 Å (Table 3) and data given by MacDonald & Hawthorne (1995) on the Si ↔ Al substitution in tourmaline. There are no indications for significant amounts of [4]B (the <T–O> distance is not reduced). Also the amount of YAl is too small (~0.1 apfu) for significant amounts of [4]B to occur through short-range order (Fig. 1 in Ertl et al., 2008).

5.8. V- and W-site occupancy

An H atom (H3) at the site associated with O3 was easily located in the refinement. Ertl et al. (2002) showed that the bond-angle distortion (σoct2) of the ZO6 octahedron in a tourmaline is largely a function of the <Y–O> distance in that tourmaline, although the occupant of the O(3) site (V site in the general formula) also affects that distortion. The correlation coefficient, r, of <Y–O> and σoct2 of the ZO6 octahedron is −0.99 (Fig. 2 in Ertl et al., 2005) for all investigated tourmalines whose V site is occupied by 3 (OH) groups. Bosiite (with Zσoct2 = 42.47 and <Y–O> 2.054 Å; Table 3) lies exactly on the V site = 3 (OH) line. Hence, the V site of bosiite is completely occupied by (OH)3. No H associated with the O(1) site (W site) was found by refinement. Hence, this site is clearly dominated by oxygen, as indicated by the chemical composition. This is also consistent with the infrared spectrum, which shows only a weak band for the hydroxyl groups at the W site (see above).

According to the IMA-CNMNC guidelines (Henry et al., 2011), the X site is dominated by Na, the Y site is dominated by Fe3+, the Z site is dominated by Al along with Mg (valency-imposed double-site occupancy), the W site by O2−. Accordingly, this tourmaline is a new species.

In addition, to accommodate the Mg-Al disorder over Y and Z, we used the procedure recommended by Henry et al. (2013), which leads to the formula:
The derived formula is consistent with the end-member formula:

Bosiite is a Na-, Fe3+- and Mg-bearing tourmaline, and is isotypic with povondraite and other trigonal (R3m) members of the tourmaline group such as oxy-dravite (Bosi & Skogby, 2013). Their ideal chemical compositions are given in Table 8. For comparison, chemical analyses and selected physical properties of the valid minerals of the tourmaline supergroup (Henry et al., 2011) are given in Table 9. A triangular plot (Fe–Al–Mg) of tourmalines from the bosiite type-locality is shown in Fig. 5 and highlights the main compositional changes.

The refined formula is X(Na0.75Ca0.25)Σ1.00Y(Fe2.01Al0.99)Σ3.00Z(Al4.15Fe1.85)Σ6.00T(Si6O18)(BO3)3V(OH)3WO, in reasonable agreement with the optimized formula (1). Slight differences may reflect minor chemical zonation.

The name follows the tourmaline nomenclature recommendations approved by the IMA-CNMNC (Henry et al., 2011), which necessitate a new root name for a mineral of the alkali-subgroup 3 with the formula NaFe3+3(Al4Mg2)(Si6O18)(BO3)3(OH)3O. The new mineral is named after Dr. Ferdinando Bosi (b. 1967), researcher at the University of Rome La Sapienza, Italy, and expert on the crystallography and mineralogy of the tourmaline-supergroup minerals and the spinels.

Bosiite, ideally NaFe3+3(Al4Mg2)(Si6O18)(BO3)3(OH)3O, is related to end-member povondraite, ideally NaFe3+3(Fe3+4Mg2)(Si6O18)(BO3)3(OH)3O, by the substitution ZAl4ZFe3+4. Further, bosiite is related to oxy-dravite, ideally Na(Al2Mg)(Al5Mg)(Si6O18)(BO3)3(OH)3O, by the substitutions [6]Fe3+3[6]Al3. It can be distinguished from schorl, povondraite, fluor-buergerite, dravite, oxy-dravite or other dark tourmalines only by accurate determination of the chemical composition, preferably by a combination of single-crystal structure refinement and chemical analysis including Mössbauer spectroscopy. It may also be confused with luinaite-(OH) (IMA 2009–046) and its currently unnamed F-analogue, both representing monoclinically distorted variants (Cm) of the tourmaline structure. Relatively large lattice parameters (a ≥ 16.1, c ≥ 7.3 Å) and the occurrence in a hydrothermal environment could be preliminary indications for the occurrence of bosiite.

We thank V.O. Yapaskurt (Department of Petrology, Lomonosov Moscow State University, Moscow, Russia), S.A. Gorbacheva (IGEM, Russian Academy of Sciences, Moscow, Russia) and A. B. Mironov (Institute of Mineralogy, Ural Division, Russian Academy of Sciences, Miass, Russia) for analytical work. We are grateful to Andreas Wagner (Institut für Mineralogie und Kristallographie, Vienna, Austria) for taking photographs of bosiite. We sincerely thank Frank C. Hawthorne and an anonymous reviewer for constructive comments on this manuscript. This work was supported in part by Austrian Science Fund (FWF) project no. P23012-N19 and no. P26903-N19 to AE, by IGCP 540 and the Russian Foundation for Basic Research (project nos. 09–05–00697a and 09–05–12037ofi–m). The purchase of the electron microscope was financially supported by the Program for the Moscow State University Development.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial Share Alike License (, which permits unrestricted non-commercial use, distribution and reproduction in any medium provided that the original work is properly cited and all further distributions of the work or adaptation are subject to the same Creative Commons License terms.

Supplementary data